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

Detailed explorations of the model organisms Caenorhabditis elegans (elegant worm) and Drosophila melanogaster (vinegar fly) have substantially improved our knowledge and understanding of biological processes and pathways in metazoan organisms. Extensive functional genomic and multi-omic data sets have enabled the discovery and characterisation of ‘essential’ genes that are critical for the survival of these organisms. Recently, we showed that a machine learning (ML)-based pipeline could be utilised to predict essential genes in both C. elegans and D. melanogaster using features from DNA, RNA, protein and/or cellular data or associated information. As these distantly-related species are within the Ecdysozoa, we hypothesised that this approach could be suited for non-model organisms within the same group (phylum) of protostome animals. In the present investigation, we cross-predicted essential genes within the phylum Nematoda – between C. elegans and the parasitic filarial nematodes Brugia malayi and Onchocerca volvulus, and then ranked and prioritised these genes. Highly ranked genes were linked to key biological pathways or processes, such as ribosome biogenesis, translation and RNA processing, and were expressed at relatively high levels in the germline, gonad, hypodermis and/or nerves. The present in silico workflow is hoped to expedite the identification of drug targets in parasitic organisms for subsequent experimental validation in the laboratory.

Advances in single cell technologies are allowing investigations of a wide range of biological processes and pathways in animals, such as the multicellular model organism Caenorhabditis elegans – a free-living nematode. However, there has been limited application of such technology to related parasitic nematodes which cause major diseases of humans and animals worldwide. With no vaccines against the vast majority of parasitic nematodes and treatment failures due to drug resistance or inefficacy, new intervention targets are urgently needed, preferably informed by a deep understanding of these nematodes’ cellular and molecular biology – which is presently lacking for most worms. Here, we created the first single cell atlas for an early developmental stage of Haemonchus contortus – a highly pathogenic, C. elegans-related parasitic nematode. We obtained and curated RNA sequence (snRNA-seq) data from single nuclei from embryonating eggs of H. contortus (150,000 droplets), and selected high-quality transcriptomic data for > 14,000 single nuclei for analysis, and identified 19 distinct clusters of cells. Guided by comparative analyses with C. elegans, we were able to reproducibly assign seven cell clusters to body wall muscle, hypodermis, neuronal, intestinal or seam cells, and identified eight genes that were transcribed in all cell clusters/types, three of which were inferred to be essential in H. contortus. Two of these genes (i.e. Hc-eef-1A and Hc-eef1G), coding for eukaryotic elongation factors (called Hc-eEF1A and Hc-eEF1G), were also demonstrated to be transcribed and expressed in all key developmental stages of H. contortus. Together with these findings, sequence- and structure-based comparative analyses indicated the potential of Hc-eEF1A and/or Hc-eEF1G as intervention targets within the protein biosynthesis machinery of H. contortus. Future work will focus on single cell studies of all key developmental stages and tissues of H. contortus, and on evaluating the suitability of the two elongation factor proteins as drug targets in H. contortus and related nematodes, with a view to finding new nematocidal drug candidates.

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.

