Improving helminth genome resources in the post-genomic era

Systems biology of Haemonchus contortus - Advancing biotechnology for parasitic nematode control

Parasitic nematodes represent a substantial global burden, impacting animal health, agriculture and economies worldwide. Of these worms, Haemonchus contortus - a blood-feeding nematode of ruminants - is a major pathogen and a model for molecular and applied parasitology research. This review synthesises some key advances in understanding the molecular biology, genetic diversity and host-parasite interactions of H. contortus, highlighting its value for comparative studies with the free-living nematode Caenorhabditis elegans. Key themes include recent developments in genomic, transcriptomic and proteomic technologies and resources, which are illuminating critical molecular pathways, including the ubiquitination pathway, protease/protease inhibitor systems and the secretome of H. contortus. Some of these insights are providing a foundation for identifying essential genes and exploring their potential as targets for novel anthelmintics or vaccines, particularly in the face of widespread anthelmintic resistance. Advanced bioinformatic tools, such as machine learning (ML) algorithms and artificial intelligence (AI)-driven protein structure prediction, are enhancing annotation capabilities, facilitating and accelerating analyses of gene functions, and biological pathways and processes. This review also discusses the integration of these tools with cutting-edge single-cell sequencing and spatial transcriptomics to dissect host-parasite interactions at the cellular level. The discussion emphasises the importance of curated databases, improved culture systems and functional genomics platforms to translate molecular discoveries into practical outcomes, such as novel interventions. New research findings and resources not only advance research on H. contortus and related nematodes but may also pave the way for innovative solutions to the global challenges with anthelmintic resistance.

Antimicrobial peptides in nematode secretions - Unveiling biotechnological opportunities for therapeutics and beyond

Gastrointestinal (GI) parasitic nematodes threaten food security and affect human health and animal welfare globally. Current anthelmintics for use in humans and livestock are challenged by continuous re-infections and the emergence and spread of multidrug resistance, underscoring an urgent need to identify novel control targets for therapeutic exploitation. Recent evidence has highlighted the occurrence of complex interplay between GI parasitic nematodes of humans and livestock and the resident host gut microbiota. Antimicrobial peptides (AMPs) found within nematode biofluids have emerged as potential effectors of these interactions. This review delves into the occurrence, structure, and function of nematode AMPs, highlighting their potential as targets for drug discovery and development. We argue that an integrated approach combining advanced analytical techniques, scalable production methods, and innovative experimental models is needed to unlock the full potential of nematode AMPs and pave the way for the discovery and development of sustainable parasite control strategies.

Gastrointestinal nematode parasites of grazing ruminants: a comprehensive literature review of diagnostic methods for quantifying parasitism, larval differentiation and measuring anthelmintic resistance

This review summarises up-to-date research on the diagnosis of gastrointestinal nematode (GIN) infection in livestock and anthelmintic resistance in GIN. It was commissioned to assist funding bodies to prioritise and guide research and extension efforts to improve the health, welfare and productivity of grazing ruminants in the face of challenge with GIN. A comprehensive review of published articles from journals, books and websites was undertaken, with a focus on peer-reviewed articles published between 2000 and 2024 involving genera of GIN in grazing sheep and cattle with economic importance to New Zealand. Suggestions for articles to include were received from 14 experts in GIN diagnostics. This review is a summary of a longer report submitted to the sponsoring organisation. Clinical signs of GIN infection in grazing ruminants in temperate grazing systems are inadequate as triggers for management interventions including anthelmintic treatment as they are visible only after economically significant pathological changes have occurred. Livestock producers benefit from monitoring GIN burdens using faecal egg counts (FEC) or associated signals such as weight gain. In future, they may use remote monitoring devices for activity in animals, as well as estimating pasture larval contamination. Methods of diagnosing GIN infections using automated FEC devices have improved the convenience of monitoring parasite burdens compared with traditional laboratory methods. However, a lack of quality control measures and a gap in training of skilled technicians for larval differentiation may lead to a shortage of diagnostic capability. Current methods of diagnosing anthelmintic resistance, particularly FEC reduction tests, are not likely to be replaced by laboratory assays in the near future and attention should be focused on facilitating application of new FEC technologies for both animal monitoring and resistance diagnosis. Extension and application of currently available methods and technology will improve animal health and productivity in ruminant grazing systems in the short term. Adoption of novel technologies for remote animal monitoring, practical tools for estimating pasture larval contamination and promoting genetic selection for immunity and resilience to GIN in both sheep and cattle will further enhance productivity in the long term.

Genome assembly variation and its implications for gene discovery in nematodes

