Search

Background The parasitic flatworm Clonorchis sinensis inhabits the biliary tree of humans and other piscivorous mammals. This parasite can survive and thrive in the bile duct, despite exposure to bile constituents and host immune attack. Although the precise biological mechanisms underlying this adaptation are unknown, previous work indicated that Niemann-pick type C2 (NPC2)-like sterol-binding proteins might be integral in the host-parasite interplay. Expansions of this family in some invertebrates, such as arthropods, have shown functional diversification, including novel forms of chemoreception. Thus, here we curated the NPC2-like protein gene complement in C. sinensis, and predicted their conserved and/or divergent functional roles. Methods We used an established comparative genomic-bioinformatic approach to curate NPC2-like proteins encoded in published genomes of Korean and Chinese isolates of C. sinensis. Protein sequence and structural homology, presence of conserved domains and phylogeny were used to group and functionally classify NPC2-like proteins. Furthermore, transcription levels of NPC2-like protein-encoding genes were explored in different developmental stages and tissues. Results Totals of 35 and 32 C. sinensis NPC2-like proteins were predicted to be encoded in the genomes of the Korean and Chinese isolates, respectively. Overall, these proteins had low sequence homology and high variability of sequence alignment coverage when compared with curated NPC2s. Most C. sinensis proteins were predicted to retain a conserved ML domain and a conserved fold conformation, with a large cavity within the protein. Only one protein sequence retained the conserved amino acid residues required in bovine NPC2 to bind cholesterol. Non-canonical C. sinensis NPC2-like protein-coding domains clustered into four distinct phylogenetic groups with members of a group frequently encoded on the same genome scaffolds. Interestingly, NPC2-like protein-encoding genes were predicted to be variably transcribed in different developmental stages and adult tissues, with most being transcribed in the metacercarial stage. Conclusions The results of the present investigation confirms an expansion of NPC2-like proteins in C. sinensis, suggesting a diverse array of functions beyond sterol binding and transport. Functional explorations of this protein family should elucidate the mechanisms enabling the establishment and survival of C. sinensis and related flukes in the biliary systems of mammalian hosts.

Background Echinococcus tapeworms cause a severe helminthic zoonosis called echinococcosis. The genus comprises various species and genotypes, of which E. granulosus (sensu stricto) represents a significant global public health and socioeconomic burden. Mitochondrial (mt) genomes have provided useful genetic markers to explore the nature and extent of genetic diversity within Echinococcus and have underpinned phylogenetic and population structure analyses of this genus. Our recent work indicated a sequence gap (> 1 kb) in the mt genomes of E. granulosus genotype G1, which could not be determined by PCR-based Sanger sequencing. The aim of the present study was to define the complete mt genome, irrespective of structural complexities, using a long-read sequencing method. Methods We extracted high molecular weight genomic DNA from protoscoleces from a single cyst of E. granulosus genotype G1 from a sheep from Australia using a conventional method and sequenced it using PacBio Sequel (long-read) technology, complemented by BGISEQ-500 short-read sequencing. Sequence data obtained were assembled using a recently-developed workflow. Results We assembled a complete mt genome sequence of 17,675 bp, which is > 4 kb larger than the complete mt genomes known for E. granulosus genotype G1. This assembly includes a previously-elusive tandem repeat region, which is 4417 bp long and consists of ten near-identical 441–445 bp repeat units, each harbouring a 184 bp non-coding region and adjacent regions. We also identified a short non-coding region of 183 bp, which includes an inverted repeat. Conclusions We report what we consider to be the first complete mt genome of E. granulosus genotype G1 and characterise all repeat regions in this genome. The numbers, sizes, sequences and functions of tandem repeat regions remain to be studied in different isolates of genotype G1 and in other genotypes and species. The discovery of such ‘new’ repeat elements in the mt genome of genotype G1 by PacBio sequencing raises a question about the completeness of some published genomes of taeniid cestodes assembled from conventional or short-read sequence datasets. This study shows that long-read sequencing readily overcomes the challenges of assembling repeat elements to achieve improved genomes.

Experimental studies of Caenorhabditis elegans and Drosophila melanogaster have contributed substantially to our understanding of molecular and cellular processes in metazoans at large. Since the publication of their genomes, functional genomic investigations have identified genes that are essential or non-essential for survival in each species. Recently, a range of features linked to gene essentiality have been inferred using a machine learning (ML)-based approach, allowing essentiality predictions within a species. Nevertheless, predictions between species are still elusive. Here, we undertake a comprehensive study using ML to discover and validate features of essential genes common to both C. elegans and D. melanogaster. We demonstrate that the cross-species prediction of gene essentiality is possible using a subset of features linked to nucleotide/protein sequences, protein orthology and subcellular localisation, single-cell RNA-seq, and histone methylation markers. Complementary analyses showed that essential genes are enriched for transcription and translation functions and are preferentially located away from heterochromatin regions of C. elegans and D. melanogaster chromosomes. The present work should enable the cross-prediction of essential genes between model and non-model metazoans.

Long non-coding, tandem-repetitive regions in mitochondrial (mt) genomes of many metazoans have been notoriously difficult to characterise accurately using conventional sequencing methods. Here, we show how the use of a third-generation (long-read) sequencing and informatic approach can overcome this problem. We employed Oxford Nanopore technology to sequence genomic DNAs from a pool of adult worms of the carcinogenic parasite, Schistosoma haematobium, and used an informatic workflow to define the complete mt non-coding region(s). Using long-read data of high coverage, we defined six dominant mt genomes of 33.4 kb to 22.6 kb. Although no variation was detected in the order or lengths of the protein-coding genes, there was marked length (18.5 kb to 7.6 kb) and structural variation in the non-coding region, raising questions about the evolution and function of what might be a control region that regulates mt transcription and/or replication. The discovery here of the largest tandem-repetitive, non-coding region (18.5 kb) in a metazoan organism also raises a question about the completeness of some of the mt genomes of animals reported to date, and stimulates further explorations using a Nanopore-informatic workflow.

Live attenuated vaccines are commonly used to control Mycoplasma gallisepticum infections in chickens. M. gallisepticum ts-304 is a novel live attenuated vaccine strain that has been shown to be safe and effective. In this study, the transcriptional profiles of genes in the tracheal mucosa in chickens challenged with the M. gallisepticum wild-type strain Ap3AS at 57 weeks after vaccination with ts-304 were explored and compared with the profiles of unvaccinated chickens that had been challenged with strain Ap3AS, unvaccinated and unchallenged chickens, and vaccinated but unchallenged chickens. At two weeks after challenge, pair-wise comparisons of transcription in vaccinated-only, vaccinated-and-challenged and unvaccinated and unchallenged birds detected no differences. However, the challenged-only birds had significant up-regulation in the transcription of genes and enrichment of gene ontologies, pathways and protein classes involved in infiltration and proliferation of inflammatory cells and immune responses mediated through enhanced cytokine and chemokine production and signaling, while those predicted to be involved in formation and motor movement of cilia and formation of the cellular cytoskeleton were significantly down-regulated. The transcriptional changes associated with the inflammatory response were less severe in these mature birds than in the relatively young birds examined in a previous study. The findings of this study demonstrated that vaccination with the attenuated M. gallisepticum strain ts-304 protects against the transcriptional changes associated with the inflammatory response and pathological changes in the tracheal mucosa caused by infection with M. gallisepticum in chickens for at least 57 weeks after vaccination.

Background Larvae of the Australian sheep blowfly, Lucilia cuprina, parasitise sheep by feeding on skin excretions, dermal tissue and blood, causing severe damage known as flystrike or myiasis. Recent advances in -omic technologies and bioinformatic data analyses have led to a greater understanding of blowfly biology and should allow the identification of protein families involved in host-parasite interactions and disease. Current literature suggests that proteins of the SCP (Sperm-Coating Protein)/TAPS (Tpx-1/Ag5/PR-1/Sc7) (SCP/TAPS) superfamily play key roles in immune modulation, cross-talk between parasite and host as well as developmental and reproductive processes in parasites. Methods Here, we employed a bioinformatics workflow to curate the SCP/TAPS protein gene family in L. cuprina. Protein sequence, the presence and number of conserved CAP-domains and phylogeny were used to group identified SCP/TAPS proteins; these were compared to those found in Drosophila melanogaster to make functional predictions. In addition, transcription levels of SCP/TAPS protein-encoding genes were explored in different developmental stages. Results A total of 27 genes were identified as belonging to the SCP/TAPS gene family: encoding 26 single-domain proteins each with a single CAP domain and a solitary double-domain protein containing two conserved cysteine-rich secretory protein/antigen 5/pathogenesis related-1 (CAP) domains. Surprisingly, 16 SCP/TAPS predicted proteins formed an extended tandem array spanning a 53 kb region of one genomic region, which was confirmed by MinION long-read sequencing. RNA-seq data indicated that these 16 genes are highly transcribed in all developmental stages (excluding the embryo). Conclusions Future work should assess the potential of selected SCP/TAPS proteins as novel targets for the control of L. cuprina and related parasitic flies of major socioeconomic importance. Graphical Abstract

