Thanasuwan et al. Parasites & Vectors 2014, 7:224 http://www.parasitesandvectors.com/content/7/1/224

RESEARCH

Open Access

Suppression of aquaporin, a mediator of water channel control in the carcinogenic liver fluke, Opisthorchis viverrini Sirikanda Thanasuwan1, Supawadee Piratae1, Paul J Brindley2, Alex Loukas3, Sasithorn Kaewkes1 and Thewarach Laha1,4*

Abstract Background: Opisthorchiasis and Opisthorchis viverrini-associated bile duct cancer represent major public health threats in Thailand and Laos. The tegument of this food borne fluke plays pivotal roles in parasite metabolism, homeostasis and osmoregulation. Excretory/secretory products also pass from the fluke to the biliary environment, products that likely underlie pathogenesis of liver fluke infection. Aquaporins (AQPs), belong to the major intrinsic protein superfamily of integral plasma membrane channel proteins that selectively transport water across cell membranes. AQPs play key roles as water and ion transport channels through the tegument of helminth parasites. Methods: Here, two forms of AQP mRNAs from the adult developmental stage of O. viverrini, termed O. viverrini aquaporin-1 and −2 (Ov-aqp-1 and −2) were investigated. Roles of Ov-aqp-1 and −2 in the movement of water across the tegument of this carcinogenic liver fluke were investigated using RNA interference. Results: Ov-AQP-1 and Ov-AQP-2 contain unique characteristic asparagine-proline-alanine (NPA) motifs of AQP transmembrane proteins. Phylogenetic analysis indicated that Ov-AQPs belong to an expanding group of aquaglyceroporin-like water channel proteins characterized from helminth and protozoan parasites, which is pivotal to the specialized requirements of water and solute control during parasitism. Elevated transcription of Ov-aqp-1 was evident in the egg, cercaria, metacercaria and adult stages of O. viverrini, whereas Ov-aqp-2 transcripts were detected at higher level in egg, metacercaria, cercaria and adult stage, respectively. RNA interference using electroporated dsRNA suppressed transcript levels of Ov-aqp-1 and Ov-aqp-2 in adult worms by 58-99% over periods of up to 16 days in vitro. Suppression of Ov-aqp-1 and Ov-aqp-2 in vitro disabled water transport in adult flukes. Conclusion: The apparently pivotal roles of Ov-AQP in solute homeostasis at the fluke surface suggest that deeper investigation will be informative for the pathophysiology of O. viverrini, and may uncover intervention targets, particularly in view of the singularly notable predilection of this pathogen for residence within ducts of the biliary tree. Keywords: O. viverrini, Cholangiocarcinoma, RNA interference, Aquaporin, Liver fluke, Tegument

* Correspondence: [email protected] 1 Department of Parasitology, Faculty of Medicine, Khon Kaen University, Khon Kaen 40002, Thailand 4 Liver Fluke and Cholangiocarcinoma Research Center, Faculty of Medicine, Khon Kaen University, Khon Kaen 40002, Thailand Full list of author information is available at the end of the article © 2014 Thanasuwan et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Thanasuwan et al. Parasites & Vectors 2014, 7:224 http://www.parasitesandvectors.com/content/7/1/224

Background Opisthorchis viverrini is highly endemic in Southeast Asia, where more than 10 million people are infected [1,2]. Studies with hamster models of both chronic O. viverrini infection and exposure to nitrosamines have determined that host responses are key factors in carcinogenesis of liver fluke-induced cholangiocarcinoma [3]. Membrane proteins of O. viverrini without orthologues/paralogues in mammals are of particular interest for the development of vaccines and drugs, because the absence of host homologues raises therapeutic margins of safety [4-6]. AQPs are the major intrinsic protein (MIP) of integral plasma membrane channel proteins that are passively permeated by water and small, uncharged solutes [7,8]. AQPs have been identified based on their highly conserved dual asparagine-proline-alanine (NPA) boxes, critical for the formation of a water-permeating pore. The sequences of conventional AQPs have only minimal or modest identity to one another but they share conserved six transmembrane domains and hydrophobic NPA box-like repeats [9,10]. Moreover, some AQP-like sequences exhibit only poor sequence conservation for the NPA motifs [11]. The NPA motifs of AQPs play crucial roles for movement of water across cell membranes [8]. AQPs have been investigated in several parasites and observed to provide key functions in transport of water and other small solutes. Moreover, AQPs facilitate and inhibit uptake of lactate and other anthelmintics [12-14]. In schistosomes, AQP is a major tegument protein with functions in trans-tegumental water movement, absorption of nutrients and other metabolites [14,15]. AQPs are highly expressed in the transcriptomes and tegumental proteome of O. viverrini [16,17]. In this study, discrete mRNA sequences encoding two AQP-like transporters were isolated from a cDNA library established from the adult stage of O. viverrini. The properties of these water channel transporters were investigated using bioinformatics approaches, phylogenetic analyses and gene silencing via RNA interference to highlight the physiological roles of AQPs in liver fluke biology and parasitism. Methods Opisthorchis viverrini

Metacercariae of O. viverrini were collected from the flesh of naturally infected cyprinid fish from Khon Kaen province, Thailand by digestion with 0.25% pepsin, as described [18]. Syrian golden hamsters, Mesocricetus auratus, purchased from the Animal Unit, Faculty of Medicine, Khon Kaen University were infected with 50 metacercariae by orogastric gavage [19]. Hamsters were maintained at the animal facility, Faculty of Medicine, Khon Kaen University using protocols approved for

Page 2 of 11

animal experimentation by the Animal Ethics Committee of Khon Kaen University, based on the Ethics of Animal Experimentation of the National Research Council of Thailand (Approval number AEKKU43/2555). Hamsters were euthanized six weeks after infection, when adult flukes were recovered from bile ducts and the gall bladder. The worms were rinsed in sterile 0.9% NaCl to remove residual host cells and debris. Eggs of O. viverrini were collected from worms recovered from euthanized hamsters; flukes were maintained in RPMI media containing antibiotics (streptomycin/penicillin, 100 μg/ml) at 37°C, in an atmosphere of 5% CO2 and incubation for 18 h [20]. Eggs were collected by centrifugation at 5,241 g for 10 min and stored at −70°C. Cercariae of O. viverrini were shed from naturally infected Bithynia sp. snails collected in farmlands in Khon Kaen province [21]. Snails were placed into plastic containers filled with de-chlorinated water, 4–5 snails per container, and exposed to the light for 2 h, after which cercariae were collected by centrifugation of the supernatant water at 5,241 g, 10 min. Isolation of aquaporin genes from O. viverrini