Genome assemblers are a critical component of genome science, but the choice of assembly software and protocols can be daunting. Here, we investigate genome assembly variation and its implications for gene discovery across three nematode species-Caenorhabditis bovis, Haemonchus contortus, and Heligmosomoides bakeri-highlighting the critical interplay between assembly choice and downstream genomic analysis. Selecting commonly used genome assemblers, we generated multiple assemblies for each species, analyzing their structure, completeness, and effect on gene family analysis. Our findings demonstrate that assembly variations can significantly affect gene family composition, with notable differences in gene families important in anthelmintic discovery and immunomodulation. Despite broadly similar performance using various assembly metrics, comparisons of assemblies with a single species revealed underlying structural rearrangements and inconsistencies in gene content, which would affect downstream analyses. This emphasizes the need for continuous refinement of genome assemblies and their annotations.

Parasite spillover rather than niche expansion explains infection of host brain by diplostomid eye flukes

Parasites often occupy specific sites within their host, which has important implications for host performance and parasite transmission. Nonetheless, parasitic infections can occur beyond their typical location within a host, significantly altering host-parasite interactions. Yet, the causes behind the atypical tissue tropism are poorly understood. Here, we focus on a ubiquitous group of diplostomid parasites that form diverse communities in fish eyes. We used targeted DNA metabarcoding (cytochrome c oxydase subunit 1, COX1, 250 bp) to evaluate potential mechanisms underlying eye parasite atypical tissue tropism to the brain of two widespread fish species (Eurasian perch and common roach). We found that the most common eye-infecting species (Tylodelphys clavata, Diplostomum baeri) are present in the brains of perch but not in roach. The bipartite network comprising 5 species and 24 mitochondrial haplotypes revealed no brain-specific haplotypes, indicating an apparent lack of genetic divergence between brain- and eye-infecting parasites. Instead, the prevalence, intensity and diversity of eye infections were positively correlated with brain infections. Thus, our results suggest that the most parsimonious mechanism underlying brain infection is density-dependent spillover rather than parasite divergence-driven niche expansion. We anticipate that 'off-target' infections are likely to be severely underestimated in nature with important ecological, evolutionary and medical implications.

Chromosomal genome assembly resolves drug resistance loci in the parasitic nematode Teladorsagia circumcincta

The parasitic nematode Teladorsagia circumcincta is one of the most important pathogens of sheep and goats in temperate climates worldwide and can rapidly evolve resistance to drugs used to control it. To understand the genetics of drug resistance, we have generated a highly contiguous genome assembly for the UK T. circumcincta isolate, MTci2. Assembly using PacBio long-reads and Hi-C long-molecule scaffolding together with manual curation resulted in a 573 Mb assembly (N50 = 84 Mb, total scaffolds = 1,286) with five autosomal and one sex-linked chromosomal-scale scaffolds consistent with its karyotype. The genome resource was further improved via annotation of 22,948 genes, with manual curation of over 3,200 of these, resulting in a robust and near complete resource (96.3% complete protein BUSCOs) to support basic and applied research on this important veterinary pathogen. Genome-wide analyses of drug resistance, combining evidence from three distinct experiments, identified selection around known candidate genes for benzimidazole, levamisole and ivermectin resistance, as well as novel regions associated with ivermectin and moxidectin resistance. These insights into contemporary and historic genetic selection further emphasise the importance of contiguous genome assemblies in interpreting genome-wide genetic variation associated with drug resistance and identifying key loci to prioritise in developing diagnostic markers of anthelmintic resistance to support parasite control.

Entering the spatial age of parasite genomics

Reconciling organism-scale biology at cellular-scale resolution has been a monumental challenge. Recently, Gramberg et al. uncovered the spatial organisation of gene expression in the common liver fluke, Fasciola hepatica. These data provide new insights and opportunities for understanding the fundamental developmental and functional biology of parasites and new targets for control.

Analysis of Haemonchus embryos at single cell resolution identifies two eukaryotic elongation factors as intervention target candidates

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.

Deep-amplicon sequencing of the complete beta-tubulin gene in Trichuris trichiura before and after albendazole treatment

