Search

Abstract Background Lymnaeid snails of the genus Austropeplea are an important vector of the liver fluke (Fasciola hepatica), contributing to livestock production losses in Australia and New Zealand. However, the species status within Austropeplea is ambiguous due to heavy reliance on morphological analysis and a relative lack of genetic data. This study aimed to characterise the mitochondrial genome of A. cf. brazieri, an intermediate host of liver fluke in eastern Victoria. Methods The mitochondrial genome was assembled and annotated from a combination of second- and third-generation sequencing data. For comparative purposes, we performed phylogenetic analyses of the concatenated nucleotide sequences of the mitochondrial protein-coding genes, cytochrome c oxidase subunit 1 and 16S genes. Results The assembled mt genome was 13,757 base pairs and comprised 37 genes, including 13 protein-coding genes, 22 transfer RNA genes and 2 ribosomal RNA genes. The mt genome length, gene order and nucleotide compositions were similar to related species of lymnaeids. Phylogenetic analyses of the mt nucleotide sequences placed A. cf. brazieri within the same clade as Orientogalba ollula with strong statistical supports. Phylogenies of the cox1 and 16S mt sequences were constructed due to the wide availability of these sequences representing the lymnaeid taxa. As expected in both these phylogenies, A. cf. brazieri clustered with other Austropeplea sequences, but the nodal supports were low. Conclusions The representative mt genome of A. cf. brazieri should provide a useful resource for future molecular, epidemiology and parasitological studies of this socio-economically important lymnaeid species. Graphical abstract

Abstract Background Filarial worms are important vector-borne pathogens of a large range of animal hosts, including humans, and are responsible for numerous debilitating neglected tropical diseases such as, lymphatic filariasis caused by Wuchereria bancrofti and Brugia spp., as well as loiasis caused by Loa loa. Moreover, some emerging or difficult-to-eliminate filarioid pathogens are zoonotic using animals like canines as reservoir hosts, for example Dirofilaria sp. ‘hongkongensis’. Diagnosis of filariasis through commonly available methods, like microscopy, can be challenging as microfilaremia may wane below the limit of detection. In contrast, conventional PCR methods are more sensitive and specific but may show limited ability to detect coinfections as well as emerging and/or novel pathogens. Use of deep-sequencing technologies obviate these challenges, providing sensitive detection of entire parasite communities, whilst also being better suited for the characterisation of rare or novel pathogens. Therefore, we developed a novel long-read metabarcoding assay for deep-sequencing the filarial nematode cytochrome c oxidase subunit I gene on Oxford Nanopore Technologies’ (ONT) MinION™ sequencer. We assessed the overall performance of our assay using kappa statistics to compare it to commonly used diagnostic methods for filarial worm detection, such as conventional PCR (cPCR) with Sanger sequencing and the microscopy-based modified Knott’s test (MKT). Results We confirmed our metabarcoding assay can characterise filarial parasites from a diverse range of genera, including, Breinlia, Brugia, Cercopithifilaria, Dipetalonema, Dirofilaria, Onchocerca, Setaria, Stephanofilaria and Wuchereria. We demonstrated proof-of-concept for this assay by using blood samples from Sri Lankan dogs, whereby we identified infections with the filarioids Acanthocheilonema reconditum, Brugia sp. Sri Lanka genotype and zoonotic Dirofilaria sp. ‘hongkongensis’. When compared to traditionally used diagnostics, such as the MKT and cPCR with Sanger sequencing, we identified an additional filarioid species and over 15% more mono- and coinfections. Conclusions Our developed metabarcoding assay may show broad applicability for the metabarcoding and diagnosis of the full spectrum of filarioids from a wide range of animal hosts, including mammals and vectors, whilst the utilisation of ONT’ small and portable MinION™ means that such methods could be deployed for field use.

Abstract Background Lucilia cuprina and L. sericata (family Calliphoridae) are globally significant ectoparasites of sheep. Current literature suggests that only one of these blowfly subspecies, L. cuprina dorsalis, is a primary parasite causing myiasis (flystrike) in sheep in Australia. These species and subspecies are difficult to distinguish using morphological features. Hence, being able to accurately identify blowflies is critical for diagnosis and for understanding their relationships with their hosts and environment. Methods In this study, adult blowflies (5 pools of 17 flies; n = 85) were collected from five locations in different states [New South Wales (NSW), Queensland (QLD), Tasmania (TAS), Victoria (VIC) and Western Australia (WA)] of Australia and their mitochondrial (mt) genomes were assembled. Results Each mt genome assembled was ~ 15 kb in size and encoded 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs and a control region. The Lucilia species mt genomes were conserved in structure, and the genes retained the same order and direction. The overall nucleotide composition was heavily biased towards As and Ts—77.7% of the whole genomes. Pairwise nucleotide diversity suggested divergence between Lucilia cuprina cuprina, L. c. dorsalis and L. sericata. Comparative analyses of these mt genomes with published data demonstrated that the blowflies collected from sheep farm in TAS clustered within a clade with L. sericata. The flies collected from an urban location in QLD were more closely related to L. sericata and represented the subspecies L. c. cuprina, whereas the flies collected from sheep farms in NSW, VIC and WA represented the subspecies L. c. dorsalis. Conclusions Phylogenetic analyses of the mt genomes representing Lucilia from the five geographic locations in Australia supported the previously demonstrated paraphyly of L. cuprina with respect to L. sericata and revealed that L. c. cuprina is distinct from L. c. dorsalis and that L. c. cuprina is more closely related to L. sericata than L. c. dorsalis. The mt genomes reported here provide an important molecular resource to develop tools for species- and subspecies-level identification of Lucilia from different geographical regions across Australia. Graphical Abstract

