Genomic landscape of drug response reveals mediators of anthelmintic resistance

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

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

The parasitic nematode Haemonchus contortus lacks molybdenum cofactor synthesis, leading to sulphite sensitivity and lethality in vitro

Sulphite oxidase has an essential role in detoxifying environmental and endogenously generated sulphite into sulphate and requires the molybdenum cofactor (Moco) to function. Until recently it was believed that the synthesis pathway for Moco was so important for survival that it was conserved in all multicellular animals. Here we report the use of comparative genomics to identify the absence of the first enzyme involved in Moco synthesis in Haemonchus contortus, a highly pathogenic and economically important helminth of livestock that, similar to many parasitic nematode species, has proved difficult to maintain in vitro. We show that Moco deficiency in Haemonchus leads to a high sensitivity to environmental sulphite and limits the ability to maintain the early parasitic larval stages in vitro. Analogous losses in Moco synthesis in other recently sequenced nematode species are also identified. These findings may lead to improved culture methods for parasitic nematodes and to novel approaches for their control.

Precise Lineage Tracking Using Molecular Barcodes Demonstrates Fitness Trade-offs for Ivermectin Resistance in Nematodes

A fundamental tenet of evolutionary genetics is that the direction and strength of selection on individual loci varies with the environment. Barcoded evolutionary lineage tracking is a powerful approach for high-throughput measurement of selection within experimental evolution that to date has largely been restricted to studies within microbial systems, largely because the random integration of barcodes within animals is limited by physical and molecular protection of the germline. Here, we use the recently developed TARDIS barcoding system in Caenorhabditis elegans (Stevenson et al., 2023) to implement the first randomly inserted genomic-barcode experimental evolution animal model and use this system to precisely measure the influence of the concentration of the anthelmintic compound ivermectin on the strength of selection on an ivermectin resistance cassette. The combination of the trio of knockouts in neuronally expressed GluCl channels, avr-14, avr-15, and glc-1, has been previously demonstrated to provide resistance to ivermectin at high concentrations. Varying the concentration of ivermectin in liquid culture allows the strength of selection on these genes to be precisely controlled within populations of millions of individuals, yielding the largest animal experimental evolution study to date. The frequency of each barcode was determined at multiple time points via sequencing at deep coverage and then used to estimate the fitness of the individual lineages in the population. The mutations display a high cost to resistance at low concentrations, rapidly losing out to wildtype genotypes, but the balance tips in their favor when the ivermectin concentration exceeds 2nM. This trade-off in resistance is likely generated by a hindered rate of development in resistant individuals. Our results demonstrate that C. elegans can be used to generate high precision estimates of fitness using a high-throughput barcoding approach to yield novel insights into evolutionarily and economically important traits.

P-glycoproteins in anthelmintic safety, efficacy, and resistance

P-glycoprotein (PGP) is a pivotal transmembrane transporter governing the cellular flux of diverse substances shielding mammals from toxics. It can thwart the effectiveness of medicines such as ivermectin (IVM) and other macrocyclic lactone (ML) anthelmintics, undermining therapeutic efforts. We analyze the role of PGPs in limiting the toxicity of these drugs in hosts, and their potential contribution to anthelmintic resistance in nematodes. Targeting nematode PGPs to increase drug sensitivity to MLs seems interesting, but is hampered by the lack of selective inhibitors. The nuclear hormone receptor (NHR)-8 should be seriously considered as a target because it upregulates multiple PGPs involved in anthelmintic resistance and it is specific to nematodes. This would advance our understanding of host-pathogen dynamics and foster innovative therapeutic strategies.

Towards precision parasite management for livestock gastrointestinal nematodes in 2030

The management of parasitic nematodes calls for a shift from conventional, indiscriminate, anthelmintic use to a more precise approach, directed by diagnostics. We should accept those parasite infection intensities that have minimal impact on production and welfare rather than attempt to eliminate them. The diagnostic toolbox for gastrointestinal nematodes (GINs) faces challenges due to anthelmintic resistance (AR), which is species-specific, drug-class-specific, and varies by region. We discuss which traditional tools may become obsolete and which tools need development to gain widespread use. Social science research highlights the need for dialogue between farmers and veterinarians that emphasises effective parasite management and upskilling the veterinary workforce for more sustainable practices centred on diagnostics to be adopted in practice by 2030.