Concerns about the emergence of benzimidazole resistance in soil-transmitted helminths (STH) infections, particularly against Trichuris trichiura, have arisen. Previous studies of veterinary nematodes have linked benzimidazole resistance to single-nucleotide polymorphisms (SNPs) at three specific codons in the beta-tubulin gene, but similar associations in STH have not been consistently observed. In this work, we screened the complete beta-tubulin gene previously linked to benzimidazole resistance in T. trichiura by deep-amplicon sequencing to identify genetic variants and associate levels of diversity with drug response to albendazole. We used 99 DNA samples extracted from T. trichiura pooled eggs, previously semi-purified from human stool samples collected in Manhiça district, Mozambique. We obtained a set of 39 amplicons of the complete gene by subjecting the pooled eggs to long-read PCR and subsequently sequencing them. Of those amplicons, 22 and 17 were obtained from stool samples collected before, and 21 days after albendazole treatment, respectively. We observed genetic variation across the whole gene sequence, in both exons and introns; however, none were associated with the previously proposed resistance-associated SNPs, and none were predicted to significantly affect protein function. No significant differences in genetic diversity were observed between pre- and post-treatment samples. Using publicly available genome-wide data, we also analysed a second beta-tubulin isotype in the T. trichiura genome. We focused on detecting the canonical SNPs and assessing for signatures of genetic selection around this second isotype gene. This analysis did not reveal evidence supporting this second isotype's role in anthelmintic resistance. Despite the limitations of our study, such as a small sample size, particularly paired pre- and post-treatment samples (n = 6), or a restricted geographical area, we found no evidence linking either of the two beta-tubulin genes to benzimidazole resistance in T. trichiura, suggesting that genetic markers of drug resistance likely exist outside the beta-tubulin genes.

Chromosome-contiguous genome for the Haecon-5 strain of Haemonchus contortus reveals marked genetic variability and enables the discovery of essential gene candidates

Millions of livestock animals worldwide are infected with the haematophagous barber's pole worm, Haemonchus contortus, the aetiological agent of haemonchosis. Despite the major significance of this parasite worldwide and its widespread resistance to current treatments, the lack of a high-quality genome for the well-defined strain of this parasite from Australia, called Haecon-5, has constrained research in a number of areas including host-parasite interactions, drug discovery and population genetics. To enable research in these areas, we report here a chromosome-contiguous genome (∼280 Mb) for Haecon-5 with high-quality models for 19,234 protein-coding genes. Comparative genomic analyses show significant genomic similarity (synteny) with a UK strain of H. contortus, called MHco3(ISE).N1 (abbreviated as "ISE"), but we also discover marked differences in genomic structure/gene arrangements, distribution of nucleotide variability (single nucleotide polymorphisms (SNPs) and indels) and orthology between Haecon-5 and ISE. We used the genome and extensive transcriptomic resources for Haecon-5 to predict a subset of essential single-copy genes employing a "cross-species" machine learning (ML) approach using a range of features from nucleotide/protein sequences, protein orthology, subcellular localisation, single-cell RNA-seq and/or histone methylation data available for the model organisms Caenorhabditis elegans and Drosophila melanogaster. From a set of 1,464 conserved single copy genes, transcribed in key life-cycle stages of H. contortus, we identified 232 genes whose homologs have critical functions in C. elegans and/or D. melanogaster, and prioritised 10 of them for further characterisation; nine of the 10 genes likely play roles in neurophysiological processes, germline, hypodermis and/or respiration, and one is an unknown (orphan) gene for which no detailed functional information exists. Future studies of these genes/gene products are warranted to elucidate their roles in parasite biology, host-parasite interplay and/or disease. Clearly, the present Haecon-5 reference genome and associated resources now underpin a broad range of fundamental investigations of H. contortus and could assist in accelerating the discovery of novel intervention targets and drug candidates to combat haemonchosis.

Assessing de novo parasite genomes assembled using only Oxford Nanopore Technologies MinION data

In this study, we assessed the quality of de novo genome assemblies for three species of parasitic nematodes (Brugia malayi, Trichuris trichiura, and Ancylostoma caninum) generated using only Oxford Nanopore Technologies MinION data. Assemblies were compared to current reference genomes and against additional assemblies that were supplemented with short-read Illumina data through polishing or hybrid assembly approaches. For each species, assemblies generated using only MinION data had similar or superior measures of contiguity, completeness, and gene content. In terms of gene composition, depending on the species, between 88.9 and 97.6% of complete coding sequences predicted in MinION data only assemblies were identical to those predicted in assemblies polished with Illumina data. Polishing MinION data only assemblies with Illumina data therefore improved gene-level accuracy to a degree. Furthermore, modified DNA extraction and library preparation protocols produced sufficient genomic DNA from B. malayi and T. trichiura to generate de novo assemblies from individual specimens.

Comparison of multi-parallel quantitative real-time PCRs targeting different DNA regions and detecting soil-transmitted helminths in stool

BACKGROUND

