WormBase ParaSite - a comprehensive resource for helminth genomics

Trematode Genomics and Proteomics

Trematode infections are among the most neglected tropical diseases despite their worldwide distribution and extraordinary ability to parasitise many different host species and host tissues. Furthermore, these parasites are of great socioeconomic, medical, veterinary and agricultural importance. During the last 10 years, there have been increasing efforts to overcome the lack of information on different "omic" resources such as proteomics and genomics. Herein, we focus on the recent advances in genomics and proteomics from trematodes of human importance, including liver, blood, intestinal and lung flukes. We also provide information on the latest technologies applied to study the biology of trematodes as well as on the resources available for the study of the molecular aspects of this group of helminths.

SCNBase: a genomics portal for the soybean cyst nematode (Heterodera glycines)

Soybean is an important worldwide crop, and farmers continue to experience significant yield loss due to the soybean cyst nematode (SCN), Heterodera glycines. This soil-borne roundworm parasite is rated the most important pathogen problem in soybean production. The infective nematodes enter into complex interactions with their host plant by inducing the development of specialized plant feeding cells that provide the parasites with nourishment. Addressing the SCN problem will require the development of genomic resources and a global collaboration of scientists to analyze and use these resources. SCNBase.org was designed as a collaborative hub for the SCN genome. All data and analyses are downloadable and can be analyzed with three integrated genomic tools: JBrowse, Feature Search and BLAST. At the time of this writing, a number of genomic and transcriptomic data sets are already available, with 43 JBrowse tracks and 21 category pages describing SCN genomic analyses on gene predictions, transcriptome and read alignments, effector-like genes, expansion and contraction of genomic repeats, orthology and synteny with related nematode species, Single Nucleotide Polymorphism (SNPs) from 15 SCN populations and novel splice sites. Standard functional gene annotations were supplemented with orthologous gene annotations using a comparison to nine related plant-parasitic nematodes, thereby enabling functional annotations for 85% of genes. These annotations led to a greater grasp on the SCN effectorome, which include over 3324 putative effector genes. By designing SCNBase as a hub, future research findings and genomic resources can easily be uploaded and made available for use by others with minimal needs for further curation. By providing these resources to nematode research community, scientists will be empowered to develop novel, more effective SCN management tools.

Comparative genomics of the major parasitic worms

Parasitic nematodes (roundworms) and platyhelminths (flatworms) cause debilitating chronic infections of humans and animals, decimate crop production and are a major impediment to socioeconomic development. Here we report a broad comparative study of 81 genomes of parasitic and non-parasitic worms. We have identified gene family births and hundreds of expanded gene families at key nodes in the phylogeny that are relevant to parasitism. Examples include gene families that modulate host immune responses, enable parasite migration though host tissues or allow the parasite to feed. We reveal extensive lineage-specific differences in core metabolism and protein families historically targeted for drug development. From an in silico screen, we have identified and prioritized new potential drug targets and compounds for testing. This comparative genomics resource provides a much-needed boost for the research community to understand and combat parasitic worms.

Genome aware CRISPR gRNA target prediction for parasitic nematodes

The pace of research towards a genetic model to understand the unique molecular biology of parasitic nematodes has increased recently. This research has developed a diverse suite of genetic tools for a variety of parasitic nematodes. CRISPR/Cas9 technology in particular offers much promise as a game changing tool for researchers studying parasitic nematodes. Unlike RNAi, which depends on diverse nematode effectors to silence gene expression, the effectors for CRISPR/Cas9 mutations are typically supplied by the experimenter, making gene editing via CRISPR/Cas9 ideal for testing on genetically intractable nematode systems. To facilitate the development of CRISPR/Cas9 technology for parasitic nematodes, I here describe a tool for identifying gRNA targets and diagnostic primers to a user supplied sequence. The software attempts to minimize non-specific targets by interrogating the genomes of parasitic nematodes. This software is freely available online and features an intuitive interface to help researchers design effective CRISPR experiments for parasitic nematodes.

Identification and partial characterization of a novel serpin from Eudiplozoon nipponicum (Monogenea, Polyopisthocotylea)

Background: Serpins are a superfamily of serine peptidase inhibitors that participate in the regulation of many physiological and cell peptidase-mediated processes in all organisms (e.g. in blood clotting, complement activation, fibrinolysis, inflammation, and programmed cell death). It was postulated that in the blood-feeding members of the monogenean family Diplozoidae, serpins could play an important role in the prevention of thrombus formation, activation of complement, inflammation in the host, and/or in the endogenous regulation of protein degradation.

Results: In silico analysis showed that the DNA and primary protein structures of serpin from Eudiplozoon nipponicum (EnSerp1) are similar to other members of the serpin superfamily. The inhibitory potential of EnSerp1 on four physiologically-relevant serine peptidases (trypsin, factor Xa, kallikrein, and plasmin) was demonstrated and its presence in the worm’s excretory-secretory products (ESPs) was confirmed.

Conclusion: EnSerp1 influences the activity of peptidases that play a role in blood coagulation, fibrinolysis, and complement activation. This inhibitory potential, together with the serpin’s presence in ESPs, suggests that it is likely involved in host-parasite interactions and could be one of the molecules involved in the control of feeding and prevention of inflammatory responses.

Terror in the dirt: Sensory determinants of host seeking in soil-transmitted mammalian-parasitic nematodes

Infection with gastrointestinal parasitic nematodes is a major cause of chronic morbidity and economic burden around the world, particularly in low-resource settings. Some parasitic nematode species, including the human-parasitic threadworm Strongyloides stercoralis and human-parasitic hookworms in the genera Ancylostoma and Necator, feature a soil-dwelling infective larval stage that seeks out hosts for infection using a variety of host-emitted sensory cues. Here, we review our current understanding of the behavioral responses of soil-dwelling infective larvae to host-emitted sensory cues, and the molecular and cellular mechanisms that mediate these responses. We also discuss the development of methods for transgenesis and CRISPR/Cas9-mediated targeted mutagenesis in Strongyloides stercoralis and the closely related rat parasite Strongyloides ratti. These methods have established S. stercoralis and S. ratti as genetic model systems for gastrointestinal parasitic nematodes and are enabling more detailed investigations into the neural mechanisms that underlie the sensory-driven behaviors of this medically and economically important class of parasites.

Temperature-dependent behaviors of parasitic helminths

Parasitic helminth infections are the most common source of neglected tropical disease among impoverished global communities. Many helminths infect their hosts via an active, sensory-driven process in which environmentally motile infective larvae position themselves near potential hosts. For these helminths, host seeking and host invasion can be divided into several discrete behaviors that are regulated by both host-emitted and environmental sensory cues, including heat. Thermosensation is a critical sensory modality for helminths that infect warm-blooded hosts, driving multiple behaviors necessary for host seeking and host invasion. Furthermore, thermosensory cues influence the host-seeking behaviors of both helminths that parasitize endothermic hosts and helminths that parasitize insect hosts. Here, we discuss the role of thermosensation in guiding the host-seeking and host-infection behaviors of a diverse group of helminths, including mammalian-parasitic nematodes, entomopathogenic nematodes, and schistosomes. We also discuss the neural circuitry and molecular pathways that underlie thermosensory responses in these species.

Combining bioinformatics, cheminformatics, functional genomics and whole organism approaches for identifying epigenetic drug targets in Schistosoma mansoni

Schistosomiasis endangers the lives of greater than 200 million people every year and is predominantly controlled by a single class chemotherapy, praziquantel (PZQ). Development of PZQ replacement (to combat the threat of PZQ insensitivity/resistance arising) or combinatorial (to facilitate the killing of PZQ-insensitive juvenile schistosomes) chemotherapies would help sustain this control strategy into the future. Here, we re-categorise two families of druggable epigenetic targets in Schistosoma mansoni, the histone methyltransferases (HMTs) and the histone demethylases (HDMs). Amongst these, a S. mansoni Lysine Specific Demethylase 1 (SmLSD1, Smp_150560) homolog was selected for further analyses. Homology modelling of SmLSD1 and in silico docking of greater than four thousand putative inhibitors identified seven (L1 – L7) showing more favourable binding to the target pocket of SmLSD1 vs Homo sapiens HsLSD1; six of these seven (L1 – L6) plus three structural analogues of L7 (L8 – L10) were subsequently screened against schistosomula using the Roboworm anthelmintic discovery platform. The most active compounds (L10 - pirarubicin > L8 – danunorubicin hydrochloride) were subsequently tested against juvenile (3 wk old) and mature (7 wk old) schistosome stages and found to impede motility, suppress egg production and affect tegumental surfaces. When compared to a surrogate human cell line (HepG2), a moderate window of selectivity was observed for the most active compound L10 (selectivity indices - 11 for schistosomula, 9 for juveniles, 1.5 for adults). Finally, RNA interference of SmLSD1 recapitulated the egg-laying defect of schistosomes co-cultivated in the presence of L10 and L8. These preliminary results suggest that SmLSD1 represents an attractive new target for schistosomiasis; identification of more potent and selective SmLSD1 compounds, however, is essential. Nevertheless, the approaches described herein highlight an interdisciplinary strategy for selecting and screening novel/repositioned anti-schistosomals, which can be applied to any druggable (epigenetic) target encoded by the parasite's genome.

