A systematically improved high quality genome and transcriptome of the human blood fluke Schistosoma mansoni

Vesicle-based secretion in schistosomes: Analysis of protein and microRNA (miRNA) content of exosome-like vesicles derived from Schistosoma mansoni

Exosomes are small vesicles of endocytic origin, which are released into the extracellular environment and mediate a variety of physiological and pathological conditions. Here we show that Schistosoma mansoni releases exosome-like vesicles in vitro. Vesicles were purified from culture medium by sucrose gradient fractionation and fractions containing vesicles verified by western blot analyses and electron microscopy. Proteomic analyses of exosomal contents unveiled 130 schistosome proteins. Among these proteins are common exosomal markers such as heat shock proteins, energy-generating enzymes, cytoskeletal proteins, and others. In addition, the schistosome extracellular vesicles contain proteins of potential importance for host-parasite interaction, notably peptidases, signaling proteins, cell adhesion proteins (e.g., integrins) and previously described vaccine candidates, including glutathione-S-transferase (GST), tetraspanin (TSP-2) and calpain. S. mansoni exosomes also contain 143 microRNAs (miRNA), of which 25 are present at high levels, including miRNAs detected in sera of infected hosts. Quantitative PCR analysis confirmed the presence of schistosome-derived miRNAs in exosomes purified from infected mouse sera. The results provide evidence of vesicle-mediated secretion in these parasites and suggest that schistosome-derived exosomes could play important roles in host-parasite interactions and could be a useful tool in the development of vaccines and therapeutics.

Evolutionary fates of universal stress protein paralogs in Platyhelminthes

Background

Universal stress proteins (USPs) are present in all domains of life. Their expression is upregulated in response to a large variety of stress conditions. The functional diversity found in this protein family, paired with the sequence degeneration of the characteristic ATP-binding motif, suggests a complex evolutionary pattern for the paralogous USP-encoding genes. In this work, we investigated the origin, genomic organization, expression patterns and evolutionary history of the USP gene family in species of the phylum Platyhelminthes.

Results

Our data showed a cluster organization, a lineage-specific distribution, and the presence of several pseudogenes among the USP gene copies identified. The absence of a well conserved -CCAATCA- motif in the promoter region was positively correlated with low or null levels of gene expression, and with amino acid changes within the ligand binding motifs. Despite evidence of the pseudogenization of various USP genes, we detected an important functional divergence at several residues, mostly located near sites that are critical for ligand interaction.

Conclusions

Our results provide a broad framework for the evolution of the USP gene family, based on the emergence of new paralogs that face very contrasting fates, including pseudogenization, subfunctionalization or neofunctionalization. This framework aims to explain the sequence and functional diversity of this gene family, providing a foundation for future studies in other taxa in which USPs occur.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1129-x) contains supplementary material, which is available to authorized users.

Functional characterization of the Ca2+-ATPase SMA1 from Schistosoma mansoni

Schistosoma mansoni is a parasite that causes bilharzia, a neglected tropical disease affecting hundreds of millions of people each year worldwide. In 2012, S. mansoni had been identified as the only invertebrate possessing two SERCA-type Ca2+-ATPases, SMA1 and SMA2. However, our analysis of recent genomic data shows that the presence of two SERCA pumps is rather frequent in parasitic flatworms. To understand the reasons of this redundancy in S. mansoni, we compared SMA1 and SMA2 at different levels. In terms of sequence and organization, the genes SMA1 and SMA2 are similar, suggesting that they might be the result of a duplication event. At the protein level, SMA1 and SMA2 only slightly differ in length and in the sequence of the nucleotide-binding domain. To get functional information on SMA1, we produced it in an active form in Saccharomyces cerevisiae, as previously done for SMA2. Using phosphorylation assays from ATP, we demonstrated that like SMA2, SMA1 bound calcium in a cooperative mode with an apparent affinity in the micromolar range. We also showed that SMA1 and SMA2 had close sensitivities to cyclopiazonic acid but different sensitivities to thapsigargin, two specific inhibitors of SERCA pumps. On the basis of transcriptomic data available in GeneDB, we hypothesize that SMA1 is a housekeeping Ca2+-ATPase, whereas SMA2 might be required in particular striated-like muscles like those present the tail of the cercariae, the infecting form of the parasite.

Tissue-specific transcriptome analyses provide new insights into GPCR signalling in adult Schistosoma mansoni

Schistosomes are blood-dwelling trematodes with global impact on human and animal health. Because medical treatment is currently based on a single drug, praziquantel, there is urgent need for the development of alternative control strategies. The Schistosoma mansoni genome project provides a platform to study and connect the genetic repertoire of schistosomes to specific biological functions essential for successful parasitism. G protein-coupled receptors (GPCRs) form the largest superfamily of transmembrane receptors throughout the Eumetazoan phyla, including platyhelminths. Due to their involvement in diverse biological processes, their pharmacological importance, and proven druggability, GPCRs are promising targets for new anthelmintics. However, to identify candidate receptors, a more detailed understanding of the roles of GPCR signalling in schistosome biology is essential. An updated phylogenetic analysis of the S. mansoni GPCR genome (GPCRome) is presented, facilitated by updated genome data that allowed a more precise annotation of GPCRs. Additionally, we review the current knowledge on GPCR signalling in this parasite and provide new insights into the potential roles of GPCRs in schistosome reproduction based on the findings of a recent tissue-specific transcriptomic study in paired and unpaired S. mansoni. According to the current analysis, GPCRs contribute to gonad-specific functions but also to nongonad, pairing-dependent processes. The latter may regulate gonad-unrelated functions during the multifaceted male-female interaction. Finally, we compare the schistosome GPCRome to that of another parasitic trematode, Fasciola, and discuss the importance of GPCRs to basic and applied research. Phylogenetic analyses display GPCR diversity in free-living and parasitic platyhelminths and suggest diverse functions in schistosomes. Although their roles need to be substantiated by functional studies in the future, the data support the selection of GPCR candidates for basic and applied studies, invigorating the exploitation of this important receptor class for drug discovery against schistosomes but also other trematodes.

Gene Duplication Analysis Reveals No Ancient Whole Genome Duplication but Extensive Small-Scale Duplications during Genome Evolution and Adaptation of Schistosoma mansoni

Gene duplication (GD), thought to facilitate evolutionary innovation and adaptation, has been studied in many phylogenetic lineages. However, it remains poorly investigated in trematodes, a medically important parasite group that has been evolutionarily specialized during long-term host-parasite interaction. In this study, we conducted a genome-wide study of GD modes and contributions in Schistosoma mansoni, a pathogen causing human schistosomiasis. We combined several lines of evidence provided by duplicate age distributions, genomic sequence similarity, depth-of-coverage and gene synteny to identify the dominant drivers that contribute to the origins of new genes in this parasite. The gene divergences following duplication events (gene structure, expression and function retention) were also analyzed. Our results reveal that the genome lacks whole genome duplication (WGD) in a long evolutionary time and has few large segmental duplications, but is extensively shaped by the continuous small-scale gene duplications (SSGDs) (i.e., dispersed, tandem and proximal GDs) that may be derived from (retro-) transposition and unequal crossing over. Additionally, our study shows that the genes generated by tandem duplications have the smallest divergence during the evolution. Finally, we demonstrate that SSGDs, especially the tandem duplications, greatly contribute to the expansions of some preferentially retained pathogenesis-associated gene families that are associated with the parasite's survival during infection. This study is the first to systematically summarize the landscape of GDs in trematodes and provides new insights of adaptations to parasitism linked to GD events for these parasites.

The Schistosoma mansoni genome encodes thousands of long non-coding RNAs predicted to be functional at different parasite life-cycle stages

Next Generation Sequencing (NGS) strategies, like RNA-Seq, have revealed the transcription of a wide variety of long non-coding RNAs (lncRNAs) in the genomes of several organisms. In the present work we assessed the lncRNAs complement of Schistosoma mansoni, the blood fluke that causes schistosomiasis, ranked among the most prevalent parasitic diseases worldwide. We focused on the long intergenic/intervening ncRNAs (lincRNAs), hidden within the large amount of information obtained through RNA-Seq in S. mansoni (88 libraries). Our computational pipeline identified 7029 canonically-spliced putative lincRNA genes on 2596 genomic loci (at an average 2.7 isoforms per lincRNA locus), as well as 402 spliced lncRNAs that are antisense to protein-coding (PC) genes. Hundreds of lincRNAs showed traits for being functional, such as the presence of epigenetic marks at their transcription start sites, evolutionary conservation among other schistosome species and differential expression across five different life-cycle stages of the parasite. Real-time qPCR has confirmed the differential life-cycle stage expression of a set of selected lincRNAs. We have built PC gene and lincRNA co-expression networks, unraveling key biological processes where lincRNAs might be involved during parasite development. This is the first report of a large-scale identification and structural annotation of lncRNAs in the S. mansoni genome.

Phenotypic, chemical and functional characterization of cyclic nucleotide phosphodiesterase 4 (PDE4) as a potential anthelmintic drug target

BACKGROUND

Reliance on just one drug to treat the prevalent tropical disease, schistosomiasis, spurs the search for new drugs and drug targets. Inhibitors of human cyclic nucleotide phosphodiesterases (huPDEs), including PDE4, are under development as novel drugs to treat a range of chronic indications including asthma, chronic obstructive pulmonary disease and Alzheimer's disease. One class of huPDE4 inhibitors that has yielded marketed drugs is the benzoxaboroles (Anacor Pharmaceuticals).