Soil-transmitted helminths infect an estimated 18% of the world's population, causing a significant health burden. Microscopy has been the primary tool for diagnosing eggs from fecal samples, but its sensitivity drops in low-prevalence settings. Quantitative real-time polymerase chain reaction (qPCR) is slowly increasing in research and clinical settings. However, there is still no consensus on preferred qPCR targets.

METHODS

We aimed to compare soil-transmitted helminth (STH) DNA detection methods by testing naïve stool samples spiked with known quantities of STH eggs and larvae. DNA extracts from spiked samples were tested using independent quantitative realtime PCR (qPCR) assays targeting ribosomal or putative non-protein coding satellite sequences.

RESULTS

For Trichuris trichiura, there was a strong correlation between egg/larvae counts and qPCR results using either qPCR method (0.86 and 0.87, respectively). Strong correlations also existed for A. lumbricoides (0.60 and 0.63, respectively), but weaker correlations were found for Ancylostoma duodenale (0.41 for both assays) and Strongyloides stercoralis (0.48 and 0.65, respectively). No correlation for Necator americanus was observed when testing with either qPCR assay. Both assays had fair-to-moderate agreement across targets when using field-collected stool samples (0.28-0.45, for all STHs), except for S. stercoralis (0.12) with slight agreement.

CONCLUSIONS

There is a strong correlation between qPCR results and egg/larvae counts. Our study confirms that qPCR is an effective diagnostic tool, even with low-intensity infections, regardless of the DNA-based diagnostic marker used. However, the moderate agreement between the two different qPCR assays when testing field samples highlights the need to understand the role of these targets in the genome so that the parasite burden can be quantified more accurately and consistently by qPCR.

Spirocerca lupi draft genome, vaccine and anthelmintic targets

Spirocerca lupi is a parasitic nematode affecting predominantly domestic dogs. It causes spirocercosis, a disease that is often fatal. The assembled draft genome of S. lupi consists of 13,627 predicted protein-coding genes and is approximately 150 Mb in length. Several known anthelmintic gene targets such as for β-Tubulin, glutamate, and GABA receptors as well as known vaccine gene targets such as cysteine protease inhibitor and cytokines were identified in S. lupi by comparing orthologs of C. elegans anthelmintic gene targets as well as orthologs to known vaccine candidates. New anthelmintic targets were predicted through an inclusion-exclusion strategy and new vaccine targets were predicted through an immunoinformatics approach. New anthelminthic targets include DNA-directed RNA polymerases, chitin synthase, polymerases, and other enzymes. New vaccine targets include cuticle collagens. These gene targets provide a starting platform for new drug identification and vaccine design.

Inference of Essential Genes of the Parasite Haemonchus contortus via Machine Learning

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.

Hookworm genomics: dusk or dawn?

Hookworms are parasites, closely related to the model nematode Caenorhabditis elegans, that are a major economic and health burden worldwide. Primarily three hookworm species (Necator americanus, Ancylostoma duodenale, and Ancylostoma ceylanicum) infect humans. Another 100 hookworm species from 19 genera infect primates, ruminants, and carnivores. Genetic data exist for only seven of these species. Genome sequences are available from only four of these species in two genera, leaving 96 others (particularly those parasitizing wildlife) without any genomic data. The most recent hookworm genomes were published 5 years ago, leaving the field in a dusk. However, assembling genomes from single hookworms may bring a new dawn. Here we summarize advances, challenges, and opportunities for studying these neglected but important parasitic nematodes.

A single locus determines praziquantel response in Schistosoma mansoni

We previously performed a genome-wide association study (GWAS) to identify the genetic basis of praziquantel (PZQ) response in schistosomes, identifying two quantitative trait loci situated on chromosomes 2 and 3. We reanalyzed this GWAS using the latest (version 10) genome assembly showing that a single locus on chromosome 3, rather than two independent loci, determines drug response. These results reveal that PZQ response is monogenic and demonstrates the importance of high-quality genomic information.

Subtelomeric plasticity contributes to gene family expansion in the human parasitic flatworm Schistosoma mansoni

BACKGROUND

The genomic region that lies between the telomere and chromosome body, termed the subtelomere, is heterochromatic, repeat-rich, and frequently undergoes rearrangement. Within this region, large-scale structural changes enable gene diversification, and, as such, large multicopy gene families are often found at the subtelomere. In some parasites, genes associated with proliferation, invasion, and survival are often found in these regions, where they benefit from the subtelomere's highly plastic, rapidly changing nature. The increasing availability of complete (or near complete) parasite genomes provides an opportunity to investigate these typically poorly defined and overlooked genomic regions and potentially reveal relevant gene families necessary for the parasite's lifestyle.

