Population genomic and evolutionary modelling analyses reveal a single major QTL for ivermectin drug resistance in the pathogenic nematode, Haemonchus 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.

Anti-inflammatory and immunomodulatory effects of mesenchymal stem cell therapy on parasitic drug resistance

INTRODUCTION

The emergence of antiparasitic drug resistance poses a concerning threat to animals and humans. Mesenchymal Stem Cells (MSCs) have been widely used to treat infections in humans, pets, and livestock. Although this is an emerging field of study, the current review outlines possible mechanisms and examines potential synergism in combination therapies and the possible harmful effects of such an approach.

AREAS COVERED

The present study delved into the latest pre-clinical research on utilizing MSCs to treat parasitic infections. As per investigations, the introduction of MSCs to patients grappling with parasitic diseases like schistosomiasis, malaria, cystic echinococcosis, toxoplasmosis, leishmaniasis, and trypanosomiasis has shown a reduction in parasite prevalence. This intervention also alters the levels of both pro- and anti-inflammatory cytokines. Furthermore, the combined administration of MSCs and antiparasitic drugs has demonstrated enhanced efficacy in combating parasites and modulating the immune response.

EXPERT OPINION

Mesenchymal stem cells are a potential solution for addressing parasitic drug resistance. This is mainly because of their remarkable immunomodulatory abilities, which can potentially help combat parasites' resistance to drugs.

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.

Genome structure and population genomics of the canine heartworm Dirofilaria immitis

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

Advancing Strongyloides omics data: bridging the gap with Caenorhabditis elegans

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

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.

Hookworms dynamically respond to loss of Type 2 immune pressure

The impact of the host immune environment on parasite transcription and fitness is currently unknown. It is widely held that hookworm infections have an immunomodulatory impact on the host, but whether the converse is true remains unclear. Immunity against adult-stage hookworms is largely mediated by Type 2 immune responses driven by the transcription factor Signal Transducer and Activator of Transcription 6 (STAT6). This study investigated whether serial passage of the rodent hookworm Nippostrongylus brasiliensis in STAT6-deficient mice (STAT6 KO) caused changes in parasites over time. After adaptation to STAT6 KO hosts, N. brasiliensis increased their reproductive output, feeding capacity, energy content, and body size. Using an improved N. brasiliensis genome, we found that these physiological changes corresponded with a dramatic shift in the transcriptional landscape, including increased expression of gene pathways associated with egg production, but a decrease in genes encoding neuropeptides, proteases, SCP/TAPS proteins, and transthyretin-like proteins; the latter three categories have been repeatedly observed in hookworm excreted/secreted proteins (ESPs) implicated in immunosuppression. Although transcriptional changes started to appear in the first generation of passage in STAT6 KO hosts for both immature and mature adult stages, downregulation of the genes putatively involved in immunosuppression was only observed after multiple generations in this immunodeficient environment. When STAT6 KO-adapted N. brasiliensis were reintroduced to a naive WT host after up to 26 generations, this progressive change in host-adaptation corresponded to increased production of inflammatory cytokines by the WT host. Surprisingly, however, this single exposure of STAT6 KO-adapted N. brasiliensis to WT hosts resulted in worms that were morphologically and transcriptionally indistinguishable from WT-adapted parasites. This work uncovers remarkable plasticity in the ability of hookworms to adapt to their hosts, which may present a general feature of parasitic nematodes.

Mapping resistance-associated anthelmintic interactions in the model nematode Caenorhabditis elegans

Parasitic nematodes infect billions of people and are mainly controlled by anthelmintic mass drug administration (MDA). While there are growing efforts to better understand mechanisms of anthelmintic resistance in human and animal populations, it is unclear how resistance mechanisms that alter susceptibility to one drug affect the interactions and efficacy of drugs used in combination. Mutations that alter drug permeability across primary nematode barriers have been identified as potential resistance mechanisms using the model nematode Caenorhabditis elegans. We leveraged high-throughput assays in this model system to measure altered anthelmintic susceptibility in response to genetic perturbations of potential cuticular, amphidial, and alimentary routes of drug entry. Mutations in genes associated with these tissue barriers differentially altered susceptibility to the major anthelmintic classes (macrocyclic lactones, benzimidazoles, and nicotinic acetylcholine receptor agonists) as measured by animal development. We investigated two-way anthelmintic interactions across C. elegans genetic backgrounds that confer resistance or hypersensitivity to one or more drugs. We observe that genetic perturbations that alter susceptibility to a single drug can shift the drug interaction landscape and lead to the appearance of novel synergistic and antagonistic interactions. This work establishes a framework for investigating combinatorial therapies in model nematodes that can potentially be translated to amenable parasite species.

Mapping resistance-associated anthelmintic interactions in the model nematode Caenorhabditis elegans

Parasitic nematodes infect billions of people and are mainly controlled by anthelmintic mass drug administration (MDA). While there are growing efforts to better understand mechanisms of anthelmintic resistance in human and animal populations, it is unclear how resistance mechanisms that alter susceptibility to one drug affect the interactions and efficacy of drugs used in combination. Mutations that alter drug permeability across primary nematode barriers have been identified as potential resistance mechanisms using the model nematode Caenorhabditis elegans. We leveraged high-throughput assays in this model system to measure altered anthelmintic susceptibility in response to genetic perturbations of potential cuticular, amphidial, and alimentary routes of drug entry. Mutations in genes associated with these tissue barriers differentially altered susceptibility to the major anthelmintic classes (macrocyclic lactones, benzimidazoles, and nicotinic acetylcholine receptor agonists) as measured by animal development. We investigated two-way anthelmintic interactions across C. elegans genetic backgrounds that confer resistance or hypersensitivity to one or more drugs. We observe that genetic perturbations that alter susceptibility to a single drug can shift the drug interaction landscape and lead to the appearance of novel synergistic and antagonistic interactions. This work establishes a framework for investigating combinatorial therapies in model nematodes that can potentially be translated to amenable parasite species.

Variation in anthelmintic responses are driven by genetic differences among diverse C. elegans wild strains

