Genomic and transcriptomic variation defines the chromosome-scale assembly of Haemonchus contortus, a model gastrointestinal worm

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.

Visualizing genomic evolution in Caenorhabditis through WormSynteny

Understanding the syntenic relationships among genomes is crucial to elucidate the genomic mechanisms that drive the evolution of species. The nematode Caenorhabditis is a good model for studying genomic evolution due to the well-established biology of Caenorhabditis elegans and the availability of > 50 genomes in the genus. However, effective alignment of more than ten species in Caenorhabditis has not been conducted before, and there is currently no tool to visualize the synteny of more than two species. In this study, we used Progressive Cactus, a recently developed multigenome aligner, to align the genomes of eleven Caenorhabditis species. Through the progressive alignment, we reconstructed nine ancestral genomes, analyzed the mutational types that cause genomic rearrangement during speciation, and found that insertion and duplication are the major driving forces for genome expansion. Dioecious species appear to expand their genomes more than androdioecious species. We then built an online interactive app called WormSynteny to visualize the syntenic relationship among the eleven species. Users can search the alignment dataset using C. elegans query sequences, construct synteny plots at different genomic scales, and use a set of options to control alignment output and plot presentation. We showcased the use of WormSynteny to visualize the syntenic conservation of one-to-one orthologues among species, tandem and dispersed gene duplication in C. elegans, and the evolution of exon and intron structures. Importantly, the integration of orthogroup information with synteny linkage in WormSynteny allows the easy visualization of conserved genomic blocks and disruptive rearrangement. In conclusion, WormSynteny provides immediate access to the syntenic relationships among the most widely used Caenorhabditis species and can facilitate numerous comparative genomics studies. This pilot study with eleven species also serves as a proof-of-concept to a more comprehensive larger-scale analysis using hundreds of nematode genomes, which is expected to reveal mechanisms that drive genomic evolution in the Nematoda phylum. Finally, the WormSynteny software provides a generalizable solution for visualizing the output of Progressive Cactus with interactive graphics, which would be useful for a broad community of genome researchers.

Molecular study of the transcription factor SKN-1 and its putative relationship with genes that encode GST and antioxidant enzymes in Haemonchus contortus

Haemonchus contortus is a parasitic nematode of ruminants. Once inside its host, it is exposed to reactive oxidative species and responds by synthesising antioxidant enzymes as a defence. In Caenorhabditis elegans, antioxidant genes are regulated by the transcription factor skinhead-1 (Cel-SKN-1). However, there is little information about the function of SKN-1 in H. contortus (Hco-SKN-1). Therefore, we performed a molecular investigation to characterise Hco-SKN-1 and its putative relationship with genes encode antioxidant enzymes, namely glutathione S-transferases (Hco-GSTs, n = 3), superoxide dismutase (Hco-SOD) and catalase (Hco-CAT), which are involved in haematophagy and defence against the host. We used in silico sequence analysis of Hco-SKN-1 and Hco-GSTs to design and perform relative expression assays involving H. contortus eggs, infective larvae (L3) and adults. Furthermore, we exposed H. contortus transitional infective larvae (xL3) to erythrocytes or hydrogen peroxide (H2O2) and evaluated the relative expression of antioxidant genes at 24 or 48 h. Gene Ontology (GO) analysis revealed 31 functions associated with Hco-SKN-1 and Hco-GSTs, including stress resistance, larval development and the active immune response. Hco-GST-5957 and Hco-SOD showed the highest expression in adults, indicating a relationship with specific functions at this mature stage. xL3 exposed to erythrocytes or H2O2 showed significant upregulation of Hco-SKN-1, but it occurred after upregulation of the antioxidant genes, indicating that these genes are not regulated by Hco-SKN-1 during the blood-feeding stage. Additional investigation is necessary to understand the putative regulation of antioxidant genes by Hco-SKN-1 during the blood-feeding stage.

The mrp-3 gene is involved in haem efflux and detoxification in a blood-feeding nematode

BACKGROUND

Haem is essential but toxic for metazoan organisms. Auxotrophic nematodes can acquire sufficient haem from the environment or their hosts in the meanwhile eliminate or detoxify excessive haem through tightly controlled machinery. In previous work, we reported a role of the unique transporter protein HRG-1 in the haem acquisition and homeostasis of parasitic nematodes. However, little is known about the haem efflux and detoxification via ABC transporters, particularly the multiple drug resistance proteins (MRPs).

RESULTS

Here, we further elucidate that a member of the mrp family (mrp-3) is involved in haem efflux and detoxification in a blood-feeding model gastrointestinal parasite, Haemonchus contortus. This gene is haem-responsive and dominantly expressed in the intestine and inner membrane of the hypodermis of this parasite. RNA interference of mrp-3 resulted in a disturbance of genes (e.g. hrg-1, hrg-2 and gst-1) that are known to be involved in haem homeostasis and an increased formation of haemozoin in the treated larvae and lethality in vitro, particularly when exposed to exogenous haem. Notably, the nuclear hormone receptor NHR-14 appears to be associated the regulation of mrp-3 expression for haem homeostasis and detoxification. Gene knockdown of nhr-14 and/or mrp-3 increases the sensitivity of treated larvae to exogenous haem and consequently a high death rate (> 80%).

CONCLUSIONS

These findings demonstrate that MRP-3 and the associated molecules are essential for haematophagous nematodes, suggesting novel intervention targets for these pathogens in humans and animals.

Multiple UDP glycosyltransferases modulate benzimidazole drug sensitivity in the nematode Caenorhabditis elegans in an additive manner

Xenobiotic biotransformation is an important modulator of anthelmintic drug potency and a potential mechanism of anthelmintic resistance. Both the free-living nematode Caenorhabditis elegans and the ruminant parasite Haemonchus contortus biotransform benzimidazole drugs by glucose conjugation, likely catalysed by UDP-glycosyltransferase (UGT) enzymes. To identify C. elegans genes involved in benzimidazole drug detoxification, we first used a comparative phylogenetic analysis of UGTs from humans, C. elegans and H. contortus, combined with available RNAseq datasets to identify which of the 63 C. elegans ugt genes are most likely to be involved in benzimidazole drug biotransformation. RNA interference knockdown of 15 prioritized C. elegans genes identified those that sensitized animals to the benzimidazole derivative albendazole (ABZ). Genetic mutations subsequently revealed that loss of ugt-9 and ugt-11 had the strongest effects. The "ugt-9 cluster" includes these genes, together with six other closely related ugts. A CRISPR-Cas-9 deletion that removed seven of the eight ugt-9 cluster genes had greater ABZ sensitivity than the single largest-effect mutation. Furthermore, a double mutant of ugt-22 (which is not a member of the ugt-9 cluster) with the ugt-9 cluster deletion further increased ABZ sensitivity. This additivity of mutant phenotypes suggest that ugt genes act in parallel, which could have several, not mutually exclusive, explanations. ugt mutations have different effects with different benzimidazole derivatives, suggesting that enzymes with different specificities could together more efficiently detoxify drugs. Expression patterns of ugt-9, ugt-11 and ugt-22 gfp reporters differ and so likely act in different tissues which may, at least in part, explain their additive effects on drug potency. Overexpression of ugt-9 alone was sufficient to confer partial ABZ resistance, indicating increasing total UGT activity protects animals. In summary, our results suggest that the multiple UGT enzymes have overlapping but not completely redundant functions in benzimidazole drug detoxification and may represent "druggable" targets to improve benzimidazole drug potency.

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.