RESULTS

Using the latest chromosome-scale genome assembly and hallmark repeat richness observed at chromosome termini, we have identified and characterised the subtelomeres of Schistosoma mansoni, a metazoan parasitic flatworm that infects over 250 million people worldwide. Approximately 12% of the S. mansoni genome is classified as subtelomeric, and, in line with other organisms, we find these regions to be gene-poor but rich in transposable elements. We find that S. mansoni subtelomeres have undergone extensive interchromosomal recombination and that these sites disproportionately contribute to the 2.3% of the genome derived from segmental duplications. This recombination has led to the expansion of subtelomeric gene clusters containing 103 genes, including the immunomodulatory annexins and other gene families with unknown roles. The largest of these is a 49-copy plexin domain-containing protein cluster, exclusively expressed in the tegument-the tissue located at the host-parasite physical interface-of intramolluscan life stages.

CONCLUSIONS

We propose that subtelomeric regions act as a genomic playground for trial-and-error of gene duplication and subsequent divergence. Owing to the importance of subtelomeric genes in other parasites, gene families implicated in this subtelomeric expansion within S. mansoni warrant further characterisation for a potential role in parasitism.

Quantifying the fitness effects of resistance alleles with and without anthelmintic selection pressure using Caenorhabditis elegans

Albendazole and ivermectin are the two most commonly co-administered anthelmintic drugs in mass-drug administration programs worldwide. Despite emerging resistance, we do not fully understand the mechanisms of resistance to these drugs nor the consequences of delivering them in combination. Albendazole resistance has primarily been attributed to variation in the drug target, a beta-tubulin gene. Ivermectin targets glutamate-gated chloride channel (GluCl) genes, but it is unknown whether these genes are involved in ivermectin resistance in nature. Using Caenorhabditis elegans, we defined the fitness costs associated with loss of the drug target genes singly or in combinations of the genes that encode GluCl subunits. We quantified the loss-of function effects on three traits: (i) multi-generational competitive fitness, (ii) fecundity, and (iii) development. In competitive fitness and development assays, we found that a deletion of the beta-tubulin gene ben-1 conferred albendazole resistance, but ivermectin resistance required loss of two GluCl genes (avr-14 and avr-15) or loss of three GluCl genes (avr-14, avr-15, and glc-1). The fecundity assays revealed that loss of ben-1 did not provide any fitness benefit in albendazole and that no GluCl deletion mutants were resistant to ivermectin. Next, we searched for evidence of multi-drug resistance across the three traits. Loss of ben-1 did not confer resistance to ivermectin, nor did loss of any single GluCl subunit or combination confer resistance to albendazole. Finally, we assessed the development of 124 C. elegans wild strains across six benzimidazoles and seven macrocyclic lactones to identify evidence of multi-drug resistance between the two drug classes and found a strong phenotypic correlation within a drug class but not across drug classes. Because each gene affects various aspects of nematode physiology, these results suggest that it is necessary to assess multiple fitness traits to evaluate how each gene contributes to anthelmintic resistance.

Genome structure and population genomics of the canine heartworm Dirofilaria immitis

The heartworm, Dirofilaria immitis, is a filarial parasitic nematode responsible for significant morbidity and mortality in wild and domesticated canids. Resistance to macrocyclic lactone drug prevention represents a significant threat to parasite control and has prompted investigations to understand the genetic determinants of resistance. This study aimed to improve the genomic resources of D. immitis to enable a more precise understanding of how genetic variation is distributed within and between parasite populations worldwide, which will inform the likelihood and rate by which parasites, and in turn, resistant alleles, might spread. We have guided the scaffolding of a recently published genome assembly for D. immitis (ICBAS_JMDir_1.0) using the chromosomal-scale reference genomes of Brugia malayi and Onchocerca volvulus, resulting in an 89.5 Mb assembly composed of four autosomal- and one sex-linked chromosomal-scale scaffolds representing 99.7% of the genome. Publicly available and new whole-genome sequencing data from 32 D. immitis samples from Australia, Italy and the USA were assessed using principal component analysis, nucleotide diversity (Pi) and absolute genetic divergence (Dxy) to characterise the global genetic structure and measure within- and between-population diversity. These population genetic analyses revealed broad-scale genetic structure among globally diverse samples and differences in genetic diversity between populations; however, fine-scale subpopulation analysis was limited and biased by differences between sample types. Finally, we mapped single nucleotide polymorphisms previously associated with macrocyclic lactone resistance in the new genome assembly, revealing the physical linkage of high-priority variants on chromosome 3, and determined their frequency in the studied populations. This new chromosomal assembly for D. immitis now allows for a more precise investigation of selection on genome-wide genetic variation and will enhance our understanding of parasite transmission and the spread of genetic variants responsible for resistance to treatment.