Extensive transmission and variation in a functional receptor for praziquantel resistance in endemic Schistosoma mansoni

Mass-drug administration (MDA) of human populations using praziquantel monotherapy has become the primary strategy for controlling and potentially eliminating the major neglected tropical disease schistosomiasis. To understand how long-term MDA impacts schistosome populations, we analysed whole-genome sequence data of 570 Schistosoma mansoni samples (and the closely related outgroup species, S. rodhaini) from eight countries incorporating both publicly-available sequence data and new parasite material. This revealed broad-scale genetic structure across countries but with extensive transmission over hundreds of kilometres. We characterised variation across the transient receptor potential melastatin ion channel, TRPMPZQ, a target of praziquantel, which has recently been found to influence praziquantel susceptibility. Functional profiling of TRPMPZQ variants found in endemic populations identified four mutations that reduced channel sensitivity to praziquantel, indicating standing variation for resistance. Analysis of parasite infrapopulations sampled from individuals pre- and post-treatment identified instances of treatment failure, further indicative of potential praziquantel resistance. As schistosomiasis is targeted for elimination as a public health problem by 2030 in all currently endemic countries, and even interruption of transmission in selected African regions, we provide an in-depth genomic characterisation of endemic populations and an approach to identify emerging praziquantel resistance alleles.

Spirocerca lupi draft genome, vaccine and anthelmintic targets

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

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.

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.

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.

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.

The proof is in the poo-ding: Benefits of the longitudinal molecular surveillance of drug resistance demonstrated in a New South Wales cattle herd

Our understanding of anthelmintic resistance in the gastrointestinal nematodes of Australian cattle relies exclusively on small-scale phenotypic reports utilising traditional faecal egg count reduction tests. This approach is not readily scalable to establish the national prevalence of resistance, nor is it conducive of routine longitudinal surveillance for the emergence of resistance in its early stages. This study introduces the benefits of applying mixed amplicon metabarcoding longitudinally for timely and cost-efficient molecular surveillance of multiple anthelmintic resistance mutations, as they emerge on farms. Using opportunistically collected faecal samples from a cattle herd in central west New South Wales (2019-2023), we detected the early emergence of Haemonchus spp. levamisole-resistant S168T shortly after levamisole introduction, while benzimidazole-resistant allele frequencies remained constant. Additionally, we observed the possible spill-over of resistant Haemonchus contortus from sheep, along with variations in faecal burdens and species diversity influenced by climate stochasticity and host immunity. This study emphasises the power of molecular diagnostics for farm-level anthelmintic resistance management, providing essential evidence to support its integration into routine surveillance programmes.

Functional validation of novel levamisole resistance marker S168T in Haemonchus contortus

Recently, a S168T variant in the acetylcholine receptor subunit ACR-8 was associated with levamisole resistance in the parasitic helminth Haemonchus contortus. Here, we used the Xenopus laevis oocyte expression system and two-electrode voltage-clamp electrophysiology to measure the functional impact of this S168T variant on the H. contortus levamisole-sensitive acetylcholine receptor, L-AChR-1.1. Expression of the ACR-8 S168T variant significantly reduced the current amplitude elicited by levamisole compared to acetylcholine, with levamisole changing from a full to partial agonist on the recombinant L-AChR. Functional validation of the S168T mutation on modulating levamisole activity at the receptor level highlights its critical importance as both a mechanism and a marker of levamisole resistance.

Altered larval activation response associated with multidrug resistance in the canine hookworm Ancylostoma caninum