•

Schistosoma mansoni contains 27 histone methyltransferases (HMTs) and 14 histone demethylases (HDMs).

•

S. mansoni lysine specific demethylase 1 (SmLSD1) is a druggable target.

•

Schistosomes treated with the putative SmLSD1 inhibitor pirarubicin or siRNAs targeting SmLSD1 are less fecund.

Genome-wide transcriptome profiling and spatial expression analyses identify signals and switches of development in tapeworms

BACKGROUND

Tapeworms are agents of neglected tropical diseases responsible for significant health problems and economic loss. They also exhibit adaptations to a parasitic lifestyle that confound comparisons of their development with other animals. Identifying the genetic factors regulating their complex ontogeny is essential to understanding unique aspects of their biology and for advancing novel therapeutics. Here we use RNA sequencing to identify up-regulated signalling components, transcription factors and post-transcriptional/translational regulators (genes of interest, GOI) in the transcriptomes of Larvae and different regions of segmented worms in the tapeworm Hymenolepis microstoma and combine this with spatial gene expression analyses of a selection of genes.

RESULTS

RNA-seq reads collectively mapped to 90% of the > 12,000 gene models in the H. microstoma v.2 genome assembly, demonstrating that the transcriptome profiles captured a high percentage of predicted genes. Contrasts made between the transcriptomes of Larvae and whole, adult worms, and between the Scolex-Neck, mature strobila and gravid strobila, resulted in 4.5-30% of the genes determined to be differentially expressed. Among these, we identified 190 unique GOI up-regulated in one or more contrasts, including a large range of zinc finger, homeobox and other transcription factors, components of Wnt, Notch, Hedgehog and TGF-β/BMP signalling, and post-transcriptional regulators (e.g. Boule, Pumilio). Heatmap clusterings based on overall expression and on select groups of genes representing 'signals' and 'switches' showed that expression in the Scolex-Neck region is more similar to that of Larvae than to the mature or gravid regions of the adult worm, which was further reflected in large overlap of up-regulated GOI.

CONCLUSIONS

Spatial expression analyses in Larvae and adult worms corroborated inferences made from quantitative RNA-seq data and in most cases indicated consistency with canonical roles of the genes in other animals, including free-living flatworms. Recapitulation of developmental factors up-regulated during larval metamorphosis suggests that strobilar growth involves many of the same underlying gene regulatory networks despite the significant disparity in developmental outcomes. The majority of genes identified were investigated in tapeworms for the first time, setting the stage for advancing our understanding of developmental genetics in an important group of flatworm parasites.

Preliminary evaluation of neoblast-like stem cell factor and transcript expression profiles in Schistosoma japonicum

Neoblast-like stem cell factors and transcripts are essential for cell proliferation, self-renewal, and differentiation. Recent studies have demonstrated that nanos, sox, and vasa-like transcription factors are associated with neoblast-like stem cells in Schistosoma mansoni and play crucial roles in the regulation of worm development. However, these neoblast-like stem cell factors and transcripts and their expression profiles remain unknown in Schistosoma japonicum. In this study, we identified orthologs of 11 neoblast-like stem cell factors and transcripts in S. japonicum using bioinformatics and confirmed them by PCR. The expression profiles of neoblast-like stem cell factors and transcripts revealed that some of them were highly expressed in certain stages. Sex-based expression analysis revealed that nanos, polo-like kinase, PCNA, cyclin B, and H2A showed significantly higher expression in female worms, whereas ago and bruli showed higher expression in male worms. In addition, we noted that ago, bruli, and pp32 exhibited higher expression in the testes, while nanos, polo-like kinase, cyclin B, H2A, and H2B showed notable higher expression in both isolated ovaries and testes. Our preliminary results are expected to provide important information about the regulatory roles of these stem cell factors in parasite development and sexual maturation.

Extreme allelic heterogeneity at a Caenorhabditis elegans beta-tubulin locus explains natural resistance to benzimidazoles

Benzimidazoles (BZ) are essential components of the limited chemotherapeutic arsenal available to control the global burden of parasitic nematodes. The emerging threat of BZ resistance among multiple nematode species necessitates the development of novel strategies to identify genetic and molecular mechanisms underlying this resistance. All detection of parasitic helminth resistance to BZ is focused on the genotyping of three variant sites in the orthologs of the β-tubulin gene found to confer resistance in the free-living nematode Caenorhabditis elegans. Because of the limitations of laboratory and field experiments in parasitic nematodes, it is difficult to look beyond these three sites to identify additional mechanisms that might contribute to BZ resistance in the field. Here, we took an unbiased genome-wide mapping approach in the free-living nematode species C. elegans to identify the genetic underpinnings of natural resistance to the commonly used BZ, albendazole (ABZ). We found a wide range of natural variation in ABZ resistance in natural C. elegans populations. In agreement with known mechanisms of BZ resistance in parasites, we found that a majority of the variation in ABZ resistance among wild C. elegans strains is caused by variation in the β-tubulin gene ben-1. This result shows empirically that resistance to ABZ naturally exists and segregates within the C. elegans population, suggesting that selection in natural niches could enrich for resistant alleles. We identified 25 distinct ben-1 alleles that are segregating at low frequencies within the C. elegans population, including many novel molecular variants. Population genetic analyses indicate that ben-1 variation arose multiple times during the evolutionary history of C. elegans and provide evidence that these alleles likely occurred recently because of local selective pressures. Additionally, we find purifying selection at all five β-tubulin genes, despite predicted loss-of-function variants in ben-1, indicating that BZ resistance in natural niches is a stronger selective pressure than loss of one β-tubulin gene. Furthermore, we used genome-editing to show that the most common parasitic nematode β-tubulin allele that confers BZ resistance, F200Y, confers resistance in C. elegans. Importantly, we identified a novel genomic region that is correlated with ABZ resistance in the C. elegans population but independent of ben-1 and the other β-tubulin loci, suggesting that there are multiple mechanisms underlying BZ resistance. Taken together, our results establish a population-level resource of nematode natural diversity as an important model for the study of mechanisms that give rise to BZ resistance.

Schistosoma "Eggs-Iting" the Host: Granuloma Formation and Egg Excretion

Schistosomiasis is a major cause of morbidity in humans invoked by chronic infection with parasitic trematodes of the genus Schistosoma. Schistosomes have a complex life-cycle involving infections of an aquatic snail intermediate host and a definitive mammalian host. In humans, adult male and female worms lie within the vasculature. Here, they propagate and eggs are laid. These eggs must then be released from the host to continue the life cycle. Schistosoma mansoni and Schistosoma japonicum reside in the mesenteric circulation of the intestines with egg excreted in the feces. In contrast, S. haematobium are present in the venus plexus of the bladder, expelling eggs in the urine. In an impressive case of exploitation of the host immune system, this process of Schistosome “eggs-iting” the host is immune dependent. In this article, we review the formation of the egg granuloma and explore how S. mansoni eggs laid in vasculature must usurp immunity to induce regulated inflammation, to facilitate extravasation through the intestinal wall and to be expelled in the feces. We highlight the roles of immune cell populations, stromal factors, and egg secretions in the process of egg excretion to provide a comprehensive overview of the current state of knowledge regarding a vastly unexplored mechanism.

Host- and Helminth-Derived Endocannabinoids That Have Effects on Host Immunity Are Generated during Infection

