Skip to main content
Advertisement
  • Loading metrics

RNA-Seq Reveals Infection-Induced Gene Expression Changes in the Snail Intermediate Host of the Carcinogenic Liver Fluke, Opisthorchis viverrini

  • Sattrachai Prasopdee,

    Affiliations Food-borne Parasite Research Group, Department of Parasitology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand, Centre for Biodiscovery and Molecular Development of Therapeutics, Australian Institute of Tropical Health and Medicine, James Cook University, Cairns, Queensland, Australia

  • Javier Sotillo,

    Affiliation Centre for Biodiscovery and Molecular Development of Therapeutics, Australian Institute of Tropical Health and Medicine, James Cook University, Cairns, Queensland, Australia

  • Smarn Tesana,

    Affiliation Food-borne Parasite Research Group, Department of Parasitology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand

  • Thewarach Laha,

    Affiliation Food-borne Parasite Research Group, Department of Parasitology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand

  • Jutharat Kulsantiwong,

    Affiliation Department of Biology, Faculty of Science, Udon Thani Rajabhat University, Udon Thani, Thailand

  • Matthew J. Nolan,

    Affiliation Centre for Biodiscovery and Molecular Development of Therapeutics, Australian Institute of Tropical Health and Medicine, James Cook University, Cairns, Queensland, Australia

  • Alex Loukas ,

    Contributed equally to this work with: Alex Loukas, Cinzia Cantacessi

    Affiliation Centre for Biodiscovery and Molecular Development of Therapeutics, Australian Institute of Tropical Health and Medicine, James Cook University, Cairns, Queensland, Australia

  • Cinzia Cantacessi

    Contributed equally to this work with: Alex Loukas, Cinzia Cantacessi

    cinzia.cantacessi@jcu.edu.au

    Affiliations Centre for Biodiscovery and Molecular Development of Therapeutics, Australian Institute of Tropical Health and Medicine, James Cook University, Cairns, Queensland, Australia, Department of Veterinary Medicine, University of Cambridge, Cambridge, United Kingdom

Abstract

Background

Bithynia siamensis goniomphalos is the snail intermediate host of the liver fluke, Opisthorchis viverrini, the leading cause of cholangiocarcinoma (CCA) in the Greater Mekong sub-region of Thailand. Despite the severe public health impact of Opisthorchis-induced CCA, knowledge of the molecular interactions occurring between the parasite and its snail intermediate host is scant. The examination of differences in gene expression profiling between uninfected and O. viverrini-infected B. siamensis goniomphalos could provide clues on fundamental pathways involved in the regulation of snail-parasite interplay.

Methodology/Principal Findings

Using high-throughput (Illumina) sequencing and extensive bioinformatic analyses, we characterized the transcriptomes of uninfected and O. viverrini-infected B. siamensis goniomphalos. Comparative analyses of gene expression profiling allowed the identification of 7,655 differentially expressed genes (DEGs), associated to 43 distinct biological pathways, including pathways associated with immune defense mechanisms against parasites. Amongst the DEGs with immune functions, transcripts encoding distinct proteases displayed the highest down-regulation in Bithynia specimens infected by O. viverrini; conversely, transcription of genes encoding heat-shock proteins and actins was significantly up-regulated in parasite-infected snails when compared to the uninfected counterparts.

Conclusions/Significance

The present study lays the foundation for functional studies of genes and gene products potentially involved in immune-molecular mechanisms implicated in the ability of the parasite to successfully colonize its snail intermediate host. The annotated dataset provided herein represents a ready-to-use molecular resource for the discovery of molecular pathways underlying susceptibility and resistance mechanisms of B. siamensis goniomphalos to O. viverrini and for comparative analyses with pulmonate snail intermediate hosts of other platyhelminths including schistosomes.

Author Summary

Despite recent significant advances in knowledge of the fundamental biology of the carcinogenic liver fluke Opisthorchis viverrini, little is known of the complement of molecular interactions occurring between this parasite and its prosobranch snail intermediate host, Bithynia siamensis goniomphalos. The determination of such interactions is a key, necessary component of the development of future integrated control strategies for liver fluke infection and associated bile duct cancer. Here, we use cutting-edge high-throughput sequencing technologies and advanced bioinformatic analyses to characterize, for the first time, qualitative and quantitative differences in gene expression between uninfected and O. viverrini-infected B. siamensis goniomphalos collected from an endemic region of Northeast Thailand. The analyses led to the identification of a number of molecules putatively involved in immune defense pathways against invading O. viverrini, and of key biological mechanisms potentially implicated in the ability of the parasite to successfully colonize its snail intermediate host. We believe that this ready-to-use molecular resource will provide the scientific community with new tools for the development of strategies to control the spread of liver fluke infection and the resulting bile duct cancer.

Introduction

Cholangiocarcinoma (CCA, or bile-duct cancer) is one of the most common malignant cancers in Southeast Asia, with incidence rates ranging from 93.8 to 317.6 per 100,000 people/year in some districts of Northeast Thailand alone [1], [2]. Chronic infections by the food-borne human liver fluke Opisthorchis viverrini have been associated with the occurrence of CCA in several Southeast Asian countries (including Thailand and Laos) [2], [3], which has led to the inclusion of this parasite amongst the Group 1 carcinogens by the International Agency for Research on Cancer [4]. Typically, CCA caused by O. viverrini is diagnosed 30–40 years after primary infection, and death occurs within 3–6 months post-diagnosis [5]. The life cycle of O. viverrini is complex, involving two intermediate hosts and multiple developmental stages [6]. Briefly, piscivorous definitive hosts, including humans, shed the embryonated eggs through the faeces; in aquatic environments, eggs are ingested by freshwater prosobranch snails of the genus Bithynia, in which the motile embryos ( = miracidia) hatch and develop into sporocysts. In the body of the first intermediate hosts, the sporocysts undergo asexual reproduction through the stages of rediae and cercariae. Subsequently, the motile cercariae are released into the aquatic environment where they actively seek the second intermediate hosts, i.e. cyprinid fishes (e.g. Puntius spp.). Upon penetrating the skin of these fishes, the parasites encyst as metacercariae under the skin and/or in muscular tissue. Humans acquire the infection by ingesting raw or undercooked infected fishes [5], [6]. Following ingestion, the metacercariae excyst in the duodenum and migrate as juvenile flukes to the intra-hepatic bile ducts, where they develop to adult hermaphrodite flukes within ∼4 weeks. The sexually mature adults release eggs that, via the bile, are excreted with the faeces into the environment [6].