Parasitic gastrointestinal nematodes pose significant health risks to humans, livestock, and companion animals, and their control relies heavily on the use of anthelmintic drugs. Overuse of these drugs has led to the emergence of resistant nematode populations. Herein, a naturally occurring isolate (referred to as BCR) of the dog hookworm, Ancylostoma caninum, that is resistant to 3 major classes of anthelmintics is characterized. Various drug assays were used to determine the resistance of BCR to thiabendazole, ivermectin, moxidectin and pyrantel pamoate. When compared to a drug-susceptible isolate of A. caninum, BCR was shown to be significantly resistant to all 4 of the drugs tested. Multiple single nucleotide polymorphisms have been shown to impart benzimidazole resistance, including the F167Y mutation in the β-tubulin isotype 1 gene, which was confirmed to be present in BCR through molecular analysis. The frequency of the resistant allele in BCR was 76.3% following its first passage in the lab, which represented an increase from approximately 50% in the founding hookworm population. A second, recently described mutation in codon 134 (Q134H) was also detected at lower frequency in the BCR population. Additionally, BCR exhibits an altered larval activation phenotype compared to the susceptible isolate, suggesting differences in the signalling pathways involved in the activation process which may be associated with resistance. Further characterization of this isolate will provide insights into the mechanisms of resistance to macrocyclic lactones and tetrahydropyrimidine anthelmintics.

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.

Gene co-expression network analysis reveal core responsive genes in Parascaris univalens tissues following ivermectin exposure

Anthelmintic resistance in equine parasite Parascaris univalens, compromises ivermectin (IVM) effectiveness and necessitates an in-depth understanding of its resistance mechanisms. Most research, primarily focused on holistic gene expression analyses, may overlook vital tissue-specific responses and often limit the scope of novel genes. This study leveraged gene co-expression network analysis to elucidate tissue-specific transcriptional responses and to identify core genes implicated in the IVM response in P. univalens. Adult worms (n = 28) were exposed to 10-11 M and 10-9 M IVM in vitro for 24 hours. RNA-sequencing examined transcriptional changes in the anterior end and intestine. Differential expression analysis revealed pronounced tissue differences, with the intestine exhibiting substantially more IVM-induced transcriptional activity. Gene co-expression network analysis identified seven modules significantly associated with the response to IVM. Within these, 219 core genes were detected, largely expressed in the intestinal tissue and spanning diverse biological processes with unspecific patterns. After 10-11 M IVM, intestinal tissue core genes showed transcriptional suppression, cell cycle inhibition, and ribosomal alterations. Interestingly, genes PgR028_g047 (sorb-1), PgB01_g200 (gmap-1) and PgR046_g017 (col-37 & col-102) switched from downregulation at 10-11 M to upregulation at 10-9 M IVM. The 10-9 M concentration induced expression of cuticle and membrane integrity core genes in the intestinal tissue. No clear core gene patterns were visible in the anterior end after 10-11 M IVM. However, after 10-9 M IVM, the anterior end mostly displayed downregulation, indicating disrupted transcriptional regulation. One interesting finding was the non-modular calcium-signaling gene, PgR047_g066 (gegf-1), which uniquely connected 71 genes across four modules. These genes were enriched for transmembrane signaling activity, suggesting that PgR047_g066 (gegf-1) could have a key signaling role. By unveiling tissue-specific expression patterns and highlighting biological processes through unbiased core gene detection, this study reveals intricate IVM responses in P. univalens. These findings suggest alternative drug uptake of IVM and can guide functional validations to further IVM resistance mechanism understanding.

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

Ascaridia galli - An old problem that requires new solutions

Reports of Ascaridia galli in laying hens in Europe have increased since the ban on conventional battery cages in 2012. As this parasite is transmitted directly via the faecal-oral route by parasite eggs containing a larva, it is reasonable to assume that the escalating problem is related to the increased exposure now occurring in modern welfare-friendly cage-free housing systems. On many farms, A. galli reappears in subsequent flocks, even though the birds have no access to the outdoors, biosecurity is high and empty houses are cleaned and disinfected during downtime. Since the egg production cycle lasts only ≈80 weeks and recombinant antigen production for helminth vaccines has not yet been solved, the development of a vaccine seems to be an unrealistic option. Therefore, disrupting the life cycle of the parasite by other means, including the strategic use of dewormers, appears to be the key to controlling infection. Of concern is that only one class of anthelmintics is licenced for poultry in Europe and that are usually administered indiscriminately through the birds' drinking water and often too late when the parasite is already established. If current calendar-based parasite control strategies are not changed, there is a risk that resistance to anthelmintics may develop, as has already been demonstrated with nematodes in livestock. We insist that treatments can be more effective and the risk of developing drug resistance can be mitigated if we invest in a better understanding of A. galli responses to more prudent and judicious use of anthelmintics. This review identifies knowledge gaps and highlights aspects of sustainable parasite control that require further research to support commercial egg producers.