METHODOLOGY/PRINCIPAL FINDINGS

A phenotypic screen involving Schistosoma mansoni and 1,085 benzoxaboroles identified a subset of huPDE4 inhibitors that induced parasite hypermotility and degeneration. To uncover the putative schistosome PDE4 target, we characterized four PDE4 sequences (SmPDE4A-D) in the parasite's genome and transcriptome, and cloned and recombinantly expressed the catalytic domain of SmPDE4A. Among a set of benzoxaboroles and catechol inhibitors that differentially inhibit huPDE4, a relationship between the inhibition of SmPDE4A, and parasite hypermotility and degeneration, was measured. To validate SmPDE4A as the benzoxaborole molecular target, we first generated Caenorhabditis elegans lines that express a cDNA for smpde4a on a pde4(ce268) mutant (hypermotile) background: the smpde4a transgene restored mutant worm motility to that of the wild type. We then showed that benzoxaborole inhibitors of SmPDE4A that induce hypermotility in the schistosome also elicit a hypermotile response in the C. elegans lines that express the smpde4a transgene, thereby confirming SmPDE4A as the relevant target.

CONCLUSIONS/SIGNIFICANCE

The orthogonal chemical, biological and genetic strategies employed identify SmPDE4A's contribution to parasite motility and degeneration, and its potential as a drug target. Transgenic C. elegans is highlighted as a potential screening tool to optimize small molecule chemistries to flatworm molecular drug targets.

Characterization of the complete nuclear ribosomal DNA sequences of Eurytrema pancreaticum

Eurytrema pancreaticum is one of the most common trematodes of cattle and sheep, and also infects humans occasionally, causing great economic losses and medical costs. In this study, the sequences of the complete nuclear ribosomal DNA (rDNA) repeat units of five E. pancreaticum individuals were determined for the first time. They were 8306–8310 bp in length, including the small subunit (18S) rDNA, internal transcribed spacer 1 (ITS1), 5.8S rDNA, internal transcribed spacer 2 (ITS2), large subunit (28S) rDNA and intergenic spacer (IGS). There were no length variations in any of the investigated 18S (1996 bp), ITS1 (1103 bp), 5.8S (160 bp), ITS2 (231 bp) or 28S (3669 bp) rDNA sequences, whereas the IGS rDNA sequences of E. pancreaticum had a 4-bp length variation, ranging from 1147 to 1151 bp. The intraspecific variations within E. pancreaticum were 0–0.2% for 18S rDNA, 0–0.5% for ITS1, 0% for 5.8S rDNA and ITS2, 0–0.2% for 28S rDNA and 2.9–20.2% for IGS. There were nine types of repeat sequences in ITS1, two types in 28S rDNA, but none in IGS. A phylogenetic analysis based on the 18S rDNA sequences classified E. pancreaticum in the family Dicrocoeliidae of Plagiorchiata, closely related to the suborder Opisthorchiata. These results provide useful information for the further study of Dicrocoeliidae trematodes.

Comparative sequence analysis reveals regulation of genes in developing schistosomula of Schistosoma mansoni exposed to host portal serum

Once inside a vertebrate host after infection, individual schistosomula of the parasite Schistosoma mansoni find a new and complex environment, which requires quick adjustments for survival, such as those that allow it to avoid the innate immune response of the host. Thus, it is very important for the parasite to remain within the skin after entering the host for a period of about 3 days, at which time it can then reach the venous system, migrate to the lungs and, by the end of eighth day post-infection, it reach the portal venous system, while undergoing minimal changes in morphology. However, after just a few days in the portal blood system, the parasite experiences an extraordinary increase in biomass and significant morphological alterations. Therefore, determining the constituents of the portal venous system that may trigger these changes that causes the parasite to consolidate its development inside the vertebrate host, thus causing the disease schistosomiasis, is essential. The present work simulated the conditions found in the portal venous system of the vertebrate host by exposing schistosomula of S. mansoni to in vitro culture in the presence of portal serum of the hamster, Mesocricetus auratus. Two different incubation periods were evaluated, one of 3 hours and one of 12 hours. These time periods were used to mimic the early contact of the parasite with portal serum during the course of natural infection. As a control, parasites were incubated in presence of hamster peripheral serum, in order to compare gene expression signatures between the two conditions. The mRNA obtained from parasites cultured under both conditions were submitted to a whole transcriptome library preparation and sequenced with a next generation platform. On average, nearly 15 million reads were produced per sample and, for the purpose of gene expression quantification, only reads mapped to one location of the transcriptome were considered. After statistical analysis, we found 103 genes differentially expressed by schistosomula cultured for 3 hours and 12 hours in the presence of hamster portal serum. After the subtraction of a second list of genes, also differentially expressed between schistosomula cultured for 3 hours and 12 hours in presence of peripheral serum, a set of 58 genes was finally established. This pattern was further validated for a subset of 17 genes, by measuring gene expression through quantitative real time polymerase chain reaction (qPCR). Processes that were activated by the portal serum stimulus include response to stress, membrane transport, protein synthesis and folding/degradation, signaling, cytoskeleton arrangement, cell adhesion and nucleotide synthesis. Additionally, a smaller number of genes down-regulated under the same condition act on cholinergic signaling, inorganic cation and organic anion membrane transport, cell adhesion and cytoskeleton arrangement. Considering the role of these genes in triggering processes that allow the parasite to quickly adapt, escape the immune response of the host and start maturation into an adult worm after contact with the portal serum, this work may point to unexplored molecular targets for drug discovery and vaccine development against schistosomiasis.

Proteomic Analysis of Trichinella spiralis Adult Worm Excretory-Secretory Proteins Recognized by Sera of Patients with Early Trichinellosis

The most commonly used serodiagnostic antigens for trichinellosis are the excretory-secretory (ES) antigens from T. spiralis muscle larvae (ML), but the specific antibodies against the ML ES antigens are usually negative during early stage of Trichinella infection. The recent studies demonstrated that T. spiralis adult worm (AW) antigens were recognized by mouse or swine infection sera on Western blot as early as 7–15 days post-infection (dpi), the AW antigens might contain the early diagnostic markers for trichinellosis. The purpose of this study was to screen early diagnostic antigens in T. spiralis AW ES proteins recognized by sera of early patients with trichinellosis. T. spiralis AW were collected at 72 h post-infection (hpi), and their ES antigens were analyzed by SDS-PAGE and Western blot. Our results showed that 5 protein bands (55, 48–50, 45, 44, and 36 kDa) were recognized by sera of early patients with trichinellosis collected at 19 dpi, and were subjected to shotgun LC–MS/MS and bioinformatics analyses. A total of 185 proteins were identified from T. spiralis protein database, of which 116 (67.2%) proteins had molecular weights of 30∼60 kDa, and 125 (67.6%) proteins with pI 4–7. Bioinformatic analyses showed that the identified proteins have a wide diversity of biological functions (binding of nucleotides, proteins, ions, carbohydrates, and lipids; hydrolase, transferase, and oxidoreductase, etc.). Several enzymes (e.g., adult-specific DNase II, serine protease and serine protease inhibitor) could be the invasion-related proteins and early diagnostic markers for trichinellosis. Moreover, recombinant T. spiralis serine protease (rTsSP-ZH68) was expressed in E. coli and its antigenicity was analyzed by Western blot with the early infection sera. The rTsSP-ZH68 was recognized by sera of infected mice at 8–10 dpi and sera of early patients with trichinellosis at 19 dpi. T. spiralis AW proteins identified in this study, especially serine protease, are the promising early diagnostic antigens and vaccine candidates for trichinellosis.

MiR-277/4989 regulate transcriptional landscape during juvenile to adult transition in the parasitic helminth Schistosoma mansoni

Schistosomes are parasitic helminths that cause schistosomiasis, a disease affecting circa 200 million people, primarily in underprivileged regions of the world. Schistosoma mansoni is the most experimentally tractable schistosome species due to its ease of propagation in the laboratory and the high quality of its genome assembly and annotation. Although there is growing interest in microRNAs (miRNAs) in trematodes, little is known about the role these molecules play in the context of developmental processes. We use the completely unaware "miRNA-blind" bioinformatics tool Sylamer to analyse the 3'-UTRs of transcripts differentially expressed between the juvenile and adult stages. We show that the miR-277/4989 family target sequence is the only one significantly enriched in the transition from juvenile to adult worms. Further, we describe a novel miRNA, sma-miR-4989 showing that its proximal genomic location to sma-miR-277 suggests that they form a miRNA cluster, and we propose hairpin folds for both miRNAs compatible with the miRNA pathway. In addition, we found that expression of sma-miR-277/4989 miRNAs are up-regulated in adults while their predicted targets are characterised by significant down-regulation in paired adult worms but remain largely undisturbed in immature "virgin" females. Finally, we show that sma-miR-4989 is expressed in tegumental cells located proximal to the oesophagus gland and also distributed throughout the male worms' body. Our results indicate that sma-miR-277/4989 might play a dominant role in post-transcriptional regulation during development of juvenile worms and suggest an important role in the sexual development of female schistosomes.

