Search

The Chinese liver fluke, Clonorchis sinensis, causes the disease clonorchiasis, affecting ~35 million people in regions of China, Vietnam, Korea and the Russian Far East. Chronic clonorchiasis causes cholangitis and can induce a malignant cancer, called cholangiocarcinoma, in the biliary system. Control in endemic regions is challenging, and often relies largely on chemotherapy with one anthelmintic, called praziquantel. Routine treatment carries a significant risk of inducing resistance to this anthelmintic in the fluke, such that the discovery of new interventions is considered important. It is hoped that the use of molecular technologies will assist this endeavour by enabling the identification of drug or vaccine targets involved in crucial biological processes and/or pathways in the parasite. Although draft genomes of C. sinensis have been published, their assemblies are fragmented. In the present study, we tackle this genome fragmentation issue by utilising, in an integrated way, advanced (second- and third-generation) DNA sequencing and informatic approaches to build a high-quality reference genome for C. sinensis, with chromosome-level contiguity and curated gene models. This substantially-enhanced genome provides a resource that could accelerate fundamental and applied molecular investigations of C. sinensis, clonorchiasis and/or cholangiocarcinoma, and assist in the discovery of new interventions against what is a highly significant, but neglected disease-complex.

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.

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.

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., Graphical Abstract Unlabelled Image

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

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

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

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

We investigated genome folding across the eukaryotic tree of life. We find two types of three-dimensional(3D) genome architectures at the chromosome scale. Each type appears and disappears repeatedlyduring eukaryotic evolution. The type of genome architecture that an organism exhibits correlates with theabsence of condensin II subunits. Moreover, condensin II depletion converts the architecture of thehuman genome to a state resembling that seen in organisms such as fungi or mosquitoes. In this state,centromeres cluster together at nucleoli, and heterochromatin domains merge. We propose a physicalmodel in which lengthwise compaction of chromosomes by condensin II during mitosis determineschromosome-scale genome architecture, with effects that are retained during the subsequent interphase.This mechanism likely has been conserved since the last common ancestor of all eukaryotes., C.H. is supported by the Boehringer Ingelheim Fonds; C.H., Á.S.C., and B.D.R. are supported by an ERC CoG (772471, “CohesinLooping”); A.M.O.E. and B.D.R. are supported by the Dutch Research Council (NWO-Echo); and J.A.R. and R.H.M. are supported by the Dutch Cancer Society (KWF). T.v.S. and B.v.S. are supported by NIH Common Fund “4D Nucleome” Program grant U54DK107965. H.T. and E.d.W. are supported by an ERC StG (637597, “HAP-PHEN”). J.A.R., T.v.S., H.T., R.H.M., B.v.S., and E.d.W. are part of the Oncode Institute, which is partly financed by the Dutch Cancer Society. Work at the Center for Theoretical Biological Physics is sponsored by the NSF (grants PHY-2019745 and CHE-1614101) and by the Welch Foundation (grant C-1792). V.G.C. is funded by FAPESP (São Paulo State Research Foundation and Higher Education Personnel) grants 2016/13998-8 and 2017/09662-7. J.N.O. is a CPRIT Scholar in Cancer Research. E.L.A. was supported by an NSF Physics Frontiers Center Award (PHY-2019745), the Welch Foundation (Q-1866), a USDA Agriculture and Food Research Initiative grant (2017-05741), the Behavioral Plasticity Research Institute (NSF DBI-2021795), and an NIH Encyclopedia of DNA Elements Mapping Center Award (UM1HG009375). Hi-C data for the 24 species were created by the DNA Zoo Consortium (www.dnazoo.org). DNA Zoo is supported by Illumina, Inc.; IBM; and the Pawsey Supercomputing Center. P.K. is supported by the University of Western Australia. L.L.M. was supported by NIH (1R01NS114491) and NSF awards (1557923, 1548121, and 1645219) and the Human Frontiers Science Program (RGP0060/2017). The draft A. californica project was supported by NHGRI. J.L.G.-S. received funding from the ERC (grant agreement no. 740041), the Spanish Ministerio de Economía y Competitividad (grant no. BFU2016-74961-P), and the institutional grant Unidad de Excelencia María de Maeztu (MDM-2016-0687). R.D.K. is supported by NIH grant RO1DK121366. V.H. is supported by NIH grant NIH1P41HD071837. K.M. is supported by a MEXT grant (20H05936). M.C.W. is supported by the NIH grants R01AG045183, R01AT009050, R01AG062257, and DP1DK113644 and by the Welch Foundation. E.F. was supported by NHGRI

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

