Beta-tubulin genes from the parasitic nematode Haemonchus contortus modulate drug resistance in Caenorhabditis elegans

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.

Quantitative tests of albendazole resistance in beta-tubulin mutants

Benzimidazole (BZ) anthelmintics are among the most important treatments for parasitic nematode infections in the developing world. Widespread BZ resistance in veterinary parasites and emerging resistance in human parasites raise major concerns for the continued use of BZs. Knowledge of the mechanisms of resistance is necessary to make informed treatment decisions and circumvent resistance. Benzimidazole resistance has traditionally been associated with mutations and natural variants in the C. elegans beta-tubulin gene ben-1 and orthologs in parasitic species. However, variants in ben-1 alone do not explain the differences in BZ responses across parasite populations. Here, we examine the roles of five C. elegans beta-tubulin genes (tbb-1, mec-7, tbb-4, ben-1, and tbb-6) to identify the role each gene plays in BZ response. We generated C. elegans strains with a loss of each beta-tubulin gene, as well as strains with a loss of tbb-1, mec-7, tbb-4, or tbb-6 in a genetic background that also lacks ben-1 to test beta-tubulin redundancy in BZ response. We found that only the individual loss of ben-1 conferred a substantial level of BZ resistance, although the loss of tbb-1 was found to confer a small benefit in the presence of albendazole (ABZ). The loss of ben-1 was found to confer an almost complete rescue of animal development in the presence of 30 μM ABZ, likely explaining why no additive effects caused by the loss of a second beta-tubulin were observed. We demonstrate that ben-1 is the only beta-tubulin gene in C. elegans where loss confers substantial BZ resistance.

Assessing Environmental Risks during the Drug Development Process for Parasitic Vector-Borne Diseases: A Critical Reflection

Parasitic vector-borne diseases (VBDs) represent nearly 20% of the global burden of infectious diseases. Moreover, the spread of VBDs is enhanced by global travel, urbanization, and climate change. Treatment of VBDs faces challenges due to limitations of existing drugs, as the potential for side effects in nontarget species raises significant environmental concerns. Consequently, considering environmental risks early in drug development processes is critically important. Here, we examine the environmental risk assessment process for veterinary medicinal products in the European Union and identify major gaps in the ecotoxicity data of these drugs. By highlighting the scarcity of ecotoxicological data for commonly used antiparasitic drugs, we stress the urgent need for considering the One Health concept. We advocate for employing predictive tools and nonanimal methodologies such as New Approach Methodologies at early stages of antiparasitic drug research and development. Furthermore, adopting progressive approaches to mitigate ecological risks requires the integration of nonstandard tests that account for real-world complexities and use environmentally relevant exposure scenarios. Such a strategy is vital for a sustainable drug development process as it adheres to the principles of One Health, ultimately contributing to a healthier and more sustainable world.

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.

Model nematodes as a practical innovation to promote high throughput screening of natural products for anthelmintics discovery in South Asia: Current challenges, proposed practical and conceptual solutions

With the increasing prevalence of anthelmintic resistance in animals recorded globally, and the threat of resistance in human helminths, the need for novel anthelmintic drugs is greater than ever. Most research aimed at discovering novel anthelmintic leads relies on high throughput screening (HTS) of large libraries of synthetic small molecules in industrial and academic settings in developed countries, even though it is the tropical countries that are most plagued by helminth infections. Tropical countries, however, have the advantage of possessing a rich flora that may yield natural products (NP) with promising anthelmintic activity. Focusing on South Asia, which produces one of the world's highest research outputs in NP and NP-based anthelmintic discovery, we find that limited basic research and funding, a lack of awareness of the utility of model organisms, poor industry-academia partnerships and lack of technological innovations greatly limit anthelmintics research in the region. Here we propose that utilizing model organisms including the free-living nematode Caenorhabditis elegans, that can potentially allow rapid target identification of novel anthelmintics, and Oscheius tipulae, a closely related, free-living nematode which is found abundantly in soil in hotter temperatures, could be a much-needed innovation that can enable cost-effective and efficient HTS of NPs for discovering compounds with anthelmintic/antiparasitic potential in South Asia and other tropical regions that historically have devoted limited funding for such research. Additionally, increased collaborations at the national, regional and international level between parasitologists and pharmacologists/ethnobotanists, setting up government-industry-academia partnerships to fund academic research, creating a centralized, regional collection of plant extracts or purified NPs as a dereplication strategy and HTS library, and holding regional C. elegans/O. tipulae-based anthelmintics workshops and conferences to share knowledge and resources regarding model organisms may collectively promote and foster a NP-based anthelmintics landscape in South Asia and beyond.

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.

Detection of SNPs and benzimidazole resistance in strongyle nematode eggs of horses by allele-specific PCR

This study was conducted to determine single nucleotide polymorphisms (SNPs) and the benzimidazole (BZ) resistance in strongyle nematode egg populations in horses using molecular techniques. A total of 200 fecal samples were collected from horses in 26 farms in two provinces (Kayseri and Nevşehir) of the Central Anatolia Region of Türkiye between May and August 2022. The flotation method was used to detect strongyle nematode eggs in the fecal samples of the horses. Afterward, strongyle nematode eggs were collected, and the allele-specific polymerase chain reaction (AS-PCR) technique was used to detect the BZ resistance. BZ-susceptible and BZ-resistant PCR products were sequenced to determine single nucleotide polymorphisms (SNPs) in the β-tubulin isotype 1 gene. The strongyle nematode eggs were determined in 85 (42.5%) out of 200 fecal samples. AS-PCR detected 50.58% (43/85) BZ-resistant (homozygous resistant) and 36.4% (31/85) BZ-susceptible (homozygous susceptible) genes in the strongyle eggs. Both BZ-resistant and BZ-susceptible genes (heterozygous) were determined in 11 samples. BZ-resistant and BZ-susceptible allele frequencies were determined as 57.0% (48.5/85) and 43.0% (36.5/85), respectively. SNPs were detected only in codon 200 of the β-tubulin isotype 1 gene in four sequenced isolates of the two resistant and two susceptible isolates. This study is the first molecular report on BZ resistance in strongyle nematode eggs in horses in Türkiye. The widespread prevalence of BZ-resistant alleles in equine strongyle nematodes shows the requirement for the immediate usage of other anthelmintics instead of the BZ group drugs for the effective management and control of equine strongyle nematodes.

Prevalence of gastrointestinal parasites and molecular identification of beta-tubulin mutations associated with benzimidazole resistance in Haemonchus contortus in goats from selected districts of Uganda

Gastrointestinal parasites are among the most economically important pathogens of small ruminants causing serious economic losses and animal welfare problems for the livestock industry worldwide. The emergence of anthelmintic resistant H. contortus in small ruminants is a serious problem because it undermines effective helminth control and results in reduced productivity. Little is known about resistance to benzimidazoles (BZ) in Haemonchus in goats and sheep in Uganda. The objective of this study was to determine the prevalence of gastrointestinal parasites and to identify the presence of benzimidazole resistance associated mutations in the β-tubulin isotype 1 gene of Haemonchus contortus in goats from selected districts of Uganda. A total of 200 goats from 10 districts of Uganda slaughtered at Kalerwe abattoir in Kampala were sampled for H. contortus adult worms. Faecal samples were also collected to detect other intestinal parasites. Faecal microscopy and analysis were performed using flotation and sedimentation techniques. DNA was extracted from adult worms and PCR and sequencing of the ITS- 2 region and β-tubulin isotype 1 gene performed to identify H. contortus species and to determine the presence of mutations associated with anthelmintic resistance respectively. Faecal microscopy showed that the most prevalent intestinal parasites were coccidia (98%), strongyles (97.5%), Strongyloides (82%), Paramphistomum (74.5%), Moniezia (46%), Fasciola (1.5%) and Trichuris (1%). Most goats had a high intestinal burden of coccidia (≥ 5000 oocyst per gram) and strongyles (≥ 1000 egg per gram), 65% and 67.5%, respectively. The prevalence of H. contortus adult worms was 63% (126/200). Sequencing of the partial β-tubulin isotype 1 gene of 54 Haemonchus contortus adult male isolates revealed the presence of mutations associated with anthelmintic resistance. The F200Y mutation was the most common mutation (13% of samples with good beta-tubulin sequences) followed by the E198A and E198K mutations, both found in 9% of sequenced samples. Mutation F167Y was not identified in any of the samples and there were no heterozygous individuals for any of the SNPS associated with BZ resistance that were identified in this study. These findings highlight the need for controlled use of anthelmintics especially benzimidazoles, to enable sustainable control of H. contortus in Uganda, and a need for further investigation to understand the resistance of other parasites identified in this study.