Abstract Background Gaining insight into molecular signalling pathways of socioeconomically important parasitic nematodes has implications for understanding their molecular biology and for developing novel anthelmintic interventions. Methods Here, we evaluated the use of a human antibody-based microarray to explore conserved elements of the signalome in the barber’s pole worm Haemonchus contortus. To do this, we prepared extracts from mixed-sex (female and male) adult worms and third-stage larvae (L3s), incubated these extracts on the antibody microarray and then measured the amounts of antibody-bound proteins (‘signal intensity’). Results In total, 878 signals were classified into two distinct categories: signals that were higher for adults than for larvae of H. contortus (n = 376), and signals that were higher for larvae than for adults of this species (n = 502). Following a data-filtering step, high confidence (‘specific’) signals were obtained for subsequent analyses. In total, 39 pan-specific signals (linked to antibodies that recognise target proteins irrespective of their phosphorylation status) and 65 phosphorylation-specific signals were higher in the adult stage, and 82 pan-specific signals and 183 phosphorylation-specific signals were higher in L3s. Thus, notably more signals were higher in L3s than in the adult worms. Using publicly available information, we then inferred H. contortus proteins that were detected (with high confidence) by specific antibodies directed against human homologues, and revealed relatively high structural conservation between the two species, with some variability for select proteins. We also in silico-matched 763 compound structures (listed in the DrugBank and Kinase SARfari public databases) to four H. contortus proteins (designated HCON_00005760, HCON_00079680, HCON_00013590 and HCON_00105100). Conclusions We conclude that the present antibody-based microarray provides a useful tool for comparative analyses of signalling pathways between/among developmental stages and/or species, as well as opportunities to explore nematocidal target candidates in H. contortus and related parasites. Graphical Abstract

Abstract 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.

Abstract 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.

Abstract Background Arthropods comprise the largest and most diverse phylum on Earth and play vital roles in nearly every ecosystem. Their diversity stems in part from variations on a conserved body plan, resulting from and recorded in adaptive changes in the genome. Dissection of the genomic record of sequence change enables broad questions regarding genome evolution to be addressed, even across hyper-diverse taxa within arthropods. Results Using 76 whole genome sequences representing 21 orders spanning more than 500 million years of arthropod evolution, we document changes in gene and protein domain content and provide temporal and phylogenetic context for interpreting these innovations. We identify many novel gene families that arose early in the evolution of arthropods and during the diversification of insects into modern orders. We reveal unexpected variation in patterns of DNA methylation across arthropods and examples of gene family and protein domain evolution coincident with the appearance of notable phenotypic and physiological adaptations such as flight, metamorphosis, sociality, and chemoperception. Conclusions These analyses demonstrate how large-scale comparative genomics can provide broad new insights into the genotype to phenotype map and generate testable hypotheses about the evolution of animal diversity.

Abstract 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.

Abstract 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.

Abstract 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.

Abstract 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.

Abstract 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.

Abstract 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.

Abstract 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.

Abstract Background Protein kinases regulate a plethora of essential signalling and other biological pathways in all eukaryotic organisms, but very little is known about them in most parasitic nematodes. Methods Here, we defined, for the first time, the entire complement of protein kinases (kinome) encoded in the barber’s pole worm (Haemonchus contortus) through an integrated analysis of transcriptomic and genomic datasets using an advanced bioinformatic workflow. Results We identified, curated and classified 432 kinases representing ten groups, 103 distinct families and 98 subfamilies. A comparison of the kinomes of H. contortus and Caenorhabditis elegans (a related, free-living nematode) revealed considerable variation in the numbers of casein kinases, tyrosine kinases and Ca2+/calmodulin-dependent protein kinases, which likely relate to differences in biology, habitat and life cycle between these worms. Moreover, a suite of kinase genes was selectively transcribed in particular developmental stages of H. contortus, indicating central roles in developmental and reproductive processes. In addition, using a ranking system, drug targets (n = 13) and associated small-molecule effectors (n = 1517) were inferred. Conclusions The H. contortus kinome will provide a useful resource for fundamental investigations of kinases and signalling pathways in this nematode, and should assist future anthelmintic discovery efforts; this is particularly important, given current drug resistance problems in parasitic nematodes.