Early studies of the free-living nematode C. elegans informed us how BCL-2-regulated apoptosis in humans is regulated. However, subsequent studies showed C. elegans apoptosis has several unique features compared with human apoptosis. To date, there has been no detailed analysis of apoptosis regulators in nematodes other than C. elegans. Here, we discovered BCL-2 orthologues in 89 free-living and parasitic nematode taxa representing four evolutionary clades (I, III, IV and V). Unlike in C. elegans, 15 species possess multiple (two to five) BCL-2-like proteins, and some do not have any recognisable BCL-2 sequences. Functional studies provided no evidence that BAX/BAK proteins have evolved in nematodes, and structural studies of a BCL-2 protein from the basal clade I revealed it lacks a functionally important feature of the C. elegans orthologue. Clade I CED-4/APAF-1 proteins also possess WD40-repeat sequences associated with apoptosome assembly, not present in C. elegans, or other nematode taxa studied., Young et al. report on the diversity in regulation of nematode apoptosis. Results indicate that the apoptotic pathways of C. elegans are not fully representative of nematodes generally.

Characterizing genes that are critical for the survival of an organism (i.e. essential) is important to gain a deep understanding of the fundamental cellular and molecular mechanisms that sustain life. Functional genomic investigations of the vinegar fly, Drosophila melanogaster, have unravelled the functions of numerous genes of this model species, but results from phenomic experiments can sometimes be ambiguous. Moreover, the features underlying gene essentiality are poorly understood, posing challenges for computational prediction. Here, we harnessed comprehensive genomic-phenomic datasets publicly available for D. melanogaster and a machine-learning-based workflow to predict essential genes of this fly. We discovered strong predictors of such genes, paving the way for computational predictions of essentiality in less-studied arthropod pests and vectors of infectious diseases.

Here, we present a draft genome of the tapeworm Dipylidium caninum (family Dipylidiidae) and compare it with other cestode genomes. This draft genome of D. caninum is 110 Mb in size, has a repeat content of ~13.4% and is predicted to encode ~10,000 protein-coding genes. We inferred excretory/secretory molecules (representing the secretome), other key groups of proteins (including peptidases, kinases, phosphatases, GTPases, receptors, transporters and ion-channels) and predicted potential intervention targets for future evaluation. Using 144 shared single-copy orthologous sequences, we investigated the genetic relationships of cestodes for which nuclear genomes are available. This study provides first insights into the molecular biology of D. caninum and a new resource for comparative genomic and genetic explorations of this and other flatworms.

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

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

Mycoplasma gallisepticum is the primary etiological agent of chronic respiratory disease in chickens. Live attenuated vaccines are most commonly used in the field to control the disease, but current vaccines have some limitations. Vaxsafe MG (strain ts-304) is a new vaccine candidate that is efficacious at a lower dose than the current commercial vaccine strain ts-11, from which it is derived. In this study, the transcriptional profiles of the trachea of unvaccinated chickens and chickens vaccinated with strain ts-304 were compared 2 weeks after challenge with M. gallisepticum strain Ap3AS during the chronic stage of infection. After challenge, genes, gene ontologies, pathways, and protein classes involved in inflammation, cytokine production and signaling, and cell proliferation were upregulated, while those involved in formation and motor movement of cilia, formation of intercellular junctional complexes, and formation of the cytoskeleton were downregulated in the unvaccinated birds compared to the vaccinated birds, reflecting immune dysregulation and the pathological changes induced in the trachea by infection with M. gallisepticum. Vaccination appears to protect the structural and functional integrity of the tracheal mucosa 2 weeks after infection with M. gallisepticum.

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

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

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

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

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

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