Clonorchiasis is a neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, and is often associated with a malignant form of bile duct cancer (cholangiocarcinoma). Although some aspects of the epidemiology of clonorchiasis are understood, little is known about the genetics of C. sinensis populations. Here, we conducted a comprehensive genetic exploration of C. sinensis from endemic geographic regions using complete mitochondrial protein gene sets. Genomic DNA samples from C. sinensis individuals (n = 183) collected from cats and dogs in China (provinces of Guangdong, Guangxi, Hunan, Heilongjiang and Jilin) as well as from rats infected with metacercariae from cyprinid fish from the Russian Far East (Primorsky Krai region) were deep sequenced using the BGISEQ-500 platform. Informatic analyses of mitochondrial protein gene data sets revealed marked genetic variation within C. sinensis; significant variation was identified within and among individual worms from distinct geographical locations. No clear affiliation with a particular location or host species was evident, suggesting a high rate of dispersal of the parasite across endemic regions. The present work provides a foundation for future biological, epidemiological and ecological studies using mitochondrial protein gene data sets, which could aid in elucidating associations between particular C. sinensis genotypes/haplotypes and the pathogenesis or severity of clonorchiasis and its complications (including cholangiocarcinoma) in humans., Author summary Clonorchiasis is an important neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, which can induce malignant bile duct cancer (cholangiocarcinoma). Little precise information is available on the biology, epidemiology and population genetics of C. sinensis. For this reason, we explored here the genetic composition of C. sinensis populations in distinct endemic areas in China and Russia. Using a deep sequencing-informatic approach, we revealed marked mitochondrial genetic variation within and between individuals and populations of C. sinensis, with no particular affiliation with geographic or host origin. These molecular findings and the methodology established should underpin future genetic studies of C. sinensis causing human clonorchiasis and associated complications (cancer) as well as transmission patterns in endemic regions.

The parasitic mite Sarcoptes scabiei is an economically highly significant parasite of the skin of humans and animals worldwide. In humans, this mite causes a neglected tropical disease (NTD), called scabies. This disease results in major morbidity, disability, stigma and poverty globally and is often associated with secondary bacterial infections. Currently, anti-scabies treatments are not sufficiently effective, resistance to them is emerging and no vaccine is available. Here, we report the first high-quality genome and transcriptomic data for S. scabiei. The genome is 56.6 Mb in size, has a a repeat content of 10.6% and codes for 9,174 proteins. We explored key molecules involved in development, reproduction, host-parasite interactions, immunity and disease. The enhanced ‘omic data sets for S. scabiei represent comprehensive and critical resources for genetic, functional genomic, metabolomic, phylogenetic, ecological and/or epidemiological investigations, and will underpin the design and development of new treatments, vaccines and/or diagnostic tests., Author summary Scabies is a highly significant parasitic disease caused by the mite S. scabiei. This NTD has a major adverse impact in disadvantaged communities around the world, particularly when associated with secondary bacterial infections and clinical complications. Here we report the first high-quality genome and transcriptomic data for S. scabiei and explore molecular aspects of S. scabiei/scabies. This genome (56.6 Mb, encoding ~ 9,200 proteins) provides a solid foundation for fundamental investigations of the molecular biology of the mite, host-parasite interactions and disease processes as well as for translational research to develop new treatments, vaccines and diagnostic tests.

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.

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.

Neglected tropical diseases cause substantial morbidity and mortality in animals and people globally. Opisthorchiasis is one such disease, caused by the carcinogenic, Asian liver fluke, Opisthorchis viverrini. This hepatobiliary disease is known to be associated with malignant cancer (cholangiocarcinoma, CCA) and affects millions of people in Asia, including Thailand, Lao People's Democratic Republic (PDR) and Cambodia. No vaccine is available, and only one drug (praziquantel) is routinely employed against the parasite. Relatively little is known about the molecular biology of the fluke itself and the disease complex that it causes in humans. With the advent of high-throughput nucleic acid sequencing and bioinformatic technologies, it has now become possible to gain global insights into the molecular biology of parasites. The purpose of this minireview is (i) to discuss recent progress on the genomics of parasitic worms, with an emphasis on the draft genome and transcriptome of O. viverrini; (ii) to use results from an integrated, global analysis of the genomic and transcriptomic data, to explain how we believe that this carcinogenic fluke establishes in the biliary system, how it feeds, survives and protects itself in such a hostile, microaerobic environment within the liver, and to propose how this parasite evades or modulates host attack; and (iii) to indicate some of the challenges, and, more importantly, the exciting opportunities that the 'omic resources for O. viverrini now provide for a plethora of fundamental and applied research areas. Looking ahead, we hope that this genomic resource stimulates vibrant and productive collaborations within a consortium context, focused on the effective control of opisthorchiasis.

As a result of limited classes of anthelmintics and an over-reliance on chemical control, there is a great need to discover new compounds to combat drug resistance in parasitic nematodes. Here, we show that deguelin, a plant-derived rotenoid, selectively and potently inhibits the motility and development of nematodes, which supports its potential as a lead candidate for drug development. Furthermore, we demonstrate that deguelin treatment significantly increases gene transcription that is associated with energy metabolism, particularly oxidative phosphorylation and mitoribosomal protein production before inhibiting motility. Mitochondrial tracking confirmed enhanced oxidative phosphorylation. In accordance, real-time measurements of oxidative phosphorylation in response to deguelin treatment demonstrated an immediate decrease in oxygen consumption in both parasitic (Haemonchus contortus) and free-living (Caenorhabditis elegans) nematodes. Consequently, we hypothesize that deguelin is exerting its toxic effect on nematodes as a modulator of oxidative phosphorylation. This study highlights the dynamic biologic response of multicellular organisms to deguelin perturbation.-Preston, S., Korhonen, P. K., Mouchiroud, L., Cornaglia, M., McGee, S. L., Young, N. D., Davis, R. A., Crawford, S., Nowell, C., Ansell, B. R. E., Fisher, G. M., Andrews, K. T., Chang, B. C. H., Gijs, M. A. M., Sternberg, P. W., Auwerx, J., Baell, J., Hofmann, A., Jabbar, A., Gasser, R. B. Deguelin exerts potent nematocidal activity via the mitochondrial respiratory chain.

Atlantic salmon Salmo salar L. farmed in south-east Tasmania, Australia, are susceptible to infection by the Tasmanian Rickettsia-like organism (TRLO), a Gram-negative bacterium. Here, we report the first isolation of TRLO from south-east Tasmania in pure culture and show that the bacterium is culturable on both specialised enriched agar and in cell culture using the CHSE-214 cell line. In vitro cultured TRLO was used to reproducibly elicit disease in Atlantic salmon parr held in fresh water. In inoculated fish, TRLO was observed intracytoplasmically in peripheral blood leucocytes, suggesting that these cells are responsible for haematogenous dispersal of the bacterium within the host. Fish with experimentally induced disease presented with gross and histopathological changes similar to TRLO-infected fish at commercial marine farms. TRLO was also isolated in culture from farmed Atlantic salmon in the Tamar River and Macquarie Harbour production areas in Tasmania, both of which have no history of TRLO-associated disease. These TRLO isolates appear to be serologically distinct from each other as well as from isolates obtained from south-east Tasmania, linking each serotype to a specific geographical location within Tasmania. Despite the lack of clinical evidence of TRLO-linked disease in fish grown in the Tamar River and Macquarie Harbour, experimental infection trials demonstrably showed the pathogenic potential of these TRLO serovars. Together, these data provide evidence that TRLO is a fastidious, facultative intracellular bacterium and confirm TRLO as a pathogen of Atlantic salmon, causing a disease designated Tasmanian salmonid rickettsiosis.