Development of a three-colour digital PCR for early and quantitative detection of benzimidazole resistance-associated single nucleotide polymorphisms in Haemonchus contortus

Haemonchus contortus is the most pathogenic nematode in small ruminants and anthelmintic resistance (AR) hampers its efficient control. Early detection of AR status is required to reduce selection for AR and cannot be achieved using phenotypic tests. For benzimidazoles (BZs), the detection of AR-associated alleles characterised by single nucleotide polymorphisms (SNPs) in the isotype 1 β-tubulin gene allows early AR detection in strongyles. The F200Y, F167Y, E198A and E198L polymorphisms have been described in BZ-resistant populations with a clear variation in frequencies between regions. A novel digital PCR (dPCR) enables the detection of all of the above-described polymorphisms in H. contortus. Assays were validated using synthetic DNA fragments containing these SNPs. Then, larvae obtained and pooled at farm level from 26 Austrian and 10 Italian sheep farms were analysed. For all assays a detection limit of 15 copies/μl of resistance alleles and a high level of accuracy were demonstrated, allowing to detect allele frequencies of 1% in most samples. In Austrian samples, elevated frequencies of F200Y resistance alleles were detected on all farms. Polymorphisms in codon 167 and codon 198 were identified in H. contortus from Austria for the first time. In Italian samples, the frequency of resistance alleles was still comparatively low, but F200Y resistance alleles were traceable. In conclusion we developed for the first time dPCR assays that target all SNPs of relevance associated with BZ-resistance in H. contortus. Future research on AR development could benefit from an early onset of SNP-based surveillance that would include the developed assays for all SNPs of relevance. Improved surveillance in the long term will include other important, though less pathogenic, nematode genera in the analyses.

Revealing novel cytb and nad5 genes-based population diversity and benzimidazole resistance in Echinococcus granulosus of bovine origin

Cystic echinococcosis (CE) is a neglected zoonotic disease caused by Echinococcus granulosus (sensu stricto). The parasite affects a wide range of livestock and wild animals. In this study, the population diversity of the Echinococcus species was investigated based on mitochondrial cytochrome b (cytb) and NADH dehydrogenase subunit 5 (nad5) genes. In addition to this, β-tubulin gene isoforms of Echinococcus granulosus were amplified to determine the resistance against benzimidazoles. For this purpose, 40 cyst samples from cattle (n = 20) and buffaloes (n = 20) were collected from the main abattoir of Sialkot. DNA extraction was performed using Qiagen Blood and Tissue Kits. Amplification was performed through PCR. Each amplicon was confirmed by GelRed™ stained agarose gel (2%). Samples were sequenced in a DNA analyzer and viewed for any misread nucleotide by using MEGA (v.11). Corrections in nucleotide sequence and multiple sequence alignment were made through the same software. NCBI-BLAST was used for sample specific sequences to identify them as belonging to a particular species. Diversity indices were estimated using DnaSP (v.6) while phylogenetic analysis was inferred using the Bayesian method using MrBayes (v.1.1). β-tubulin gene isoforms sequence analysis was performed to find out the candidate gene causing benzimidazole resistance. All 40 isolates were found positive for E. granulosus. BLAST-based searches of sequences of each isolate for each gene (nad5 and cytb) confirmed their maximum similarity with the G1 genotype. Overall, high haplotype diversity (Hd nad5 = 1.00; Hd cytb = 0.833) and low nucleotide diversity (π nad5 = 0.00560; π = cytb = 0.00763) was identified based on diversity indices. For both the genes, non-significant values of Tajima's D (nad5 = -0.81734; cytb = -0.80861) and Fu's Fs (nad5 = -1.012; cytb = 0.731) indicate recent population expansion. Bayesian phylogeny-based results of nad5 and cytb sequences confirmed their genotypic status as distinct from other Echinococcus species. This study shed light on the status of benzimidazole resistance in Echinococcus granulosus for the very first time from Pakistan. The findings of this study will significantly add in the information available on genetic diversity of Echinoccous granulosus based on cytb and nad5 genes sequences.

Functional validation of the truncated UNC-63 acetylcholine receptor subunit in levamisole resistance

Levamisole is a broad-spectrum anthelmintic which permanently activates cholinergic receptors from nematodes, inducing a spastic paralysis of the worms. Whereas this molecule is widely used to control parasitic nematodes impacting livestock, its efficacy is compromised by the emergence of drug-resistant parasites. In that respect, there is an urgent need to identify and validate molecular markers associated with resistance. Previous transcriptomic analyses revealed truncated cholinergic receptor subunits as potential levamisole resistance markers in the trichostrongylid nematodes Haemonchus contortus, Telodorsagia circumcincta and Trichostrongylus colubriformis. In the present study we used the Xenopus oocyte, as well as the free-living model nematode Caenorhabditis elegans, as heterologous expression systems to functionally investigate truncated isoforms of the levamisole-sensitive acetylcholine receptor (L-AChR) UNC-63 subunit. In the Xenopus oocyte, we report that truncated UNC-63 from C. elegans has a strong dominant negative effect on the expression of the recombinant C. elegans L-AChRs. In addition, we show that when expressed in C. elegans muscle cells, truncated UNC-63 induces a drastic reduction in levamisole susceptibility in transgenic worms, thus providing the first known functional validation for this molecular marker in vivo.

Pharmacological Profiling of a Brugia malayi Muscarinic Acetylcholine Receptor as a Putative Antiparasitic Target

The diversification of anthelmintic targets and mechanisms of action will help ensure the sustainable control of nematode infections in response to the growing threat of drug resistance. G protein-coupled receptors (GPCRs) are established drug targets in human medicine but remain unexploited as anthelmintic substrates despite their important roles in nematode neuromuscular and physiological processes. Bottlenecks in exploring the druggability of parasitic nematode GPCRs include a limited helminth genetic toolkit and difficulties establishing functional heterologous expression. In an effort to address some of these challenges, we profile the function and pharmacology of muscarinic acetylcholine receptors in the human parasite Brugia malayi, an etiological agent of human lymphatic filariasis. While acetylcholine-gated ion channels are intensely studied as targets of existing anthelmintics, comparatively little is known about metabotropic receptor contributions to parasite cholinergic signaling. Using multivariate phenotypic assays in microfilariae and adults, we show that nicotinic and muscarinic compounds disparately affect parasite fitness traits. We identify a putative G protein-linked acetylcholine receptor of B. malayi (Bma-GAR-3) that is highly expressed across intramammalian life stages and adapt spatial RNA in situ hybridization to map receptor transcripts to critical parasite tissues. Tissue-specific expression of Bma-gar-3 in Caenorhabditis elegans (body wall muscle, sensory neurons, and pharynx) enabled receptor deorphanization and pharmacological profiling in a nematode physiological context. Finally, we developed an image-based feeding assay as a reporter of pharyngeal activity to facilitate GPCR screening in parasitized strains. We expect that these receptor characterization approaches and improved knowledge of GARs as putative drug targets will further advance the study of GPCR biology across medically important nematodes.