A novel ex vivo immunoproteomic approach characterising Fasciola hepatica tegumental antigens identified using immune antibody from resistant sheep

A more thorough understanding of the immunological interactions between Fasciola spp. and their hosts is required if we are to develop new immunotherapies to control fasciolosis. Deeper knowledge of the antigens that are the target of the acquired immune responses of definitive hosts against both Fasciola hepatica and Fasciola gigantica will potentially identify candidate vaccine antigens. Indonesian Thin Tail sheep express a high level of acquired immunity to infection by F. gigantica within 4weeks of infection and antibodies in Indonesian Thin Tail sera can promote antibody-dependent cell-mediated cytotoxicity against the surface tegument of juvenile F. gigantica in vitro. Given the high protein sequence similarity between F. hepatica and F. gigantica, we hypothesised that antibody from F. gigantica-infected sheep could be used to identify the orthologous proteins in the tegument of F. hepatica. Purified IgG from the sera of F. gigantica-infected Indonesian Thin Tail sheep collected pre-infection and 4weeks p.i. were incubated with live adult F. hepatica ex vivo and the immunosloughate (immunoprecipitate) formed was isolated and analysed via liquid chromatography-electrospray ionisation-tandem mass spectrometry to identify proteins involved in the immune response. A total of 38 proteins were identified at a significantly higher abundance in the immunosloughate using week 4 IgG, including eight predicted membrane proteins, 20 secreted proteins, nine proteins predicted to be associated with either the lysosomes, the cytoplasm or the cytoskeleton and one protein with an unknown cellular localization. Three of the membrane proteins are transporters including a multidrug resistance protein, an amino acid permease and a glucose transporter. Interestingly, a total of 21 of the 38 proteins matched with proteins recently reported to be associated with the proposed small exosome-like extracellular vesicles of adult F. hepatica, suggesting that the Indonesian Thin Tail week 4 IgG is either recognising individual proteins released from extracellular vesicles or is immunoprecipitating intact exosome-like extracellular vesicles. Five extracellular vesicle membrane proteins were identified including two proteins predicted to be associated with vesicle transport/ exocytosis (VPS4, vacuolar protein sorting-associated protein 4b and the Niemann-Pick C1 protein). RNAseq analysis of the developmental transcription of the 38 immunosloughate proteins showed that the sequences are expressed over a wide abundance range with 21/38 transcripts expressed at a relatively high level from metacercariae to the adult life cycle stage. A notable feature of the immunosloughates was the absence of cytosolic proteins which have been reported to be secreted markers for damage to adult flukes incubated in vitro, suggesting that the proteins observed are not inadvertent contaminants leaking from damaged flukes ex vivo. The identification of tegument protein antigens shared between F. gigantica and F. hepatica is beneficial in terms of the possible development of a dual purpose vaccine effective against both fluke species.

Histone deacetylase inhibition modulates histone acetylation at gene promoter regions and affects genome-wide gene transcription in Schistosoma mansoni

BACKGROUND

Schistosomiasis is a parasitic disease infecting hundreds of millions of people worldwide. Treatment depends on a single drug, praziquantel, which kills the Schistosoma spp. parasite only at the adult stage. HDAC inhibitors (HDACi) such as Trichostatin A (TSA) induce parasite mortality in vitro (schistosomula and adult worms), however the downstream effects of histone hyperacetylation on the parasite are not known.

METHODOLOGY/PRINCIPAL FINDINGS

TSA treatment of adult worms in vitro increased histone acetylation at H3K9ac and H3K14ac, which are transcription activation marks, not affecting the unrelated transcription repression mark H3K27me3. We investigated the effect of TSA HDACi on schistosomula gene expression at three different time points, finding a marked genome-wide change in the transcriptome profile. Gene transcription activity was correlated with changes on the chromatin acetylation mark at gene promoter regions. Moreover, combining expression data with ChIP-Seq public data for schistosomula, we found that differentially expressed genes having the H3K4me3 mark at their promoter region in general showed transcription activation upon HDACi treatment, compared with those without the mark, which showed transcription down-regulation. Affected genes are enriched for DNA replication processes, most of them being up-regulated. Twenty out of 22 genes encoding proteins involved in reducing reactive oxygen species accumulation were down-regulated. Dozens of genes encoding proteins with histone reader motifs were changed, including SmEED from the PRC2 complex. We targeted SmEZH2 methyltransferase PRC2 component with a new EZH2 inhibitor (GSK343) and showed a synergistic effect with TSA, significantly increasing schistosomula mortality.

CONCLUSIONS/SIGNIFICANCE

Genome-wide gene expression analyses have identified important pathways and cellular functions that were affected and may explain the schistosomicidal effect of TSA HDACi. The change in expression of dozens of histone reader genes involved in regulation of the epigenetic program in S. mansoni can be used as a starting point to look for possible novel schistosomicidal targets.

(Epi)genetic Inheritance in Schistosoma mansoni: A Systems Approach

The G×E concept, in which genotype × environment interactions bring about the phenotype, is widely used to describe biological phenomena. We propose to extend the initial notion of the concept, replacing G by 'inheritance system'. This system, comprised of both genome and epigenome components, collectively interacts with the environment to shape the development of a phenotype. In the case of the human blood fluke Schistosoma mansoni, responsible for intestinal bilharzia, the phenotypic trait that is most relevant to global health is infection success. Taking a systems biology view we show how genetic and epigenetic interactions result in ephemeral, but also heritable, phenotypic variations that are important for infection success.

De Novo Transcriptome Characterization of a Sterilizing Trematode Parasite (Microphallus sp.) from Two Species of New Zealand Snails

Snail-borne trematodes represent a large, diverse, and evolutionarily, ecologically, and medically important group of parasites, often imposing strong selection on their hosts and causing host morbidity and mortality. Even so, there are very few genomic and transcriptomic resources available for this important animal group. We help to fill this gap by providing transcriptome resources from trematode metacercariae infecting two congeneric snail species, Potamopyrgus antipodarum and P. estuarinus. This genus of New Zealand snails has gained prominence in large part through the development of P. antipodarum and its sterilizing trematode parasite Microphallus livelyi into a textbook model for host–parasite coevolutionary interactions in nature. By contrast, the interactions between Microphallus trematodes and P. estuarinus, an estuary-inhabiting species closely related to the freshwater P. antipodarum, are relatively unstudied. Here, we provide the first annotated transcriptome assemblies from Microphallus isolated from P. antipodarum and P. estuarinus. We also use these transcriptomes to produce genomic resources that will be broadly useful to those interested in host–parasite coevolution, local adaption, and molecular evolution and phylogenetics of this and other snail–trematode systems. Analyses of the two Microphallus transcriptomes revealed that the two trematode types are more genetically differentiated from one another than are the M. livelyi infecting different populations of P. antipodarum, suggesting that the Microphallus infecting P. estuarinus represent a distinct lineage. We also provide a promising set of candidate genes likely involved in parasitic infection and response to salinity stress.

Functional diversity of secreted cestode Kunitz proteins: Inhibition of serine peptidases and blockade of cation channels

We previously reported a multigene family of monodomain Kunitz proteins from Echinococcus granulosus (EgKU-1-EgKU-8), and provided evidence that some EgKUs are secreted by larval worms to the host interface. In addition, functional studies and homology modeling suggested that, similar to monodomain Kunitz families present in animal venoms, the E. granulosus family could include peptidase inhibitors as well as channel blockers. Using enzyme kinetics and whole-cell patch-clamp, we now demonstrate that the EgKUs are indeed functionally diverse. In fact, most of them behaved as high affinity inhibitors of either chymotrypsin (EgKU-2-EgKU-3) or trypsin (EgKU-5-EgKU-8). In contrast, the close paralogs EgKU-1 and EgKU-4 blocked voltage-dependent potassium channels (Kv); and also pH-dependent sodium channels (ASICs), while showing null (EgKU-1) or marginal (EgKU-4) peptidase inhibitory activity. We also confirmed the presence of EgKUs in secretions from other parasite stages, notably from adult worms and metacestodes. Interestingly, data from genome projects reveal that at least eight additional monodomain Kunitz proteins are encoded in the genome; that particular EgKUs are up-regulated in various stages; and that analogous Kunitz families exist in other medically important cestodes, but not in trematodes. Members of this expanded family of secreted cestode proteins thus have the potential to block, through high affinity interactions, the function of host counterparts (either peptidases or cation channels) and contribute to the establishment and persistence of infection. From a more general perspective, our results confirm that multigene families of Kunitz inhibitors from parasite secretions and animal venoms display a similar functional diversity and thus, that host-parasite co-evolution may also drive the emergence of a new function associated with the Kunitz scaffold.

Evidence for Integrin - Venus Kinase Receptor 1 Alliance in the Ovary of Schistosoma mansoni Females Controlling Cell Survival