Tetrahymena rostrata is a ciliated protist which can parasitize the gray field slug, Deroceras reticulatum. Here, we report the sequence of the mitochondrial genome (mt genome) of a strain of T. rostrata that was isolated from the egg of D. reticulatum. Whole cell genomic DNA was sequenced using Illumina® MiSeq and the mitochondrial DNA sequence reads were extracted and assembled. The resulting 47,235 bp assembly contained rRNAs, tRNAs, and 45 protein coding DNA sequences of which 21 encoded proteins of unknown function. Phylogenetic analysis showed T. rostrata clustered with Tetrahymena thermophila, Tetrahymena pigmentosa, Tetrahymena pyriformis, Tetrahymena paravorax, and Tetrahymena malaccensis.

Protein phosphorylation plays essential roles in many cellular processes. Despite recent progress in the genomics, transcriptomics and proteomics of socioeconomically important parasitic nematodes, there is scant phosphoproteomic data to underpin molecular biological discovery. Here, using the phosphopeptide enrichment-based LC-MS/MS and data-independent acquisition (DIA) quantitation, we characterised the first developmental phosphoproteome of the parasitic nematode Haemonchus contortus - one of the most pathogenic parasites of ruminant livestock. Totally, 1804 phosphorylated proteins with 4406 phosphorylation sites ('phosphosites') from different developmental stages/sexes were identified. Bioinformatic analyses of quantified 'phosphosites' exhibited distinctive stage- and sex-specific patterns during development, and identified a subset of phosphoproteins proposed to play crucial roles in processes such as spindle positioning, signal transduction and kinase activity. A sequence-based comparison of the phosphoproteome of H. contortus with those of two free-living nematode species (Caenorhabditis elegans and Pristionchus pacificus) suggested a limited number of common protein phosphorylation events among these species. Our findings infer active roles for protein phosphorylation in the adaptation of a parasitic nematode to a constantly changing external environment. The phosphoproteomic data set for H. contortus provides a basis to better understand phosphorylation and associated biological processes (e.g., regulation of signal transduction), and might enable the discovery of novel anthelmintic targets. SIGNIFICANCE: Here, we report the first phosphoproteome for a socioeconomically parasitic nematode (Haemonchus contortus). This phosphoproteome exhibits distinctive patterns during development, suggesting active roles of post-translational modification in the parasite's adaptation to changing environments within and outside of the host animal. This work sheds a light on the developmental phosphorylation in a parasitic nematode, and could enable the discovery of novel interventions against major pathogens.

This study aimed to elucidate virus, host and environmental dynamics of Vietnamese H5 highly pathogenic avian influenza viruses (HPAIVs) during 2014–2017. Epidemiologically, H5 HPAIVs were frequently detected in apparently healthy domestic and Muscovy ducks and therefore these are preferred species for H5 HPAIV detection in active surveillance. Virologically, clade 2.3.2.1c and 2.3.4.4 H5 HPAIVs were predominant and exhibited distinct phylogeographic evolution. Clade 2.3.2.1c viruses clustered phylogenetically in North, Central and South regions, whilst clade 2.3.4.4 viruses only detected in North and Central regions formed small groups. These viruses underwent diverse reassortment with existence of at least 12 genotypes and retained typical avian-specific motifs. These H5 HPAIVs exhibited large antigenic distance from progenitor viruses and commercial vaccines currently used in poultry. Bayesian phylodynamic analysis inferred that clade 2.3.2.1c viruses detected during 2014–2017 were likely descended from homologous clade viruses imported to Vietnam previously and/or preexisting Chinese viruses during 2012–2013. Vietnamese clade 2.3.4.4 viruses closely shared genetic traits with contemporary foreign spillovers, suggesting that there existed multiple transboundary virus dispersals to Vietnam. This study provides insights into the evolution of Vietnamese H5 HPAIVs and highlights the necessity of strengthening control measures such as, preventive surveillance and poultry vaccination.

Background Natural compounds from plants are known to provide a source of anthelmintic molecules. In previous studies, we have shown that plant extracts from the plant Picria fel-terrae Lour. and particular fractions thereof have activity against the free-living nematode Caenorhabditis elegans, causing quite pronounced stress responses in this nematode. We have also shown that a fraction, designated Pf-fraction 5, derived from this plant has a substantial adverse effect on this worm; however, nothing is known about the molecular processes affected in the worm. In the present study, we explored this aspect. Results Key biological processes linked to upregulated genes (n = 214) included ‘response to endoplasmic reticulum stress’ and ‘lipid metabolism’, and processes representing downregulated genes (n = 357) included ‘DNA-conformation change’ and ‘cellular lipid metabolism’. Conclusions Exposure of C. elegans to Pf-fraction 5 induces significant changes in the transcriptome. Gene ontology analysis suggests that Pf-fraction 5 induces endoplasmic reticulum and mitochondrial stress, and the changes in gene expression are either a direct or indirect consequence of this. Further work is required to assess specific responses to sub-fractions of Pf-fraction 5 in time-course experiments in C. elegans, to define the chemical(s) with potent anthelmintic properties, to attempt to unravel their mode(s) of action and to assess their selectivity against nematodes.

Background Signalling pathways have been extensively investigated in the free-living nematode Caenorhabditis elegans, but very little is known about these pathways in parasitic nematodes. Here, we constructed a model for the dauer-associated signalling pathways in an economically highly significant parasitic worm, Haemonchus contortus. Methods Guided by data and information available for C. elegans, we used extensive genomic and transcriptomic datasets to infer gene homologues in the dauer-associated pathways, explore developmental transcriptomic, proteomic and phosphoproteomic profiles in H. contortus and study selected molecular structures. Results The canonical cyclic guanosine monophosphate (cGMP), transforming growth factor-β (TGF-β), insulin-like growth factor 1 (IGF-1) and steroid hormone signalling pathways of H. contortus were inferred to represent a total of 61 gene homologues. Compared with C. elegans, H. contortus has a reduced set of genes encoding insulin-like peptides, implying evolutionary and biological divergences between the parasitic and free-living nematodes. Similar transcription profiles were found for all gene homologues between the infective stage of H. contortus and dauer stage of C. elegans. High transcriptional levels for genes encoding G protein-coupled receptors (GPCRs), TGF-β, insulin-like ligands (e.g. ins-1, ins-17 and ins-18) and transcriptional factors (e.g. daf-16) in the infective L3 stage of H. contortus were suggestive of critical functional roles in this stage. Conspicuous protein expression patterns and extensive phosphorylation of some components of these pathways suggested marked post-translational modifications also in the L3 stage. The high structural similarity in the DAF-12 ligand binding domain among nematodes indicated functional conservation in steroid (i.e. dafachronic acid) signalling linked to worm development. Conclusions Taken together, this pathway model provides a basis to explore hypotheses regarding biological processes and regulatory mechanisms (via particular microRNAs, phosphorylation events and/or lipids) associated with the development of H. contortus and related nematodes as well as parasite-host cross talk, which could aid the discovery of new therapeutic targets. Electronic supplementary material The online version of this article (10.1186/s13071-019-3419-6) contains supplementary material, which is available to authorized users.

Infection with the food-borne liver fluke Opisthorchis viverrini is the principal risk factor (IARC Working Group on the Evaluation of Carcinogenic Risks to Humans, 2012) for cholangiocarcinoma (CCA) in the Lower Mekong River Basin countries including Thailand, Lao PDR, Vietnam and Cambodia. We exploited this link to explore the role of the secreted growth factor termed liver fluke granulin (Ov-GRN-1) in pre-malignant lesions by undertaking programmed CRISPR/Cas9 knockout of the Ov-GRN-1 gene from the liver fluke genome. Deep sequencing of amplicon libraries from genomic DNA of gene-edited parasites revealed Cas9-catalyzed mutations within Ov-GRN-1. Gene editing resulted in rapid depletion of Ov-GRN-1 transcripts and the encoded Ov-GRN-1 protein. Gene-edited parasites colonized the biliary tract of hamsters and developed into adult flukes, but the infection resulted in reduced pathology as evidenced by attenuated biliary hyperplasia and fibrosis. Not only does this report pioneer programmed gene-editing in parasitic flatworms, but also the striking, clinically-relevant pathophysiological phenotype confirms the role for Ov-GRN-1 in virulence morbidity during opisthorchiasis., eLife digest In the rural regions alongside the Mekong River in South East Asia, traditional cuisines often use uncooked or under cooked fish, many of which carry a worm known as Opisthorchis viverrini. Once inside the body, this parasite settles in the human liver, causing a tropical disease known as liver fluke infection. Out of the 10 million people affected by O. viverrini, thousands will also develop a type of liver cancer that is triggered by the presence of the worm. In particular, the parasite secretes a protein known as granulin that may encourage certain liver cells to multiply, potentially raising the risk for cancer. A gene editing technique called CRISPR/Cas9 allows scientist to precisely target and then deactivate the genetic information a cell needs to produce a given protein. While the tool has been used in other species before, it was unknown if it could be applied to O. viverrini. Here, Arunsan et al. harnessed CRISPR/Cas9 to deactivate the gene that codes for granulin and create parasites that can only produce very little of the protein. Hamsters infected with the gene-edited worms had fewer symptoms of liver fluke infection compared to those carrying normal O. viverrini. The animals with parasites that cannot produce granulin also had fewer changes to the liver that are associated with cancer. These findings confirm that granulin has a role in promoting liver fluke infection and liver cancer. Alongside this work, Ittiprasert et al. used CRISPR/Cas9 to inactivate a gene in a species of worm that causes a human disease called schistosomiasis. Together, these findings demonstrate for the first time that the gene editing method can be adapted for use in parasitic worms, which are a major public health problem in tropical climates. This tool should help scientists understand how the parasites invade and damage our bodies, and provide new ideas for treatment and disease control.