The Mitogenome of the Haecon-5 Strain of Haemonchus contortus and a Comparative Analysis of Its Nucleotide Variation with Other Laboratory Strains

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

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

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

Haemonchus contortus, an obligatory haematophagus worm infection in small ruminants: Population genetics and genetic diversity

Haemonchus contortus, a stomach worm, is prevalent in ruminants worldwide. They particularly hamper profitable small ruminant production. Here, we estimate the genetic variation of H. contortus collected from slaughtered goats and sheep from various geographic zones of Bangladesh using multiple genes. To perform this, adult parasites were isolated from the abomasum of slaughtered animals (sheep and goats). Among them, 79 male H. contortus were identified by microscopy. Following the extraction of DNA, ITS-2 and cox1 genes were amplified and subsequently considered for sequencing. After alignment and editing, sequences were analyzed to find out sequence variation, diversity pattern of genes, and population genetics of isolates. Among the sequence data, the analyses identified 19 genotypes of ITS-2 and 77 haplotypes of cox1 genes. The diversity of nucleotides was 0.0103 for ITS-2 and 0.029 for cox1 gene. The dendogram constructed by the genotype and haplotype sequences of H. contortus revealed that two populations were circulating in Bangladesh without any demarcation of host and geographic regions. Analysis of population genetics demonstrated a high flow of genes (89.2 %) within the population of the worm in Bangladesh. The Fst value showed very little amount of genetic difference among the worm populations of Bangladesh but marked genetic variation between different continents. The findings are expected to help explain the risks of anthelmintic resistance and the transmission pattern of the parasite, and also provide a control strategy against H. contortus.

Phylogenetic and transcriptomic study of aldo-keto reductases in Haemonchus contortus and their inducibility by flubendazole

Aldo-keto reductases (AKRs), a superfamily of NADP(H)-dependent oxidoreductases, catalyze the oxidoreduction of a wide variety of eobiotic and xenobiotic aldehydes and ketones. In mammals, AKRs play essential roles in hormone and xenobiotic metabolism, oxidative stress, and drug resistance, but little is known about these enzymes in the parasitic nematode Haemonchus contortus. In the present study, 22 AKR genes existing in the H. contortus genome were investigated and a phylogenetic analysis with comparison to AKRs in Caenorhabditis elegans, sheep and humans was conducted. The constitutive transcription levels of all AKRs were measured in eggs, larvae, and adults of H. contortus, and their expression was compared in a drug-sensitive strain (ISE) and a benzimidazole-resistant strain (IRE) previously derived from the sensitive strain by imposing benzimidazole selection pressure. In addition, the inducibility of AKRs by exposure of H. contortus adults to benzimidazole anthelmintic flubendazole in vitro was tested. Phylogenetic analysis demonstrated that the majority of AKR genes in H. contortus lack orthologues in the sheep genome, which is a favorable finding for considering AKRs as potential drug targets. Large differences in the expression levels of individual AKRs were observed, with AKR1, AKR3, AKR8, and AKR10 being the most highly expressed at most developmental stages. Significant changes in the expression of AKRs during the life cycle and pronounced sex differences were found. Comparing the IRE and ISE strains, three AKRs were upregulated, and seven AKRs were downregulated in adults. In addition, the expression of three AKRs was induced by flubendazole exposure in adults of the ISE strain. Based on these results, AKR1, AKR2, AKR3, AKR5, AKR10 and AKR19 in particular merit further investigation and functional characterization with respect to their potential involvement in drug biotransformation and anthelmintic resistance in H. contortus.

Comparative structure activity and target exploration of 1,2-diphenylethynes in Haemonchus contortus and Caenorhabditis elegans

Infections and diseases caused by parasitic nematodes have a major adverse impact on the health and productivity of animals and humans worldwide. The control of these parasites often relies heavily on the treatment with commercially available chemical compounds (anthelmintics). However, the excessive or uncontrolled use of these compounds in livestock animals has led to major challenges linked to drug resistance in nematodes. Therefore, there is a need to develop new anthelmintics with novel mechanism(s) of action. Recently, we identified a small molecule, designated UMW-9729, with nematocidal activity against the free-living model organism Caenorhabditis elegans. Here, we evaluated UMW-9729's potential as an anthelmintic in a structure-activity relationship (SAR) study in C. elegans and the highly pathogenic, blood-feeding Haemonchus contortus (barber's pole worm), and explored the compound-target relationship using thermal proteome profiling (TPP). First, we synthesised and tested 25 analogues of UMW-9729 for their nematocidal activity in both H. contortus (larvae and adults) and C. elegans (young adults), establishing a preliminary nematocidal pharmacophore for both species. We identified several compounds with marked activity against either H. contortus or C. elegans which had greater efficacy than UMW-9729, and found a significant divergence in compound bioactivity between these two nematode species. We also identified a UMW-9729 analogue, designated 25, that moderately inhibited the motility of adult female H. contortus in vitro. Subsequently, we inferred three H. contortus proteins (HCON_00134350, HCON_00021470 and HCON_00099760) and five C. elegans proteins (F30A10.9, F15B9.8, B0361.6, DNC-4 and UNC-11) that interacted directly with UMW-9729; however, no conserved protein target was shared between the two nematode species. Future work aims to extend the SAR investigation in these and other parasitic nematode species, and validate individual proteins identified here as possible targets of UMW-9729. Overall, the present study evaluates this anthelmintic candidate and highlights some challenges associated with early anthelmintic investigation.

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.

Characterisation of protective vaccine antigens from the thiol-containing components of excretory/secretory material of Ostertagia ostertagi

Previous vaccination trials have demonstrated that thiol proteins affinity purified from Ostertagia ostertagi excretory-secretory products (O. ostertagi ES-thiol) are protective against homologous challenge. Here we have shown that protection induced by this vaccine was consistent across four independent vaccine-challenge experiments. Protection is associated with reduced cumulative faecal egg counts across the duration of the trials, relative to control animals. To better understand the diversity of antigens in O. ostertagi ES-thiol we used high-resolution shotgun proteomics to identify 490 unique proteins in the vaccine preparation. The most numerous ES-thiol proteins, with 91 proteins identified, belong to the sperm-coating protein/Tpx/antigen 5/pathogenesis-related protein 1 (SCP/TAPS) family. This family includes previously identified O. ostertagi vaccine antigens O. ostertagi ASP-1 and ASP-2. The ES-thiol fraction also has numerous proteinases, representing three distinct classes, including: metallo-; aspartyl- and cysteine proteinases. In terms of number of family members, the M12 astacin-like metalloproteinases, with 33 proteins, are the most abundant proteinase family in O. ostertagi ES-thiol. The O. ostertagi ES-thiol proteome provides a comprehensive database of proteins present in this vaccine preparation and will guide future vaccine antigen discovery projects.

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.

Next-generation sequencing technologies for helminth diagnostics and surveillance in ruminants: shifting diagnostic barriers