In metazoan integrin signaling is an important process of mediating extracellular and intracellular communication processes. This can be achieved by cooperation of integrins with growth factor receptors (GFRs). Schistosoma mansoni is a helminth parasite inducing schistosomiasis, an infectious disease of worldwide significance for humans and animals. First studies on schistosome integrins revealed their role in reproductive processes, being involved in spermatogenesis and oogenesis. With respect to the roles of eggs for maintaining the parasite´s life cycle and for inducing the pathology of schistosomiasis, elucidating reproductive processes is of high importance. Here we studied the interaction of the integrin receptor Smβ-Int1 with the venus kinase receptor SmVKR1 in S. mansoni. To this end we cloned and characterized SmILK, SmPINCH, and SmNck2, three putative bridging molecules for their role in mediating Smβ-Int1/SmVKR1 cooperation. Phylogenetic analyses showed that these molecules form clusters that are specific for parasitic platyhelminths as it was shown for integrins before. Transcripts of all genes colocalized in the ovary. In Xenopus oocytes germinal vesicle breakdown (GVBD) was only induced if all members were simultaneously expressed. Coimmunoprecipitation results suggest that a Smβ-Int1-SmILK-SmPINCH-SmNck2-SmVKR1 complex can be formed leading to the phosphorylation and activation of SmVKR1. These results indicate that SmVKR1 can be activated in a ligand-independent manner by receptor-complex interaction. RNAi and inhibitor studies to knock-down SmILK as a representative complex member concurrently revealed effects on the extracellular matrix surrounding the ovary and oocyte localization within the ovary, oocyte survival, and egg production. By TUNEL assays, confocal laser scanning microscopy (CLSM), Caspase-3 assay, and transcript profiling of the pro-apoptotic BCL-2 family members BAK/BAX we obtained first evidence for roles of this signaling complex in mediating cell death in immature and primary oocytes. These results suggest that the Smβ-Int1/SmVKR1 signaling complex is important for differentiation and survival in oocytes of paired schistosomes.

Parasites of the genus Schistosoma cause schistosomiasis, a life-threatening infectious disease for humans and animals worldwide. Among the remarkable biological features of schistosomes is the differentiation of the female gonads which is controlled by pairing with the male and a prerequisite for egg production. Eggs, however, are not only important for the maintenance of the life-cycle; they also cause the pathological consequences of schistosomiasis. Part of the eggs gets trapped in host tissues such as liver and spleen and trigger inflammatory processes, finally leading to liver cirrhosis. Research activities of the last decade have indicated that different families of cellular and receptor-type kinases but also integrins contribute to the control of mitogenic activity and differentiation the female goands. In this context an unusual class of receptor tyrosine kinases (RTKs) has been identified, the venus kinase receptors (SmVKRs). By biochemical and molecular approaches we demonstrate that SmVKR1 activation can be achieved by cooperation with a signaling complex consisting of the beta integrin receptor Smβ-Int1 and the bridging molecules SmILK, SmPINCH, SmNck2. Besides unravelling a novel way of SmVKR1 activation, we provide evidence that this complex controls the differentiation status of oocytes by regulating cell death-associated processes.

Whole Genome Amplification and Reduced-Representation Genome Sequencing of Schistosoma japonicum Miracidia

Background

In areas where schistosomiasis control programs have been implemented, morbidity and prevalence have been greatly reduced. However, to sustain these reductions and move towards interruption of transmission, new tools for disease surveillance are needed. Genomic methods have the potential to help trace the sources of new infections, and allow us to monitor drug resistance. Large-scale genotyping efforts for schistosome species have been hindered by cost, limited numbers of established target loci, and the small amount of DNA obtained from miracidia, the life stage most readily acquired from humans. Here, we present a method using next generation sequencing to provide high-resolution genomic data from S. japonicum for population-based studies.

Methodology/Principal Findings

We applied whole genome amplification followed by double digest restriction site associated DNA sequencing (ddRADseq) to individual S. japonicum miracidia preserved on Whatman FTA cards. We found that we could effectively and consistently survey hundreds of thousands of variants from 10,000 to 30,000 loci from archived miracidia as old as six years. An analysis of variation from eight miracidia obtained from three hosts in two villages in Sichuan showed clear population structuring by village and host even within this limited sample.

Conclusions/Significance

This high-resolution sequencing approach yields three orders of magnitude more information than microsatellite genotyping methods that have been employed over the last decade, creating the potential to answer detailed questions about the sources of human infections and to monitor drug resistance. Costs per sample range from $50-$200, depending on the amount of sequence information desired, and we expect these costs can be reduced further given continued reductions in sequencing costs, improvement of protocols, and parallelization. This approach provides new promise for using modern genome-scale sampling to S. japonicum surveillance, and could be applied to other schistosome species and other parasitic helminthes.

Schistosomiasis is a devastating tropical disease that affects more than 200 million people worldwide. Over the past several decades, transmission control strategies implemented in China have reduced the prevalence and morbidity of Schistosoma japonicum in many areas. Infections still persist, however, and it is therefore of great interest to determine the sources of recurring infections. Genetic analysis is a promising means to achieve this. Towards this aim, we conducted a pilot study to assess the feasibility of using high-throughput sequencing to assess the geographic distribution of schistosome genetic variants. Because DNA yields from miracidia, the most easily accessible life stage, are insufficient for high throughput sequencing, we first employed whole genome amplification to obtain sufficient quantities of DNA. We then employed a technique that reproducibly sequences the same fraction of a genome across numerous samples. We successfully sequenced 6-year old S. japonicum samples from Sichuan Province, China, easily and economically identifying tens of thousands of variable loci, a sufficient number to discriminate fine-scale population structure. Further population sampling will help answer important questions concerning the persistence of infections, the sources of new infections, and whether parasite populations have undergone incipient evolution of drug resistance.

Protease Inhibitors of Parasitic Flukes: Emerging Roles in Parasite Survival and Immune Defence

Protease inhibitors play crucial roles in parasite development and survival, counteracting the potentially damaging immune responses of their vertebrate hosts. However, limited information is currently available on protease inhibitors from schistosomes and food-borne trematodes. Future characterization of these molecules is important not only to expand knowledge on parasitic fluke biology but also to determine whether they represent novel vaccine and/or drug targets. Moreover, protease inhibitors from flukes may represent lead compounds for the development of a new range of therapeutic agents against inflammatory disorders and cancer. This review discusses already identified protease inhibitors of fluke origin, emphasizing their biological function and their possible future development as new intervention targets.

Echinococcus-Host Interactions at Cellular and Molecular Levels

The potentially lethal zoonotic diseases alveolar and cystic echinococcosis are caused by the metacestode larval stages of the tapeworms Echinococcus multilocularis and Echinococcus granulosus, respectively. In both cases, metacestode growth and proliferation occurs within the inner organs of mammalian hosts, which is associated with complex molecular host-parasite interactions that regulate nutrient uptake by the parasite as well as metacestode persistence and development. Using in vitro cultivation systems for parasite larvae, and informed by recently released, comprehensive genome and transcriptome data for both parasites, these molecular host-parasite interactions have been subject to significant research during recent years. In this review, we discuss progress in this field, with emphasis on parasite development and proliferation. We review host-parasite interaction mechanisms that occur early during an infection, when the invading oncosphere stage undergoes a metamorphosis towards the metacestode, and outline the decisive role of parasite stem cells during this process. We also discuss special features of metacestode morphology, and how this parasite stage takes up nutrients from the host, utilizing newly evolved or expanded gene families. We comprehensively review mechanisms of host-parasite cross-communication via evolutionarily conserved signalling systems and how the parasite signalling systems might be exploited for the development of novel chemotherapeutics. Finally, we point to an urgent need for the development of functional genomic techniques in this parasite, which will be imperative for hypothesis-driven analyses into Echinococcus stem cell biology, developmental mechanisms and immunomodulatory activities, which are all highly relevant for the development of anti-infective measures.

Genomes of Fasciola hepatica from the Americas Reveal Colonization with Neorickettsia Endobacteria Related to the Agents of Potomac Horse and Human Sennetsu Fevers

Food borne trematodes (FBTs) are an assemblage of platyhelminth parasites transmitted through the food chain, four of which are recognized as neglected tropical diseases (NTDs). Fascioliasis stands out among the other NTDs due to its broad and significant impact on both human and animal health, as Fasciola sp., are also considered major pathogens of domesticated ruminants. Here we present a reference genome sequence of the common liver fluke, Fasciola hepatica isolated from sheep, complementing previously reported isolate from cattle. A total of 14,642 genes were predicted from the 1.14 GB genome of the liver fluke. Comparative genomics indicated that F. hepatica Oregon and related food-borne trematodes are metabolically less constrained than schistosomes and cestodes, taking advantage of the richer millieux offered by the hepatobiliary organs. Protease families differentially expanded between diverse trematodes may facilitate migration and survival within the heterogeneous environments and niches within the mammalian host. Surprisingly, the sequencing of Oregon and Uruguay F. hepatica isolates led to the first discovery of an endobacteria in this species. Two contigs from the F. hepatica Oregon assembly were joined to complete the 859,205 bp genome of a novel Neorickettsia endobacterium (nFh) closely related to the etiological agents of human Sennetsu and Potomac horse fevers. Immunohistochemical studies targeting a Neorickettsia surface protein found nFh in specific organs and tissues of the adult trematode including the female reproductive tract, eggs, the Mehlis' gland, seminal vesicle, and oral suckers, suggesting putative routes for fluke-to-fluke and fluke-to-host transmission. The genomes of F. hepatica and nFh will serve as a resource for further exploration of the biology of F. hepatica, and specifically its newly discovered trans-kingdom interaction with nFh and the impact of both species on disease in ruminants and humans.