BackgroundSchistosoma haematobium causes urogenital schistosomiasis, a neglected tropical disease affecting >100 million people worldwide. Chronic infection with this parasitic trematode can lead to urogenital conditions including female genital schistosomiasis and bladder cancer. At the molecular level, little is known about this blood fluke and the pathogenesis of the disease that it causes. To support molecular studies of this carcinogenic worm, we reported a draft genome for S. haematobium in 2012. Although a useful resource, its utility has been somewhat limited by its fragmentation.FindingsHere, we systematically enhanced the draft genome of S. haematobium using a single-molecule and long-range DNA-sequencing approach. We achieved a major improvement in the accuracy and contiguity of the genome assembly, making it superior or comparable to assemblies for other schistosome species. We transferred curated gene models to this assembly and, using enhanced gene annotation pipelines, inferred a gene set with as many or more complete gene models as those of other well-studied schistosomes. Using conserved, single-copy orthologs, we assessed the phylogenetic position of S. haematobium in relation to other parasitic flatworms for which draft genomes were available.ConclusionsWe report a substantially enhanced genomic resource that represents a solid foundation for molecular research on S. haematobium and is poised to better underpin population and functional genomic investigations and to accelerate the search for new disease interventions.

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

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

Colletotrichum tanaceti is an emerging foliar fungal pathogen of commercially grown pyrethrum (Tanacetum cinerariifolium). Despite being reported consistently from field surveys in Australia, the molecular basis of pathogenicity of C. tanaceti on pyrethrum is unknown. Herein, the genome of C. tanaceti (isolate BRIP57314) was assembled de novo and annotated using transcriptomic evidence. The inferred putative pathogenicity gene suite of C. tanaceti comprised a large array of genes encoding secreted effectors, proteases, CAZymes and secondary metabolites. Comparative analysis of its putative pathogenicity gene profiles with those of closely related species suggested that C. tanaceti likely has additional hosts to pyrethrum. The genome of C. tanaceti had a high repeat content and repetitive elements were located significantly closer to genes inferred to influence pathogenicity than other genes. These repeats are likely to have accelerated mutational and transposition rates in the genome, resulting in a rapid evolution of certain CAZyme families in this species. The C. tanaceti genome showed strong signals of Repeat Induced Point (RIP) mutation which likely caused its bipartite nature consisting of distinct gene-sparse, repeat and A-T rich regions. Pathogenicity genes within these RIP affected regions were likely to have a higher evolutionary rate than the rest of the genome. This "two-speed" genome phenomenon in certain Colletotrichum spp. was hypothesized to have caused the clustering of species based on the pathogenicity genes, to deviate from taxonomic relationships. The large repertoire of pathogenicity factors that potentially evolve rapidly due to the plasticity of the genome, indicated that C. tanaceti has a high evolutionary potential. Therefore, C. tanaceti poses a high-risk to the pyrethrum industry. Knowledge of the evolution and diversity of the putative pathogenicity genes will facilitate future research in disease management of C. tanaceti and other Colletotrichum spp.

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