Advancing Strongyloides omics data: bridging the gap with Caenorhabditis elegans

Among nematodes, the free-living model organism Caenorhabditis elegans boasts the most advanced portfolio of high-quality omics data. The resources available for parasitic nematodes, including Strongyloides spp., however, are lagging behind. While C. elegans remains the most tractable nematode and has significantly advanced our understanding of many facets of nematode biology, C. elegans is not suitable as a surrogate system for the study of parasitism and it is important that we improve the omics resources available for parasitic nematode species. Here, we review the omics data available for Strongyloides spp. and compare the available resources to those for C. elegans and other parasitic nematodes. The advancements in C. elegans omics offer a blueprint for improving omics-led research in Strongyloides. We suggest areas of priority for future research that will pave the way for expansions in omics resources and technologies. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

A standard workflow for community-driven manual curation of Strongyloides genome annotations

Advances in the functional genomics and bioinformatics toolkits for Strongyloides species have positioned these species as genetically tractable model systems for gastrointestinal parasitic nematodes. As community interest in mechanistic studies of Strongyloides species continues to grow, publicly accessible reference genomes and associated genome annotations are critical resources for researchers. Genome annotations for multiple Strongyloides species are broadly available via the WormBase and WormBase ParaSite online repositories. However, a recent phylogenetic analysis of the receptor-type guanylate cyclase (rGC) gene family in two Strongyloides species highlights the potential for errors in a large percentage of current Strongyloides gene models. Here, we present three examples of gene annotation updates within the Strongyloides rGC gene family; each example illustrates a type of error that may occur frequently within the annotation data for Strongyloides genomes. We also extend our analysis to 405 previously curated Strongyloides genes to confirm that gene model errors are found at high rates across gene families. Finally, we introduce a standard manual curation workflow for assessing gene annotation quality and generating corrections, and we discuss how it may be used to facilitate community-driven curation of parasitic nematode biodata. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

Trematode Genomics and Proteomics

Trematode infections stand out as one of the frequently overlooked tropical diseases, despite their wide global prevalence and remarkable capacity to parasitize diverse host species and tissues. Furthermore, these parasites hold significant socio-economic, medical, veterinary and agricultural implications. Over the past decades, substantial strides have been taken to bridge the information gap concerning various "omic" tools, such as proteomics and genomics, in this field. In this edition of the book, we highlight recent progress in genomics and proteomics concerning trematodes with a particular focus on the advances made in the past 5 years. Additionally, we present insights into cutting-edge technologies employed in studying trematode biology and shed light on the available resources for exploring the molecular facets of this particular group of parasitic helminths.

A single locus determines praziquantel response in Schistosoma mansoni

We previously performed a genome-wide association study (GWAS) to identify the genetic basis of praziquantel (PZQ) response in schistosomes, identifying two quantitative trait loci (QTL) situated on chromosome 2 and chromosome 3. We reanalyzed this GWAS using the latest (v10) genome assembly showing that a single locus on chromosome 3, rather than two independent loci, determines drug response. These results reveal that praziquantel response is monogenic and demonstrates the importance of high-quality genomic information.

Fecundity, in vitro early larval development and karyotype of the zoonotic nematode Anisakis pegreffii

The in vitro life cycle of zoonotic helminths is an essential tool for -omic translational studies focused on disease control and treatment. Anisakiosis is an emerging zoonosis contracted by the ingestion of raw or undercooked fish infected with the third stage larvae (L3) of two sibling species Anisakis simplex sensu stricto (s.s.) and Anisakis pegreffii, the latter being the predominant species in the Mediterranean basin. Recently, in vitro culture of A. pegreffii has been developed to enable fast and large-scale production of fertile adults. However, the conditions for larval development from hatching to infective L3 were not fulfilled to complete the cycle. Herein, we used a Drosophila medium supplemented with chicken serum and adjusted different osmolarities to maintain the culture of L3 hatched from eggs for up to 17 weeks. The highest survival rate was observed in the medium with the highest osmolarities, which also allowed the highest larval exsheathment rate. Key morphological features of embryogenesis and postembryogenesis studied by transmission electron microscopy revealed that the excretory gland cell is differentiated already up to 48 h post-hatching. Extracellular vesicles and cell-free mitochondria are discharged between the two cuticle sheets of the second stage larvae (L2). Contemporarly cultivated, two populations of adult A. simplex s.s. and A. pegreffii reached an average production of 29,914.05 (± 27,629.36) and 24,370.96 (± 12,564.86) eggs/day/female, respectively. The chromosome spreads of A. pegreffii obtained from mature gonads suggests a diploid karyotype formula of 2n = 18. The development of a reliable protocol for the in vitro culture of a polyxenous nematode such as Anisakis spp. will serve to screen for much needed novel drug targets, but also to study the intricated and unknown ecological and physiological traits of these trophically transmitted marine nematodes.