Abstract Background Although fascioliasis has been relatively well studied, little is known about the molecular basis of this disease. This is particularly relevant, considering the very different response that sheep have to Fasciola hepatica relative to cattle. The acute phase of this disease is severe in sheep, whereas chronic fascioliasis is more common in cattle. Methods To begin to explore the host-response to Fasciola in sheep and improve the understanding of the host-pathogen interactions during the parasite’s migration through liver parenchyma to the bile duct, we used RNA sequencing (RNA-seq) to investigate livers from sheep infected for eight weeks compared with those from uninfected controls. Results This study identified 572 and 42 genes that were up- and down-regulated, respectively, in infected livers relative to uninfected controls. Our molecular findings provide significant new insights into the mechanisms linked to metabolism, fibrosis and tissue-repair in sheep, and highlight the relative importance of specific components of immune response pathways, which appear to be driven toward a suppression of inflammation. Conclusions This study is, to our knowledge, the first detailed investigation of the transcriptomic responses in the liver tissue of any host to F. hepatica infection. It defines the involvement of specific genes associated with the host’s metabolism, immune response and tissue repair/regeneration, and highlights an apparent overlapping function of many genes involved in these processes.

Abstract Background Schistosomiasis is a parasitic disease affecting ~200 million people worldwide. Schistosoma haematobium and S. mansoni are two relatively closely related schistosomes (blood flukes), and the causative agents of urogenital and hepatointestinal schistosomiasis, respectively. The availability of genomic, transcriptomic and proteomic data sets for these two schistosomes now provides unprecedented opportunities to explore their biology, host interactions and schistosomiasis at the molecular level. A particularly important group of molecules involved in a range of biological and developmental processes in schistosomes and other parasites are the G protein-coupled receptors (GPCRs). Although GPCRs have been studied in schistosomes, there has been no detailed comparison of these receptors between closely related species. Here, using a genomic-bioinformatic approach, we identified and characterised key GPCRs in S. haematobium and S. mansoni (two closely related species of schistosome). Methods Using a Hidden Markov Model (HMM) and Support Vector Machine (SVM)-based pipeline, we classified and sub-classified GPCRs of S. haematobium and S. mansoni, combined with phylogenetic and transcription analyses. Results We identified and classified classes A, B, C and F as well as an unclassified group of GPCRs encoded in the genomes of S. haematobium and S. mansoni. In addition, we characterised ligand-specific subclasses (i.e. amine, peptide, opsin and orphan) within class A (rhodopsin-like). Conclusions Most GPCRs shared a high degree of similarity and conservation, except for members of a particular clade (designated SmGPR), which appear to have diverged between S. haematobium and S. mansoni and might explain, to some extent, some of the underlying biological differences between these two schistosomes. The present set of annotated GPCRs provides a basis for future functional genomic studies of cellular GPCR-mediated signal transduction and a resource for future drug discovery efforts in schistosomes.

Abstract Background Fascioliasis is an important and neglected disease of humans and other mammals, caused by trematodes of the genus Fasciola. Fasciola hepatica and F. gigantica are valid species that infect humans and animals, but the specific status of Fasciola sp. (‘intermediate form’) is unclear. Methods Single specimens inferred to represent Fasciola sp. (‘intermediate form’; Heilongjiang) and F. gigantica (Guangxi) from China were genetically identified and characterized using PCR-based sequencing of the first and second internal transcribed spacer regions of nuclear ribosomal DNA. The complete mitochondrial (mt) genomes of these representative specimens were then sequenced. The relationships of these specimens with selected members of the Trematoda were assessed by phylogenetic analysis of concatenated amino acid sequence datasets by Bayesian inference (BI). Results The complete mt genomes of representatives of Fasciola sp. and F. gigantica were 14,453 bp and 14,478 bp in size, respectively. Both mt genomes contain 12 protein-coding genes, 22 transfer RNA genes and two ribosomal RNA genes, but lack an atp 8 gene. All protein-coding genes are transcribed in the same direction, and the gene order in both mt genomes is the same as that published for F. hepatica. Phylogenetic analysis of the concatenated amino acid sequence data for all 12 protein-coding genes showed that the specimen of Fasciola sp. was more closely related to F. gigantica than to F. hepatica. Conclusions The mt genomes characterized here provide a rich source of markers, which can be used in combination with nuclear markers and imaging techniques, for future comparative studies of the biology of Fasciola sp. from China and other countries.