Helminths have coevolved with their hosts, resulting in the development of specialized host immune mechanisms and parasite-specific regulatory products. Identification of new pathways that regulate helminth infection could provide a better understanding of host-helminth interaction and may identify new therapeutic targets for helminth infection. Here we identify the endocannabinoid system as a new mechanism that influences host immunity to helminths. Endocannabinoids are lipid-derived signaling molecules that control important physiologic processes, such as feeding behavior and metabolism. Following murine infection with Nippostrongylus brasiliensis, an intestinal nematode with a life cycle similar to that of hookworms, we observed increased levels of endocannabinoids (2-arachidonoylglycerol [2-AG] or anandamide [AEA]) and the endocannabinoid-like molecule oleoylethanolamine (OEA) in infected lung and intestine. To investigate endocannabinoid function in helminth infection, we employed pharmacological inhibitors of cannabinoid subtype receptors 1 and 2 (CB₁R and CB₂R). Compared to findings for vehicle-treated mice, inhibition of CB₁R but not CB₂R resulted in increased N. brasiliensis worm burden and egg output, associated with significantly decreased expression of the T helper type 2 cytokine interleukin 5 (IL-5) in intestinal tissue and splenocyte cultures. Strikingly, bioinformatic analysis of genomic and transcriptome sequencing (RNA-seq) data sets identified putative genes encoding endocannabinoid biosynthetic and degradative enzymes in many parasitic nematodes. To test the novel hypothesis that helminth parasites produce their own endocannabinoids, we measured endocannabinoid levels in N. brasiliensis by mass spectrometry and quantitative PCR and found that N. brasiliensis parasites produced endocannabinoids, especially at the infectious larval stage. To our knowledge, this is the first report of helminth- and host-derived endocannabinoids that promote host immune responses and reduce parasite burden.

Virus discovery in all three major lineages of terrestrial arthropods highlights the diversity of single-stranded DNA viruses associated with invertebrates

Viruses encoding a replication-associated protein (Rep) within a covalently closed, single-stranded (ss)DNA genome are among the smallest viruses known to infect eukaryotic organisms, including economically valuable agricultural crops and livestock. Although circular Rep-encoding ssDNA (CRESS DNA) viruses are a widespread group for which our knowledge is rapidly expanding, biased sampling toward vertebrates and land plants has limited our understanding of their diversity and evolution. Here, we screened terrestrial arthropods for CRESS DNA viruses and report the identification of 44 viral genomes and replicons associated with specimens representing all three major terrestrial arthropod lineages, namely Euchelicerata (spiders), Hexapoda (insects), and Myriapoda (millipedes). We identified virus genomes belonging to three established CRESS DNA viral families (Circoviridae, Genomoviridae, and Smacoviridae); however, over half of the arthropod-associated viral genomes are only distantly related to currently classified CRESS DNA viral sequences. Although members of viral and satellite families known to infect plants (Geminiviridae, Nanoviridae, Alphasatellitidae) were not identified in this study, these plant-infecting CRESS DNA viruses and replicons are transmitted by hemipterans. Therefore, members from six out of the seven established CRESS DNA viral families circulate among arthropods. Furthermore, a phylogenetic analysis of Reps, including endogenous viral sequences, reported to date from a wide array of organisms revealed that most of the known CRESS DNA viral diversity circulates among invertebrates. Our results highlight the vast and unexplored diversity of CRESS DNA viruses among invertebrates and parallel findings from RNA viral discovery efforts in undersampled taxa.

Parasitic Nematodes Exert Antimicrobial Activity and Benefit From Microbiota-Driven Support for Host Immune Regulation

Intestinal parasitic nematodes live in intimate contact with the host microbiota. Changes in the microbiome composition during nematode infection affect immune control of the parasites and shifts in the abundance of bacterial groups have been linked to the immunoregulatory potential of nematodes. Here we asked if the small intestinal parasite Heligmosomoides polygyrus produces factors with antimicrobial activity, senses its microbial environment and if the anti-nematode immune and regulatory responses are altered in mice devoid of gut microbes. We found that H. polygyrus excretory/secretory products exhibited antimicrobial activity against gram^+/− bacteria. Parasites from germ-free mice displayed alterations in gene expression, comprising factors with putative antimicrobial functions such as chitinase and lysozyme. Infected germ-free mice developed increased small intestinal Th2 responses coinciding with a reduction in local Foxp3⁺RORγt⁺ regulatory T cells and decreased parasite fecundity. Our data suggest that nematodes sense their microbial surrounding and have evolved factors that limit the outgrowth of certain microbes. Moreover, the parasites benefit from microbiota-driven immune regulatory circuits, as an increased ratio of intestinal Th2 effector to regulatory T cells coincides with reduced parasite fitness in germ-free mice.

Population genetic analysis of Chadian Guinea worms reveals that human and non-human hosts share common parasite populations

Following almost 10 years of no reported cases, Guinea worm disease (GWD or dracunculiasis) reemerged in Chad in 2010 with peculiar epidemiological patterns and unprecedented prevalence of infection among non-human hosts, particularly domestic dogs. Since 2014, animal infections with Guinea worms have also been observed in the other three countries with endemic transmission (Ethiopia, Mali, and South Sudan), causing concern and generating interest in the parasites' true taxonomic identity and population genetics. We present the first extensive population genetic data for Guinea worm, investigating mitochondrial and microsatellite variation in adult female worms from both human and non-human hosts in the four endemic countries to elucidate the origins of Chad's current outbreak and possible host-specific differences between parasites. Genetic diversity of Chadian Guinea worms was considerably higher than that of the other three countries, even after controlling for sample size through rarefaction, and demographic analyses are consistent with a large, stable parasite population. Genealogical analyses eliminate the other three countries as possible sources of parasite reintroduction into Chad, and sequence divergence and distribution of genetic variation provide no evidence that parasites in human and non-human hosts are separate species or maintain isolated transmission cycles. Both among and within countries, geographic origin appears to have more influence on parasite population structure than host species. Guinea worm infection in non-human hosts has been occasionally reported throughout the history of the disease, particularly when elimination programs appear to be reaching their end goals. However, no previous reports have evaluated molecular support of the parasite species identity. Our data confirm that Guinea worms collected from non-human hosts in the remaining endemic countries of Africa are Dracunculus medinensis and that the same population of worms infects both humans and dogs in Chad. Our genetic data and the epidemiological evidence suggest that transmission in the Chadian context is currently being maintained by canine hosts.

Silent Witness: Dual-Species Transcriptomics Reveals Epithelial Immunological Quiescence to Helminth Larval Encounter and Fostered Larval Development

Gastrointestinal nematodes are among the most prevalent parasites infecting humans and livestock worldwide. Infective larvae of the soil-transmitted nematode Ascaris spp. enter the host and start tissue migration by crossing the intestinal epithelial barrier. The initial interaction of the intestinal epithelium with the parasite, however, has received little attention. In a time-resolved interaction model of porcine intestinal epithelial cells (IPEC-J2) and infective Ascaris suum larvae, we addressed the early transcriptional changes occurring simultaneously in both organisms using dual-species RNA-Seq. Functional analysis of the host response revealed an overall induction of metabolic activity, without induction of immune responsive genes or immune signaling pathways and showing suppression of chemotactic genes like CXCL8/IL-8 or CHI3L1. Ascaris larvae, when getting in contact with the epithelium, showed induction of genes that orchestrate motor activity and larval development, such as myosin, troponin, myoglobin, and protein disulfide isomerase 2 (PDI-2). In addition, excretory-secretory products that likely facilitate parasite invasion were increased, among them, aspartic protease 6 or hyaluronidase. Integration of host and pathogen data in an interspecies gene co-expression network indicated links between nematode fatty acid biosynthesis and host ribosome assembly/protein synthesis. In summary, our study provides new molecular insights into the early factors of parasite invasion, while at the same time revealing host immunological unresponsiveness. Reproducible software for dual RNA-Seq analysis of non-model organisms is available at https://gitlab.com/mkuhring/project_asuum and can be applied to similar studies.

In silico analyses of protein glycosylating genes in the helminth Fasciola hepatica (liver fluke) predict protein-linked glycan simplicity and reveal temporally-dynamic expression profiles