Crystal Structure of Schistosoma mansoni Adenosine Phosphorylase/5'-Methylthioadenosine Phosphorylase and Its Importance on Adenosine Salvage Pathway

Schistosoma mansoni do not have de novo purine pathways and rely on purine salvage for their purine supply. It has been demonstrated that, unlike humans, the S. mansoni is able to produce adenine directly from adenosine, although the enzyme responsible for this activity was unknown. In the present work we show that S. mansoni 5´-deoxy-5´-methylthioadenosine phosphorylase (MTAP, E.C. 2.4.2.28) is capable of use adenosine as a substrate to the production of adenine. Through kinetics assays, we show that the Schistosoma mansoni MTAP (SmMTAP), unlike the mammalian MTAP, uses adenosine substrate with the same efficiency as MTA phosphorolysis, which suggests that this enzyme is part of the purine pathway salvage in S. mansoni and could be a promising target for anti-schistosoma therapies. Here, we present 13 SmMTAP structures from the wild type (WT), including three single and one double mutant, and generate a solid structural framework for structure description. These crystal structures of SmMTAP reveal that the active site contains three substitutions within and near the active site when compared to it mammalian counterpart, thus opening up the possibility of developing specific inhibitors to the parasite MTAP. The structural and kinetic data for 5 substrates reveal the structural basis for this interaction, providing substract for inteligent design of new compounds for block this enzyme activity.

The huge challenge in parasitic chemotherapy development is finding a specific compound to attack the parasite organisms without damaging their host. Schistosoma mansoni, which is the causative agent of schistosomiasis, is one of the major health concerns in the developing world. Purine bases are essential for organisms that make DNA, RNA and energetic molecules during parasitic growth and egg laying. The parasites depend entirely on re-utilizing existing purines, not being able to synthesize them from more simple molecules. The adenosine phosphorylase is an important activity for this process and kinetic assays we performed with this S. mansoni MTAP confirm that it displays this specific activity that is not present in the human metabolism. Therefore, understanding the properties of this enzyme is an important step in achieving an efficient anti-schistosomiasis drug with minimal collateral effects to humans.

WormBase ParaSite - a comprehensive resource for helminth genomics

•

WormBase ParaSite is a new resource for helminth genomics.

•

The resource provides access to over 100 nematode and platyhelminth genomes.

•

The genomes are consistently annotated, organised and presented.

•

A variety of views and tools for exploring and querying the data are provided.

The number of publicly available parasitic worm genome sequences has increased dramatically in the past three years, and research interest in helminth functional genomics is now quickly gathering pace in response to the foundation that has been laid by these collective efforts. A systematic approach to the organisation, curation, analysis and presentation of these data is clearly vital for maximising the utility of these data to researchers. We have developed a portal called WormBase ParaSite (http://parasite.wormbase.org) for interrogating helminth genomes on a large scale. Data from over 100 nematode and platyhelminth species are integrated, adding value by way of systematic and consistent functional annotation (e.g. protein domains and Gene Ontology terms), gene expression analysis (e.g. alignment of life-stage specific transcriptome data sets), and comparative analysis (e.g. orthologues and paralogues). We provide several ways of exploring the data, including genome browsers, genome and gene summary pages, text search, sequence search, a query wizard, bulk downloads, and programmatic interfaces. In this review, we provide an overview of the back-end infrastructure and analysis behind WormBase ParaSite, and the displays and tools available to users for interrogating helminth genomic data.

Curcumin Generates Oxidative Stress and Induces Apoptosis in Adult Schistosoma mansoni Worms

Inducing apoptosis is an interesting therapeutic approach to develop drugs that act against helminthic parasites. Researchers have investigated how curcumin (CUR), a biologically active compound extracted from rhizomes of Curcuma longa, affects Schistosoma mansoni and several cancer cell lines. This study evaluates how CUR influences the induction of apoptosis and oxidative stress in couples of adult S. mansoni worms. CUR decreased the viability of adult worms and killed them. The tegument of the parasite suffered morphological changes, the mitochondria underwent alterations, and chromatin condensed. Different apoptotic parameters were determined in an attempt to understand how CUR affected adult S. mansoni worms. CUR induced DNA damage and fragmentation and increased the expression of SmCASP3/7 transcripts and the activity of Caspase 3 in female and male worms. However, CUR did not intensify the activity of Caspase 8 in female or male worms. Evaluation of the superoxide anion and different antioxidant enzymes helped to explore the mechanism of parasite death further. The level of superoxide anion and the activity of Superoxide Dismutase (SOD) increased, whereas the activity of Glutathione-S-Transferase (GST), Glutathione reductase (GR), and Glutathione peroxidase (GPX) decreased, which culminated in the oxidation of proteins in adult female and male worms incubated with CUR. In conclusion, CUR generated oxidative stress followed by apoptotic-like-events in both adult female and male S. mansoni worms, ultimately killing them.

GPCR and IR genes in Schistosoma mansoni miracidia

Background

Schistosoma species are responsible for the disease schistosomiasis, a highly prevalent helminthic disease that requires a freshwater snail as intermediate host. The S. mansoni free-living miracidium must utilize olfaction to find a suitable snail host, and certain types of rhodopsin G protein-coupled receptors (GPCRs) and ionotropic receptors (IRs) have been identified as olfactory receptors in other animal phyla. The Schistosoma genome project, together with the recent availability of proteomic databases, allowed for studies to explore receptors within S. mansoni, some of which may contribute to host finding.

Results

We have identified 17 rhodopsin-type GPCR sequences in S. mansoni belonging to four subclasses, including ligand-specific GPCRs (i.e. neuropeptide and opsin). RT-PCR demonstrated the expression of nine out of the 17 GPCRs in the free-living miracidia, each of which have been characterized for homology to S. haematobium. Among the nine GPCRs, two are predicted as Gq-opsins. We also describe the characterization of a Schistosoma-encoded IR based on similarity with other species IR and conservation of IR-like domains. Schistosoma mansoni IR is expressed in miracidia at 3 and 6 h post-hatch.

Conclusions

The identification of receptors in S. mansoni miracidia, presented here, contributes not only to further understanding of Schistosoma biology and signal transduction but also provides a basis for approaches that may modify parasite behaviour.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1837-2) contains supplementary material, which is available to authorized users.

The NIH-NIAID Schistosomiasis Resource Center at the Biomedical Research Institute: Molecular Redux

Schistosomiasis remains a health burden in many parts of the world. The complex life cycle of Schistosoma parasites and the economic and societal conditions present in endemic areas make the prospect of eradication unlikely in the foreseeable future. Continued and vigorous research efforts must therefore be directed at this disease, particularly since only a single World Health Organization (WHO)-approved drug is available for treatment. The National Institutes of Health (NIH)–National Institute of Allergy and Infectious Diseases (NIAID) Schistosomiasis Resource Center (SRC) at the Biomedical Research Institute provides investigators with the critical raw materials needed to carry out this important research. The SRC makes available, free of charge (including international shipping costs), not only infected host organisms but also a wide array of molecular reagents derived from all life stages of each of the three main human schistosome parasites. As the field of schistosomiasis research rapidly advances, it is likely to become increasingly reliant on omics, transgenics, epigenetics, and microbiome-related research approaches. The SRC has and will continue to monitor and contribute to advances in the field in order to support these research efforts with an expanding array of molecular reagents. In addition to providing investigators with source materials, the SRC has expanded its educational mission by offering a molecular techniques training course and has recently organized an international schistosomiasis-focused meeting. This review provides an overview of the materials and services that are available at the SRC for schistosomiasis researchers, with a focus on updates that have occurred since the original overview in 2008.

HIV-1 Integrates Widely throughout the Genome of the Human Blood Fluke Schistosoma mansoni

Schistosomiasis is the most important helminthic disease of humanity in terms of morbidity and mortality. Facile manipulation of schistosomes using lentiviruses would enable advances in functional genomics in these and related neglected tropical diseases pathogens including tapeworms, and including their non-dividing cells. Such approaches have hitherto been unavailable. Blood stream forms of the human blood fluke, Schistosoma mansoni, the causative agent of the hepatointestinal schistosomiasis, were infected with the human HIV-1 isolate NL4-3 pseudotyped with vesicular stomatitis virus glycoprotein. The appearance of strong stop and positive strand cDNAs indicated that virions fused to schistosome cells, the nucleocapsid internalized and the RNA genome reverse transcribed. Anchored PCR analysis, sequencing HIV-1-specific anchored Illumina libraries and Whole Genome Sequencing (WGS) of schistosomes confirmed chromosomal integration; >8,000 integrations were mapped, distributed throughout the eight pairs of chromosomes including the sex chromosomes. The rate of integrations in the genome exceeded five per 1,000 kb and HIV-1 integrated into protein-encoding loci and elsewhere with integration bias dissimilar to that of human T cells. We estimated ~ 2,100 integrations per schistosomulum based on WGS, i.e. about two or three events per cell, comparable to integration rates in human cells. Accomplishment in schistosomes of post-entry processes essential for HIV-1replication, including integrase-catalyzed integration, was remarkable given the phylogenetic distance between schistosomes and primates, the natural hosts of the genus Lentivirus. These enigmatic findings revealed that HIV-1 was active within cells of S. mansoni, and provided the first demonstration that HIV-1 can integrate into the genome of an invertebrate.

