Analysis of putative resistance gene loci in UK field populations of Haemonchus contortus after 6years of macrocyclic lactone use

The equine ascarids: resuscitating historic model organisms for modern purposes

The equine ascarids, Parascaris spp., are important nematode parasites of juvenile horses and were historically model organisms in the field of cell biology, leading to many important discoveries, and are used for the study of chromatin diminution. In veterinary parasitology, Parascaris spp. are important not only because they can cause clinical disease in young horses but also because they are the only ascarid parasites to have developed widespread anthelmintic resistance. Despite this, much of the general biology and mechanisms of anthelmintic resistance are poorly understood. This review condenses known basic biological information and knowledge on the mechanisms of anthelmintic resistance in Parascaris spp., highlighting the importance of foundational research programs. Although two variants of this parasite were recognized based on the number of chromosomes in the 1870s and suggested to be two species in 1890, one of these, P. univalens, appears to have been largely forgotten in the veterinary scientific literature over the past 100 years. We describe how this omission has had a century-long effect on nomenclature and data analysis in the field, highlighting the importance of proper specimen identification in public repositories. A summary of important basic biology, including life cycle, in vitro maintenance, and immunology, is given, and areas of future research for the improvement of knowledge and development of new systems are given. Finally, the limited knowledge regarding anthelmintic resistance in Parascaris spp. is summarized, along with caution regarding assumptions that resistance mechanisms can be applied across clades.

Comparative analysis on transcriptomics of ivermectin resistant and susceptible strains of Haemonchus contortus

Background

Ivermectin (IVM) is one of the most important and widely used anthelmintics in veterinary medicine. However, its efficacy is increasingly compromised by widespread resistance, and the exact mechanism of IVM resistance remains unclear for most parasitic nematodes, including Haemonchus contortus, a blood-sucking parasitic nematode of small ruminants.

Methods

In this study, an H. contortus IVM-resistant strain from Zhaosu, Xinjiang, China, was isolated and assessed by the control test, faecal egg count reduction test (FECRT) and the larval development assay (LDA). Subsequently, comparative analyses on the transcriptomics of IVM-susceptible and IVM-resistant adult worms of this parasite were carried out using RNA sequencing (RNA-seq) and bioinformatics.

Results

In total, 543 (416 known, 127 novel) and 359 (309 known, 50 novel) differentially expressed genes (DEGs) were identified in male and female adult worms of the resistant strain compared with those of the susceptible strain, respectively. In addition to several previously known candidate genes which were supposed to be associated with IVM resistance and whose functions were involved in receptor activity, transport, and detoxification, we found some new potential target genes, including those related to lipid metabolism, structural constituent of cuticle, and important pathways such as antigen processing and presentation, lysosome, autophagy, apoptosis, and NOD1-like receptor signalling pathways. Finally, the results of quantitative real-time polymerase chain reaction confirmed that the transcriptional profiles of selected DEGs (male: 8 genes, female: 10 genes) were consistent with those obtained by the RNA-seq.

Conclusions

Our results indicate that IVM has multiple effects, including both neuromuscular and non-neuromuscular targets, and provide valuable information for further studies on the IVM resistance mechanism in H. contortus.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05274-y.

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

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

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

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

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

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

Transcriptional responses in Parascaris univalens after in vitro exposure to ivermectin, pyrantel citrate and thiabendazole

Background

Parascaris univalens is a pathogenic parasite of foals and yearlings worldwide. In recent years, Parascaris spp. worms have developed resistance to several of the commonly used anthelmintics, though currently the mechanisms behind this development are unknown. The aim of this study was to investigate the transcriptional responses in adult P. univalens worms after in vitro exposure to different concentrations of three anthelmintic drugs, focusing on drug targets and drug metabolising pathways.

Methods

Adult worms were collected from the intestines of two foals at slaughter. The foals were naturally infected and had never been treated with anthelmintics. Worms were incubated in cell culture media containing different concentrations of either ivermectin (10⁻⁹ M, 10⁻¹¹ M, 10⁻¹³ M), pyrantel citrate (10⁻⁶ M, 10⁻⁸ M, 10⁻¹⁰ M), thiabendazole (10⁻⁵ M, 10⁻⁷ M, 10⁻⁹ M) or without anthelmintics (control) at 37 °C for 24 h. After incubation, the viability of the worms was assessed and RNA extracted from the anterior region of 36 worms and sequenced on an Illumina NovaSeq 6000 system.