Background Toxocara canis is quite closely related to Ascaris suum but its biology is more complex, involving a phase of arrested development (diapause or hypobiosis) in tissues as well as transplacental and transmammary transmission routes. In the present study, we explored and compared dauer-like signalling pathways of T. canis and A. suum to infer which components in these pathways might associate with, or regulate, this added complexity in T. canis. Methods Guided by information for Caenorhabditis elegans, we bioinformatically inferred and compared components of dauer-like signalling pathways in T. canis and A. suum using genomic and transcriptomic data sets. In these two ascaridoids, we also explored endogenous dafachronic acids (DAs), which are known to be critical in regulating larval developmental processes in C. elegans and other nematodes, by liquid chromatography-mass spectrometry (LC-MS). Results Orthologues of C. elegans dauer signalling genes were identified in T. canis (n = 55) and A. suum (n = 51), inferring the presence of a dauer-like signalling pathway in both species. Comparisons showed clear differences between C. elegans and these ascaridoids as well as between T. canis and A. suum, particularly in the transforming growth factor-β (TGF-β) and insulin-like signalling pathways. Specifically, in both A. suum and T. canis, there was a paucity of genes encoding SMAD transcription factor-related protein (daf-3, daf-5, daf-8 and daf-14) and insulin/insulin-like peptide (daf-28, ins-4, ins-6 and ins-7) homologues, suggesting an evolution and adaptation of the signalling pathway in these parasites. In T. canis, there were more orthologues coding for homologues of antagonist insulin-like peptides (Tc-ins-1 and Tc-ins-18), an insulin receptor substrate (Tc-ist-1) and a serine/threonine kinase (Tc-akt-1) than in A. suum, suggesting potentiated functional roles for these molecules in regulating larval diapause and reactivation. A relatively conserved machinery was proposed for DA synthesis in the two ascaridoids, and endogenous Δ4- and Δ7-DAs were detected in them by LC-MS analysis. Differential transcription analysis between T. canis and A. suum suggests that ins-17 and ins-18 homologues are specifically involved in regulating development and migration in T. canis larvae in host tissues. Conclusion The findings of this study provide a basis for functional explorations of insulin-like peptides, signalling hormones (i.e. DAs) and related nuclear receptors, proposed to link to development and/or parasite-host interactions in T. canis. Elucidating the functional roles of these molecules might contribute to the discovery of novel anthelmintic targets in ascaridoids. Electronic supplementary material The online version of this article (10.1186/s13071-018-3265-y) contains supplementary material, which is available to authorized users.

Background Babesiosis is a socioeconomically important tick-borne disease of animals (including humans) caused by haemoprotozoan parasites. The severity of babesiosis relates to host and parasite factors, particularly virulence/pathogenicity. Although Babesia bovis is a particularly pathogenic species of cattle, there are species of Babesia of ruminants that have limited pathogenicity. For instance, the operational taxonomic unit Babesia sp. Xinjiang (abbreviated here as Bx) of sheep from China is substantially less virulent/pathogenic than B. bovis is in cattle. Although the reason for this distinctiveness is presently unknown, it is possible that Bx has a reduced ability to adhere to cells or evade/suppress immune responses, which might relate to particular proteins, such as the variant erythrocyte surface antigens (VESAs). Results We sequenced and annotated the 8.4 Mb nuclear draft genome of Bx and compared it with those of B. bovis and B. bigemina by synteny analysis; we also investigated the genetic relationship of Bx with selected Babesia species and related apicomplexans for which genomic datasets are available, and explored the VESA complement in Bx. Conclusions The availability of the Bx genome now provides unique opportunities to elucidate aspects of the molecular biology, biochemistry and physiology of Bx, and to explore the reason(s) for its limited virulence and/or apparent ability to evade immune attack by the host animal. Moreover, the present genomic resource and an in vitro culture system for Bx raises the prospect of establishing a functional genomic platform to explore essential genes as new intervention targets against babesiosis. Electronic supplementary material The online version of this article (doi:10.1186/s13071-016-1846-1) contains supplementary material, which is available to authorized users.

Billions of people and animals are infected with parasitic worms (helminths). Many of these worms cause diseases that have a major socioeconomic impact worldwide, and are challenging to control because existing treatment methods are often inadequate. There is, therefore, a need to work toward developing new intervention methods, built on a sound understanding of parasitic worms at molecular level, the relationships that they have with their animal hosts and/or the diseases that they cause. Decoding the genomes and transcriptomes of these parasites brings us a step closer to this goal. The key focus of this article is to critically review and discuss bioinformatic tools used for the assembly and annotation of these genomes and transcriptomes, as well as various post-genomic analyses of transcription profiles, biological pathways, synteny, phylogeny, biogeography and the prediction and prioritisation of drug target candidates. Bioinformatic pipelines implemented and established recently provide practical and efficient tools for the assembly and annotation of genomes of parasitic worms, and will be applicable to a wide range of other parasites and eukaryotic organisms. Future research will need to assess the utility of long-read sequence data sets for enhanced genomic assemblies, and develop improved algorithms for gene prediction and post-genomic analyses, to enable comprehensive systems biology explorations of parasitic organisms.

Parasitic worms of the genus Trichinella (phylum Nematoda; class Enoplea) represent a complex of at least twelve taxa that infect a range of different host animals, including humans, around the world. They are foodborne, intracellular nematodes, and their life cycles differ substantially from those of other nematodes. The recent characterization of the genomes and transcriptomes of all twelve recognized taxa of Trichinella now allows, for the first time, detailed studies of their molecular biology. In the present study, we defined, curated, and compared the protein kinase complements (kinomes) of Trichinella spiralis and T. pseudospiralis using an integrated bioinformatic workflow employing transcriptomic and genomic data sets. We examined how variation in the kinome might link to unique aspects of Trichinella morphology, biology, and evolution. Furthermore, we utilized in silico structural modeling to discover and characterize a novel, MOS-like kinase with an unusual, previously undescribed N-terminal domain. Taken together, the present findings provide a basis for comparative investigations of nematode kinomes, and might facilitate the identification of Enoplea-specific intervention and diagnostic targets. Importantly, the in silico modeling approach assessed here provides an exciting prospect of being able to identify and classify currently unknown (orphan) kinases, as a foundation for their subsequent structural and functional investigation.