Sex-Biased Transcriptome of Schistosoma mansoni: Host-Parasite Interaction, Genetic Determinants and Epigenetic Regulators Are Associated with Sexual Differentiation

Background

Among more than 20,000 species of hermaphroditic trematodes, Schistosomatidae are unusual since they have evolved gonochorism. In schistosomes, sex is determined by a female heterogametic system, but phenotypic sexual dimorphism appears only after infection of the vertebrate definitive host. The completion of gonad maturation occurs even later, after pairing. To date, the molecular mechanisms that trigger the sexual differentiation in these species remain unknown, and in vivo studies on the developing schistosomulum stages are lacking. To study the molecular basis of sex determination and sexual differentiation in schistosomes, we investigated the whole transcriptome of the human parasite Schistosoma mansoni in a stage- and sex-comparative manner.

Methodology/ Principal Findings

We performed a RNA-seq on males and females for five developmental stages: cercariae larvae, three in vivo schistosomulum stages and adults. We detected 7,168 genes differentially expressed between sexes in at least one of the developmental stages, and 4,065 of them were functionally annotated. Transcriptome data were completed with H3K27me3 histone modification analysis using ChIP-Seq before (in cercariae) and after (in adults) the phenotypic sexual dimorphism appearance. In this paper we present (i) candidate determinants of the sexual differentiation, (ii) sex-biased players of the interaction with the vertebrate host, and (iii) different dynamic of the H3K27me3 histone mark between sexes as an illustration of sex-biased epigenetic landscapes.

Conclusions/ Significance

Our work presents evidence that sexual differentiation in S. mansoni is accompanied by distinct male and female transcriptional landscapes of known players of the host-parasite crosstalk, genetic determinants and epigenetic regulators. Our results suggest that such combination could lead to the optimized sexual dimorphism of this parasitic species. As S. mansoni is pathogenic for humans, this study represents a promising source of therapeutic targets, providing not only data on the parasite development in interaction with its vertebrate host, but also new insights on its reproductive function.

Parasitic flatworms include more than 20,000 species that are classically hermaphrodites. Among them, the roughly hundred species of Schistosomatidae are intriguing because they are gonochoric. Schistosomes are responsible of the second most important parasitic disease worldwide, and eggs are the main cause of the inflammatory symptoms. Thus, studying the sexual reproduction mechanisms of schistosomes is of particular interest for drug development. Schistosome’s sex is genetically determined by the presence of sex chromosomes: ZZ in males or ZW in females. There is, however, no phenotypic dimorphism in the larval stages: sexual dimorphism appears only in the vertebrate host. In order to understand the molecular mechanisms underlying phenotypic sexual dimorphism, we performed a transcriptome analysis (RNA-Seq) in five different stages of the parasite lifecycle as well as a chromatin status analysis (ChIP-Seq) in the non-differentiated stage cercariae and in the adult differentiated stage, for males and females separately. Our work presents evidence that sexual differentiation in S. mansoni is accompanied by distinct male and female transcriptional landscapes of known players of the host-parasite crosstalk, developmental pathways and epigenetic regulators. Our sex-comparative approach provides therefore new potential therapeutic targets to affect development and sexual reproduction of parasite.

Hsp70 May Be a Molecular Regulator of Schistosome Host Invasion

Schistosomiasis is a debilitating disease that affects over 240 million people worldwide and is considered the most important neglected tropical disease following malaria. Free-swimming freshwater cercariae, one of the six morphologically distinct schistosome life stages, infect humans by directly penetrating through the skin. Cercariae identify and seek the host by sensing chemicals released from human skin. When they reach the host, they burrow into the skin with the help of proteases and other contents released from their acetabular glands and transform into schistosomula, the subsequent larval worm stage upon skin infection. Relative to host invasion, studies have primarily focused on the nature of the acetabular gland secretions, immune response of the host upon exposure to cercariae, and cercaria-schistosomulum transformation methods. However, the molecular signaling pathways involved from host-seeking through the decision to penetrate skin are not well understood. We recently observed that heat shock factor 1 (Hsf1) is localized to the acetabular glands of infectious schistosome cercariae, prompting us to investigate a potential role for heat shock proteins (HSPs) in cercarial invasion. In this study, we report that cercarial invasion behavior, similar to the behavior of cercariae exposed to human skin lipid, is regulated through an Hsp70-dependent process, which we show by using chemical agents that target Hsp70. The observation that biologically active protein activity modulators can elicit a direct and clear behavioral change in parasitic schistosome larvae is itself interesting and has not been previously observed. This finding suggests a novel role for Hsp70 to act as a switch in the cercaria-schistosomulum transformation, and it allows us to begin elucidating the pathways associated with cercarial host invasion. In addition, because the Hsp70 protein and its structure/function is highly conserved, the model that Hsp70 acts as a behavior transitional switch could be relevant to other parasites that also undergo an invasion process and can apply more broadly to other organisms during morphological transitions. Finally, it points to a new function for HSPs in parasite/host interactions.

Parasitic schistosome worms cause morbid disease in over 240 million individuals worldwide. Acute infections with these worms can lead to Katayama fever, while chronic infections can lead to portal hypertension, enlarged abdomen, and liver damage. The infective larval stage, called cercariae, are free-swimming and can detect, seek, and penetrate human skin to enter the human host circulatory system, eventually developing into egg-laying adult worms that cause schistosomiasis. Molecular pathways associated with the initial cercarial invasion of the host, however, are largely unknown, especially with respect to the parasite-specific signals involved in host detection and subsequent decision to invade. Here, we describe a role for Hsp70 in cercarial invasion behavior. To date, only generic stimulation with skin lipid, linoleic acid or L-arginine are known to induce cercarial invasion behavior; thus, we can begin an initial investigation of molecular requirements for host invasion and environment transition for schistosomes and possibly other parasitic organisms.

Evolutionary characterization of Ty3/gypsy-like LTR retrotransposons in the parasitic cestode Echinococcus granulosus

Cyclophyllidean cestodes including Echinococcus granulosus have a smaller genome and show characteristics such as loss of the gut, a segmented body plan, and accelerated growth rate in hosts compared with other tissue-invading helminths. In an effort to address the molecular mechanism relevant to genome shrinkage, the evolutionary status of long-terminal-repeat (LTR) retrotransposons, which are known as the most potent genomic modulators, was investigated in the E. granulosus draft genome. A majority of the E. granulosus LTR retrotransposons were classified into a novel characteristic clade, named Saci-2, of the Ty3/gypsy family, while the remaining elements belonged to the CsRn1 clade of identical family. Their nucleotide sequences were heavily corrupted by frequent base substitutions and segmental losses. The ceased mobile activity of the major retrotransposons and the following intrinsic DNA loss in their inactive progenies might have contributed to decrease in genome size. Apart from the degenerate copies, a gag gene originating from a CsRn1-like element exhibited substantial evidences suggesting its domestication including a preserved coding profile and transcriptional activity, the presence of syntenic orthologues in cestodes, and selective pressure acting on the gene. To my knowledge, the endogenized gag gene is reported for the first time in invertebrates, though its biological function remains elusive.

An Isochore-Like Structure in the Genome of the Flatworm Schistosoma mansoni

Eukaryotic genomes are compositionally heterogeneous, that is, composed by regions that differ in guanine-cytosine (GC) content (isochores). The most well documented case is that of vertebrates (mainly mammals) although it has been also noted among unicellular eukaryotes and invertebrates. In the human genome, regarded as a typical mammal, this heterogeneity is associated with several features. Specifically, genes located in GC-richest regions are the GC3-richest, display CpG islands and have shorter introns. Furthermore, these genes are more heavily expressed and tend to be located at the extremes of the chromosomes. Although the compositional heterogeneity seems to be widespread among eukaryotes, the associated properties noted in the human genome and other mammals have not been investigated in depth in other taxa Here we provide evidence that the genome of the parasitic flatworm Schistosoma mansoni is compositionally heterogeneous and exhibits an isochore-like structure, displaying some features associated, until now, only with the human and other vertebrate genomes, with the exception of gene concentration.

The influence of trematode parasite burden on gene expression in a mammalian host

Background

Parasites can profoundly impact their hosts and are responsible for a plethora of debilitating diseases. To identify global changes in host gene expression related to parasite infection, we sequenced, assembled, and annotated the liver transcriptomes of Balb/cj mice infected with the trematode parasite Schistosoma mansoni and compared the results to uninfected mice. We used two different methodologies (i.e. de novo and reference guided) to evaluate the influence of parasite sequences on host transcriptome assembly.

Results

Our results demonstrate that the choice of assembly methodology significantly impacted the proportion of parasitic reads detected from the host library, yet the presence of non-target (xenobiotic) sequences did not create significant structural errors in the assembly. After removing parasite sequences from the mouse transcriptomes, we analyzed host gene expression under different parasite infection levels and observed significant differences in the associated immunologic and metabolic responses based on infection level. In particular, genes associated with T–helper type 1 (Th–1) and T–helper type 2 (Th–2) were up-regulated in infected mice whereas genes related to amino acid and carbohydrate metabolism were down-regulated in infected mice. These changes in gene expression scale with infection status and likely impact the evolutionary fitness of hosts.

Conclusions

