Inefficacy of ivermectin and moxidectin treatments against Dictyocaulus viviparus in dairy calves.

BACKGROUND: The bovine lungworm Dictyocaulus viviparus negatively impacts bovine health and leads to substantial economic losses. Lungworm infections can be difficult to manage due to the unpredictable and severe nature of clinical outbreaks. Despite the widespread use of macrocyclic lactones (MLs) in grazing cattle in the UK, there have been no confirmed reports of resistant lungworms to date, with only one case of anthelmintic-resistant (ML) lungworm confirmed worldwide. METHODS: Lungworm Baermann filtrations were conducted on first-season grazing dairy calves as part of a wider study investigating anthelmintic resistance in gastrointestinal nematodes in Scotland using the faecal egg count reduction test. RESULTS: Clinical signs and significant numbers of lungworm larvae in faeces were observed after treatment with either ivermectin or moxidectin. LIMITATIONS: There are no established guidelines for the diagnosis of resistant lungworms in the field. Currently, resistance can only be diagnosed after a controlled efficacy test has been conducted. This limits the conclusions that can be drawn; however, they are highly suggestive of resistance. CONCLUSION: This short report describes the inefficacy of ivermectin and moxidectin against D. viviparus and is highly suggestive of ML resistance.

Functional validation of novel levamisole resistance marker S168T in Haemonchus contortus.

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

A mixed amplicon metabarcoding and sequencing approach for surveillance of drug resistance to levamisole and benzimidazole in Haemonchus spp.

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

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

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

New technologies to study helminth development and host-parasite interactions.

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

Moxidectin use in Scottish sheep flocks suggests a need for clearer product labelling and communication of updated SCOPS guidelines.

BACKGROUND: Guidelines for sustainable use of moxidectin were established in 2020. This study aimed to identify how Scottish sheep farmers are using this key endectocide and estimate its effectiveness against gastrointestinal nematodes. METHODS: Questionnaires were distributed to sheep farmers across Scotland, and analysis focused on moxidectin use in relation to Sustainable Control of Parasites in Sheep (SCOPS) guidelines. Farmers using moxidectin in their flock volunteered to submit post-treatment sheep faecal samples, which were analysed for the presence of gastrointestinal nematodes using faecal egg counts with polymerase chain reaction to determine species. RESULTS: Despite 70% of farmers using moxidectin in 2020, knowledge levels varied: 24% of farmers included other anthelmintics when asked about moxidectin use. Moxidectin was used for a wide variety of reasons, and most farmers did not consistently follow SCOPS guidelines. Despite only 2 of 76 farmers reporting failure of moxidectin treatment, gastrointestinal nematodes were found following moxidectin treatment on five out of six farms tested and included Teladorsagia circumcincta, Cooperia curticei, Haemonchus contortus and Nematodirus sp. CONCLUSION: Findings from this project indicate the need for improved anthelmintic product labelling and farmer support to encourage sustainable use. The presence of nematodes in treated animals is suggestive of anthelmintic resistance.

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.

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

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

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

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

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

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

Small RNAs in parasitic nematodes - forms and functions.

Small RNAs are important regulators of gene expression. They were first identified in Caenorhabditis elegans, but it is now apparent that the main small RNA silencing pathways are functionally conserved across diverse organisms. Availability of genome data for an increasing number of parasitic nematodes has enabled bioinformatic identification of small RNA sequences. Expression of these in different lifecycle stages is revealed by small RNA sequencing and microarray analysis. In this review we describe what is known of the three main small RNA classes in parasitic nematodes - microRNAs (miRNAs), Piwi-interacting RNAs (piRNAs) and small interfering RNAs (siRNAs) - and their proposed functions. miRNAs regulate development in C. elegans and the temporal expression of parasitic nematode miRNAs suggest modulation of target gene levels as parasites develop within the host. miRNAs are also present in extracellular vesicles released by nematodes in vitro, and in plasma from infected hosts, suggesting potential regulation of host gene expression. Roles of piRNAs and siRNAs in suppressing target genes, including transposable elements, are also reviewed. Recent successes in RNAi-mediated gene silencing, and application of small RNA inhibitors and mimics will continue to advance understanding of small RNA functions within the parasite and at the host-parasite interface.