Anthelmintic Agents from African Medicinal Plants: Review and Prospects

Soil-transmitted helminthiasis affects more than 1.5 billion people globally and largely remains a sanitary problem in Africa. These infections place a huge economic burden on poor countries and affect livestock production, causing substantial economic losses and poor animal health. The emergence of anthelmintic resistance, especially in livestock, and the potential for its widespread in humans create a need for the development of alternative therapies. Medicinal plants play a significant role in the management of parasitic diseases in humans and livestock, especially in Africa. This report reviews anthelmintic studies that have been conducted on medicinal plants growing in Africa and published within the past two decades. A search was made in various electronic databases, and only full articles in English were included in the review. Reports show that aqueous and hydroalcoholic extracts and polar fractions obtained from these crude extracts form the predominant (80%) form of the extracts studied. Medicinal plants, extracts, and compounds with different chemical groups have been studied for their anthelmintic potential. Polyphenols and terpenoids are the most reported groups. More than 64% of the studies employed in vitro assays against parasitic and nonparasitic nematode models. Egg hatch inhibition, larval migration inhibition, and paralysis are the common parameters assessed in vitro. About 72% of in vivo models involved small ruminants, 15% rodents, and 5% chicken. Egg and worm burden are the main factors assessed in vivo. There were no reports on interventions in humans cited within the period under consideration. Also, few reports have investigated the potential of combining plant extracts with common anthelmintic drugs. This review reveals the huge potential of African medicinal plants as sources of anthelmintic agents and the dire need for in-depth clinical studies of extracts, fractions, and compounds from African plants as anthelmintic agents in livestock, companion animals, and humans.

A Subset Screen of the Compounds Australia Scaffold Library Identifies 7-Acylaminodibenzoxazepinones as Potent and Selective Hits for Anti-Giardia Drug Discovery

On an annual basis the flagellate protozoan, Giardia duodenalis, is responsible for an estimated one billion human infections of which approximately two hundred million cause disease. However, the treatment of Giardia infections is reliant on a small group of chemotherapeutic classes that have a broad spectrum of antimicrobial activity and increasing treatment failure rates. To improve this situation, we need new drugs. In this study we screened the Compounds Australia Scaffolds Library for compounds with potent and selective activity against these parasites. Unlike previous drug discovery efforts that have focused on drug repurposing, this library is comprised of commercially available synthetic compounds arranged into lead-like scaffolds to facilitate structure activity relationship assessments and de novo drug discovery. A screen of 2451 compounds in this library identified 40 hits (>50% inhibitory activity at 10 µM, over 48 h). Secondary testing identified three compounds with IC50 values <1 μM and >50-fold selectivity for parasites over mammalian cells and a hit series, CL9406, comprising compounds with potent (lowest IC50 180 nM) and selective activity for Giardia parasites. The most promising compound in this series, SN00797640, displayed selective activity against assemblage A, B, and metronidazole resistant parasites which was parasiticidal (minimum lethal concentration 625 nM) and synergistic with albendazole. SN00797640 was well-tolerated when administered to mice at doses of 50 mg/kg daily for three days paving the way for pre-clinical in vivo activity assessment.

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.

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.

Advances in diagnosis of gastrointestinal nematodes in livestock and companion animals

Diagnosis of gastrointestinal nematodes in livestock and companion animals has been neglected for years and there has been an historical underinvestment in the development and improvement of diagnostic tools, undermining the undoubted utility of surveillance and control programmes. However, a new impetus by the scientific community and the quickening pace of technological innovations, are promoting a renaissance of interest in developing diagnostic capacity for nematode infections in veterinary parasitology. A cross-cutting priority for diagnostic tools is the development of pen-side tests and associated decision support tools that rapidly inform on the levels of infection and morbidity. This includes development of scalable, parasite detection using artificial intelligence for automated counting of parasitic elements and research towards establishing biomarkers using innovative molecular and proteomic methods. The aim of this review is to assess the state-of-the-art in the diagnosis of helminth infections in livestock and companion animals and presents the current advances of diagnostic methods for intestinal parasites harnessing (i) automated methods for copromicroscopy based on artificial intelligence, (ii) immunodiagnosis, and (iii) molecular- and proteome-based approaches. Regardless of the method used, multiple factors need to be considered before diagnostics test results can be interpreted in terms of control decisions. Guidelines on how to apply diagnostics and how to interpret test results in different animal species are increasingly requested and some were recently made available in veterinary parasitology for the different domestic species.

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.

Interactions of C. elegans β-tubulins with the microtubule inhibitor and anthelmintic drug albendazole

Parasitic nematodes are major human and agricultural pests, and benzimidazoles are amongst the most important broad spectrum anthelmintic drug class used for their control. Benzimidazole resistance is now widespread in many species of parasitic nematodes in livestock globally and an emerging concern for the sustainable control of human soil transmitted helminths. β-tubulin is the major benzimidazole target, although other genes may influence resistance. Among the six C. elegans β-tubulin genes, loss of ben-1 causes resistance without other apparent defects. Here, we explored the genetics of C. elegans β-tubulin genes in relation to the response to the benzimidazole derivative albendazole. The most highly expressed β-tubulin isotypes, encoded by tbb-1 and tbb-2, were known to be redundant with each other for viability, and their products are predicted not to bind benzimidazoles. We found that tbb-2 mutants, and to a lesser extent tbb-1 mutants, were hypersensitive to albendazole. The double mutant tbb-2 ben-1 is uncoordinated and short, resembling the wild type exposed to albendazole, but the tbb-1 ben-1 double mutant did not show the same phenotypes. These results suggest that tbb-2 is a modifier of ABZ sensitivity. To better understand how BEN-1 mutates to cause benzimidazole resistance, we isolated mutants resistant to albendazole and found that 15 of 16 mutations occurred in the ben-1 coding region. Mutations ranged from likely nulls to hypomorphs, and several corresponded to residues that cause resistance in other organisms. Null alleles of ben-1 are albendazole-resistant and BEN-1 shows high sequence identity with tubulins from other organisms, suggesting that many amino acid changes could cause resistance. However, our results suggest that missense mutations conferring resistance are not evenly distributed across all possible conserved sites. Independent of their roles in benzimidazole resistance, tbb-1 and tbb-2 may have specialized functions as null mutants of tbb-1 or tbb-2 were cold or heat sensitive, respectively.

Absence of Polymorphisms in Codons 167, 198 and 200 of All Seven β-tubulin Isotypes of Benzimidazole Susceptible and Resistant Parascaris spp. Specimens from Australia

Benzimidazoles resistance is widespread in strongyle parasitic nematodes and associated with polym orphisms in the codons 167, 198 and 200 of isotype 1 β-tubulin (tbb-1). In ascarids, benzimidazole (BZ) resistance has rarely been reported and in none of these cases were any of these polymorphisms detected. Here, available genome and transcriptome data from WormBase ParaSite were used to compare the complete β-tubulin reservoirs of Parascaris univalens, Ascaris suum and Ascaris lumbricoides. Adult Parascaris spp. specimens collected in Australia from horses after BZ treatment (susceptible, n = 13) or surviving BZ treatment and collected after ivermectin treatment (resistant, n = 10) were genotyped regarding codons 167, 198 and 200 using Sanger sequencing. Phylogenetic analyses clearly showed that there are no one-to-one ascarid orthologs of strongyle tbb-1 genes. In the reference genomes, as well as phenotypically susceptible and resistant Parascaris spp. from Australia, six out of seven β-tubulin genes showed a BZ-susceptible genotype (F167, E198, F200). The only exception were the testis-specific β-tubulin D genes from all three ascarid species that encode tyrosine at codon 200. This was observed independently of the BZ-susceptibility phenotype of Parascaris spp. These data suggest that different mechanisms lead to BZ resistance in ascarid and strongyle nematodes.

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.