Helminth infections in grazing ruminants are a major issue for livestock farming globally, but are unavoidable in outdoor grazing systems and must be effectively managed to avoid deleterious effects to animal health, and productivity. Next-generation sequencing (NGS) technologies are transforming our understanding of the genetic basis of anthelmintic resistance (AR) and epidemiological studies of ruminant gastrointestinal parasites. They also have the potential to not only help develop and validate molecular diagnostic tests but to be directly used in routine diagnostics integrating species-specific identification and AR into a single test. Here, we review how these developments have opened the pathway for the development of multi-AR and multispecies identification in a single test, with widespread implications for sustainable livestock farming for the future.

Occurrence of gastrointestinal nematodes in lambs in Norway, as assessed by copromicroscopy and droplet digital polymerase chain reaction

BACKGROUND

Gastrointestinal nematodes (GINs) have a major impact on sheep production, health, and welfare worldwide. Norway is no exception, but there are only a few studies on the prevalence of GINs in Norwegian sheep. The aim of this study was to investigate the current occurrence of the most important nematodes in sheep flocks in Norway. Faecal samples were collected from flocks in 2021/2022, mainly from three geographical regions in Norway, i.e., northern, eastern, and western. In each of 134 flocks included, individual samples from 10 lambs (autumn) were pooled. Third stage larvae (L3) were cultivated and harvested (Baermann method) from the pooled samples. The DNA was then extracted and further analysed using droplet digital PCR (ddPCR). This enables assessment of the proportions of the three most important nematode species/genera, i.e., H. contortus, T. circumcincta, and Trichostrongylus. The fractional abundance/relative proportion of each species/genus was assessed by performing duplex assays with universal strongyle and species/genus-specific primers and probe sets. In addition, the occurrence of Nematodirus eggs was assessed by standard faecal egg counts (i.e., McMaster method).

RESULTS

Of the 134 flocks sampled, 24 were from the northern region, 31 from eastern, and 71 from western Norway. In addition, some flocks from central (n = 7), and southern (n = 1) Norway were included. Among the sampled flocks, T. circumcincta occurred most commonly (94%), followed by H. contortus (60%) and Trichostrongylus (55%), and Nematodirus (51%). In general, mixed infections were observed, with 38% and 18% of flocks infected with three or all four genera, respectively.

CONCLUSIONS

The results of this study indicate that GINs are widespread in Norway. Teladorsagia circumcincta seems to be present in most flocks based on this screening. Moreover, the results show that Nematodirus spp. infect lambs throughout the country, predominantly N. battus, and indicate that this nematode has become more abundant, which could lead to an increase in nematodirosis.

SMC-mediated dosage compensation in C. elegans evolved in the presence of an ancestral nematode mechanism

Mechanisms of X chromosome dosage compensation have been studied extensively in a few model species representing clades of shared sex chromosome ancestry. However, the diversity within each clade as a function of sex chromosome evolution is largely unknown. Here, we anchor ourselves to the nematode Caenorhabditis elegans, for which a well-studied mechanism of dosage compensation occurs through a specialized structural maintenance of chromosomes (SMC) complex, and explore the diversity of dosage compensation in the surrounding phylogeny of nematodes. Through phylogenetic analysis of the C. elegans dosage compensation complex and a survey of its epigenetic signatures, including X-specific topologically associating domains (TADs) and X-enrichment of H4K20me1, we found that the condensin-mediated mechanism evolved recently in the lineage leading to Caenorhabditis through an SMC-4 duplication. Intriguingly, an independent duplication of SMC-4 and the presence of X-specific TADs in Pristionchus pacificus suggest that condensin-mediated dosage compensation arose more than once. mRNA-seq analyses of gene expression in several nematode species indicate that dosage compensation itself is ancestral, as expected from the ancient XO sex determination system. Indicative of the ancestral mechanism, H4K20me1 is enriched on the X chromosomes in Oscheius tipulae, which does not contain X-specific TADs or SMC-4 paralogs. Together, our results indicate that the dosage compensation system in C. elegans is surprisingly new, and condensin may have been co-opted repeatedly in nematodes, suggesting that the process of evolving a chromosome-wide gene regulatory mechanism for dosage compensation is constrained.

Significance statement

X chromosome dosage compensation mechanisms evolved in response to Y chromosome degeneration during sex chromosome evolution. However, establishment of dosage compensation is not an endpoint. As sex chromosomes change, dosage compensation strategies may have also changed. In this study, we performed phylogenetic and epigenomic analyses surrounding Caenorhabditis elegans and found that the condensin-mediated dosage compensation mechanism in C. elegans is surprisingly new, and has evolved in the presence of an ancestral mechanism. Intriguingly, condensin-based dosage compensation may have evolved more than once in the nematode lineage, the other time in Pristionchus. Together, our work highlights a previously unappreciated diversity of dosage compensation mechanisms within a clade, and suggests constraints in evolving new mechanisms in the presence of an existing one.

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

Albendazole (a benzimidazole) and ivermectin (a macrocyclic lactone) 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 channels (GluCls), but it is unknown whether GluCl 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 the loss of two GluCl genes (avr-14 and avr-15). The fecundity assays revealed that loss of ben-1 did not provide any fitness benefit in albendazole conditions 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.

Structure-activity relationship and target investigation of 2-aryl quinolines with nematocidal activity

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

A phenotypic screen of the Global Health Priority Box identifies an insecticide with anthelmintic activity

BACKGROUND

Infection with parasitic nematodes (helminths), particularly those of the order Strongylida (such as Haemonchus contortus), can cause significant and burdensome diseases in humans and animals. Widespread drug (anthelmintic) resistance in livestock parasites, the absence of vaccines against most of these nematodes, and a lack of new and effective chemical entities on the commercial market demands the discovery of new anthelmintics. In the present study, we searched the Global Health Priority Box (Medicines for Malaria Venture) for new candidates for anthelmintic development.

METHODS

We employed a whole-organism, motility-based phenotypic screening assay to identify compounds from the Global Health Priority Box with activity against larvae of the model parasite H. contortus, and the free-living comparator nematode Caenorhabditis elegans. Hit compounds were further validated via dose-response assays, with lead candidates then assessed for nematocidal activity against H. contortus adult worms, and additionally, for cytotoxic and mitotoxic effects on human hepatoma (HepG2) cells.

RESULTS

The primary screen against H. contortus and C. elegans revealed or reidentified 16 hit compounds; further validation established MMV1794206, otherwise known as 'flufenerim', as a significant inhibitor of H. contortus larval motility (half-maximal inhibitory concentration [IC50] = 18 μM) and development (IC50 = 1.2 μM), H. contortus adult female motility (100% after 12 h of incubation) and C. elegans larval motility (IC50 = 0.22 μM). Further testing on a mammalian cell line (human hepatoma HepG2 cells), however, identified flufenerim to be both cytotoxic (half-maximal cytotoxic concentration [CC50] < 0.7 μM) and mitotoxic (half-maximal mitotoxic concentration [MC50] < 0.7 μM).

CONCLUSIONS