Results

All worms were alive at the end of the incubation but showed varying degrees of viability depending on the drug and concentration used. Differential expression (Padj < 0.05 and log2 fold change ≥ 1 or ≤ − 1) analysis showed similarities and differences in the transcriptional response after exposure to the different drug classes. Candidate genes upregulated or downregulated in drug exposed worms include members of the phase I metabolic pathway short-chain dehydrogenase/reductase superfamily (SDR), flavin containing monooxygenase superfamily (FMO) and cytochrome P450-family (CYP), as well as members of the membrane transporters major facilitator superfamily (MFS) and solute carrier superfamily (SLC). Generally, different targets of the anthelmintics used were found to be upregulated and downregulated in an unspecific pattern after drug exposure, apart from the GABA receptor subunit lgc-37, which was upregulated only in worms exposed to 10⁻⁹ M of ivermectin.

Conclusions

To our knowledge, this is the first time the expression of lgc-37 and members of the FMO, SDR, MFS and SLC superfamilies have been described in P. univalens and future work should be focused on characterising these candidate genes to further explore their potential involvement in drug metabolism and anthelmintic resistance.

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

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

Review of the Eprinomectin effective doses required for dairy goats: Where do we go from here?

Eprinomectin (EPM) has been recently granted a marketing authorisation in the European Union for use in goats, with a zero-day milk withdrawal period. Considering the high prevalence of benzimidazole resistance worldwide and the economic implications of managing milk residues, EPM may today be considered the main (or even the only) affordable treatment option, at least in dairy goats in the EU. However, the chosen dose (1 mg/kg) seems to be suboptimal, especially for lactating goats, and the chosen route of administration (Pour-on) highly subject to inter-individual variability. Considering the scarcity of anthelmintic resources, such a dosage regimen might threat the sustainability of this crucial drug in goat milk production and needs to be urgently discussed and reassessed.

Refugia and anthelmintic resistance: Concepts and challenges

Anthelmintic resistance is a threat to global food security. In order to alleviate the selection pressure for resistance and maintain drug efficacy, management strategies increasingly aim to preserve a proportion of the parasite population in 'refugia', unexposed to treatment. While persuasive in its logic, and widely advocated as best practice, evidence for the ability of refugia-based approaches to slow the development of drug resistance in parasitic helminths is currently limited. Moreover, the conditions needed for refugia to work, or how transferable those are between parasite-host systems, are not known. This review, born of an international workshop, seeks to deconstruct the concept of refugia and examine its assumptions and applicability in different situations. We conclude that factors potentially important to refugia, such as the fitness cost of drug resistance, the degree of mixing between parasite sub-populations selected through treatment or not, and the impact of parasite life-history, genetics and environment on the population dynamics of resistance, vary widely between systems. The success of attempts to generate refugia to limit anthelmintic drug resistance are therefore likely to be highly dependent on the system in hand. Additional research is needed on the concept of refugia and the underlying principles for its application across systems, as well as empirical studies within systems that prove and optimise its usefulness.

Population genomic and evolutionary modelling analyses reveal a single major QTL for ivermectin drug resistance in the pathogenic nematode, Haemonchus contortus

BACKGROUND

Infections with helminths cause an enormous disease burden in billions of animals and plants worldwide. Large scale use of anthelmintics has driven the evolution of resistance in a number of species that infect livestock and companion animals, and there are growing concerns regarding the reduced efficacy in some human-infective helminths. Understanding the mechanisms by which resistance evolves is the focus of increasing interest; robust genetic analysis of helminths is challenging, and although many candidate genes have been proposed, the genetic basis of resistance remains poorly resolved.

RESULTS

Here, we present a genome-wide analysis of two genetic crosses between ivermectin resistant and sensitive isolates of the parasitic nematode Haemonchus contortus, an economically important gastrointestinal parasite of small ruminants and a model for anthelmintic research. Whole genome sequencing of parental populations, and key stages throughout the crosses, identified extensive genomic diversity that differentiates populations, but after backcrossing and selection, a single genomic quantitative trait locus (QTL) localised on chromosome V was revealed to be associated with ivermectin resistance. This QTL was common between the two geographically and genetically divergent resistant populations and did not include any leading candidate genes, suggestive of a previously uncharacterised mechanism and/or driver of resistance. Despite limited resolution due to low recombination in this region, population genetic analyses and novel evolutionary models supported strong selection at this QTL, driven by at least partial dominance of the resistant allele, and that large resistance-associated haplotype blocks were enriched in response to selection.