Treatment of parasitic nematode infections in humans and livestock relies on a limited arsenal of anthelmintic drugs that have historically reduced parasite burdens. However, anthelmintic resistance (AR) is increasing, and little is known about the molecular and genetic causes of resistance for most drugs. The free-living roundworm Caenorhabditis elegans has proven to be a tractable model to understand AR, where studies have led to the identification of molecular targets of all major anthelmintic drug classes. Here, we used genetically diverse C. elegans strains to perform dose-response analyses across 26 anthelmintic drugs that represent the three major anthelmintic drug classes (benzimidazoles, macrocyclic lactones, and nicotinic acetylcholine receptor agonists) in addition to seven other anthelmintic classes. First, we found that C. elegans strains displayed similar anthelmintic responses within drug classes and significant variation across drug classes. Next, we compared the effective concentration estimates to induce a 10% maximal response (EC10) and slope estimates of each dose-response curve of each strain to the laboratory reference strain, which enabled the identification of anthelmintics with population-wide differences to understand how genetics contribute to AR. Because genetically diverse strains displayed differential susceptibilities within and across anthelmintics, we show that C. elegans is a useful model for screening potential nematicides before applications to helminths. Third, we quantified the levels of anthelmintic response variation caused by genetic differences among individuals (heritability) to each drug and observed a significant correlation between exposure closest to the EC10 and the exposure that exhibited the most heritable responses. These results suggest drugs to prioritize in genome-wide association studies, which will enable the identification of AR genes.

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.

A major locus confers triclabendazole resistance in Fasciola hepatica and shows dominant inheritance

Fasciola hepatica infection is responsible for substantial economic losses in livestock worldwide and poses a threat to human health in endemic areas. The mainstay of control in livestock and the only drug licenced for use in humans is triclabendazole (TCBZ). TCBZ resistance has been reported on every continent and threatens effective control of fasciolosis in many parts of the world. To date, understanding the genetic mechanisms underlying TCBZ resistance has been limited to studies of candidate genes, based on assumptions of their role in drug action. Taking an alternative approach, we combined a genetic cross with whole-genome sequencing to localise a ~3.2Mbp locus within the 1.2Gbp F. hepatica genome that confers TCBZ resistance. We validated this locus independently using bulk segregant analysis of F. hepatica populations and showed that it is the target of drug selection in the field. We genotyped individual parasites and tracked segregation and reassortment of SNPs to show that TCBZ resistance exhibits Mendelian inheritance and is conferred by a dominant allele. We defined gene content within this locus to pinpoint genes involved in membrane transport, (e.g. ATP-binding cassette family B, ABCB1), transmembrane signalling and signal transduction (e.g. GTP-Ras-adenylyl cyclase and EGF-like protein), DNA/RNA binding and transcriptional regulation (e.g. SANT/Myb-like DNA-binding domain protein) and drug storage and sequestration (e.g. fatty acid binding protein, FABP) as prime candidates for conferring TCBZ resistance. This study constitutes the first experimental cross and genome-wide approach for any heritable trait in F. hepatica and is key to understanding the evolution of drug resistance in Fasciola spp. to inform deployment of efficacious anthelmintic treatments in the field.

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.

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.

Ivermectin and albendazole coadministration: opportunities for strongyloidiasis control

In 2020, WHO recognised the importance of strongyloidiasis alongside soil-transmitted helminths (STH) in their 2021-30 roadmap, which aspires to target Strongyloides stercoralis with preventive chemotherapy by use of ivermectin. Combination treatment with both albendazole, the primary drug used to treat STH, and ivermectin, would improve the efficiency of mass drug administration targeting both STH and S stercoralis. In this Personal View, we discuss the challenges and opportunities towards the development of an efficient control programme for strongyloidiasis, particularly if it is to run concurrently with STH control. We argue the need to define the prevalence threshold to implement preventive chemotherapy for S stercoralis, the target populations and optimal dosing schedules, and discuss the added benefits of a fixed-dose coformulation of ivermectin and albendazole. Implementation of an efficient control programme will require improvements to current diagnostics, and validation of new diagnostics, to target and monitor S stercoralis infections, and consideration of the challenges of multispecies diagnostics for S stercoralis and STH control. Finally, the evolution of ivermectin resistance represents a credible risk to control S stercoralis; we argue that genome-wide approaches, together with improved genome resources, are needed to characterise and prevent the emergence of resistance. Overcoming these challenges will help to reduce strongyloidiasis burden and enhance the feasibility of controlling it worldwide.

Genomic landscape of drug response reveals mediators of anthelmintic resistance

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

Polymorphisms in exon 11 of the mptl-1 gene and monepantel resistance in Haemonchus contortus

Chemical treatments are the main strategy to control gastrointestinal nematodes in sheep, and the emergence of anthelmintic resistance, as consequence, results in control failures and leads to economic losses. Thus, molecular tests may constitute an excellent tool for the early detection of anthelmintic resistance-related mutations. Thus, a polymerase chain reaction (PCR)-based genotyping assay followed by polyacrylamide gel electrophoresis (PAGE) was developed to detect polymorphisms in exon 11 of the acetylcholine receptor monepantel-1 gene (mptl-1) that were previously associated with monepantel resistance through a genome-wide study in Haemonchus contortus. DNA samples recovered from individual and pooled third-stage larvae from two susceptible field-derived isolates and five (three in vivo-derived and two field-derived) resistant populations were used. New polymorphisms, including a 6-bp deletion and a 3-bp insertion, were detected in resistant individuals. These indels, confirmed using sequencing of cloned PCR products, are predicted to result in amino acid changes in transmembrane domain 2 (TMD2) of the MPTL-1 protein. The two susceptible isolates showed only the presence of the wild-type allele (100%), whereas lower frequencies of the wild-type allele were detected in monepantel-resistant populations (11.1 to 66.7%). These findings report new polymorphisms in the mptl-1 gene, validate the results obtained through genomic mapping for monepantel resistance, and provide a PCR-based assay to genotype indels located in exon 11 of mptl-1 in H. contortus.

Improving helminth genome resources in the post-genomic era

Rapid advancement in high-throughput sequencing and analytical approaches has seen a steady increase in the generation of genomic resources for helminth parasites. Now, helminth genomes and their annotations are a cornerstone of numerous efforts to compare genetic and transcriptomic variation, from single cells to populations of globally distributed parasites, to genome modifications to understand gene function. Our understanding of helminths is increasingly reliant on these genomic resources, which are primarily static once published and vary widely in quality and completeness between species. This article seeks to highlight the cause and effect of this variation and argues for the continued improvement of these genomic resources - even after their publication - which is necessary to provide a more accurate and complete understanding of the biology of these important pathogens.

Genomic signatures of selection associated with benzimidazole drug treatments in Haemonchus contortus field populations