Revolutionizing Medical Microbiology: How Molecular and Genomic Approaches Are Changing Diagnostic Techniques

Molecular and genomic approaches have revolutionized medical microbiology by offering faster and more accurate diagnostic techniques for infectious diseases. Traditional methods, which include culturing microbes and biochemical testing, are time-consuming and may not detect antibiotic-resistant strains. In contrast, molecular and genomic methods, including polymerase chain reaction (PCR)-based techniques and whole-genome sequencing, provide rapid and precise detection of pathogens, early-stage diseases, and antibiotic-resistant strains. These approaches have advantages such as high sensitivity and specificity, the potential for targeted therapies, and personalized medicine. However, implementing molecular and genomic techniques faces challenges related to cost, equipment, expertise, and data analysis. Ethical and legal considerations regarding patient privacy and genetic data usage also arise. Nonetheless, the future of medical microbiology lies in the widespread adoption of molecular and genomic approaches, which can lead to improved patient outcomes and the identification of antibiotic-resistant strains. Continued advancements, education, and exploration of ethical implications are necessary to fully harness the potential of molecular and genomic techniques in medical microbiology.

Exploring the antimicrobial peptidome of nematodes through phylum-spanning in silico analyses highlights novel opportunities for pathogen control

Antimicrobial Peptides (AMPs) are key constituents of the invertebrate innate immune system and provide critical protection against microbial threat. Nematodes display diverse life strategies where they are exposed to heterogenous, microbe rich, environments highlighting their need for an innate immune system. Within the Ecdysozoa, arthropod AMPs have been well characterised, however nematode-derived AMP knowledge is limited. In this study the distribution and abundance of putative AMP-encoding genes was examined in 134 nematode genomes providing the most comprehensive profile of AMP candidates within phylum Nematoda. Through genome and transcriptome analyses we reveal that phylum Nematoda is a rich source of putative AMP diversity and demonstrate (i) putative AMP group profiles that are influenced by nematode lifestyle where free-living nematodes appear to display enriched putative AMP profiles relative to parasitic species; (ii) major differences in the putative AMP profiles between nematode clades where Clade 9/V and 10/IV species possess expanded putative AMP repertoires; (iii) AMP groups with highly restricted profiles (e.g. Cecropins and Diapausins) and others [e.g. Nemapores and Glycine Rich Secreted Peptides (GRSPs)] which are more widely distributed; (iv) complexity in the distribution and abundance of CSαβ subgroup members; and (v) that putative AMPs are expressed in host-facing life stages and biofluids of key nematode parasites. These data indicate that phylum Nematoda displays diversity in putative AMPs and underscores the need for functional characterisation to reveal their role and importance to nematode biology and host-nematode-microbiome interactions.

An insight into the functional genomics and species classification of Eudiplozoon nipponicum (Monogenea, Diplozoidae), a haematophagous parasite of the common carp Cyprinus carpio

BACKGROUND

Monogenea (Platyhelminthes, Neodermata) are the most species-rich class within the Neodermata superclass of primarily fish parasites. Despite their economic and ecological importance, monogenean research tends to focus on their morphological, phylogenetic, and population characteristics, while comprehensive omics analyses aimed at describing functionally important molecules are few and far between. We present a molecular characterisation of monogenean representative Eudiplozoon nipponicum, an obligate haematophagous parasite infecting the gills of the common carp. We report its nuclear and mitochondrial genomes, present a functional annotation of protein molecules relevant to the molecular and biochemical aspect of physiological processes involved in interactions with the fish hosts, and re-examinate the taxonomic position of Eudiplozoon species within the Diplozoidae family.

RESULTS