CONCLUSIONS

We have described the genetic architecture and mode of ivermectin selection, revealing a major genomic locus associated with ivermectin resistance, the most conclusive evidence to date in any parasitic nematode. This study highlights a novel genome-wide approach to the analysis of a genetic cross in non-model organisms with extreme genetic diversity, and the importance of a high-quality reference genome in interpreting the signals of selection so identified.

GluClR-mediated inhibitory postsynaptic currents reveal targets for ivermectin and potential mechanisms of ivermectin resistance

Glutamate-gated chloride channel receptors (GluClRs) mediate inhibitory neurotransmission at invertebrate synapses and are primary targets of parasites that impact drastically on agriculture and human health. Ivermectin (IVM) is a broad-spectrum pesticide that binds and potentiates GluClR activity. Resistance to IVM is a major economic and health concern, but the molecular and synaptic mechanisms of resistance are ill-defined. Here we focus on GluClRs of the agricultural endoparasite, Haemonchus contortus. We demonstrate that IVM potentiates inhibitory input by inducing a tonic current that plateaus over 15 minutes and by enhancing post-synaptic current peak amplitude and decay times. We further demonstrate that IVM greatly enhances the active durations of single receptors. These effects are greatly attenuated when endogenous IVM-insensitive subunits are incorporated into GluClRs, suggesting a mechanism of IVM resistance that does not affect glutamate sensitivity. We discovered functional groups of IVM that contribute to tuning its potency at different isoforms and show that the dominant mode of access of IVM is via the cell membrane to the receptor.

Glutamate-gated chloride channel receptors (GluClRs) mediate chemoelectric inhibition in invertebrate animals and are targets for broad-spectrum pesticides such as ivermectin. However, resistance to ivermectin threatens the effective control of invertebrates that cause a range of agricultural and human diseases. This study investigates different isoforms of GluClR expressed by the major agricultural endoparasite, Haemonchus contortus, on a synaptic and single receptor level. We discovered that ivermectin enhances synaptic current amplitude and decay and lengthens single receptor activity. Furthermore, ivermectin is less efficacious at GluClRs that incorporate a naturally ivermectin-resistant subunit, suggesting a potential resistance mechanism. Finally, we identify two chemical interactions between the GluClR and ivermectin that determine its potency and show that ivermectin binds to GluClRs via cell membrane interactions.

Bioassay-guided isolation of active substances from Semen Torreyae identifies two new anthelmintic compounds with novel mechanism of action

ETHNOPHARMACOLOGICAL RELEVANCE

Semen Torreyae, the seeds of Torreya grandis Fortune ex Lindley (Cephalotaxaceae) is a well-known traditional Chinese medicinal plant recorded in the Chinese Pharmacopeia (2010 version). It is widely used for treating intestinal parasites in China, owing to its desirable efficacy and safety. However, the anthelmintic compounds in Semen Torreyae have not yet been identified.

AIM OF THE STUDY

This study aims to identify the compounds active against helminths from Semen Torreyae. In addition, we tested whether C. elegans strains resistant to currently-used anthelmintic drugs showed cross-resistance to these compounds.

METHODS

A bioassay-guided isolation of anthelmintic compounds from Semen Torreyae was performed using a Caenorhabditis elegans (C. elegans) testing model. The structures of active compounds were elucidated by a combination of GC-MS, high resolution MS, and NMR. The median-effect method was employed to generate a combination index (CI) to evaluate the synergistic effect of the anthelmintic compounds. A panel of C. elegans mutant strains resistant against the major anthelmintic drug classes was used to study the cross-resistance to currently-used anthelmintic drugs. A panel of transient receptor potential (TRP) channel mutant strains was also tested to explore the possible mechanisms of action of the anthelmintic compounds.

RESULTS

The bioassay-guided isolation led to two active compounds, i.e. galangal acetate (IC₅₀: 58.5 ± 8.9 μM) and miogadial (IC₅₀: 25.1 ± 5.4 μM). The combination of galangal acetate and miogadial resulted in a synergistic effect at IC₅₀, IC₇₀, and IC₉₀ levels (CIs < 1). Galangal acetate and miogadial demonstrated similar activity against drug-resistant C. elegans strains compared to the wild-type strain. In addition, none of the TRP mutants was significantly resistant to galangal acetate or miogadial compared to wild type worms.

