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

Genome-wide RNA interference of the nhr gene family in barber's pole worm identified members crucial for larval viability in vitro

Nuclear hormone receptors (NHRs) are emerging target candidates against nematode infection and resistance. However, there is a lack of comprehensive information on NHR-coding genes in parasitic nematodes. In this study, we curated the nhr gene family for 60 major parasitic nematodes from humans and animals. Compared with the free-living model organism Caenorhabditis elegans, a remarkable contraction of the nhr family was revealed in parasitic species, with genetic diversification and conservation unveiled among nematode Clades I (10-13), III (16-42), IV (33-35) and V (25-64). Using an in vitro biosystem, we demonstrated that 40 nhr genes in a blood-feeding nematode Haemonchus contortus (clade V; barber's pole worm) were responsive to host serum and one nhr gene (i.e., nhr-64) was consistently stimulated by anthelmintics (i.e., ivermectin, thiabendazole and levamisole); Using a high-throughput RNA interference platform, we knocked down 43 nhr genes of H. contortus and identified at least two genes that are required for the viability (i.e., nhr-105) and development (i.e., nhr-17) of the infective larvae of this parasitic nematode in vitro. Harnessing this preliminary functional atlas of nhr genes for H. contortus will prime the biological studies of this gene family in nematode genetics, infection, and anthelmintic metabolism within host animals, as well as the promising discovery of novel intervention targets.

Fatty acid- and retinol-binding protein 6 does not control worm fatty acid content in Caenorhabditis elegans but might play a role in Haemonchus contortus parasitism

BACKGROUND

Nematodes have lost the ability to synthesise necessary lipids de novo and have complementally evolved the capacity to acquire fatty acids and their derivatives from a diet or host animal. Nematode-specific fatty acid- and retinol-binding protein (FAR) family is one approach that facilitates lipid acquisition, representing an Achilles heel and potential target against roundworms of socioeconomic significance. However, little is known about their detailed functional roles in either free-living or parasitic nematodes.

METHODS

A genome-wide identification and curation were performed to screen the FAR family members of Haemonchus contortus. Their transcription patterns in worms were also analysed to identify the targets. Ligand binding assay and molecular docking were conducted to verify the fatty acid binding activities of FAR proteins of interest. RNA interference (RNAi) and heterologous expression (rescuing) experiments were designed to explore the potential roles of the selected FAR protein in nematodes. Localisation of the protein was shown in sections of paraffin-embedded worms after an immunohistochemistry (IHC) assay.

RESULTS

Here, an orthologue of far-6 in the model organism Caenorhabditis elegans (Ce-far-6) was functionally characterised in a parasitic nematode, H. contortus (Hc-far-6). It is demonstrated that knockdown of Ce-far-6 gene did not affect worm fat content, reproduction, or lifespan, but decreased worm body length at an early life stage of C. elegans. In particular, the Ce-far-6 mutant associated phenotype was completely rescued by Hc-far-6, suggesting a conserved functional role. Surprisingly, there were distinct tissue expression patterns of FAR-6 in the free-living C. elegans and parasitic H. contortus. High transcriptional level of Hc-far-6 and dominant expression of FAR-6 in the intestine of the parasitic stage of H. contortus link this gene/protein to nematode parasitism.

CONCLUSIONS

These findings substantially enhance our understanding of far genes and the associated lipid biology of this important parasitic nematode at a molecular level, and the approaches established are readily applicable to the studies of far genes in a broad range of parasites.

Ubiquitination pathway model for the barber's pole worm - Haemonchus contortus

The ubiquitin-mediated pathway has been comprehensively explored in the free-living nematode Caenorhabditis elegans, but very little is known about this pathway in parasitic nematodes. Here, we inferred the ubiquitination pathway for an economically significant and pathogenic nematode - Haemonchus contortus - using abundant resources available for C. elegans. We identified 215 genes encoding ubiquitin (Ub; n = 3 genes), ubiquitin-activating enzyme (E1; one), -conjugating enzymes (E2s; 21), ligases (E3s; 157) and deubiquitinating enzymes (DUBs; 33). With reference to C. elegans, Ub, E1 and E2 were relatively conserved in sequence and structure, and E3s and DUBs were divergent, likely reflecting functional and biological uniqueness in H. contortus. Most genes encoding ubiquitination pathway components exhibit high transcription in the egg compared with other stages, indicating marked protein homeostasis in this early developmental stage. The ubiquitination pathway model constructed for H. contortus provides a foundation to explore the ubiquitin-proteasome system, crosstalk between autophagy and the proteasome system, and the parasite-host interactions. Selected E3 and DUB proteins which are very divergent in sequence and structure from host homologues or entirely unique to H. contortus and related parasitic nematodes may represent possible anthelmintic targets.

The transcriptome of Balamuthia mandrillaris trophozoites for structure-guided drug design

Balamuthia mandrillaris, a pathogenic free-living amoeba, causes cutaneous skin lesions as well as granulomatous amoebic encephalitis, a 'brain-eating' disease. As with the other known pathogenic free-living amoebas (Naegleria fowleri and Acanthamoeba species), drug discovery efforts to combat Balamuthia infections of the central nervous system are sparse; few targets have been validated or characterized at the molecular level, and little is known about the biochemical pathways necessary for parasite survival. Current treatments of encephalitis due to B. mandrillaris lack efficacy, leading to case fatality rates above 90%. Using our recently published methodology to discover potential drugs against pathogenic amoebas, we screened a collection of 85 compounds with known antiparasitic activity and identified 59 compounds that impacted the growth of Balamuthia trophozoites at concentrations below 220 µM. Since there is no fully annotated genome or proteome of B. mandrillaris, we sequenced and assembled its transcriptome from a high-throughput RNA-sequencing (RNA-Seq) experiment and located the coding sequences of the genes potentially targeted by the growth inhibitors from our compound screens. We determined the sequence of 17 of these target genes and obtained expression clones for 15 that we validated by direct sequencing. These will be used in the future in combination with the identified hits in structure guided drug discovery campaigns to develop new approaches for the treatment of Balamuthia infections.

High-quality nuclear genome for Sarcoptes scabiei-A critical resource for a neglected parasite

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.

Expanded complement of Niemann-Pick type C2-like protein genes in Clonorchis sinensis suggests functions beyond sterol binding and transport

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.

High-quality Schistosoma haematobium genome achieved by single-molecule and long-range sequencing

BACKGROUND

Schistosoma haematobium causes urogenital schistosomiasis, a neglected tropical disease affecting >100 million people worldwide. Chronic infection with this parasitic trematode can lead to urogenital conditions including female genital schistosomiasis and bladder cancer. At the molecular level, little is known about this blood fluke and the pathogenesis of the disease that it causes. To support molecular studies of this carcinogenic worm, we reported a draft genome for S. haematobium in 2012. Although a useful resource, its utility has been somewhat limited by its fragmentation.

FINDINGS

Here, we systematically enhanced the draft genome of S. haematobium using a single-molecule and long-range DNA-sequencing approach. We achieved a major improvement in the accuracy and contiguity of the genome assembly, making it superior or comparable to assemblies for other schistosome species. We transferred curated gene models to this assembly and, using enhanced gene annotation pipelines, inferred a gene set with as many or more complete gene models as those of other well-studied schistosomes. Using conserved, single-copy orthologs, we assessed the phylogenetic position of S. haematobium in relation to other parasitic flatworms for which draft genomes were available.

CONCLUSIONS

We report a substantially enhanced genomic resource that represents a solid foundation for molecular research on S. haematobium and is poised to better underpin population and functional genomic investigations and to accelerate the search for new disease interventions.

Dauer signalling pathway model for Haemonchus contortus

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.

Comparative bioinformatic analysis suggests that specific dauer-like signalling pathway components regulate Toxocara canis development and migration in the mammalian host

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.

Interactive online application for the prediction, ranking and prioritisation of drug targets in Schistosoma haematobium

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.