Effect of Clinically Used Microtubule Targeting Drugs on Viral Infection and Transport Function

Microtubule targeting agents (MTAs) have been exploited mainly as anti-cancer drugs because of their impact on cellular division and angiogenesis. Additionally, microtubules (MTs) are key structures for intracellular transport, which is frequently hijacked during viral infection. We have analyzed the antiviral activity of clinically used MTAs in the infection of DNA and RNA viruses, including SARS-CoV-2, to find that MT destabilizer agents show a higher impact than stabilizers in the viral infections tested, and FDA-approved anti-helminthic benzimidazoles were among the most active compounds. In order to understand the reasons for the observed antiviral activity, we studied the impact of these compounds in motor proteins-mediated intracellular transport. To do so, we used labeled peptide tools, finding that clinically available MTAs impaired the movement linked to MT motors in living cells. However, their effect on viral infection lacked a clear correlation to their effect in motor-mediated transport, denoting the complex use of the cytoskeleton by viruses. Finally, we further delved into the molecular mechanism of action of Mebendazole by combining biochemical and structural studies to obtain crystallographic high-resolution information of the Mebendazole-tubulin complex, which provided insights into the mechanisms of differential toxicity between helminths and mammalians.

Molecular modelling and de novo fragment-based design of potential inhibitors of beta-tubulin gene of Necator americanus from natural products

The emergence of drug resistance against the known hookworm drugs namely albendazole and mebendazole and their reduced efficacies necessitate the need for new drugs. Chemically diverse natural products present plausible templates to augment hookworm drug discovery. The present work utilized pharmacoinformatics techniques to predict African natural compounds ZINC95486082, ZINC95486052 and euphohelionon as potential inhibitory molecules of the hookworm Necator americanus β tubulin gene. A library of 3390 compounds was screened against a homology-modelled structure of β tubulin. The docking results obtained from AutoDock Vina was validated with an acceptable area under the curve (AUC) of 0.714 computed from the receiver operating characteristic (ROC) curve. The three selected compounds had favourable binding affinities and were predicted to form no interactions with the resistance-associated mutations Phe167, Glu198 and Phe200. The compounds were predicted as anthelmintics using a Bayesian-based technique and were pharmacologically profiled to be druglike. Further molecular dynamics simulations and MM-PBSA calculations showed the compounds as promising anthelmintic drug leads. Novel critical residues comprising Leu246, Asn247 and Asn256 were also predicted for binding. Euphohelionon was selected as a template for the de novo fragment-based design of five compounds labelled A1, A2, A3, A4 and A5; with four of them having SAscore values below 6, denoting easy synthesis. All the five de novo molecules docked firmly in the binding pocket of the β tubulin with no binding interactions with the three known resistance mutation residues. Binding energies of −8.2, −7.6, −7.3, −7.2 and −6.8 kcal/mol were obtained for A1, A2, A3, A4 and A5, respectively. The identified compounds can serve as treasure troves from which future potent anthelmintics can be designed. The current study strives to assuage the hookworm disease burden, especially making available molecules with the potential to circumvent the chemoresistance.

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.

Newly identified parasitic nematode beta-tubulin alleles confer resistance to benzimidazoles

Infections by parasitic nematodes cause large health and economic burdens worldwide. We use anthelmintic drugs to reduce these infections. However, resistance to anthelmintic drugs is extremely common and increasing worldwide. It is essential to understand the mechanisms of resistance to slow its spread. Recently, four new parasitic nematode beta-tubulin alleles have been identified in benzimidazole (BZ) resistant parasite populations: E198I, E198K, E198T, and E198stop. These alleles have not been tested for the ability to confer resistance or for any effects that they might have on organismal fitness. We introduced these four new alleles into the sensitive C. elegans laboratory-adapted N2 strain and exposed these genome-edited strains to both albendazole and fenbendazole. We found that all four alleles conferred resistance to both BZ drugs. Additionally, we tested for fitness consequences in both control and albendazole conditions over seven generations in competitive fitness assays. We found that none of the edited alleles had deleterious effects on fitness in control conditions and that all four alleles conferred strong and equivalent fitness benefits in BZ drug conditions. Because it is unknown if previously validated alleles confer a dominant or recessive BZ resistance phenotype, we tested the phenotypes caused by five of these alleles and found that none of them conferred a dominant BZ resistance phenotype. Accurate measurements of resistance, fitness effects, and dominance caused by the resistance alleles allow for the generation of better models of population dynamics and facilitate control practices that maximize the efficacy of this critical anthelmintic drug class.

•

Four newly identified parasitic nematode beta-tubulin alleles confer benzimidazole resistance.

•

The four newly identified alleles do not confer deleterious fitness consequences.

•

Five beta-tubulin alleles confer recessive benzimidazole resistance.

Development of a proteochemometric-based support vector machine model for predicting bioactive molecules of tubulin receptors

Microtubules are receiving enormous interest in drug discovery due to the important roles they play in cellular functions. Targeting tubulin polymerization presents an excellent opportunity for the development of anti-tubulin drugs. Drug resistance and high toxicity of currently used tubulin-binding agents have necessitated the pursuit of novel drug candidates with increased therapeutic potency. The design of novel drug candidates can be achieved using efficient computational techniques to support existing efforts. Proteochemometric (PCM) modeling is a computational technique that can be employed to elucidate the bioactivity relations between related targets and multiple ligands. We have developed a PCM-based Support Vector Machine (SVM) approach for predicting the bioactivity between tubulin receptors and small, drug-like molecules. The bioactivity datasets used for training the SVM algorithm were obtained from the Binding DB database. The SVM-based PCM model yielded a good overall predictive performance with an area under the curve (AUC) of 87%, Matthews correlation coefficient (MCC) of 72%, overall accuracy of 93%, and a classification error of 7%. The algorithm allows the prediction of the likelihood of new interactions based on confidence scores between the query datasets, comprising ligands in SMILES format and protein sequences of tubulin targets. The algorithm has been implemented as a web server known as TubPred, accessible via http://35.167.90.225:5000/ .

From Magic Bullet to Magic Bomb: Reductive Bioactivation of Antiparasitic Agents

Paul Ehrlich coined the term "magic bullet" to describe how a drug kills the parasite inside its human host without harming the host itself. Ehrlich concluded that the drug must have a greater affinity to the parasite than to human cells. Today, the specificity of drug action is understood in terms of the drug target. An ideal target is a protein that is essential for the proliferation of the pathogen but absent in human cells. Examples are the enzymes of folate synthesis or of the nonmevalonate pathway in the malaria parasites. However, there are other ways how a drug can kill selectively. Of particular relevance is the specific activation of a prodrug inside the pathogen but not in the host, as this is how the current frontrunners of parasite chemotherapy work. Artemisinins for malaria, fexinidazole for human African trypanosomiasis, benznidazole for Chagas' disease, metronidazole for intestinal protozoa: these molecules are "magic bombs" that are triggered selectively. They are prodrugs that need to be activated by chemical reduction, i.e., the acquisition of an electron, which occurs in the parasite. Such a mode of action is shared by the novel antimalarial peroxides arterolane and artefenomel, which are activated by reduction of the endoperoxide bond with ferrous heme as the likely electron donor, a metabolic end-product of Plasmodium falciparum. Here we provide an overview on the molecular basis of selectivity of antiparasitic drug action with particular reference to the ozonides, the new generation of antimalarial peroxides designed by Jonathan Vennerstrom.

Susceptible trichostrongyloid species mask presence of benzimidazole-resistant Haemonchus contortus in cattle

Background