CONCLUSIONS

We identified the bioactive compounds from Semen Torreyae responsible for its anthelmintic activity: galangal acetate and miogadial. The two anthelmintic compounds demonstrated a synergistic effect against C. elegans. Galangal acetate and miogadial are unlikely to act on the targets of currently-used anthelmintics (ivermectin, levamisole, benomyl and aldicarb), and an action on TRP channels appears to be ruled out as well. In summary, galangal acetate and miogadial are promising anthelmintic hits worth further investigation.

Genetic variants in dyf-7 validated by droplet digital PCR are not drivers for ivermectin resistance in Haemonchus contortus

Resistance to ivermectin (IVM) in the nematode Haemonchus contortus in small ruminants is an increasing problem throughout the world. Access to molecular diagnostics will allow early detection of IVM resistance, which in turn can limit the spread of resistant isolates. One candidate gene which has recently been suggested as a marker for IVM resistance is that for dye-filling protein (dyf-7). In this study, we critically investigated the suitability of A141G and G153T single nucleotide polymorphisms (SNP) of dyf-7 as a marker in larval cultures collected from sheep farms in Sweden, involving several isolates for which resistance status had been characterised by the faecal egg count reduction test (FECRT). Initially, we designed dyf-7 primers from a worldwide collection of adult Haemonchus contortus DNA. With the sequence data, we created a haplotype network. We then optimised and used the same sets of primers and probes in a droplet digital PCR (ddPCR) assay for precise quantification of dyf-7 allele frequencies in pre- and post-anthelmintic treatment faecal larval cultures. The fractional abundance (FA) of the mutant SNP was within the range 7.8 and 31%. However, the FA was generally stable in samples collected from the same farms, even though they were obtained on different occasions up to 25 months apart. There was also no indication that the level of IVM resistance as measured by the faecal egg count reduction test was higher on farms with high FA. Furthermore, by comparing FA in samples from the same farms pre- and post-IVM treatment, we found no evidence of a correlation between dyf-7 and level of IVM resistance. Based on these results, dyf-7 is not a suitable marker for field testing of IVM resistance in H. contortus.

Bioassay-guided isolation of three anthelmintic compounds from Warburgia ugandensis Sprague subspecies ugandensis, and the mechanism of action of polygodial

Parasitic helminths continue to pose problems in human and veterinary medicine, as well as in agriculture. Resistance to current anthelmintics has prompted the search for new drugs. Anthelmintic metabolites from medicinal plants could be good anthelmintic drug candidates. However, the compounds active against nematodes have not been identified in most medicinal plants with anthelmintic activity. In this study, we aimed to identify the active compounds against helminths in Warburgia ugandensis Sprague subspecies ugandensis (Canellaceae) and study the underlying mechanism of action. A bioassay-guided isolation of anthelmintic compounds from the plant was performed using a Caenorhabditis elegans (C. elegans) test model with a WMicrotracker instrument to monitor motility. Three active compounds were purified and identified by nuclear magnetic resonance and high resolution MS: warburganal (IC₅₀: 28.2 ± 8.6 μM), polygodial (IC₅₀: 13.1 ± 5.3 μM) and alpha-linolenic acid (ALA, IC₅₀: 70.1 ± 17.5 μM). A checkerboard assay for warburganal and ALA as well as polygodial and ALA showed a fractional inhibitory concentration index of 0.41 and 0.37, respectively, suggesting that polygodial and ALA, as well as warburganal and ALA, have a synergistic effect against nematodes. A preliminary structure-activity relationship study for polygodial showed that the α,β-unsaturated 1,4-dialdehyde structural motif is essential for the potent activity. None of a panel of C. elegans mutant strains, resistant against major anthelmintic drug classes, showed significant resistance to polygodial, implying that polygodial may block C. elegans motility through a mechanism which differs from that of currently marketed drugs. Further measurements showed that polygodial inhibits mitochondrial ATP synthesis of C. elegans in a dose-dependent manner (IC₅₀: 1.8 ± 1.0 μM). Therefore, we believe that the underlying mechanism of action of polygodial is probably inhibition of mitochondrial ATP synthesis. In conclusion, polygodial could be a promising anthelmintic drug candidate worth considering for further development.