Over the last decade, much effort has been directed towards the elucidation of the molecular mechanisms that determine the occurrence of CCA in O. viverrini-infected individuals [2], [7][9], and advances in high-throughput sequencing technologies and bioinformatics are contributing towards a better understanding of the systems biology of this parasite [10], [11]. While these advancements undoubtedly provide a foundation for the development of novel strategies for the diagnosis, treatment and prevention of O. viverrini-induced CCA [3], long-term control of this devastating disease strictly relies on the development of integrated approaches, targeting the parasite as well as its intermediate hosts [12], [13]. Thus far, large-scale molecular investigations of snail species acting as intermediate hosts of human parasites have mainly involved pulmonate species of the genus Biomphalaria, which harbor Schistosoma blood flukes (see Biomphalaria glabrata Genome Initiative at http://biology.unm.edu/biomphalaria-genome/index.html; [14], [15]) while, until recently, knowledge of the fundamental molecular biology of prosobranch snails has remained limited [12], [13]. Our group has recently utilized high-throughput (Illumina) sequencing and advanced bioinformatic analyses of sequence data to generate and characterize the first reference transcriptome of Bithynia siamensis goniomphalos (Gastropoda, Bithyniidae), the snail intermediate host of O. viverrini in areas of Northeast Thailand, Lao PDR, Cambodia and South Vietnam [13]. This sequence dataset represents a valid resource for molecular studies of these organisms and provides a solid basis for investigations of snail-parasite interactions at the molecular level. In the present study, we employ high-throughput (Illumina) sequencing and bioinformatic technologies to characterize differences in levels of gene transcription between uninfected and O. viverrini-infected B. siamensis goniomphalos and identify key molecules and biological pathways inferred to be associated with the processes of invasion and survival of the parasite in the snail intermediate host.

Materials and Methods

Procurement of snail materials and assessment of infection status

Adult B. siamensis goniomphalos snails were collected from natural freshwater bodies in Amphoe Muang, Khon Kaen Province, Thailand. The soft tissues of each snail were removed from the shells and washed three times with sterile distilled water. A portion of the hepatopancreas was isolated from each snail and subjected to DNA extraction using the Direct PCR kit (Thermo Scientific) according to manufacturer's instructions, while the remaining material was preserved in RNAlater (Sigma) for subsequent extraction of total RNA. A PCR targeting the O. viverrini pOV-A6 gene was employed to screen individual snails for infection by the liver fluke following the amplification protocol described by Sermswan et al. [16]. Amplicons were visualized on a 1.5% agarose gel stained in ethidium bromide and the presence of bands of approximately 330 bp in size was considered diagnostic for O. viverrini infection [17].

RNA isolation and Illumina sequencing

Two distinct pools of soft tissues from uninfected (UB; n = 4) and O. viverrini-infected snails (IB; n = 4) were subjected to isolation of total RNA using the Trizol reagent (Invitrogen), followed by DNAse I treatment (Promega) according to the manufacturer's instructions. The amounts and integrity of total RNA were determined using a 2100 BioAnalyzer (Agilent Technologies). Polyadenylated (PolyA+) RNA was purified from 10 µg of total RNA from each pool of uninfected and O. viverrini-infected snails using Sera-Mag oligo(dT) beads, fragmented to a length of 100–500 nucleotides and reverse transcribed to cDNA using random hexamers. The size-fractionated cDNA was end-repaired and adaptor-ligated according to the manufacturer's protocol (Illumina). Ligated products of ∼300 bp were excised from agarose gels and PCR-amplified (15 cycles) [18]. Products were cleaned using a MinElute column (Qiagen, Netherlands) and paired-end sequenced on an Illumina HiSeq 2000 [19] according to the manufacturer's protocol.

Bioinformatic analyses

Following removal of adapter sequences and sequences with suboptimal read quality (i.e., PHRED score of 32.0), the remaining 90 bp single-read sequences generated from the non-normalized cDNA libraries from UB and IB, as well as the raw sequence data generated previously [13] were pooled and assembled using the program Trinity [20] (http://trinityrnaseq.sourceforge.net/). Briefly, a representative set of transcripts was assembled using the software Inchworm; next, the Chrysalis software was used to cluster portions of alternatively spliced transcripts and/or unique portions of paralogous genes and construct a de Brujin graph for each cluster of transcripts. Finally, each de Brujin graph was analysed and alternatively spliced isoforms and transcripts derived from paralogous genes were resolved using the Butterfly software [20] (http://trinityrnaseq.sourceforge.net/). The non-redundant, assembled dataset was subsequently compared to transcriptomic sequence data available for O. viverrini [10], [11] as well as the entire genome sequences of Clonorchis sinensis, S. mansoni, S. japonicum and S. haematobium [10], [21][24]; sequences with high similarity to known O. viverrini and/or other molecules from parasitic trematodes (e-value cut-off: 1e-15) were excluded [25]. The non-redundant transcriptomic dataset for B. siamensis goniomphalos was then annotated using an established approach [23], [25]. Briefly, assembled contigs were compared using the BLASTn and BLASTx algorithms [26] to sequences available in the nucleotide sequence collection (nt) of NCBI (www.ncbi.nlm.nih.gov), and in the non-redundant (nr) (www.ncbi.nlm.nih.gov), Swiss-Prot (http://expasy.org/), Kyoto Encyclopedia of Genes and Genomes (KEGG; http://www.genome.jp/kegg/) and Clusters of Orthologous Groups (COG; http://www.ncbi.nlm.nih.gov/COG/) databases, respectively (e-value<0.00001), to identify putative homologues in other mollusks and organisms other than mollusks. In addition, the non-redundant B. siamensis goniomphalos transcriptome was compared, using a BLASTn algorithm, to an in-house built nucleotide sequence database containing 86,837 expressed sequence tags (ESTs) available for B. glabrata (http://www.ncbi.nlm.nih.gov/nucest/?term=biomphalaria+glabrata).

Proteins with the highest sequence similarity to conceptually translated B. siamensis goniomphalos transcripts were assigned Gene Ontology terms (GO; http://www.geneontology.org/) [27], according to the categories ‘Biological Process’, ‘Cellular Component’ and ‘Molecular Function’, using the Blast2GO software [28]. B. siamensis goniomphalos contigs with no homologues in available protein databases were conceptually translated using the ESTScan software and annotated using the same approach [29].

The short reads from each UB and IB dataset were then separately mapped to the non-redundant sequence data using SOAPAligner, such that each read was mapped to a unique transcript with a minimum alignment overlap of 40 nucleotides and a maximum of three mismatches per read. Reads that mapped to more than one transcript (“multi-reads”) were randomly allocated to a unique transcript, such that they were mapped only once. Comparative analyses of gene transcription levels between UB and IB were performed using the fragments per kilobase of exon per million fragments method (i.e., FPKM) [30]. The analysis of statistical difference of transcription was determined using a method developed for serial analysis of gene expression (SAGE) and applied to RNA-seq data [31]. Statistical significance was set at a p value of ≤0.01 and, to control for errors associated with multiple pairwise comparisons, a false-discovery rate (FDR) correction was applied to the data set [32]. GO and KEGG functional enrichment analyses for B. siamensis goniomphalos genes whose transcription levels were altered by infection with the parasite, were performed using the software GO::TermFinder with multiple testing correction using the Benjamini & Hochberg calculations of FDR (cut-off: FDR<0.01) [33], [34].

Results

A total number of 60,210,044 and 56,426,082 raw reads were generated from UB and IB, respectively; of these, 55,137,440 (92%) and 52,476,468 (92%), respectively, were retained for further analyses (Table 1). All single reads generated in the present study have been deposited in the Sequence Read Archive (SRA) database of NCBI (http://www.ncbi.nlm.nih.gov/sra) under accession numbers SRR1046838 (UB) and SRR1046839 (IB), respectively. The assembly yielded 43,259 contigs, of which 20,171 displayed significant sequence similarity (e-value cut-off: 1e-15) to publicly available sequences from O. viverrini, C. sinensis, S. mansoni, S. japonicum and/or S. haematobium (data not shown) and were therefore excluded from subsequent analyses. All remaining 23,088 contigs (average length = 775 nt±679; maximum length = 12,119 nt; minimum length = 200; N50 = 671; cf. Table 1), matched previously assembled sequence data for B. siamensis goniomphalos, and displayed significant sequence similarity (e-value cut-off: 1e-05) with available ESTs from B. glabrata (not shown). Of the assembled sequence dataset, 5,461 (23.6%) and 13,283 (57.5%) contigs had homologues in the non-redundant nucleotide (nt) and protein (nr) sequence databases of NCBI, respectively (e-value cut-off: 1e-05; Table 1). Of the conceptually translated contigs with significant similarity to sequences in the nr database, 18.4% matched known proteins from the bivalve mollusk Crassostrea gigas (Table S1). This Transcriptome Shotgun Assembly Project has been deposited at DDJB/EMBL/GenBank under the accession GAQQ00000000.

thumbnail
Table 1. Summary of the RNA-Seq data from uninfected and Opisthorchis viverrini-infected Bithynia siamensis goniomphalos snails (UB and IB, respectively), prior and following assembly, and results of bioinformatic analyses and annotation.

https://doi.org/10.1371/journal.pntd.0002765.t001

Of the 17,849 contigs that encoded a predicted protein (cf. Table 1), 10,761 (60.3%) had homologues in the Swiss-Prot database (cf. Table 1); these included mucins (n = 378), histones (n = 318), collagens (n = 251) and serine/threonine protein phosphatases (STPs; n = 181) and kinases (STKs; n = 119) (Table S1). A total of 5,950 predicted proteins (33.3%) could be assigned GO terms, while 9,497 (53.2%) matched homologous proteins in the KEGG database associated to 259 distinct biological pathways (Table 1). The 5,287 (29.6%) predicted proteins with matches in the COG database could be assigned to at least one of 25 functional categories, of which ‘general function prediction’ (84%), ‘translation, ribosomal structure and biogenesis’ (50.2%) and ‘post-translational modification, protein turnover and chaperones’ (41%) were the most represented (Tables 1 and S1).

Comparative analyses of levels of gene expression revealed 6,920 and 735 DEGs in UB and IB, respectively (FPKM≥2.0; Table 1, Figure 1). Of the entire complement of molecules up-regulated in UB, 4,488 (64%) did not match known proteins in the nr database (Table S2); of the remaining 2,462, 7% were mucins, followed by uncharacterized proteins (2%) and immunoglobulin receptor homologues (0.6%). Of the 208 DEGs in IB with known homologues in the nr database (cf. Table S2), 38% were uncharacterized proteins, followed by actins (1%) and transcription factors (0.7%) (Table S2). Two haemocyanins and a Kunitz-type protease inhibitor were amongst the most significantly up-regulated DEGs in UB, whereas a histone and dynein displayed the highest expression levels in IB (Table 2). GO-enrichment analysis of DEGs in UB and IB revealed significant enrichment of the GO terms ‘cellular process’ (80.8% of DEGs), ‘organic substance metabolic process’ (53.3%) and ‘cellular metabolic process’ (52.1%) within ‘biological process’, ‘cell’ (90.4%), ‘intracellular part’ (85.4%) and ‘organelle’ (68.7%) within ‘cellular component’ and ‘protein binding’ (28.3%), ‘hydrolase activity’ (25.1%) and ‘ribonucleoside binding’ within ‘molecular function’ (Figure 2a). KEGG pathway enrichment analysis identified a total of 43 pathways linked to DEGs, including ‘RNA transport’ (8.6%), ‘Huntington’s disease’ (7.1%) and ‘spliceosome’ (8.5%) (Figure 2b). Out of 375 DEGs predicted to be involved in biological pathways linked to immune responses (cf. Figure 2b), homologues of metalloproteases, and actins and heat shock proteins displayed the largest up-regulation in UB and IB, respectively (Figure 3).

thumbnail
Figure 1. Scatter plot comparing log2 ratios of FPKM expression values for uninfected and Opisthorchis viverrini-infected Bithynia siamensis goniomphalos snails (UB and IB, respectively).

https://doi.org/10.1371/journal.pntd.0002765.g001

thumbnail
Figure 2. (a) Enriched Gene Ontology (GO) terms assigned to differentially expressed genes (DEGs) in uninfected and Opisthorchis viverrini-infected Bithynia siamensis goniomphalos snails (UB and IB, respectively), according to the categories ‘biological process’, ‘cellular component’ and molecular function; (b) enriched KEGG biological pathways assigned to DEGs in UB and IB, respectively.

https://doi.org/10.1371/journal.pntd.0002765.g002

thumbnail
Figure 3. The top ten up-regulated genes linked to biological pathways associated to immune responses [based on orthology to molecules in the Kyoto Encyclopedia of Genes and Genomes (KEGG)] in each uninfected and Opisthorchis viverrini-infected Bithynia siamensis goniomphalos snails (UB and IB, respectively).

https://doi.org/10.1371/journal.pntd.0002765.g003

thumbnail
Table 2. The top-ten up-regulated genes in uninfected and Opisthorchis viverrini-infected Bithynia siamensis goniomphalos snails (UB and IB, respectively).

https://doi.org/10.1371/journal.pntd.0002765.t002

Discussion

Despite the high prevalence of infection with O. viverrini in humans in the Greater Mekong sub-region (i.e. >90% in lowland districts of Lao PDR), and the large number of cyprinid fishes which may act as potential second intermediate hosts for this parasite [35][37], snail infection rates are extremely low, i.e. typically <1% [37][41] with peaks of ∼3% recorded in wetlands of Thailand and Lao PDR [42]. These observations lead to the hypothesis that infections of Bithynia snails by O. viverrini may stimulate a cascade of immune-molecular events similar to those induced in Biomphalaria spp. in response to the invasion by schistosomes, which are aimed at eliminating and/or preventing the dissemination of the infection in the snail intermediate hosts [14], [15], [43][45]. While knowledge of the complex mechanisms underlying the activation of immune defense pathways in Biomphalaria snails is expanding [46], [47], studies of the immune-molecular interactions in the Bithynia-Opisthorchis system are lacking. By undertaking the first large-scale comparative analysis of gene expression in UB and IB, we hoped to provide the community with a ready-to-use molecular infrastructure for in-depth studies of biological pathways specifically involved in snail-parasite interactions. An obvious limitation of this study was represented by the fact that the snail specimens subjected to RNA-Seq were collected from natural freshwater bodies, and it was therefore impossible to establish whether the uninfected Bithynia snails had been previously exposed to infection by O. viverrini and/or other parasites, or whether a previous infection had successfully been cleared. However, unlike the monoecious Biomphalaria spp., the dioecious nature of Bithynia poses some challenges for the rearing and maintenance of this snail species under experimental conditions. Indeed, it has been reported that Bithynia spp. specimens reared in the laboratory are smaller and morphologically altered when compared to the wild counterparts [40], [48]; since changes in body morphology are likely to correspond to alterations in gene regulation and expression, both in the steady-state and following infection by O. viverrini, we preferred to investigate transcriptional differences between UB and IB collected in the field.

Despite displaying significant sequence similarities to ESTs from B. glabrata, a significant proportion of B. siamensis goniomphalos DEGs (65%) could not be functionally annotated based on homology to known protein sequence data available in public databases. This finding is likely to be related to the limited number of annotated proteins available from mollusks [49], which impairs meaningful functional analysis of large-scale transcriptomic datasets. For DEGs matching available protein sequences, GO and KEGG enrichment analyses identified a range of biological pathways linked to immune responses, e.g. T and B cell signaling pathways (cf. Figure 2b), consistent with the results of previous investigations of levels of gene transcription in Biomphalaria snails exposed to infection by schistosomes [14], [43][45], [50], [51]. Amongst the immunity-associated DEGs, two proteases, a metalloprotease and a cysteine protease (cathepsin L) displayed the largest degree of up-regulation in UB when compared with IB (cf. Figure 3). The role of proteolytic enzymes, such as cathepsins, in the defense of mollusks against invading trematodes (e.g. in the breakdown of components of invading parasites) has been well documented [52], [53]. For instance, enhanced proteolytic activity has been described from extracts from S. mansoni-resistant compared to susceptible B. glabrata [53], while other studies have consistently detected transcripts encoding proteases amongst the molecules up-regulated in schistosome-resistant snail lines [14], [44]. In particular, two distinct expressed sequence tags (ESTs) encoding cathepsins L had been previously identified as highly expressed in B. glabrata snails susceptible to infection by the trematode parasite Echinostoma caproni [54] and resistant to S. mansoni [14], respectively. Based on the results of bioinformatic analyses of these two ESTs, Lockyer et al. [14] hypothesized that these two sequences might have represented two non-overlapping fragments of the same gene, and that the discrepancies in gene expression which had been observed followed infection by E. caproni and S. mansoni might have been due to “the different response elicited by two different parasite species” [14]. Comparative analyses between the abovementioned Biomphalaria ESTs and the cathepsin L-encoding sequence identified in the present study suggest that the latter transcript is indeed distinct from the two previously identified (data available from the corresponding author upon request). Moreover, in addition to the transcript encoding a ‘novel’ cathepsin L, a molecule encoding a ficolin was identified amongst the most highly expressed transcripts in UB (cf. Figure 3). Ficolins are a group of lectin-like defense proteins that, in vertebrates, are known to recognize carbohydrate molecules on pathogens and/or apoptotic and necrotic cells [55]. Structurally, ficolins consist of a small N-terminal domain, a collagen-like domain and a fibrinogen-like domain similar to the C-terminal domains of the beta and gamma chains of fibrinogen [56]. These molecules, which belong to a large protein family known as Fibrinogen-Related Proteins (FREPs), play key roles in the innate immunity of both vertebrate and invertebrate organisms [57]. In particular, in mollusks, FREPs are believed to exert ‘antibody-like’ properties in biological mechanisms of recognition and binding of invading pathogens [58], [59]. In B. glabrata, transcription of FREP-encoding genes is enhanced after exposure of the snail to both S. mansoni and E. paraensei [43], [60], [61]; furthermore, FREPs have been shown to bind sporocysts of both S. mansoni and E. paraensei, thus actively participating in molecular events linked to defense against pathogens [62]. A mucin and an immunoglobulin receptor homologue were also amongst highly transcribed molecules in UB (cf. Figure 3). While the up-regulation of transcripts encoding mucins could be attributed to a role of these proteins in the establishment of a physical defense barrier against invading pathogens [63], the detection of an immune receptor-like molecule amongst the DEGs is interesting, since mollusks lack adaptive immunity. However, in a previous study, Bayne and co-workers [64] demonstrated that hemocytes of B. glabrata were unable to recognize IgG-bound antigens of S. mansoni sporocysts, thus providing evidence for the existence of specific antigen-binding sites on snail hemocytes [64]. Results from the present study support this hypothesis.

While it is tempting to speculate that the identification of homologues of a protease, a FREP and an immunoglobulin-receptor amongst the most highly expressed transcripts in UB (collected from the same natural water bodies as IB) may be associated with a ‘successful’ defense response against O. viverrini, which ultimately resulted in prevention or clearance of the infection (similarly to resistant strains of B. glabrata), such a hypothesis requires rigorous testing.

Amongst the DEGs that could be linked to pathways associated with immune responses, actins and heat-shock proteins (HSPs) were identified amongst the molecules with the largest degree of up-regulation in IB (cf. Figure 3). The function/s of HSPs in stress responses of hemocytes of B. glabrata in response to schistosome infection has been documented [65], [66]; conversely; the role/s played by actins in snail-parasite relationships is unclear. Since defense against pathogens in mollusks mainly relies on the migration of circulating hemocytes to the sites of infection and tissue damage, it could be hypothesized that the relative increase in transcription of actin-encoding genes might reflect a general increase in cellular motility, as well as phagocytosis and cell trafficking. Besides actins and heat-shock proteins, the absence of key immune-regulatory molecules amongst the most-significantly up-regulated genes in IB remains to be investigated; one possible explanation could be that, similarly to schistosome parasites infecting susceptible Biomphalaria spp. [14], O. viverrini manipulates the innate defense system of susceptible Bithynia, effectively resulting into the inability of the hemocytes to recognize the parasite and/or into the suppression of the snail humoral responses against the parasite invasion. This hypothesis is also supported by the fact that, out of the DEGs identified in this study, the vast majority (i.e. ∼90%) were molecules whose levels of transcription were up-regulated in UB, thus reinforcing a potential role of O. viverrini in mechanisms of immune suppression in susceptible snails [14].

Taken together, the similarities in gene expression profiling between parasite-free specimens of wild-caught B. siamensis goniomphalos examined in the present study and schistosome-resistant populations of B. glabrata, and between infected Bithynia and schistosome-susceptible Biomphalaria, lead to the speculation that an O. viverrini-resistant population of Bithynia snails may occur in natural freshwater bodies of Northeast Thailand. However, it must be pointed out that the snail specimens herein examined were screened for the presence of O. viverrini DNA only; therefore, the occurrence of underlying infections by other snail pathogens in both UB and IB cannot be excluded. Ultimately, in order to test the hypothesis that O. viverrini-resistant B. siamensis goniomphalos may occur in Thailand, the development of adequate and efficient methods to successfully rear and maintain populations of this snail species under laboratory conditions is necessary. Once an optimal rearing protocol is established, future studies could investigate changes in gene expression of key genes at different time points following experimental infection of Bithynia specimens with O. viverrini; such studies will allow to further explore the role/s of selected molecules, such as proteases and actins, in the cascades of molecular events that determine susceptibility or resistance of snails to trematode parasites.

A large collaborative project aimed at generating the first reference whole genome sequence of B. glabrata is currently underway (http://biology.unm.edu/biomphalaria-genome/). Once completed, the B. glabrata genome will represent an unprecedented resource for in-depth studies of snail-trematode interactions, as well as for comparative analyses between and among species acting as intermediate hosts of parasitic flukes, and for the discovery of the molecular and evolutionary mechanisms associated to the selection of parasite-resistant snail populations [67], [68]. Similarly, the wealth of molecular data that is currently available for O. viverrini [10], [11], [69], together with future efforts to define its entire genome sequence, will set the basis for the identification of genes and gene products essential for the recognition and infection of susceptible snail hosts. Taken together, these advances will provide the scientific community with new tools for the development of strategies to control the spread of opisthorchiasis and Opisthorchis induced-CCA.

Supporting Information

Table S1.

The annotated transcriptome of Bithynia siamensis goniomphalos and corresponding sequence data. Summary of annotation data linked to individual assembled transcripts. Abbreviations: non-redundant protein database (Nr), non-redundant nucleotide database (Nt), SwissProt database (SwissProt), Clusters of Orthologous Groups of proteins (COG), Kyoto Encyclopaedia of Genes and Genomes (KEGG); Gene Ontology (GO).

https://doi.org/10.1371/journal.pntd.0002765.s001

(ZIP)

Table S2.

Differentially expressed genes (DEGs) inferred from the transcriptomes of uninfected and Opisthorchis viverrini-infected Bithynia siamensis goniomphalos snails (UB and IB, respectively). Abbreviations: fragments per kilobase of exon per million fragments (FPKM), false discovery rate (FDR), non-redundant protein database (Nr), non-redundant nucleotide database (Nt), SwissProt database (SwissProt), Clusters of Orthologous Groups of proteins (COG), Kyoto Encyclopaedia of Genes and Genomes (KEGG); Gene Ontology (GO).

https://doi.org/10.1371/journal.pntd.0002765.s002

(XLS)

Author Contributions

Conceived and designed the experiments: SP AL CC. Performed the experiments: SP JS. Analyzed the data: SP JS CC TL JK MJN. Contributed reagents/materials/analysis tools: ST AL. Wrote the paper: SP AL CC.

References

  1. 1. Sriamporn S, Pisani P, Pipitgool V, Suwanrungruang K, Kamsa-ard S, et al. (2004) Prevalence of Opisthorchis viverrini infection and incidence of cholangiocarcinoma in Khon Kaen, Northeast Thailand. Trop Med Int Health 9: 588–594.
  2. 2. Sithithaworn P, Andrews RH, Petney TN, Saijuntha W, Laoprom N (2012) The systematics and population genetics of Opisthorchis viverrini sensu lato: implications in parasite epidemiology and bile duct cancer. Parasitol Int 61: 32–37.
  3. 3. Smout MJ, Sripa B, Laha T, Mulvenna J, Gasser RB, et al. (2011) Infection with the carcinogenic human liver fluke, Opisthorchis viverrini. Mol Biosyst 7: 1367–1375.
  4. 4. Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, et al. (2009) A review of human carcinogens – Part B: biological agents. Lacent Oncol 10: 321–322.
  5. 5. Sripa B, Kaewkes S, Sithithaworn P, Mairiang E, Laha T, et al. (2007) Liver fluke induces cholangiocarcinoma. PLoS Med 4: e201.
  6. 6. Kaewkes S (2003) Taxonomy and biology of liver flukes. Acta Trop 88: 177–186.
  7. 7. Smout MJ, Laha T, Mulvenna J, Sripa B, Suttiprapa S, et al. (2009) A granulin-growth factor secreted by the carcinogenic liver fluke, Opisthorchis viverrini, promotes proliferation of host cells. PLoS Pathog 5: e1000611.
  8. 8. Daorueang D, Thuwajit P, Roitrakul S, Laha T, Kaewkes S, et al. (2012) Secreted Opisthorchis viverrini glutathione S-transferase regulates cell proliferation through AKT and ERK pathways in cholangiocarcinoma. Parasitol Int 61: 155–161.
  9. 9. Chusorn P, Namwat N, Loilome W, Techasen A, Pairojkul C, et al. (2013) Overexpression of microRNA-21 regulating PDCD4 during tumorigenesis of liver-fluke associated cholangiocarcinoma contributes to tumor growth and metastasis. Tumor Biol 34: 1579–1588.
  10. 10. Young ND, Campbell BE, Hall RS, Jex AR, Cantacessi C, et al. (2010) Unlocking the transcriptomes of two carcinogenic parasites, Clonorchis sinensis and Opisthorchis viverrini. PLoS Negl Trop Dis 4: e719.
  11. 11. Jex AR, Young ND, Sripa J, Hall RS, Scheerlinck JP, et al. (2012) Molecular changes in Opisthorchis viverrini (Southeast Asian liver fluke) during the transition from the juvenile to the adult stage. PLoS Negl Trop Dis 6: e1916.
  12. 12. Adema CM, Bayne CJ, Bridger JM, Knight M, Loker ES, et al. (2012) Will all scientists working on snails and the diseases they transmit please stand up? PLoS Negl Trop Dis 6: e1835.
  13. 13. Cantacessi C, Prasopdee S, Sotillo J, Mulvenna J, Tesana S, et al. (2013) Coming out of the shell: building the molecular infrastructure for research on parasite-harbouring snails. PLoS Negl Trop Dis 7: e2284.
  14. 14. Lockyer AE, Spinks J, Kane RA, Hoffmann KF, Fitzpatrick JM, et al. (2008) Biomphalaria glabrata transcriptome: cDNA microarray profiling identifies resistant- and susceptible-specific gene expression in haemocytes from snail strains exposed to Schistosoma mansoni. BMC Genomics 9: 634.
  15. 15. Loker ES (2010) Gastropod immunobiology. Adv Exp Med Biol 708: 17–43.
  16. 16. Sermswan R, Mongkolsuk S, Panyim S, Sirisinha S (1991) Isolation and characterization of Opisthorchis viverrini specific DNA probe. Mol Cell Probes 5: 399–407.
  17. 17. Wongratanacheewin S, Pumidonming W, Sermswan RW, Maleewong W (2001) Development of a PCR based method for the detection of Opisthorchis viverrini in experimentally infected hamsters. Parasitology 122: 175–180.
  18. 18. Cantacessi C, Mulvenna J, Young ND, Kasny M, Horak P, et al. (2012) A deep exploration of the transcriptome and the “excretory/secretory” proteome of adult Fascioloides magna. Mol Cell Proteomics 11: 1340–1353.
  19. 19. Bentley DR, Balasubramanian S, Swerdlow HP, Smith GP, Milton J, et al. (2008) Accurate whole human genome sequencing using reversible terminator chemistry. Nature 456: 53–59.
  20. 20. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, et al. (2011) Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 29: 644–52.
  21. 21. Berriman M, Haas BJ, LoVerde PT, Wilson RA, Dillon GP, et al. (2009) The genome of the blood fluke Schistosoma mansoni. Nature 460: 352–358.
  22. 22. The Schistosoma japonicum Genome Sequencing and Functional Analysis Consortium (2009) The Schisosoma japonicum genome reveals features of host-parasite interplay. Nature 460: 345–351.
  23. 23. Young ND, Jex AR, Li B, Liu S, Yang L, et al. (2012) Whole-genome sequence of Schistosoma haematobium. Nat Genet 44: 221–225.
  24. 24. Huang Y, Chen W, Wang X, Liu H, Chen Y, et al. (2013) The carcinogenic liver fluke, Clonorchis sinensis: new assembly, reannotation and analysis of the genome and characterization of tissue transcriptomes. PLoS One 8: e54732.
  25. 25. Cantacessi C, Jex AR, Hall RS, Young ND, Campbell BE, et al. (2010) A practical, bioinformatic workflow system for large data sets generated by next generation sequencing. Nucleic Acids Res 38: e171.
  26. 26. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402.
  27. 27. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, et al. (2000) Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 25: 25–29.
  28. 28. Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, et al. (2008) High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res 36: 3420–3435.
  29. 29. Iseli C, Jongeneel CV, Bucher P (1999) ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc Int Conf Intell Syst Mol Biol 138–148.
  30. 30. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-seq. Nat Methods 5: 621–628.
  31. 31. Audic S, Claverie JM (1997) The significance of digital gene expression profiles. Genome Res 7: 986–995.
  32. 32. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I (2001) Controlling the false discovery rate in behavior genetics research. Behav Brain Res 125: 279–284.
  33. 33. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B 57: 289–300.
  34. 34. Boyle EI, Weng S, Gollub J, Jin H, Botstein D, et al. (2004) GO: TermFinder – open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology term associated with a list of genes. Bioinformatics 20: 3710–3715.
  35. 35. Sithithaworn P, Haswell-Elkins M (2003) Epidemiology of Opisthorchis viverrini. Acta Trop 88: 187–194.
  36. 36. Sayasone C, Odermatt P, Phoumindr N, Vongsaravane X, Sensombath V, et al. (2007) Epidemiology of Opisthorchis viverrini in a rural district of southern Lao PDR. Trans R Soc Trop Med Hyg 101: 40–47.
  37. 37. Wang YC (2012) Examining landscape determinants of Opisthorchis viverrini transmission. Ecohealth 9: 328–341.
  38. 38. Vichasri S, Viyanant V, Upatham ES (1982) Opisthorchis viverrini: intensity and rates of infection in cyprinoid fish from an endemic focus in northeast Thailand. Southeast Asian J Trop Med Public Health 13: 138–141.
  39. 39. Upatham ES, Viyanant V, Kurathong S, Rojborwonwitaya J, Brockelman WY, et al. (1984) Relationship between prevalence and intensity of Opisthorchis viverrini infection and clinical symptoms and signs in a rural community in northeast Thailand. Bull WHO 62: 451–461.
  40. 40. Brockelman WY, Upatham ES, Viyanant V, Ardsungnoen S, Chantanawat R (1986) Field studies on the transmission of the human liver fluke, Opisthorchis viverrini, in northeast Thailand: population changes of the snail intermediate host. Int J Parasitol 16: 545–552.
  41. 41. Prasopdee S, Kulsantiwong J, Piratae S, Khampoosa P, Thammasiri C, et al. (2013) Temperature dependence of Opisthorchis viverrini infection in first intermediate host snail, Bithynia siamensis goniomphalos. Acta Trop In press.
  42. 42. Kiatsopit N, Sithithaworn P, Saijuntha W, Boonmars T, Tesana S, et al. (2012) Exceptionally high prevalence of infection of Bithynia siamensis goniomphalos with Opisthorchis viverrini cercariae in different wetlands in Thailand and Lao PDR. Am J Trop Med Hyg 86: 464–469.
  43. 43. Hanington PC, Lun CM, Adema CM, Loker ES (2010) Time series analysis of the transcriptional responses of Biomphalaria glabrata throughout the course of intramolluscan development of Schistosoma mansoni and Echinostoma paraensei. Int J Parasitol 40: 819–831.
  44. 44. Ittiprasert W, Miller A, Myers J, Nene V, El-Sayed N, et al. (2010) Identification of immediate response genes dominantly expressed in juvenile resistant and susceptible Biomphalaria glabrata snails upon exposure to Schistosoma mansoni. Mol Biochem Parasitol 169: 27–39.
  45. 45. Galinier R, Portela J, Mone Y, Allienne JF, Henri H, et al. (2013) Biomphalysis, a new β pore-forming toxin involved in Biomphalaria glabrata immune defense against Schistosoma mansoni. PLoS Pathog 9: e1003216.
  46. 46. Lockyer AE, Emery AM, Kane RA, Walker AJ, Mayer CD, et al. (2012) Early differential gene expression in haemocytes from resistant and susceptible Biomphalaria glabrata snails in response to Schistosoma mansoni. PLoS One 7: e51102.
  47. 47. Yoshino TP, Wu XJ, Gonzalez LA, Hokke CH (2013) Circulating Biomphalaria glabrata hemocyte subpopulations possess shared schistosome glycans and receptors capable of binding larval glycoconjugates. Exp parasitol 133: 28–36.
  48. 48. Kruatrachue M, Jantataeme S, Ratanatham S, Vichasri S, Upatham ES (1982) A culture method for Bithynia (Prosobranchia: Bithyniidae), snail hosts for the trematode Opisthorchis viverrini. Malacol Rev 15: 63–67.
  49. 49. Hathaway JJ, Adema CM, Stout BA, Mobarak CD, Loker ES (2010) Identification of protein components of egg masses indicates parental investment in immunoprotection of offspring by Biomphalaria glabrata (Gastropoda, Mollusca). Dev Comp Immunol 34: 425–435.
  50. 50. Lockyer AE, Spinks JN, Walker AJ, Kane RA, Noble LR, et al. (2007) Biomphalaria glabrata transcriptome: identification of cell-signalling, transcriptional control and immune-related genes from open reading frame expressed sequence tags. Dev Comp Immunol 31: 763–782.
  51. 51. Baeza Garcia A, Pierce RJ, Gourbal B, Werkmeister E, Colinet D, et al. (2010) Involvement of the cytokine MIF in the snail host immune response to the parasite Schistosoma mansoni. PLoS Pathog 6: e1001115.
  52. 52. Yoshino TP, Lodes MJ, Rege AA, Chappell CL (1993) Proteinase activity in miracidia, transformation excretory-secretory products, and primary sporocysts of Schistosoma mansoni. J Parasitol 79: 23–31.
  53. 53. Myers J, Ittiprasert W, Raghavan N, Miller A, Knight M (2008) Differences in cysteine protease activity in Schistosoma mansoni-resistant and -susceptible Biomphalaria glabrata and characterization of the hepatopancreas cathepsin B full-length cDNA. J Parasitol 94: 659–668.
  54. 54. Bouchut A, Coustau C, Gourbal B, Mitta G (2007) Compatibility in the Biomphalaria glabrata/Echinostoma caproni model: new candidate genes evidenced by suppressive-subtractive hybridization. Parasitology 134: 575–588.
  55. 55. Zhang XL, Ali MA (2008) Ficolins: structure, function and associated diseases. Adv exp Mol Biol 632: 105–115.
  56. 56. Endo Y, Liu Y, Fujita T (2006) Structure and function of ficolins. Current topics in complement. Adv Exp Med Biol 586: 265–279.
  57. 57. Doolittle RF, McNamara K, Lin K (2012) Correlating structure and function during the evolution of fibrinogen-related domains. Protein Sci 21: 1808–1823.
  58. 58. Zhang SM, Zeng Y, Loker ES (2007) Characterization of immune genes from the schistosome host snail Biomphalaria glabrata that encode peptidoglycan recognition proteins and gram-negative bacteria binding protein. Immunogenetics 59: 883–898.
  59. 59. Yoshino TP, Dinguirard N, Kunert J, Hokke CH (2008) Molecular and functional characterization of a tandem-repeat galectin from the freshwater snail Biomphalaria glabrata, intermediate host of the human blood fluke Schistosoma mansoni. Gene 411: 46–58.
  60. 60. Hertel LA, Adema CM, Loker ES (2005) Differential expression of FREP genes in two strains of Biomphalaria glabrata following exposure to the digenetic trematodes Schistosoma mansoni and Echinostoma paraensei. Dev Comp Immunol 29: 295–303.
  61. 61. Zhang SM, Zeng Y, Loker ES (2008) Expression profiling and binding properties of fibrinogen-related proteins (FREPs), plasma proteins from the schistosome snail host Biomphalaria glabrata. Innate Inmmun 14: 175–189.
  62. 62. Adema CM, Hertel LA, Miller RD, Loker ES (1997) A family of fibrinogen-related proteins that precipitates parasite-derived molecules is produced by an invertebrate after infection. Proc Natl Acad Sci USA 94: 8691–8696.
  63. 63. Hasnain SZ, Gallagher AL, Grencis RK, Thornton DJ (2013) A new role for mucins in immunity: insights from gastrointestinal nematode infections. Int J Biochem Cell Biol 45: 364–374.
  64. 64. Bayne CJ, Loker ES, Yui MA, Stephens JA (1984) Immune-recognition of Schistosoma mansoni primary sporocysts may require specific receptors on Biomphalaria glabrata hemocytes. Parasite Immunol 6: 519–528.
  65. 65. Ittiprasert W, Nene R, Miller A, Raghavan N, Lewis F, et al. (2009) Schistosoma mansoni infection of juvenile Biomphalaria glabrata induces a differential stress response between resistant and susceptible snails. Exp Parasitol 123: 203–211.
  66. 66. Zahoor Z, Davies AJ, Kirk RS, Rollinson D, Walker AJ (2010) Larval excretory-secretory products from the parasite Schistosoma mansoni modulate HSP70 protein expression in defence cells of its snail host, Biomphalaria glabrata. Cell Stress Chaperones 15: 639–650.
  67. 67. Bayne CJ (2009) Successful parasitism of vector snail Biomphalaria glabrata by the human blood fluke (trematode) Schistosoma mansoni: a 2009 assessment. Mol Biochem Parasitol 165: 8–18.
  68. 68. Blouin MS, Bonner KM, Cooper B, Amarasinghe V, O'Donnell RP, et al. (2013) Three genes involved in the oxidative burst are closely linked in the genome of the snail, Biomphalaria glabrata. Int J Parasitol 43: 51–55.
  69. 69. Laha T, Pinlaor P, Mulvenna J, Sripa B, Sripa M, et al. (2007) Gene discovery for the carcinogenic human liver fluke, Opisthorchis viverrini. BMC Genomics 8: 189.