Genome-wide methods offer a powerful approach to detect signatures of drug selection. However, limited availability of suitable reference genomes and the difficulty of obtaining field populations with well-defined, distinct drug treatment histories mean there is little information on the signatures of selection in parasitic nematodes and on how best to detect them. This study addresses these knowledge gaps by using field populations of Haemonchus contortus with well-defined benzimidazole treatment histories, leveraging a recently completed chromosomal-scale reference genome assembly. We generated a panel of 49,393 genomic markers to genotype 20 individual adult worms from each of four H. contortus populations: two from closed sheep flocks with an approximate 20 year history of frequent benzimidazole treatment, and two populations with a history of little or no treatment. Sampling occurred in the same geographical region to limit genetic differentiation and maximise the detection sensitivity. A clear signature of selection was detected on chromosome I, centred on the isotype-1 β-tubulin gene. Two additional, but weaker, signatures of selection were detected; one near the middle of chromosome I spanning 3.75 Mbp and 259 annotated genes, and one on chromosome II spanning a region of 3.3 Mbp and 206 annotated genes, including the isotype-2 β-tubulin locus. We also assessed how sensitivity was impacted by sequencing depth, worm number, and pooled versus individual worm sequence data. This study provides the first known direct genome-wide evidence for any parasitic nematode, that the isotype-1 β-tubulin gene is quantitatively the single most important benzimidazole resistance locus. It also identified two additional genomic regions that likely contain benzimidazole resistance loci of secondary importance. This study provides an experimental framework to maximise the power of genome-wide approaches to detect signatures of selection driven by anthelmintic drug treatments in field populations of parasitic nematodes.

Transgenic Expression of Haemonchus contortus Cytochrome P450 Hco-cyp-13A11 Decreases Susceptibility to Particular but Not All Macrocyclic Lactones in the Model Organism Caenorhabditis elegans

The number of reported macrocyclic lactones (ML) resistance cases across all livestock hosts is steadily increasing. Different studies in the parasitic nematode Haemonchus contortus assume the participation of cytochrome P450s (Cyps) enzymes in ML resistance. Still, functional data about their individual contribution to resistance or substrate specificity is missing. Via microinjection, transgenic Caenorhabditis elegans expressing HCON_00141052 (transgene-Hco-cyp-13A11) from extrachromosomal arrays were generated. After 24 h of exposure to different concentrations of ivermectin (IVM), ivermectin aglycone (IVMa), selamectin (SEL), doramectin (DRM), eprinomectin (EPR), and moxidectin (MOX), motility assays were performed to determine the impact of the H. contortus Cyp to the susceptibility of the worms against each ML. While transgene-Hco-cyp-13A11 significantly decreased susceptibility to IVM (four-fold), IVMa (2-fold), and SEL (3-fold), a slight effect for DRM and no effect for MOX, and EPR was observed. This substrate specificity of Hco-cyp-13A11 could not be explained by molecular modeling and docking studies. Hco-Cyp-13A11 molecular models were obtained for alleles from isolates with different resistance statuses. Although 14 amino acid polymorphisms were detected, none was resistance specific. In conclusion, Hco-cyp-13A11 decreased IVM, IVMa, and SEL susceptibility to a different extent, but its potential impact on ML resistance is not driven by polymorphisms.

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

Genome-wide analysis of Schistosoma mansoni reveals limited population structure and possible praziquantel drug selection pressure within Ugandan hot-spot communities

Populations within schistosomiasis control areas, especially those in Africa, are recommended to receive regular mass drug administration (MDA) with praziquantel (PZQ) as the main strategy for controlling the disease. The impact of PZQ treatment on schistosome genetics remains poorly understood, and is limited by a lack of high-resolution genetic data on the population structure of parasites within these control areas. We generated whole-genome sequence data from 174 individual miracidia collected from both children and adults from fishing communities on islands in Lake Victoria in Uganda that had received either annual or quarterly MDA with PZQ over four years, including samples collected immediately before and four weeks after treatment. Genome variation within and between samples was characterised and we investigated genomic signatures of natural selection acting on these populations that could be due to PZQ treatment. The parasite population on these islands was more diverse than found in nearby villages on the lake shore. We saw little or no genetic differentiation between villages, or between the groups of villages with different treatment intensity, but slightly higher genetic diversity within the pre-treatment compared to post-treatment parasite populations. We identified classes of genes significantly enriched within regions of the genome with evidence of recent positive selection among post-treatment and intensively treated parasite populations. The differential selection observed in post-treatment and pre-treatment parasite populations could be linked to any reduced susceptibility of parasites to praziquantel treatment.

Vaccine-induced time- and age-dependent mucosal immunity to gastrointestinal parasite infection

Individuals vary broadly in their response to vaccination and subsequent challenge infection, with poor vaccine responders causing persistence of both infection and transmission in populations. Yet despite having substantial economic and societal impact, the immune mechanisms that underlie such variability, especially in infected tissues, remain poorly understood. Here, to characterise how antihelminthic immunity at the mucosal site of infection developed in vaccinated lambs, we inserted gastric cannulae into the abomasa of three-month- and six-month-old lambs and longitudinally analysed their local immune response during subsequent challenge infection. The vaccine induced broad changes in pre-challenge abomasal immune profiles and reduced parasite burden and egg output post-challenge, regardless of age. However, age affected how vaccinated lambs responded to infection across multiple immune pathways: adaptive immune pathways were typically age-dependent. Identification of age-dependent and age-independent protective immune pathways may help refine the formulation of vaccines, and indicate specificities of pathogen-specific immunity more generally.

Transcriptomic analyses implicate neuronal plasticity and chloride homeostasis in ivermectin resistance and response to treatment in a parasitic nematode

The antiparasitic drug ivermectin plays an essential role in human and animal health globally. However, ivermectin resistance is widespread in veterinary helminths and there are growing concerns of sub-optimal responses to treatment in related helminths of humans. Despite decades of research, the genetic mechanisms underlying ivermectin resistance are poorly understood in parasitic helminths. This reflects significant uncertainty regarding the mode of action of ivermectin in parasitic helminths, and the genetic complexity of these organisms; parasitic helminths have large, rapidly evolving genomes and differences in evolutionary history and genetic background can confound comparisons between resistant and susceptible populations. We undertook a controlled genetic cross of a multi-drug resistant and a susceptible reference isolate of Haemonchus contortus, an economically important gastrointestinal nematode of sheep, and ivermectin-selected the F2 population for comparison with an untreated F2 control. RNA-seq analyses of male and female adults of all populations identified high transcriptomic differentiation between parental isolates, which was significantly reduced in the F2, allowing differences associated specifically with ivermectin resistance to be identified. In all resistant populations, there was constitutive upregulation of a single gene, HCON_00155390:cky-1, a putative pharyngeal-expressed transcription factor, in a narrow locus on chromosome V previously shown to be under ivermectin selection. In addition, we detected sex-specific differences in gene expression between resistant and susceptible populations, including constitutive upregulation of a P-glycoprotein, HCON_00162780:pgp-11, in resistant males only. After ivermectin selection, we identified differential expression of genes with roles in neuronal function and chloride homeostasis, which is consistent with an adaptive response to ivermectin-induced hyperpolarisation of neuromuscular cells. Overall, we show the utility of a genetic cross to identify differences in gene expression that are specific to ivermectin selection and provide a framework to better understand ivermectin resistance and response to treatment in parasitic helminths.