Understanding anthelmintic resistance in livestock using "omics" approaches

Widespread and improper use of various anthelmintics, genetic, and epidemiological factors has resulted in anthelmintic-resistant (AR) helminth populations in livestock. This is currently quite common globally in different livestock animals including sheep, goats, and cattle to gastrointestinal nematode (GIN) infections. Therefore, the mechanisms underlying AR in parasitic worm species have been the subject of ample research to tackle this challenge. Current and emerging technologies in the disciplines of genomics, transcriptomics, metabolomics, and proteomics in livestock species have advanced the understanding of the intricate molecular AR mechanisms in many major parasites. The technologies have improved the identification of possible biomarkers of resistant parasites, the ability to find actual causative genes, regulatory networks, and pathways of parasites governing the AR development including the dynamics of helminth infection and host-parasite infections. In this review, various "omics"-driven technologies including genome scan, candidate gene, quantitative trait loci, transcriptomic, proteomic, and metabolomic approaches have been described to understand AR of parasites of veterinary importance. Also, challenges and future prospects of these "omics" approaches are also discussed.

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.

Recent Advances in the Control of Endoparasites in Ruminants from a Sustainable Perspective

Consumer awareness of animal welfare and environmental health has led to a plateau level of global consumption putting serious pressure on the livestock industry [...].

Identification of cytochrome P450 gene family and functional analysis of HgCYP33E1 from Heterodera glycines

The cytochrome P450 (CYP) genes of nematode play a crucial role in the metabolic detoxification of xenobiotics including pesticides. Heterodera glycines, also known as the soybean cyst nematode, is a sedentary endoparasite that infests plant roots, causing high annual economic losses in soybean production regions globally. In this study, we identified 36 CYP genes at a genome-wide level of the H. glycines isolate TN10 using all CYPs from Caenorhabditis elegans as queries. Subsequently, a full-length cDNA of HgCYP33E1 which was significantly up-regulated by the conventional nematicide abamectin was initially cloned from H. glycines. It presented significantly higher expressions in the second-stage juvenile (J2) compared to other parasitic stages of H. glycines. qRT-PCR analysis suggested that the expression of HgCYP33E1 was also xenobiotically induced by soybean root exudate and the metabolites of biocontrol agents. Using RNA interference (RNAi), we investigated the function of HgCYP33E1 in H. glycines parasitism and nematicide selectivity. Compared to the control and dsGFP-treated group, silencing of HgCYP33E1 did not affect the J2 behaviors and the early invasion ability, while it decreased the number of J4s in soybean roots after 18-d inoculation with the dsHgCYP33E1-treated nematodes. In addition, knockdown of HgCYP33E1 in H. glycines resulted in an increase in J2 mortality after 24-h incubation with abamectin compared to the GFP dsRNA-soaked and the control group. These findings revealed the potential role of HgCYP33E1 in the xenobiotic detoxification pathway of H. glycines. Moreover, our data also provided valuable gene information for studying the functions of the CYP family in H. glycines host adaption.

New paradigms in research on Dirofilaria immitis

BACKGROUND

Since the advent of ivermectin (along with melarsomine and doxycycline), heartworm has come to be viewed as a solved problem in veterinary medicine, diminishing investment into non-clinical research on Dirofilaria immitis. However, heartworm infections continue to pose problems for practitioners and their patients and seem to be increasing in frequency and geographic distribution. Resistance to preventative therapies (macrocyclic lactones) complicates the picture. The use of chemotherapy for other kinds of pathogens has benefitted enormously from research into the basic biology of the pathogen and on the host-pathogen interface. A lack of basic information on heartworms as parasites and how they interact with permissive and non-permissive hosts greatly limits the ability to discover new ways to prevent and treat heartworm disease. Recent advances in technical platforms will help overcome the intrinsic barriers that hamper research on D. immitis, most notably, the need for experimentally infected dogs to maintain the life cycle and provide material for experiments. Impressive achievements in the development of laboratory animal models for D. immitis will enhance efforts to discover new drugs for prevention or treatment, to characterize new diagnostic biomarkers and to identify key parasite-derived molecules that are essential for survival in permissive hosts, providing a rational basis for vaccine discovery. A 'genomics toolbox' for D. immitis could enable unprecedented insight into the negotiations between host and parasite that enable survival in a permissive host. The more we know about the pathogen and how it manipulates its host, the better able we will be to protect companion animals far into the future.