Glycoproteins secreted by helminth parasites are immunogenic and represent appealing components of vaccine preparations. Our poor knowledge of the pathways that mediate protein glycosylation in parasitic flatworms hinders our understanding of how proteins are synthesised and modified, and our ability to target these pathways for parasite control. Here we provide the first detailed description of genes associated with protein glycosylation in a parasitic flatworm, focusing on the genome of the liver fluke (Fasciola hepatica), which is a globally important trematode parasite of humans and their livestock. Using 190 human sequences as search queries against currently available F. hepatica genomes, we identified 149 orthologues with putative roles in sugar uptake or nucleotide sugar synthesis, and an array of glycosyltransferase and glycosidase activities required for protein N- and O-glycosylation. We found appreciable duplication within these orthologues, describing just 87 non-redundant genes when paralogues were excluded. F. hepatica lacks many of the enzymes required to produce complex N- and O-linked glycans, which explains the genomic basis for the structurally simple glycans described by F. hepatica glycomic datasets, and predicts pervasive structural simplicity in the wider glycome. These data provide a foundation for functional genomic interrogation of these pathways with the view towards novel parasite intervention strategies.

Dramatic evolution of body length due to postembryonic changes in cell size in a newly discovered close relative of Caenorhabditis elegans

Understanding morphological diversity-and morphological constraint-has been a central question in evolutionary biology since its inception. Nematodes of the genus Caenorhabditis, which contains the well-studied model organism C. elegans, display remarkable morphological consistency in the face of extensive genetic divergence. Here, we provide a description of the broad developmental patterns of a newly discovered species, C. sp. 34, which was isolated from fresh figs in Okinawa and which is among the closest known relatives of C. elegans. C. sp. 34 displays an extremely large body size; it can grow to be nearly twice as long as C. elegans and all other known members of the genus. Observations of the timing of developmental milestones reveal that C. sp. 34 develops about twice as slowly as C. elegans. Measurements of embryonic and larval size show that the size difference between C. sp. 34 and C. elegans is largely due to postembryonic events, particularly during the transition from larval to adult stages. This difference in size is not attributable to differences in germ line chromosome number or the number of somatic cells. The overall difference in body size is therefore largely attributable to changes in cell size via increased cytoplasmic volume. Because of its close relationship to C. elegans, the distinctness of C. sp. 34 provides an ideal system for the detailed analysis of evolutionary diversification. The context of over 40 years of C. elegans developmental genetics also reveals clues into how natural selection and developmental constraint act jointly to promote patterns of morphological stasis and divergence in this group.

Trichuris muris whey acidic protein induces type 2 protective immunity against whipworm

Human whipworm (Trichuris trichiura) infects approximately 1 in 15 people worldwide, representing the leading infectious cause of colitis and subsequent, inflammatory bowel disease (IBD). Current control measures focused on mass deworming have had limited success due to low drug efficacies. Vaccination would be an ideal, cost-effective strategy to induce protective immunity, leading to control of infection and transmission. Here we report the identification of whey acidic protein, a whipworm secretory protein, as a strong immunogen for inducing protective efficacy in a surrogate mouse T. muris infection model. The recombinant WAP protein (rTm-WAP49), as well as a single, highly conserved repeat within WAP (fragment 8) expressed as an Na-GST-1 fusion protein (rTm-WAP-F8+Na-GST-1), generate a strong T helper type 2 (Th2) immune response when delivered as subcutaneous vaccines formulated with Montanide ISA 720. Oral challenge with T. muris infective eggs following vaccination led to a significant reduction in worm burden of 48% by rTm-WAP49 and 33% by rTm-WAP-F8+Na-GST-1. The cellular immune correlates of protection included significant antigen-specific production of Th2 cytokines IL-4, IL-9, and IL-13 by cells isolated from the vaccine-draining inguinal lymph nodes, parasite-draining mesenteric lymph nodes, and spleen in mice vaccinated with either rTm-WAP49 or rTm-WAP-F8+Na-GST-1. The humoral immune correlates included a high antigen-specific ratio of IgG1 to IgG2a, without eliciting an IgE-mediated allergic response. Immunofluorescent staining of adult T. muris with WAP antisera identified the worm's pathogenic stichosome organ as the site of secretion of native Tm-WAP protein into the colonic mucosa. Given the high sequence conservation for the WAP proteins from T. muris and T. trichiura, the results presented here support the WAP protein to be further evaluated as a potential human whipworm vaccine candidate.

Evolution of gene dosage on the Z-chromosome of schistosome parasites

XY systems usually show chromosome-wide compensation of X-linked genes, while in many ZW systems, compensation is restricted to a minority of dosage-sensitive genes. Why such differences arose is still unclear. Here, we combine comparative genomics, transcriptomics and proteomics to obtain a complete overview of the evolution of gene dosage on the Z-chromosome of Schistosoma parasites. We compare the Z-chromosome gene content of African (Schistosoma mansoni and S. haematobium) and Asian (S. japonicum) schistosomes and describe lineage-specific evolutionary strata. We use these to assess gene expression evolution following sex-linkage. The resulting patterns suggest a reduction in expression of Z-linked genes in females, combined with upregulation of the Z in both sexes, in line with the first step of Ohno's classic model of dosage compensation evolution. Quantitative proteomics suggest that post-transcriptional mechanisms do not play a major role in balancing the expression of Z-linked genes.

A galactokinase-like protein from the liver fluke Fasciola hepatica

Galactokinase catalyses the ATP-dependent phosphorylation of galactose. A galactokinase-like sequence was identified in a Fasciola hepatica EST library. Recombinant expression of the corresponding protein in Escherichia coli resulted in a protein of approximately 50 kDa. The protein is monomeric, like galactokinases from higher animals, yeasts and some bacteria. The protein has no detectable enzymatic activity with galactose or N-acetylgalactosamine as a substrate. However, it does bind to ATP. Molecular modelling predicted that the protein adopts a similar fold to galactokinase and other GHMP kinases. However, a key loop in the active site was identified which may influence the lack of activity. Sequence analysis strongly suggested that this protein (and other proteins annotated as "galactokinase" in the trematodes Schistosoma mansoni and Clonorchis sinensis) are closer to N-acetylgalactosamine kinases. No other galactokinase-like sequences appear to be present in the genomes of these three species. This raises the intriguing possibility that these (and possibly other) trematodes are unable to catabolise galactose through the Leloir pathway due to the lack of a functional galactokinase.

Hypothesis: Multiple sclerosis is caused by three-hits, strictly in order, in genetically susceptible persons

Multiple Sclerosis is a chronic, progressive and debilitating neurological disease which, despite extensive study for over 100 years, remains of enigmatic aetiology. Drawn from the epidemiological evidence, there exists a consensus that there are environmental (possibly infectious) factors that contribute to disease pathogenesis that have not yet been fully elucidated. Here we propose a three-tiered hypothesis: 1) a clinic-epidemiological model of multiple sclerosis as a rare late complication of two sequential infections (with the temporal sequence of infections being important); 2) a proposal that the first event is helminthic infection with Enterobius Vermicularis, and the second is Epstein Barr Virus infection; and 3) a proposal for a testable biological mechanism, involving T-Cell exhaustion for Epstein-Barr Virus protein LMP2A. We believe that this model satisfies some of the as-yet unexplained features of multiple sclerosis epidemiology, is consistent with the clinical and neuropathological features of the disease and is potentially testable by experiment. This model may be generalizable to other autoimmune diseases.

Refined ab initio gene predictions of Heterorhabditis bacteriophora using RNA-seq

Interest has recently grown in developing the entomopathogenic nematode Heterorhabditis bacteriophora as a model to genetically dissect the process of parasitic infection. Despite the availability of a full genome assembly, there is substantial variation in gene model accuracy. Here, a methodology is presented for leveraging RNA-seq evidence to generate improved annotations using ab initio gene prediction software. After alignment of reads and subsequent generation of a RNA-seq supported annotation, the new gene prediction models were verified on a selection of genes by comparison with sequenced 5' and 3' rapid amplification of cDNA ends products. By utilising a whole transcriptome for genome annotation, the current reference annotation was enriched, demonstrating the importance of coupling transcriptional data with genome assemblies.

Nanomedicine Approaches Against Parasitic Worm Infections

Nanomedicine approaches have the potential to transform the battle against parasitic worm (helminth) infections, a major global health scourge from which billions are currently suffering. It is anticipated that the intersection of two currently disparate fields, nanomedicine and helminth biology, will constitute a new frontier in science and technology. This progress report surveys current innovations in these research fields and discusses research opportunities. In particular, the focus is on: (1) major challenges that helminth infections impose on mankind; (2) key aspects of helminth biology that inform future research directions; (3) efforts to construct nanodelivery platforms to target drugs and genes to helminths hidden in their hosts; (4) attempts in applying nanotechnology to enable vaccination against helminth infections; (5) outlooks in utilizing nanoparticles to enhance immunomodulatory activities of worm-derived factors to cure allergy and autoimmune diseases. In each section, achievements are summarized, limitations are explored, and future directions are assessed.