Profiling microRNAs through development of the parasitic nematode Haemonchus identifies nematode-specific miRNAs that suppress larval development.

Parasitic nematodes transition between dramatically different free-living and parasitic stages, with correctly timed development and migration crucial to successful completion of their lifecycle. However little is known of the mechanisms controlling these transitions. microRNAs (miRNAs) negatively regulate gene expression post-transcriptionally and regulate development of diverse organisms. Here we used microarrays to determine the expression profile of miRNAs through development and in gut tissue of the pathogenic nematode Haemonchus contortus. Two miRNAs, mir-228 and mir-235, were enriched in infective L3 larvae, an arrested stage analogous to Caenorhabditis elegans dauer larvae. We hypothesized that these miRNAs may suppress development and maintain arrest. Consistent with this, inhibitors of these miRNAs promoted H. contortus development from L3 to L4 stage, while genetic deletion of C. elegans homologous miRNAs reduced dauer arrest. Epistasis studies with C. elegans daf-2 mutants showed that mir-228 and mir-235 synergise with FOXO transcription factor DAF-16 in the insulin signaling pathway. Target prediction suggests that these miRNAs suppress metabolic and transcription factor activity required for development. Our results provide novel insight into the expression and functions of specific miRNAs in regulating nematode development and identify miRNAs and their target genes as potential therapeutic targets to limit parasite survival within the host.

Evaluation of DNA Extraction Methods on Individual Helminth Egg and Larval Stages for Whole-Genome Sequencing.

Whole-genome sequencing is being rapidly applied to the study of helminth genomes, including de novo genome assembly, population genetics, and diagnostic applications. Although late-stage juvenile and adult parasites typically produce sufficient DNA for molecular analyses, these parasitic stages are almost always inaccessible in the live host; immature life stages found in the environment for which samples can be collected non-invasively offer a potential alternative; however, these samples typically yield very low quantities of DNA, can be environmentally resistant, and are susceptible to contamination, often from bacterial or host DNA. Here, we have tested five low-input DNA extraction protocols together with a low-input sequencing library protocol to assess the feasibility of whole-genome sequencing of individual immature helminth samples. These approaches do not use whole-genome amplification, a common but costly approach to increase the yield of low-input samples. We first tested individual parasites from two species spotted onto FTA cards-egg and L1 stages of Haemonchus contortus and miracidia of Schistosoma mansoni-before further testing on an additional five species-Ancylostoma caninum, Ascaridia dissimilis, Dirofilaria immitis, Strongyloides stercoralis, and Trichuris muris-with an optimal protocol. A sixth species-Dracunculus medinensis-was included for comparison. Whole-genome sequencing followed by analyses to determine the proportion of on- and off-target mapping revealed successful sample preparations for six of the eight species tested with variation both between species and between different life stages from some species described. These results demonstrate the feasibility of whole-genome sequencing of individual parasites, and highlight a new avenue toward generating sensitive, specific, and information-rich data for the diagnosis and surveillance of helminths.

Genotypic characterisation of monepantel resistance in historical and newly derived field strains of Teladorsagia circumcincta.

Recent reports of monepantel (MPTL) resistance in UK field isolates of Teladorsagia circumcincta has highlighted the need for a better understanding of the mechanism of MPTL-resistance in order to preserve its anthelmintic efficacy in this economically important species. Nine discrete populations of T. circumcincta were genotypically characterised; three MPTL-susceptible isolates, three experimentally selected MPTL-resistant strains and three field derived populations. Full-length Tci-mptl-1 gene sequences were generated and comparisons between the MPTL-susceptible isolates, MPTL-resistant strains and one field isolate, showed that different putative MPTL-resistance conferring mutations were present in different resistant isolates. Truncated forms of the Tci-mptl-1 gene were also observed. The genetic variability of individual larvae, within and between populations, was examined using microsatellite analyses at 10 'neutral' loci (presumed to be unaffected by MPTL). Results confirmed that there was little background genetic variation between the populations, global FST <0.038. Polymorphisms present in exons 7 and 8 of Tci-mptl-1 enabled genotyping of individual larvae. A reduction in the number of genotypes was observed in all MPTL-resistant strains compared to the MPTL-susceptible strains that they were derived from, suggesting there was purifying selection at Tci-mptl-1 as a result of MPTL-treatment. The potential link between benzimidazole (BZ)-resistance and MPTL-resistance was examined by screening individual larvae for the presence of three SNPs associated with BZ-resistance in the ß-tubulin isotype-1 gene. The majority of larvae were BZ-susceptible homozygotes at positions 167 and 198. Increased heterozygosity at position 200 was observed in the MPTL-resistant strains compared to their respective MPTL-susceptible population. There was no decrease in the occurrence of BZ-resistant genotypes in larvae from each population. These differences, in light of the purifying selection at this locus in all MPTL-resistant isolates, suggests that Tci-mptl-1 confers MPTL-resistance in T. circumcincta, as in Haemonchus contortus, but that different mutations in Tci-mptl-1 can confer resistance in different populations.

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

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

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.