The turkey ascarid, Ascaridia dissimilis, as a model genetic system

Parasitic nematodes cause significant effects on humans each year, with the most prevalent being Ascaris lumbricoides. Benzimidazoles (BZ) are the most widely used anthelmintic drug in humans, and although the biology of resistance to this drug class is understood in some species, resistance is poorly characterized in ascarids. Models such as Caenorhabditis elegans were essential in developing our current understanding of BZ resistance, but more closely related model nematodes are needed to understand resistance in ascarids. Here, we propose a new ascarid model species that infects turkeys, Ascaridia dissimilis, to develop a better understanding of BZ resistance.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Intraspecific diversity in the mechanisms underlying abamectin resistance in a cosmopolitan pest

Pesticide resistance relies on a myriad of mechanisms, ranging from single mutations to a complex and polygenic architecture, and it involves mechanisms such as target-site insensitivity, metabolic detoxification, or a combination of these, with either additive or synergistic effects. Several resistance mechanisms against abamectin, a macrocyclic lactone widely used in crop protection, have been reported in the cosmopolitan pest Tetranychus urticae. However, it has been shown that a single mechanism cannot account for the high levels of abamectin resistance found across different mite populations. Here, we used experimental evolution combined with bulked segregant analyses to map quantitative trait loci (QTL) associated with abamectin resistance in two genetically unrelated populations of T. urticae. In these two independent QTL mapping experiments, three and four QTLs were identified, of which three were shared between experiments. Shared QTLs contained genes encoding subunits of the glutamate-gated chloride channel (GluCl) and harboured previously reported mutations, including G314D in GluCl1 and G326E in GluCl3, but also novel resistance candidate loci, including DNA helicases and chemosensory receptors. Surprisingly, the fourth QTL, present only in only one of the experiments and thus unique for one resistant parental line, revealed a non-functional variant of GluCl2, suggesting gene knock-out as resistance mechanism. Our study uncovers the complex basis of abamectin resistance, and it highlights the intraspecific diversity of genetic mechanisms underlying resistance in a cosmopolitan pest.

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.

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.

Tanacetum species: Bridging empirical knowledge, phytochemistry, nutritional value, health benefits and clinical evidence

Introduction: The Tanacetum genus consists of 160 accepted flowering species thriving throughout temperate regions, mainly in the Mediterranean Basin, Northern America, and southwestern and eastern Asia. Tanacetum species bear a long-standing record of use in the folk medicine of indigenous tribes and communities worldwide, along with multitudinous applications in traditional cuisines, cosmeceuticals, and agricultural fields. Methods: Up-to-date data related to traditional uses, phytochemistry, biological activities, toxicity and clinical trials of the genus Tanacetum were systematically reviewed from several online scientific engines, including PubMed, Web of Science, Scopus, SciFinder, Wiley Online, Science Direct, and Cochrane library. Results and discussion: Over the past three decades, 241 metabolites have been isolated from nearly twenty species, including phenolic acids, flavonoids, coumarins, fatty acids and alkanes, aldehydes, volatile compounds, and naphthoquinones. Some unique metabolites have also been identified, such as the ceramides tanacetamide (A-D) from T. artemisioides, pyrethrins from T. cinerariifolium, and sesquiterpene lactones from several species. However, these secondary metabolites are still poorly studied despite in vitro clues highlighting their colossal pharmacological properties, especially as hypotensive, neuroprotective, anticancer, and antimicrobial agents. Scientific studies have validated some traditional claims of the plant, such as antidiabetic, anticancer, anthelmintic, insecticide, antioxidant, and hepatoprotective activities, as well as against festering wounds, skin ulcers, urinary tract infections, and sexually transmitted diseases. Other ethnomedicinal uses for arthritis, gout, rheumatism, anemia, and as a litholytic, antivenom and diaphoretic have not yet been supported and would constitute the subject of further research.