Parasitic helminths (worms) infect people and animals throughout the world and are largely controlled with mass administration of anthelmintic drugs. There are a very limited number of anthelmintics available and parasitic helminths can rapidly develop resistance to these drugs. Ivermectin is a widely used anthelmintic in both humans and animals, but resistance is now widespread in the veterinary field. We crossed ivermectin resistant and ivermectin susceptible parasitic helminths and treated them with ivermectin or left them as untreated controls. This provided resistant and susceptible populations with a similar genetic background with which to study differences in gene expression associated with ivermectin resistance. We identified upregulation of a gene with no previous association with drug resistance (HCON_00155390:cky-1) in male and female worms in all resistant populations. This gene is thought to be expressed in the helminth pharynx (mouthpart) and, in mammals, plays a role in controlling nerve function and protecting nerves from damage. This is consistent with the known effects of ivermectin in inhibiting helminth feeding through pharyngeal paralysis and may implicate a novel mechanism that allows resistant worms to survive treatment.

Genome-wide analysis of the response to ivermectin treatment by a Swedish field population of Haemonchus contortus

Haemonchus contortus is a pathogenic gastrointestinal nematode of small ruminants and, in part due to its capacity to develop resistance to drugs, contributes to significant losses in the animal production sector worldwide. Despite decades of research, comparatively little is known about the specific mechanism(s) driving resistance to drugs such as ivermectin in this species. Here we describe a genome-wide approach to detect evidence of selection by ivermectin treatment in a field population of H. contortus from Sweden, using parasites sampled from the same animals before and seven days after ivermectin exposure followed by whole-genome sequencing. Despite an 89% reduction in parasites recovered after treatment measured by the fecal egg count reduction test, the surviving population was highly genetically similar to the population before treatment, suggesting that resistance has likely evolved over time and that resistance alleles are present on diverse haplotypes. Pairwise gene and SNP frequency comparisons indicated the highest degree of differentiation was found at the terminal end of chromosome 4, whereas the most striking difference in nucleotide diversity was observed in a region on chromosome 5 previously reported to harbor a major quantitative trait locus involved in ivermectin resistance. These data provide novel insight into the genome-wide effect of ivermectin selection in a field population as well as confirm the importance of the previously established quantitative trait locus in the development of resistance to ivermectin.

Phenotypic and genotypic analysis of benzimidazole resistance in reciprocal genetic crosses of Haemonchus contortus

Haemonchus contortus is arguably one of the most economically important and ubiquitous parasites of livestock globally and commonly involved in cases of anthelmintic resistance. Here, we performed reciprocal genetic crosses using susceptible (MHco3(ISE)) and multiple anthelmintic resistant (MHco18(UGA2004)) H. contortus isolates. Resultant admixed populations were designated MHco3/18 or MHco18/3, where the lead isolate reflects the origin of the females. Three independent filial generations were generated for each cross, which were subjected to bioassays, molecular approaches and population genetic analyses to investigate the phenotypic and genotypic inheritance of benzimidazole (BZ) resistance at each stage. A panel of microsatellite markers confirmed the success of the genetic cross as markers from both parents were seen in the F1 crosses. Egg hatch tests revealed a stark difference between the two F1 crosses with ED50 estimates for MHco18/3 being 9 times greater than those for MHco3/18. Resistance factors based on ED50 estimates ranged from 6 to 57 fold in the filial progeny compared to MHco3(ISE) parents. Molecular analysis of the F167Y and F200Y SNP markers associated with BZ resistance were analysed by pyrosequencing and MiSeq deep amplicon sequencing, which showed that MHco3/18.F1 and MHco18/3.F1 both had similar frequencies of the F200Y resistant allele (45.3% and 44.3%, respectively), whereas for F167Y, MHco18/3.F1 had a two-fold greater frequency of the resistant-allele compared to MHco3/18.F1 (18.2% and 8.8%, respectively). Comparison between pyrosequencing and MiSeq amplicon sequencing revealed that the allele frequencies derived from both methods were concordant at codon 200 (rc = 0.97), but were less comparable for codon 167 (rc = 0.55). The use of controlled reciprocal genetic crosses have revealed a potential difference in BZ resistance phenotype dependent on whether the resistant allele is paternally or maternally inherited. These findings provide new insight and prompt further investigation into the inheritance of BZ resistance in H. contortus.

Anthelmintic resistance in ruminants: challenges and solutions

Anthelmintic resistance (AR) is a growing concern for effective parasite control in farmed ruminants globally. Combatting AR will require intensified and integrated research efforts in the development of innovative diagnostic tests to detect helminth infections and AR, sustainable anthelmintic treatment strategies and the development of complementary control approaches such as vaccination and plant-based control. It will also require a better understanding of socio-economic drivers of anthelmintic treatment decisions, in order to support a behavioural shift and develop targeted communication strategies that promote the uptake of evidence-based sustainable solutions. Here, we review the state-of-the-art in these different fields of research activity related to AR in helminths of livestock ruminants in Europe and beyond. We conclude that in the advent of new challenges and solutions emerging from continuing spread of AR and intensified research efforts, respectively, there is a strong need for transnational multi-actor initiatives. These should involve all key stakeholders to develop indicators of infection and sustainable control, set targets and promote good practices to achieve them.

Point of care colourimetric and lateral flow LAMP assay for the detection of Haemonchus contortus in ruminant faecal samples

In this study, we present an optimised colourimetric and a lateral flow LAMP assay for the detection of Haemonchus contortus in small ruminant faecal samples. Using a previously published LAMP primer set, we made use of commercially available colourimetric LAMP and lateral flow kits and combined this into an optimised diagnostic assay which was then tested on field faecal samples from Eastern and South-Eastern Hungary as well as a pure H. contortus egg faecal sample from Košice, Slovakia. Both assays showed no conflicts in visual detection of the results. Additionally, we modified and tested several centrifuge-free DNA extraction methods and one bead-beating egg lysis DNA extraction method to develop a true point of care protocol, as the source of the starting DNA is the main rate-limiting step in farm-level molecular diagnosis. Out of the various methods trialed, promising results were obtained with the magnetic bead extraction method. Sample solutions from the Fill-FLOTAC^® technique were also utilised, which demonstrated that it could be efficiently adapted for field-level egg concentration to extract DNA. This proof of concept study showed that isothermal amplification technologies with a colourimetric detection or when combined with a lateral flow assay could be an important step for a true point of care molecular diagnostic assay for H. contortus.