The in vitro efficacy of MMV1794206 against the most pathogenic stages of H. contortus, as well as the free-living C. elegans, suggests the potential for development as a broad-spectrum anthelmintic compound; however, the high toxicity towards mammalian cells presents a significant hindrance. Further work should seek to establish the protein-drug interactions of MMV1794206 in a nematode model, to unravel the mechanism of action, in addition to an advanced structure-activity relationship investigation to optimise anthelmintic activity and eliminate mammalian cell toxicity.

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.

Multi-omics data elucidate parasite-host-microbiota interactions and resistance to Haemonchus contortus in sheep

BACKGROUND

The integration of molecular data from hosts, parasites, and microbiota can enhance our understanding of the complex biological interactions underlying the resistance of hosts to parasites. Haemonchus contortus, the predominant sheep gastrointestinal parasite species in the tropics, causes significant production and economic losses, which are further compounded by the diminishing efficiency of chemical control owing to anthelmintic resistance. Knowledge of how the host responds to infection and how the parasite, in combination with microbiota, modulates host immunity can guide selection decisions to breed animals with improved parasite resistance. This understanding will help refine management practices and advance the development of new therapeutics for long-term helminth control.

METHODS

Eggs per gram (EPG) of feces were obtained from Morada Nova sheep subjected to two artificial infections with H. contortus and used as a proxy to select animals with high resistance or susceptibility for transcriptome sequencing (RNA-seq) of the abomasum and 50 K single-nucleotide genotyping. Additionally, RNA-seq data for H. contortus were generated, and amplicon sequence variants (ASV) were obtained using polymerase chain reaction amplification and sequencing of bacterial and archaeal 16S ribosomal RNA genes from sheep feces and rumen content.

RESULTS

The heritability estimate for EPG was 0.12. GAST, GNLY, IL13, MGRN1, FGF14, and RORC genes and transcripts were differentially expressed between resistant and susceptible animals. A genome-wide association study identified regions on chromosomes 2 and 11 that harbor candidate genes for resistance, immune response, body weight, and adaptation. Trans-expression quantitative trait loci were found between significant variants and differentially expressed transcripts. Functional co-expression modules based on sheep genes and ASVs correlated with resistance to H. contortus, showing enrichment in pathways of response to bacteria, immune and inflammatory responses, and hub features of the Christensenellaceae, Bacteroides, and Methanobrevibacter genera; Prevotellaceae family; and Verrucomicrobiota phylum. In H. contortus, some mitochondrial, collagen-, and cuticle-related genes were expressed only in parasites isolated from susceptible sheep.

CONCLUSIONS

The present study identified chromosome regions, genes, transcripts, and pathways involved in the elaborate interactions between the sheep host, its gastrointestinal microbiota, and the H. contortus parasite. These findings will assist in the development of animal selection strategies for parasite resistance and interdisciplinary approaches to control H. contortus infection in sheep.

Getting around the roundworms: Identifying knowledge gaps and research priorities for the ascarids

The ascarids are a large group of parasitic nematodes that infect a wide range of animal species. In humans, they cause neglected diseases of poverty; many animal parasites also cause zoonotic infections in people. Control measures include hygiene and anthelmintic treatments, but they are not always appropriate or effective and this creates a continuing need to search for better ways to reduce the human, welfare and economic costs of these infections. To this end, Le Studium Institute of Advanced Studies organized a two-day conference to identify major gaps in our understanding of ascarid parasites with a view to setting research priorities that would allow for improved control. The participants identified several key areas for future focus, comprising of advances in genomic analysis and the use of model organisms, especially Caenorhabditis elegans, a more thorough appreciation of the complexity of host-parasite (and parasite-parasite) communications, a search for novel anthelmintic drugs and the development of effective vaccines. The participants agreed to try and maintain informal links in the future that could form the basis for collaborative projects, and to co-operate to organize future meetings and workshops to promote ascarid research.

Analysis of P-gp genes relative expression associated to ivermectin resistance in Haemonchus contortus larval stages from in vitro cultures (L3 and xL3) and from gerbils (Meriones unguiculatus) (L4) as models of study

The aim of the study was to compare the relative gene expression of Haemonchus contortus P-glycoprotein genes (Hco-pgp) between fourth (L4), infective (L3), and transitory infective (xL3) larval stages as laboratory models to study ivermectin (IVM) resistance. The H. contortus resistant to IVM (IVMr) and susceptible to IVM (IVMs) strains were used to develop xL3in vitro culture and to infect Meriones unguiculatus (gerbils) to collect L4 stages. Morphometric differences were evaluated from 25 individuals of H. contortus from each strain. Relative gene expression from xL3 and L4 was determined between comparison of IVMr stages and from IVMr vs IVMs stages. Seven Hco-pgp genes (1, 2, 3, 4, 9, 10, and 16) were analysed by RT-qPCR using L3 stage as control group, per strain, and GAPDH and β-tubulin as constitutive genes. Morphological changes were confirmed between xL3 and L4 developing oral shape, oesophagus, and intestinal tube. In addition, the body length and width showed statistical differences (p < 0.05). The Hco-pgp1, 2, 3, and 4 genes (p < 0.05) were upregulated from 7.1- to 463.82-fold changes between IVMr stages, and Hco-pgp9 (13.12-fold) and Hco-pgp10 (13.56-fold) genes showed differences between L4 and xL3, respectively. The comparative study between IVMr vs IVMs strains associated to xL3 and L4 displayed significant upregulation for most of the Hco-pgp genes among 4.89-188.71 fold-change. In conclusion, these results suggest the use of H. contortus xL3 and L4 as suitable laboratory models to study IVMr associated with Hco-pgp genes to contribute to the understanding of anthelmintic resistance.

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.

Real-time single-base specific detection of the Haemonchus contortus S168T variant associated with levamisole resistance using loop-primer endonuclease cleavage loop-mediated isothermal amplification

Haemonchus contortus is a parasitic haematophagous nematode that primarily affects small ruminants and causes significant economic loss to the global livestock industry. Treatment of haemonchosis typically relies on broad-spectrum anthelmintics, resistance to which is an important cause of treatment failure. Resistance to levamisole remains less widespread than to other major anthelmintic classes, prompting the need for more effective and accurate surveillance to maintain its efficacy. Loop-primer endonuclease cleavage loop-mediated isothermal amplification (LEC-LAMP) is a recently developed diagnostic method that facilitates multiplex target detection with single nucleotide polymorphism (SNP) specificity and portable onsite testing. In this study, we designed a new LEC-LAMP assay and applied it to detect the levamisole resistance marker S168T in H. contortus. We explored multiplexing probes for both the resistant S168T and the susceptible S168 alleles in a single-tube assay. We then included a generic probe to detect the acr-8 gene in the multiplex assay, which could facilitate the quantification of both resistance markers and overall genetic material from H. contortus in a single step. Our results showed promising application of these technologies, demonstrating a proof-of-concept assay which is amenable to detection of resistance alleles within the parasite population, with the potential for multiplex detection, and point-of-care application enabled by lateral flow end-point detection. However, further optimisation and validation is necessary.

A haplotype-like, chromosome-level assembled and annotated genome of Biomphalaria glabrata, an important intermediate host of schistosomiasis and the best studied model of schistosomiasis vector snails