Overall, our data indicate that a) infected mice reduce the expression of key metabolic genes in direct proportion to their infection level; b) infected mice similarly increase the expression of key immune genes in response to infection; c) patterns of gene expression correspond to the pathological symptoms of schistosomiasis; and d) identifying and filtering out non-target sequences (xenobiotics) improves differential expression prediction. Our findings identify parasite targets for RNAi or other therapies and provide a better understanding of the pathology and host immune repertoire involved in response to S. mansoni infections.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2950-5) contains supplementary material, which is available to authorized users.

Discovery of New Anti-Schistosomal Hits by Integration of QSAR-Based Virtual Screening and High Content Screening

Schistosomiasis is a debilitating neglected tropical disease, caused by flatworms of Schistosoma genus. The treatment relies on a single drug, praziquantel (PZQ), making the discovery of new compounds extremely urgent. In this work, we integrated QSAR-based virtual screening (VS) of Schistosoma mansoni thioredoxin glutathione reductase (SmTGR) inhibitors and high content screening (HCS) aiming to discover new antischistosomal agents. Initially, binary QSAR models for inhibition of SmTGR were developed and validated using the Organization for Economic Co-operation and Development (OECD) guidance. Using these models, we prioritized 29 compounds for further testing in two HCS platforms based on image analysis of assay plates. Among them, 2-[2-(3-methyl-4-nitro-5-isoxazolyl)vinyl]pyridine and 2-(benzylsulfonyl)-1,3-benzothiazole, two compounds representing new chemical scaffolds have activity against schistosomula and adult worms at low micromolar concentrations and therefore represent promising antischistosomal hits for further hit-to-lead optimization.

TRP channels in schistosomes

Praziquantel (PZQ) is effectively the only drug currently available for treatment and control of schistosomiasis, a disease affecting hundreds of millions of people worldwide. Many anthelmintics, likely including PZQ, target ion channels, membrane protein complexes essential for normal functioning of the neuromusculature and other tissues. Despite this fact, only a few classes of parasitic helminth ion channels have been assessed for their pharmacological properties or for their roles in parasite physiology. One such overlooked group of ion channels is the transient receptor potential (TRP) channel superfamily. TRP channels share a common core structure, but are widely diverse in their activation mechanisms and ion selectivity. They are critical to transducing sensory signals, responding to a wide range of external stimuli. They are also involved in other functions, such as regulating intracellular calcium and organellar ion homeostasis and trafficking. Here, we review current literature on parasitic helminth TRP channels, focusing on those in schistosomes. We discuss the likely roles of these channels in sensory and locomotor activity, including the possible significance of a class of TRP channels (TRPV) that is absent in schistosomes. We also focus on evidence indicating that at least one schistosome TRP channel (SmTRPA) has atypical, TRPV1-like pharmacological sensitivities that could potentially be exploited for future therapeutic targeting.

•

We provide an overview of transient receptor potential (TRP) channels in schistosomes and other parasitic helminths.

•

TRP channels are important for sensory signaling, ion homeostasis, organellar trafficking, and a host of other functions.

•

Very little work has been done on TRP channels in parasitic helminths.

•

TRPV channels, found throughout the Metazoa, appear not to be present in parasitic platyhelminths.

•

TRP channels in schistosomes appear to have atypical pharmacology, perhaps an entrée for therapeutic targeting.

Uncovering Notch pathway in the parasitic flatworm Schistosoma mansoni

Several signaling molecules that govern development in higher animals have been identified in the parasite Schistosoma mansoni, including the transforming growth factor β, protein tyrosine kinases, nuclear hormone receptors, among others. The Notch pathway is a highly conserved signaling mechanism which is involved in a wide variety of developmental processes including embryogenesis and oogenesis in worms and flies. Here we aimed to provide the molecular reconstitution of the Notch pathway in S. mansoni using the available transcriptome and genome databases. Our results also revealed the presence of the transcripts coded for SmNotch, SmSu(H), SmHes, and the gamma-secretase complex (SmNicastrin, SmAph-1, and SmPen-2), throughout all the life stages analyzed. Besides, it was observed that the viability and separation of adult worm pairs were not affected by treatment with N-[N(3,5)-difluorophenacetyl)-L-Alanyl]-S-phenylglycine t-butyl ester (DAPT), a Notch pathway inhibitor. Moreover, DAPT treatment decreased the production of phenotypically normal eggs and arrested their development in culture. Our results also showed a significant decrease in SmHes transcript levels in both adult worms and eggs treated with DAPT. These results provide, for the first time, functional validation of the Notch pathway in S. mansoni and suggest its involvement in parasite oogenesis and embryogenesis. Given the complexity of the Notch pathway, further experiments shall highlight the full repertoire of Notch-mediated cellular processes throughout the S. mansoni life cycle.

Pharmacological profiling an abundantly expressed schistosome serotonergic GPCR identifies nuciferine as a potent antagonist

5-hydroxytryptamine (5-HT) is a key regulator of muscle contraction in parasitic flatworms. In Schistosoma mansoni, the myoexcitatory action of 5-HT is effected through activation of a serotonergic GPCR (Sm.5HTR_L), prioritizing pharmacological characterization of this target for anthelmintic drug discovery. Here, we have examined the effects of several aporphine alkaloids on the signaling activity of a heterologously expressed Sm.5HTR_L construct using a cAMP biosensor assay. Four structurally related natural products - nuciferine, D-glaucine, boldine and bulbocapnine - were demonstrated to block Sm.5HTR_L evoked cAMP generation with the potency of GPCR blockade correlating well with the ability of each drug to inhibit contractility of schistosomule larvae. Nuciferine was also effective at inhibiting both basal and 5-HT evoked motility of adult schistosomes. These data advance our understanding of structure-affinity relationships at Sm.5HTR_L, and demonstrate the effectiveness of Sm.5HTR_L antagonists as hypomotility-evoking drugs across different parasite life cycle stages.

QSAR-Driven Discovery of Novel Chemical Scaffolds Active against Schistosoma mansoni

Schistosomiasis is a neglected tropical disease that affects millions of people worldwide. Thioredoxin glutathione reductase of Schistosoma mansoni (SmTGR) is a validated drug target that plays a crucial role in the redox homeostasis of the parasite. We report the discovery of new chemical scaffolds against S. mansoni using a combi-QSAR approach followed by virtual screening of a commercial database and confirmation of top ranking compounds by in vitro experimental evaluation with automated imaging of schistosomula and adult worms. We constructed 2D and 3D quantitative structure-activity relationship (QSAR) models using a series of oxadiazoles-2-oxides reported in the literature as SmTGR inhibitors and combined the best models in a consensus QSAR model. This model was used for a virtual screening of Hit2Lead set of ChemBridge database and allowed the identification of ten new potential SmTGR inhibitors. Further experimental testing on both shistosomula and adult worms showed that 4-nitro-3,5-bis(1-nitro-1H-pyrazol-4-yl)-1H-pyrazole (LabMol-17) and 3-nitro-4-{[(4-nitro-1,2,5-oxadiazol-3-yl)oxy]methyl}-1,2,5-oxadiazole (LabMol-19), two compounds representing new chemical scaffolds, have high activity in both systems. These compounds will be the subjects for additional testing and, if necessary, modification to serve as new schistosomicidal agents.

A Miniaturized Screen of a Schistosoma mansoni Serotonergic G Protein-Coupled Receptor Identifies Novel Classes of Parasite-Selective Inhibitors

Schistosomiasis is a tropical parasitic disease afflicting ~200 million people worldwide and current therapy depends on a single drug (praziquantel) which exhibits several non-optimal features. These shortcomings underpin the need for next generation anthelmintics, but the process of validating physiologically relevant targets (‘target selection’) and pharmacologically profiling them is challenging. Remarkably, even though over a quarter of current human therapeutics target rhodopsin-like G protein coupled receptors (GPCRs), no library screen of a flatworm GPCR has yet been reported. Here, we have pharmacologically profiled a schistosome serotonergic GPCR (Sm.5HTR) implicated as a downstream modulator of PZQ efficacy, in a miniaturized screening assay compatible with high content screening. This approach employs a split luciferase based biosensor sensitive to cellular cAMP levels that resolves the proximal kinetics of GPCR modulation in intact cells. Data evidence a divergent pharmacological signature between the parasitic serotonergic receptor and the closest human GPCR homolog (Hs.5HTR7), supporting the feasibility of optimizing parasitic selective pharmacophores. New ligands, and chemical series, with potency and selectivity for Sm.5HTR over Hs.5HTR7 are identified in vitro and validated for in vivo efficacy against schistosomules and adult worms. Sm.5HTR also displayed a property resembling irreversible inactivation, a phenomenon discovered at Hs.5HTR7, which enhances the appeal of this abundantly expressed parasite GPCR as a target for anthelmintic ligand design. Overall, these data underscore the feasibility of profiling flatworm GPCRs in a high throughput screening format competent to resolve different classes of GPCR modulators. Further, these data underscore the promise of Sm.5HTR as a chemotherapeutically vulnerable node for development of next generation anthelmintics.

