Novel cholinesterase paralogs of Schistosoma mansoni have perceived roles in cholinergic signalling and drug detoxification and are essential for parasite survival

Cholinesterase (ChE) function in schistosomes is essential for orchestration of parasite neurotransmission but has been poorly defined with respect to the molecules responsible. Interrogation of the S. mansoni genome has revealed the presence of three ChE domain-containing genes (Smche)s, which we have shown to encode two functional acetylcholinesterases (AChE)s (Smache1 –smp_154600 and Smache2 –smp_136690) and a butyrylcholinesterase (BChE) (Smbche1 –smp_125350). Antibodies to recombinant forms of each SmChE localized the proteins to the tegument of adults and schistosomula and developmental expression profiling differed among the three molecules, suggestive of functions extending beyond traditional cholinergic signaling. For the first time in schistosomes, we identified ChE enzymatic activity in fluke excretory/secretory (ES) products and, using proteomic approaches, attributed this activity to the presence of SmAChE1 and SmBChE1. Parasite survival in vitro and in vivo was significantly impaired by silencing of each smche, either individually or in combination, attesting to the essential roles of these molecules. Lastly, in the first characterization study of a BChE from helminths, evidence is provided that SmBChE1 may act as a bio-scavenger of AChE inhibitors as the addition of recombinant SmBChE1 to parasite cultures mitigated the effect of the anti-schistosome AChE inhibitor 2,2- dichlorovinyl dimethyl phosphate—dichlorvos (DDVP), whereas smbche1-silenced parasites displayed increased sensitivity to DDVP.


Introduction
The functioning of the nervous system is a tightly regulated process controlled through multiple catalytic and non-catalytic signaling proteins. Among the catalytic molecules, cholinesterases (ChEs) play a pivotal role in regulating the signaling activity of the nervous system. There are two major types of ChEs, acetylcholinesterase (AChE) and pseudocholinesterase, or butyrylcholinesterase (BChE), and they can be distinguished both kinetically and pharmacologically [1]. AChE selectively hydrolyzes the neurotransmitter acetylcholine (ACh) to maintain neurotransmitter homeostasis [2] while the main role of BChE is widely accepted to be the detoxification of organophosphorus esters which are inhibitors of AChE [3]. ChEs are generally believed to be functionally redundant in cholinergic signaling with the main differences between paralogs lying in their spatial and temporal expression as well as non-cholinergic functionality [4,5].
The nervous system of helminths has long been a potential target for therapeutic agents as it plays several crucial roles in parasite biology that are fundamental to survival, including coordinating motility within and outside of the host, feeding and reproduction [6][7][8][9][10]. The Schistosoma nervous system is particularly important in this respect as this parasite lacks a body cavity and circulating body fluid [11,12] and, as a result, its signaling functions are chiefly achieved via neurotransmission. The primary neurotransmitter that schistosomes utilize is acetylcholine (ACh), which allows muscle contraction. The physiological concentration of ACh, however, must be maintained otherwise it triggers paralysis and this is achieved primarily through the action of AChE [6][7][8].
While AChE activity has been documented extensively in S. mansoni (reviewed in [13]), most of the work has involved studies on parasite extracts or native SmChE purified by inhibitor-affinity chromatography, making it difficult to attribute function to any one particular SmChE molecule. More recent studies have documented two AChE-encoding genes in S. mansoni and other species [14][15][16]. In 2016, You et al. characterized AChE activity in S. japonicum extracts and at a molecular level, but only through the expression of one recombinant AChE [17]. Moreover, to the best of our knowledge, genes encoding proteins with BChE activity have not been previously described in schistosomes or any other helminth. Interrogation of the now fully annotated S. mansoni genome [18] has revealed three different SmChE paralogs; however, their individual contributions to ChE function remain unknown.
Traditionally, ChEs have been regarded solely as neurotransmitter terminators; however, there is increasing data to suggest that these enzymes play a variety of roles that extend beyond this cholinergic function due to their presence in multiple cell types and subcellular locations [4,5,19,20]. In schistosomes, AChE has been localized to the tegument as well as the neuromusculature [17,21,22] and a proposed function for tegumental AChE has been mediation of glucose uptake by the parasite from the external environment [23]. The exact mechanism for this process is still unclear but the proposed initiating step is by limiting the interaction of host ACh with tegumental nicotinic ACh receptors (nAChRs), a hypothesis bolstered by the observation that glucose uptake is ablated by the use of membrane-impermeable AChE and nAChR inhibitors in S. mansoni [23,24] and RNAi-mediated AChE silencing in S. japonicum [25]. The nAChRs are also associated both spatially and temporally with surface AChE expression and are concentrated on the tegument [26], the major site of glucose uptake [27].
Many intestinal nematodes secrete AChE [28][29][30][31], which, where studied, orchestrate exogenous cholinergic activities. It has also been indirectly shown that the nematode N. brasiliensis employs parasite-derived AChE to alter the host cytokine environment to inhibit M2 macrophage recruitment, a condition favorable to worm survival [32]. Despite this breadth of literature in nematodes, there has been no documentation of secreted AChE activity from schistosomes.
Herein we describe and functionally characterize using gene silencing and enzymatic approaches, a novel AChE and BChE from S. mansoni and further characterize the only previously identified AChE-encoding gene from the parasite. Importantly, we show through gene knockdown that each smche is essential to S. mansoni development and survival, highlighting them as targets for novel anti-schistosomal intervention strategies.