Schistosomiasis is one of the world's most devastating parasitic diseases, afflicting 251 million people globally. The Neotropical snail Biomphalaria glabrata is an important intermediate host of the human blood fluke Schistosoma mansoni and a predominant model for schistosomiasis research. To fully exploit this model snail for biomedical research, here we report a haplotype-like, chromosome-level assembled and annotated genome of the homozygous iM line of B. glabrata that we developed at the University of New Mexico. Using multiple sequencing platforms, including Illumina, PacBio, and Omni-C sequencing, 18 sequence contact matrices representing 18 haploid chromosomes (2n = 36) were generated (337x genome coverage), and 96.5% of the scaffold sequences were anchored to the 18 chromosomes. Protein-coding genes (n = 34,559), non-coding RNAs (n = 2,406), and repetitive elements (42.52% of the genome) were predicted for the whole genome, and detailed annotations for individual chromosomes were also provided. Using this genomic resource, we have investigated the genomic structure and organization of the Toll-like receptor (TLR) and fibrinogen-domain containing protein (FReD) genes, the two important immune-related gene families. Notably, TLR-like genes are scattered on 13 chromosomes. In contrast, almost all (39 of 40) fibrinogen-related genes (FREPs) (immunoglobulin superfamily (IgSF) + fibrinogen (FBG)) are clustered within a 5-million nucleotide region on chromosome 13, yielding insight into mechanisms involved in the diversification of FREPs. This is the first genome of schistosomiasis vector snails that has been assembled at the chromosome level, annotated, and analyzed. It serves as a valuable resource for a deeper understanding of the biology of vector snails, especially Biomphalaria snails.

Structure activity relationship and target prediction for ABX464 analogues in Caenorhabditis elegans

Global challenges with treatment failures and/or widespread resistance in parasitic worms against commercially available anthelmintics lend impetus to the development of new anthelmintics with novel mechanism(s) of action. The free-living nematode Caenorhabditis elegans is an important model organism used for drug discovery, including the screening and structure-activity investigation of new compounds, and target deconvolution. Previously, we conducted a whole-organism phenotypic screen of the 'Pandemic Response Box' (from Medicines for Malaria Venture, MMV) and identified a hit compound, called ABX464, with activity against C. elegans and a related, parasitic nematode, Haemonchus contortus. Here, we tested a series of 44 synthesized analogues to explore the pharmacophore of activity on C. elegans and revealed five compounds whose potency was similar or greater than that of ABX464, but which were not toxic to human hepatoma (HepG2) cells. Subsequently, we employed thermal proteome profiling (TPP), protein structure prediction and an in silico-docking algorithm to predict ABX464-target candidates. Taken together, the findings from this study contribute significantly to the early-stage drug discovery of a new nematocide based on ABX464. Future work is aimed at validating the ABX464-protein interactions identified here, and at assessing ABX464 and associated analogues against a panel of parasitic nematodes, towards developing a new anthelmintic with a mechanism of action that is distinct from any of the compounds currently-available commercially.

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 mixed amplicon metabarcoding and sequencing approach for surveillance of drug resistance to levamisole and benzimidazole in Haemonchus spp

Anthelmintic-resistant parasitic nematodes present a significant threat to sustainable livestock production worldwide. The ability to detect the emergence of anthelmintic resistance at an early stage, and therefore determine which drugs remain most effective, is crucial for minimising production losses. Despite many years of research into the molecular basis of anthelmintic resistance, no molecular-based tools are commercially available for the diagnosis of resistance as it emerges in field settings. We describe a mixed deep amplicon sequencing approach to determine the frequency of the levamisole (LEV)-resistant single nucleotide polymorphism (SNP) within arc-8 exon 4 (S168T) in Haemonchus spp., coupled with benzimidazole (BZ)-resistant SNPs within β-tubulin isotype-1 and the internal transcribed spacer-2 (ITS-2) nemabiome. This constitutes the first known multi-drug and multi-species molecular diagnostic developed for helminths of veterinary importance. Of the ovine, bovine, caprine and camelid Australian field isolates we tested, S168T was detected in the majority of Haemonchus spp. populations from sheep and goats, but rarely at a frequency greater than 16%; an arbitrary threshold we set based on whole genome sequencing (WGS) of LEV-resistant Haemonchus contortus GWBII. Overall, BZ resistance was far more prevalent in Haemonchus spp. than LEV resistance, confirming that LEV is still an effective anthelmintic class for small ruminants in New South Wales, Australia. The mixed amplicon metabarcoding approach described herein paves the way towards the use of large scale sequencing as a surveillance technology in the field, the results of which can be translated into evidence-based recommendations for the livestock sector.

Comparative Genomics of Sex, Chromosomes, and Sex Chromosomes in Caenorhabditis elegans and Other Nematodes

The nematode phylum has evolved a remarkable diversity of reproductive modes, including the repeated emergence of asexuality and hermaphroditism across divergent clades. The species-richness and small genome size of nematodes make them ideal systems for investigating the genome-wide causes and consequences of such major transitions. The availability of functional annotations for most Caenorhabditis elegans genes further allows the linking of patterns of gene content evolution with biological processes. Such gene-centric studies were recently complemented by investigations of chromosome evolution that made use of the first chromosome-scale genome assemblies outside the Caenorhabditis genus. This review highlights recent comparative genomic studies of reproductive mode evolution addressing the hybrid origin of asexuality and the parallel gene loss following the emergence of hermaphroditism. It further summarizes ongoing efforts to characterize ancient linkage blocks called Nigon elements, which form central units of chromosome evolution. Fusions between Nigon elements have been demonstrated to impact recombination and speciation. Finally, multiple recent fusions between autosomal and the sex-linked Nigon element reveal insights into the dynamic evolution of sex chromosomes across various timescales.

Ancient diversity in host-parasite interaction genes in a model parasitic nematode