cDNA sequences encoding full open reading frames (ORFs) of Ov-aqp-1 and Ov-aqp-2 were amplified by PCR from a cDNA library of transcripts from the adult developmental stage of the fluke [16]. The specific primers for PCR amplification of the Ov-aqp genes were designed based on expressed sequence tags (ESTs) and EST contigs. cDNA sequences encoding full length ORFs of Ov-aqp-1 (GenBank accession EL618688) [16] and Ov-aqp-2 (OV_ contiq1681, available at http://bioinfosecond.vet.unimelb. edu.au/) [22] were identified from our previous transcriptomics study. The primers for Ov-AQP-1 F were 5′-AG CATGGCTGGTAGTCTCTCATC and Ov-AQP-1R5′-AG CGGATCCTCAGTTTTTTTTCTGGCG. The primers of Ov-AQP-2 F were 5′-AGCCATATGATGAGTTTGGAA TGCGAAAC and Ov-AQP-2R5′-AGCGGATCCCTAGG CAAGCAGTTCAGTTC. PCR reaction mixes included the O. viverrini cDNA library (100 ng), 0.2 mM dNTP, 1.5 mM MgCl2 with 1 unit Taq polymerase (Invitrogen, USA). Amplification was accomplished with 35 cycles of denaturation at 95°C for 1 min, annealing at 60°C for 1 min, extension at 72°C for 2 min and a final extension at 72°C for 10 min. The amplicons were separated and sized by electrophoresis through agarose followed by staining with ethidium bromide. Products of interest were isolated from the gel using a kit (GeneJET™ gel extraction, Fermentas, EU), ligated into the pGEM-T Easy vector (Promega, USA), after which ligation products were used to transform Escherichia coli strain JM109 competent cells (Promega). Plasmids isolated from the resulting colonies were sequenced using BigDye terminator method (1st BASE, Singapore); sequences were analyzed using Blast search

Thanasuwan et al. Parasites & Vectors 2014, 7:224 http://www.parasitesandvectors.com/content/7/1/224

against GenBank databases [23] and compared to consensus sequences of Ov-AQPs from ESTs of O. viverrini [16,22]. Recombinant plasmids containing cDNA sequences encoding the entire ORFs of Ov-aqp-1 and Ovaqp-2 were termed pOv-aqp-1 and pOv-aqp-2. Sequence and phylogenetic analyses

Nucleotide sequences and chromatograms were evaluated using BioEdit V7.0.5 [24]. The edited sequences were translated to ORFs with assistance of software at http:// bio.lundberg.gu.se/edu/translat.html. Signal peptides from the deduced amino acid sequences were predicted and analyzed by SignalP 3.0 Server at http://www.cbs.dtu.dk/ services/SignalP/. Transmembrane helices were predicted by using TMHMM Server v. 2.0 (http://www.cbs.dtu.dk/ services/TMHMM-2.0/). N-Glycosylation sites were predicted by using the NetNGlyc 1.0 Server at http://www. cbs.dtu.dk/services/NetNGlyc/. Other divergent sequences were compared and multiple alignments constructed by ClustalW in the BioEdit program [24]. ORFs of Ov-AQPs were aligned with well characterized AQP protein sequences from informative species using ClustalW in the BioEdit program [24]. A phylogenetic tree was constructed with p-distance matrix using the neighbor-joining method [25] with 1,000 bootstrap samplings in the MEGA version 6.06 [26]. RNA extraction and quantitative real time reverse transcription PCR

Real time, reverse transcription PCR was undertaken to monitor the expression of the Ov-aqp genes during the developmental cycle of the liver fluke - egg, cercaria, metacercaria and adult stages were examined, and at intervals after exposure by electroporation to dsRNA. Briefly, total RNA was extracted from the developmental stages of the parasite using TRIZOL (Invitrogen). Any residual DNA remaining in the RNA preparations was removed by DNase digestion. Double stranded cDNA was synthesized from equal amounts of total RNA template (1 μg) using a cDNA synthesis kit (Fermentas). Quantitative real-time PCR was performed using custom SYBR Green Assays. The primers to detect Ov-aqp-1 (spanning coding DNA positions 3–280), were AQP1_EXF, 5′-GG CTGGTAGTCTCTCATC-3′ and AQP1_EXR, 5′-CGTA TCCCATAGTACCGCTG-3′. Ov-aqp-2 transcripts (spanning nt positions 16 to 256) were designated Ov-AQP2_ EXF: 5′-GAAACCCGATTTCGAAGAGG and Ov-AQ P2_EXR: 5′-TGATCCCGGAGAAGAATACG. PCRs were performed in triplicate using SYBR Green reagents and a thermal cycler with a real time detector (ABI 7500); SYBR Green reactions were prepared by adding 12.5 μl of SYBR Green Master Mix (TAKARA Perfect Real-time Kit, Japan), 0.5 μl (10 mM) of forward primer and reverse primers, 0.5 μl of reference dye (ROX), 1 μl (equivalent to

Page 3 of 11

50 ng of total RNA) of first-stand cDNA and water to a final volume of 25 μl. The thermal cycling conditions used were: initiation pre-heat for one cycle at 95°C, 10 min; 40 cycles of denaturation at 95°C, 30 sec; annealing at 55°C, 30 sec; extension at 72°C, 45 sec. Expression levels of the Ov-aqp-1, Ov-aqp-2and actin mRNAs (OvAE1657, GenBank EL620339) were determined as described [20]. To determine the extent of gene silencing induced by dsRNAs, the mRNA expression levels of Ov-aqp-1 and Ov-aqp-2 (or firefly luciferase as an irrelevant control, below) were normalized with actin mRNA and presented as the unit value of 2-ΔΔCt where ΔΔCt = ΔCt (treated worms) - ΔCt (non-treated worms) [20,27]. Data are presented as the mean ± 1 standard deviation. Differences between groups were assessed using Student’s t-test (GraphPad Prism Software); p values of ≤ 0.05 were considered statistically significant. Preparation and delivery of dsRNA

dsRNAs were designed to span 605 nt of Ov-aqp-1 (full ORF spanning nt193-798) and 575 nt of Ov-aqp-2 (full ORF spanning nt 201–776). Target sequences were amplified from plasmids pOv-aqp-1 and pOv-aqp-2 (above) using primers flanked with a T7 RNA polymerase promoter sequence, indicated in underlined italic bold faced, at the 5′ end. Ov-aqp1 was generated using primers ds-aqp1_T7-F, 5′-TAATACGACTCACTATAGGGGGT AGCAACGTCTCGGCTand ds-aqp1_T7-R, 5′-TAATAC GACTCACTATAGGGGTAGAGTAACACTCCGAG. Ovaqp-2 was generated using primers ds-aqp2_T7-F, 5′TAATACGACTCACTATAGGGAGTGAGTCTTGGCTG GGGTA and ds-aqp2_T7-R, 5′-TAATACGACTCACTA TAGGGTACGGTCCGACGATTGGTAT. dsRNAs were synthesized using a MEGAscript RNAi Kit (Ambion, USA). The irrelevant negative control, luciferase (LUC) dsRNA was constructed from plasmid pGL3-basic (Promega), because this sequence does not match any targets in the O. viverrini genome. It was amplified using primers ds-LUC_T7-F5′-TAATACGACTCACTATAGGGTGCGC CCGCGAACGACATTTA and ds-LUC_T7-R5′-TAATA CGACTCACTATAGGGGCAACCGCTTCCCCGACTTC CTTA [28]. The thermal cycling conditions were 35 cycles of denaturation at 94°C, 30 sec; annealing at 55°C, 30 sec; and extension at 72°C, 60 sec. Amplicons were sized by electrophoresis and purified, as above. Concentrations of dsRNA were determined by spectrophotometer (NanoVue, GE Healthcare, USA). To deliver Ov-aqp dsRNAs, adult worms (30 per treatment group) were washed with sterile phosphate buffer, then transferred to a cuvette, 4 mm gap (Bio-Rad) in 100 μl electroporation buffer (RPMI-1640, 1x antibiotic/antimycotic, 1% glucose, 1 mM trans-epoxysuccinyl-L-leucylamido (4-guanidino) butane (E64) (Sigma) containing 50 μg Ov-aqp-1, Ov-aqp-2 or luc dsRNAs. Each group of

Thanasuwan et al. Parasites & Vectors 2014, 7:224 http://www.parasitesandvectors.com/content/7/1/224

worms was subjected to square wave electroporation using a single 20 ms impulse at 125 V (Gene PulserXcell, Bio-Rad, USA). Subsequently, worms were maintained in RPMI culture medium (one ml culture medium in each well of a 24 well plate) supplemented with 2 mg/ml of dsRNA at 37°C under 5% CO2 in air. Similarly treated worms in electroporation buffer without dsRNA served as controls. Worms were soaked in 2 mg dsRNA for 16 days (Ov-aqp-1) and 7 days (Ov-aqp-2) with changes of media containing dsRNA every second day. Worms were collected for analysis on days 1, 3, 6, 10 and 16 for Ov-aqp-1 and on days 1, 3 and 7 for Ov-aqp-2 following electroporation. Water transport assay and size measurements of parasites

Water transportation in normal worms was investigated to obtain determine optimal times for observation of the transformed worms. The volume of water transport via the body wall of adult flukes was investigated by monitoring swelling of the flukes after exposure to water. Wild type adults were incubated in distilled water for 5, 10, 30, 60 and 120 min. Photomicrographs were captured on these cultures, after which surface area of the flukes was measured using the NIS-Elements D 3.22.00 software (Build 710) – Driver selection (Nikon, Japan). The function of the O. viverrini AQPs as water channels was investigated using RNAi to knock down the genes after treatment with Ov-aqp-1, Ov-aqp-2 and Ov-aqp-1 + Ov-aqp-2 dsRNAs. Water transportation through the tegument of control and dsRNA-treated worms was determined by measuring the surface area of the fluke after placing the worms in water for 10 min. Briefly, worms electroporated with dsRNAs of Ov-aqp or luc were cultured in complete RPMI medium for 24 h prior to transfer to hypo-osmotic medium (distilled water) for 10 min. Photomicrographs were recorded and the surface area of individual worms ascertained before and after exposure to the hypo-osmotic conditions. To compare surface area of the liver flukes before and after transfer, images were recorded so that the surface area of individual adult worms could be determined. The liver flukes were photographed before and 10 min after transfer into water, and the sizes of a sample of individual parasites from control and Ov-aqp suppressed groups were ascertained and compared. Histological examination

To observe swelling of the worms, sections of fixed flukes after dsRNA treatment for 1, 3 and 7 days were examined. Briefly, the whole worms were fixed in hot 10% buffered formalin overnight then processed conventionally. Adult worms in paraffin sections were de-paraffinized in xylene three times, five minutes each. The 4 μm thick sections were stained with hematoxylin and eosin. Thereafter, the

Page 4 of 11

sections were dehydrated, cleared and mounted. Histology of adult worms was compared and contrasted using light microscopy [29]. Statistical analysis

Means ± standard deviations were calculated and presented using Microsoft Excel or GraphPad Prism software. Differences between groups were assessed using Student’s t-test (GraphPad Prism Software, USA); p ≤ 0.05 was considered statistically significant.

Results Sequences and structure of O. viverrini aquaporin

The full length cDNA sequences encoding two isoforms of O. viverrini AQPs (Ov-AQPs) were designated Ov-aqp1and Ov-aqp-2; nucleotide sequences of Ov-aqp-1 and Ov-aqp-2 have been assigned GenBank accession numbers KF697690 and KF697691, respectively. Ov-aqp-1and Ov-aqp-2 are 864 and 921 nt and encode 287 and 306 amino acid ORFs, respectively. Both Ov-AQP-1 and OvAQP-2 are predicted to contain six- transmembrane domains (TM1-6) connected by five loops (A-E) (Figure 1). Potential N-glycosylation sites are present at Asn-67 in Ov-AQP-1 and Asn-293 in Ov-AQP-2 (Figure 1A). OvAQP-1 and Ov-AQP-2 sequences contain conserved NPA motifs as found in other AQPs. Both the amino and carboxyl termini of Ov-AQP-1 and Ov-AQP-2 are predicted to be intracellular Ov-AQP-1 consists of two hemipores, each hemipore containing three transmembrane domains (Figure 1B). The pore that forms in the AQPs is composed of two hemipores where loops B and E form the poreforming signature motif, NPA (Asn-Pro-Ala) via an overlap to form the characteristic pore between the membrane bilayer (Figure 1B). It is noteworthy that both of these two new water channel molecules include the aspartic acid residue (D) in the second NPA box that is a diagnostic sequence for aquaglyceroporins, a category of aquaporinlike water channel proteins discrete from the ‘classical’ aquaporins. Aquaporins sensu stricto are water selective or specific water channels, whereas the aquaglyceroporins are permeable to other small, uncharged molecules, in particular glycerol, in addition to water (Benga 2012). Phylogenetic analysis of Ov-AQP

The phylogenetic relationships of Ov-AQP-1 and OvAQP-2 protein sequences were investigated by comparison with well-characterized AQPs from various animals [9]. A neighbor-joining tree indicated that Ov-AQPs formed a clade with AQPs from several invertebrates, including other trematodes, and also nematode and protozoan parasites such asaquaporin-3 from Clonorchis sinensis (GAA31414) [30] (Figure 2). Ov-AQP-2 formed a robust clade with aquaporin-9 from C. sinensis (GAA55320 and GAA55323) [30] and a group of aquaporins from the

Thanasuwan et al. Parasites & Vectors 2014, 7:224 http://www.parasitesandvectors.com/content/7/1/224

Page 5 of 11

A B1

TM1 TM2 N B2 Ov-AQP-1 Ov_AQP1 : M AG S L S S Y E Y ER V SC Q S K Y E S T I H R WAD KMR L T R P L L L R V F I G E L MG T M I L V I F GD A V V AQ S V F G - - - - - SN V S A L S V S F GWG L A V I F G I 85 Cs-AQP-3 : M PG S L S S Y E Y ER V SC Q S K Y EN T I H R WAD K I R L T R L P L L R V F I G E L AG T M I L I I F GN AM V AQ S V F G - - - - - SN V S A F S I S F GWG L A V A F G V 85 Cs_AQP3 B3 B4

TM1

Ov-AQP-2 Ov_AQP2 : Cs-AQP-9 : Cs-AQP9 Sm-AQP : Sm_AQP Human-AQP Human_AQP: Fg-AQP-1 : Fg_AQP1

-

-

-

-

-

-

-

- - - M S L EC E T R F R R G L R T L AD K I R L R SC P L I R AC I - - - M S L EC E T R F R KG L R T L AD K I R L R SC P L I R AC I - MGR T D T F EGR L ER R QH S L AD K I R L N R C P F I R AC L - - - - - - - - - MGR Q K E L V SR C G EM L H I R YR L L R Q A L - - - - - - - - M SQ Y YD D VN Y K A T T A T L PR C D WY - - H I

TM2

A E F C G T A I L I I F GC G V L AQ V F L GD H G K K SH GG F L A E F C G T A V L I I F GC G V L AQ V F L GD H G K K SH GG F L A E F C G T A I L M I F G SG V L AQ T I L GD H G KQ AH GG F L A EC L G T L I L VM F GC G S V AQ V V L SR G - - - T H GG F L K F MC R V F A A E A L G L GM I V F I V AC H G A EGN L A SN I

S V S L GWGM A V T F G V S V S L GWGM A V T F G V S V S L GWG F A V YMG V T I N L A FG F A V T L G I SG P V T A S A A F A V A V

80 80 82 71 73

B1

TM3 TM2 TM4 NPA B2 Ov-AQP-1 Ov_AQP1 : F T SG T MG YG L L N P A V A L A F A F N G K V PW I A V P V A T I A E L I G A F L GG L I V YC N YQ VN I Q L - - - - L D PG - - L T T K T L G I F C T A P S - N AGN L AG 168 Cs-AQP-3 : F I SG T L G YG L L N P A V A L A F A F D A K V PW I V V P V A T I A E V I G A F L GG L I V YC N YQ VN I Q L - - - - L D PG - - L T T K T L G I F C T A P S - N AGN L AG 168 Cs_AQP3 B3

TM3 TM2 NPA B4 Ov-AQP-2 Ov_AQP2 : F F SG I T G SG N I N P A I T L A F A L I G K V P L KR I P F Y T I SQ V L G A Y VG A A T V Y V I YR E K I I ER VN T VD EG L Y L VD T T G N I F V T N P - - WA SH T T C Cs-AQP-9 : F F SG I AG SG N I N P A I T L A F A L I G K V P L KR I P F Y T V SQ V L G A Y VG A S T V Y L I YR E K I L ER VN T VD EG S Y L VD T T G N I F V T N P - - WA SH T T C Cs-AQP9 Sm-AQP : I F SG C GG SG H I N P A V T L A F S I I G K F P L R R V P F Y T I SQ I L G A F VG S L A I YG V Y Y E K I H E Y A L T R D N G S L I V K T T G G I F V T N P - - G A SH L T C Sm_AQP Human-AQP Human_AQP: L I AG Q V SG AH L N P A V T F AMC F L AR E PW I K L P I Y T L AQ T L G A F L G AG I V F G L Y YD A I WH F AD N Q L F A S - G PN G T AG I F A T Y P SGH L D M I N G Fg-AQP-1 : WT VG P I N G PQ I T A A L S V A L L M T R R I N F V Y F I L G I AG Q L C G A L VG I S L A SR L V PG - - - - - - - - - - - - - L L D K AN L N L H I PG PG - - V T D GQ A Fg_AQP1

168 168 170 160 148

B1

TM4 TM6 TM5 NPA B2 Ov-AQP-1 Ov_AQP1 : F YD Q I I G T A L L T L M V L AMG D C - - F N V P T Y - - - - V F P I L VG L L I T A L VG A F G VN AG A A L N P AR D L G PR L A A F A VG Y - KD A F T A YD S YC WV P 251 Cs-AQP-3 : F YD Q I I G T A L L A L M V L AMG D Y - - F D V P T Y - - - - V F P I F VG L L I T A L VG A F G VN AG A A L N P AR D L G PR L A A L A VG Y - KD A F T A YN S YC WV P 251 Cs_AQP3 B3 B4

TM4

Ov-AQP-2 Ov_AQP2 : Cs-AQP-9 : Cs-AQP9 Sm-AQP : Sm_AQP Human-AQP Human_AQP: Fg-AQP-1 Fg_AQP1 :

F I DQ V I G T A L L F I DQ V I G T A L L F LDQ V LG T A L L F FDQ F I G T A S L F GM EC I C S F L L

TM5 S AG A L S AG A L A AG A L I VC V L I VC C L

TM6

NPA

A V VD Y KGWKM PR Y - - - - L H P L Y L G L L V Y T L VGC F AMN AGC A L N P AR D L G PR L M L L MC G WG T S A F R D KD Y F F W I P A V VD Y KGWKM PG Y - - - - L H P L Y L G L L V Y T L VGC F AMN AGC A L N P AR D L G PR L M L L MC G WG T S A F R D KN Y F F W I P A I T D R N GWK L PD Y - - - - L H P L H L A F L V Y A L VGC F A L N AG A A L N P AR D L G PR L M I L MC G WG G S A F T SGN Y F F WV P A I VD P YN N P V PR G - - - - L E A F T VG L V V L V I G T SMG F N SG Y A VN P AR D F G PR L F T A L AG WG S A V F T T GQH WWWV P S T VD E F R N PH WAQGH V T V F S V I V Y F L I L M L A A I L V K V T GC GMN P A A S L S - - - - - - - - - - - A A I F N N N Y A K I W I Y

254 254 256 246 227

B1

TM6 B2 Ov-AQP-1 Ov_AQP1 : I VG P Y I G A I L G V L L Y E L I I G VH VR G - - - - - - - L K E E E ER Q K KN - - - - - - - - - - - - - - - - - - - - - - - - - - - - Cs-AQP-3 : V AG P Y I G A I VG V L L Y E L L I G VH VR G - - - - - - - L K E E E EH Q K K L KD SD D Y - - - - - - - - - - - - - - - - - - - - - - Cs_AQP3 B3 B4 Ov-AQP-2 Ov_AQP2 : Cs-AQP-9 : Cs-AQP9 Sm-AQP : Sm_AQP Human-AQP Human_AQP: Fg-AQP-1 : Fg_AQP1

TM6 I I I I I

287 293

N

VG P Y VG A V I G A V L Y E L T I G I H L D T - - - AG ER L A I E T PG K AN D T K KD - - - E SG T E L L A - - - - - VG P Y VG A V I G A V L Y E L T I G I H L D T - - - AG ER L A I E T PG K SN D T K KD - - - E SG T E L L A - - - - - L G P Y I G A V I G A I L Y E L T I G I H L D R K S I S S S P AQ T T T P SG EQR F D L D S PG K SG SH L L T R G L - - V S P L L G S I AG V F V YQ L M I GC H L EQ P P P SN E E EN V K L AH V KH K EQ I - - - - - - - - - - - - - - - - - V A P V VG S I I A V L F WEMM L SD G A S A ER I KH WWT D PH F D R T KD Y KR M ER E S KG V YD T E E A I A F S S

B

- - - - - - - - - - - - - - - - - - - - - - - - - - - - I R S K ED L R

306 306 317 292 299

C

A B

EXTRACELLULAR

1

2

N P A

3

4

5

N P A

6

E

CYTOSOL

Hemipore 2

Hemipore 1 NH2

D HOOC

Figure 1 Deduced amino acid sequences, multiple sequence alignment, and predicted secondary structure of aquaporins from Opisthorchis viverrini. Multiple sequence alignment of aquaporins from Opisthorchis viverrini and other eukaryotic organisms (panel A). Homologous protein sequences are derived from GenBank: Ov-AQP-1 (KF697690), Ov-AQP-2 (KF697691), Cs-AQP-3 and Cs-AQP-9 from Clonorchis sinensis (GAA31414 and GAA55320), Sm-AQP from Schistosoma mansoni (XP_002572046), humanAQP from Homo sapiens (CAG46822) and Fg-AQP-1 from Fasciola gigantica (AD032834). Schematic representation of the predicted architecture of Ov-AQPs in the cell membrane (panel B). Ov-AQPs consist of 2 hemipores, each hemipore containing 3 transmembrane domains. Hemipore-1 is composed of transmembrane domains 1, 2, 3 and hemipore-2 is composed of transmembrane domains 4, 5, and 6. The pore-forming signature motif NPA (Asn-Pro-Ala) is located in both loops B and E.

schistosomes (blood flukes) [12]. Ov-AQP-1 was less closely related to AQPs from other trematodes, including Ov-AQP-2, branching at amore distant position in the phylogram along with a predicted AQP from C. sinensis, GAA31414). However, both new Ov-AQP water channel proteins branched with the aquaglyceroporins rather than the ‘classical’ (strict) aquaporins.

Developmental expression of Ov-aqp genes

Levels of expression of Ov-aqp-1 and -aqp-2 mRNAs relative to expression of the actin control gene were examined by qRT-PCR among developmental stages of O. viverrini. Expression of Ov-aqp-1 expression was seen in all developmental stages examined but was highest in metacercariae (Figure 3). Ov-aqp-2 expression was also

Thanasuwan et al. Parasites & Vectors 2014, 7:224 http://www.parasitesandvectors.com/content/7/1/224

Page 6 of 11

100

Sm-AQP, Schistosoma mansoni XP_002572048 Schistosoma japonicum AAW24850

83 82 68

S. mansoni XP_002572046 S. japonicum CAX 76488

91

100

OV-AQP-2 , Opisthorchis viverrini Clonorchis sinensis GAA55320

100 100

61

Aquaglyceroporins

Clonorchis sinensis GAA55323 OV-AQP-1, O. viverrini Clonorchis sinensis GAA31414

95

Ascaris suum ADY46722 Toxocara canis AAC32826 Pf- AQP, Plasmodium falciparum XP_001348009

60

AQP10 ( Danio rerio )

89

AQP9 ( Homo sapiens) AQP9 (Gallus gallus)

58

AQP9 ( Tetraodon nigroviridis ) CAF99586 AQP10 ( Homo sapiens)

99 100

AQP3 (Homo sapiens)CAG46822 AQP3 ( Mus musculus) AAD20305 AQP3 ( Xenopus laevis ) AAH70543 AQP7 (Homo sapiens ) BAC05693

100 95

AQP7 ( Rattus norvegicus )AAN52930 AQP7 ( Mus musculus)AAH22223

Trypanosoma brucei XP_827897

100

Toxoplasma gondii EEB03724 Toxoplasma gondii CAE46485

62 99

Leishmania donovani XP_003860873 Trypanosoma cruzi AAM76680

100 100

80

Fasciola gigantica ADO32834 Fasciola gigantica ADO32835

77

S. mansoni XP_002572574 C. sinensis GAA33659

55

AQP1 ( Mus musculus) AAH07125 AQP1 ( Rattus norvegicus ) NP_036910 AQP1 (Homo sapiens) P29972

95

99

AQP6 (H. sapiens) NP_001643 AQP6(R. norvegicus ) NP_071517

77 100

AQP2 (H. sapiens ) NP_000477 AQP2(R. norvegicus ) P34080

68 100

AQP5 (H. sapiens ) NP_001642

Strict Aquaporins

100

AQP5(R. norvegicus ) NP_036911

Figure 2 Phylogram constructed using the neighbor joining method to compare the relationship between the open reading frames of Ov-AQPs of Opisthorchis viverrini and homologues within the aquaporin superfamily. The accession numbers for sequences included in the phylogenetic analysis are shown in each operational taxonomy unit. The numbers above the branches refer to bootstrap values (where greater than 50% support was obtained). Note the separate branches for strict aquaporins and aquaglyceroporins. The novel O. viverrini aquaporins belong to the latter assemblage, indicating that they can transport additional permants beyond water.

detected in all stages; however, highest expression was evident in the egg stage (Figure 3). Suppression of Ov-aqp gene expression

Ov-aqp-1 and Ov-aqp-2 gene expression was suppressed in O. viverrini in vitro by introducing target-specific

dsRNAs using electroporation followed by soaking to mediate knockdown. Figure 4A shows the specific suppression of Ov-aqp-1 by 99%, 57%, 72%, 83% and 91% measured on days 1, 3 6, 10 and 16 after treatment with the dsRNAs, respectively. Figure 4B shows the specific suppression of Ov-aqp-2 by 58%, 92% and 77% measured

Thanasuwan et al. Parasites & Vectors 2014, 7:224 http://www.parasitesandvectors.com/content/7/1/224

Page 7 of 11

A

Ov-aqp-1

Relative expression level

Ov-aqp-1

Egg

Cercaria

Metacercaria

Adult Day1

Day3

Day6

Day10

Ov-aqp-2

Ov-aqp-2

Egg

Cercaria

Metacercaria

Adult

Relative expression level

B

Figure 3 Levels of expression of mRNAs of Ov-aqp-1 and Ov-aqp2 among developmental stages of the liver fluke, Opisthorchis viverrini. Expression levels relative to Ov-actin were analyzed by real-time qRT-PCR. Two biological samples plus technical duplicates were used to assess expression levels. Data are presented as mean ±1SD.

Day16

dsRNA treatment

Day1

Day3

Day7

dsRNA treatment

Reduced Ov-AQP water transport in Ov-aqp dsRNA-treated flukes

Figure 4 Suppression of aquaporin Ov-aqp mRNA in adult O. viverrini by RNA interference (RNAi). Relative Ov-aqp-1 and Ov-aqp-2 expression in adult Opisthorchis viverrini flukes following electroporation with dsRNAs corresponding to Ov-aqp-1 (panel A) and Ov-aqp-2 (panel B). Transcript levels were determined using SYBR-green real time qRT-PCR. Ov-actin was used as the control gene for normalization. The control group (no dsRNA treatment) was compared to the negative control group that received luc dsRNA. Three biological samples with technical duplicates were used for analysis. Means ±1SD are shown.

At the outset, to assess the capacity of adult O.viverrini to swell following prolonged incubation in distilled water, parasites were measured at intervals after transfer from physiological conditions in vitro (i.e. in RPMI) to hypotonic conditions (distilled water). This led to marked swelling of the worms by 10 min, with the increase stabilizing by 30 min of incubation (Figure 5). The average size of worms at time 0 min was 10.5 mm2; after 5 min of incubation in water the average size of worms was 10.2 mm2 (p ≤ 0.001) (Figure 5A). However, after 10 min of incubation in water, average worm size increased to 14.6 mm2 (p ≤ 0.001) and 20.45 mm2 (p ≤ 0.001) by 30 min of incubation. After 30 min in water, worms did not continue to swell further for the duration of the experiment (120 min). Based on these findings, further studies were carried out using the 10 min time point.

To determine whether silencing of Ov-aqp-1, Ov-aqp-2 and Ov-aqp-1 + 2 expression affected water movement in O. viverrini, adult flukes treated with Ov-aqp dsRNAs were transferred from RPMI into hypotonic medium as described above (distilled water). Parasites were photographed before and 10 min after transfer into water, and the sizes of a sample of individual parasites from control (mock control or luc) and Ov-aqp-suppressed groups were compared. When flukes were not subjected to osmotic shock by culture in RPMI medium, no significant difference in body size was detected between any of the control or dsRNA-treated flukes (range = 9.2 – 11.8 mm2; Figure 6A). When flukes were transferred into water, control flukes swelled to roughly twice their original size (in RPMI) with average sizes of 20.2 mm2 (mock control)

on days 1, 3 and 7 after treatment with the dsRNAs, respectively, compared to the negative control group that received luc dsRNA.

Thanasuwan et al. Parasites & Vectors 2014, 7:224 http://www.parasitesandvectors.com/content/7/1/224

Worm size, mm2

30

***

A

*** *** ***

30

Worm size, mm2

A

Page 8 of 11

20

10

0

5

10

30

60

20

10 5

120

Time (min)

Control

luc

Ov-aqp-1 Ov-aqp-2 Ov-aqp-1+2

dsRNA treatment

B

B

5 min

1.5 mm

10 min

1.5 mm

Figure 5 Swelling of adult Opisthorchis viverrini when incubated in distilled water. Each point represents an individual parasite; the horizontal line indicates the mean for that group; n = 18 per treatment. p ≤ 0.001 denoted with *** (panel A). Morphology of a typical set of worms for each group is shown (panel B).

Worm size, mm2

30

20

***

***

***

10

Control

luc

Ov-aqp-1 Ov-aqp-2 Ov-aqp-1+2

dsRNA treatment

and 19.2 mm2 (luc), however, flukes treated with Ov-aqp1, Ov-aqp-2 and Ov-aqp-1 + 2 dsRNAs had significantly reduced average body sizes of15.55 (p ≤ 0.001), 16.05 (p ≤ 0.001) and 16.3 mm2 (p ≤ 0.001), respectively (Figure 6B). The mean body sizes of mock and luc control groups were not significantly different. Tissue sections of Ov-aqp-2 dsRNA-treated worms maintained in isotonic culture medium (RPMI) for 7 days showed that the parenchyma was more condensed than control worms (Figure 7), and parenchyma pore sizes were reduced in Ov-aqp-2 dsRNA-treated worms (Figure 7).

Discussion O. viverrini has a complex life cycle with each developmental stage inhabiting dissimilar environments. Freeliving stages including miracidia and cercariae live in fresh water for short periods before invading their respective snail and fish hosts. Parasitic stages including sporocysts, rediae, metacercariae and adult flukes live in the snail (sporocyst, redia), fish (metacercaria) and mammalian (adult fluke) hosts [31]. Like other eukaryotes, trematodes need to regulate water movement into and out of their cells. Transcriptomic and proteomic surveys of O. viverrini indicate abundant expression of a small number of AQPs, transporters that have been well characterized as water channels in other species [16,17,22].

Figure 6 Suppression of Ov-aqp-1, Ov-aqp-2 and Ov-aqp-1 + 2 gene expression results in an inability to maintain osmotic flux. Size of cultured adult parasites after treatment with control (luc), Ov-aqp-1, Ov-aqp-2 and Ov-aqp-1 + 2dsRNAs when incubated in RPMI medium (A) or hypotonic medium (distilled water; B). Each point represents an individual parasite; horizontal line indicates the mean for that group; n = 8 per treatment. p ≤ 0.001 is indicated with ***.

In this study, cDNAs encoding aqp genes from O. viverrini were identified and cloned. AQP is a membrane integral protein that forms pores that facilitate passage of water and other small solutes through cell membranes. They function in cells where the membrane is intact and there is a need for a specific amount of water in the cell to maintain osmotic pressure. The structure of AQPs is well known [10]. Ov-AQP-1 and AQP-2 exhibited the typical AQP structure of six transmembrane domains with two NPA boxes that allow pore formation [9,32]. The AQP family is comprised of three groups based on permeability properties. The first group contains the strict AQPs that allow only water to permeate the channel. Members of the second group conduct other small solutes, such as glycerol and urea; these transporters are termed aquaglyceroporins [9,33]. The third family of AQPs has a unique group of conserved residues surrounding the NPA

Thanasuwan et al. Parasites & Vectors 2014, 7:224 http://www.parasitesandvectors.com/content/7/1/224

Page 9 of 11

Day 1

Parenchyma Parenchyma Parenchyma

Ov-AQP-2

LUC

Control

Day 3 Parenchyma Parenchyma Parenchyma

Control

LUC

Ov-AQP-2 Parenchyma

Parenchyma

Parenchyma

Control

LUC

Day 7

Ov-AQP-2

Figure 7 Cross-section of adult Opisthorchis viverrini after RNAi-mediated suppression of Ov-aqp-2; stained with hematoxylin and eosin. Adult O. viverrini were electroporated with Ov-aqp-2 dsRNA and were maintained thereafter in isotonic medium (RPMI), and observed on days 1, 3 and 7. Regions of trematode parenchyma are indicated by the boxes.

boxes; this group has numerous names including the S-aquaporins (superaquaporins), aquaporins with unusual (or deviated) NPA boxes, and subcellular (or sip-like) aquaporins [11,33,34]. To date, the structures of nematode and trematode AQPs have generally consisted of six transmembrane domains and exhibited two NPA domaincontaining hemipores that overlap to form a pore between the membrane bilayer,a defining characteristic of AQPs [13,35]. Some aquaporins from the liver fluke Fasciola gigantica do not show conserved two NPA motifs, where the first NPA box is replaced by TAA although the second NPA box is conserved [36]. It is notable, therefore, that the O.viverrini AQPs reported here showed higher sequence similarity with the aquaglyceroporin group (human AQP3, -7, −9 and −10) [9,10] as well as SmAQP from the human blood fluke, S. mansoni [12] and PfAQP from Plasmodium falciparum which exhibit a functional channel facilitating water and glycerol [37] compared to the traditional AQP group (human AQP-1, −2, −5 and −6) [9,10],

suggesting that permeants beyond water or glycerol might be transported into/from cells of O. viverrini. Transcripts encoding aqp were detected in each of the developmental stages examined here, suggesting a requirement for these proteins throughout the entire life cycle of the parasite. The highest expression levels of Ov-aqp-1 and Ov-aqp-2 were detected in metacercariae and eggs, respectively, implying key roles in the physiology of the infectious process for AQP-1 in particular. To further address the functional roles of Ov-AQPs, we used RNAi to suppress their expression and observed that both Ov-aqp-1 and aqp-2 were essential for controlling osmosis in the adult fluke. Likewise, SmAQP is important for both entry and exit of water and other solutes in schistosomes [12]. Accordingly, AQP is being considered as a potential vaccine candidate for combating Schistosoma japonicum infections due to the presence of numerous B cell epitopes that are predicted to be accessible to host antibodies [38].

Thanasuwan et al. Parasites & Vectors 2014, 7:224 http://www.parasitesandvectors.com/content/7/1/224

Histological analysis of Ov-aqp-2 dsRNA-treated flukes showed that the average pore size in the parenchymal tissue was reduced compared to control worms. This finding, coupled with the inability of aqp-2-dsRNA treated worms to maintain osmotic pressure, implies that OvAQP-2 is integral to entry and egress of water from the tissues of O. viverrini. Phylogenetic analysis revealed that the O. viverrini AQPs belong to the aquaglyceroporin group, suggesting that their roles extend to substrates beyond other small solutes [37,39,40]. However, in this study we investigated only regulation of water into the body wall. Several authors have investigated other functions of AQPs in schistosomes, including the entry/exit of drugs, highlighting the importance of the trematode tegument in osmoregulation, and the uptake and excretion of metabolic waste and drugs [12,14]. In like fashion, Ov-AQPs might play similar roles in drug uptake and waste excretion, and we propose to explore these potential functions in the future. In conclusion, AQPs likely perform essential roles in the movement of water between the external environment and the liver fluke. Drugs that counteract mammalian AQPs are in development [41,42], and accordingly liver fluke AQPs represent potential new targets for interventions to control and treat this carcinogenic pathogen.

Conclusion Aquaporins (AQPs) belong to the major intrinsic protein superfamily of integral plasma membrane channel proteins that selectively transport water across cell membranes. AQPs play key roles as water and ion transport channels through the tegument of trematodes. Here we investigated two discrete forms of AQP mRNAs from the adult O. viverrini. Phylogenetic analysis indicated that Ov-AQPs belong to aquaglyceroporin-like water channel proteins characterized from parasites, which may be integral to the specialized requirements of water and solute control during parasitism. Ov-aqps transcripts were detected in all developmental stages of O. viverrini. RNA interference using electroporated dsRNA suppressed transcript levels of Ov-aqps in adult worms by 58-99%, and disabled water transport in adult flukes.The apparently pivotal roles of Ov-AQP in solute homeostasis at the fluke surface suggest that deeper investigation will be informative for the pathophysiology of O. viverrini, and may uncover intervention targets, particularly in view of the singularly notable predilection of this pathogen for residence within the biliary tree. Competing interests The authors declare that they have no competing interests. Authors’ contributions ST, PJB, AL and TL conceived and designed the study. ST and SP performed experiments. ST, SK and TL performed data analysis and interpretation. ST, PJB, AL and TL wrote the manuscript. All authors reviewed the draft and

Page 10 of 11

accepted it for submission. All authors read and approved the final version of the manuscript. Acknowledgements This work was supported by the Higher Education Research Promotion and National Research University Project of Thailand, Office of the Higher Education Commission, Grant under the program Strategic Scholarships for Frontier Research Network for the Ph.D. Program Thai Doctoral degree from the Office of the Higher Education Commission, Thailand and from the National Institute of Allergy and Infectious Disease (NIAID) of the United States National Institutes of Health (NIH) award number P50AI098639 (Tropical Medicine Research Center). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIAID or the NIH.AL is supported by a NHMRC principal research fellowship. Author details 1 Department of Parasitology, Faculty of Medicine, Khon Kaen University, Khon Kaen 40002, Thailand. 2Department of Microbiology, Immunology and Tropical Medicine, Research Center for Neglected Diseases of Poverty, School of Medicine & Health Sciences, George Washington University, Washington, DC 20037, USA. 3Centre for Biodiscovery and Molecular Development of Therapeutics, Australian Institute of Tropical Health and Medicine, James Cook University, Cairns, Queensland 4878, Australia. 4Liver Fluke and Cholangiocarcinoma Research Center, Faculty of Medicine, Khon Kaen University, Khon Kaen 40002, Thailand. Received: 18 March 2014 Accepted: 9 May 2014 Published: 14 May 2014 References 1. Sithithaworn P, Andrews RH, Nguyen VD, Wongsaroj T, Sinuon M, Odermatt P, Nawa Y, Liang S, Brindley PJ, Sripa B: The current status of opisthorchiasis and clonorchiasis in the Mekong Basin. Parasitol Int 2012, 61(1):10–16. 2. Sripa B: Global burden of food-borne trematodiasis. Lancet Infect Dis 2012, 12(3):171–172. 3. Sripa B, Brindley PJ, Mulvenna J, Laha T, Smout MJ, Mairiang E, Bethony JM, Loukas A: The tumorigenic liver fluke Opisthorchis viverrini–multiple pathways to cancer. Trends Parasitol 2012, 28(10):395–407. 4. Geary TG, Woods DJ, Williams T, Nwaka S: Target identification and mechanism-based screening for anthelmintics: Application of veterinary antiparasitic research programmes to search for new antiparasitic drugs for human indications. In Antiparasitic and Antibacterial Drug Discovery: From Molecular Targets to Drug Candidates. Edited by Selzer PM. Weinheim: Wiley-Blackwell; 2009:3–15. 5. Taylor CM, Wang Q, Rosa BA, Huang SC, Powell K, Schedl T, Pearce EJ, Abubucker S, Mitreva M: Discovery of anthelmintic drug targets and drugs using chokepoints in nematode metabolic pathways. PLoS Pathog 2013, 9(8):e1003505. 6. Kaminsky R, Ducray P, Jung M, Clover R, Rufener L, Bouvier J, Weber SS, Wenger A, Wieland-Berghausen S, Goebel T, Gauvry N, Pautrat F, Skripsky T, Froelich O, Komoin-Oka C, Westlund B, Sluder A, Maser P: A new class of anthelmintics effective against drug-resistant nematodes. Nature 2008, 452(7184):176–180. 7. Beitz E: Aquaporins from pathogenic protozoan parasites: structure, function and potential for chemotherapy. Biol Cell 2005, 97(6):373–383. 8. Zeuthen T: How water molecules pass through aquaporins. Trends Biochem Sci 2001, 26(2):77–79. 9. Zardoya R: Phylogeny and evolution of the major intrinsic protein family. Biol Cell 2005, 97(6):397–414. 10. Verkman AS, Mitra AK: Structure and function of aquaporin water channels. Am J Physiol Renal Physiol 2000, 278(1):F13–F28. 11. Ishibashi K: Aquaporin subfamily with unusual NPA boxes. Biochim Biophys Acta 2006, 1758(8):989–993. 12. Faghiri Z, Skelly PJ: The role of tegumental aquaporin from the human parasitic worm, Schistosoma mansoni, in osmoregulation and drug uptake. Faseb J 2009, 23(8):2780–2789. 13. Song J, Mak E, Wu B, Beitz E: Parasite aquaporins: Current developments in drug facilitation and resistance. Biochim Biophys Acta 2014, 1840(5):1566–1573.

Thanasuwan et al. Parasites & Vectors 2014, 7:224 http://www.parasitesandvectors.com/content/7/1/224

14. Faghiri Z, Camargo SM, Huggel K, Forster IC, Ndegwa D, Verrey F, Skelly PJ: The tegument of the human parasitic worm Schistosoma mansoni as an excretory organ: the surface aquaporin SmAQP is a lactate transporter. PLoS One 2010, 5(5):e10451. 15. Da’dara A, Krautz-Peterson G, Faghiri Z, Skelly PJ: Metabolite movement across the schistosome surface. J Helminthol 2012, 86(2):141–147. 16. Laha T, Pinlaor P, Mulvenna J, Sripa B, Sripa M, Smout MJ, Gasser RB, Brindley PJ, Loukas A: Gene discovery for the carcinogenic human liver fluke, Opisthorchis viverrini. BMC Genomics 2007, 8:189. 17. Mulvenna J, Sripa B, Brindley PJ, Gorman J, Jones MK, Colgrave ML, Jones A, Nawaratna S, Laha T, Suttiprapa S, Smout MJ, Loukas A: The secreted and surface proteomes of the adult stage of the carcinogenic human liver fluke Opisthorchis viverrini. Proteomics 2010, 10(5):1063–1078. 18. Sithithaworn P, Pipitgool V, Srisawangwong T, Elkins DB, Haswell-Elkins MR: Seasonal variation of Opisthorchis viverrini infection in cyprinoid fish in north-east Thailand: implications for parasite control and food safety. Bull World Health Organ 1997, 75(2):125–131. 19. Sripa B, Kaewkes S: Gall bladder and extrahepatic bile duct changes in Opisthorchis viverrini-infected hamsters. Acta Trop 2002, 83(1):29–36. 20. Piratae S, Tesana S, Jones MK, Brindley PJ, Loukas A, Lovas E, Eursitthichai V, Sripa B, Thanasuwan S, Laha T: Molecular characterization of a tetraspanin from the human liver fluke, Opisthorchis viverrini. PLoS Negl Trop Dis 2012, 6(12):e1939. 21. Kaewkes S, Kaewkes W, Boonmars T, Sripa B: Effect of light intensity on Opisthorchis viverrini cercarial shedding levels from Bithynia snails–a preliminary study. Parasitol Int 2012, 61(1):46–48. 22. Young ND, Campbell BE, Hall RS, Jex AR, Cantacessi C, Laha T, Sohn WM, Sripa B, Loukas A, Brindley PJ, Gasser RB: Unlocking the transcriptomes of two carcinogenic parasites, Clonorchis sinensis and Opisthorchis viverrini. PLoS Negl Trop Dis 2010, 4(6):e719. 23. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25(17):3389–3402. 24. Hall T: BioEdit is a biological sequence alignment editor written for Windows 95/98/NT/2000/XP. Carlsbad, CA: Nucl Acids Symp Ser; 1999. 25. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 1987, 4(4):406–425. 26. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S: MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol 2013, 30(12):2725–2729. 27. Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc 2008, 3(6):1101–1108. 28. Rinaldi G, Morales ME, Cancela M, Castillo E, Brindley PJ, Tort JF: Development of functional genomic tools in trematodes: RNA interference and luciferase reporter gene activity in Fasciola hepatica. PLoS Negl Trop Dis 2008, 2(7):e260. 29. Sripa B, Kaewkes S: Localisation of parasite antigens and inflammatory responses in experimental opisthorchiasis. Int J Parasitol 2000, 30(6):735–740. 30. Wang X, Chen W, Huang Y, Sun J, Men J, Liu H, Luo F, Guo L, Lv X, Deng C, Zhou C, Fan Y, Li X, Huang L, Hu Y, Liang C, Hu X, Xu J, Yu X: The draft genome of the carcinogenic human liver fluke Clonorchis sinensis. Genome Biol 2011, 12(10):R107. 31. Wykoff DE, Harinasuta C, Juttijudata P, Winn MM: Opisthorchis viverrini in Thailand–the Life Cycle and Comparison with O. Felineus. J Parasitol 1965, 51:207–214. 32. Benga G: The first discovered water channel protein, later called aquaporin 1: molecular characteristics, functions and medical implications. Mol Aspects Med 2012, 33(5–6):518–534. 33. Benga G: On the definition, nomenclature and classification of water channel proteins (aquaporins and relatives). Mol Aspects Med 2012, 33(5–6):514–517. 34. Ishibashi K: Aquaporin superfamily with unusual npa boxes: S-aquaporins (superfamily, sip-like and subcellular-aquaporins). Cell Mol Biol 2006, 52(7):20–27. 35. Loukas A, Hunt P, Maizels RM: Cloning and expression of an aquaporinlike gene from a parasitic nematode. Mol Biochem Parasitol 1999, 99(2):287–293. 36. Geadkaew A, Von Bulow J, Beitz E, Grams SV, Viyanant V, Grams R: Functional analysis of novel aquaporins from Fasciola gigantica. Mol Biochem Parasitol 2011, 175(2):144–153.

Page 11 of 11

37. Hansen M, Kun JF, Schultz JE, Beitz E: A single, bi-functional aquaglyceroporin in blood-stage Plasmodium falciparum malaria parasites. J Biol Chem 2002, 277(7):4874–4882. 38. Song J, He QF: Bioinformatics analysis of the structure and linear B-cell epitopes of aquaporin-3 from Schistosoma japonicum. Asian Pac J Trop Med 2012, 5(2):107–109. 39. Liu F, Lu J, Hu W, Wang SY, Cui SJ, Chi M, Yan Q, Wang XR, Song HD, Xu XN, Wang JJ, Zhang XL, Zhang X, Wang ZQ, Xue CL, Brindley PJ, McManus DP, Yang PY, Feng Z, Chen Z, Han ZG: New perspectives on host-parasite interplay by comparative transcriptomic and proteomic analyses of Schistosoma japonicum. PLoS Pathog 2006, 2(4):e29. 40. Bassarak B, Uzcategui NL, Schonfeld C, Duszenko M: Functional characterization of three aquaglyceroporins from Trypanosoma brucei in osmoregulation and glycerol transport. Cell Physiol Biochem 2011, 27(3–4):411–420. 41. Wacker SJ, Aponte-Santamaria C, Kjellbom P, Nielsen S, De Groot BL, Rutzler M: The identification of novel, high affinity AQP9 inhibitors in an intracellular binding site. Mol Membr Biol 2013, 30(3):246–260. 42. Martins AP, Marrone A, Ciancetta A, Galan Cobo A, Echevarria M, Moura TF, Re N, Casini A, Soveral G: Targeting aquaporin function: potent inhibition of aquaglyceroporin-3 by a gold-based compound. PLoS One 2012, 7(5):e37435. doi:10.1186/1756-3305-7-224 Cite this article as: Thanasuwan et al.: Suppression of aquaporin, a mediator of water channel control in the carcinogenic liver fluke, Opisthorchis viverrini. Parasites & Vectors 2014 7:224.

Submit your next manuscript to BioMed Central and take full advantage of: • Convenient online submission • Thorough peer review • No space constraints or color figure charges • Immediate publication on acceptance • Inclusion in PubMed, CAS, Scopus and Google Scholar • Research which is freely available for redistribution Submit your manuscript at www.biomedcentral.com/submit