Benzimidazole (BZ) anthelmintics are widely used to control infections with parasitic nematodes, but BZ resistance is an emerging threat among several nematode species infecting humans and animals. In Sudan, BZ-resistant Haemonchus contortus populations were recently reported in goats in South Darfur State. The objective of this study was to collect data regarding the situation of BZ resistance in cattle parasitic nematodes in South Darfur using phenotypic and molecular approaches, besides providing some epidemiological data on nematodes in cattle.

Methods

The faecal egg count reduction test and the egg hatch test (EHT) were used to evaluate benzimidazole efficacy in cattle nematodes in five South Darfur study areas: Beleil, Kass, Nyala, Rehed Al-Birdi and Tulus. Genomic DNA was extracted from pools of third-stage larvae (L3) (n = 40) during trials, before and after treatment, and pools of adult male Haemonchus spp. (n = 18) from abattoirs. The polymorphisms F167Y, E198A and F200Y in isotype 1 β-tubulin genes of H. contortus and H. placei were analysed using Sanger and pyrosequencing.

Results

Prevalence of gastro-intestinal helminths in cattle was 71% (313/443). Reduced albendazole faecal egg count reduction efficacy was detected in three study areas: Nyala (93.7%), Rehed Al-Birdi (89.7%) and Tulus (88.2%). In the EHT, EC₅₀ values of these study areas ranged between 0.032 and 0.037 µg/ml thiabendazole. Genus-specific PCRs detected the genera Haemonchus, Trichostrongylus and Cooperia in L3 samples collected after albendazole treatment. Sanger sequencing followed by pyrosequencing assays did not detect elevated frequencies of known BZ resistance-associated alleles in codon F167Y, E198A and F200Y in isotype 1 β-tubulin gene of H. placei (≤ 11.38%). However, polymorphisms were detected in H. contortus and in samples with mixed infections with H. contortus and H. placei at codon 198, including E198L (16/58), E198V (2/58) and potentially E198Stop (1/58). All pooled L3 samples post-albendazole treatment (n = 13) were identified as H. contortus with an E198L substitution at codon 198.

Conclusions

To the knowledge of the authors, this is the first report of reduced albendazole efficacy in cattle in Sudan and is the first study describing an E198L substitution in phenotypically BZ-resistant nematodes collected from cattle.

Haemonchosis: A Challenging Parasitic Infection of Sheep and Goats

The paper reviews the challenges about haemonchosis-a significant and common parasitic infection of small ruminants. Haemonchus contortus is a highly pathogenic parasite that localises in the abomasum of affected animals and exerts its pathogenicity by blood-sucking activity, adversely affecting the health and productivity of animals. The first challenge is the uneven distribution of the infection globally, this being more prevalent in tropical and subtropical and warm temperate and summer rainfall regions than in cool and cold temperate and arid regions; hence, this leads in differences in the approaches required for its control. Another challenge is the widespread presence of Haemonchus strains resistant to the various anthelmintics available: Benzimidazoles, imidazothiazoles, macrocyclic lactones, closantel and monepantel, which makes the control of the infection difficult. The third challenge refers to the difficulty of diagnosing the disease, given that field evidence can provide suspicion about the infection, which needs to be subsequently confirmed by laboratory tests through parasitological or molecular techniques. The final challenge relates to the difficulties in the control of the infection and the necessity to use pharmaceutical products cautiously and with a planned approach, to avoid further development of anthelmintic resistance, also given that use of a recently licenced vaccine is not widespread. In conclusion, at the moment, we should be concerned, but not worried, about this infection, and apply correctly the appropriate health management plans.

A 4 year observation of gastrointestinal nematode egg counts, nemabiomes and the benzimidazole resistance genotypes of Teladorsagia circumcincta on a Scottish sheep farm

Anthelmintic resistance threatens the sustainability of sheep production globally. Advice regarding strategies to reduce the development of anthelmintic resistance incorporates the outcomes of modelling exercises. Further understanding of gastrointestinal nematode species diversity, and population dynamics and genetics (which may vary between species) is required to refine these models; and field studies combining faecal egg outputs, species composition and resistance genetics are needed to calibrate them. In this study, faecal samples were taken from ewes and lambs on a commercial farm in south-eastern Scotland at approximately 3 t-4 week intervals between spring and autumn over a period of 4 years. Faecal egg counts were performed on these samples, and L₃ were collected from pooled coprocultures. Deep amplicon sequencing was used to determine both the species composition of these L₃ and the proportions of benzimidazole-resistant single nucleotide polymorphisms in the isotype-1 β-tubulin locus of the predominant species, Teladorsagia circumcincta L₃. Despite consistent management throughout the study, the results show variation in gastrointestinal nematode species composition with time and between age groups, that was potentially associated with weather conditions. The F200Y benzimidazole resistance mutation is close to genetic fixation in the T. circumcincta population on this farm. There was no evidence of variation in isotype-1 β-tubulin single nucleotide polymorphisms frequency between age groups, and no genetic evidence of reversion to benzimidazole susceptibility, despite targeted benzimidazole usage. This study highlights the need to include speciation when investigating gastrointestinal nematode epidemiology and anthelmintic resistance, and serves as an example of how genetic data may be analysed alongside species diversity and faecal egg counts, when markers for other anthelmintic classes are identified.

Considerations for anthelmintic resistance emergence in hookworm at a single locus

Over 800 million people are infected with hookworms around the world. Hookworms of the genus Ancylostoma and Necator are examples of nematodes that harbor the ability to enter a host by penetrating the skin, and after entry the infective larvae migrate to the small intestine where they encounter host-specific signals that initiate developmental pathways and culminate in maturation to the adult stage. Currently no vaccine is available for the treatment of hookworm infection. The control strategy is limited to anthelmintic drugs, which run the risk of losing efficacy as resistance grows. Genetic resistance has developed against all classes of anthelmintic drugs against livestock parasites, and recently markers of anthelmintic resistance in human hookworm populations have been reported. As anthelmintic resistance develops in human populations of hookworm, new drugs and novel control methods like vaccines will be required in the future to control hookworm transmission. This review outlines how population genetics and anthelmintic resistance could interact at a single locus to influence current control strategies.

•

This review outlines how population genetics and anthelmintic resistance may interact to influence current control strategies.

•

Altering mutation rate modifies rate of emergence of anthelmintic resistance.

•

Genotypic fitness changes the efficacy of anthelmintic administration.

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.

Molecular and phenotypic characterisation of fenbendazole resistance in a field-derived isolate of Ostertagia ostertagi

The prevalence of anthelmintic resistance in the bovine nematode Cooperia oncophora has been well documented globally but lack of efficacy against the more pathogenic nematode species Ostertagia ostertagi is less common. The sensitivity of an O. ostertagi isolate to the benzimidazole class of anthelmintic was investigated using classical parasitological techniques following apparent clinical failure of controlled release fenbendazole capsule administration in first season grazers at pasture. A controlled efficacy test (CET) was conducted in conjunction with sequencing of the β-tubulin isotype 1 gene of larvae pre- and post-fenbendazole administration. Twelve helminth-naïve calves were infected experimentally with 20,000 third stage larvae; six received oral fenbendazole (7.5 mg/kg bodyweight) 28 days post infection. Total abomasal nematode burdens were compared between treatment and control groups to determine efficacy. Fenbendazole resistance in O. ostertagi was confirmed with a total treatment failure in reducing worm burden: efficacy of 0%. Sequence analysis of the β-tubulin isotype-1 gene from forty-five infective larvae from both control and treated groups was performed. The three commonest single nucleotide polymorphisms (SNPs) associated with benzimidazole resistance, namely F167Y, E198A and F200Y, were examined. The predominant resistance-associated SNPs were F200Y (78 % control and 79 % treated groups) and F167Y (remaining genotypes) and emphasises the importance of these SNPs in clinical disease in this isolate. The development of diagnostic molecular tools based on a characterised field-derived isolate of benzimidazole-resistant Ostertagia will enable future prevalence surveys to be undertaken to assess the possible risk posed by resistance in this economically important species.