Introduction: Overview of Fungal Genomics

Genome sequences and postgenomic tools have had a major impact on fungal research. When the first fungal genome sequences became available it became clear how much more complex fungal biology was than had been previously assumed. Since then an increasing number of genomes have become available enabling detailed comparative studies, especially when combined with postgenomic tools such as transcriptomics and proteomics. This chapter provides an overview and current state of fungal genomics.

Methyl-CpG-binding (SmMBD2/3) and chromobox (SmCBX) proteins are required for neoblast proliferation and oviposition in the parasitic blood fluke Schistosoma mansoni

While schistosomiasis remains a significant health problem in low to middle income countries, it also represents a recently recognised threat to more economically-developed regions. Until a vaccine is developed, this neglected infectious disease is primarily controlled by praziquantel, a drug with a currently unknown mechanism of action. By further elucidating how Schistosoma molecular components cooperate to regulate parasite developmental processes, next generation targets will be identified. Here, we continue our studies on schistosome epigenetic participants and characterise the function of a DNA methylation reader, the Schistosoma mansoni methyl-CpG-binding domain protein (SmMBD2/3). Firstly, we demonstrate that SmMBD2/3 contains amino acid features essential for 5-methyl cytosine (5mC) binding and illustrate that adult schistosome nuclear extracts (females > males) contain this activity. We subsequently show that SmMBD2/3 translocates into nuclear compartments of transfected murine NIH-3T3 fibroblasts and recombinant SmMBD2/3 exhibits 5mC binding activity. Secondly, using a yeast-two hybrid (Y2H) screen, we show that SmMBD2/3 interacts with the chromo shadow domain (CSD) of an epigenetic adaptor, S. mansoni chromobox protein (SmCBX). Moreover, fluorescent in situ hybridisation (FISH) mediated co-localisation of Smmbd2/3 and Smcbx to mesenchymal cells as well as somatic- and reproductive- stem cells confirms the Y2H results and demonstrates that these interacting partners are ubiquitously expressed and found within both differentiated as well as proliferating cells. Finally, using RNA interference, we reveal that depletion of Smmbd2/3 or Smcbx in adult females leads to significant reductions (46-58%) in the number of proliferating somatic stem cells (PSCs or neoblasts) as well as in the quantity of in vitro laid eggs. Collectively, these results further expand upon the schistosome components involved in epigenetic processes and suggest that pharmacological inhibition of SmMBD2/3 and/or SmCBX biology could prove useful in the development of future schistosomiasis control strategies.

Molecular alterations during larval development of Haemonchus contortus in vitro are under tight post-transcriptional control

In this study, we explored the molecular alterations in the developmental switch from the L3 to the exsheathed L3 (xL3) and to the L4 stage of Haemonchus contortus in vitro using an integrated transcriptomic, proteomic and bioinformatic approach. Totals of 9,754 mRNAs, 88 microRNAs (miRNAs) and 1,591 proteins were identified, and 6,686 miRNA-mRNA pairs inferred in all larval stages studied. Approximately 16% of transcripts in the combined transcriptome (representing all three larval stages) were expressed as proteins, and there were positive correlations (r = 0.39-0.44) between mRNA transcription and protein expression in the three distinct developmental stages of the parasite. Of the predicted targets, 1,019 (27.0%) mRNA transcripts were expressed as proteins, and there was a negative correlation (r = -0.60 to -0.50) in the differential mRNA transcription and protein expression between developmental stages upon pairwise comparison. The changes in transcription (mRNA and miRNA) and protein expression from the free-living to the parasitic life cycle phase of H. contortus related to enrichments in biological pathways associated with metabolism (e.g., carbohydrate and lipid degradation, and amino acid metabolism), environmental information processing (e.g., signal transduction, signalling molecules and interactions) and/or genetic information processing (e.g., transcription and translation). Specifically, fatty acid degradation, steroid hormone biosynthesis and the Rap1 signalling pathway were suppressed, whereas transcription, translation and protein processing in the endoplasmic reticulum were upregulated during the transition from the free-living L3 to the parasitic xL3 and L4 stages of the nematode in vitro. Dominant post-transcriptional regulation was inferred to elicit these changes, and particular miRNAs (e.g., hco-miR-34 and hco-miR-252) appear to play roles in stress responses and/or environmental adaptations during developmental transitions of H. contortus. Taken together, these integrated results provide a comprehensive insight into the developmental biology of this important parasite at the molecular level in vitro. The approach applied here to H. contortus can be readily applied to other parasitic nematodes.

Assessment of the behaviour and survival of nematodes under low oxygen concentrations

Oxygen is required for the completion of almost all known metazoan lifecycles, but many metazoans harbour abilities to withstand varying degrees and periods of hypoxia. Caenorhabditis elegans, one of the most popular model organism is extensively used as a model for the study of hypoxia and anoxia biology and it has been found that this nematode is capable of tolerance to varying degrees of hypoxia. Considering the extremely high diversity of nematodes, the effects of low oxygen concentration and mechanisms of adaptation to oxygen depletion differ among species. In this study, we used a simple assay to examine anoxia tolerance in four nematode species, including three free-living and one plant parasitic nematode. We found that the plant parasitic nematode Bursaphelenchus xylophilus can survive more than 14 days under anoxic conditions. Comparisons of behaviour during anoxia induction and the repertoire of oxygen sensation genes among the tested species suggested the existence of different oxygen sensation systems between B. xylophilus and C. elegans, which quickly introduce suspended animation in response to oxygen depletion to survive long-term anoxia.

Genome-wide expert annotation of the epigenetic machinery of the plant-parasitic nematodes Meloidogyne spp., with a focus on the asexually reproducing species

Background

The renewed interest in epigenetics has led to the understanding that both the environment and individual lifestyle can directly interact with the epigenome to influence its dynamics. Epigenetic phenomena are mediated by DNA methylation, stable chromatin modifications and non-coding RNA-associated gene silencing involving specific proteins called epigenetic factors.

Multiple organisms, ranging from plants to yeast and mammals, have been used as model systems to study epigenetics. The interactions between parasites and their hosts are models of choice to study these mechanisms because the selective pressures are strong and the evolution is fast.

The asexually reproducing root-knot nematodes (RKN) offer different advantages to study the processes and mechanisms involved in epigenetic regulation. RKN genomes sequencing and annotation have identified numerous genes, however, which of those are involved in the adaption to an environment and potentially relevant to the evolution of plant-parasitism is yet to be discovered.

Results

Here, we used a functional comparative annotation strategy combining orthology data, mining of curated genomics as well as protein domain databases and phylogenetic reconstructions.

Overall, we show that (i) neither RKN, nor the model nematode Caenorhabditis elegans possess any DNA methyltransferases (DNMT) (ii) RKN do not possess the complete machinery for DNA methylation on the 6th position of adenine (6mA) (iii) histone (de)acetylation and (de)methylation pathways are conserved between C. elegans and RKN, and the corresponding genes are amplified in asexually reproducing RKN (iv) some specific non-coding RNA families found in plant-parasitic nematodes are dissimilar from those in C. elegans. In the asexually reproducing RKN Meloidogyne incognita, expression data from various developmental stages supported the putative role of these proteins in epigenetic regulations.

Conclusions

Our results refine previous predictions on the epigenetic machinery of model species and constitute the most comprehensive description of epigenetic factors relevant to the plant-parasitic lifestyle and/or asexual mode of reproduction of RKN. Providing an atlas of epigenetic factors in RKN is an informative resource that will enable researchers to explore their potential role in adaptation of these parasites to their environment.

Electronic supplementary material

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

Improved strategy for the curation and classification of kinases, with broad applicability to other eukaryotic protein groups

Despite the substantial amount of genomic and transcriptomic data available for a wide range of eukaryotic organisms, most genomes are still in a draft state and can have inaccurate gene predictions. To gain a sound understanding of the biology of an organism, it is crucial that inferred protein sequences are accurately identified and annotated. However, this can be challenging to achieve, particularly for organisms such as parasitic worms (helminths), as most gene prediction approaches do not account for substantial phylogenetic divergence from model organisms, such as Caenorhabditis elegans and Drosophila melanogaster, whose genomes are well-curated. In this paper, we describe a bioinformatic strategy for the curation of gene families and subsequent annotation of encoded proteins. This strategy relies on pairwise gene curation between at least two closely related species using genomic and transcriptomic data sets, and is built on recent work on kinase complements of parasitic worms. Here, we discuss salient technical aspects of this strategy and its implications for the curation of protein families more generally.

