Calmodulin disruption impacts growth and motility in juvenile liver fluke
© McCammick et al. 2016
Received: 29 October 2015
Accepted: 19 January 2016
Published: 27 January 2016
Deficiencies in effective flukicide options and growing issues with drug resistance make current strategies for liver fluke control unsustainable, thereby promoting the need to identify and validate new control targets in Fasciola spp. parasites. Calmodulins (CaMs) are small calcium-sensing proteins with ubiquitous expression in all eukaryotic organisms and generally use fluctuations in intracellular calcium levels to modulate cell signalling events. CaMs are essential for fundamental processes including the phosphorylation of protein kinases, gene transcription, calcium transport and smooth muscle contraction. In the blood fluke Schistosoma mansoni, calmodulins have been implicated in egg hatching, miracidial transformation and larval development. Previously, CaMs have been identified amongst liver fluke excretory-secretory products and three CaM-like proteins have been characterised biochemically from adult Fasciola hepatica, although their functions remain unknown.
In this study, we set out to investigate the biological function and control target potential of F. hepatica CaMs (FhCaMs) using RNAi methodology alongside novel in vitro bioassays.
Our results reveal that: (i) FhCaMs are widely expressed in parenchymal cells throughout the forebody region of juvenile fluke; (ii) significant transcriptional knockdown of FhCaM1-3 was inducible by exposure to either long (~200 nt) double stranded (ds) RNAs or 27 nt short interfering (si) RNAs, although siRNAs were less effective than long dsRNAs; (iii) transient long dsRNA exposure-induced RNA interference (RNAi) of FhCaMs triggered transcript knockdown that persisted for ≥ 21 days, and led to detectable suppression of FhCaM proteins; (iv) FhCaM RNAi significantly reduced the growth of juvenile flukes maintained in vitro; (v) FhCaM RNAi juveniles also displayed hyperactivity encompassing significantly increased migration; (vi) both the reduced growth and increased motility phenotypes were recapitulated in juvenile fluke using the CaM inhibitor trifluoperazine hydrochloride, supporting phenotype specificity.
These data indicate that the Ca2+-modulating functions of FhCaMs are important for juvenile fluke growth and movement and provide the first functional genomics-based example of a growth-defect resulting from gene silencing in liver fluke. Whilst the phenotypic impacts of FhCaM silencing on fluke behaviour do not strongly support their candidature as new flukicide targets, the growth impacts encourage further consideration, especially in light of the speed of juvenile fluke growth in vivo.
Fasciola spp. liver flukes cause fasciolosis, a disease that impacts heavily on farmed animal health and food security, and is of growing importance as a human neglected tropical disease [1–3]. There are only a handful of available flukicides with which to combat these infections, and their utility is being eroded by the increasing occurrence of drug resistance in both veterinary and human infections [4, 5]. Triclabendazole is the only flukicide with significant activity against juvenile liver fluke , the life-cycle stage that causes host tissue damage during its migration from the gut lumen, through the intestinal mucosa and hepatic parenchyma to the bile ducts. Novel flukicides, especially those with activity against juvenile fluke, are a pressing need.
Calcium ions (Ca2+) are ubiquitous messengers with an array of roles in physiological and metabolic pathways. External stimuli or signals often lead to variations in intracellular Ca2+ levels that regulate processes such as transcription, neurotransmitter release and muscle function . The range of responses that can result from Ca2+ signalling are a direct result of the temporal and spatial organisation of the signal . In many cases, particularly longer acting responses, effector proteins cannot bind Ca2+ directly and require an intermediate Ca2+ sensing protein in order to initiate a biological response.
Calmodulin (CaM) is the most studied Ca2+ sensing protein in mammals and can act as a monomer, or as a subunit in multimeric proteins . Within its dumbbell shape it possesses two globular domains, each with the ability to bind two Ca2+ ions (four binding sites per CaM), connected by a flexible linker region. Normally, cytosolic Ca2+ levels are low, but even a small rise in cytosolic [Ca2+] can trigger Ca2+ binding and a conformational change in CaM that facilitates the interaction between the Ca2+/CaM complex and target proteins . The cooperative binding properties of CaM mean that small increases in cytosolic Ca2+ can be translated, via CaM and other calcium sensors, to a multitude of biological processes such as transcription, muscle activity, neurotransmission, cytoskeletal assembly/reorganisation and cell death [10, 11]. The key role played by CaM in diverse biological processes makes it an appealing target for the control of pathogens. Indeed, the disruption of Ca2+ regulation is believed to be key to the mode of action of praziquantel, used for the control of Schistosoma mansoni and various cestodiases .
Despite this interest, CaMs in flatworms are still quite poorly understood. CaM-like proteins have been identified in several important fluke species. Two CaMs have been identified and characterised from S. mansoni , showing high sequence identity (>98 %) to CaMs from mammals. Silencing these schistosome CaMs using RNA interference (RNAi) resulted in shorter worms in larval stages and somatic muscle contraction/dilation in adults, while CaM antagonists inhibited miracidial transformation in a concentration dependent manner [13, 14]. Previously, a CaM inhibitor was shown to inhibit egg hatching and miracidial transformation in S. mansoni [15, 16]. A Clonorchis sinensis CaM was shown to bind Ca2+ and to be expressed in diverse developmental stages including the egg, metacercaria and adult where it was localised in the intestine, pharynx and tegument . A divergent CaM-like protein has been reported in Schistosoma japonicum  and three CaM-like proteins have been identified and biochemically described in Fasciola hepatica [19, 20]. The first, FhCaM1, displays high identity to mammalian CaM with only two amino acids differing. FhCaM2 and FhCaM3 are more divergent, displaying lower sequence similarity, as well as distinct biochemical and ion binding properties. This prompted the hypothesis that they are involved in different pathways or perform roles distinct from that of FhCaM1 in the organism. All three fluke CaMs were predicted to adopt a calmodulin-like fold and shown to bind Ca2+. Immunocytochemistry failed to detect FhCaM1 but localised FhCaM2 and FhCaM3 to eggs and vitelline cells in adult worms [19, 20].
This study builds upon our description of RNAi dynamics in juvenile fluke , employing RNAi methods to investigate the biological functions of FhCaM1, FhCaM2 and FhCaM3 in juvenile F. hepatica. We show that FhCaM2 and FhCaM3 are expressed widely in parenchymal cells of juvenile fluke and that FhCaM RNAi inhibits growth and modulates motility, phenotypes that can be recapitulated using a CaM antagonist drug. These data provide the first functional data for liver fluke calmodulins, including the first description of an RNAi-induced growth phenotype in liver fluke, suggesting that drugs capable of targeting CaM’s role in growth/development might represent appealing new additions to our flukicidal arsenal.
Parasites and maintenance
F. hepatica metacercariae (Oregon strain) were purchased from Baldwin Aquatics Inc. (Monmouth, Or, USA) and stored in distilled water at 4 °C. These metacercariae were supplied without their outer cyst walls, with excystment performed as described previously . Worms were transferred, as they excysted, into RPMI 1640 supplemented with antibiotic/antimycotic (100 U/ml penicillin, 100 μg/ml streptomycin and 25 μg/ml amphotericin B, Life Technologies, Carlsbad, USA), before further processing as described below. At all stages of these experiments, excysted juvenile fluke were handled under aseptic conditions and maintained at 37 °C in a humidified 5 % CO2 atmosphere.
Genomic organisation of CaM encoding genes
The F. hepatica genome  was queried using tBLASTn analysis for the translated sequences of FhCaM1, FhCaM2, and FhCaM3. Returned hits were assessed for bit score and E-value and manually assembled from the parent scaffold using Artemis to decipher exon organisation . Contigs were aligned pairwise with the corresponding query sequence to ensure identity.
Freshly excysted fluke were fixed for 4 h at room temperature in 4 % paraformaldehyde in 0.1 M Phosphate Buffered Saline (PBS prepared from tablets; Sigma Aldrich, St Louis, USA). Primary antisera were raised in rabbits against each of the following: recombinant human-CaM (structurally similar to FhCaM1 and differs by two amino acids; antiserum shown previously to cross-react with FhCaM1; Sigma-Aldrich); FhCaM2; and, FhCaM3 [19, 20]. Incubations in primary antisera were performed at 1/200 dilution (diluted in antibody diluent, AbD; 0.1 M PBS, 0.1 % (v/v) TritonX-100, 0.1 % (w/v) bovine serum albumin, 0.1 % (w/v) NaN3) for two days. Samples were then incubated consecutively in secondary antibody (fluorescein isothiocyanate (FITC)-labelled goat anti-rabbit; 1:100) overnight, and then overnight in tetramethylrhodamine isothiocyanate (TRITC)-labelled phalloidin (200 ng/ml), both diluted in AbD. Incubations were performed at 4 °C unless stated otherwise, and at least three washes in 1 M NaCl in PBS, 0.1 % TritonX-100, 0.2 % Tween-20, were performed between incubations. Samples were mounted in Vectashield (Vector Laboratories, Peterborough, UK) and viewed on a Leica TCS SP5 confocal scanning laser microscope.
CaM RNAi methodology
CaM transcript dsRNA templates were amplified from PCR amplicons representing the coding sequence of each CaM gene. These coding sequences were amplified using the following primers, designed against the GenBank accession numbers provided: FhCaM1 [AM412546], forward 5’-atggctgatcaactcacagaaga-3’, reverse 5’-tcattttgcagtcatcattttcac-3’; FhCaM2 [AM412547], forward 5’-ggtgcaccatctctcacaag-3’, reverse 5’-ttgcgcagaagagtcaagaa-3’; FhCaM3 [JQ792169], forward 5’-atggctgactttgatatggatca-3’, reverse 5’-ctacaggtattcgtatgcaaattcc-3’. Templates for dsRNA constructs were amplified from these primary amplicons in a nested PCR reaction, using primers labelled (on 5’ end) with a T7 RNA polymerase promoter sequence (5’-taatacgactcactatagggt-3’), to create templates labelled with T7 in either sense or antisense directions (FhCaM1, forward 5’- cagaagagcagatcgcagaa-3’, reverse 5’-aattcggggaaatcaatggt-3’; FhCaM2, forward 5’- agatgccacggaactgaaag-3’, reverse 5’-aaatgtacccgtcgccatta-3’; FhCaM3, forward 5’- atgacgctgtacgacaccaa-3’, reverse 5’-agtggtgggcacaacgatac-3’). All PCR amplicons were sequence confirmed . This approach generated dsRNA templates for FhCaM1, 2, and 3 of 190 bp, 221 bp and 210 bp, respectively. Exogenously derived (i.e. non-worm) negative control dsRNA templates (bacterial neomycin phosphotransferase; [U55762]) were generated from the cloning vector pEGFP-NI (Clontech Laboratories Inc., Saint-Germain-en-Laye, France), using the following primers: forward, 5’-ggtggagaggctattcggct-3’; reverse: 5’-ccttcccggttcagtgacaa-3’. These templates were used to generate double stranded (ds)RNAs as described in McVeigh et al. .
Short interfering (si)RNAs (27mer) for each FhCaM were designed and synthesised by ; a randomised siRNA control was also used. Target sequences for these siRNAs were as follows: FhCaM1, 5’-augacugucccaagcuccuucguggua-3’; FhCaM2, 5’-uauguaagacugaucucaauucuuccg-3’; FhCaM3, 5'-ucccugucacuaacccuccucuuuguu-3’; and negative control siRNA sequence siCTRL, 5’-cttcctctctttctctcccttgtga-3’.
In order to expose juvenile fluke to RNAi trigger molecules, fluke were soaked (20 NEJs per replicate for transcriptional analysis, 10 NEJs per replicate for phenotype assays and 50 NEJs per replicate for western blot analysis) in 50 μl of 100 ng/μl (dsRNA), or 50 ng/μl (siRNA) solution of FhCaM ds/siRNA dissolved in RPMI 1640, alongside both untreated (no ds/siRNA) and negative control (dsCTRL/siCTRL) treatments, with all soaks performed at least in triplicate. Soaks were performed for 4 h at 37 °C, in 5 % CO2. After this 4 h incubation, 250 μl maintenance media was added (RPMI 1640, containing antibiotic/antimycotic and supplemented with 20 % foetal bovine serum (FBS)). Worms were maintained in this medium for up to 21 days with a change of media once every 2–3 days. During media changes worms were observed using light microscopy for obvious phenotypic changes (aberrant morphology or motility), or imaged for measurements as described below.
Transcript abundance was assessed in NEJs maintained in the presence (+) and absence (−) of serum in the maintenance media at 0-, 1-, 3-, 5- and 7-days. Note that 97 % of juvenile fluke are viable by day 7 in RPMI and 100 % are viable in RPMI supplemented with serum. Parasites remain viable and motile for several weeks under these conditions (data not shown). For RNAi experiments, transcript abundance was assessed following maintenance at 3-, 14- and 21-days post dsRNA (or siRNA) exposure, by quantitative PCR (qPCR). Poly-adenylated mRNA, extracted from treatment groups (20 NEJs per replicate) using Dynabeads mRNA direct Kit (Life Technologies, Carlsbad, USA), was used to synthesise cDNA (High Capacity RNA-to-cDNA kit, Life Technologies) following DNase treatment (Turbo DNA-free, Life Technologies). qPCRs were performed on a Rotor-Gene Q 5-plex HRM PCR system (Qiagen), using the following gene-specific primers: FhCaM1, forward 5’- tggtaacgggaccattgatt-3’, reverse 5’- ttatcgaacacacggaatgc-3’; FhCaM2, forward 5’- atcgcggcaatcttgataat-3’, reverse 5’-ttgcgcagaagagtcaagaa-3’; FhCaM3, forward 5’- agcttcggaactgcacaaat-3’, reverse 5’-cgttgtcaccatctttgtcg-3’. Target primers were amplified alongside the glyceraldehyde 3-phosphate dehydrogenase (GAPDH) reference transcript (FhGAPDH [GenBank; AY005475]), forward 5’-ggctgtgggcaaagtcattc-3’, reverse 5’- agatccacgacggaaacatca-3’). Primers were added to 8 μl PCR reactions at a final concentration of 1.25 μM, including cDNA diluted 1:1 with water, and FastStart SYBR Green Master (Roche Applied Science). Cycling parameters employed were a 10 min “hot-start” at 95 °C, followed by 40 cycles of 95 °C 10 s, 60 °C 15 s, 72 °C 30 s. All PCRs were performed in triplicate, and included no-template controls and melt-curve analyses as standard. Relative expression analysis employed Pfaffl’s Augmented ΔΔCt method  (which normalises expression in each sample relative to the untreated control, standardised to a GAPDH reference amplicon), with amplification efficiencies of individual reactions calculated using Real Time PCR Miner software [28, 29]. This analysis method produces a ratio of target transcript abundance relative to the untreated control, which was then converted into percentage change (i.e. where a value of 100 % represents no change). These data (mean % transcript abundance relative to FhGAPDH reference ± SEM) are plotted in bar graphs in the qPCR figures. All datasets were normally distributed according to the Kolmogorov-Smirnov test and were, therefore, analysed using ANOVA with Tukey's or Dunnett’s multiple comparison post hoc test. Statistical significance was determined relative to the effects of negative control treatments (dsCTRL or siCTRL) on target gene expression.
Analysis of CaM protein
Protein abundance was assessed after RNAi treatments and maintenance in the presence of serum at days −7 and −14 post dsRNA, using western blot methods. Protein was extracted from tissue matched (50 NEJs) homogenised treatment groups in 50 μl radioimmunoprecipitation assay (RIPA) buffer (Sigma-Aldrich), with blots performed as described previously . Following protein transfer onto nitrocellulose membrane, the top half of each membrane (relatively high molecular weight proteins) was probed with rabbit anti-actin (1:400) (Sigma Aldrich) as a loading control. The bottom half (relatively low molecular weight proteins) was probed with rabbit anti-HsCaM1 (1:1000), rabbit anti-FhCaM2 (1:1000) or rabbit anti-FhCaM3 (1:1000); these incubations were performed in parallel. Secondary antibody was an alkaline phosphatase conjugated goat-anti-rabbit (Sigma Aldrich), used at 1:1000 dilution. Development was performed in nitro blue tetrazolium chloride/5-bromo-4-chloro-3-indolyl-phosphate, toluidine-salt (NBT/BCIP) solution.
We used the CaM antagonists, trifluoperazine hydrochloride (TFP) and N-(6-aminohexyl)-5-chloro-1-naphthalenesulfonamide hydrochloride (W7) (both Sigma Aldrich) . Drugs were dissolved in 100 % dimethylsulfoxide (DMSO) and added to maintenance media such that the final DMSO concentration was 0.1 %. Two methods of administration were employed: (i) short term exposure performed in RPMI 1640 for 18 h at final concentrations 10, 5, 2 and 0.5 μM, including vector (DMSO) only control and untreated control; and (ii) long term exposure performed in RPMI + 20 % FBS over 7 days in constant presence of TFP with final concentrations at 5, 1, 0.5, 0.1, 0.05 μM plus controls as outlined above.
Following RNAi and drug treatments, we assayed for changes in both motility/migration and growth phenotypes. All analyses were performed in an operator-blinded fashion, on coded samples. Analysis of worm motility was performed on drug-treated (18 h) worms by capturing videos (30 s), and manually counting the numbers of worms displaying full, “normal” motility (i.e. waves of circular muscle contraction travelling from anterior to posterior) during the analysis period. Analysis of migration was performed on both drug treated and RNAi-treated worms. Migratory capacity under these experimental conditions was assessed using an agar dispersal assay, in which worms were placed into an agar substrate (0.2 % cell culture grade agar (Sigma Aldrich) in HEPES-buffered Ringer’s (123 mM NaCl, 5 mM KCl, 1.6 mM CaCl2, 20 mM HEPES, 11.1 mM d-glucose; pH 7.4)) at the centre of a petri dish. The proportion of worms migrating beyond an empirically determined radius (5 mm) was then recorded over a 3 h time-course using a measuring grid under a dissecting microscope.
Growth analyses were performed on both drug and RNAi treated worms (and controls). Images captured using a Leica MZ 12.5 stereomicroscope and Unibrain Fire-i digital camera were subject to measurements of length and area using the measurement tools in ImageJ software . Growth and motility assay data were analysed using one-way ANOVA with Dunnett’s post hoc test, analysing significance of treatments relative to the negative control ds/siRNA treated worms (RNAi experiments), or vehicle control treated worms (drug experiments).
Results and Discussion
This study reports functional characterisation of three CaM-like proteins in F. hepatica juveniles. Although the in vitro biochemistry of the three FhCaM recombinants, and their adult localisations have been reported previously [19, 20], functional genomics analyses have not been performed on CaMs in any F. hepatica life-cycle stage. Motivated by reports of CaM RNAi phenotypes in other helminths [13, 33] we focused here on assessing CaM function in juvenile fluke, given the amenability of this life-cycle stage to our target validation platform for fluke that encompasses RNAi methods [21, 22, 34] and diverse functional assay platforms. A primary motivation for focusing on CaM, a protein intimately involved in cellular Ca2+ signalling processes, is the fact that Ca2+ signalling and regulation play key roles in neuromuscular function, which represents a well-established repository of targets for empirically-discovered anthelmintics . Indeed, the disruption of Ca2+ regulation by praziquantel, the drug of choice for the treatment of schistosome and tapeworm infections, directly highlights the validity of interventions disrupting Ca2+ signalling/regulation. Whilst previous work on liver fluke CaMs focussed on in vitro function and adult stage expression, our focus was CaM expression and biological importance as a measure of control target validity in juvenile fluke. It is well established that effective flukicides need to target both the migrating juvenile and adult fluke stages. In this way our data inform CaM expression, function, and potential as a therapeutic target for liver fluke control.
Analysis of hits collected from querying the genomic dataset with pre-confirmed FhCaM1, FhCaM2 and FhCaM3 protein sequences revealed that the FhCaM1 gene exists in at least two intronless copies on scaffold 932: one full length with an S18F substitution, and one missing the initial 18 amino acids due to an undefined region on the scaffold (Additional file 1: Figure S1). A third copy of FhCaM1 is present on scaffold 210 (exists as 4 exons) with a stop codon (TAG) in place of the initiating methionine. This could represent a sequencing error or the premature stop codon may indicate that the third FhCaM1 copy is a pseudogene that would not be translated but could play alternative regulatory roles. FhCaM2 is present as a single copy on scaffold 2277. Four exons have been identified with possible variants for the amino acid at position 60 (N60S or N60D). The initiating methionine is also absent from the genome sequence of FhCaM2, likely due to an artefact of genome assembly as this sequence has been characterised and the initiating methionine confirmed to be present experimentally . FhCaM3 is also represented by a single copy on scaffold 1560 where the full gene presents as three exons with an apparent duplication of exon 3 (amino acids 71–131) within the scaffold. This duplication likely represents a sequence assembly artefact, as the surrounding non-CaM coding sequence is identical.
These data show that distinct genes code for the CaM proteins. Due to the current annotation status of the genome, we have no data on the relative chromosomal positions of these genes. Schistosoma mansoni has also been shown to possess two distinct CaM genes, although they exhibit high (>99 %) sequence identity relative to each other .
FhCaMs are abundantly expressed in parenchymal cells of juvenile liver fluke
FhCaM RNAi triggers specific suppression of target transcript and protein
We next analysed the longevity of CaM transcript knockdown in vitro over a 3-week period in order to begin to assess the time frame over which functional assays could be run. CaM knockdown persisted at significant levels during maintenance (measurements taken at 7, 14 and 21 days post dsRNA exposure; Fig. 3d–f), with transcript levels at day 21 remaining significantly lower than controls (FhCaM1 = 43.1 ± 4.1 %, n = 4, P < 0.01; FhCaM2 = 1.5 ± 1.5 %, n = 5, P < 0.001; FhCaM3 = 4.0 ± 1.3 %, n = 5, P < 0.001; Fig. 3f). These data showed that CaM transcript knockdown persists throughout this maintenance period; this persistence occurs even in the absence of dsRNA during maintenance, since worms were only exposed to dsRNA during the initial 4 h exposure soak. Such long-term maintenance of knockdown is consistent with our previous observations on RNAi of liver fluke cathepsin B, L and glutathione transferase genes, in which we reported the persistence of knockdown for 21 days following a similar RNAi protocol . Similarly, electro-soaking methods have also shown the persistence of silencing effects up to 21 days  and reports from schistosomes illustrate that electroporation-induced RNAi can be maintained for up to 40 days . Such persistent silencing in trematodes, even in the absence of prolonged dsRNA supplementation, indicates the presence of an efficient dsRNA uptake and amplification mechanism and bodes well for the success of putative in vivo RNAi experiments.
FhCaM RNAi induces motility and growth phenotypes in vitro
Assay selection was guided by our hypotheses that CaMs are involved in worm growth and motility. Previous work has linked CaM function with growth and development of helminths, including reduced growth associated with CaM silencing in both S. mansoni sporocysts , and in the free-living nematode C. elegans . Recently, the silencing of adult S. mansoni CaM has implicated the protein in muscular function and resulted in a characteristic somatic contraction/dilation phenotype . Given the persistence of CaM transcript knockdown (Fig. 3) and significant protein suppression from day 7 onwards (Fig. 4), we hypothesised that any associated phenotypes resulting from this suppression would be similarly detectable from day 7 onwards.
The rationale for investigating the impact of FhCaM RNAi on fluke motility stemmed partly from our localisation data (see Fig. 2), where we suspected that some of the stained cells in the parenchyma would include muscle myocytons, combined with the well-documented role of CaM in regulating smooth muscle contraction . If that was the case, we hypothesised that FhCaM RNAi should affect the ability of worms to move and/or migrate (although CaM suppression may also have been expected to affect motility via disruption of Ca2+-induced synaptic vesicle release, or through more general induction of cellular stress associated with widespread Ca2+ disruption). Certainly, data from other systems link disruption of CaM function with hindered movement/motility, particularly in ciliated or flagellated organisms/cells . Impacts of FhCaM RNAi on fluke motility were quantified based on the ability of worms to migrate through an agar substrate, across a 5 mm radius from the point of origin. Performed at 10 days post dsRNA, this assay illustrated that combinatorially-treated dsFhCaM1-3 (P < 0.01), but not individual dsFhCaM treatments, migrated significantly faster than untreated controls through this assay (proportion of worms escaping 5 mm radius: untreated, 38.6 ± 5.2 %, n = 13; dsCTRL, 45.2 ± 5.3 %, n = 13; dsFhCaM1, 33.6 ± 7.1 %, n = 7; dsFhCaM2, 35.3 ± 5.1 %, n = 13; dsFhCaM3, 50.2 ± 6.1 %, n = 13; dsFhCaM1-3, 64.3 ± 5.9 %, n = 12, P < 0.01; Fig. 5), i.e. worms exposed to a cocktail of three FhCaM dsRNAs (dsFhCaM1-3) were hyperactive and migrated more quickly over a 3 h period than untreated controls. None of the individual dsFhCaM treatments induced motor phenotypes that were significantly different from controls, suggesting that this phenotype may be the result of removing functional redundancy that may exist between the three fluke CaM genes.
FhCaM RNAi phenotypes are recapitulated by the CaM antagonist, trifluoperazine
In order to further validate our RNAi-induced growth and motility phenotypes, we incubated juvenile worms in the presence of TFP, a compound with CaM antagonising activity. This compound (and another calmodulin antagonist, W7) has been shown to bind to FhCaM1 and FhCaM3 . Additionally, CaM antagonists have been shown to partially block transformation in S. mansoni miracidia  and inhibit hatching of schistosome eggs . Growth inhibition associated with F. hepatica CaM-RNAi was also seen in worms maintained in the presence of the CaM antagonist, TFP, consistent with the effects of this drug as reported in schistosome sporocysts . Juvenile fluke were assayed after 7 days in the presence of TFP, with those treated with 0.5 and 1 μM showing significant reductions in size (P < 0.0001) compared to vehicle controls (untreated (no drug), 28,702 ± 890 μm2, n = 38; vehicle (DMSO) control, 26,775 ± 580 μm2, n = 50; 0.5 μM TFP, 20,276 ± 801 μm2, n = 50; 1 μM TFP, 18,776 ± 737 μm2, n = 60; Fig. 5).
As with the inhibited growth phenotype discussed above, this accelerated motility phenotype was recapitulated in the presence of the CaM-inhibiting drug, TFP. CaM RNAi has not, to our knowledge, been linked to increased motility in any other system to date, although TFP has been reported to increase spontaneous muscle contractions in the filarial nematode Acanthocheilonema vitae . This increased motility phenotype might occur due to the disruption of intracellular Ca2+ dynamics in neuromuscular cells, either in muscle cells, by triggering aberrant contractility, or in neurones, by impacting on Ca2+-evoked synaptic release. We do not yet know how increased motility would impact on worm infectivity under in vivo conditions, although these experiments are currently in progress. While increased motility may not make an obvious therapeutic strategy, it is possible that the increased motility triggered by CaM RNAi might impact worm survival through increased energy expenditure. Alternatively, stunted growth and development as observed when FhCaM function was disrupted by both molecular genetic and pharmacological means, does represent a desirable therapeutic outcome. This might reflect either: (i) FhCaM’s direct involvement in controlling Ca2+ dynamics in dividing cells (Ca2+ fluxes have a well-established role in the cell cycle; see [45, 46]) where disrupting CaM function could directly inhibit cell division; or (ii) FhCaMs are important housekeeping proteins required for normal cellular function , such that FhCaM’s impact on growth occurs via a more general impact on cellular health/function in FhCaM-expressing cells. Although the data presented in this study cannot further delineate between these possibilities, growth suppression could represent an appealing therapeutic outcome following FhCaM-selective drug intervention.
This work represents only the second published study to use gene-silencing methods to identify aberrant phenotypes that inform gene function in liver fluke parasites. We have demonstrated that CaMs are expressed widely in the parenchyma of the invasive juvenile life stage of F. hepatica. The development and application of novel F. hepatica assays demonstrated that CaMs are important for the normal growth and motility of NEJs. Additionally, the growth phenotype observed was only evident upon the addition of serum-stimulated development and we recommend the addition of serum to any further RNAi experiments in F. hepatica. It is worth noting that juvenile fluke growth in vitro is slower than that seen in vivo such that the growth phenotype associated with FhCaM dysregulation could be markedly enhanced in vivo. While the increased motility effect does not implicate FhCaM as an obvious therapeutic target, the growth defects observed encourage further investigation. Future work will focus on deciphering functional differences in the roles of FhCaM1, FhCaM2 and FhCaM3, particularly during development in vivo. The role of FhCaMs in growth supports their consideration as control targets in liver fluke parasites.
The authors acknowledge the following: BBSRC research grants BB/H009477/1 and BB/K009583/1; a Department of Agriculture and Rural Development PhD studentship award to EMC; and, the Glover family for the Professor John Glover Memorial Studentship award to PMC. We also thank Norman Baldwin for providing a continuous supply of liver fluke metacercariae.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Keiser J, Utzinger J. Emerging foodborne trematodiasis. Emerging Infect Dis. 2005;11(10):1507–14.PubMed CentralView ArticlePubMedGoogle Scholar
- Hotez PJ, Brindley PJ, Bethony JM, King CH, Pearce EJ, Jacobson J. Helminth infections: the great neglected tropical diseases. J Clin Invest. 2008;118(4):1311–21.PubMed CentralView ArticlePubMedGoogle Scholar
- Fairweather I. Reducing the future threat from (liver) fluke: realistic prospect or quixotic fantasy? Vet Parasitol. 2011;180(1–2):133–43.View ArticlePubMedGoogle Scholar
- Brennan GP, Fairweather I, Trudgett A, Hoey E, McCoy E, McConville M, et al. Understanding triclabendazole resistance. Exp Mol Pathol. 2007;82:104–9.View ArticlePubMedGoogle Scholar
- Winkelhagen AJ, Mank T, de Vries PJ, Soetekouw R. Apparent triclabendazole-resistant human Fasciola hepatica infection, the Netherlands. Emerging Infect Dis. 2012;18(6):1028–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Fairweather I. Triclabendazole: new skills to unravel an old(ish) enigma. J Helminthol. 2005;79(3):227–34.View ArticlePubMedGoogle Scholar
- Bootman MD, Collins TJ, Peppiatt CM, Prothero LS, MacKenzie L, De Smet P, et al. Calcium signalling--an overview. Semin Cell Dev Biol. 2001;12(1):3–10.View ArticlePubMedGoogle Scholar
- Berridge MJ, Bootman MD, Roderick HL. Calcium signalling: dynamics, homeostasis and remodelling. Nat Rev Mol Cell Biol. 2003;4(7):517–29.View ArticlePubMedGoogle Scholar
- Saimi Y, Kung C. Calmodulin as an ion channel subunit. Annu Rev Physiol. 2002;64:289–311.View ArticlePubMedGoogle Scholar
- O’Day DH. CaMBOT: profiling and characterizing calmodulin-binding proteins. Cell Signal. 2003;15(4):347–54.View ArticlePubMedGoogle Scholar
- Hoeflich KP, Ikura M. Calmodulin in action: diversity in target recognition and activation mechanisms. Cell. 2002;108(6):739–42.View ArticlePubMedGoogle Scholar
- Day TA, Bennett JL, Pax RA. Praziquantel: The enigmatic antiparasitic. Parasitol Today. 1992;8(10):342–4.View ArticlePubMedGoogle Scholar
- Taft AS, Yoshino TP. Cloning and functional characterization of two calmodulin genes during larval development in the parasitic flatworm Schistosoma mansoni. J Parasitol. 2011;97(1):72–81.PubMed CentralView ArticlePubMedGoogle Scholar
- Guidi A, Mansour NR, Paveley RA, Carruthers IM, Besnard J, Hopkins AL, et al. Application of RNAi to genomic drug target validation in schistosomes. PLoS Negl Trop Dis. 2015;9(5):e0003801.PubMed CentralView ArticlePubMedGoogle Scholar
- Katsumata T, Kohno S, Yamaguchi K, Hara K, Aoki Y. Hatching of Schistosoma mansoni eggs is a Ca2+/calmodulin-dependent process. Parasitol Res. 1989;76(1):90–1.View ArticlePubMedGoogle Scholar
- Kawamoto F, Shozawa A, Kumada N, Kojima K. Possible roles of cAMP and Ca2+ in the regulation of miracidial transformation in Schistosoma mansoni. Parasitol Res. 1989;75(5):368–74.View ArticlePubMedGoogle Scholar
- Zhou J, Sun J, Huang Y, Zhou C, Liang P, Zheng M, et al. Molecular identification, immunolocalization, and characterization of Clonorchis sinensis calmodulin. Parasitol Res. 2013;112(4):1709–17.View ArticlePubMedGoogle Scholar
- Hu S, Law PK, Lv Z, Wu Z, Fung MC. Molecular characterization of a calcium-binding protein SjCa8 from Schistosoma japonicum. Parasitol Res. 2008;103(5):1047–53.View ArticlePubMedGoogle Scholar
- Russell SL, McFerran NV, Hoey EM, Trudgett A, Timson DJ. Characterisation of two calmodulin-like proteins from the liver fluke, Fasciola hepatica. Biol Chem. 2007;388(6):593–9.View ArticlePubMedGoogle Scholar
- Russell SL, McFerran NV, Moore CM, Tsang Y, Glass P, Hoey EM, et al. A novel calmodulin-like protein from the liver fluke, Fasciola hepatica. Biochimie. 2012;94(11):2398–406.View ArticlePubMedGoogle Scholar
- McVeigh P, McCammick EM, McCusker P, Morphew RM, Mousley A, Abidi A, et al. RNAi dynamics in juvenile Fasciola spp. liver flukes reveals the persistence of gene silencing. PLoS Negl Trop Dis. 2014;8(9):e3185.PubMed CentralView ArticlePubMedGoogle Scholar
- McGonigle L, Mousley A, Marks NJ, Brennan GP, Dalton JP, Sprithill TW, et al. The silencing of cysteine proteases in Fasciola hepatica newly excysted juveniles using RNA interference reduces gut penetration. Int J Parasitol. 2008;38:149–55.View ArticlePubMedGoogle Scholar
- Cwiklinski K, Dalton JP, Dufresne PJ, La Course J, Williams DJL, Hodgkinson J, et al. The Fasciola hepatica genome: gene duplication and polymorphism reveals adaption to the host environment and the capacity for rapid evolution. Genome Biol. 2015;16(1):71.PubMed CentralView ArticlePubMedGoogle Scholar
- Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, et al. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16(10):944–5.View ArticlePubMedGoogle Scholar
- GATC Biotech. www.gatc-biotech.com. Accessed April 2013.
- Integrated DNA technologies. www.idtdna.com. Accessed May 2013.
- Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9):e45.PubMed CentralView ArticlePubMedGoogle Scholar
- Real Time PCR Miner software. http://miner.ewindup.info. Accessed May-August 2013.
- Zhao S, Fernald RD. Comprehensive algorithm for quantitative real-time polymerase chain reaction. J Comput Biol. 2005;12(8):1047–64.PubMed CentralView ArticlePubMedGoogle Scholar
- Pierson L, Mousley A, Devine L, Marks NJ, Day TA, Maule AG. RNA interference in a cestode reveals specific silencing of selected highly expressed gene transcripts. Int J Parasitol. 2010;40:605–15.View ArticlePubMedGoogle Scholar
- Bar-Sagi D, Prives J. Trifluoperazine, a calmodulin antagonist, inhibits muscle cell fusion. J Cell Biol. 1983;97(5 Pt 1):1375–80.View ArticlePubMedGoogle Scholar
- Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–5.View ArticlePubMedGoogle Scholar
- Karabinos A, Büssing I, Schulze E, Wang J, Weber J, Schnabel R. Functional analysis of the single calmodulin gene in the nematode Caenorhabditis elegans by RNA interference and 4-D microscopy. Eur J Cell Biol. 2003;82(11):557–63.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- McVeigh P, Mair GR, Novozhilova E, Day A, Zamanian M, Marks NJ, et al. Schistosome I/Lamides--a new family of bioactive helminth neuropeptides. Int J Parasitol. 2011;41(8):905–13.PubMed CentralView ArticlePubMedGoogle Scholar
- Davies C, Smyth JD. In vitro cultivation of Fasciola hepatica metacercariae and of partially developed flukes recovered from mice. Int J Parasitol. 1978;8(2):125–31.View ArticlePubMedGoogle Scholar
- Threadgold LT, Gallagher SS. Electron microscope studies of Fasciola hepatica. I. The ultrastructure and interrelationship of the parenchymal cells. Parasitology. 1966;56(2):299–304.View ArticlePubMedGoogle Scholar
- Gallagher SS, Threadgold LT. Electron-microscope studies of Fasciola hepatica. II. The interrelationship of the parenchyma with other organ systems. Parasitology. 1967;57(4):627–32.View ArticleGoogle Scholar
- Pax RA, Day TA, Miller CL, Bennett JL. Neuromuscular physiology and pharmacology of parasitic flatworms. Parasitology. 1996;113(Suppl):S83–96.View ArticlePubMedGoogle Scholar
- Dell'Oca N, Basika T, Corvo I, Castillo E, Brindley PJ, Rinaldi G, et al. RNA interference in Fasciola hepatica newly excysted juveniles: long dsRNA induces more persistent silencing than siRNA. Mol Biochem Parasitol. 2014;197(1–2):28–35.View ArticlePubMedGoogle Scholar
- Krautz-Peterson G, Radwanska M, Ndegwa D, Shoemaker CB, Skelly PJ. Optimizing gene suppression in schistosomes using RNA interference. Mol Biochem Parasitol. 2007;153:194–202.View ArticlePubMedGoogle Scholar
- Walsh MP. Calmodulin and the regulation of smooth muscle contraction. Mol Cell Biochem. 1994;135(1):21–41.View ArticlePubMedGoogle Scholar
- Ashizawa K, Tomonaga H, Tsuzuki Y. Regulation of flagellar motility of fowl spermatozoa: evidence for the involvement of intracellular free Ca2+ and calmodulin. J Reprod Fertil. 1994;101(2):265–72.View ArticlePubMedGoogle Scholar
- Minardi AJ, Christ D, Saz HJ. Effects of calmodulin and protein kinase C antagonists on muscle in the filariid, Acanthocheilonema viteae. J Parasitol. 1995;81(6):989–96.View ArticlePubMedGoogle Scholar
- Poenie M, Alderton J, Tsien RY, Steinhardt RA. Changes of free calcium levels with stages of the cell division cycle. Nature. 1985;315(6015):147–9.View ArticlePubMedGoogle Scholar
- Hepler P. The role of calcium in cell division. Cell Calcium. 1994;16(4):322–30.View ArticlePubMedGoogle Scholar
- Berchtold MW, Villalobo A. The many faces of calmodulin in cell proliferation, programmed cell death, autophagy, and cancer. Biochim Biophys Acta. 2014;1843(2):398–435.View ArticlePubMedGoogle Scholar