Haemonchus sp. in beef cattle in Brazil: species composition and frequency of benzimidazole resistance alleles

The aim of the present study was to investigate the occurrence of Haemonchus contortus and Haemonchus placei in beef cattle and the frequency of single nucleotide polymorphisms associated with benzimidazole (BZ)-resistance in Haemonchus spp. in Brazil. For such, fecal samples were collected from 61 beef cattle ranches in 11 Brazilian states. Third-stage larvae (L3) were produced for morphological species identification and DNA extraction. PCR was performed for the analysis of the isotype 1 β-tubulin gene and the products were sequenced to confirm the presence of H. placei and H. contortus. For each field population, pyrosequencing assays were performed to quantify the frequency of the F167Y, E198A and F200Y polymorphisms in the isotype-1 β-tubulin gene. The results of the morphometric analysis of 2345 larvae showed that H. placei was present on all ranches. The analysis of the isotype 1 β-tubulin gene confirmed 100% prevalence for H. placei and 23.7% for H. contortus. Pyrosequencing assays demonstrated single-nucleotide polymorphisms (SNPs) associated with BZ-resistance in all three codons (F167Y, E198A and F200Y) of the isotype 1 β-tubulin gene in H. placei field populations. Frequencies of resistance-associated alleles above background (≥ 15%) were found for at least one codon in 11.4% of the field isolates and maximum frequencies of 30, 21 and 29% were found for codons 167, 198 and 200, respectively, on individual ranches. This study confirms the presence of H. contortus in beef cattle in the major livestock farming states in Brazil and demonstrates that genotypes associated with BZ resistance are present in field populations of Haemonchus spp..

Caenorhabditis elegans in anthelmintic research - Old model, new perspectives

For more than four decades, the free-living nematode Caenorhabditis elegans has been extensively used in anthelmintic research. Classic genetic screens and heterologous expression in the C. elegans model enormously contributed to the identification and characterization of molecular targets of all major anthelmintic drug classes. Although these findings provided substantial insights into common anthelmintic mechanisms, a breakthrough in the treatment and control of parasitic nematodes is still not in sight. Instead, we are facing increasing evidence that the enormous diversity within the phylum Nematoda cannot be recapitulated by any single free-living or parasitic species and the development of novel broad-spectrum anthelmintics is not be a simple goal. In the present review, we summarize certain milestones and challenges of the C. elegans model with focus on drug target identification, anthelmintic drug discovery and identification of resistance mechanisms. Furthermore, we present new perspectives and strategies on how current progress in C. elegans research will support future anthelmintic research.

Quantitative benzimidazole resistance and fitness effects of parasitic nematode beta-tubulin alleles

Infections by parasitic nematodes inflict a huge burden on the health of humans and livestock throughout the world. Anthelmintic drugs are the first line of defense against these infections. Unfortunately, resistance to these drugs is rampant and continues to spread. To improve treatment strategies, we must understand the genetics and molecular mechanisms that underlie resistance. Studies of the fungus Aspergillus nidulans and the free-living nematode Caenorhabditis elegans discovered that a beta-tubulin gene is mutated in benzimidazole (BZ) resistant strains. In parasitic nematode populations, three beta-tubulin alleles, F167Y, E198A, and F200Y, have long been correlated with resistance. Additionally, improvements in sequencing technologies have identified new alleles - E198V, E198L, E198K, E198I, and E198Stop - also correlated with BZ resistance. However, none of these alleles have been proven to cause resistance. To empirically demonstrate this point, we independently introduced the F167Y, E198A, and F200Y alleles as well as two of the newly identified alleles, E198V and E198L, into the BZ susceptible C. elegans N2 genetic background using the CRISPR-Cas9 system. These genome-edited strains were exposed to both albendazole and fenbendazole to quantitatively measure animal responses to BZs. We used a range of concentrations for each BZ compound to define response curves and found that all five of the alleles conferred resistance to BZ compounds equal to a loss of the entire beta-tubulin gene. These results prove that the parasite beta-tubulin alleles cause resistance. The E198V allele is found at low frequencies along with the E198L allele in natural parasite populations, suggesting that it could affect fitness. We performed competitive fitness assays and demonstrated that the E198V allele reduces animal health, supporting the hypothesis that this allele might be less fit in field populations. Overall, we present a powerful platform to quantitatively assess anthelmintic resistance and effects of specific resistance alleles on organismal fitness in the presence or absence of the drug.

•

Beta-tubulin alleles from parasitic nematodes were introduced into drug-susceptible C. elegans using CRISPR-Cas9.

•

Beta-tubulin variants are strongly selected in albendazole competitive fitness assays.

•

The E198V allele confers a fitness cost in control conditions as compared to other alleles.

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.

High frequency of benzimidazole resistance alleles in trichostrongyloids from Austrian sheep flocks in an alpine transhumance management system

Background

Infections of small ruminants with trichostrongyloid nematodes often result in reduced productivity and may be detrimental to the host. Anthelmintic resistance (AR) against most anthelmintic drug classes is now widespread amongst the trichostrongyloids. Baseline establishment, followed by regular monitoring of the level of AR, is necessary for farmers and veterinarians to make informed decisions about parasite management. The detection of single nucleotide polymorphisms (SNPs) is a sensitive method to detect AR against benzimidazoles (BZs), one of the most widely used anthelmintic classes. Alpine transhumance constitutes a special type of pasturing of sheep from many different farms, the aim of this study was to investigate the prevalence of benzimidazole resistance alleles in this particular management system.

Results

Sixteen sheep flocks in Styria and Salzburg in Austria were examined by pyrosequencing for SNPs at codons 167, 198 and 200 of the isotype-1 β-tubulin gene. The frequency of the resistance-associated exchange F200Y was 87–100% for H. contortus, 77–100% for T. colubriformis and <  5–66% for T. circumcincta. Additionally, the F167Y polymorphism was detected in T. colubriformis from two farms at a frequency of 19 and 23% respectively.

Conclusions

The high resistance allele frequency in H. contortus and T. colubriformis in the examined sheep population urgently calls for the development of new treatment strategies to sustainably control trichostrongyloid infections for this kind of pasturing, since the frequent mixing of flocks during the alpine summer grazing must be considered an important risk factor for the spread of resistant nematodes to a large number of farms.

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.

Assessing anthelmintic resistance risk in the post-genomic era: a proof-of-concept study assessing the potential for widespread benzimidazole-resistant gastrointestinal nematodes in North American cattle and bison

As genomic research continues to improve our understanding of the genetics of anthelmintic drug resistance, the revolution in DNA sequencing technologies will provide increasing opportunities for large-scale surveillance for the emergence of drug resistance. In most countries, parasite control in cattle and bison has mainly depended on pour-on macrocyclic lactone formulations resulting in widespread ivermectin resistance. Consequently, there is an increased interest in using benzimidazole drugs which have been used comparatively little in cattle and bison in recent years. This situation, together with our understanding of benzimidazole resistance genetics, provides a practical opportunity to use deep-amplicon sequencing to assess the risk of drug resistance emergence. In this paper, we use deep-amplicon sequencing to scan for those mutations in the isotype-1 β-tubulin gene previously associated with benzimidazole resistance in many trichostrongylid nematode species. We found that several of these mutations occur at low frequency in many cattle and bison parasite populations in North America, suggesting increased use of benzimidazole drugs in cattle has the potential to result in widespread emergence of resistance in multiple parasite species. This work illustrates a post-genomic approach to large-scale surveillance of early emergence of anthelmintic resistance in the field.

Large scale screening for benzimidazole resistance mutations in Nematodirus battus, using both pyrosequence genotyping and deep amplicon sequencing, indicates the early emergence of resistance on UK sheep farms