Genomic introgression mapping of field-derived multiple-anthelmintic resistance in Teladorsagia circumcincta

Preventive chemotherapy has long been practiced against nematode parasites of livestock, leading to widespread drug resistance, and is increasingly being adopted for eradication of human parasitic nematodes even though it is similarly likely to lead to drug resistance. Given that the genetic architecture of resistance is poorly understood for any nematode, we have analyzed multidrug resistant Teladorsagia circumcincta, a major parasite of sheep, as a model for analysis of resistance selection. We introgressed a field-derived multiresistant genotype into a partially inbred susceptible genetic background (through repeated backcrossing and drug selection) and performed genome-wide scans in the backcross progeny and drug-selected F2 populations to identify the major genes responsible for the multidrug resistance. We identified variation linking candidate resistance genes to each drug class. Putative mechanisms included target site polymorphism, changes in likely regulatory regions and copy number variation in efflux transporters. This work elucidates the genetic architecture of multiple anthelmintic resistance in a parasitic nematode for the first time and establishes a framework for future studies of anthelmintic resistance in nematode parasites of humans.

Teladorsagia circumcincta is an economically significant nematode (roundworm) pathogen affecting sheep and goats in temperate regions of the world. The widespread use of prophylactic treatment has resulted in rapid selection for anthelmintic (anti-worm drug) resistance in this and other species of livestock parasites. The mechanism of resistance is not well understood because most studies have focused on the role of candidate genes using simplistic models of single gene selection, despite evidence that the evolution of resistance is more complex. Here, we report on a comprehensive whole-genome analysis that elucidated resistance-associated genes, which was facilitated by developing a pair of T. circumcincta strains sharing a largely common genetic background but differing markedly in their susceptibility to anthelmintic drugs. The results show that multiple genetic factors contribute to anthelmintic resistance in a variety of ways, including possible reduction/modulation in target site sensitivity, reduced target site expression, and increased drug efflux, to name a few. This suggests that drug resistance in these parasites is a multifactorial quantitative trait rather than a simple discrete Mendelian character. With this study, we established a genomics-based experimental paradigm for investigating anthelmintic resistance, at a time when its medical importance is rapidly increasing.

Evidence from two independent backcross experiments supports genetic linkage of microsatellite Hcms8a20, but not other candidate loci, to a major ivermectin resistance locus in Haemonchus contortus

Haemonchus contortus is the leading parasitic nematode species used to study anthelmintic drug resistance. A variety of candidate loci have been implicated as being associated with ivermectin resistance in this parasite but definitive evidence of their importance is still lacking. We have previously performed two independent serial backcross experiments to introgress ivermectin resistance loci from two H. contortus ivermectin-resistant strains - MHco4(WRS) and MHco10(CAVR) - into the genetic background of the ivermectin-susceptible genome reference strain MHco3(ISE). We have interrogated a number of candidate ivermectin resistance loci in the resulting backcross populations and assessed the evidence for their genetic linkage to an ivermectin resistance locus. These include the microsatellite marker Hcms8a20 and six candidate genes Hco-glc-5, Hco-avr-14, Hco-lgc-37 (previously designated Hco-hg-1), Hco-pgp-9 (previously designated Hco-pgp-1), Hco-pgp-2 and Hco-dyf-7. We have sampled the haplotype diversity of amplicon markers within, or adjacent to, each of these loci in the parental strains and fourth generation backcross populations to assess the evidence for haplotype introgression from the resistant parental strain into the genomic background of the susceptible parental strain in each backcross. The microsatellite Hcms8a20 locus showed strong evidence of such introgression in both independent backcrosses, suggesting it is linked to an important ivermectin resistance mutation in both the MHco4(WRS) and MHco10(CAVR) strains. In contrast, Hco-glc-5, Hco-avr-14, Hco-pgp-9 and Hco-dyf-7 showed no evidence of introgression in either backcross. Hco-lgc-37 and Hco-pgp-2 showed only weak evidence of introgression in the MHco3/4 backcross but not in the MHco3/10 backcross. Overall, these results suggest that microsatellite marker Hcms8a20, but not the other candidate genes tested, is linked to a major ivermectin resistance locus in the MHco4(WRS) and MHco10(CAVR) strains. This work also emphasises the need for genome-wide approaches to identify mutations responsible for the ivermectin resistance in this parasite.