UDP-glycosyltransferase family in Haemonchus contortus: Phylogenetic analysis, constitutive expression, sex-differences and resistance-related differences.

UDP-glycosyltransferases (UGT), catalysing conjugation of UDP-activated sugar donors to small lipophilic chemicals, are widespread in living organisms from bacteria to fungi, plant, or animals. The progress of genome sequencing has enabled an assessment of the UGT multigene family in Haemonchus contortus (family Trichostrongylidae, Nematoda), a hematophagous gastrointestinal parasite of small ruminants. Here we report 32 putative UGT genes divided into 15 UGT families. Phylogenetic analysis in comparison with UGTs from Caenorhabditis elegans, a free-living model nematode, revealed several single member homologues, a lack of the dramatic gene expansion seen in C. elegans, but also several families (UGT365, UGT366, UGT368) expanded in H. contortus only. The assessment of constitutive UGT mRNA expression in H. contortus adults identified significant differences between females and males. In addition, we compared the expression of selected UGTs in the drug-sensitive ISE strain to two benzimidazole-resistant strains, IRE and WR, with different genetic backgrounds. Constitutive expression of UGT368B2 was significantly higher in both resistant strains than in the sensitive strain. As resistant strains were able to deactivate benzimidazole anthelmintics via glycosylation more effectively then the sensitive strain, UGT368B2 enhanced constitutive expression might contribute to drug resistance in H. contortus.

Hidden in plain sight - Multiple resistant species within a strongyle community.

Ovine parasitic gastroenteritis is a complex disease routinely treated using anthelmintics. Although many different strongyle species may contribute to parasitic gastroenteritis, not all are equally pathogenic: in temperate regions, the primary pathogen is Teladorsagia circumcincta. In this study we investigated benzimidazole and ivermectin resistance on a commercial sheep farm in southeast Scotland. We assessed the impact of species diversity on the diagnosis of resistance using the faecal egg count reduction test and in vitro bioassays, and correlated the results with the frequency of benzimidazole resistance-associated genotypes measured in the T. circumcincta population by pyrosequencing of the ß-tubulin isotype-1 gene. Faecal egg count reduction test results showed efficacies of 65% for albendazole and 77% for ivermectin, indicating moderate resistance levels on the farm. However, PCR speciation of the same populations pre- and post-treatment revealed that removal of susceptible species had masked the presence of a highly resistant population of T. circumcincta. Less than 25% of individuals in the pre-treatment populations were T. circumcincta, the remainder consisting of Cooperia curticei, Chabertia ovina, Oesophagostomum venulosum and Trichostrongylus spp. In contrast, post-treatment with albendazole or ivermectin, the majority (88% and 100% respectively) of the populations consisted of T. circumcincta. The egg hatch test for benzimidazole resistance and the larval development test for ivermectin resistance were carried out using eggs obtained from the same populations and the results were broadly consistent with the faecal egg count reduction test. Thirty individual T. circumcincta from each sampling time point were assessed for benzimidazole resistance by pyrosequencing, revealing a high frequency and diversity of resistance-associated mutations, including within the population sampled post-ivermectin treatment. These results highlight the potential diversity of parasite species present on UK farms, and their importance in the diagnosis of anthelmintic resistance. On this particular farm, we demonstrate the presence of a highly dual-resistant population of T. circumcincta, which was strongly selected by treatment with either benzimidazoles or ivermectin, while other potentially less pathogenic species were removed.