Benzimidazoles (BZ) have been the anthelmintic of choice for controlling Nematodirus battus infections since their release in the 1950s. Despite heavy reliance on this single anthelmintic drug class, resistance was not identified in this nematode until 2010 (Mitchell et al., 2011). The study aimed to explore the prevalence of BZ-resistance mutations in N. battus from UK sheep flocks using deep amplicon sequencing and pyrosequencing platforms. Based on evidence from other gastrointestinal nematodes, resistance in N. battus is likely to be conferred by single nucleotide polymorphisms (SNP) within the β-tubulin isotype 1 locus at codons 167, 198 and 200. Pyrosequencing and deep amplicon sequencing assays were designed to identify the F167Y (TTC to TAC), E198A (GAA to GCA) and F200Y (TTC to TAC) SNPs. Nematodirus battus populations from 253 independent farms were analysed by pyrosequencing; 174 farm populations were included in deep amplicon sequencing and 170 were analysed using both technologies. F200Y was the most prevalent SNP identified throughout the UK, in 12-27% of the populations tested depending on assay, at a low overall individual frequency of 2.2 ± 0.6% (mean ± SEM, based on pyrosequencing results). Four out of the five populations with high frequencies (>20%) of the F200Y mutation were located in NW England. The F167Y SNP was identified, for the first time in this species, in four of the populations tested at a low frequency (1.2% ± 0.01), indicating the early emergence of the mutation. E198A or E198L were not identified in any of the isolates. Results obtained were comparable between both techniques for F200Y (Lins' CCC, r_c = 0.96) with discrepancies being limited to populations with low frequencies. The recent emergence of resistance in this species will provide a unique opportunity to study the early stages of anthelmintic resistance within a natural setting and track its progress in the future.

A perspective on the discovery of selected compounds with anthelmintic activity against the barber's pole worm-Where to from here?

Parasitic roundworms (nematodes) cause substantial morbidity and mortality in animals worldwide. Anthelmintic treatment is central to controlling these worms, but widespread resistance to most of the commercially available anthelmintics for veterinary and agricultural use is compromising control, such that there is an urgency to discover new and effective drugs. The purpose of this article is to review information on parasitic nematodes, the treatment and control of parasitic nematode infections and aspects of discovering new anthelmintics in the context of anthelmintic resistance problems, and then to discuss some progress that our group has made in identifying selected compounds with activity against nematodes. The focus of our recent work has been on discovering new chemical entities and known drugs with anthelmintic activities against Haemonchus contortus as well as other socioeconomically important parasitic nematodes for subsequent development. Using whole worm-based phenotypic assays, we have been screening compound collections obtained via product-development-partnerships and/or collaborators, and active compounds have been assessed for their potential as anthelmintic candidates. Following the screening of 15,333 chemicals from five distinct compound collections against H. contortus, we have discovered one new chemical entity (designated SN00797439), two human kinase inhibitors (SNS-032 and AG-1295), 14 tetrahydroquinoxaline analogues, one insecticide (tolfenpyrad) and two tolfenpyrad (pyrazole-5-carboxamide) derivatives (a-15 and a-17) with anthelmintic activity in vitro. Some of these 20 'hit' compounds have selectivity against H. contortus in vitro when compared to particular human cell lines. In our opinion, some of these compounds could represent starting points for 'lead' development. Accordingly, the next research steps to be pursued include: (i) chemical optimisation of representative chemicals via structure-activity relationship (SAR) evaluations; (ii) assessment of the breadth of spectrum of anthelmintic activity on a range of other parasitic nematodes, such as strongyloids, ascaridoids, enoplids and filarioids; (iii) detailed investigations of the absorption, distribution, metabolism, excretion and toxicity (ADMET) of optimised chemicals with broad nematocidal or nematostatic activity; and (iv) establishment of the modes of action of lead candidates.

New codon 198 β-tubulin polymorphisms in highly benzimidazole resistant Haemonchus contortus from goats in three different states in Sudan

Background

Benzimidazole (BZ) resistance in gastrointestinal nematodes is a worldwide problem for livestock production, particularly in small ruminants. Assignment of the emergence of resistance using sensitive and reliable methods is required to adopt the correct strategies for control. In Sudan, BZ resistant Haemonchus contortus populations were recently reported in goats in South Darfur. This study aimed to provide additional data regarding albendazole efficacy and to describe the prevailing molecular BZ resistance mechanisms.

Methods

Faecal egg count reduction and egg hatch tests (EHT) were used to evaluate albendazole efficacy in three different areas of South Darfur using naturally (Rehed Al-Birdi and Tulus) and experimentally infected (Tulus and Um Dafuq) goats. Using samples from Central, East and South Darfur, pyro- and Sanger sequencing were used to detect the polymorphisms F167Y, E198A and F200Y in H. contortus isotype 1 β-tubulin in DNA extracted from pooled third-stage larval (L3) samples (n = 36) on days 0 and 10 during trials, and from pooled adult male H. contortus (treated goats, n = 14; abattoirs, n = 83) including samples from populations previously found to be resistant in South Darfur.

Results

Albendazole efficacies at 5, 7.5 and 10 mg/kg doses were 73.5–90.2% on day 14 in natural and experimental infections while 12.5 mg/kg showed > 96.6% efficacy. EC₅₀ in the EHT were 0.8 and 0.11 µg/ml thiabendazole in natural and experimental infection trials, respectively. PCRs detected Haemonchus, Trichostrongylus and Cooperia in L3 samples from albendazole-treated goats. Haemonchus contortus allele frequencies in codons 167 and 200 using pyrosequencing assays were ≤ 7.4% while codon 198 assays failed. Sanger sequencing revealed five novel polymorphisms at codon 198. Noteworthy, an E198L substitution was present in 82% of the samples (L3 and adults) including all post-treatment samples. Moreover, E198V, E198K and potentially E198I, and E198Stop were identified in a few samples.

Conclusions

To our knowledge, this is the first report of E198L in BZ resistant H. contortus and the second where this is the predominant genotype associated with resistance in any strongyle species. Since this variant cannot be quantified using pyrosequencing, the results highlight important limitations in the general applicability of pyrosequencing to quantify BZ resistance genotypes.

The critical importance of planned small ruminant livestock health and production in addressing global challenges surrounding food production and poverty alleviation

Small ruminants are particularly well suited to meet United Nations Sustainable Development Goals surrounding food security, human wellbeing and poverty alleviation in different environmental and climatic settings. However the current efficiency of food production from small ruminants in both developed agricultural regions and in lower and middle income countries is woefully inadequate to meet predicted global needs over the forthcoming decades. Most global research to address this challenge is focussed on the genetics of animal growth, conformation and disease tolerance or resistance traits, albeit the practical consequences of such selection and strategies to use genetically improved animals in the field are uncertain. Any long-term benefits derived from small ruminant genetic selection will only be impactful if steps are first taken to keep animals alive, healthy and productive through iterative planned health management. Parasites are the foremost global infectious disease constraints to efficient small ruminant production. Their genetic adaptability to exploit opportunities afforded by effects of climatic or management changes on free-living stages, or exposure of parasitic stages to drugs, presents specific challenges to their sustainable control. Hence, parasite control provides a relevant means of engagement with livestock keepers and farmers on the topic of planned animal health management. The value of parasitology in this regard is enhanced by the availability of simple to use and accessible diagnostic tools.

Drug resistance in Giardia: Mechanisms and alternative treatments for Giardiasis

The use of chemotherapeutic drugs is the main resource against clinical giardiasis due to the lack of approved vaccines. Resistance of G. duodenalis to the most used drugs to treat giardiasis, metronidazole and albendazole, is a clinical issue of growing concern and yet unknown impact, respectively. In the search of new drugs, the completion of the Giardia genome project and the use of biochemical, molecular and bioinformatics tools allowed the identification of ligands/inhibitors for about one tenth of ≈150 potential drug targets in this parasite. Further, the synthesis of second generation nitroimidazoles and benzimidazoles along with high-throughput technologies have allowed not only to define overall mechanisms of resistance to metronidazole but to screen libraries of repurposed drugs and new pharmacophores, thereby increasing the known arsenal of anti-giardial compounds to some hundreds, with most demonstrating activity against metronidazole or albendazole-resistant Giardia. In particular, cysteine-modifying agents which include omeprazole, disulfiram, allicin and auranofin outstand due to their pleiotropic activity based on the extensive repertoire of thiol-containing proteins and the microaerophilic metabolism of this parasite. Other promising agents derived from higher organisms including phytochemicals, lactoferrin and propolis as well as probiotic bacteria/fungi have also demonstrated significant potential for therapeutic and prophylactic purposes in giardiasis. In this context the present chapter offers a comprehensive review of the current knowledge, including commonly prescribed drugs, causes of therapeutic failures, drug resistance mechanisms, strategies for the discovery of new agents and alternative drug therapies.

Similarities and differences in the biotransformation and transcriptomic responses of Caenorhabditis elegans and Haemonchus contortus to five different benzimidazole drugs

We have undertaken a detailed analysis of the biotransformation of five of the most therapeutically important benzimidazole anthelmintics - albendazole (ABZ), mebendazole (MBZ), thiabendazole (TBZ), oxfendazole (OxBZ) and fenbendazole (FBZ) - in Caenorhabditis elegans and the ruminant parasite Haemonchus contortus. Drug metabolites were detected by LC-MS/MS analysis in supernatants of C. elegans cultures with a hexose conjugate, most likely glucose, dominating for all five drugs. This work adds to a growing body of evidence that glucose conjugation is a major pathway of xenobiotic metabolism in nematodes and may be a target for enhancement of anthelmintic potency. Consistent with this, we found that biotransformation of albendazole by C. elegans reduced drug potency. Glucose metabolite production by C. elegans was reduced in the presence of the pharmacological inhibitor chrysin suggesting that UDP-glucuronosyl/glucosyl transferase (UGT) enzymes may catalyze benzimidazole glucosidation. Similar glucoside metabolites were detected following ex vivo culture of adult Haemonchus contortus. As a step towards identifying nematode enzymes potentially responsible for benzimidazole biotransformation, we characterised the transcriptomic response to each of the benzimidazole drugs using the C. elegans resistant strain CB3474 ben-1(e1880)III. In the case of albendazole, mebendazole, thiabendazole, and oxfendazole the shared transcriptomic response was dominated by the up-regulation of classical xenobiotic response genes including a shared group of UGT enzymes (ugt-14/25/33/34/37/41/8/9). In the case of fenbendazole, a much greater number of genes were up-regulated, as well as developmental and brood size effects suggesting the presence of secondary drug targets in addition to BEN-1. The transcriptional xenobiotic response of a multiply resistant H. contortus strain UGA/2004 was essentially undetectable in the adult stage but present in the L3 infective stage, albeit more muted than C. elegans. This suggests that xenobiotic responses may be less efficient in stages of parasitic nematodes that reside in the host compared with the free-living stages.

•

C. e. & H. c. display hexose conjugation (likely glucose) and excretion of 5 BZs.

•

C. elegans (C.e.) biotransformation of ABZ reduces drug potency.

•

UGT inhibitor chrysin reduces ABZ biotransformation by C. elegans.

•

Transcriptomic response of C. e. (ben-1) to 5 BZs dominated by xenobiotic response and additional targets for FBZ.

•

Minimal transcriptomic response of H. contortus to ABZ exposure.

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.

A detection of benzimidazole resistance-associated SNPs in the isotype-1 β-tubulin gene in Haemonchus contortus from wild blue sheep (Pseudois nayaur) sympatric with sheep in Helan Mountains, China

Background

Benzimidazole (BZ) resistance is an increasingly serious problem due to the excessive use of this anthelmintic for controlling Haemonchus contortus, which is one of the major gastrointestinal nematodes infecting small ruminants worldwide. Three known single nucleotide polymorphisms (SNPs), F167Y (TAC), E198A (GCA) and F200Y (TAC), in the isotype-1 β-tubulin gene of H. contortus are associated with BZ resistance. Comprehending the spread and origins of BZ resistance-associated SNPs has important implications for the control of this nematode.

Results

Twenty-seven adult H. contortus were harvested from wild blue sheep (Pseudois nayaur), small wild ruminants sympatric with domestic ruminants, inhabiting the Helan Mountains, China, to monitor the status of BZ resistance. In addition, 20 adult H. contortus from domestic sheep sympatric with this wild ruminant and 36 isotype-1 β-tubulin haplotype sequences of H. contortus (two of these haplotypes, E198A3 and E198A4, possessed resistance-associated SNP E198A (GCA) from domestic ruminants in eight other geographical regions of China were used to further define the origins of BZ resistance-associated SNPs within the worms collected from blue sheep. The BZ resistance-associated SNP E198A was detected, whereas SNPs F167Y (TAC) and F200Y (TAC) were not found within the worms collected from blue sheep, and the frequency of homozygous resistant E198A (GCA) was 7.40%. The evolutionary tree and network showed consistent topologies for which there was no obvious boundary among the worms from the wild and domestic hosts, and two haplotypes (E198A1 and E198A2) possessing E198A from the wild blue sheep had two different independent origins. E198A1 had the same origin with E198A3 but E198A2 had a different origin with them. Population genetic analyses revealed a low level of Fst values (ranging from 0 to 0.19749) between all H. contortus worm groups in China.

Conclusions

Results of the current study of the three BZ resistance-associated SNPs of H. contortus from wild blue sheep suggested that only E198A (GCA) was present within the worms collected from the wild ruminants and had multiple independent origins.

Electronic supplementary material

The online version of this article (10.1186/s12917-019-1838-4) contains supplementary material, which is available to authorized users.

Isolation and characterization of a naturally occurring multidrug-resistant strain of the canine hookworm, Ancylostoma caninum

Soil-transmitted nematodes infect over a billion people and place several billion more at risk of infection. Hookworm disease is the most significant of these soil-transmitted nematodes, with over 500 million people infected. Hookworm infection can result in debilitating and sometimes fatal iron-deficiency anemia, which is particularly devastating in children and pregnant women. Currently, hookworms and other soil-transmitted nematodes are controlled by administration of a single dose of a benzimidazole to targeted populations in endemic areas. While effective, people are quickly re-infected, necessitating frequent treatment. Widespread exposure to anthelmintic drugs can place significant selective pressure on parasitic nematodes to generate resistance, which has severely compromised benzimidazole anthelmintics for control of livestock nematodes in many areas of the world. Here we report, to our knowledge, the first naturally occurring multidrug-resistant strain of the canine hookworm Ancylostoma caninum. We reveal that this isolate is resistant to fenbendazole at the clinical dosage of 50 mg/kg for 3 days. Our data shows that this strain harbors a fixed, single base pair mutation at amino acid 167 of the β-tubulin isotype 1 gene, and by using CRISPR/Cas9 we demonstrate that introduction of this mutation into the corresponding amino acid in the orthologous β-tubulin gene of Caenorhabditis elegans confers a similar level of resistance to thiabendazole. We also show that the isolate is resistant to the macrocyclic lactone anthelmintic ivermectin. Understanding the mechanism of anthelmintic resistance is important for rational design of control strategies to maintain the usefulness of current drugs, and to monitor the emergence of resistance. The isolate we describe represents the first multidrug-resistant strain of A. caninum reported, and our data reveal a resistance marker that can emerge naturally in response to heavy anthelminthic treatment.