A journey through 50 years of research relevant to the control of gastrointestinal nematodes in ruminant livestock and thoughts on future directions

This review article provides an historical perspective on some of the major research advances of relevance to ruminant livestock gastrointestinal nematode control over the last 50 years. Over this period, gastrointestinal nematode control has been dominated by the use of broad-spectrum anthelmintic drugs. Whilst this has provided unprecedented levels of successful control for many years, this approach has been gradually breaking down for more than two decades and is increasingly unsustainable which is due, at least in part, to the emergence of anthelmintic drug resistance and a number of other factors discussed in this article. We first cover the remarkable success story of the discovery and development of broad-spectrum anthelmintic drugs, the changing face of anthelmintic drug discovery research and the emergence of anthelmintic resistance. This is followed by a review of some of the major advances in the increasingly important area of non-pharmaceutical gastrointestinal nematode control including immunology and vaccine development, epidemiological modelling and some of the alternative control strategies such as breeding for host resistance, refugia-based methods and biological control. The last 50 years have witnessed remarkable innovation and success in research aiming to improve ruminant livestock gastrointestinal nematode control, particularly given the relatively small size of the research community and limited funding. In spite of this, the growing global demand for livestock products, together with the need to maximise production efficiencies, reduce environmental impacts and safeguard animal welfare - as well as specific challenges such as anthelmintic drug resistance and climate change- mean that gastrointestinal nematode researchers will need to be as innovative in the next 50 years as in the last.

Anthelmintic resistance: is a solution possible?

More than 50 years after anthelmintic resistance was first identified, its prevalence and impact on the animal production industry continues to increase across the world. The term "anthelmintic resistance" (AR) can be briefly defined as the reduction in efficacy of a certain dose of anthelmintic drugs (AH) in eliminating the presence of a parasite population that was previously susceptible. The main aim of this study is to examine anthelmintic resistance in domestic herbivores. There are numerous factors playing a role in the development of AR, but the most important is livestock management. The price of AH and the need to treat a high number of animals mean that farmers face significant costs in this regard, yet, since 1981, little progress has been made in the discovery of new molecules and the time and cost required to bring a new AH to market has increased dramatically in recent decades. Furthermore, resistance has also emerged for new AH, such as monepantel or derquantel. Consequently, ruminant parasitism cannot be controlled solely by using synthetic chemicals. A change in approach is needed, using a range of preventive measures in order to achieve a sustainable control programme. The use of nematophagous fungi or of plant extracts rich in compounds with anthelmintic properties, such as terpenes, condensed tannins, or flavonoids, represent potential alternatives. Nevertheless, although new approaches are showing promising results, there is still much to do. More research focused on the control of AR is needed.

The Potential Role of Clinical Metagenomics in Infectious Diseases: Therapeutic Perspectives

Clinical metagenomics (CMg) is the process of sequencing nucleic acid of clinical samples to obtain clinically relevant information such as the identification of microorganisms and their susceptibility to antimicrobials. Over the last decades, sequencing and bioinformatic solutions supporting CMg have much evolved and an increasing number of case reports and series covering various infectious diseases have been published. Metagenomics is a new approach to infectious disease diagnosis that is currently being developed and is certainly one of the most promising for the coming years. However, most CMg studies are retrospective, and few address the potential impact CMg could have on patient management, including initiation, adaptation, or cessation of antimicrobials. In this narrative review, we have discussed the potential role of CMg in bacteriology, virology, mycology, and parasitology. Several reports and case-series confirm that CMg is an innovative tool with which one can (i) identify more microorganisms than with conventional methods in a single test, (ii) obtain results within hours, and (iii) tailor the antimicrobial regimen of patients. However, the cost-efficiency of CMg and its real impact on patient management are still to be determined.

A telomere-to-telomere assembly of Oscheius tipulae and the evolution of rhabditid nematode chromosomes

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, each chromosome has a single functional centromere with essential roles in spindle attachment and segregation. Fusion and fission can generate chromosomes with no or multiple centromeres, leading to genome instability. Groups with holocentric chromosomes (where centromeric function is distributed along each chromosome) might be expected to show karyotypic instability. This is generally not the case, and in Caenorhabditis elegans, it has been proposed that the role of maintenance of a stable karyotype has been transferred to the meiotic pairing centers, which are found at one end of each chromosome. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60-Mb O. tipulae genome is resolved into six chromosomal molecules. We find the evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes, we identify seven ancestral chromosomal elements (Nigon elements) and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex chromosome-associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

Two novel loci underlie natural differences in Caenorhabditis elegans abamectin responses

Parasitic nematodes cause a massive worldwide burden on human health along with a loss of livestock and agriculture productivity. Anthelmintics have been widely successful in treating parasitic nematodes. However, resistance is increasing, and little is known about the molecular and genetic causes of resistance for most of these drugs. The free-living roundworm Caenorhabditis elegans provides a tractable model to identify genes that underlie resistance. Unlike parasitic nematodes, C. elegans is easy to maintain in the laboratory, has a complete and well annotated genome, and has many genetic tools. Using a combination of wild isolates and a panel of recombinant inbred lines constructed from crosses of two genetically and phenotypically divergent strains, we identified three genomic regions on chromosome V that underlie natural differences in response to the macrocyclic lactone (ML) abamectin. One locus was identified previously and encodes an alpha subunit of a glutamate-gated chloride channel (glc-1). Here, we validate and narrow two novel loci using near-isogenic lines. Additionally, we generate a list of prioritized candidate genes identified in C. elegans and in the parasite Haemonchus contortus by comparison of ML resistance loci. These genes could represent previously unidentified resistance genes shared across nematode species and should be evaluated in the future. Our work highlights the advantages of using C. elegans as a model to better understand ML resistance in parasitic nematodes.

Complementary Approaches with Free-living and Parasitic Nematodes to Understanding Anthelmintic Resistance

Anthelmintic drugs are the major line of defense against parasitic nematode infections, but the arsenal is limited and resistance threatens sustained efficacy of the available drugs. Discoveries of the modes of action of these drugs and mechanisms of resistance have predominantly come from studies of a related nonparasitic nematode species, Caenorhabditis elegans, and the parasitic nematode Haemonchus contortus. Here, we discuss how our understanding of anthelmintic resistance and modes of action came from the interplay of results from each of these species. We argue that this 'cycle of discovery', where results from one species inform the design of experiments in the other, can use the complementary strengths of both to understand anthelmintic modes of action and mechanisms of resistance.