Host-parasite interactions exert strong selection pressures on the genomes of both host and parasite. These interactions can lead to negative frequency-dependent selection, a form of balancing selection that is hypothesised to explain the high levels of polymorphism seen in many host immune and parasite antigen loci. Here, we sequence the genomes of several individuals of Heligmosomoides bakeri, a model parasite of house mice, and Heligmosomoides polygyrus, a closely related parasite of wood mice. Although H. bakeri is commonly referred to as H. polygyrus in the literature, their genomes show levels of divergence that are consistent with at least a million years of independent evolution. The genomes of both species contain hyper-divergent haplotypes that are enriched for proteins that interact with the host immune response. Many of these haplotypes originated prior to the divergence between H. bakeri and H. polygyrus, suggesting that they have been maintained by long-term balancing selection. Together, our results suggest that the selection pressures exerted by the host immune response have played a key role in shaping patterns of genetic diversity in the genomes of 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 (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.

Benzimidazole-resistance associated mutation in Haemonchus contortus in Norwegian sheep, as detected by droplet digital PCR

The aim of this study was to investigate the occurrence of benzimidazole-resistant Haemonchus contortus in Norwegian sheep flocks. Screening was based on detection of one of the resistance-conferring mutations in the β tubulin isotype 1 gene (F200Y, TAC) in larvae (L3) cultivated from H. contortus eggs from naturally infected sheep. Faecal samples were collected in 2021/2022 from flocks in the northern (n = 34), central (n = 7), eastern (n = 40), southern (n = 1), and western (n = 87) areas of Norway. In total, samples were taken from 169 flocks (spring-ewes samples: 167, autumn-lambs samples: 134). Individual faecal samples were collected from 10 randomly selected ewes (spring) and 10 randomly selected lambs (autumn) in each flock. Faecal samples collected from each flock on each occasion were pooled (lamb and ewe samples pooled separately) and cultured for L3 development. After harvest of larvae (Baermann method), DNA was extracted and then analysed using droplet digital PCR with primer/probe sets targeting the BZ-associated F200Y (TAC) mutation. Haemonchus was found in 60% (80/134) of samples from lambs, and in 63% (106/167) from ewes. Among these, the F200Y mutation was detected in 73% (58/80) of larval samples from lambs and 69% (73/106) of larval samples from ewes, respectively. Although regional differences were evident, the mutation was detected in all areas indicating a widespread distribution and a strong potential for an increasing problem with treatment-resistant haemonchosis in Norwegian sheep flocks.

Novel and improved Caenorhabditis briggsae gene models generated by community curation

BACKGROUND

The nematode Caenorhabditis briggsae has been used as a model in comparative genomics studies with Caenorhabditis elegans because of their striking morphological and behavioral similarities. However, the potential of C. briggsae for comparative studies is limited by the quality of its genome resources. The genome resources for the C. briggsae laboratory strain AF16 have not been developed to the same extent as C. elegans. The recent publication of a new chromosome-level reference genome for QX1410, a C. briggsae wild strain closely related to AF16, has provided the first step to bridge the gap between C. elegans and C. briggsae genome resources. Currently, the QX1410 gene models consist of software-derived gene predictions that contain numerous errors in their structure and coding sequences. In this study, a team of researchers manually inspected over 21,000 gene models and underlying transcriptomic data to repair software-derived errors.

RESULTS

We designed a detailed workflow to train a team of nine students to manually curate gene models using RNA read alignments. We manually inspected the gene models, proposed corrections to the coding sequences of over 8,000 genes, and modeled thousands of putative isoforms and untranslated regions. We exploited the conservation of protein sequence length between C. briggsae and C. elegans to quantify the improvement in protein-coding gene model quality and showed that manual curation led to substantial improvements in the protein sequence length accuracy of QX1410 genes. Additionally, collinear alignment analysis between the QX1410 and AF16 genomes revealed over 1,800 genes affected by spurious duplications and inversions in the AF16 genome that are now resolved in the QX1410 genome.

CONCLUSIONS

Community-based, manual curation using transcriptome data is an effective approach to improve the quality of software-derived protein-coding genes. The detailed protocols provided in this work can be useful for future large-scale manual curation projects in other species. Our manual curation efforts have brought the QX1410 gene models to a comparable level of quality as the extensively curated AF16 gene models. The improved genome resources for C. briggsae provide reliable tools for the study of Caenorhabditis biology and other related nematodes.

A member of the CAP protein superfamily, Hc-CAP-15, is important for the parasitic-stage development of Haemonchus contortus

BACKGROUND

The CAP superfamily proteins are distributed widely in eukaryotes and play crucial roles in various biological processes. However, very little is known about their functions in parasitic nematodes, including Haemonchus contortus, a socioeconomically important parasitic nematode. We have therefore studied a member of the CAP protein family of H. contortus, named Hc-CAP-15, with the aim to explore its roles in regulating the parasitic developmental process.

METHODS

The conservation and phylogenetic relationships, spatial expression and temporal transcription profiles of Hc-CAP/cap-15, as well its biological function during parasite development were investigated using bioinformatics, immunofluorescence, real-time PCR and RNA interference (RNAi).

RESULTS

Hc-CAP-15 was found to be a single-domain CAP protein consisting of four conserved motifs that is localized in the cuticle, intestine and oocyte of adult worms. Hc-cap-15 was transcribed at all developmental stages of H. contortus, with the highest transcription level in parasitic fourth-stage larvae (L4s). Silencing of Hc-cap-15 resulted in a significant increase in the body length of L4s.

CONCLUSIONS

The results suggested that Hc-CAP-15 is important for the development of H. contortus. Our findings provide a basis for further study of the functions of the CAP family proteins in H. contortus and related parasitic nematodes.

Genome-Wide Analysis of Haemonchus contortus Proteases and Protease Inhibitors Using Advanced Informatics Provides Insights into Parasite Biology and Host-Parasite Interactions

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

New technologies to study helminth development and host-parasite interactions

How parasites develop and survive, and how they stimulate or modulate host immune responses are important in understanding disease pathology and for the design of new control strategies. Microarray analysis and bulk RNA sequencing have provided a wealth of data on gene expression as parasites develop through different life-cycle stages and on host cell responses to infection. These techniques have enabled gene expression in the whole organism or host tissue to be detailed, but do not take account of the heterogeneity between cells of different types or developmental stages, nor the spatial organisation of these cells. Single-cell RNA-seq (scRNA-seq) adds a new dimension to studying parasite biology and host immunity by enabling gene profiling at the individual cell level. Here we review the application of scRNA-seq to establish gene expression cell atlases for multicellular helminths and to explore the expansion and molecular profile of individual host cell types involved in parasite immunity and tissue repair. Studying host-parasite interactions in vivo is challenging and we conclude this review by briefly discussing the applications of organoids (stem-cell derived mini-tissues) to examine host-parasite interactions at the local level, and as a potential system to study parasite development in vitro. Organoid technology and its applications have developed rapidly, and the elegant studies performed to date support the use of organoids as an alternative in vitro system for research on helminth parasites.

The Proteome and Lipidome of Extracellular Vesicles from Haemonchus contortus to Underpin Explorations of Host-Parasite Cross-Talk

Many parasitic worms have a major adverse impact on human and animal populations worldwide due to the chronicity of their infections. There is a growing body of evidence indicating that extracellular vesicles (EVs) are intimately involved in modulating (suppressing) inflammatory/immune host responses and parasitism. As one of the most pathogenic nematodes of livestock animals, Haemonchus contortus is an ideal model system for EV exploration. Here, employing a multi-step enrichment process (in vitro culture, followed by ultracentrifugation, size exclusion and filtration), we enriched EVs from H. contortus and undertook the first comprehensive (qualitative and quantitative) multi-omic investigation of EV proteins and lipids using advanced liquid chromatography-mass spectrometry and informatics methods. We identified and quantified 561 proteins and 446 lipids in EVs and compared these molecules with those of adult worms. We identified unique molecules in EVs, such as proteins linked to lipid transportation and lipid species (i.e., sphingolipids) associated with signalling, indicating the involvement of these molecules in parasite-host cross-talk. This work provides a solid starting point to explore the functional roles of EV-specific proteins and lipids in modulating parasite-host cross-talk, and the prospect of finding ways of disrupting or interrupting this relationship to suppress or eliminate parasite infection.

Novel and improved Caenorhabditis briggsae gene models generated by community curation

Background

The nematode Caenorhabditis briggsae has been used as a model for genomics studies compared to Caenorhabditis elegans because of its striking morphological and behavioral similarities. These studies yielded numerous findings that have expanded our understanding of nematode development and evolution. However, the potential of C. briggsae to study nematode biology is limited by the quality of its genome resources. The reference genome and gene models for the C. briggsae laboratory strain AF16 have not been developed to the same extent as C. elegans . The recent publication of a new chromosome-level reference genome for QX1410, a C. briggsae wild strain closely related to AF16, has provided the first step to bridge the gap between C. elegans and C. briggsae genome resources. Currently, the QX1410 gene models consist of protein-coding gene predictions generated from short- and long-read transcriptomic data. Because of the limitations of gene prediction software, the existing gene models for QX1410 contain numerous errors in their structure and coding sequences. In this study, a team of researchers manually inspected over 21,000 software-derived gene models and underlying transcriptomic data to improve the protein-coding gene models of the C. briggsae QX1410 genome.

Results

We designed a detailed workflow to train a team of nine students to manually curate genes using RNA read alignments and predicted gene models. We manually inspected the gene models using the genome annotation editor, Apollo, and proposed corrections to the coding sequences of over 8,000 genes. Additionally, we modeled thousands of putative isoforms and untranslated regions. We exploited the conservation of protein sequence length between C. briggsae and C. elegans to quantify the improvement in protein-coding gene model quality before and after curation. Manual curation led to a substantial improvement in the protein sequence length accuracy of QX1410 genes. We also compared the curated QX1410 gene models against the existing AF16 gene models. The manual curation efforts yielded QX1410 gene models that are similar in quality to the extensively curated AF16 gene models in terms of protein-length accuracy and biological completeness scores. Collinear alignment analysis between the QX1410 and AF16 genomes revealed over 1,800 genes affected by spurious duplications and inversions in the AF16 genome that are now resolved in the QX1410 genome.

Conclusions

Community-based, manual curation using transcriptome data is an effective approach to improve the quality of software-derived protein-coding genes. Comparative genomic analysis using a related species with high-quality reference genome(s) and gene models can be used to quantify improvements in gene model quality in a newly sequenced genome. The detailed protocols provided in this work can be useful for future large-scale manual curation projects in other species. The chromosome-level reference genome for the C. briggsae strain QX1410 far surpasses the quality of the genome of the laboratory strain AF16, and our manual curation efforts have brought the QX1410 gene models to a comparable level of quality to the previous reference, AF16. The improved genome resources for C. briggsae provide reliable tools for the study of Caenorhabditis biology and other related nematodes.

Immunization with thiol-binding proteins from Haemonchus contortus adult worms partially protects goats against infection during prepatency

To contribute of the knowledge of the immune mechanisms underlying the response to the immunization of goats with thiol-binding proteins fractions (PBS-TSBP) from Haemonchus contortus (H. contortus) adult worms, this study analyzed the degree of protection and the immune responses developed against the parasite after vaccination with this antigenic complex during the time-elapsing between challenge with L3 of the parasite and the development of adult worms, evidenced by the appearance of first faecal eggs (prepatent period or prepatency). Goat kids immunized with PBS-TBSP generated an immune response during the prepatency which translates into a reduction in the number of worms, as well as a lower reduction on packed cell volume and plasma protein levels in relation to the non-vaccinated animals. As previously described in other studies carried out after the prepatent period, this protection was associated with a systemic humoral response. At the local level, a specific humoral response was also observed, together with an immune-inflammatory infiltrate in the gastric mucosa of MCH-II + cells and CD4⁺ lymphocytes, whose number was associated with a reduction in the number of worms and an increase in plasma proteins. A high peripheral eosinophilia was detected, but no corresponding increased infiltration of the gastric mucosa by eosinophils or globular leukocytes was observed. In agreement with previous data on the immunolocalization of the antigens used here, the results obtained contribute to the idea that these may be excretion/secretion (E/S) products necessary for parasite survival, whose inactivation during the larval and/or pre-adult stages may have contributed to immunoprotection.

Assembly and Analysis of Haemonchus contortus Transcriptome as a Tool for the Knowledge of Ivermectin Resistance Mechanisms

Haemonchus contortus (Hc) is an important parasitic nematode of small ruminants. In this study we assembled the transcriptome of Hc as a model to contribute to the knowledge about the profile of the differential gene expression between two Mexican Hc strains under different anthelmintic resistance statuses, one susceptible and the other resistant to ivermectin (IVMs and IVMr, respectively), in order to improve and/or to have new strategies of control and diagnosis. The transcript sequence reads were assembled and annotated. Overall, ~127 Mbp were assembled and distributed into 77,422 transcript sequences, and 4394 transcripts of the de novo transcriptome were matched base on at least one of the following criteria: (1) Phylum Nemathelminthes and Platyhelminthes, important for animal health care, and (2) ≥55% of sequence identity with other organisms. The gene ontology (GO) enrichment analysis (GOEA) was performed to study the level of gene regulation to IVMr and IVMs strains using Log Fold Change (LFC) filtering values ≥ 1 and ≥ 2. The upregulated-displayed genes obtained via GOEA were: 1993 (for LFC ≥ 1) and 1241 (for LFC ≥ 2) in IVMr and 1929 (for LFC ≥ 1) and 835 (for LFC ≥ 2) in IVMs. The enriched GO terms upregulated per category identified the intracellular structure, intracellular membrane-bounded organelle and integral component of the cell membrane as some principal cellular components. Meanwhile, efflux transmembrane transporter activity, ABC-type xenobiotic transporter activity and ATPase-coupled transmembrane transporter activity were associated with molecular function. Responses to nematicide activity, pharyngeal pumping and positive regulation of synaptic assembly were classified as biological processes that might be involved in events related to the anthelmintic resistance (AR) and nematode biology. The filtering analysis of both LFC values showed similar genes related to AR. This study deepens our knowledge about the mechanisms behind the processes of H. contortus in order to help in tool production and to facilitate the reduction of AR and promote the development of other control strategies, such as anthelmintic drug targets and vaccines.

Population genomics of helminth parasites

Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.

Chromosome fusions repatterned recombination rate and facilitated reproductive isolation during Pristionchus nematode speciation

Large-scale genome-structural evolution is common in various organisms. Recent developments in speciation genomics revealed the importance of inversions, whereas the role of other genome-structural rearrangements, including chromosome fusions, have not been well characterized. We study genomic divergence and reproductive isolation of closely related nematodes: the androdioecious (hermaphroditic) model Pristionchus pacificus and its dioecious sister species Pristionchus exspectatus. A chromosome-level genome assembly of P. exspectatus using single-molecule and Hi-C sequencing revealed a chromosome-wide rearrangement relative to P. pacificus. Strikingly, genomic characterization and cytogenetic studies including outgroup species Pristionchus occultus indicated two independent fusions involving the same chromosome, ChrIR, between these related species. Genetic linkage analysis indicated that these fusions altered the chromosome-wide pattern of recombination, resulting in large low-recombination regions that probably facilitated the coevolution between some of the ~14.8% of genes across the entire genomes. Quantitative trait locus analyses for hybrid sterility in all three sexes revealed that major quantitative trait loci mapped to the fused chromosome ChrIR. While abnormal chromosome segregations of the fused chromosome partially explain hybrid female sterility, hybrid-specific recombination that breaks linkage of genes in the low-recombination region was associated with hybrid male sterility. Thus, recent chromosome fusions repatterned recombination rate and drove reproductive isolation during Pristionchus speciation.

An informatic workflow for the enhanced annotation of excretory/secretory proteins of Haemonchus contortus

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

Chromosome-length genome assembly of Teladorsagia circumcincta - a globally important helminth parasite in livestock

BACKGROUND

Gastrointestinal (GIT) helminthiasis is a global problem that affects livestock health, especially in small ruminants. One of the major helminth parasites of sheep and goats, Teladorsagia circumcincta, infects the abomasum and causes production losses, reductions in weight gain, diarrhoea and, in some cases, death in young animals. Control strategies have relied heavily on the use of anthelmintic medication but, unfortunately, T. circumcincta has developed resistance, as have many helminths. Vaccination offers a sustainable and practical solution, but there is no commercially available vaccine to prevent Teladorsagiosis. The discovery of new strategies for controlling T. circumcincta, such as novel vaccine targets and drug candidates, would be greatly accelerated by the availability of better quality, chromosome-length, genome assembly because it would allow the identification of key genetic determinants of the pathophysiology of infection and host-parasite interaction. The available draft genome assembly of T. circumcincta (GCA_002352805.1) is highly fragmented and thus impedes large-scale investigations of population and functional genomics.

RESULTS

We have constructed a high-quality reference genome, with chromosome-length scaffolds, by purging alternative haplotypes from the existing draft genome assembly and scaffolding the result using chromosome conformation, capture-based, in situ Hi-C technique. The improved (Hi-C) assembly resulted in six chromosome-length scaffolds with length ranging from 66.6 Mbp to 49.6 Mbp, 35% fewer sequences and reduction in size. Substantial improvements were also achieved in both the values for N50 (57.1 Mbp) and L50 (5 Mbp). A higher and comparable level of genome and proteome completeness was achieved for Hi-C assembly on BUSCO parameters. The Hi-C assembly had a greater synteny and number of orthologs with a closely related nematode, Haemonchus contortus.

CONCLUSION

This improved genomic resource is suitable as a foundation for the identification of potential targets for vaccine and drug development.

The presence and relative frequency detection of the levamisole-resistance-associated S168T substitution in hco-acr-8 in Haemonchus contortus

Parasitic sheep nematodes, among which Haemonchus contortus is often considered to be the most clinically important, exact a significant toll on the animals, not least because of their capacity to evolve drug resistance. Despite decades of research, our understanding of the mechanism of resistance to compounds such as levamisole is fairly limited, which therefore constrains our ability to develop sensitive and efficient molecular diagnostic tools for rapid and accurate resistance detection in field settings. Herein, we investigated the presence and frequency of the newly reported, levamisole-resistance-associated, mutation, yielding a S168T substitution in exon 4 of hco-acr-8, in six different phenotypically described isolates (three susceptible and three resistant), three Swedish field isolates and eight larvae culture samples, the latter two of which originated on farms where levamisole showed complete parasite elimination. For this purpose, we created both an allele-specific and droplet digital PCR approaches and found the mutated allele to be present only in the Kokstad isolate, whereas the other five as well as both the Swedish isolates and larvae cultures displayed only the non-mutated, serine-encoding, allele. While the finding of only the non-mutated allele in the phenotypically susceptible and Swedish isolate and larvae culture samples seemed sensible, we speculate that for the other two phenotypically resistant isolates, different (perhaps secondary) variants are responsible for conferring the resistance to levamisole phenotype, given the polygenic nature of levamisole resistance. All in all, despite the limited number of samples tested here, the mutation causing the S168T substitution in hco-acr-8 represents a plausible levamisole resistance-associated variant in, at least, some isolates of H. contortus.

Transcriptomics of ivermectin response in Caenorhabditis elegans: Integrating abamectin quantitative trait loci and comparison to the Ivermectin-exposed DA1316 strain

Parasitic nematodes pose a significant threat to human and animal health, as well as cause economic losses in the agricultural sector. The use of anthelmintic drugs, such as Ivermectin (IVM), to control these parasites has led to widespread drug resistance. Identifying genetic markers of resistance in parasitic nematodes can be challenging, but the free-living nematode Caenorhabditis elegans provides a suitable model. In this study, we aimed to analyze the transcriptomes of adult C. elegans worms of the N2 strain exposed to the anthelmintic drug Ivermectin (IVM), and compare them to those of the resistant strain DA1316 and the recently identified Abamectin Quantitative Trait Loci (QTL) on chromosome V. We exposed pools of 300 adult N2 worms to IVM (10-7 and 10-8 M) for 4 hours at 20°C, extracted total RNA and sequenced it on the Illumina NovaSeq6000 platform. Differentially expressed genes (DEGs) were determined using an in-house pipeline. The DEGs were compared to genes from a previous microarray study on IVM-resistant C. elegans and Abamectin-QTL. Our results revealed 615 DEGs (183 up-regulated and 432 down-regulated genes) from diverse gene families in the N2 C. elegans strain. Of these DEGs, 31 overlapped with genes from IVM-exposed adult worms of the DA1316 strain. We identified 19 genes, including the folate transporter (folt-2) and the transmembrane transporter (T22F3.11), which exhibited an opposite expression in N2 and the DA1316 strain and were deemed potential candidates. Additionally, we compiled a list of potential candidates for further research including T-type calcium channel (cca-1), potassium chloride cotransporter (kcc-2), as well as other genes such as glutamate-gated channel (glc-1) that mapped to the Abamectin-QTL.

Allele specific PCR for a major marker of levamisole resistance in Haemonchus contortus

Haemonchus contortus is a haematophagous parasitic nematode that infects small ruminants and causes significant animal health concerns and economic losses within the livestock industry on a global scale. Treatment primarily depends on broad-spectrum anthelmintics, however, resistance is established or rapidly emerging against all major drug classes. Levamisole (LEV) remains an important treatment option for parasite control, as resistance to LEV is less prevalent than to members of other major classes of anthelmintics. LEV is an acetylcholine receptor (AChR) agonist that, when bound, results in paralysis of the worm. Numerous studies implicated the AChR sub-unit, ACR-8, in LEV sensitivity and in particular, the presence of a truncated acr-8 transcript or a deletion in the acr-8 locus in some resistant isolates. Recently, a single non-synonymous SNP in acr-8 conferring a serine-to-threonine substitution (S168T) was identified that was strongly associated with LEV resistance. Here, we investigate the role of genetic variation at the acr-8 locus in a controlled genetic cross between the LEV susceptible MHco3(ISE) and LEV resistant MHco18(UGA2004) isolates of H. contortus. Using single worm PCR assays, we found that the presence of S168T was strongly associated with LEV resistance in the parental isolates and F3 progeny of the genetic cross surviving LEV treatment. We developed and optimised an allele-specific PCR assay for the detection of S168T and validated the assay using laboratory isolates and field samples that were phenotyped for LEV resistance. In the LEV-resistant field population, a high proportion (>75%) of L3 encoded the S168T variant, whereas the variant was absent in the susceptible isolates studied. These data further support the potential role of acr-8 S168T in LEV resistance, with the allele-specific PCR providing an important step towards establishing a sensitive molecular diagnostic test for LEV resistance.