Background Here, we created an automated pipeline for the de novoassembly of genomes from Pacific Biosciences long-read and Illumina short-read data using common workflow language (CWL). To evaluate the performance of this pipeline, we assembled the nuclear genomes of the eukaryotes Caenorhabditis elegans (∼100 Mb), Drosophila melanogaster (∼138 Mb), and Plasmodium falciparum (∼23 Mb) directly from publicly accessible nucleotide sequence datasets and assessed the quality of the assemblies against curated reference genomes. Findings We showed a dependency of the accuracy of assembly on sequencing technology and GC content and repeatedly achieved assemblies that meet the high standards set by the National Human Genome Research Institute, being applicable to gene prediction and subsequent genomic analyses. Conclusions This CWL pipeline overcomes current challenges of achieving repeatability and reproducibility of assembly results and offers a platform for the re-use of the workflow and the integration of diverse datasets. This workflow is publicly available via GitHub (https://github.com/vetscience/Assemblosis) and is currently applicable to the assembly of haploid and diploid genomes of eukaryotes.

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

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

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

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.

16 páginas, 1 tabla, 3 figuras., Mesenteric infection by the parasitic blood fluke Schistosoma bovis is a common veterinary problem in Africa and the Middle East and occasionally in the Mediterranean Region. The species also has the ability to form interspecific hybrids with the human parasite S. haematobium with natural hybridisation observed in West Africa, presenting possible zoonotic transmission. Additionally, this exchange of alleles between species may dramatically influence disease dynamics and parasite evolution. We have generated a 374 Mb assembly of the S. bovis genome using Illumina and PacBio-based technologies. Despite infecting different hosts and organs, the genome sequences of S. bovis and S. haematobium appeared strikingly similar with 97% sequence identity. The two species share 98% of protein-coding genes, with an average sequence identity of 97.3% at the amino acid level. Genome comparison identified large continuous parts of the genome (up to several 100 kb) showing almost 100% sequence identity between S. bovis and S. haematobium. It is unlikely that this is a result of genome conservation and provides further evidence of natural interspecific hybridization between S. bovis and S. haematobium. Our results suggest that foreign DNA obtained by interspecific hybridization was maintained in the population through multiple meiosis cycles and that hybrids were sexually reproductive, producing viable offspring. The S. bovis genome assembly forms a highly valuable resource for studying schistosome evolution and exploring genetic regions that are associated with species-specific phenotypic traits., This work was supported by grants from the QIMR Berghofer Medical Research Institute (Chenhall Estate) and the Australian Infectious Diseases Research Centre. NDY was supported by a National Health and Medical Research Council (NHMRC) CDF1, PKK by an NHMRC ECRF and DPM by an NHMRC SPRF (APP1102926). The work has been supported by NHMRC project grant APP1102322.

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

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

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

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

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

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

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

Infections with several flatworm parasites represent group 1 biological carcinogens, i.e. definite causes of cancer. Infection with the food-borne liver fluke Opisthorchis viverrini causes cholangiocarcinoma (CCA). Whereas the causative agent for most cancers, including CCA in the West, remains obscure, the principal risk factor for CCA in Thailand is opisthorchiasis. We exploited this established link to explore the role of the secreted parasite growth factor termed liver fluke granulin (Ov-GRN-1) in pre-malignant lesions of the biliary tract. We targeted the Ov-grn-1 gene for programmed knockout and investigated gene-edited parasites in vitro and in experimentally infected hamsters. Both adult and juvenile stages of the liver fluke were transfected with a plasmid encoding a guide RNA sequence specific for exon 1 of Ov-grn-1 and the Cas9 nuclease. Deep sequencing of amplicon libraries from genomic DNA from gene-edited parasites exhibited programmed, Cas9-catalyzed mutations within the Ov-grn-1 locus, and tandem analyses by RT-PCR and western blot revealed rapid depletion of Ov-grn-1 transcripts and protein. Newly excysted juvenile flukes that had undergone editing of Ov-grn-1 colonized the biliary tract, grew and developed over a period of 60 days, were active and motile, and induced a clinically relevant pathophysiological tissue phenotype of attenuated biliary hyperplasia and fibrosis in comparison to infection with wild type flukes. This is the first report of gene knock-out using CRISPR/Cas9 in a parasitic flatworm, demonstrating the activity and utility of the process for functional genomics in these pathogens. The striking clinical phenotype highlights the role in virulence that liver fluke growth factors play in biliary tract morbidity during chronic opisthorchiasis.

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

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

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

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

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

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

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

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

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