Challenges and opportunities for the adoption of molecular diagnostics for anthelmintic resistance

Anthelmintic resistance is a significant threat to livestock production systems worldwide and is emerging as an important issue in companion animal parasite management. It is also an emerging concern for the control of human soil-transmitted helminths and filaria. An important aspect of managing anthelmintic resistance is the ability to utilise diagnostic tests to detect its emergence at an early stage. In host-parasite systems where resistance is already widespread, diagnostics have a potentially important role in determining those drugs that remain the most effective. The development of molecular diagnostics for anthelmintic resistance is one focus of the Consortium for Anthelmintic Resistance and Susceptibility (CARS) group. The present paper reflects discussions of this issue that occurred at the most recent meeting of the group in Wisconsin, USA, in July 2019. We compare molecular resistance diagnostics with in vivo and in vitro phenotypic methods, and highlight the advantages and disadvantages of each. We assess whether our knowledge on the identity of molecular markers for resistance towards the different drug classes is sufficient to provide some expectation that molecular tests for field use may be available in the short-to-medium term. We describe some practical aspects of such tests and how our current capabilities compare to the requirements of an 'ideal' test. Finally, we describe examples of drug class/parasite species interactions that provide the best opportunity for commercial use of molecular tests in the near future. We argue that while such prototype tests may not satisfy the requirements of an 'ideal' test, their potential to provide significant advances over currently-used phenotypic methods warrants their development as field diagnostics.

Ivermectin: An Anthelmintic, an Insecticide, and Much More

Here we tell the story of ivermectin, describing its anthelmintic and insecticidal actions and recent studies that have sought to reposition ivermectin for the treatment of other diseases that are not caused by helminth and insect parasites. The standard theory of its anthelmintic and insecticidal mode of action is that it is a selective positive allosteric modulator of glutamate-gated chloride channels found in nematodes and insects. At higher concentrations, ivermectin also acts as an allosteric modulator of ion channels found in host central nervous systems. In addition, in tissue culture, at concentrations higher than anthelmintic concentrations, ivermectin shows antiviral, antimalarial, antimetabolic, and anticancer effects. Caution is required before extrapolating from these preliminary repositioning experiments to clinical use, particularly for Covid-19 treatment, because of the high concentrations of ivermectin used in tissue-culture experiments.

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

Haemonchus contortus is a globally distributed and economically important gastrointestinal pathogen of small ruminants and has become a key nematode model for studying anthelmintic resistance and other parasite-specific traits among a wider group of parasites including major human pathogens. Here, we report using PacBio long-read and OpGen and 10X Genomics long-molecule methods to generate a highly contiguous 283.4 Mbp chromosome-scale genome assembly including a resolved sex chromosome for the MHco3(ISE).N1 isolate. We show a remarkable pattern of conservation of chromosome content with Caenorhabditis elegans, but almost no conservation of gene order. Short and long-read transcriptome sequencing allowed us to define coordinated transcriptional regulation throughout the parasite's life cycle and refine our understanding of cis- and trans-splicing. Finally, we provide a comprehensive picture of chromosome-wide genetic diversity both within a single isolate and globally. These data provide a high-quality comparison for understanding the evolution and genomics of Caenorhabditis and other nematodes and extend the experimental tractability of this model parasitic nematode in understanding helminth biology, drug discovery and vaccine development, as well as important adaptive traits such as drug resistance.

The use of high resolution melting analysis of ITS-1 for rapid differentiation of parasitic nematodes Haemonchus contortus and Ashworthius sidemi

Among gastrointestinal nematodes, haematophagous strongylids Haemonchus contortus and Ashworthius sidemi belong to the most pathogenic parasites of both domestic and wild ruminants. Correct identification of parasitic taxa is of crucial importance in many areas of parasite research, including monitoring of occurrence, epidemiological studies, or testing of effectiveness of therapy. In this study, we identified H. contortus and A. sidemi in a broad range of ruminant hosts that occur in the Czech Republic using morphological/morphometric and molecular approaches. As an advanced molecular method, we employed qPCR followed by High Resolution Melting analysis, specifically targeting the internal transcribed spacer 1 (ITS-1) sequence to distinguish the two nematode species. We demonstrate that High Resolution Melting curves allow for taxonomic affiliation, making it a convenient, rapid, and reliable identification tool.

Toward integrative 'omics of the barber's pole worm and related parasitic nematodes

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

Moxidectin: an oral treatment for human onchocerciasis

INTRODUCTION

Moxidectin is a milbemycin endectocide recently approved for the treatment of human onchocerciasis. Onchocerciasis, earmarked for elimination of transmission, is a filarial infection endemic in Africa, Yemen, and the Amazonian focus straddling Venezuela and Brazil. Concerns over whether the predominant treatment strategy (yearly mass drug administration (MDA) of ivermectin) is sufficient to achieve elimination in all endemic foci have refocussed attention upon alternative treatments. Moxidectin's stronger and longer microfilarial suppression compared to ivermectin in both phase II and III clinical trials indicates its potential as a novel powerful drug for onchocerciasis elimination.

AREAS COVERED

This work summarizes the chemistry and pharmacology of moxidectin, reviews the phase II and III clinical trials evidence on tolerability, safety, and efficacy of moxidectin versus ivermectin, and discusses the implications of moxidectin's current regulatory status.

EXPERT OPINION

Moxidectin's superior clinical performance has the potential to substantially reduce times to elimination compared to ivermectin. If donated, moxidectin could mitigate the additional programmatic costs of biannual ivermectin distribution because, unlike other alternatives, it can use the existing community-directed treatment infrastructure. A pediatric indication (for children <12 years) and determination of its usefulness in onchocerciasis-loiasis co-endemic areas will greatly help fulfill the potential of moxidectin for the treatment and elimination of onchocerciasis.

Post-genomic progress in helminth parasitology

Helminth parasitology is an important discipline, which poses often unique technical challenges. One challenge is that helminth parasites, particularly those in humans, are often difficult to obtain alive and in sufficient quantities for study; another is the challenge of studying these organisms in vitro – no helminth parasite life cycle has been fully recapitulated outside of a host. Arguably, the key issue retarding progress in helminth parasitology has been a lack of experimental tools and resources, certainly relative to the riches that have driven many parasitologists to adopt free-living model organisms as surrogate systems. In response to these needs, the past 10–12 years have seen the beginnings of helminth parasitology's journey into the ‘omics’ era, with the release of abundant sequencing resources, and the functional genomics tools with which to test biological hypotheses. To reflect this progress, the 2019 Autumn Symposium of the British Society for Parasitology was held in Queen's University Belfast on the topic of ‘post-genomic progress in helminth parasitology’. This issue presents examples of the current state of play in the field, while this editorial summarizes how genomic datasets and functional genomic tools have stimulated impressive recent progress in our understanding of parasite biology.