Identification and characterization of the Fasciola hepatica sodium- and chloride-dependent taurine transporter

The parasitic liver fluke Fasciola hepatica infests mainly ruminants, but it can also cause fasciolosis in people, who ingest the metacercariae encysted on plants. The drug of choice to treat fasciolosis is triclabendazole (TBZ), which has been on the market for several decades. This is also true for the other available drugs. Accordingly, drug-resistant flukes have been emerging at an increasing rate making it desirable to identify alternative drug targets. Here, we focused on the fact that adult F. hepatica persists in the hostile environment of the bile ducts of infected organisms. A common way to render bile acids less toxic is to conjugate them to taurine (2-aminoethanesulfonic acid). We cloned a transporter from the solute carrier-6 (SLC6) family, which was most closely related to the GABA-transporter-2 of other organisms. When heterologously expressed, this F. hepatica transporter supported the high-affinity cellular uptake of taurine (KM = 12.0 ± 0.5 μM) but not of GABA. Substrate uptake was dependent on Na+- and Cl- (calculated stoichiometry 2:1). Consistent with the low chloride concentration in mammalian bile, the F. hepatica transporter had a higher apparent affinity for Cl- (EC50 = 14±3 mM) than the human taurine transporter (EC50 = 55±7 mM). We incubated flukes with unconjugated bile acids in the presence and absence of taurine: taurine promoted survival of flukes; the taurine transporter inhibitor guanidinoethansulfonic acid abolished this protective effect of taurine. Based on these observations, we conclude that the taurine transporter is critical for the survival of liver flukes in the bile. Thus, the taurine transporter represents a candidate drug target.

Profiling G protein-coupled receptors of Fasciola hepatica identifies orphan rhodopsins unique to phylum Platyhelminthes

G protein-coupled receptors (GPCRs) are established drug targets. Despite their considerable appeal as targets for next-generation anthelmintics, poor understanding of their diversity and function in parasitic helminths has thwarted progress towards GPCR-targeted anti-parasite drugs. This study facilitates GPCR research in the liver fluke, Fasciola hepatica, by generating the first profile of GPCRs from the F. hepatica genome. Our dataset describes 147 high confidence GPCRs, representing the largest cohort of GPCRs, and the largest set of in silico ligand-receptor predictions, yet reported in any parasitic helminth. All GPCRs fall within the established GRAFS nomenclature; comprising three glutamate, 135 rhodopsin, two adhesion, five frizzled, one smoothened, and one secretin GPCR. Stringent annotation pipelines identified 18 highly diverged rhodopsins in F. hepatica that maintained core rhodopsin signatures, but lacked significant similarity with non-flatworm sequences, providing a new sub-group of potential flukicide targets. These facilitated identification of a larger cohort of 76 related sequences from available flatworm genomes, representing new members of existing groups (PROF1/Srfb, Rho-L, Rho-R, Srfa, Srfc) of flatworm-specific rhodopsins. These receptors imply flatworm specific GPCR functions, and/or co-evolution with unique flatworm ligands, and could facilitate the development of exquisitely selective anthelmintics. Ligand binding domain sequence conservation relative to deorphanised rhodopsins enabled high confidence ligand-receptor matching of seventeen receptors activated by acetylcholine, neuropeptide F/Y, octopamine or serotonin. RNA-Seq analyses showed expression of 101 GPCRs across various developmental stages, with the majority expressed most highly in the pathogenic intra-mammalian juvenile parasites. These data identify a broad complement of GPCRs in F. hepatica, including rhodopsins likely to have key functions in neuromuscular control and sensory perception, as well as frizzled and adhesion/secretin families implicated, in other species, in growth, development and reproduction. This catalogue of liver fluke GPCRs provides a platform for new avenues into our understanding of flatworm biology and anthelmintic discovery.

The anti-fecundity effect of 5-azacytidine (5-AzaC) on Schistosoma mansoni is linked to dis-regulated transcription, translation and stem cell activities

Uncontrolled host immunological reactions directed against tissue-trapped eggs precipitate a potentially lethal, pathological cascade responsible for schistosomiasis. Blocking schistosome egg production, therefore, presents a strategy for simultaneously reducing immunopathology as well as limiting disease transmission in endemic or emerging areas. We recently demonstrated that the ribonucleoside analogue 5-azacytidine (5-AzaC) inhibited Schistosoma mansoni oviposition, egg maturation and ovarian development. While these anti-fecundity effects were associated with a loss of DNA methylation, other molecular processes affected by 5-AzaC were not examined at the time. By comparing the transcriptomes of 5-AzaC-treated females to controls, we provide evidence that this ribonucleoside analogue also modulates other crucial aspects of schistosome egg-laying biology. For example, S. mansoni gene products associated with amino acid-, carbohydrate-, fatty acid-, nucleotide- and tricarboxylic acid (TCA)- homeostasis are all dysregulated in 5-AzaC treated females. To validate the metabolic pathway most significantly affected by 5-AzaC, amino acid metabolism, nascent protein synthesis was subsequently quantified in adult schistosomes. Here, 5-AzaC inhibited this process by 68% ±16.7% (SEM) in male- and 81% ±4.8% (SEM) in female-schistosomes. Furthermore, the transcriptome data indicated that adult female stem cells were also affected by 5-AzaC. For instance, 40% of transcripts associated with proliferating schistosome cells were significantly down-regulated by 5-AzaC. This finding correlated with a considerable reduction (95%) in the number of 5-ethynyl-2'-deoxyuridine (EdU) positive cells found in 5-AzaC-treated females. In addition to protein coding genes, the effect that 5-AzaC had on repetitive element expression was also assessed. Here, 46 repeats were found differentially transcribed between 5-AzaC-treated and control females with long terminal repeat (LTR) and DNA transposon classes being amongst the most significant. This study demonstrates that the anti-fecundity activity of 5-AzaC affects more than just DNA methylation in schistosome parasites. Further characterisation of these processes may reveal novel targets for schistosomiasis control.

Comparison of PCR Methods for Onchocerca volvulus Detection in Skin Snip Biopsies from the Tshopo Province, Democratic Republic of the Congo

Defining the optimal diagnostic tools for evaluating onchocerciasis elimination efforts in areas co-endemic for other filarial nematodes is imperative. This study compared three published polymerase chain reaction (PCR) methods: the Onchocerca volvulus-specific qPCR-O150, the pan-filarial qPCR melt curve analysis (MCA), and the O150-PCR enzyme-linked immunosorbent assay (ELISA) currently used for vector surveillance in skin snip biopsies (skin snips) collected from the Democratic Republic of the Congo. The pan-filarial qPCR-MCA was compared with species-specific qPCRs for Loa loa and Mansonella perstans. Among the 471 skin snips, 47.5%, 43.5%, and 27.0% were O. volvulus positive by qPCR-O150, qPCR-MCA, and O150-PCR ELISA, respectively. Using qPCR-O150 as the comparator, the sensitivity and specificity of qPCR-MCA were 89.3% and 98.0%, respectively, whereas for O150-PCR ELISA, they were 56.7% and 100%, respectively. Although qPCR-MCA identified the presence of L. loa and Mansonella spp. in skin snips, species-specific qPCRs had greater sensitivity and were needed to identify M. perstans. Most of the qPCR-MCA misclassifications occurred in mixed infections. The reduced sensitivity of O150-PCR ELISA was associated with lower microfilaria burden and with lower amounts of O. volvulus DNA. Although qPCR-MCA identified most of the O. volvulus-positive skin snips, it is not sufficiently robust to be used for stop-mass drug administration (MDA) evaluations in areas co-endemic for other filariae. Because O150-PCR ELISA missed 43.3% of qPCR-O150-positive skin snips, the qPCR-O150 assay is more appropriate for evaluating skin snips of OV-16 + children in stop-MDA assessments. Although improving the sensitivity of the O150-PCR ELISA as an alternative to qPCR might be possible, qPCR-O150 offers distinct advantages aside from increased sensitivity.

Flatworm-specific transcriptional regulators promote the specification of tegumental progenitors in Schistosoma mansoni