Identification of novel genes encoding ChE proteins in S. mansoni
Three putative ChE paralogs were identified from interrogation of the S. mansoni genome: smache1 (Smp_154600), smbche1 (Smp_125350) and smache2 (Smp_136690). The predicted SmChEs were then aligned with characterized AChE enzymes from Homo sapiens, the electric eel Torpedo californica, and the nematode Caenorhabditis elegans (Fig 1). Homology analysis of amino acid sequences revealed that SmAChE1, SmBChE1, and SmAChE2 share (32-35%) sequence identity and (49-52%) sequence similarity. Further, all SmChEs have 36-40% amino acid identity with H. sapiens and T. californica AChE. All identified SmChEs had ChE-specific characteristics, including a catalytic triad with an active site serine, which is required for ester hydrolysis [33]. Interestingly, the His residue of the catalytic triad of SmBChE1 appears to have been substituted for Gln, a change consistent among all the BChE1 homologs shown for other Platyhelminthes, but not nematode or model organism BChE1 sequences (S1 Fig). A 3D model of the three SmChEs was constructed by homology modeling with AChE from model organisms (H. sapiens and T. californica (S2 Fig). All three SmChEs exhibited predicted folding characteristics of the functional globular enzymes as most of the α-helical and β-stranded sheets were tightly aligned. Each predicted SmChE structure consisted of a ChE catalytic domain but, although the core architecture of the catalytic gorges was well aligned, regions that are associated with substrate specificity and catalytic efficiency were disparate. In particular, and in agreement with the sequence alignment, the catalytic triad of SmBChE1 was predicted to be Ser-Gln-Glu instead of the canonical Ser-His-Glu present in the other two paralogs. A phylogenetic tree of the alignment (S3 Fig) shows that SmChEs were clustered into three distinct branches, with SmChE1 being phylogenetically distinct from SmBChE1 and SmAChE2. In addition, each SmChE was grouped together with closely related flatworms, including other Schistosoma species. Importantly, as shown in the sequence alignment, SmChEs are divergent from the human homolog. Reflective of the catalytic triad residue difference (S1 expression patterns were variable, the transcript levels of all three genes were relatively lower in cercariae compared to the other developmental stages. Overall, the transcript levels of smache1 and smache2 genes in most life stages were higher than that of smbche1. In adult worms, smache1 was expressed at higher levels, specifically in male parasites, followed by sporocysts.

Immunolocalization of SmChEs
To gain insight into the anatomical sites of expression of ChE proteins in S. mansoni, SmChEs were immunolocalized in whole juvenile and sectioned adult parasites. In adults, and consistent with their predicted cholinergic function, all SmChEs were expressed throughout the worms' internal structures (presumably localizing to the neuromusculature) and on the tegument surface. SmAChE1 was the least uniformly distributed of all SmChEs, localizing mostly to the tegument ( Fig 3A). Additionally, anti-SmChE antibodies were able to detect homologous ChEs in adult S. haematobium sections. SmChE proteins were detected in all stages of larval development tested and, as was the case with adult worms, localized to the tegument (Fig 3B).

Expression and ChE activity of SmChEs
Soluble, functionally active proteins were expressed in P. pastoris, purified via IMAC and tested for ChE activity. Both SmAChE1 and SmAChE2 demonstrated significantly stronger hydrolase activity when AcSCh was used as a substrate, compared to SmBChE1 and, conversely, SmBChE1 hydrolyzed BcSCh to significantly higher levels compared to SmAChE1 and SmAChE2 ( Fig 4A). All paralogs exhibited Michaelis-Menten kinetics (Table 1) when hydrolyzing their designated substrate, with SmAChE1 having a substrate affinity approximately twice that of SmAChE2. In addition, preferred substrate activity of both SmAChE1 and SmAChE2 was inhibited by DDVP, an AChE inhibitor, whereas iso-OMPA, a specific inhibitor of BChE, only inhibited SmBChE1 activity (Fig 4B).

BChE and secreted ChE activity in schistosomes
Although the presence of nonspecific ChE activity has long been known in schistosomes [34], the identity of the gene product and its function remain unknown. Prompted by the identification of SmBChE1 as a BChE, based on its substrate preference and enzymatic inhibition by iso-OMPA, we sought to investigate the distribution of BChE activity in juvenile and adult schistosomes. Extracts from S. mansoni schistosomula had higher BChE activity compared to S. mansoni adult worms (Fig 5A), and that activity was significantly greater in S. mansoni compared with S. haematobium adults (Fig 5B). Varied amounts of AChE activity were detected in ES from all developmental stages tested. ES products from adult males had double the AChE activity of adult female ES products (P < 0.001), while cercariae ES exhibited the highest activity (at least ten-fold more than male ES products (P < 0.0001)) and egg ES had the lowest ( Fig  5C). Availability of ES precluded the measurement of secreted BChE activity from all developmental stages but, of those tested, activity in schistosomula ES products was the highest-twice as high as that of adult (P < 0.01) and cercariae (P < 0.01) ES ( Fig 5D). SmChEs were purified from ES products of S. mansoni adult worms using edrophonium-sepharose affinity chromatography. Purification resulted in an activity increase of more than 200-fold relative to crude ES ( Fig 5E). Resolution of the purified proteins by SDS-PAGE resulted in a doublet with a major band migrating at 70 kDa under denaturing and reducing conditions (Fig 5F). The identity of purified, secreted SmChEs was substantiated by in-gel LC-MS/MS analysis with the peptide data generated used to interrogate the S. mansoni proteome (predicted from the S. mansoni genome -http://www.genedb.org/Homepage/Smansoni)). The false discovery rate was set at <1% and only proteins with at least two unique peptides having significant Mascot identification scores (P < 0.05) were considered. The top protein hits were identified as SmAChE1 (Smp_154600) and SmBChE1 (Smp_125350); SmAChE1 had a relative abundance of more than 40-fold that of SmBChE1 (S3 Table).

RNAi-mediated smche transcript and SmChE protein reduction
Schistosomula electroporated with smache1 siRNA showed respective decreases in SmAChE1 mRNA levels of 43% (P � 0.05) and 84% (P � 0.001) at three and seven days post-treatment, respectively, compared to the luc control (Fig 6A), while treating parasites with smbche1 siRNA caused 51.0% (P � 0.001) and 78% (P � 0.001) suppression of smbche1 mRNA expression at day 3 and 7 after electroporation, respectively, compared to the luc control ( Fig 6B). Treatment of schistosomula with smache2 siRNA resulted in respective decreases in smache2  mRNA levels of 60% (P � 0.001) and 47% (P � 0.01) five and seven days after electroporation, compared to the luc control ( Fig 6C). Schistosomula electroporated with a cocktail of all three smche siRNAs showed decreases of all three transcript levels over time, with smache2 mRNA levels decreasing by an average of 90% (P � 0.001) by day 3 after treatment, compared to the luc control ( Fig 6D).

Suppression of SmChE activity
Suppression of AChE activity was seen in smache1 and smache2 siRNA-treated parasites from 5 and 3 days after electroporation, respectively (Fig 7A), compared to the luc control, while schistosomula treated with smbche1 siRNA did not show any significant reduction in BChE activity, even 7 days after electroporation ( Fig 7B). Parasites electroporated with a cocktail of all three smche siRNAs showed significant decreases in AChE activity at 3 days (62% reduction, P � 0.001), 5 days (67% reduction, P � 0.001) and 7 days (71% reduction, P � 0.001) after treatment ( Fig 7C). BChE activity was not measured in the cocktail siRNA treatment group.
Individual silencing of smache1 or smache2 genes and combined silencing of all three smche genes reduced glucose uptake in schistosomula by 24.9% (P � 0.001), 32.34% (P � 0.001) and 38.61% (P � 0.001) at 48 h post-treatment, respectively, relative to the luc control. However, smbche1-silenced parasites showed no significant changes in glucose uptake at the same timepoint and there was no difference in the glucose consumed by the smache1 or smache2 siRNA-treated groups compared with the cocktail siRNA-treated group ( Fig 7D). Transcript levels of the glucose transporters, sgtp1 and sgtp4, were neither decreased nor significantly increased in individual or cocktail smche-silenced parasites (S5 Fig).

Effects of smche silencing on schistosomula viability in vitro and development in vivo
Parasites treated with smache1, smbche1 or smache2 siRNAs showed significant decreases in viability at days 3, 5 and 1 after treatment, respectively, compared to luc controls. At days 5 and 7 post-treatment, the most significant decrease in parasite viability was seen in the group which received the cocktail siRNA treatment, compared to luc controls. Furthermore, viability in this group was also significantly lower than it was for any individual treatment at these two time points (Fig 8A).
To examine whether RNAi-mediated smche suppression reduced parasite viability in vivo, mice were infected with smche-silenced parasites and worm burdens were measured after three weeks. From two independent experiments, there was an average 88.15%, 55.15%, 75.95% and 88.35% decrease in adult fluke burdens from mice injected with smache1-, smbche1-, smache2-and smche cocktail-silenced schistosomula, respectively, compared to mice infected with luc-treated parasites (Fig 8B and 8C). All worm burden decreases were significant and there was no significant difference in fluke burdens between mice injected with luc-treated parasites and non-electroporated control parasites. All mice had been successfully infected with parasites, as serum from necropsied mice contained parasite-specific antibodies (S6A Fig

Bio-scavenging of carboxylic esters by SmBChE1
The hypothesis that SmBChE may act as a molecular decoy in schistosomes to detoxify the effects of organophosphorus AChE inhibitors was examined by testing whether, firstly, addition of DDVP (an organophosphorous AChE inhibitor) to parasite extracts or recombinant SmBChE1 inhibited BChE activity, secondly, inhibition of parasite-derived BChE potentiated the effects of DDVP and, thirdly, addition of exogenous BChE (SmBChE1) mitigated the effects of DDVP. Firstly, BChE activity of schistosome extracts and recombinant SmBChE1 significantly decreased in the presence of increasing concentrations of DDVP (Fig 9A and 9B, respectively). Secondly, DDVP-mediated killing of schistosomula was significantly increased in the presence of iso-OMPA (60.9% compared with 21.8%; P < 0.0001) ( Fig 9C). AChE activity in live parasites was not affected by the addition of iso-OMPA (S7A and S7B Fig). Further, smbche1-silenced schistosomula were significantly more susceptible to DDVP-mediated killing than luc-treated controls (67.8% compared with 27.5%; P < 0.0001) ( Fig 9D). The iso-OMPA addition experiment was attempted using extracts and, while AChE activity in DDVPtreated schistosome extracts significantly decreased in the presence of increasing amounts of the BChE inhibitor, iso-OMPA, it also decreased, albeit to a lesser degree, when iso-OMPA was added to extracts without DDVP (S8A Fig). Thirdly, schistosomula were increasingly resistant to DDVP-mediated killing with the addition of increasing amounts of recombinant protein to the culture media ( Fig 9E). The SmBChE1 addition experiment was attempted using extracts and, while AChE activity in DDVP-treated extracts was significantly increased

Discussion
Cholinesterase (ChE) activity in S. mansoni was first described by Bueding in 1952 [34] and was well characterized biochemically in the four decades succeeding this discovery. The technological limitations of this time period meant that most of the evidence for SmChEs came from whole worm studies and analyses of crude parasite extracts (reviewed in [13]), which could not ascribe ChE activity to any particular protein. Several studies in the early 2000's characterized a single AChE from S. mansoni (Smp_154600 in the current gene annotation nomenclature) and its direct homolog in other species of schistosomes [14,15,22], but lack of smche siRNAs or luc siRNA were cultured for 7 days in complete Basch medium with viability determined at day 1, 3, 5 and 7 after treatment by Trypan Blue exclusion (mean ± SEM of duplicate assays from 2 biological replicates of each treatment). (B and C) One-day-old schistosomula treated with individual or a cocktail of all three smche siRNAs or luc siRNA were intramuscularly injected (2,000 parasites) into mice. After 3 weeks, adult worms were recovered and counted. Data from two independent experiments are shown. Differences between smche-and luc-treated groups were measured by the student's t test. � P � 0.05, �� P � 0.01, ��� P � 0.001.
https://doi.org/10.1371/journal.ppat.1008213.g008 a comprehensive schistosome genome annotation at the time precluded identification of more ChE family members. Interrogation of the most recent iteration of the S. mansoni genome assembly has identified two additional ChE-encoding genes that are paralogs to Smp_154600 (which we have termed SmAChE1); Smp_125350 (SmBChE1) and Smp_136690 (SmAChE2).
In this current study, we have provided a more in-depth characterization of the previously documented SmAChE1 and two additional ChEs from S. mansoni: SmAChE2 -an AChE which has recently been characterized in terms of its gene expression in adult S. mansoni [16]; and SmBChE1 -a BChE which, to the best of our knowledge, has never been documented in the helminth literature. All SmChEs share a modest level of identity which is consistent with their divergence over evolutionary time, an occurrence that is possibly due to a series of gene duplications due to the phylogenetic distance between the relative clades. This divergence between SmChEs and, also, ChEs of other organisms, provides evidence for the increasing reports of non-cholinergic functions of ChEs in the literature. Additionally, the relative lack of sequence identity between SmChEs and human ChEs suggests potential scope for the development of intervention strategies targeting schistosome ChEs that will not affect the host. Despite the diversity between ChEs, all enzymes analyzed herein would appear to be enzymatically active as they possessed a catalytic triad with an active site serine, the amino acid responsible for ester hydrolysis [33]. It is interesting to note, however, the catalytic triad His-Gln substitution in SmBChE1 (and the other platyhelminth BChE1 homologs); while this change is not a hallmark of model BChEs, that it occurs within an entire parasite lineage is noteworthy and will be investigated further.
The transcript levels of each smche varied among parasite developmental stages and this is likely a response to the differing cholinergic and cholinesterase-independent requirements of the parasite throughout its lifecycle. For example smache1 is expressed at a higher level in adult males than females, probably due to the more "muscular" roles of attachment and movement orchestrated by the male compared to the female, which remains sedentary once inside the gynecophoric canal of the male [35]. This hypothesis is supported by the additional observation that smache1 is downregulated in females after pairing [16]. Expression of smbche1 was highest in the egg stage; there is evidence for BChE involvement in chicken embryo neurogenesis and development, independent of its enzymatic function [36], which suggests that SmBChE1 could play a role in parasite embryogenesis. The miracidium and sporocyst stages had the highest levels of smache2 expression, in agreement with Parker-Manuel et al [37].
Immunolocalization of the SmChEs revealed apparent expression in the neuromusculature and tegument to varying degrees, depending on the paralog, and is consistent with early localization experiments [21], although the antibodies used in those studies were raised against AChEs purified from parasite extracts and so the localization could not be attributed to a specific family member. Tegumental distribution is suggestive of non-neuronal cholinergic and/ or non-cholinergic roles. Indeed, surface-expressed SmAChE has been implicated in mediating glucose scavenging by the parasite, as this process can be ablated by membrane-impermeable AChE inhibitors [23,24]. Tegumental SmAChE may also act to hydrolyze exogenous ACh, neutralizing its immune-mediating function to create an environment more conducive to parasite establishment [32]. The localization of these proteins to the tegument of schistosomula should also be noted since early developing schistosomula are considered most vulnerable to immune attack [38], and so SmAChE-targeted immunotherapeutics could be used effectively to vaccinate against schistosomiasis. Indeed, antibodies against SmAChEs have been shown to interact with the surface of schistosomula, resulting in complement-dependent killing of the parasite [39].
Full-length and functional SmChEs were expressed in P. pastoris. SmAChE1 had preferred substrate specificity for AcSCh over BcSCh, albeit at a three-fold lower affinity than previously reported for SmAChE1 expressed in Xenopus laevis oocytes [15]. SmAChE2 also had a substrate preference for AcSCh and an affinity twice that of SmAChE1. Extremely low enzyme activity was observed with SmBChE1 when AcSCh was used as a substrate, but enzymatic activity significantly increased with the use of BcSCh as the substrate. Although sequence alignment of SmBChE1 with the other two SmChEs revealed a single amino acid substitution in the peripheral anionic site (Glu-Trp), acyl binding pocket (Val-Leu) and catalytic triad (His-Gln), it was unclear whether these changes alone were enough to classify SmBChE1 as a BChE; based on the significant difference in substrate preference, however, this classification would appear valid. Cloning of a recombinant BChE from S. mansoni is consistent with our observations of BChE activity in parasite extracts, and S. mansoni schistosomula exhibited significantly more activity than adults, as did S. mansoni compared to S. haematobium adults. It has been reported that S. mansoni is more sensitive to the BChE inhibitor, iso-OMPA, than S. haematobium [40] and it may be due to the increased BChE activity in S. mansoni. Indeed, this relationship has been documented between AChE and metrifonate (precursor of DDVP used in this study); S. haematobium is more sensitive to the inhibitor than S. mansoni because of the greater amount of AChE on the worm's surface [40].
For the first time, we document the presence of secreted SmChE activity in schistosomes and AChE activity was highest in cercarial ES products. Of the intra-mammalian stages tested, AChE activity was highest in schistosomula and adults and may be acting to bind and neutralize exogenous AChE inhibitors (thus protecting tegumental and somatic AChE) or hostderived ACh to mitigate the immunomodulatory effects of this molecule. Extending this hypothesis, ES products from cultured female worms had lower AChE activity than males and could be due to females having less of a requirement for this defensive mechanism as they reside in the relative shelter of the gynecophoric canal. BChE activity was present in the ES products of adults, schistosomula and cercariae and was significantly higher in the intra-mammalian larval stage than the other two stages. The SmChE molecules present in ES were isolated by purification on edrophonium (a reversible ChE inhibitor) sepharose and, consistent with the class of activity observed in ES, identified by mass spectrometry as SmAChE1 and SmBChE1; the former being forty-fold more abundant than the latter.
RNAi-mediated silencing of smache1 and smache2 in schistosomula showed decreases in AChE activity, consistent with reductions in transcript and protein expression levels. Moreover, inhibition of this biochemical activity was greater in schistosomula treated with the smche siRNA cocktail than parasites receiving any of the individual treatments, further evidence suggestive of simultaneous silencing of all smache paralogs. AChE activity inhibition in smache2-silenced parasites was more pronounced than in smache1-silenced parasites, which was inconsistent with protein level reductions and this may be due to the increased AChE activity reported for SmAChE1 (Vmax = 5.57 nmol/min/mg, Km = 5.83 mM) compared to SmAChE2 (Vmax = 5.59 nmol/min/mg, Km = 10.87 mM).
Previous studies have documented the involvement of SmChEs in the uptake of exogenous glucose by schistosomes through the ablation of the glucose uptake pathway by organophosphorus [40] and large molecule [24] AChE inhibitors, so we sought to identify the SmChE paralog(s) responsible for this mediation through the use of RNAi targeting smche genes. Individual gene knockdown of smache1 and smache2 suppressed glucose uptake in schistosomula, implying that both genes were involved in regulation of this mechanism. Glucose uptake in these experiments was normalized to the amount of viable parasites in each treatment to control for any observed loss of uptake due to parasite death. Tegumental AChE is speculated to mediate glucose uptake by limiting the interaction of ACh with tegumental nicotinic ACh receptors which is thought to decrease the amount of glucose uptake through surface glucose transporters. The fact that both molecules are localized to the tegument and can hydrolyze ACh therefore provides evidence for their role in this pathway. Silencing of smbche1 in schistosomula did not show any difference in glucose uptake and is probably reflective of the molecule's limited role in ACh hydrolysis. Transcript levels of sgtp1 and sgtp4 were not significantly changed in smache1-and smache2-silenced parasites, suggesting that SmAChEs may facilitate glucose uptake in a manner which does not directly involve glucose transporters. Indeed, at least in nematodes, AChEs have been proposed to be involved in altering the permeability of surrounding host cells, allowing nutrients (such as glucose) to leak into the parasite niche and be uptaken [41].
Individual smche silencing in schistosomula resulted in significant decreases in parasite viability at various timepoints after treatment, with smache2-silenced parasites showing the most rapid and significant decrease in viability after treatment. Of the smche paralogs studied, smche2 is the only one whose expression is significantly upregulated between S. mansoni cercariae and schistosomula [37], an observation consistent with qPCR data, and so silencing this relatively highly expressed gene may have the most profound effects of all smche silencing on parasite viability. The viability of parasites treated with all three smche siRNAs was significantly decreased compared to parasites treated with an individual siRNA, suggesting functional overlap exists between the paralogs. This redundancy has been documented in AChE-knockout mice where BChE has the ability to hydrolyze ACh in the absence of AChE [42,43]. Moreover, AChE deletion is found to be lethal in Drosophila only because there is no alternative BChE paralog to compensate for the lack of ACh hydrolysis [44,45]. Similar to the observations in this study, simultaneous knockdown of multiple ChE genes has been reported to have deleterious effects on their target organisms including the insects Plutella xylostella [46], Chilo suppressalis [47] and Tribolium castaneum [48] and the nematodes Nippostrongylus brasiliensis [49] and Caenorhabditis elegans [50]. Similarly, chemotherapy with "broad spectrum" ChE inhibitors has shown to be effective against a range of organisms, including pest insects [51], schistosomes [24,[52][53][54] and parasitic nematodes [55]. It is likely that the simultaneous silencing of the smche genes in this study has a profound effect on parasite viability due to the knockdown of cholinergic signaling, a process to which all paralogs contribute, as they have all been shown to hydrolyze ChE substrates. Also possible is that knockdown of these three genes might have resulted in the ablation of multiple other functions that have been suggested for these molecules [5,48,56]. Reflective of in vitro silencing experiments, worm recovery from mice infected with all groups of smche-silenced parasites was significantly less than controls, indicating that suppression of smache1, smbche1 or smache2 could inhibit schistosome establishment and/or development in the host. An immunomodulatory function that results in a host environment favorable to parasite survival has been suggested for AChE secreted by N. brasiliensis [32] and so it may be that, if schistosome ChEs have similar non-neuronal roles, impairment of ChE function in these parasites may lead to more efficient immune-mediated worm expulsion. There is now a growing body of evidence that AChEs play non-classical roles as adhesion molecules [4,57,58] due to these two protein families sharing significant domain homology and so SmChE-silenced worms may be unable to properly establish in their site of predilection due to impaired adherence to host vasculature. Schistosomula silenced for all three smche genes exhibited the highest mortality in vivo when worm recovery was averaged across the two independent trials, suggesting that, not only is simultaneous knockdown of SmChEs required to overcome any functional redundancy between the molecules, but that this treatment has the largest impact on parasite pathogenesis due to the inhibition of multiple biological functions collectively orchestrated by these proteins.
The smche transcript levels of silenced worms recovered from mice were no different from control parasites and it is likely that the surviving worms received less siRNA than those worms that died in the host. Alternatively, surviving worms might have recovered from the transient effects of RNAi, highlighting the advantage of targeted gene knockout techniques such as CRISPR/Cas9 which has recently been reported for the first time in schistosomes [59].
It is generally accepted that vertebrate BChE has a predominant role in the detoxification of ingested or inhaled drugs and poisons such as the AChE-inhibitory organophosphorous esters that constitute nerve agents and pesticides due to the binding of the enzyme to these molecules, resulting in blockade of their active site (reviewed in [3]). The addition of the organophophorous AChE inhibitor, DDVP, to schistosome extracts decreased extract BChE activity, suggesting a similar detoxification role exists for schistosome BChE as for the vertebrate enzyme. Moreover, inactivation of SmBChE1 in live schistosomula by the BChE inhibitor iso-OMPA or through RNAi-mediated silencing potentiated the parasite-killing effects of DDVP (presumably less BChE available for DDVP binding and inactivation) whereas addition of exogenous SmBChE1 mitigated the results (presumably more BChE available for DDVP binding and inactivation), further supporting this hypothesis. It should be noted, however, that when AChE activity was used as a readout for similar experiments performed on extracts, the data was less compelling and could be due to the fact that BChE has a role (albeit limited) in ACh hydrolysis and so modification of BChE availability could have impacted on AcSCh conversion, confounding assay results. Numerous plant species in the Solanaceae used for their nutritional value (potatoes, for example) and others employed in traditional medicine for their anthelmintic properties contain naturally-occurring AChE-inhibitory compounds [60,61], and so it may be that the evolution of this dietary behavior in schistosomiasis endemic populations has resulted in selective pressure on the parasite to produce these particular ChE molecules. The localization of SmBChE1 to the tegument and its presence in ES products may further support this hypothesis, as the enzyme would be spatially available to interact with toxins present in the host environment, thus safeguarding parasite AChE against AChE inhibitors. Moreover, BChE activity is higher in S. mansoni than S. haematobium, which is more sensitive to the effects of the organophosphorous AChE inhibitor metrifonate. It has been reported that this sensitivity is due to the larger amount of tegumental AChE present in S. haematobium [40], but it may also be due to the reduced amount of BChE available to detoxify the inhibitor, as a similar relationship has been reported in studies which use human BChE to counter organophosphate toxicity [62]. Plasma-derived human BChE is currently in a phase I clinical trial as a nerve agent detoxifier, and a recombinant human BChE mutant is being used to prevent relapse in cocaine addicts due to the enzyme's ability to hydrolyze the drug into inactive by-products [3]. One of the major limitations of these approaches, however, is the catalytic turnover of human BChE [3] and so there is emphasis on the identification of BChE homologues from other organisms, such as SmBChE1, that might offer improved detoxification activity in this regard.
Inhibition of BChE in the absence of DDVP results in parasite death so a bio-scavenging role is possibly not the only function of this enzyme. Indeed, vertebrate BChE has also been shown to have roles in (1) ACh hydrolysis in situations of AChE deficiency [63], (2) fat metabolism by hydrolyzing the feeding stimulant peptide octanoyl ghrelin [64], and (3) scavenging polyproline-rich peptides to regulate protein-protein and protein-DNA interactions [65].
In summary, the work herein has identified multiple ChE paralogs in the genome of S. mansoni where previous studies, making use of the technology available at the time, attributed ChE activity to a single AChE. Consistent with previous observations that ChEs are multi-faceted enzymes, we posit that the three ChE paralogs described herein may fulfil distinct neuronal and non-neuronal functions based on their anatomical and temporal expression in the parasite and its ES products and the enzymatic activity of recombinant molecules. In addition to providing valuable insight into the functionality of individual ChE molecules, the study herein documents the essentiality of these proteins, providing a compelling evidence base for their use as intervention targets against schistosomiasis.

Ethics statement
All experimental procedures reported in the study were approved by the James Cook University (JCU) animal ethics committee (ethics approval numbers A2271 and A2391). Mice were maintained in cages in the university's quarantine facility (Q2152) for the duration of the experiments. The study protocols were in accordance with the 2007 Australian Code of Practice for the Care and Use of Animals for Scientific Purposes and the 2001 Queensland Animal Care and Protection Act.

Parasite maintenance, culture and ES collection
Biomphalaria glabrata snails infected with S. mansoni (NMRI strain) were obtained from the Biomedical Research Institute (BRI), MD, USA. Cercariae were shed from infected snails through exposure to light at 28˚C for 1.5 h and were mechanically transformed into schistosomula [66]. To obtain adult worms, 6-8 week old male BALB/c mice (Animal Resource Centre, WA, Australia) were infected with 180 cercariae via tail penetration and adults were harvested by vascular perfusion at 7-8 weeks post-infection [67]. Both adult worms and schistosomula were cultured (10 adult pairs/ml and 2000 schistosomula/ml) at 37˚C and 5% CO 2 in serumfree modified Basch medium [68] supplemented with 4 × antibiotic/antimycotic (AA-200 units/ml penicillin, 200 μg/ml streptomycin and 0.5 μg/ml amphotericin B) (SFB) in 6 well plates. Media containing ES products was initially collected after 3 h for schistosomula or 24 h for adults, and replenished daily thereafter. ES products were stored at -80˚C. Media was thawed when needed, concentrated through Amicon centrifugation filters (Sigma) with a 3 kDa molecular weight cutoff (MWCO), buffer exchanged into phosphate buffered saline, pH 7.4 (PBS) and aliquoted. Protein concentration of ES products was determined using the Pierce BCA™ Protein Assay kit (Thermofisher). To collect cercarial ES products, freshly-shed cercariae were incubated in H 2 O (4000/ml) at 25˚C for 3 h. H 2 O was filtered through Whatman filter paper (11 μm) to remove cercariae and associated debris, and ES products were concentrated, quantified and stored as described for adult and schistosomula ES products.

Parasite extract preparation
To make PBS-soluble extracts, worms were homogenized in PBS (50 μl/adult worm pair or 50 μl/1000 schistosomula) at 4˚C using a TissueLyser II (Qiagen), homogenates were incubated overnight with mixing at 4˚C and the supernatants collected by centrifugation at 15,000 g for 1 h at 4˚C. Triton X-100-soluble extracts were made from the PBS-insoluble pellets by resuspension in 1% Triton X-100, 40 mM Tris-HCl, pH 7.4, mixing overnight at 4˚C and the supernatant collected by centrifugation at 15,000 g for 1 h at 4˚C. Tegument extraction was achieved using a combination of freeze/thaw/vortex [69]. In brief, parasites were slowly thawed on ice, washed in TBS (10 mM Tris/HCl, 0.84% NaCl, pH 7.4) and incubated for 5 min on ice in 10 mM Tris/HCl, pH 7.4 followed by vortexing (5 × 1 s bursts). Subsequently, the tegumental extract was pelleted at 1000 g for 30 min and solubilized (3×) in 200 μl of 0.1% (w/v) SDS, 1% (v/v) Triton X-100 in 40 mM Tris, pH 7.4 with pelleting at 15,000 g between each wash. Protein concentration was determined using the Pierce BCA Protein Assay kit, aliquoted and stored at -80˚C until use.

Real-time qPCR
Real-time qPCR was used to assess developmental expression of smche genes and to determine smche transcript suppression resulting from RNAi experiments. RNA from miracidia, sporocysts, cercariae, adult male worms, adult female worms, and eggs were obtained from BRI. Schistosomula were cultured as described above, harvested (1,000 parasites) after either 3 h, 24 h, 3 or 5 days, washed three times in PBS and stored at -80˚C until use. Schistosomula from RNAi experiments were similarly processed. Total RNA extraction was performed using the Trizol (Thermofisher) reagent according to manufacturer's instructions. After air-drying, RNA pellets were re-suspended in 12 μl diethylpyrocarbonate (DEPC)-treated water. Concentration and purity of RNA was determined using an ND2000 Nanodrop spectrophotometer (Thermofisher). Synthesis of cDNA was carried out with 1 μg of total RNA using Superscript-III-Reverse Transcriptase (Invitrogen) according to the manufacturer's instructions. Finally, cDNA was quantified, diluted to 50 ng/μl, aliquoted and stored at -20˚C.
Real-time qPCR primers for each smche (S1 Table) were designed using Primer3 (http:// frodo.wi.mit.edu/). The housekeeping gene smcox1 was selected as an internal control to normalize relative smche gene expression [71]. Each qPCR (1 μl (50 ng) of cDNA, 5 μl of 2x SYBR green master mix (Bioline), 1μl (5 pmol/μl) each of forward and reverse primers and 2 μl of nuclease-free water) was run in a Rotor-Gene Q thermal cycler (Qiagen) using 40 cycles of 95˚C for 10 seconds, 50-55˚C for 15 seconds and 72˚C for 20 seconds. Stage-specific smche gene expression levels were normalized against smcox1 (smp_900000) expression using the comparative 2 -ΔΔCT method [72]. All results represent the average of 5 independent experiments with data presented as mean ± SEM.

Cloning, expression and purification of smche gene fragments in E. coli
Complete ORFs for smache1, smbche1 and smache2 were synthesized by Genewiz. Attempts to express full-length sequences in E. coli were unsuccessful, so primer sets incorporating NdeI (forward primer) and XhoI restriction enzyme sites (reverse primer) were designed (S1 Table) to amplify partial, non-conserved regions of each smche, which might prove more amenable to expression. Sequences (containing NdeI/ XhoI sites) for each SmChE were amplified from each full-length template by PCR and cloned into the pET41a expression vector (Novagen) such that the N-terminal GST tag was removed. Protein expression was induced for 24 h in E. coli BL21(DE3) by addition of 1 mM Isopropyl beta-D-1-thiogalactopyranoside (IPTG) using standard methods. Cultures were harvested by centrifugation (8,000 g for 20 min at 4˚C), resuspended in 50 ml lysis buffer (50 mM sodium phosphate, pH 8.0, 300 mM NaCl, 40 mM imidazole) and stored at -80˚C. Cell pellets were lysed by three freeze-thaw cycles at -80˚C and 42˚C followed by sonication on ice (10 × 5 s pulses [70% amplitude] with 30 s rest periods between each pulse) with a Qsonica Sonicator. Triton X-100 was added to each lysate at a final concentration of 3% and incubated for 1 h at 4˚C with end-over-end mixing. Insoluble material (containing SmChEs) was pelleted by centrifugation at 20,000 g for 20 min at 4˚C. The supernatant was discarded, and inclusion bodies (IBs) were washed twice by resuspension in 30 ml of lysis buffer followed by centrifugation at 20,000 g for 20 min at 4˚C. IBs were then solubilized sequentially by resuspension in 25 ml lysis buffers containing either 2 M, 4 M or 8 M urea, end-over-end mixing overnight at 4˚C and centrifugation at 20,000 g for 20 min at 4˚C. Finally, supernatant containing solubilized IBs was diluted 1:4 in lysis buffer containing 8M urea and filtered through a 0.22 μm membrane (Millipore). Solubilized IBs were purified by immobilized metal affinity chromatography (IMAC) by loading onto a prepacked 1 ml His-Trap HP column (GE Healthcare) equilibrated with lysis buffer containing 8M urea at a flow rate of 1 ml/min using an AKTA-pure-25 FPLC (GE Healthcare). After washing with 20 ml lysis buffer containing 8M urea, bound His-tagged proteins were eluted using the same buffer with a stepwise gradient of 50-250 mM imidazole (50 mM steps). Fractions containing SmChEs (as determined by SDS-PAGE) were pooled and concentrated using Amicon Ultra-15 centrifugal devices with a 3 kDa MWCO and quantified using the Pierce BCA Protein Assay kit. The final concentration of each SmChE was adjusted to 1 mg/ml and proteins were aliquoted and stored at -80˚C.

Generation of anti-rSmChE antisera and purification of IgG
Three groups of five male BALB/c mice (6-week-old) were intraperitoneally immunized with either SmAChE1, SmBChE1 or SmAChE2 subunits (50 μg/mouse). Antigens were mixed with an equal volume of Imject alum adjuvant (Thermofisher) and administered three times, two weeks apart. Two weeks after the final immunization, mice were sacrificed and blood was collected via cardiac puncture. Blood from all mice in each group was pooled and serum was separated by centrifugation after clotting and stored at −20˚C. Polyclonal antibodies were purified from mouse sera using Protein A Sepharose-4B (Thermofisher) according to the manufacturer's instructions. Serum from naïve mice was similarly processed.

Immunolocalization using anti-SmChE antisera
Adult worm sections. Freshly perfused adult S. mansoni and S. haematobium worms were fixed in 4% paraformaldehyde, embedded in paraffin and sections (7 μm thick) were cut in a cryostat. Following deparaffinization in xylene and rehydration in an ethanol series, antigen retrieval was performed by boiling the slides in 10 mM sodium citrate, pH 6.0, for 40 min followed by a solution of 10 mM Tris, 1 mM EDTA, 0.05% Tween, pH 9.0, for 20 min. All sections were then blocked with 10% heat-inactivated goat serum for 1 h RT. After washing 3 times with PBST, sections were incubated with anti-SmAChE1, anti-SmBChE1, anti-SmAChE2, naïve sera (negative control), S. mansoni or S. haematobium infected mouse sera (positive controls) (1:50 in PBST) overnight at 4˚C and then washed again (3 × 5 min each). Finally, the sections were incubated with goat-anti-mouse IgG-alexafluor647 (Sigma) (1:200 in PBST) for 1 h in the dark at RT. After a final washing step, slides were mounted with coverslips in Entellan mounting medium (Millipore). Fluorescence and bright-field microscopy were performed with an AxioImager M1 fluorescence microscope (Zeiss) using 10× and 20× objectives.
Live schistosomula. In vitro cultured living cercariae and schistosomula (3 h, 24 h, 3 and 5 days old) were harvested, washed with PBS and then blocked with PBST/10% heat-inactivated goat serum for 30 min at RT. Following three washes, the larvae were incubated with anti-SmAChE1, anti-SmBChE1, anti-SmAChE2 or naïve serum (negative control) (1:100 in PBST) overnight at 4˚C. Parasites were washed again before incubation with goat-anti-mouse IgG-alexafluor647 (Sigma) (1: 200 in PBST) for 1 h in the dark at RT, followed by 3 washes. Finally, schistosomes were fixed in 4% paraformaldehyde and transferred to a microscope slide for fluorescence microscopy using an AxioImager M1 fluorescence microscope.

Cloning and expression of full-length SmChEs in P. pastoris
Full-length sequences (minus the signal peptide) of SmAChE1, SmBChE1 and SmAChE2 were EcoRI/XbaI cloned into the C-terminal 6-His-tagged pPICZαA expression vector (Invitrogen) to facilitate secretory expression. Recombinant plasmids (20 μg) were linearized with Pmel (SmAChE1 & SmBChE1) and SacI (SmAChE2), purified by ethanol precipitation and resuspended in 15 μl of H 2 O. Linearized vectors were electroporated according to the manufacturer's instructions into P. pastoris X-33 cells (Thermofisher) in 2 mm cuvettes (2 ms, 2000V, 25 μF, 200 O, square wave pulse), using a Gene Pulser Xcell (Bio-Rad), plated onto YPDS agar plates containing 100 μg/ml zeocin and incubated for 3 days at 30˚C. Resultant colonies were then picked and patched onto YPDS agar containing 2 mg/ml zeocin and plates incubated at 30˚C until colonies were visible. A high-expressing clone of each SmChE (determined from pilot expression experiments) was used to inoculate 5 ml BMGY media supplemented with 50 ug/ml zeocin and grown overnight at 30˚C with rotation at 250 rpm. The entire culture was then used to inoculate 250 ml of BMGY in a 2L baffled flask and incubation was continued for 24 h at 30˚C. Cells were pelleted at 5000 g for 20 min at RT, re-suspended in 1L of BMMY (to induce protein expression) and split between 2 × 2L baffled flasks, which were incubated with shaking (250 rpm) at 30˚C for 72 h. Methanol was added to a final concentration of 0.5% (2.5 ml/flask) every 24 h to maintain induction of protein expression. Culture medium containing the secreted SmChE proteins was harvested by centrifugation (5000 g for 20 min at RT) and filtered through a 0.22 μm membrane filter (Millipore). Recombinant proteins were purified by IMAC using an AKTA-pure-25 FPLC (GE Healthcare). Briefly, culture medium was loaded onto a 5 ml His-excel column, pre-equilibrated with binding buffer (50 mM PBS pH 7.4, 300 mM NaCl), washed with 20 column volumes of binding buffer and then eluted with binding buffer containing a linear imidazole gradient (20 to 500 mM). The purity of fractions within the main peak was analyzed by SDS-PAGE and fractions of appropriate purity were pooled, concentrated and buffer exchanged into PBS using Amicon Ultra-15 centrifugal devices with a 3 kDa MWCO and quantified using the Pierce BCA Protein Assay kit. The final concentration of each SmChE was adjusted to 1 mg/ml and proteins were aliquoted and stored at -80˚C.

SmChE enzyme assays
Activity of SmChEs (10 μg), extracts (20 μg) and ES samples (20 μg) was determined by the Ellman method [73]; modified for use with 96 well microplates. Samples (parasite extracts, ES and SmChEs) were diluted in assay buffer (0.1M sodium phosphate, pH 7.4), and 2 mM acetylthiocholine (AcSCh) or butyrylthiocholine (BcSCh) (Sigma) and 0.5 mM 5, 5'-dithio-bis (2-nitrobenzoic acid) (DTNB) (Sigma) was added. The absorbance increase was monitored every 5 min at 405 nm in a Polarstar Omega microplate reader (BMG Labtech). Specific activity was calculated using the initial velocity of the reaction and extinction coefficient of 13,260 M −1 cm −1 for TNB. To investigate sensitivity of parasite ES products to AChE inhibitors, 25 μg of adult ES was pre-treated with 1 μM DDVP-active metabolite of the organophosphorous AChE inhibitor metrifonate-for 20 min at RT before measuring activity. Kinetic parameters of SmChEs were characterized by measuring enzyme activity at differing substrate concentrations and plotting enzyme activity [V] vs. substrate concentration [S]. The Km ([S] at 1/2 Vmax) was calculated using the Michaelis Menton equation. Enzyme assays with inhibitors were performed as above except that SmChEs in assay buffer were pre-treated with 1 μM DDVP, in the case of SmAChE1 and SmAChE2, or 1 mM iso-OMPA-a membrane-impermeable specific BChE inhibitor-in the case of SmBChE1, for another 20 min at RT. Experiments were performed in triplicate with data presented as the mean ± SEM.

Purification of secreted SmChEs from adult S. mansoni ES products
Affinity chromatography using edrophonium chloride-sepharose was used to purify SmAChE from S. mansoni based on the method of Hodgson and Chubb [74]. Briefly, 1 g of epoxy-activated sepharose 6B beads was washed with distilled H 2 O, the slurry centrifuged at 814 g for 5 min and the pellet gently resuspended in 50 mM sodium phosphate, pH 8.0, containing 200 mM edrophonium chloride (1:2 ratio of sepharose:edrophonium chloride). The pH of the solution was adjusted to 10.0 and coupling of edrophonium with the sepharose was facilitated by incubating the mixture overnight with shaking at 50˚C. The gel was then washed sequentially with 10 volumes each of 100 mM sodium acetate, pH 4.5, 12 mM sodium borate, pH 10.0, and distilled H 2 O and finally resuspended in distilled H 2 O to generate a 1 ml gel slurry. The gel slurry was packed into a chromatography column (10 cm long, 1 cm diameter) and equilibrated by gravity flow at 4˚C with 10 column volumes (CV) of equilibration buffer (50 mM phosphate buffer, pH 8.0). Approximately 20 ml of ES from adult S. mansoni (concentrated through a 10 kDa MWCO centrifugal filter from a starting volume of 500 ml of media, harvested each day for 7 days from 100 pairs of adult worms and buffer exchanged into equilibration buffer) was added to the column followed by washes with 20 CV of equilibration buffer and 20 CV of equilibration buffer containing 500 mM NaCl. Bound SmChE was then eluted with 10 CV of equilibration buffer containing 500 mM NaCl and 20 mM edrophonium chloride. The eluate was concentrated and buffer exchanged into PBS using a 10 kDa MWCO centrifugal filter (edrophonium chloride is an AChE inhibitor and would interfere with subsequent activity assays) and resolved by 10% SDS-PAGE to check purity and facilitate identification by mass spectrometry.

Mass spectrometric analysis of purified, secreted SmChE
Bands of interest were manually excised from the SDS polyacrylamide gel, washed with 50% acetonitrile and dried under vacuum at 30˚C. Cysteine residues were reduced with 20 mM DTT for 1 h at 65˚C followed by alkylation with 50 mM iodoacetamide for 40 min at 37˚C in the dark. In-gel trypsin digestion was performed at 37˚C overnight with 0.8 ng of trypsin in trypsin reaction buffer (40 mM ammonium bicarbonate, 9% acetonitrile). The supernatant was removed to a fresh microfuge tube and stored at 4˚C, and the remaining peptides were further extracted from the gel pieces by incubation with 0.1% trifluoracetic acid (TFA) at 37˚C for 45 min. The newly extracted supernatant was combined with the previously collected supernatant, then dried under vacuum. Prior to the matrix-assisted laser desorption/ionization-timeof-flight mass spectrometry (MALDI-TOF MS) analysis, peptides were concentrated and desalted using ZipTips (Millipore) following the manufacturer's instructions. Tryptic peptides were re-dissolved in 10 μl 5% formic acid and 6 μl was injected onto a 50 mm 300 μm C18 trap column (Agilent Technologies) followed by an initial wash step with Buffer A (5% (v/v) ACN, 0.1% (v/v) formic acid) for 5 min at 30 μl/min. Peptides were eluted at a flow rate of 0.3 μl/min onto an analytical nano HPLC column (15 cm × 75 μm 300SBC18, 3.5 μm, Agilent Technologies). The eluted peptides were then separated by a 55-min gradient of buffer B (90/10 acetonitrile/ 0.1% formic acid) 1-40% followed by a 5 min steeper gradient from 40-80%. The mass spectrometer (ABSciex 5600 Triple Tof) was operated in data-dependent acquisition mode, in which full scan TOF-MS data was acquired over the range of 350-1400 m/z, and over the range of 80-1400 m/z for product-ion observed in the TOF-MS scan exceeding a threshold of 100 counts and a charge state of +2 to +5. Analyst 1.6.1 (ABSCIEX) software was used for data acquisition and analysis.
For protein identification, a database was built using the S. mansoni genome v5.0 [http:// www.genedb.org/Homepage/Smansoni] with the common repository of adventitious proteins (cRAP, http://www.thegpm.org/crap/) appended to it. Mascot v.2.5.1 (Matrix Science) was used for database search. Carbamidomethylation of Cys was set as a fixed modification and oxidation of Met and deamidation of Asn and Gln were set as variable modifications. MS and MS/MS tolerance were set at 10 ppm and 0.1 Da, respectively and only proteins with at least two unique peptides (each composed of at least seven amino acid residues) identified were considered reliably identified and used for analysis.

siRNA design and synthesis
Three short interfering RNA duplexes (siRNAs) targeting each of the three identified smche paralogs were designed (S2 Table) and checked to avoid off-target silencing by BLAST search using the S. mansoni genome. An irrelevant siRNA from firefly luciferase (luc) was selected as a negative control [75]. All siRNAs were commercially synthesized (Integrated DNA Technologies) and oligonucleoitdes were suspended to a concentration of 1 μg/μl in DEPC-treated water.

Electroporation of schistosomula with siRNA
Prior to electroporation, mechanically transformed schistosomula were cultured for 24 h (2,000 schistosomula/ml), at 37˚C and 5% CO 2 in SFB in 6 well plates. After 3 washes with PBS, schistosomula were re-suspended in modified Basch medium (3,000 schistosomula/ 100 μl) and 3,000 schistosomula were transferred into a Genepulser 4 mm electroporation cuvette (Bio-Rad) for every siRNA treatment (four) and timepoint (four for each siRNA treatment-1, 3, 5 and 7 days). Schistosomula were electroporated with 10 μg of either luc, smache1, smbche1 or smache2 siRNA (for individual knockdown experiments) or an equal combination of all three smche siRNAs (30 μg total) or 30 μg luc siRNA (for combinatorial knockdown experiments) using a Bio-Rad Gene Pulser Xcell (single 20 ms pulse-125 V, 25 μF capacitance, 200 O resistance, square wave electroporation) at RT, added to 24 well plates containing 1 ml pre-warmed SFB and incubated (37˚C, 5% CO 2 ) for 7 days. Schistosomula were harvested at each timepoint and approximately 1,000 parasites were used for qPCR analysis (to assess transcript knockdown), 1,700 parasites were used for protein extract preparation (to examine phenotypic knockdown) and 300 parasites were used for Trypan Blue exclusion assays (to determine parasite viability). All parasite material was generated and separately analyzed from 2 independent experiments.

Determination of viability in siRNA-treated schistosomula
Schistosomula (100 parasites/replicate) were harvested at each timepoint and viability was determined by Trypan Blue exclusion staining [76]. Briefly, schistosomula were stained with 0.16% Trypan Blue in PBS with gentle shaking for 30 min at RT and then excess stain was removed by multiple washes in PBS before fixing in 10% formalin. Parasites were counted under 10× objective and live parasites (which had not taken up stain) were expressed as a percentage of total worms. Each assay was performed in triplicate.

Evaluation of protein expression in siRNA-treated schistosomula
Western blots were performed with day-7 parasite extracts (20 μg) following standard procedures. The blots were probed with the anti-SmChE antibodies (1:1000 in PBST) generated herein. A polyclonal anti-Sm-paramyosin antibody [75] was used as a loading control. Expression levels were quantified by comparative densitometry of bands resulting from luc-and smche-treated parasite extracts (n = 1) using ImageJ software.

Glucose uptake in schistosomula treated with siRNA
In a separate RNAi experiment, newly transformed schistosomula (5,000/treatment) were incubated for 5 days in SFB. Parasites were then electroporated with siRNAs as described above and finally transferred to serum-free DMEM (1 mg/ml glucose) supplemented with 4×AA. Media (50 μl) from each experiment was collected 48 h post-treatment and the amount of glucose was quantified using a colorimetric glucose assay kit (Sigma) following the manufacturer's instructions. Parasite viability at this timepoint was determined by Trypan Blue exclusion and transcript levels of each smche, as well as the glucose transporters sgtp1 (smp_012440) and sgtp4 (smp_105410), were also measured. Glucose levels were normalized according to the number of viable parasites and expressed relative to the luc group. Data is the average of 2 biological and 3 technical replicates ± SEM.

Infection of mice with SmChE siRNA-treated schistosomula
One-day-old schistosomula (10,000) were electroporated as above in 500 μl of SFB with 50 μg of either luc, smache1, smbche1 or smache2 siRNA or a combination of all three smche siRNAs (150 μg total). Parasites were injected intramuscularly into both thighs (1,000 per thigh) of male 6-8 week BALB/c mouse (5 mice per treatment group) using a 23-gauge needle. A control group of mice were similarly injected with non-electroporated schistosomula. Adult worms were perfused 20 days later to assess the number of worms that had matured and reached the mesenteries. Experiments were performed independently in duplicate. After trial 1, transcript levels of each smche from surviving worms were assessed using real-time qPCR. Terminal bleeds were taken from all mice from each trial and serum (1:500 in PBST) was analysed by ELISA for the presence of anti-schistosome antibodies using ELISA plates coated with cercarial transformation fluid [46].

Bio-scavenging of carboxylic esters by SmBChE1
To test the hypothesis that SmBChE1 may play a role in the bio-scavenging of AChE-inhibitory molecules, we first sought to determine whether BChE activity of schistosome extracts or recombinant SmBChE1 would be inhibited by DDVP. Schistosomula extracts (20 μg) or SmBChE1 (10 μg) were diluted in assay buffer (final reaction volume was 200 μl), DDVP was added to a final concentration of either 1, 2 or 5 μM (extract reactions) or 1, 5 or 10 μM (SmBChE1 reactions) and the reactions incubated for 20 mins at RT. BcSCh (final concentration 2 mM) and DTNB (final concentration 0.5 mM) were then added and the absorbance was monitored every 5 min at 405 nm in a Polarstar Omega microplate reader (BMG Labtech).
Next, we investigated whether inhibition of BChE activity would potentiate the AChEinhibitory and anti-schistosome effects of organophosphates (OP)s. Schistosomula extracts (20 μg) were diluted in assay buffer (final reaction volume was 200 μl), then iso-OMPA was added to a final concentration of either 1 or 2 mM and incubated for 20 min at RT. DDVP was then added to a final concentration of 1 μM and the samples were further incubated for 20 min at RT. AcSCh (final concentration 2 mM) and DTNB (final concentration 0.5 mM) were then added and the absorbance was monitored every 5 min at 405 nm. Extracts that were not treated with iso-OMPA with or without DDVP treatment were used as controls. Experiments were performed in triplicate with data presented as the mean ± SEM.
The same experiments were performed on live schistosomula using either an inhibitor-or RNAi-based approach. For the inhibitor-based experiment, 24 h schistosomula (1,000/treatment in 1 ml SFB) were pretreated with iso-OMPA at the non-lethal concentration of 100 μM and, 1 h after iso-OMPA treatment, schistosomula were treated with 1 μM DDVP and cultured for 5 h at 37˚C in 5% CO 2 . Parasites that were untreated, not treated with iso-OMPA but treated with DDVP and treated with iso-OMPA alone were used as controls. The parasites treated with iso-OMPA alone were assayed for AChE activity as described above to investigate any effect of iso-OMPA on AChE activity. For the RNAi-based experiment, 24 h schistosomula (1,500/100 μl SFB) were electroporated with 10 μg of either smbche1 or luc siRNA as described above, added to 24 well plates containing 1 ml pre-warmed SFB and incubated (37˚C, 5% CO 2 ) for 3 days before being treated with 1 μM DDVP and cultured for a further 5 h. Identical experiments in the absence of DDVP were performed to control for the effects of RNAi alone. For both inhibitor-and RNAi-based experiments, schistosomula viability was determined using Trypan Blue staining. Data is presented as the mean ± SEM of 2 biological and 3 technical replicates.
In a reverse testing of the bio-scavenging hypothesis, we sought to determine whether addition of SmBChE could mitigate the effects of DDVP. Ten micrograms of SmBChE1 was preincubated with 1 μM final concetration DDVP in AChE assay buffer (170 μl final volume) for 20 min at RT. Schistosomula extracts (20 μg), ACh (final concentration 2 mM) and DTNB (final concentration 0.5 mM) were then added and the absorbance was monitored every 5 min at 405 nm in a Polarstar Omega microplate reader. Reactions without SmBChE1 or without DDVP were used as controls. Experiments were performed in triplicate with data presented as the mean ± SEM. Again, the same experiments were performed on live schistosomula. After the pre-treatment of different amounts of SmBChE1 (10, 5, and 2.5 μg) with 1 μM DDVP in 500 μl SFB, 24 h schistosomula (1,000/treatment in 500 μl SFB) were added, incubated at 37˚C and 5% CO 2 for 24 h and then parasite viability was measured by Trypan Blue staining. Experiments where a similarly expressed and purified, but irrelevant, protein (SmTSP2 -smp_181530) was used instead of SmBChE1, and schistosomula cultured in media alone, were used as controls. Data is presented as the mean ± SEM of 2 biological and 3 technical replicates.

Statistical analyses
Data were reported as the means ± SEM. Statistical differences were assessed using the student's t test. P values less than 0.05 were considered statistically significant. and then DDVP, before assaying AsCSh conversion activity. Identical experiments were performed in the absence of DDVP to control for AsCSh conversion activity by BChE. (B) Schistosomula extracts were pre-incubated with fSmBChE1 (10 μg) then treated with DDVP (1 μM), or treated with DDVP alone, before assaying AcSCh conversion activity. Identical experiments were performed in the absence of DDVP to control for AsCSh conversion activity by BChE. For all assays, data are the average of triplicate biological and technical experiments ± SEM and differences were measured by the student's t test. �� P � 0.01, ��� P � 0.001. (TIF) S1