A Whole Genome Re-Sequencing Based GWA Analysis Reveals Candidate Genes Associated with Ivermectin Resistance in Haemonchus contortus

The most important and broad-spectrum drug used to control the parasitic worms to date is ivermectin (IVM). Resistance against IVM has emerged in parasites, and preserving its efficacy is now becoming a serious issue. The parasitic nematode Haemonchus contortus (Rudolphi, 1803) is economically an important parasite of small ruminants across the globe, which has a successful track record in IVM resistance. There are growing evidences regarding the multigenic nature of IVM resistance, and although some genes have been proposed as candidates of IVM resistance using lower magnification of genome, the genetic basis of IVM resistance still remains poorly resolved. Using the full magnification of genome, we herein applied a population genomics approach to characterize genome-wide signatures of selection among pooled worms from two susceptible and six ivermectin-resistant isolates of H. contortus, and revealed candidate genes under selection in relation to IVM resistance. These candidates also included a previously known IVM-resistance-associated candidate gene HCON_00148840, glc-3. Finally, an RNA-interference-based functional validation assay revealed the HCON_00143950 as IVM-tolerance-associated gene in H. contortus. The possible role of this gene in IVM resistance could be detoxification of xenobiotic in phase I of xenobiotic metabolism. The results of this study further enhance our understanding on the IVM resistance and continue to provide further evidence in favor of multigenic nature of IVM resistance.

Genomic Epidemiology in Filarial Nematodes: Transforming the Basis for Elimination Program Decisions

Onchocerciasis and lymphatic filariasis are targeted for elimination, primarily using mass drug administration at the country and community levels. Elimination of transmission is the onchocerciasis target and global elimination as a public health problem is the end point for lymphatic filariasis. Where program duration, treatment coverage, and compliance are sufficiently high, elimination is achievable for both parasites within defined geographic areas. However, transmission has re-emerged after apparent elimination in some areas, and in others has continued despite years of mass drug treatment. A critical question is whether this re-emergence and/or persistence of transmission is due to persistence of local parasites-i.e., the result of insufficient duration or drug coverage, poor parasite response to the drugs, or inadequate methods of assessment and/or criteria for determining when to stop treatment-or due to re-introduction of parasites via human or vector movement from another endemic area. We review recent genetics-based research exploring these questions in Onchocerca volvulus, the filarial nematode that causes onchocerciasis, and Wuchereria bancrofti, the major pathogen for lymphatic filariasis. We focus in particular on the combination of genomic epidemiology and genome-wide associations to delineate transmission zones and distinguish between local and introduced parasites as the source of resurgence or continuing transmission, and to identify genetic markers associated with parasite response to chemotherapy. Our ultimate goal is to assist elimination efforts by developing easy-to-use tools that incorporate genetic information about transmission and drug response for more effective mass drug distribution, surveillance strategies, and decisions on when to stop interventions to improve sustainability of elimination.

Zoonotic transmission of intestinal helminths in southeast Asia: Implications for control and elimination

Intestinal helminths are extremely widespread and highly prevalent infections of humans, particularly in rural and poor urban areas of low and middle-income countries. These parasites have chronic and often insidious effects on human health and child development including abdominal problems, anaemia, stunting and wasting. Certain animals play a fundamental role in the transmission of many intestinal helminths to humans. However, the contribution of zoonotic transmission to the overall burden of human intestinal helminth infection and the relative importance of different animal reservoirs remains incomplete. Moreover, control programmes and transmission models for intestinal helminths often do not consider the role of zoonotic reservoirs of infection. Such reservoirs will become increasingly important as control is scaled up and there is a move towards interruption and even elimination of parasite transmission. With a focus on southeast Asia, and the Philippines in particular, this review summarises the major zoonotic intestinal helminths, risk factors for infection and highlights knowledge gaps related to their epidemiology and transmission. Various methodologies are discussed, including parasite genomics, mathematical modelling and socio-economic analysis, that could be employed to improve understanding of intestinal helminth spread, reservoir attribution and the burden associated with infection, as well as assess effectiveness of interventions. For sustainable control and ultimately elimination of intestinal helminths, there is a need to move beyond scheduled mass deworming and to consider animal and environmental reservoirs. A One Health approach to control of intestinal helminths is proposed, integrating interventions targeting humans, animals and the environment, including improved access to water, hygiene and sanitation. This will require coordination and collaboration across different sectors to achieve best health outcomes for all.

The global diversity of Haemonchus contortus is shaped by human intervention and climate

Haemonchus contortus is a haematophagous parasitic nematode of veterinary interest. We have performed a survey of its genome-wide diversity using single-worm whole genome sequencing of 223 individuals sampled from 19 isolates spanning five continents. We find an African origin for the species, together with evidence for parasites spreading during the transatlantic slave trade and colonisation of Australia. Strong selective sweeps surrounding the β-tubulin locus, a target of benzimidazole anthelmintic drug, are identified in independent populations. These sweeps are further supported by signals of diversifying selection enriched in genes involved in response to drugs and other anthelmintic-associated biological functions. We also identify some candidate genes that may play a role in ivermectin resistance. Finally, genetic signatures of climate-driven adaptation are described, revealing a gene acting as an epigenetic regulator and components of the dauer pathway. These results begin to define genetic adaptation to climate in a parasitic nematode.

Based on single worm whole genome sequencing, the authors here characterise the global evolution of the gastrointestinal parasite Haemonchus contortus and identify genes that play a role in drug resistance as well as climate-driven adaptations involving an epigenetic regulator.

The confounding effects of high genetic diversity on the determination and interpretation of differential gene expression analysis in the parasitic nematode Haemonchus contortus