Schistosomes infect more than 200 million people. These parasitic flatworms rely on a syncytial outer coat called the tegument to survive within the vasculature of their host. Although the tegument is pivotal for their survival, little is known about maintenance of this tissue during the decades schistosomes survive in the bloodstream. Here, we demonstrate that the tegument relies on stem cells (neoblasts) to specify fusogenic progenitors that replace tegumental cells lost to turnover. Molecular characterization of neoblasts and tegumental progenitors led to the discovery of two flatworm-specific zinc finger proteins that are essential for tegumental cell specification. These proteins are homologous to a protein essential for neoblast-driven epidermal maintenance in free-living flatworms. Therefore, we speculate that related parasites (i.e., tapeworms and flukes) employ similar strategies to control tegumental maintenance. Since parasitic flatworms infect every vertebrate species, understanding neoblast-driven tegumental maintenance could identify broad-spectrum therapeutics to fight diseases caused by these parasites.

Comparative proteomics of the larval and adult stages of the model cestode parasite Mesocestoides corti

Mesocestoides corti is a widely used model for the study of cestode biology, and its transition from the larval tetrathyridium (TT) stage to the strobilated, adult worm (ST) stage can be induced and followed in vitro. Here, a proteomic approach was used to describe and compare M. corti TT and ST protein repertories. Overall, 571 proteins were identified, 238 proteins in TT samples and 333 proteins in ST samples. Among the identified proteins, 207 proteins were shared by TTs and STs, while 157 were stage-specific, being 31 exclusive from TTs, and 126 from STs. Functional annotation revealed fundamental metabolic differences between the TT and the ST stages. TTs perform functions related mainly to basic metabolism, responsible for growth and vegetative development by asexual reproduction. STs, in contrast, perform a wider range of functions, including macromolecule biosynthetic processes, gene expression and control pathways, which may be associated to its proglottization/segmentation, sexual differentiation and more complex physiology. Furthermore, the generated results provided an extensive list of cestode proteins of interest for functional studies in M. corti. Many of these proteins are novel candidate diagnostic antigens, and/or potential targets for the development of new and more effective antihelminthic drugs.

BIOLOGICAL SIGNIFICANCE

Cestodiases are parasitic diseases with serious impact on human and animal health. Efforts to develop more effective strategies for diagnosis, treatment or control of cestodiases are impaired by the still limited knowledge on many aspects of cestode biology, including the complex developmental processes that occur in the life cycles of these parasites. Mesocestoides corti is a good experimental model to study the transition from the larval to the adult stage, called strobilation, which occur in typical cestode life-cycles. The performed proteomics approach provided large-scale identification and quantification of M. corti proteins. Many stage-specific or differentially expressed proteins were detected in the larval tetrathyridium (TT) stage and in the strobilated, adult worm (ST) stage. Functional comparative analyses of the described protein repertoires shed light on function and processes associated to specific features of both stages, such as less differentiation and asexual reproduction in TTs, and proglottization/segmentation and sexual differentiation in ST. Moreover, many of the identified stage-specific proteins are useful as cestode developmental markers, and are potential targets for development of novel diagnostic methods and therapeutic drugs for cestodiases.

Advances in kinome research of parasitic worms - implications for fundamental research and applied biotechnological outcomes

Protein kinases are enzymes that play essential roles in the regulation of many cellular processes. Despite expansions in the fields of genomics, transcriptomics and bioinformatics, there is limited information on the kinase complements (kinomes) of most eukaryotic organisms, including parasitic worms that cause serious diseases of humans and animals. The biological uniqueness of these worms and the draft status of their genomes pose challenges for the identification and classification of protein kinases using established tools. In this article, we provide an account of kinase biology, the roles of kinases in diseases and their importance as drug targets, and drug discovery efforts in key socioeconomically important parasitic worms. In this context, we summarise methods and resources commonly used for the curation, identification, classification and functional annotation of protein kinase sequences from draft genomes; review recent advances made in the characterisation of the worm kinomes; and discuss the implications of these advances for investigating kinase signalling and developing small-molecule inhibitors as new anti-parasitic drugs.

Transcriptomic profiling of nematode parasites surviving vaccine exposure

•

Surviving Haemonchus contortus from vaccinated sheep were compared with control worms.

•

There is no evidence for changes in expression of genes encoding Barbervax® antigens.

•

There was increased expression of other proteases and regulators of lysosome trafficking.

•

Surviving worms displayed up-regulated lipid storage and defecation abilities.

Some nematode species are economically important parasites of livestock, while others are important human pathogens causing some of the most important neglected tropical diseases. In both humans and animals, anthelmintic drug administration is the main control strategy, but the emergence of drug-resistant worms has stimulated the development of alternative control approaches. Among these, vaccination is considered to be a sustainable and cost effective strategy. Currently, Barbervax® for the ruminant strongylid Haemonchus contortus is the only registered subunit vaccine for a nematode parasite, although a vaccine for the human hookworm Necator americanus is undergoing clinical trials (HOOKVAC consortium). As both these vaccines comprise a limited number of proteins, there is potential for selection of nematodes with altered sequences or expression of the vaccine antigens. Here we compared the transcriptome of H. contortus populations from sheep vaccinated with Barbervax® with worms from control animals. Barbervax® antigens are native integral membrane proteins isolated from the brush border of the intestinal cells of the adult parasite and many of those are proteases. Our findings provide no evidence for changes in expression of genes encoding Barbervax® antigens in the surviving parasite populations. However, surviving parasites from vaccinated animals showed increased expression of other proteases and regulators of lysosome trafficking, and displayed up-regulated lipid storage and defecation abilities that may have circumvented the effect of the vaccine. Implications for other potential vaccines for human and veterinary nematodes are discussed.

Cysteine peptidases of Eudiplozoon nipponicum: a broad repertoire of structurally assorted cathepsins L in contrast to the scarcity of cathepsins B in an invasive species of haematophagous monogenean of common carp

Background

Cysteine peptidases of clan CA, family C1 account for a major part of proteolytic activity in the haematophagous monogenean Eudiplozoon nipponicum. The full spectrum of cysteine cathepsins is, however, unknown and their particular biochemical properties, tissue localisation, and involvement in parasite-host relationships are yet to be explored.

Methods

Sequences of cathepsins L and B (EnCL and EnCB) were mined from E. nipponicum transcriptome and analysed bioinformatically. Genes encoding two EnCLs and one EnCB were cloned and recombinant proteins produced in vitro. The enzymes were purified by chromatography and their activity towards selected substrates was characterised. Antibodies and specific RNA probes were employed for localisation of the enzymes/transcripts in tissues of E. nipponicum adults.

Results

Transcriptomic analysis revealed a set of ten distinct transcripts that encode EnCLs. The enzymes are significantly variable in their active sites, specifically the S2 subsites responsible for interaction with substrates. Some of them display unusual structural features that resemble cathepsins B and S. Two recombinant EnCLs had different pH activity profiles against both synthetic and macromolecular substrates, and were able to hydrolyse blood proteins and collagen I. They were localised in the haematin cells of the worm’s digestive tract and in gut lumen. The EnCB showed similarity with cathepsin B2 of Schistosoma mansoni. It displays molecular features typical of cathepsins B, including an occluding loop responsible for its exopeptidase activity. Although the EnCB hydrolysed haemoglobin in vitro, it was localised in the vitelline cells of the parasite and not the digestive tract.

Conclusions

To our knowledge, this study represents the first complex bioinformatic and biochemical characterisation of cysteine peptidases in a monogenean. Eudiplozoon nipponicum adults express a variety of CLs, which are the most abundant peptidases in the worms. The properties and localisation of the two heterologously expressed EnCLs indicate a central role in the (partially extracellular?) digestion of host blood proteins. High variability of substrate-binding sites in the set of EnCLs suggests specific adaptation to a range of biological processes that require proteolysis. Surprisingly, a single cathepsin B is expressed by the parasite and it is not involved in digestion, but probably in vitellogenesis.

Electronic supplementary material

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

TGF-β mimic proteins form an extended gene family in the murine parasite Heligmosomoides polygyrus