Parasitic flatworms express a diverse array of G protein coupled receptors, but our knowledge of their pharmacological profile is limited. No high throughput screen of a flatworm GPCR has been reported, even though these targets have precedent for high druggability and functionality in the chemotherapeutically vulnerable excitable cell niche. The goal of this study was to establish a method for profiling flatworm G protein coupled receptors that can be scaled to high content screening. Using a cAMP biosensor, we have performed a proof of principle miniaturized screen on a schistosome serotonergic GPCR that resolves new ligands that potently and selectivity block 5-HT receptor activity in vitro, and 5-HT evoked responses in schistosomules and adult worms. This approach evidences the pharmacological divergence of a parasitic GPCR from the closest human homolog and a capacity for high content interrogation of flatworm GPCR properties and ligand specificities.

Identification of Candidate Serum Biomarkers for Schistosoma mansoni Infected Mice Using Multiple Proteomic Platforms

BACKGROUND

Schistosomiasis is an important helminth infection of humans. There are few reliable diagnostic biomarkers for early infection, for recurrent infection or to document successful treatment. In this study, we compared serum protein profiles in uninfected and infected mice to identify disease stage-specific biomarkers.

METHODS

Serum collected from CD1 mice infected with 50-200 Schistosoma mansoni cercariae were analyzed before infection and at 3, 6 and 12 weeks post-infection using three mass spectrometric (MS) platforms.

RESULTS

Using SELDI-TOF MS, 66 discriminating m/z peaks were detected between S. mansoni infected mice and healthy controls. Used in various combinations, these peaks could 1) reliably diagnose early-stage disease, 2) distinguish between acute and chronic infection and 3) diagnose S. mansoni infection regardless the parasite burden. The most important contributors to these diagnostic algorithms were peaks at 3.7, 13 and 46 kDa. Employing sample fractionation and differential gel electrophoresis, we analyzed gel slices either by MALDI-TOF MS or Velos Orbitrap MS. The former yielded eight differentially-expressed host proteins in the serum at different disease stages including transferrin and alpha 1- antitrypsin. The latter suggested the presence of a surprising number of parasite-origin proteins in the serum during both the acute (n = 200) and chronic (n = 105) stages. The Orbitrap platform also identified many differentially-expressed host-origin serum proteins during the acute and chronic stages (296 and 220 respectively). The presence of one of the schistosome proteins, glutathione S transferase (GST: 25 KDa), was confirmed by Western Blot. This study provides proof-of-principle for an approach that can yield a large number of novel candidate biomarkers for Schistosoma infection.

Modern approaches to accelerate discovery of new antischistosomal drugs

INTRODUCTION

The almost exclusive use of only praziquantel for the treatment of schistosomiasis has raised concerns about the possible emergence of drug-resistant schistosomes. Consequently, there is an urgent need for new antischistosomal drugs. The identification of leads and the generation of high quality data are crucial steps in the early stages of schistosome drug discovery projects.

AREAS COVERED

Herein, the authors focus on the current developments in antischistosomal lead discovery, specifically referring to the use of automated in vitro target-based and whole-organism screens and virtual screening of chemical databases. They highlight the strengths and pitfalls of each of the above-mentioned approaches, and suggest possible roadmaps towards the integration of several strategies, which may contribute for optimizing research outputs and led to more successful and cost-effective drug discovery endeavors.

EXPERT OPINION

Increasing partnerships and access to funding for drug discovery have strengthened the battle against schistosomiasis in recent years. However, the authors believe this battle also includes innovative strategies to overcome scientific challenges. In this context, significant advances of in vitro screening as well as computer-aided drug discovery have contributed to increase the success rate and reduce the costs of drug discovery campaigns. Although some of these approaches were already used in current antischistosomal lead discovery pipelines, the integration of these strategies in a solid workflow should allow the production of new treatments for schistosomiasis in the near future.

A next-generation proteome array for Schistosoma mansoni

A proteome microarray consisting of 992 Schistosoma mansoni proteins was produced and screened with sera to determine antibody signatures indicative of the clinical stages of schistosomiasis and the identification of subunit vaccine candidates. Herein, we describe the methods used to derive the gene list for this array (representing approximately 10% of the predicted S. mansoni proteome). We also probed a pilot version of the microarray with sera from individuals either acutely or chronically infected with S. mansoni from endemic areas in Brazil and sera from individuals resident outside the endemic area (USA) to determine if the array is functional and informative.

Haemonchus contortus: Genome Structure, Organization and Comparative Genomics

One of the first genome sequencing projects for a parasitic nematode was that for Haemonchus contortus. The open access data from the Wellcome Trust Sanger Institute provided a valuable early resource for the research community, particularly for the identification of specific genes and genetic markers. Later, a second sequencing project was initiated by the University of Melbourne, and the two draft genome sequences for H. contortus were published back-to-back in 2013. There is a pressing need for long-range genomic information for genetic mapping, population genetics and functional genomic studies, so we are continuing to improve the Wellcome Trust Sanger Institute assembly to provide a finished reference genome for H. contortus. This review describes this process, compares the H. contortus genome assemblies with draft genomes from other members of the strongylid group and discusses future directions for parasite genomics using the H. contortus model.

Identification of new markers for the Schistosoma mansoni vitelline lineage

•

Transcriptional profiling identifies candidate factors associated with the schistosome vitellarium.

•

In situ hybridization confirms many new markers of this tissue.

•

New cell type-specific markers for various stages of vitellocyte development are reported.

Schistosomes cause significant morbidity and mortality in millions of the world’s poorest people. While parasite egg-induced inflammation is the primary driver of host pathology, relatively little is known at the molecular level about the organ systems that participate in schistosome egg production (i.e., testes, ovaries and vitellaria). Here we use transcriptional profiling and in situ hybridization to characterise the vitellarium of Schistosoma mansoni. We uncovered several previously uncharacterised vitellaria-specific factors and defined molecular markers for various stages in the vitellocyte differentiation process. These data provide the framework for future in-depth molecular studies exploring the biology of this important parasite organ.

Epithelial sodium channel (ENaC) family: Phylogeny, structure-function, tissue distribution, and associated inherited diseases

The epithelial sodium channel (ENaC) is composed of three homologous subunits and allows the flow of Na(+) ions across high resistance epithelia, maintaining body salt and water homeostasis. ENaC dependent reabsorption of Na(+) in the kidney tubules regulates extracellular fluid (ECF) volume and blood pressure by modulating osmolarity. In multi-ciliated cells, ENaC is located in cilia and plays an essential role in the regulation of epithelial surface liquid volume necessary for cilial transport of mucus and gametes in the respiratory and reproductive tracts respectively. The subunits that form ENaC (named as alpha, beta, gamma and delta, encoded by genes SCNN1A, SCNN1B, SCNN1G, and SCNN1D) are members of the ENaC/Degenerin superfamily. The earliest appearance of ENaC orthologs is in the genomes of the most ancient vertebrate taxon, Cyclostomata (jawless vertebrates) including lampreys, followed by earliest representatives of Gnathostomata (jawed vertebrates) including cartilaginous sharks. Among Euteleostomi (bony vertebrates), Actinopterygii (ray finned-fishes) branch has lost ENaC genes. Yet, most animals in the Sarcopterygii (lobe-finned fish) branch including Tetrapoda, amphibians and amniotes (lizards, crocodiles, birds, and mammals), have four ENaC paralogs. We compared the sequences of ENaC orthologs from 20 species and established criteria for the identification of ENaC orthologs and paralogs, and their distinction from other members of the ENaC/Degenerin superfamily, especially ASIC family. Differences between ENaCs and ASICs are summarized in view of their physiological functions and tissue distributions. Structural motifs that are conserved throughout vertebrate ENaCs are highlighted. We also present a comparative overview of the genotype-phenotype relationships in inherited diseases associated with ENaC mutations, including multisystem pseudohypoaldosteronism (PHA1B), Liddle syndrome, cystic fibrosis-like disease and essential hypertension.

Pipeline for the identification and classification of ion channels in parasitic flatworms

Background

Ion channels are well characterised in model organisms, principally because of the availability of functional genomic tools and datasets for these species. This contrasts the situation, for example, for parasites of humans and animals, whose genomic and biological uniqueness means that many genes and their products cannot be annotated. As ion channels are recognised as important drug targets in mammals, the accurate identification and classification of parasite channels could provide major prospects for defining unique targets for designing novel and specific anti-parasite therapies. Here, we established a reliable bioinformatic pipeline for the identification and classification of ion channels encoded in the genome of the cancer-causing liver fluke Opisthorchis viverrini, and extended its application to related flatworms affecting humans.

Methods

We built an ion channel identification + classification pipeline (called MuSICC), employing an optimised support vector machine (SVM) model and using the Kyoto Encyclopaedia of Genes and Genomes (KEGG) classification system. Ion channel proteins were first identified and grouped according to amino acid sequence similarity to classified ion channels and the presence and number of ion channel-like conserved and transmembrane domains. Predicted ion channels were then classified to sub-family using a SVM model, trained using ion channel features.

Results

Following an evaluation of this pipeline (MuSICC), which demonstrated a classification sensitivity of 95.2 % and accuracy of 70.5 % for known ion channels, we applied it to effectively identify and classify ion channels in selected parasitic flatworms.

Conclusions

MuSICC provides a practical and effective tool for the identification and classification of ion channels of parasitic flatworms, and should be applicable to a broad range of organisms that are evolutionarily distant from taxa whose ion channels are functionally characterised.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1428-2) contains supplementary material, which is available to authorized users.