We have generated 50.81 Gbp of raw sequencing data (Illumina and Oxford Nanopore reads), bioinformatically processed, and de novo assembled them into a genome draft 0.94 Gbp long, consisting of 21,044 contigs (N50 = 87 kbp). The final assembly represents 57% of the estimated total genome size (~ 1.64 Gbp), whereby repetitive and low-complexity regions account for ~ 64% of the assembled length. In total, 36,626 predicted genes encode 33,031 proteins and homology-based annotation of protein-coding genes (PCGs) and proteins characterises 14,785 (44.76%) molecules. We have detected significant representation of functional proteins and known molecular functions. The numbers of peptidases and inhibitors (579 proteins), characterised GO terms (16,016 unique assigned GO terms), and identified KEGG Orthology (4,315 proteins) acting in 378 KEGG pathways demonstrate the variety of mechanisms by which the parasite interacts with hosts on a macromolecular level (immunomodulation, feeding, and development). Comparison between the newly assembled E. nipponicum mitochondrial genome (length of 17,038 bp) and other diplozoid monogeneans confirms the existence of two distinct Eudiplozoon species infecting different fish hosts: Cyprinus carpio and Carassius spp.

CONCLUSIONS

Although the amount of sequencing data and characterised molecules of monogenean parasites has recently increased, a better insight into their molecular biology is needed. The E. nipponicum nuclear genome presented here, currently the largest described genome of any monogenean parasite, represents a milestone in the study of monogeneans and their molecules but further omics research is needed to understand these parasites' biological nature.

Helminthic host defense peptides: using the parasite to defend the host

Parasitic helminths are destined to share niches with a variety of microbiota that inevitably influence their interaction with the host. To modulate the microbiome for their benefit and defend against pathogenic isolates, helminths have developed host defense peptides (HDPs) and proteins as integral elements of their immunity. These often exert a relatively nonspecific membranolytic activity toward bacteria, sometimes with limited or no toxicity toward host cells. With a few exceptions, such as nematode cecropin-like peptides and antibacterial factors (ABFs), helminthic HDPs are largely underexplored. This review scrutinizes current knowledge on the repertoire of such peptides in helminths and promotes their research as potential leads for an anti-infective solution to the burgeoning problem of antibiotic resistance.

Evaluation of genome skimming to detect and characterise human and livestock helminths

The identification of gastrointestinal helminth infections of humans and livestock almost exclusively relies on the detection of eggs or larvae in faeces, followed by manual counting and morphological characterisation to differentiate species using microscopy-based techniques. However, molecular approaches based on the detection and quantification of parasite DNA are becoming more prevalent, increasing the sensitivity, specificity and throughput of diagnostic assays. High-throughput sequencing, from single PCR targets through to the analysis of whole genomes, offers significant promise towards providing information-rich data that may add value beyond traditional and conventional molecular approaches; however, thus far, its utility has not been fully explored to detect helminths in faecal samples. In this study, low-depth whole genome sequencing, i.e. genome skimming, has been applied to detect and characterise helminth diversity in a set of helminth-infected human and livestock faecal material. The strengths and limitations of this approach are evaluated using three methods to characterise and differentiate metagenomic sequencing data based on (i) mapping to whole mitochondrial genomes, (ii) whole genome assemblies, and (iii) a comprehensive internal transcribed spacer 2 (ITS2) database, together with validation using quantitative PCR (qPCR). Our analyses suggest that genome skimming can successfully identify most single and multi-species infections reported by qPCR and can provide sufficient coverage within some samples to resolve consensus mitochondrial genomes, thus facilitating phylogenetic analyses of selected genera, e.g. Ascaris spp. Key to this approach is both the availability and integrity of helminth reference genomes, some of which are currently contaminated with bacterial and host sequences. The success of genome skimming of faecal DNA is dependent on the availability of vouchered sequences of helminths spanning both taxonomic and geographic diversity, together with methods to detect or amplify minute quantities of parasite nucleic acids in mixed samples.

Genomic landscape of drug response reveals mediators of anthelmintic resistance

Like other pathogens, parasitic helminths can rapidly evolve resistance to drug treatment. Understanding the genetic basis of anthelmintic drug resistance in parasitic nematodes is key to tracking its spread and improving the efficacy and sustainability of parasite control. Here, we use an in vivo genetic cross between drug-susceptible and multi-drug-resistant strains of Haemonchus contortus in a natural host-parasite system to simultaneously map resistance loci for the three major classes of anthelmintics. This approach identifies new alleles for resistance to benzimidazoles and levamisole and implicates the transcription factor cky-1 in ivermectin resistance. This gene is within a locus under selection in ivermectin-resistant populations worldwide; expression analyses and functional validation using knockdown experiments support that cky-1 is associated with ivermectin survival. Our work demonstrates the feasibility of high-resolution forward genetics in a parasitic nematode and identifies variants for the development of molecular diagnostics to combat drug resistance in the field.