Within the context of our anthelmintic discovery program, we recently identified and evaluated a quinoline derivative, called ABX464 or obefazimod, as a nematocidal candidate; synthesised a series of analogues which were assessed for activity against the free-living nematode Caenorhabditis elegans; and predicted compound-target relationships by thermal proteome profiling (TPP) and in silico docking. Here, we logically extended this work and critically evaluated the anthelmintic activity of ABX464 analogues on Haemonchus contortus (barber's pole worm) – a highly pathogenic nematode of ruminant livestock. First, we tested a series of 44 analogues on H. contortus (larvae and adults) to investigate the nematocidal pharmacophore of ABX464, and identified one compound with greater potency than the parent compound and showed moderate activity against a select number of other parasitic nematodes (including Ancylostoma, Heligmosomoides and Strongyloides species). Using TPP and in silico modelling studies, we predicted protein HCON_00074590 (a predicted aldo-keto reductase) as a target candidate for ABX464 in H. contortus. Future work aims to optimise this compound as a nematocidal candidate and investigate its pharmacokinetic properties. Overall, this study presents a first step toward the development of a new nematocide.

Haemonchus contortus (the barber’s pole worm)—a highly pathogenic gastric nematode of ruminants—causes significant economic losses in the livestock industry worldwide. H. contortus has become a valuable model organism for both fundamental and applied research (e.g., drug and vaccine discovery) because of the availability of well-defined laboratory strains (e.g., MHco3(ISE).N1 in the UK and Haecon-5 in Australia) and genomic, transcriptomic and proteomic data sets. Many recent investigations have relied heavily on the use of the chromosome-contiguous genome of MHco3(ISE).N1 in the absence of a genome for Haecon-5. However, there has been no genetic comparison of these and other strains to date. Here, we assembled and characterised the mitochondrial genome (14.1 kb) of Haecon-5 and compared it with that of MHco3(ISE).N1 and two other strains (i.e., McMaster and NZ_Hco_NP) from Australasia. We detected 276 synonymous and 25 non-synonymous single nucleotide polymorphisms (SNPs) within Haecon-5. Between the Haecon-5 and MHco3(ISE).N1 strains, we recorded 345 SNPs, 31 of which were non-synonymous and linked to fixed amino acid differences in seven protein-coding genes (nad5, nad6, nad1, atp6, nad2, cytb and nad4) between these strains. Pronounced variation (344 and 435 SNPs) was seen between Haecon-5 and each of the other two strains from Australasia. The question remains as to what impact these mitogenomic mutations might have on the biology and physiology of H. contortus, which warrants exploration. The high degree of mitogenomic variability recorded here among these strains suggests that further work should be undertaken to assess the nature and extent of the nuclear genomic variation within H. contortus.

Nematodes of the genus Trichinella are important pathogens of humans and animals. This study aimed to enhance the genomic and transcriptomic resources for T. pseudospiralis (non-encapsulated phenotype) and T. spiralis (encapsulated phenotype) and to explore transcriptional profiles. First, we improved the assemblies of the genomes of T. pseudospiralis (code ISS13) and T. spiralis (code ISS534), achieving genome sizes of 56.6 Mb (320 scaffolds, and an N50 of 1.02 Mb) and 63.5 Mb (568 scaffolds, and an N50 value of 0.44 Mb), respectively. Then, for each species, we produced RNA sequence data for three key developmental stages (first-stage muscle larvae [L1s], adults, and newborn larvae [NBLs]; three replicates for each stage), analysed differential transcription between stages, and explored enriched pathways and processes between species. Stage-specific upregulation was linked to cellular processes, metabolism, and host–parasite interactions, and pathway enrichment analysis showed distinctive biological processes and cellular localisations between species. Indeed, the secreted molecules calmodulin, calreticulin, and calsyntenin—with possible roles in modulating host immune responses and facilitating parasite survival—were unique to T. pseudospiralis and not detected in T. spiralis. These insights into the molecular mechanisms of Trichinella–host interactions might offer possible avenues for developing new interventions against trichinellosis.

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

Over the years, comprehensive explorations of the model organisms Caenorhabditis elegans (elegant worm) and Drosophila melanogaster (vinegar fly) have contributed substantially to our understanding of complex biological processes and pathways in multicellular organisms generally. Extensive functional genomic–phenomic, genomic, transcriptomic, and proteomic data sets have enabled the discovery and characterisation of genes that are crucial for life, called ‘essential genes’. Recently, we investigated the feasibility of inferring essential genes from such data sets using advanced bioinformatics and showed that a machine learning (ML)-based workflow could be used to extract or engineer features from DNA, RNA, protein, and/or cellular data/information to underpin the reliable prediction of essential genes both within and between C. elegans and D. melanogaster. As these are two distantly related species within the Ecdysozoa, we proposed that this ML approach would be particularly well suited for species that are within the same phylum or evolutionary clade. In the present study, we cross-predicted essential genes within the phylum Nematoda (evolutionary clade V)—between C. elegans and the pathogenic parasitic nematode H. contortus—and then ranked and prioritised H. contortus proteins encoded by these genes as intervention (e.g., drug) target candidates. Using strong, validated predictors, we inferred essential genes of H. contortus that are involved predominantly in crucial biological processes/pathways including ribosome biogenesis, translation, RNA binding/processing, and signalling and which are highly transcribed in the germline, somatic gonad precursors, sex myoblasts, vulva cell precursors, various nerve cells, glia, or hypodermis. The findings indicate that this in silico workflow provides a promising avenue to identify and prioritise panels/groups of drug target candidates in parasitic nematodes for experimental validation in vitro and/or in vivo.

Major advances in genomic and associated technologies have demanded reliable bioinformatic tools and workflows for the annotation of genes and their products via comparative analyses using well-curated reference data sets, accessible in public repositories. However, the accurate in silico annotation of molecules (proteins) encoded in organisms (e.g., multicellular parasites) which are evolutionarily distant from those for which these extensive reference data sets are available, including invertebrate model organisms (e.g., Caenorhabditis elegans – free-living nematode, and Drosophila melanogaster – the vinegar fly) and vertebrate species (e.g., Homo sapiens and Mus musculus), remains a major challenge. Here, we constructed an informatic workflow for the enhanced annotation of biologically-important, excretory/secretory (ES) proteins (“secretome”) encoded in the genome of a parasitic roundworm, called Haemonchus contortus (commonly known as the barber’s pole worm). We critically evaluated the performance of five distinct methods, refined some of them, and then combined the use of all five methods to comprehensively annotate ES proteins, according to gene ontology, biological pathways and/or metabolic (enzymatic) processes. Then, using optimised parameter settings, we applied this workflow to comprehensively annotate 2591 of all 3353 proteins (77.3%) in the secretome of H. contortus. This result is a substantial improvement (10–25%) over previous annotations using individual, “off-the-shelf” algorithms and default settings, indicating the ready applicability of the present, refined workflow to gene/protein sequence data sets from a wide range of organisms in the Tree-of-Life.

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

A high-quality genome for the parasitic tapeworm, Echinococcus granulosus, provides further insight into the EG95 vaccine target for cystic echinococcosis.

The snail Bulinus truncatus is an intermediate host of the carcinogenic human blood fluke Schistosoma haematobium. Here the authors report the genome of Bu. truncatus, explore protein groups inferred to play a role in its interaction with the schistosome parasite, and identify expansions in gene families linked to immune response regulation.

Biodiversity within the animal kingdom is associated with extensive molecular diversity. The expansion of genomic, transcriptomic and proteomic data sets for invertebrate groups and species with unique biological traits necessitates reliable in silico tools for the accurate identification and annotation of molecules and molecular groups. However, conventional tools are inadequate for lesser-known organismal groups, such as eukaryotic pathogens (parasites), so that improved approaches are urgently needed. Here, we established a combined sequence- and structure-based workflow system to harness well-curated publicly available data sets and resources to identify, classify and annotate proteases and protease inhibitors of a highly pathogenic parasitic roundworm (nematode) of global relevance, called Haemonchus contortus (barber’s pole worm). This workflow performed markedly better than conventional, sequence-based classification and annotation alone and allowed the first genome-wide characterisation of protease and protease inhibitor genes and gene products in this worm. In total, we identified 790 genes encoding 860 proteases and protease inhibitors representing 83 gene families. The proteins inferred included 280 metallo-, 145 cysteine, 142 serine, 121 aspartic and 81 “mixed” proteases as well as 91 protease inhibitors, all of which had marked physicochemical diversity and inferred involvements in >400 biological processes or pathways. A detailed investigation revealed a remarkable expansion of some protease or inhibitor gene families, which are likely linked to parasitism (e.g., host–parasite interactions, immunomodulation and blood-feeding) and exhibit stage- or sex-specific transcription profiles. This investigation provides a solid foundation for detailed explorations of the structures and functions of proteases and protease inhibitors of H. contortus and related nematodes, and it could assist in the discovery of new drug or vaccine targets against infections or diseases.

Abstract The redlegged earth mite, Halotydeus destructor (Tucker, 1925: Trombidiformes, Eupodoidea, Penthaleidae), is an invasive mite species. In Australia, this mite has become a pest of winter pastures and grain crops. We report the complete mitogenome for H. destructor, the first to represent the family Penthaleidae, superfamily Eupodoidea. The mitogenome of H. destructor is 14,691 bp in size, and has a GC content of 27.87%, 13 protein‐coding genes, two rRNA genes, and 22 tRNA genes. We explored evolutionary relationships of H. destructor with other members of the Trombidiformes using phylogenetic analyses of nucleotide sequences and the order of protein‐coding and rRNA genes. We found strong, consistent support for the superfamily Tydeoidea being the sister taxon to the superfamily Eupodoidea based on nucleotide sequences and gene arrangements. Moreover, the gene arrangements of Eupodoidea and Tydeoidea are not only identical to each other but also identical to that of the hypothesized arthropod ancestor, showing a high level of conservatism in the mitogenomic structure of these mite superfamilies. Our study illustrates the utility of gene arrangements for providing complementary information to nucleotide sequences with respect to inferring the evolutionary relationships of species within the order Trombidiformes. The mitogenome of H. destructor provides a valuable resource for further population genetic studies of this important agricultural pest. Given the co‐occurrence of closely related, morphologically similar Penthaleidae mites with H. destructor in the field, a complete mitogenome provides new opportunities to develop metabarcoding tools to study mite diversity in agro‐ecosystems. Moreover, the H. destructor mitogenome fills an important taxonomic gap that will facilitate further study of trombidiform mite evolution.

Presently, there is a dearth of proteomic data for parasitic mites and their relationship with the host animals. Here, using a high throughput LC-MS/MS-based approach, we undertook the first comprehensive, large-scale proteomic investigation of egg and adult female stages of the scabies mite, Sarcoptes scabiei–one of the most important parasitic mites of humans and other animals worldwide. In total, 1,761 S. scabiei proteins were identified and quantified with high confidence. Bioinformatic analyses revealed differentially expressed proteins to be involved predominantly in biological pathways or processes including genetic information processing, energy (oxidative phosphorylation), nucleotide, amino acid, carbohydrate and/or lipid metabolism, and some adaptive processes. Selected, constitutively and highly expressed proteins, such as peptidases, scabies mite inactivated protease paralogues (SMIPPs) and muscle proteins (myosin and troponin), are proposed to be involved in key biological processes within S. scabiei, host-parasite interactions and/or the pathogenesis of scabies. These proteomic data will enable future molecular, biochemical and physiological investigations of early developmental stages of S. scabiei and the discovery of novel interventions, targeting the egg stage, given its non-susceptibility to acaricides currently approved for the treatment of scabies in humans. Author summary Scabies is a neglected tropical disease caused by the parasitic mite Sarcoptes scabiei. The treatment and control of scabies are challenging, as there is no vaccine and the two mostly used broad-spectrum acaricides (i.e. ivermectin and permethrin) do not kill the key developmental stage (egg) of the mite that enables the re-establishment of infection. The availability of a well-assembled genome for S. scabiei now provides a foundation to explore the molecular biology, biochemistry and physiology of this mite. Here, we characterised the first somatic proteome of key developmental stages of S. scabiei using high throughput LC-MS/MS. Bioinformatic analyses of proteomic data indicate that proteins expressed differentially between egg and female adult stages are mainly involved in biological pathways or processes, such as genetic information processing, energy (oxidative phosphorylation), nucleotide, amino acid, carbohydrate and/or lipid metabolism in the mite. These proteomic data should underpin further investigations of early developmental stages of S. scabiei with a focus on identifying novel intervention targets for scabies.

Parasitic roundworms (nematodes) cause destructive diseases, and immense suffering in humans and other animals around the world. The control of these parasites relies heavily on anthelmintic therapy, but treatment failures and resistance to these drugs are widespread. As efforts to develop vaccines against parasitic nematodes have been largely unsuccessful, there is an increased focus on discovering new anthelmintic entities to combat drug resistant worms. Here, we employed thermal proteome profiling (TPP) to explore hit pharmacology and to support optimisation of a hit compound (UMW-868), identified in a high-throughput whole-worm, phenotypic screen. Using advanced structural prediction and docking tools, we inferred an entirely novel, parasite-specific target (HCO_011565) of this anthelmintic small molecule in the highly pathogenic, blood-feeding barber’s pole worm, and in other socioeconomically important parasitic nematodes. The “hit-to-target” workflow constructed here provides a unique prospect of accelerating the simultaneous discovery of novel anthelmintics and associated parasite-specific targets.

Cytotoxic lymphocytes are essential for anti-tumor immunity, and for effective responses to cancer immunotherapy. Natural killer cell granule protein 7 (NKG7) is expressed at high levels in cytotoxic lymphocytes infiltrating tumors from patients treated with immunotherapy, but until recently, the role of this protein in cytotoxic lymphocyte function was largely unknown. Unexpectedly, we found that highly CD8+ T cell-immunogenic murine colon carcinoma (MC38-OVA) tumors grew at an equal rate in Nkg7+/+ and Nkg7-/- littermate mice, suggesting NKG7 may not be necessary for effective CD8+ T cell anti-tumor activity. Mechanistically, we found that deletion of NKG7 reduces the ability of CD8+ T cells to degranulate and kill target cells in vitro. However, as a result of inefficient cytotoxic activity, NKG7 deficient T cells form a prolonged immune synapse with tumor cells, resulting in increased secretion of inflammatory cytokines, including tumor necrosis factor alpha (TNF). By deleting the TNF receptor, TNFR1, from MC38-OVA tumors, we demonstrate that this hyper-secretion of TNF compensates for reduced synapse-mediated cytotoxic activity against MC38-OVA tumors in vivo, via increased TNF-mediated tumor cell death. Taken together, our results demonstrate that NKG7 enhances CD8+ T cell immune synapse efficiency, which may serve as a mechanism to accelerate direct cytotoxicity and limit potentially harmful inflammatory responses.

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.

Urogenital schistosomiasis is caused by the blood fluke Schistosoma haematobium and is one of the most neglected tropical diseases worldwide, afflicting > 100 million people. It is characterised by granulomata, fibrosis and calcification in urogenital tissues, and can lead to increased susceptibility to HIV/AIDS and squamous cell carcinoma of the bladder. To complement available treatment programs and break the transmission of disease, sound knowledge and understanding of the biology and ecology of S. haematobium is required. Hybridisation/introgression events and molecular variation among members of the S. haematobium-group might effect important biological and/or disease traits as well as the morbidity of disease and the effectiveness of control programs including mass drug administration. Here we report the first chromosome-contiguous genome for a well-defined laboratory line of this blood fluke. An exploration of this genome using transcriptomic data for all key developmental stages allowed us to refine gene models (including non-coding elements) and annotations, discover ‘new’ genes and transcription profiles for these stages, likely linked to development and/or pathogenesis. Molecular variation within S. haematobium among some geographical locations in Africa revealed unique genomic ‘signatures’ that matched species other than S. haematobium, indicating the occurrence of introgression events. The present reference genome (designated Shae.V3) and the findings from this study solidly underpin future functional genomic and molecular investigations of S. haematobium and accelerate systematic, large-scale population genomics investigations, with a focus on improved and sustained control of urogenital schistosomiasis. Author summary More than 100 million people are infected with the carcinogenic blood fluke Schistosoma haematobium, the aetiological agent of urogenital schistosomiasis—a neglected tropical disease (NTD). In spite of its major significance, little is known about this fluke, its interactions with the human and snail intermediate hosts and the pathogenesis of the urogenital form of schistosomiasis at the molecular and biochemical levels. To enable research in these areas, we report the first chromosome-level genome and markedly enhanced gene models for S. haematobium. Comparative genomic analyses also reveal evidence of past introgression events between or among closely related schistosome species. This present reference genome for S. haematobium and the findings from this study should underpin future functional genomic and molecular investigations of S. haematobium and accelerate systematic, large-scale population genomics investigations, with a focus on improved control of urogenital schistosomiasis.

Defining genes that are essential for life has major implications for understanding critical biological processes and mechanisms. Although essential genes have been identified and characterised experimentally using functional genomic tools, it is challenging to predict with confidence such genes from molecular and phenomic data sets using computational methods. Using extensive data sets available for the model organism Caenorhabditis elegans, we constructed here a machine-learning (ML)-based workflow for the prediction of essential genes on a genome-wide scale. We identified strong predictors for such genes and showed that trained ML models consistently achieve highly-accurate classifications. Complementary analyses revealed an association between essential genes and chromosomal location. Our findings reveal that essential genes in C. elegans tend to be located in or near the centre of autosomal chromosomes; are positively correlated with low single nucleotide polymorphim (SNP) densities and epigenetic markers in promoter regions; are involved in protein and nucleotide processing; are transcribed in most cells; are enriched in reproductive tissues or are targets for small RNAs bound to the argonaut CSR-1. Based on these results, we hypothesise an interplay between epigenetic markers and small RNA pathways in the germline, with transcription-based memory; this hypothesis warrants testing. From a technical perspective, further work is needed to evaluate whether the present ML-based approach will be applicable to other metazoans (including Drosophila melanogaster) for which comprehensive data sets (i.e. genomic, transcriptomic, proteomic, variomic, epigenetic and phenomic) are available.

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.

Here, we explored transcriptomic differences among early egg (Ee), late egg (Le) and adult female (Af) stages of the scabies mite, Sarcoptes scabiei, using an integrative bioinformatic approach. We recorded a high, negative correlation between miRNAs and genes with decreased mRNA transcription between the developmental stages, indicating substantial post-transcriptional repression; we also showed a positive correlation between miRNAs and genes with increased mRNA transcription, suggesting indirect post-transcriptional regulation. The alterations in mRNA transcription between the egg and adult female stages of S. scabiei were inferred to be linked to metabolism (including carbohydrate and lipid degradation, amino acid and energy metabolism), environmental information processing (e.g., signal transduction and signalling molecules), genetic information processing (e.g., transcription and translation) and/or organismal systems. Taken together, these results provide insight into the transcription of this socioeconomically important parasitic mite, with a particular focus on the egg stage. This work encourages further, detailed laboratory studies of miRNA regulation across all developmental stages of S. scabiei and might assist in discovering new intervention targets in the egg stage of S. scabiei.

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.

The availability of whole-genome sequences and associated multi-omics data sets, combined with advances in gene knockout and knockdown methods, has enabled large-scale annotation and exploration of gene and protein functions in eukaryotes. Knowing which genes are essential for the survival of eukaryotic organisms is paramount for an understanding of the basic mechanisms of life, and could assist in identifying intervention targets in eukaryotic pathogens and cancer. Here, we studied essential gene orthologs among selected species of eukaryotes, and then employed a systematic machine-learning approach, using protein sequence-derived features and selection procedures, to investigate essential gene predictions within and among species. We showed that the numbers of essential gene orthologs comprise small fractions when compared with the total number of orthologs among the eukaryotic species studied. In addition, we demonstrated that machine-learning models trained with subsets of essentiality-related data performed better than random guessing of gene essentiality for a particular species. Consistent with our gene ortholog analysis, the predictions of essential genes among multiple (including distantly-related) species is possible, yet challenging, suggesting that most essential genes are unique to a species. The present work provides a foundation for the expansion of genome-wide essentiality investigations in eukaryotes using machine learning approaches. Keywords: Machine-learning, Essential genes, Essentiality prediction, Eukaryotes

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.

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.

Many freshwater snails of the genus Bulinus act as intermediate hosts in the life-cycles of schistosomes in Africa and adjacent regions. Currently, 37 species of Bulinus representing four groups are recognised. The mitochondrial cytochrome c oxidase subunit 1 (cox1) gene has shown utility for identifying and differentiating Bulinus species and groups, but taxonomic relationships based on genetic data are not entirely consistent with those inferred using morphological and biological features. To underpin future systematic studies of members of the genus, we characterised here the mitochondrial genome of Bulinus truncatus (from a defined laboratory strain) using a combined second- and third-generation sequencing and informatics approach, enabling taxonomic comparisons with other planorbid snails for which mitochondrial (mt) genomes were available. Analyses showed consistency in gene order and length among mitochondrial genomes of representative planorbid snails, with the lowest and highest nucleotide diversities being in the cytochrome c oxidase and nicotinamide dehydrogenase subunit genes, respectively. This first mt genome for a representative of the genus Bulinus should provide a useful resource for future investigations of the systematics, population genetics, epidemiology and/or ecology of Bulinus and related snails. The sequencing and informatic workflow employed here should find broad applicability to a range of other snail intermediate hosts of parasitic trematodes.

Despite advances in high-throughput sequencing and bioinformatics, molecular investigations of snail intermediate hosts that transmit parasitic trematodes are scant. Here, we report the first transcriptome for Bulinus truncatus – a key intermediate host of Schistosoma haematobium – a blood fluke that causes urogenital schistosomiasis in humans. We assembled this transcriptome from short- and long-read RNA-sequence data. From this transcriptome, we predicted 12,998 proteins, 58% of which had orthologs in Biomphalaria glabrata – an intermediate host of Schistosoma mansoni – a blood fluke that causes hepato-intestinal schistosomiasis. We predicted that select protein groups are involved in signal transduction, cell growth and death, the immune system, environmental adaptation and/or the excretory/secretory system, suggesting roles in immune responses, pathogen defence and/or parasite-host interactions. The transcriptome of Bu. truncatus provides a useful resource to underpin future molecular investigations of this and related snail species, and its interactions with pathogens including S. haematobium. The present resource should enable comparative investigations of other molluscan hosts of socioeconomically important parasites in the future.

The prevalent porcine helminth, Ascaris suum, compromises pig health and reduces farm productivity worldwide. The closely related human parasite, A. lumbricoides, infects more than 800 million people representing a disease burden of 1.31 million disability-adjusted life years. The infections are often chronic in nature, and the parasites have a profound ability to modulate their hosts’ immune responses. This study provides the first in-depth characterisation of extracellular vesicles (EVs) from different developmental stages and body parts of A. suum and proposes the role of these vesicles in the host–parasite interplay. The release of EVs from the third- (L3) and fourth-stage (L4) larvae and adults was demonstrated by transmission electron microscopy (TEM), and sequencing of EV-derived RNA identified a number of microRNAs (miRNAs) and transcripts of potential host immune targets, such as IL-13, IL-25 and IL-33, were identified. Furthermore, proteomics of EVs identified several proteins with immunomodulatory properties and other proteins previously shown to be associated with parasite EVs. Taken together, these results suggest that A. suum EVs and their cargo may play a role in host–parasite interactions. This knowledge may pave the way to novel strategies for helminth infection control and knowledge of their immune modulatory potential.

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

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.

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.

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.

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.

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.

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 Trichinellaspecies and genotypes, and reconstruct their phylogeny and biogeography.

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.

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.