Advances in nucleic acid sequencing, mass spectrometry and computational biology have facilitated the identification, annotation and analysis of genes, transcripts, proteins and metabolites in model nematodes (Caenorhabditis elegans and Pristionchus pacificus) and socioeconomically important parasitic nematodes (Clades I, III, IV and V). Significant progress has been made in genomics and transcriptomics as well as in the proteomics and lipidomics of Haemonchus contortus (the barber's pole worm) - one of the most pathogenic representatives of the order Strongylida. Here, we review salient aspects of genomics, transcriptomics, proteomics, lipidomics, glycomics and functional genomics, and discuss the rise of integrative 'omics of this economically important parasite. Although our knowledge of the molecular biology, genetics and biochemistry of H. contortus and related species has progressed significantly, much remains to be explored, particularly in areas such as drug resistance, unique/unknown genes, host-parasite interactions, parasitism and the pathogenesis of disease, by integrating the use of multiple 'omics methods. This approach should lead to a better understanding of H. contortus and its relatives at a 'systems biology' level, and should assist in developing new interventions against these parasites.

The surface tegument of Fasciola hepatica is a crucial tissue due to its key role at the host-parasite interface. We characterised three novel proteins, termed Fhteg1, Fhteg5 and Fhteg8, that are found in the tegument membrane fraction of adult F. hepatica. Bioinformatic analysis of proteomic datasets identified Fhteg5 and Fhteg8 as tegumental glycoproteins and revealed that Fhteg1, Fhteg5 and Fhteg8 are associated with exosomes of adult F. hepatica. Fhteg1, Fhteg5 and Fhteg8 appear to be related to uncharacterised sequences in F. gigantica, Fasciolopsis buski, Echinostoma caproni, Clonorchis sinensis, Opisthorchis viverrini, Schistosoma japonicum and S. mansoni, although F. hepatica appears to have expanded this family. Fhteg1 and Fhteg5 were characterised in detail. The Fhteg1 and Fhteg5 gene transcripts each demonstrate significant upregulation in juvenile fluke 2-4 days post-excystment, with transcript levels maintained during development over 3 weeks in vitro. RNAseq data showed that both Fhtegs are expressed in the adult life stage, although the transcript levels were about 8 fold lower than those in juveniles (3 week post infection). Using immunocytochemistry, Fhteg1 and Fhteg5 were each shown to be expressed in cells adjacent to the muscle layer as well as on the surface of 1 week old juveniles, whilst Fhteg5 was also present in cells at the base of the pharynx. RNAi mediated knockdown of Fhteg1 and Fhteg5 transcripts in 4-10 day old juveniles had no effect on parasite survival, movement or growth in vitro. Although no IgG responses were observed for Fhteg1 or Fhteg5 during infection in sheep and cattle, both proteins elicited a low IgG response in a proportion of infected rats. Rats vaccinated with Fhteg1 and Fhteg5 showed good IgG responses to both proteins and a mean 48.2 % reduction in worm burden following parasite challenge. Although vaccination of cattle with both proteins induced a range of IgG responses, no protection was observed against parasite challenge. This is the first study to provide insights into the molecular properties of two novel, developmentally regulated surface tegument proteins in F. hepatica.

Background Mitochondrial genomes provide useful genetic markers for systematic and population genetic studies of parasitic helminths. Although many such genome sequences have been published and deposited in public databases, there is evidence that some of them are incomplete relating to an inability of conventional techniques to reliably sequence non-coding (repetitive) regions. In the present study, we characterise the complete mitochondrial genome—including the long, non-coding region—of the carcinogenic Chinese liver fluke, Clonorchis sinensis, using long-read sequencing. Methods The mitochondrial genome was sequenced from total high molecular-weight genomic DNA isolated from a pool of 100 adult worms of C. sinensis using the MinION sequencing platform (Oxford Nanopore Technologies), and assembled and annotated using an informatic approach. Results From > 93,500 long-reads, we assembled a 18,304 bp-mitochondrial genome for C. sinensis. Within this genome we identified a novel non-coding region of 4,549 bp containing six tandem-repetitive units of 719–809 bp each. Given that genomic DNA from pooled worms was used for sequencing, some variability in length/sequence in this tandem-repetitive region was detectable, reflecting population variation. Conclusions For C. sinensis, we report the complete mitochondrial genome, which includes a long (> 4.5 kb) tandem-repetitive region. The discovery of this non-coding region using a nanopore-sequencing/informatic approach now paves the way to investigating the nature and extent of length/sequence variation in this region within and among individual worms, both within and among C. sinensis populations, and to exploring whether this region has a functional role in the regulation of replication and transcription, akin to the mitochondrial control region in mammals. Although applied to C. sinensis, the technological approach established here should be broadly applicable to characterise complex tandem-repetitive or homo-polymeric regions in the mitochondrial genomes of a wide range of taxa., Author summary In the present study, we characterised the complete mitochondrial genome of Clonorchis sinensis—a carcinogenic liver fluke. To do this, we sequenced from total genomic DNA from multiple adult worms using a new method (Oxford Nanopore technology) to obtain data for long stretches of DNA, and then assembled these data to construct a mitochondrial genome of 18,304 bp, containing a > 4.5 kb-long tandem-repetitive region—not previously detected in this species. The results demonstrate that this method is effective at sequencing long and complex non-coding elements—not achievable using conventional techniques. The discovery of this long tandem-repetitive region in C. sinensis provides an opportunity to now explore its origin(s) and length/sequence diversity in populations of this species, and also to characterise its function(s). The technological approach employed here should have broad applicability to characterise previously-elusive non-coding mitochondrial genomic regions in a wide range of taxa.

Background Human schistosomiasis is a neglected tropical disease caused by parasitic worms of the genus Schistosoma that still affects some 200 million people. The mainstay of schistosomiasis control is a single drug, praziquantel. The reliance on this drug carries a risk of resistance emerging to this anthelmintic, such that research towards alternative anti-schistosomal drugs is warranted. In this context, a number of studies have employed computational approaches to prioritise proteins for investigation as drug targets, based on extensive genomic, transcriptomic and small-molecule data now available. Methods Here, we established a customisable, online application for the prioritisation of drug targets and applied it, for the first time, to the entire inferred proteome of S. haematobium. This application enables selection of weighted and ranked proteins representing potential drug targets, and integrates transcriptional data, orthology and gene essentiality information as well as drug-drug target associations and chemical properties of predicted ligands. Results Using this application, we defined 25 potential drug targets in S. haematobium that associated with approved drugs, and 3402 targets that (although they could not be linked to any compounds) are conserved among a range of socioeconomically important flatworm species and might represent targets for new trematocides. Conclusions The online application developed here represents an interactive, customisable, expandable and reproducible drug target ranking and prioritisation approach that should be useful for the prediction of drug targets in schistosomes and other species of parasitic worms in the future. We have demonstrated the utility of this online application by predicting potential drug targets in S. haematobium that can now be evaluated using functional genomics tools and/or small molecules, to establish whether they are indeed essential for parasite survival, and to assist in the discovery of novel anti-schistosomal compounds. Electronic supplementary material The online version of this article (10.1186/s13071-018-3197-6) contains supplementary material, which is available to authorized users.

An unknown member of the family Pasteurellaceae was repeatedly isolated from 20- to 24-week-old pigs with severe pulmonary lesions reared on the same farm in Victoria, Australia. The etiological diagnosis of the disease was inconclusive. The complete genome sequence analysis of one strain, 15-184, revealed some phylogenic proximity to Glaesserella ( Haemophilus ) parasuis , the cause of Glasser's disease. However, the sequences of the 16S rRNA and housekeeping genes, as well as the average nucleotide identity scores, differed from those of all other known species in the family Pasteurellaceae . The protein content of 15-184 was composite, with 60% of coding sequences matching known G. parasuis products, while more than 20% had a closer relative in the genera Actinobacillus , Mannheimia , Pasteurella , and Bibersteinia . Several putative virulence genes absent from G. parasuis but present in other Pasteurellaceae were also found, including the apxIII RTX toxin gene from Actinobacillus pleuropneumoniae , ABC transporters from Actinobacillus minor , and iron transporters from various species. Three prophages and one integrative conjugative element were present in the isolate. Horizontal gene transfers might explain the mosaic genomic structure and atypical metabolic and virulence characteristics of 15-184. This organism has not been assigned a taxonomic position in the family, but this study underlines the need for a large-scale epidemiological and clinical characterization of this novel pathogen in swine populations, as a genomic analysis suggests it could have a severe impact on pig health. IMPORTANCE Several species of Pasteurellaceae cause a range of significant diseases in pigs. A novel member of this family was recently isolated from Australian pigs suffering from severe respiratory infections. Comparative whole-genome analyses suggest that this bacterium represents a new species, which possesses a number of virulence genes horizontally acquired from a diverse range of other Pasteurellaceae . While the possible contribution of other coinfecting noncultivable agents to the disease has not been ruled out in this study, the repertoire of virulence genes found in this organism may nevertheless explain some aspects of the associated pathology observed on the farm. The prevalence of this novel pathogen within pig populations is currently unknown. This finding is of particular importance for the pig industry, as this organism can have a serious impact on the health of these animals.

Background Blood flukes of the genus Schistosoma cause schistosomiasis—a neglected tropical disease (NTD) that affects more than 200 million people worldwide. Studies of schistosome genomes have improved our understanding of the molecular biology of flatworms, but most of them have focused largely on protein-coding genes. Small non-coding RNAs (sncRNAs) have been explored in selected schistosome species and are suggested to play essential roles in the post-transcriptional regulation of genes, and in modulating flatworm-host interactions. However, genome-wide small RNA data are currently lacking for key schistosomes including Schistosoma haematobium—the causative agent of urogenital schistosomiasis of humans. Methodology MicroRNAs (miRNAs) and other sncRNAs of male and female adults of S. haematobium and small RNA transcription levels were explored by deep sequencing, genome mapping and detailed bioinformatic analyses. Principal findings In total, 89 transcribed miRNAs were identified in S. haematobium—a similar complement to those reported for the congeners S. mansoni and S. japonicum. Of these miRNAs, 34 were novel, with no homologs in other schistosomes. Most miRNAs (n = 64) exhibited sex-biased transcription, suggestive of roles in sexual differentiation, pairing of adult worms and reproductive processes. Of the sncRNAs that were not miRNAs, some related to the spliceosome (n = 21), biogenesis of other RNAs (n = 3) or ribozyme functions (n = 16), whereas most others (n = 3798) were novel (‘orphans’) with unknown functions. Conclusions This study provides the first genome-wide sncRNA resource for S. haematobium, extending earlier studies of schistosomes. The present work should facilitate the future curation and experimental validation of sncRNA functions in schistosomes to enhance our understanding of post-transcriptional gene regulation and of the roles that sncRNAs play in schistosome reproduction, development and parasite-host cross-talk., Author summary Human schistosomiasis is a chronic, neglected tropical disease (NTD) that is predominantly caused by the blood flukes Schistosoma haematobium, S. mansoni and S. japonicum. Infections by S. haematobium and/or S. mansoni are highly prevalent in Africa, affecting ~ 200 million people. The decoding of schistosome draft genomes has, to some extent, improved our understanding of the molecular biology of these parasites and now allows for non-protein-coding regions in these genomes to be characterised. Here, we explored small RNAs in adult S. haematobium by deep sequencing, reference genome mapping and detailed bioinformatic analyses. This study provides the first genome-wide miRNA and sncRNA resource for S. haematobium, extending earlier work on schistosomes and facilitating future curation efforts and functional investigations of schistosome sncRNAs. These efforts should enable a better understanding of post-transcriptional RNA modifications, gene regulation and novel aspects of parasite development, parasite-host cross-talk and disease at the molecular level.

Clonorchis sinensis (family Opisthorchiidae) is an important foodborne parasite that has a major socioeconomic impact on ~35 million people predominantly in China, Vietnam, Korea and the Russian Far East. In humans, infection with C. sinensis causes clonorchiasis, a complex hepatobiliary disease that can induce cholangiocarcinoma (CCA), a malignant cancer of the bile ducts. Central to understanding the epidemiology of this disease is knowledge of genetic variation within and among populations of this parasite. Although most published molecular studies seem to suggest that C. sinensis represents a single species, evidence of karyotypic variation within C. sinensis and cryptic species within a related opisthorchiid fluke (Opisthorchis viverrini) emphasise the importance of studying and comparing the genes and genomes of geographically distinct isolates of C. sinensis. Recently, we sequenced, assembled and characterised a draft nuclear genome of a C. sinensis isolate from Korea and compared it with a published draft genome of a Chinese isolate of this species using a bioinformatic workflow established for comparing draft genome assemblies and their gene annotations. We identified that 50.6% and 51.3% of the Korean and Chinese C. sinensis genomic scaffolds were syntenic, respectively. Within aligned syntenic blocks, the genomes had a high level of nucleotide identity (99.1%) and encoded 15 variable proteins likely to be involved in diverse biological processes. Here, we review current technical challenges of using draft genome assemblies to undertake comparative genomic analyses to quantify genetic variation between isolates of the same species. Using a workflow that overcomes these challenges, we report on a high-quality draft genome for C. sinensis from Korea and comparative genomic analyses, as a basis for future investigations of the genetic structures of C. sinensis populations, and discuss the biotechnological implications of these explorations.

Opisthorchiasis is a neglected tropical disease of major proportion, caused by the carcinogenic, Asian liver fluke, Opisthorchis viverrini. This hepatobiliary disease is known to be associated with malignant cancer (cholangiocarcinoma, CCA) and affects millions of people in Southeast Asia. No vaccine is available, and only one drug (praziquantel) is routinely employed against the parasite. Despite technological advances, little is known about the molecular biology of the fluke itself and the disease complex that it causes in humans. The advent of high-throughput nucleic acid sequencing and bioinformatic technologies is enabling researchers to gain global insights into the molecular pathways and processes in parasites. The principal aims of this chapter are to (1) review molecular research of O. viverrini and opisthorchiasis; (2) provide an account of recent advances in the sequencing and characterization of the genome and transcriptomes of O. viverrini; (3) describe the complex life of this worm in the biliary system of the definitive (human) host and how the fluke interacts with this host and causes disease at the molecular level; (4) discuss the implications of systems biological research and (5) consider how progress in genomics and informatics might enable explorations of O. viverrini and related worms and the discovery of new interventions against opisthorchiasis and CCA.

Parasitic trematodes (flukes) cause substantial mortality and morbidity in humans. The Chinese liver fluke, Clonorchis sinensis, is one of the most destructive parasitic worms in humans in China, Vietnam, Korea and the Russian Far East. Although C. sinensis infection can be controlled relatively well using anthelmintics, the worm is carcinogenic, inducing cholangiocarcinoma and causing major suffering in ~15 million people in Asia. This chapter provides an account of C. sinensis and clonorchiasis research-covering aspects of biology, epidemiology, pathogenesis and immunity, diagnosis, treatment and control, genetics and genomics. It also describes progress in the area of molecular biology (genetics, genomics, transcriptomics and proteomics) and highlights challenges associated with comparative genomics and population genetics. It then reviews recent advances in the sequencing and characterisation of the mitochondrial and nuclear genomes for a Korean isolate of C. sinensis and summarises salient comparative genomic work and the implications thereof. The chapter concludes by considering how advances in genomic and informatics will enable research on the genetics of C. sinensis and related parasites, as well as the discovery of new fluke-specific intervention targets.

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.

Background The accurate tracking of Cryptosporidium in faecal, water and/or soil samples in water catchment areas is central to developing strategies to manage the potential risk of cryptosporidiosis transmission to humans. Various PCR assays are used for this purpose. Although some assays achieve specific amplification from Cryptosporidium DNA in animal faecal samples, some do not. Indeed, we have observed non-specificity of some oligonucleotide primers in the small subunit of nuclear ribosomal RNA gene (SSU), which has presented an obstacle to the identification and classification of Cryptosporidium species and genotypes (taxa) from faecal samples. Results Using a novel bioinformatic approach, we explored all available Cryptosporidium genome sequences for new and diagnostically-informative, multi-copy regions to specifically design oligonucleotide primers in the large subunit of nuclear ribosomal RNA gene (LSU) as a basis for an effective nested PCR-based sequencing method for the identification and/or classification of Cryptosporidium taxa. Conclusion This newly established PCR, which has high analytical specificity and sensitivity, is now in routine use in our laboratory, together with other assays developed by various colleagues. Although the present bioinformatic workflow used here was for the specific design of primers in nuclear DNA of Cryptosporidium, this approach should be broadly applicable to many other microorganisms. Electronic supplementary material The online version of this article (10.1186/s13071-017-2462-4) contains supplementary material, which is available to authorized users.

Throughout history, the life sciences have been revolutionised by technological advances; in our era this is manifested by advances in instrumentation for data generation, and consequently researchers now routinely handle large amounts of heterogeneous data in digital formats. The simultaneous transitions towards biology as a data science and towards a ‘life cycle’ view of research data pose new challenges. Researchers face a bewildering landscape of data management requirements, recommendations and regulations, without necessarily being able to access data management training or possessing a clear understanding of practical approaches that can assist in data management in their particular research domain.Here we provide an overview of best practice data life cycle approaches for researchers in the life sciences/bioinformatics space with a particular focus on ‘omics’ datasets and computer-based data processing and analysis. We discuss the different stages of the data life cycle and provide practical suggestions for useful tools and resources to improve data management practices.

Roundworms belong to a diverse phylum (Nematoda) which is comprised of many parasitic species including whipworms (genus Trichuris). These worms have adapted to a biological niche within the host and exhibit unique morphological characteristics compared with other nematodes. Although these adaptations are known, the underlying molecular mechanisms remain elusive. The availability of genomes and transcriptomes of some whipworms now enables detailed studies of their molecular biology. Here, we defined and curated the full complement of an important class of enzymes, the protein kinases (kinomes) of two species of Trichuris, using an advanced and integrated bioinformatic pipeline. We investigated the transcription of Trichuris suis kinase genes across developmental stages, sexes and tissues, and reveal that selectively transcribed genes can be linked to central roles in developmental and reproductive processes. We also classified and functionally annotated the curated kinomes by integrating evidence from structural modelling and pathway analyses, and compared them with other curated kinomes of phylogenetically diverse nematode species. Our findings suggest unique adaptations in signalling processes governing worm morphology and biology, and provide an important resource that should facilitate experimental investigations of kinases and the biology of signalling pathways in nematodes.

Flystrike, or cutaneous myiasis, is caused by blow fly larvae of the genus Lucilia. This disease is a major problem in countries with large sheep populations. In Australia, Lucilia cuprina (Wiedemann, 1830) is the principal fly involved in flystrike. While much research has been conducted on L. cuprina, including physical, chemical, immunological, genetic and biological investigations, the molecular biology of this fly is still poorly understood. The recent sequencing, assembly and annotation of the draft genome and analyses of selected transcriptomes of L. cuprina have given a first global glimpse of its molecular biology and insights into host-fly interactions, insecticide resistance genes and intervention targets. The present article introduces L. cuprina, flystrike and associated issues, details past control efforts and research foci, reviews salient aspects of the L. cuprina genome project and discusses how the new genomic and transcriptomic resources for this fly might accelerate fundamental molecular research of L. cuprina towards developing new methods for the treatment and control of flystrike.

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.

Infectious endometritis is a major cause of reduced pregnancy rates in horses. The objectives of this study were to establish a timeline of the innate immune response in the uterus of healthy horses and to investigate the oestrous cycle effect on this. Endometrial biopsies were collected from five horses before and at 3, 12, 24, 48 and 72 h after inoculation of Escherichia coli, once in oestrus and once in dioestrus. They were analysed by quantitative real-time PCR, microbiology and histology. Neutrophil numbers increased from very low levels in the absence of inflammation to severe neutrophilia 3 h after inoculation. The concentrations of mRNAs for Toll-like receptor (TLR)2, TLR4, NOD-like receptor NLRC5, tissue inhibitor of metallopeptidases 1 (TIMP1) and chemokines CCL2, CXCL9, CXCL10 and CXCL11 were all increased 3 h after inoculation of E. coli compared to levels detected prior to inoculation. Chemokine mRNA levels remained elevated for 48 h. Concentrations of mRNAs for the antimicrobial peptides equine β-defensin 1 (EBD1), lysozyme, secretory leukoprotease inhibitor (SLPI), lipocalin 2 (LCN2), lactoferrin and uteroferrin were increased between 3 and 12 h post inoculation. The gene for secreted phospholipase A2 (sPLA2) was expressed constitutively. P19 uterocalin mRNA levels were higher in dioestrus than in oestrus over the first 24 h of inflammation. Neutrophils and many innate immune genes responded rapidly to the introduction of E. coli into the uterus, while the oestrous cycle stage had only a relatively minor effect on the response to E. coli. This study has delineated a useful model of innate immunity in infectious endometritis of healthy animals. Electronic supplementary material The online version of this article (doi:10.1186/s13567-016-0398-x) contains supplementary material, which is available to authorized users.

Parasitic nematodes of humans, animals and plants have a major, adverse impact on global health and agricultural production worldwide. To cope with their surrounding environment in and the immune attack from the host, excretory-secretory (ES) proteins are released by nematodes to orchestrate or regulate parasite-host interactions. In the present study, we characterised the ES products from short-term (12 h) in vitro culture of different developmental stages/sexes of Haemonchus contortus (one of the most important parasitic nematodes of livestock animals worldwide) using a high throughput tandem mass-spectrometry, underpinned by the most recent genomic dataset. In total, 878 unique proteins from key developmental stages/sexes (third-stage and fourth-stage larvae, and female and male adults) were identified and quantified with high confidence. Bioinformatic analyses showed noteworthy ES protein alterations during the transition from the free-living to the parasitic phase, especially for proteins which are likely involved in nutrient digestion and acquisition as well as parasite-host interactions, such as proteolytic cascade-related peptidases, glycoside hydrolases, C-type lectins and sperm-coating protein/Tpx/antigen 5/pathogenesis related-1/Sc7 (= SCP/TAPS) proteins. Our findings provide an avenue to better explore interactive processes between the host and this highly significant parasitic nematode, to underpin the search for novel drug and vaccine targets. SIGNIFICANCE: The present study represents a comprehensive proteomic analysis of the secretome of key developmental stages/sexes of H. contortus maintained in short-term in vitro culture. High throughput LC-MS/MS analysis of ES products allowed the identification of a large repertoire of proteins (secretome) and the establishment of a new proteomic database for H. contortus. The secretome of H. contortus undergoes substantial changes during the nematode's transition from free-living to parasitic stages, suggesting a constant adaptation to different environments outside of and within the host animal. Understanding the host-parasite relationship at the molecular level could assist significantly in the development of intervention strategies (i.e. novel drugs and vaccines) against H. contortus and related nematodes.

The trehalose biosynthetic pathway is of great interest for the development of novel therapeutics because trehalose is an essential disaccharide in many pathogens but is neither required nor synthesized in mammalian hosts. As such, trehalose-6-phosphate phosphatase (TPP), a key enzyme in trehalose biosynthesis, is likely an attractive target for novel chemotherapeutics. Based on a survey of genomes from a panel of parasitic nematodes and bacterial organisms and by way of a structure-based amino acid sequence alignment, we derive the topological structure of monoenzyme TPPs and classify them into 3 groups. Comparison of the functional roles of amino acid residues located in the active site for TPPs belonging to different groups reveal nuanced variations. Because current literature on this enzyme family shows a tendency to infer functional roles for individual amino acid residues, we investigated the roles of the strictly conserved aspartate tetrad in TPPs of the nematode Brugia malayi by using a conservative mutation approach. In contrast to aspartate-213, the residue inferred to carry out the nucleophilic attack on the substrate, we found that aspartate-215 and aspartate-428 of BmTPP are involved in the chemistry steps of enzymatic hydrolysis of the substrate. Therefore, we suggest that homology-based inference of functionally important amino acids by sequence comparison for monoenzyme TPPs should only be carried out for each of the 3 groups.-Cross, M., Lepage, R., Rajan, S., Biberacher, S., Young, N. D., Kim, B.-N., Coster, M. J., Gasser, R. B., Kim, J.-S., Hofmann, A. Probing function and structure of trehalose-6-phosphate phosphatases from pathogenic organisms suggests distinct molecular groupings.

The apicoplast (ap) is a unique, non-photosynthetic organelle found in most apicomplexan parasites. Due to the essential roles that this organelle has, it has been widely considered as target for drugs against diseases caused by apicomplexans. Exploring the ap genomes of such parasites would provide a better understanding of their systematics and their basic molecular biology for therapeutics. However, there is limited information available on the ap genomes of apicomplexan parasites. In the present study, the ap genomes of two operational taxonomic units of Babesia (known as Babesia sp. Lintan [Bl] and Babesia sp. Xinjiang [Bx]) from sheep were sequenced, assembled and annotated using a massive parallel sequencing-based approach. Then, the gene content and gene order in these ap genomes (∼30.7kb in size) were defined and compared, and the genetic differences were assessed. In addition, a phylogenetic analysis of ap genomic data sets was carried out to assess the relationships of these taxonomic units with other apicomplexan parasites for which complete ap genomic data sets were publicly available. The results showed that the ap genomes of Bl and Bx encode 59 and 57 genes, respectively, including 2 ribosomal RNA genes, 25 transfer RNA genes and 30-32 protein-encoding genes, being similar in content to those of Babesia bovis and B. orientalis. Ap gene regions that might serve as markers for future epidemiological and population genetic studies of Babesia species were identified. Using sequence data for a subset of six protein-encoding genes, a close relationship of Bl and Bx with Babesia bovis from cattle and B. orientalis from water buffalo was inferred. Although the focus of the present study was on Babesia, we propose that the present sequencing-bioinformatic approach should be applicable to organellar genomes of a wide range of apicomplexans of veterinary importance.

Here, we sequenced, assembled and annotated the mitochondrial (mt) genomes of two operational taxonomic units of Babesia from sheep from China using a deep sequencing-coupled approach. Then, we defined and compared the gene order of these mt genomes (~5.8 to 6.2kb in size), assessed sequence differences in mt genes among Babesia taxa and evaluated genetic relationships among these taxa and related apicomplexans (Theileria) for which mt genomic data sets were available. We also identified mt genetic regions that might be useful as markers for future population genetic and molecular epidemiological studies of Babesia from small ruminants. We propose that the sequencing-bioinformatic approach used here should be applicable to a wide range of protists of veterinary importance.

Background Toxocara canis is the causative agent of toxocariasis of humans and other animals. This parasitic nematode (roundworm) has a complex life cycle, in which substantial developmental changes and switches occur. As small non-coding RNAs (sRNAs) are key regulators of gene expression in a wide range of organisms, we explored these RNAs in T. canis to provide a basis for future studies of its developmental biology as well as host interactions and disease at the molecular level. Methods We conducted high-throughput RNA sequencing and bioinformatic analyses to define sRNAs in individual male and female adults of T. canis. Results Apart from snRNA and snoRNA, 560 and 619 microRNAs (miRNAs), including 5 and 2 novel miRNAs, were identified in male and female worms, respectively, without piRNAs being detected in either sex. An analysis of transcriptional profiles showed that, of 564 miRNAs predicted as being differentially transcribed between male and female individuals of T. canis, 218 miRNAs were transcribed exclusively in male and 277 in female worms. Functional enrichment analysis predicted that both male and female miRNAs were mainly involved in regulating embryonic morphogenesis, hemidesmosome assembly and genetic information processing. The miRNAs differentially transcribed between the sexes were predicted to be associated with sex determination, embryonic morphogenesis and nematode larval development. The roles of miRNAs were predicted based on gene ontology (GO) and KEGG pathway annotations. The miRNAs Tc-miR-2305 and Tc-miR-6090 are proposed to have roles in reproduction, embryo development and larval development, and Tc-let-7-5p, Tc-miR-34 and Tc-miR-100 appear to be involved in host-parasite interactions. Together with published information from previous studies, some miRNAs (such as Tc-miR-2861, Tc-miR-2881 and Tc-miR-5126) are predicted to represent drug targets and/or associated with drug resistance. Conclusions This is the first exploration of miRNAs in T. canis, which could provide a basis for fundamental investigations of the developmental biology of the parasite, parasite-host interactions and toxocariasis as well as applied areas, such as the diagnosis of infection/disease, drug target discovery and drug resistance detection. Electronic supplementary material The online version of this article (doi:10.1186/s13071-016-1508-3) contains supplementary material, which is available to authorized users.

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.

Trichinellosis is a globally important food-borne parasitic disease of humans caused by roundworms of the Trichinella complex. Extensive biological diversity is reflected in substantial ecological and genetic variability within and among Trichinella taxa, and major controversy surrounds the systematics of this complex. Here we report the sequencing and assembly of 16 draft genomes representing all 12 recognized Trichinella species and genotypes, define protein-coding gene sets and assess genetic differences among these taxa. Using thousands of shared single-copy orthologous gene sequences, we fully reconstruct, for the first time, a phylogeny and biogeography for the Trichinella complex, and show that encapsulated and non-encapsulated Trichinella taxa diverged from their most recent common ancestor ∼21 million years ago (mya), with taxon diversifications commencing ∼10−7 mya., Trichinellosis is a globally important food-borne disease caused by roundworms of the Trichinella complex. Here the authors present genomic sequences representing all 12 recognized Trichinella species and genotypes, and reconstruct their phylogeny and biogeography.

Parasitic roundworms (nematodes) cause substantial mortality and morbidity in animals globally. The barber's pole worm, Haemonchus contortus, is one of the most economically significant parasitic nematodes of small ruminants worldwide. Although this and related nematodes can be controlled relatively well using anthelmintics, resistance against most drugs in common use has become a major problem. Until recently, almost nothing was known about the molecular biology of H. contortus on a global scale. This chapter gives a brief background on H. contortus and haemonchosis, immune responses, vaccine research, chemotherapeutics and current problems associated with drug resistance. It also describes progress in transcriptomics before the availability of H. contortus genomes and the challenges associated with such work. It then reviews major progress on the two draft genomes and developmental transcriptomes of H. contortus, and summarizes their implications for the molecular biology of this worm in both the free-living and the parasitic stages of its life cycle. The chapter concludes by considering how genomics and transcriptomics can accelerate research on Haemonchus and related parasites, and can enable the development of new interventions against haemonchosis.

Parasitic worms, such as flatworms (platyhelminths) and roundworms (nematodes), cause substantial morbidity and mortality in animals and people globally. The ascaridoid nematode Toxocara canis is a zoonotic parasite of socioeconomic significance worldwide. In humans, this worm causes toxocariasis (disease) mainly in underprivileged communities in both the developed and developing worlds. While reasonably well studied from clinical and epidemiological perspectives, little is understood about the molecular biology of T. canis, its relationship with its hosts and the disease that it causes. However, a recent report of the draft genome and transcriptomes of T. canis should underpin many fundamental and applied research areas in the future. The present article gives a background on Toxocara and toxocariasis, a brief account of diagnostic approaches for specific identification and genetic analysis, and gives a perspective on the impact that the genome of T. canis and advanced molecular technologies could have on our understanding of the parasite and the diseases that it causes as well as the design of new and improved approaches for the diagnosis, treatment and control of toxocariasis.

The blow fly, Lucilia cuprina (Wiedemann, 1830) is a parasitic insect of major global economic importance. Maggots of this fly parasitize the skin of animal hosts, feed on excretions and tissues, and cause severe disease (flystrike or myiasis). Although there has been considerable research on L. cuprina over the years, little is understood about the molecular biology, biochemistry and genetics of this parasitic fly, as well as its relationship with its hosts and the disease that it causes. This situation might change with the recent report of the draft genome and transcriptome of this blow fly, which has given new and global insights into its biology, interactions with the host animal and aspects of insecticide resistance at the molecular level. This genomic resource will likely enable many fundamental and applied research areas in the future. The present article gives a background on L. cuprina and myiasis, a brief account of past and current treatment, prevention and control approaches, and provides a perspective on the impact that the L. cuprina genome should have on future research of this and related parasitic flies, and the design of new and improved interventions for myiasis.

Although substantial research has been focused on the 'hidden antigen' H11 of Haemonchus contortus as a vaccine against haemonchosis in small ruminants, little is know about this and related aminopeptidases. In the present article, we reviewed genomic and transcriptomic data sets to define, for the first time, the complement of aminopeptidases (designated Hc-AP-1 to Hc-AP-13) of the family M1 with homologues in Caenorhabditis elegans, characterised by zinc-binding (HEXXH) and exo-peptidase (GAMEN) motifs. The three previously published H11 isoforms (accession nos. X94187, FJ481146 and AJ249941) had most sequence similarity to Hc-AP-2 and Hc-AP-8, whereas unpublished isoforms (accession nos. AJ249942 and AJ311316) were both most similar to Hc-AP-3. The aminopeptidases characterised here had homologues in C. elegans. Hc-AP-1 to Hc-AP-8 were most similar in amino acid sequence (28-41%) to C. elegans T07F10.1; Hc-AP-9 and Hc-AP-10 to C. elegans PAM-1 (isoform b) (53-54% similar); Hc-AP-11 and Hc-AP-12 to C. elegans AC3.5 and Y67D8C.9 (26% and 50% similar, respectively); and Hc-AP-13 to C. elegans C42C1.11 and ZC416.6 (50-58% similar). Comparative analysis suggested that Hc-AP-1 to Hc-AP-8 play roles in digestion, metabolite excretion, neuropeptide processing and/or osmotic regulation, with Hc-AP-4 and Hc-AP-7 having male-specific functional roles. The analysis also indicated that Hc-AP-9 and Hc-AP-10 might be involved in the degradation of cyclin (B3) and required to complete meiosis. Hc-AP-11 represents a leucyl/cystinyl aminopeptidase, predicted to have metallopeptidase and zinc ion binding activity, whereas Hc-AP-12 likely encodes an aminopeptidase Q homologue also with these activities and a possible role in gonad function. Finally, Hc-AP-13 is predicted to encode an aminopeptidase AP-1 homologue of C. elegans with hydrolase activity, suggested to operate, possibly synergistically with a PEPT-1 ortholog, as an oligopeptide transporter in the gut for protein uptake and normal development and/or reproduction of the worm. An appraisal of structure-based amino acid sequence alignments revealed that all conceptually translated Hc-AP proteins, with the exception of Hc-AP-12, adopt a topology similar to those observed for the two subgroups of mammalian M1 aminopeptidases, which possess either three (I, II and IV) or four (I-IV) domains. In contrast, Hc-AP-12 lacks the N-terminal domain (I), but possesses a substantially expanded domain III. Although further work needs to be done to assess amino acid sequence conservation of the different aminopeptidases among individual worms within and among H. contortus populations, we hope that these insights will support future localisation, structural and functional studies of these molecules in H. contortus as well as facilitate future assessments of a recombinant subunit or cocktail vaccine against haemonchosis.