We recently reported the discovery of a new parasite-derived protein that functionally mimics the immunosuppressive cytokine transforming growth factor (TGF)-β. The Heligmosomoides polygyrus TGF-β Mimic (Hp-TGM) shares no homology to any TGF-β family member, however it binds the mammalian TGF-β receptor and induces expression of Foxp3, the canonical transcription factor of both mouse and human regulatory T cells. Hp-TGM consists of five atypical Complement Control Protein (CCP, Pfam 00084) domains, each lacking certain conserved residues and 12-15 amino acids longer than the 60-70 amino acids consensus domain, but with a recognizable 3-cysteine, tryptophan, cysteine motif. We now report on the identification of a family of nine related Hp-TGM homologues represented in the secreted proteome and transcriptome of H. polygyrus. Recombinant proteins from five of the nine new TGM members were tested for TGF-β activity, but only two were functionally active in an MFB-F11 reporter assay, and by the induction of T cell Foxp3 expression. Sequence comparisons reveal that proteins with functional activity are similar or identical to Hp-TGM across the first three CCP domains, but more variable in domains 4 and 5. Inactive proteins diverged in all domains, or lacked some domains entirely. Testing truncated versions of Hp-TGM confirmed that domains 1-3 are essential for full activity in vitro, while domains 4 and 5 are not required. Further studies will elucidate whether these latter domains fulfill other functions in promoting host immune regulation during infection and if the more divergent family members play other roles in immunomodulation.

Serine proteases in schistosomes and other trematodes

Trematodes, also known as flukes, are phylogenetically ancient parasitic organisms. Due to their importance as human and veterinary parasites, their proteins have been investigated extensively as drug and vaccine targets. Among those, proteases, as crucial enzymes for parasite survival, are considered candidate molecules for anti-parasitic interventions. Surprisingly however, trematode serine proteases, in comparison with other groups of proteases, are largely neglected. Genes encoding serine proteases have been identified in trematode genomes in significant abundance, but the biological roles and biochemical functions of these proteases are poorly understood. However, increasing volumes of genomic and proteomic studies, and accumulated experimental evidence, indicate that this class of proteases plays a substantial role in host-parasite interactions and parasite survival. Here, we discuss in detail serine proteases at genomic and protein levels, and their known or hypothetical functions.

Advances in Fasciola hepatica research using 'omics' technologies

The liver fluke Fasciola hepatica is an economically important pathogen of livestock worldwide, as well as being an important neglected zoonosis. Parasite control is reliant on the use of drugs, particularly triclabendazole, which is effective against multiple parasite stages. However, the spread of parasites resistant to triclabendazole has intensified the pursuit for novel control strategies. Emerging 'omics' technologies are helping advance our understanding of liver fluke biology, specifically the molecules that act at the host-parasite interface and are central to infection, virulence and long-term survival within the definitive host. This review discusses the technological sequencing advances that have facilitated the unbiased analysis of liver fluke biology, resulting in an extensive range of 'omics' datasets. In addition, we highlight the 'omics' studies of host responses to F. hepatica infection that, when combined with the parasite datasets, provide the opportunity for integrated analyses of host-parasite interactions. These extensive datasets will form the foundation for future in-depth analysis of F. hepatica biology and development, and the search for new drug or vaccine interventions.

Pleiotropic alterations in gene expression in Latin American Fasciola hepatica isolates with different susceptibility to drugs

Background

Fasciola hepatica is the main agent of fasciolosis, a zoonotic disease affecting livestock worldwide, and an emerging food-borne disease in humans. Even when effective treatments are available, drugs are costly and can result in tolerance, liver damage and normally they do not prevent reinfection. Drug-resistant strains in livestock have been reported in various countries and, more worryingly, drug resistance in human cases has emerged in South America. The present study aims to characterize the transcriptome of two South American resistant isolates, the Cajamarca isolate from Peru, resistant to both triclabendazole and albendazole (TCBZR/ABZR) and the Rubino isolate from Uruguay, resistant to ABZ (TCBZS/ABZR), and compare them to a sensitive strain (Cenapa, Mexico, TCBZS/ABZS) to reveal putative molecular mechanisms leading to drug resistance.

Results

We observed a major reduction in transcription in the Cajamarca TCBZR/ABZR isolate in comparison to the other isolates. While most of the differentially expressed genes are still unannotated, several trends could be detected. Specific reduction in the expression levels of cytoskeleton proteins was consistent with a role of tubulins as putative targets of triclabendazole (TCBZ). A marked reduction of adenylate cyclase might be underlying pleiotropic effects on diverse metabolic pathways of the parasite. Upregulation of GST mu isoforms suggests this detoxifying mechanism as one of the strategies associated with resistance.

Conclusions

Our results stress the value of transcriptomic approaches as a means of providing novel insights to advance the understanding of drug mode of action and drug resistance. The results provide evidence for pleiotropic variations in drug-resistant isolates consistent with early observations of TCBZ and ABZ effects and recent proteomic findings.

Electronic supplementary material

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

Tissue-specific transcriptomes of Anisakis simplex (sensu stricto) and Anisakis pegreffii reveal potential molecular mechanisms involved in pathogenicity

Background

Larval stages of the sibling species of parasitic nematodes Anisakis simplex (sensu stricto) (s.s.) (AS) and Anisakis pegreffii (AP) are responsible for a fish-borne zoonosis, known as anisakiasis, that humans aquire via the ingestion of raw or undercooked infected fish or fish-based products. These two species differ in geographical distribution, genetic background and peculiar traits involved in pathogenicity. However, thus far little is known of key molecules potentially involved in host-parasite interactions. Here, high-throughput RNA-Seq and bioinformatics analyses of sequence data were applied to the characterization of the whole sets of transcripts expressed by infective larvae of AS and AP, as well as of their pharyngeal tissues, in a bid to identify transcripts potentially involved in tissue invasion and host-pathogen interplay.

Results

Approximately 34,000,000 single-end reads were generated from cDNA libraries for each species. Transcripts identified in AS and AP encoded 19,403 and 10,424 putative peptides, respectively, and were classified based on homology searches, protein motifs, gene ontology and biological pathway mapping. Differential gene expression analysis yielded 226 and 339 transcripts upregulated in the pharyngeal regions of AS and AP, respectively, compared with their corresponding whole-larvae datasets. These included proteolytic enzymes, molecules encoding anesthetics, inhibitors of primary hemostasis and virulence factors, anticoagulants and immunomodulatory peptides.

Conclusions

This work provides the scientific community with a list of key transcripts expressed by AS and AP pharyngeal tissues and corresponding annotation information which represents a ready-to-use resource for future functional studies of biological pathways specifically involved in host-parasite interplay.

Electronic supplementary material

The online version of this article (10.1186/s13071-017-2585-7) contains supplementary material, which is available to authorized users.

Rapid Gene Family Evolution of a Nematode Sperm Protein Despite Sequence Hyper-conservation

Reproductive proteins are often observed to be the most rapidly evolving elements within eukaryotic genomes. The major sperm protein (MSP) is unique to the phylum Nematoda and is required for proper sperm locomotion and fertilization. Here, we annotate the MSP gene family and analyze their molecular evolution in 10 representative species across Nematoda. We show that MSPs are hyper-conserved across the phylum, having maintained an amino acid sequence identity of 83.5–97.7% for over 500 million years. This extremely slow rate of evolution makes MSPs some of the most highly conserved genes yet identified. However, at the gene family level, we show hyper-variability in both gene copy number and genomic position within species, suggesting rapid, lineage-specific gene family evolution. Additionally, we find evidence that extensive gene conversion contributes to the maintenance of sequence identity within chromosome-level clusters of MSP genes. Thus, while not conforming to the standard expectation for the evolution of reproductive proteins, our analysis of the molecular evolution of the MSP gene family is nonetheless consistent with the widely repeatable observation that reproductive proteins evolve rapidly, in this case in terms of the genomic properties of gene structure, copy number, and genomic organization. This unusual evolutionary pattern is likely generated by strong pleiotropic constraints acting on these genes at the sequence level, balanced against expansion at the level of the whole gene family.

Using WormBase ParaSite: An Integrated Platform for Exploring Helminth Genomic Data

WormBase ParaSite ( parasite.wormbase.org ) is a comprehensive resource for the genomes of parasitic nematodes and flatworms (helminths). It currently includes genomic data for over 100 helminth species, adding value by way of consistent functional annotation, gene comparative analysis and gene expression analysis. We provide several ways of exploring the data including a choice of genome browsers, genome and gene summary pages, text and sequence searching, a query wizard, bulk downloads, and programmatic interfaces. WormBase ParaSite is released three to six times per year, and is developed in collaboration with WormBase ( www.wormbase.org ) and Ensembl Genomes ( www.ensemblgenomes.org ).