Differential expression analysis between parasitic nematode strains is commonly used to implicate candidate genes in anthelmintic resistance or other biological functions. We have tested the hypothesis that the high genetic diversity of an organism such as Haemonchus contortus could complicate such analyses. First, we investigated the extent to which sequence polymorphism affects the reliability of differential expression analysis between the genetically divergent H. contortus strains MHco3(ISE), MHco4(WRS) and MHco10(CAVR). Using triplicates of 20 adult female worms from each population isolated under parallel experimental conditions, we found that high rates of sequence polymorphism in RNAseq reads were associated with lower efficiency read mapping to gene models under default TopHat2 parameters, leading to biased estimates of inter-strain differential expression. We then showed it is possible to largely compensate for this bias by optimising the read mapping single nucleotide polymorphism (SNP) allowance and filtering out genes with particularly high single nucleotide polymorphism rates. Once the sequence polymorphism biases were removed, we then assessed the genuine transcriptional diversity between the strains, finding ≥824 differentially expressed genes across all three pairwise strain comparisons. This high level of inter-strain transcriptional diversity not only suggests substantive inter-strain phenotypic variation but also highlights the difficulty in reliably associating differential expression of specific genes with phenotypic differences. To provide a practical example, we analysed two gene families of potential relevance to ivermectin drug resistance; the ABC transporters and the ligand-gated ion channels (LGICs). Over half of genes identified as differentially expressed using default TopHat2 parameters were shown to be an artifact of sequence polymorphism differences. This work illustrates the need to account for sequence polymorphism in differential expression analysis. It also demonstrates that a large number of genuine transcriptional differences can occur between H. contortus strains and these must be considered before associating the differential expression of specific genes with phenotypic differences between strains.

Extreme-QTL mapping of monepantel resistance in Haemonchus contortus

Background

Haemonchus contortus, a gastrointestinal nematode parasite of sheep, is mainly controlled by anthelmintics; the occurrence of anthelmintic resistance leads to treatment failures and increases economic burden. Because molecular mechanisms involved in drug resistance can be elucidated by genomic studies, an extreme quantitative trait locus (X-QTL) mapping approach was used to identify co-segregation of the resistance phenotype with genetic markers to detect the genome-wide variants associated with monepantel resistance in H. contortus.

Methods

A cross between H. contortus isolates using parental susceptible (Par-S) males and monepantel resistant (Par-R) females resulted in SR progeny, while reciprocal cross resulted in RS progeny. Pools (n = 30,000) of infective larvae (L3) recovered from Par-R, and from SR and RS populations in the F3 generation, collected both before (unselected group) and 7 days after (selected group) selection with monepantel treatment in sheep hosts, were subjected to genome sequencing (Pool-Seq). Pairwise comparisons of allele frequencies between unselected and selected groups were performed for each population by Fisher’s exact test (FET) and for both populations combined by a Cochran-Mantel-Haenszel (CMH) test.

Results

Mapping rates varied from 80.29 to 81.77% at a 90.4X mean coverage of aligned reads. After correction for multiple testing, significant (P < 0.05) changes in allele frequencies were detected by FET for 6 and 57 single nucleotide polymorphisms (SNPs) in the SR and RS populations, respectively, and by the CMH test for 124 SNPs in both populations. The significant variants located on chromosome 2 generated a selection signal in a genomic region harboring the mptl-1, deg-3 and des-2 genes, previously reported as candidates for monepantel resistance. In addition, three new variants were identified in the mptl-1 gene.

Conclusions

This study expands knowledge on genome-wide molecular events underlying H. contortus resistance to monepantel. The identification of a genome region harboring major genes previously associated with monepantel resistance supports the results of the employed X-QTL approach. In addition, a deletion in exon 11 of the mptl-1 gene should be further investigated as the putative causal mutation leading to monepantel resistance.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3663-9) contains supplementary material, which is available to authorized users.

Development of amplicon sequencing for the analysis of benzimidazole resistance allele frequencies in field populations of gastrointestinal nematodes

Anthelmintic resistant gastrointestinal helminths have become a major cause of poor health in sheep and goats. Sensitive and specific molecular markers are needed to monitor the genotypic frequency of resistance in field parasite populations. Gastrointestinal nematode resistance to benzimidazole is caused by a mutation in one of three positions within the isotype 1 β-tubulin gene. In the absence of markers for resistance to the other broad spectrum anthelmintic classes, these provide a relevant study example. Determination of the prevalence of these single nucleotide polymorphisms in field nematode populations can be impractical using conventional molecular methods to examine individual parasites; which can be laborious and lack sensitivity in determining low levels of resistance in parasite populations. Here, we report the development of a novel method based on an Illumina MiSeq deep amplicon sequencing platform to sequence the isotype 1 β-tubulin locus of the small ruminant gastrointestinal nematode, Teladorsagia circumcincta, and determine the frequency of the benzimidazole resistance mutations. We validated the method by assessing sequence representation bias, comparing the results of Illumina MiSeq and pyrosequencing, and applying the method to populations containing known proportions of resistant and susceptible larvae. We applied the method to field samples collected from ewes and lambs on over a period of one year on three farms, each highlighting different aspects of sheep management and approaches to parasite control. The results show opportunities to build hypotheses with reference to selection pressures leading to differences in resistance allele frequencies between sampling dates, farms and ewes or lambs, and to consider the impact of their genetic fixation or otherwise. This study provides proof of concept of a practical, accurate, sensitive and scalable method to determine frequency of anthelmintic resistance mutations in gastrointestinal nematodes in field studies and as a management tool for livestock farmers.

Current therapeutic applications and pharmacokinetic modulations of ivermectin

Ivermectin is considered to be a wonder drug due to its broad-spectrum antiparasitic activity against both ectoparasites and endoparasites (under class of endectocide) and has multiple applications in both veterinary and human medicine. In particular, ivermectin is commonly used in the treatment of different kinds of infections and infestations. By altering the vehicles used in the formulations, the pharmacokinetic properties of different ivermectin preparations can be altered. Since its development, various vehicles have been evaluated to assess the efficacy, safety, and therapeutic systemic concentrations of ivermectin in different species. A subcutaneous route of administration is preferred over a topical or an oral route for ivermectin due to superior bioavailability. Different formulations of ivermectin have been developed over the years, such as stabilized aqueous formulations, osmotic pumps, controlled release capsules, silicone carriers, zein microspheres, biodegradable microparticulate drug delivery systems, lipid nanocapsules, solid lipid nanoparticles, sustained-release ivermectin varnish, sustained-release ivermectin-loaded solid dispersion suspension, and biodegradable subcutaneous implants. However, several reports of ivermectin resistance have been identified in different parts of the world over the past few years. Continuous use of suboptimal formulations or sub-therapeutic plasma concentrations may predispose an individual to resistance toward ivermectin. The current research trend is focused toward the need for developing ivermectin formulations that are stable, effective, and safe and that reduce the number of doses required for complete clinical cure in different parasitic diseases. Therefore, single-dose long-acting preparations of ivermectin that provide effective therapeutic drug concentrations need to be developed and commercialized, which may revolutionize drug therapy and prophylaxis against various parasitic diseases in the near future. The present review highlights the current advances in pharmacokinetic modulation of ivermectin formulations and their potent therapeutic applications, issues related to emergence of ivermectin resistance, and future trends of ivermectin usage.