Introgression of ivermectin resistance genes into a susceptible Haemonchus contortus strain by multiple backcrossing

Chromosomal genome assembly resolves drug resistance loci in the parasitic nematode Teladorsagia circumcincta

The parasitic nematode Teladorsagia circumcincta is one of the most important pathogens of sheep and goats in temperate climates worldwide and can rapidly evolve resistance to drugs used to control it. To understand the genetics of drug resistance, we have generated a highly contiguous genome assembly for the UK T. circumcincta isolate, MTci2. Assembly using PacBio long-reads and Hi-C long-molecule scaffolding together with manual curation resulted in a 573 Mb assembly (N50 = 84 Mb, total scaffolds = 1,286) with five autosomal and one sex-linked chromosomal-scale scaffolds consistent with its karyotype. The genome resource was further improved via annotation of 22,948 genes, with manual curation of over 3,200 of these, resulting in a robust and near complete resource (96.3% complete protein BUSCOs) to support basic and applied research on this important veterinary pathogen. Genome-wide analyses of drug resistance, combining evidence from three distinct experiments, identified selection around known candidate genes for benzimidazole, levamisole and ivermectin resistance, as well as novel regions associated with ivermectin and moxidectin resistance. These insights into contemporary and historic genetic selection further emphasise the importance of contiguous genome assemblies in interpreting genome-wide genetic variation associated with drug resistance and identifying key loci to prioritise in developing diagnostic markers of anthelmintic resistance to support parasite control.

A Comprehensive Review on Haemonchus contortus Excretory and Secretory Proteins (HcESPs): TH-9 stimulated ESPs as a potential candidate for Vaccine Development and Diagnostic Antigen

Haemonchus contortus (Barber pole worm) is one of the dominant helminth parasitic infections in small ruminants which is economically important and causes severe losses in the livestock industry, particularly in tropical and subtropical regions. This parasite resides in the abomasum and is responsible for severe blood loss, leading to anemia, emaciation, hypoproteinemia, weight loss, and potentially death. The economic impact of H. contortus on the livestock industry necessitates effective control measures, including early diagnosis and the development of effective vaccines. H. contortus secretes a variety of excretory and secretory proteins (ESPs), which are glycoproteins that play a crucial role in modulating the host's immune response. These ESPs are not only vital for understanding the immunological interactions between the parasite and the host but also serve as potential diagnostic tools and vaccine candidates. Similar ESPs have been identified in other parasitic species such as Cooperia spp, Ostertagia ostertagia, Teladorsagia circumcincta, Ascaris sum, Schistosoma japonicum, and Echinococcus multilocularis, underscoring their importance in both detection and vaccine development. In addition, there is a lack of highly potential specific proteins which having immunogenic properties that can be used for the accurate, early diagnosis serologically and serve as a potential candidate for the vaccine development against H. contortus. Recent research highlights that TH-9 stimulated proteins from H. contortus are emerging as promising candidates for vaccine development due to their immunomodulatory effects. These proteins have been shown to induce a TH-9 immune response, characterized by increased production of interleukin-9 (IL-9), which is critical for enhancing protective immunity against helminth infections. It is suggested to investigate TH-9 stimulated protein as potential candidates for vaccine development and diagnostic antigen.

Towards precision parasite management for livestock gastrointestinal nematodes in 2030

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

Expression of transporter genes in anthelmintic resistant isolates of Haemonchus contortus

ATP-binding cassette (ABC) transporters, including P-glycoproteins (PGP), have been implicated in drug resistance in different organisms including Haemonchus contortus. This study confirmed the resistance status of H. contortus isolates selected for ivermectin (IVM) and oxfendazole (OXF) resistances using the fecal egg count reduction test and evaluated the gene expression of seven ABC transporters using RT-qPCR for two biological scenarios: the effect of selection for anthelmintic resistance and the effect of drug exposure on gene expression. Gene expression results showed that selection for IVM resistance led to the significant upregulation of Hco-pgp-9a (1.5-fold), Hco-pgp-11 (3-fold) and Hco-haf-9 (1.5-fold) (p < 0.05). Similarly, selection for OXF resistance led to the significant upregulation of Hco-pgp-9a (3-fold), Hco-pgp-11 (4-fold) and Hco-haf-9 (2-fold) when comparing with the unselected ISE isolate (p < 0.05). Exposure of selected isolates to anthelmintics lead to no significant upregulation of the studied transporter genes. We also observed instances where there was strong intragroup variation regarding samples originating from parasites obtained from different individual hosts pointing that the interactions of the animal host with the tested anthelmintics may also play a role in the expression of the studied nematode genes.

Anthelmintic resistance of gastrointestinal nematodes in cattle in Brazil and Argentina - current status and global perspectives

This review outlines the current state of anthelmintic resistance (AHR) of gastrointestinal nematodes (GINs) among cattle in Argentina and Brazil, emphasizing the economic repercussions, animal health and welfare. The analysis explores factors associated with AHR and proposes a potential solution: the use of drug combinations. Both countries are grappling with a severe AHR scenario in cattle, having progressed through incipient, established, and advanced phases, leading to extreme cases of animal mortality due to ineffective control strategies. Genera such as Cooperia and Haemonchus have the highest reports of resistance, with Oesophagostomum radiatum also posing significant problems. While oral benzimidazoles and levamisole remain effective in most herds, moxidectin is entering an advanced resistance phase, and avermectins are increasingly deemed ineffective. The review explores the impact ofclimate, mixed grazing, animal movement and other husbandry practices, and the relationship between ectoparasite control and the emergence of resistant helminths. Notably, the discussion includes the strategic use of drug combinations as a valuable approach to address resistant GINs control in livestock, highlighting its significant potential to mitigate the challenges posed by AHR in the cattle industry of these countries.

The Mitogenome of the Haecon-5 Strain of Haemonchus contortus and a Comparative Analysis of Its Nucleotide Variation with Other Laboratory Strains

Haemonchus contortus (the barber's pole worm)-a highly pathogenic gastric nematode of ruminants-causes significant economic losses in the livestock industry worldwide. H. contortus has become a valuable model organism for both fundamental and applied research (e.g., drug and vaccine discovery) because of the availability of well-defined laboratory strains (e.g., MHco3(ISE).N1 in the UK and Haecon-5 in Australia) and genomic, transcriptomic and proteomic data sets. Many recent investigations have relied heavily on the use of the chromosome-contiguous genome of MHco3(ISE).N1 in the absence of a genome for Haecon-5. However, there has been no genetic comparison of these and other strains to date. Here, we assembled and characterised the mitochondrial genome (14.1 kb) of Haecon-5 and compared it with that of MHco3(ISE).N1 and two other strains (i.e., McMaster and NZ_Hco_NP) from Australasia. We detected 276 synonymous and 25 non-synonymous single nucleotide polymorphisms (SNPs) within Haecon-5. Between the Haecon-5 and MHco3(ISE).N1 strains, we recorded 345 SNPs, 31 of which were non-synonymous and linked to fixed amino acid differences in seven protein-coding genes (nad5, nad6, nad1, atp6, nad2, cytb and nad4) between these strains. Pronounced variation (344 and 435 SNPs) was seen between Haecon-5 and each of the other two strains from Australasia. The question remains as to what impact these mitogenomic mutations might have on the biology and physiology of H. contortus, which warrants exploration. The high degree of mitogenomic variability recorded here among these strains suggests that further work should be undertaken to assess the nature and extent of the nuclear genomic variation within H. contortus.

Next-generation sequencing technologies for helminth diagnostics and surveillance in ruminants: shifting diagnostic barriers

Helminth infections in grazing ruminants are a major issue for livestock farming globally, but are unavoidable in outdoor grazing systems and must be effectively managed to avoid deleterious effects to animal health, and productivity. Next-generation sequencing (NGS) technologies are transforming our understanding of the genetic basis of anthelmintic resistance (AR) and epidemiological studies of ruminant gastrointestinal parasites. They also have the potential to not only help develop and validate molecular diagnostic tests but to be directly used in routine diagnostics integrating species-specific identification and AR into a single test. Here, we review how these developments have opened the pathway for the development of multi-AR and multispecies identification in a single test, with widespread implications for sustainable livestock farming for the future.

Quantifying the fitness effects of resistance alleles with and without anthelmintic selection pressure using Caenorhabditis elegans

Albendazole (a benzimidazole) and ivermectin (a macrocyclic lactone) are the two most commonly co-administered anthelmintic drugs in mass-drug administration programs worldwide. Despite emerging resistance, we do not fully understand the mechanisms of resistance to these drugs nor the consequences of delivering them in combination. Albendazole resistance has primarily been attributed to variation in the drug target, a beta-tubulin gene. Ivermectin targets glutamate-gated chloride channels (GluCls), but it is unknown whether GluCl genes are involved in ivermectin resistance in nature. Using Caenorhabditis elegans, we defined the fitness costs associated with loss of the drug target genes singly or in combinations of the genes that encode GluCl subunits. We quantified the loss-of-function effects on three traits: (i) multi-generational competitive fitness, (ii) fecundity, and (iii) development. In competitive fitness and development assays, we found that a deletion of the beta-tubulin gene ben-1 conferred albendazole resistance, but ivermectin resistance required the loss of two GluCl genes (avr-14 and avr-15). The fecundity assays revealed that loss of ben-1 did not provide any fitness benefit in albendazole conditions and that no GluCl deletion mutants were resistant to ivermectin. Next, we searched for evidence of multi-drug resistance across the three traits. Loss of ben-1 did not confer resistance to ivermectin, nor did loss of any single GluCl subunit or combination confer resistance to albendazole. Finally, we assessed the development of 124 C. elegans wild strains across six benzimidazoles and seven macrocyclic lactones to identify evidence of multi-drug resistance between the two drug classes and found a strong phenotypic correlation within a drug class but not across drug classes. Because each gene affects various aspects of nematode physiology, these results suggest that it is necessary to assess multiple fitness traits to evaluate how each gene contributes to anthelmintic resistance.

Population genomics of helminth parasites

Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.

In vitro culture of parasitic stages of Haemonchus contortus

Haemonchus contortus is a constraint to sheep production. Seeking to reduce the use of hosts and produce parasitic stages in large-scale, a 42-day in vitro culture protocol of H. contortus third-stage larvae was optimized using Dulbecco's modified Eagle's medium (DMEM). In cell-free culture, larvae were maintained at 39.6°C, in acidic media (pH 6.1) for 3 or 6 days with Δ4-dafachronic acid followed by DMEM pH 7.4 supplemented or not with Fildes' reagent. In DMEM pH 7.4 at 37°C, supplementation with Caco-2 cells was compared to Fildes. On Day 14, fourth-stage larvae (L4) development rates in acidic media supplemented (86.8-88.4%) or not (74.4-77.8%) with Fildes and in Caco-2 cell co-culture (92.6%) were similar, and superior to DMEM pH 7.4 with Fildes (0.0%). On Day 21, Caco-2 cell co-culture resulted in higher larvae differentiation (25.0%) and lower degeneration (13.9%) compared to acidic media (1.5-8.1% and 48.6-69.9%, respectively). This is the first report of prolonged in vitro culture of H. contortus larvae using commercial media in co-culture with Caco-2 cells. Although no progression to the adult stage, Caco-2 cell co-culture resulted in morphological differentiation of H. contortus L4 and larval viability for up to 28 days.

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.

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.

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.

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.

Assessment of Avermectins-Induced Toxicity in Animals

Macrocyclic lactones, particularly the avermectins, have completely revolutionized the approaches aimed at control of parasites. These avermectins are the most widely used anti-parasitic drugs in veterinary field with sales exceeding one billion US dollars annually. However, before clinical usage, their safety evaluation in the animals is a major critical factor that must be considered. Many studies have reported the negative effects of avermectins like ivermectin, abamectin, doramectin, and eprinomectin on the host animals. These harmful effects arise from avermectins targeting GABA and glutamate-gated chloride channels present both in the parasites and the host animals. In this review, various modes of avermectins action along with the negative effects on the host like nephrotoxicity, hepatotoxicity, neurotoxicity, reproductive toxicity, and endocrine disruption were discussed in detail. Furthermore, other important issues like ecotoxicity, drug resistance, and drug residues in milk associated with avermectins usage were also discussed, which need special attention.

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.

Evaluation of changes in drug susceptibility and population genetic structure in Haemonchus contortus following worm replacement as a means to reverse the impact of multiple-anthelmintic resistance on a sheep farm

A population of Haemonchus contortus that was highly resistant to benzimidazoles and avermectin/milbemycins with a subpopulation that was resistant to levamisole, was replaced with a susceptible laboratory isolate of H. contortus in a flock of sheep. The anthelmintic susceptibility and population genetics of the newly established population were evaluated for 3.5 years using in vivo, in vitro, and molecular methods. Successful replacement of the resistant population with a susceptible population was confirmed using phenotypic and genotypic measurements; larval development assay indicated full anthelmintic susceptibility; albendazole treatment yielded 98.7% fecal egg count reduction; pyrosequence genotyping of single nucleotide polymorphisms in positions 167 and 200 of the isotype-1 beta tubulin gene were present at 0.0 and 1.7%, respectively; microsatellite genotyping indicated the background haplotype was similar to the susceptible isolate; and haplotypes of the isotype-1 beta tubulin gene were similar to the susceptible isolate. To sustain the susceptibility of the new population, targeted selective treatment was implemented using albendazole. Surprisingly, within 1.5 years post-replacement, the population reverted to a resistant phenotype. Resistance to albendazole, ivermectin, and moxidectin was confirmed via fecal egg count reduction test, larval development assay, and pyrosequencing-based genotyping. Targeted selective treatment was then carried out using levamisole. However, within one year, resistance was detected to levamisole. Population genetics demonstrated a gradual change in the genetic structure of the population until the final population was similar to the initial resistant population. Genetic analyses showed a lack of diversity in the susceptible isolate, suggesting the susceptible isolate had reduced environmental fitness compared to the resistant population, providing a possible explanation for the rapid reversion to resistance. This work demonstrates the power of combining molecular, in vitro, and in vivo assays to study phenotypic and genotypic changes in a field population of nematodes, enabling improved insights into the epidemiology of anthelmintic resistance.

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.

Fasciola Species Introgression: Just a Fluke or Something More?

The threats posed by a range of viral and bacterial zoonotic diseases inevitably receive renewed attention in the wake of global pandemic events due to their overt and devastating impacts on human health and the economy. Parasitic zoonoses, however, many of which affect millions of people each day, are frequently ignored. In the case of fasciolosis, caused by infection with Fasciola hepatica or Fasciola gigantica, this oversight has allowed for the expansion of areas of parasite sympatry and thus increased the incidence of hybridization and possible introgression between the two species. Here we highlight how an increased demand for animal-derived protein, combined with a lack of appropriate tools for detection of these events, is changing the status quo of these zoonotic parasites.

Increased demand for animal-derived protein from Fasciola hepatica-endemic countries has led to a growing number of reports of hybridization between F. hepatica and Fasciola gigantica in Southeast Asia.

Hybridization and eventual introgression have been reported in a range of protozoan, helminth, and arthropod parasites and act as important drivers of evolutionary change and adaptation.

Introgression between Fasciola spp. remains unproven but has potentially serious human and animal health consequences as seen in other parasites.

New tools for the characterization of hybridization and introgression events between Fasciola spp. are needed.

Suppression of hyaluronidase reduces invasion and establishment of Haemonchus contortus larvae in sheep

Haemonchus contortus is a hematophagous endoparasite of small ruminants, which is responsible for huge economic losses in livestock sector. Hyaluronidase produced by infective larvae of H. contortus can degrade hyaluronic acid present in the host’s abomasal tissue. Thus, it facilitates larval tissue invasion and early establishment. We herein explored this ability of hyaluronidase in H. contortus, and tested whether hyaluronidase is utilized as a virulence factor by H. contortus while establishing the infection. We first successfully blocked the hyaluronidase gene in L3 larvae by RNA interference (RNAi), which was subsequently confirmed by qPCR, enzymatic activity, and immunohistochemistry assays. Using these larvae we then conducted in vivo and in vitro assays on sheep to assess the effects of hyaluronidase suppression on larval invasion and establishment of infection. The in vivo assay showed a significant drop in worm burden in siRNA treated group in comparison to control group. During in vitro assay we applied an ovine ex vivo model where siRNA treated group of larvae showed significantly reduced invasion of the abomasal tissue explants as compared to control group. These findings indicate that hyaluronidase plays a key role in host’s tissue invasion and larval establishment, and it is used as a virulence factor by H. contortus while establishing the infection. As an invasive virulence molecule, its functional research is thus conducive to the prevention of haemonchosis.

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.

Recombinant elongation factor 1 alpha of Haemonchus contortus affects the functions of goat PBMCs

Excretory/secretory proteins of Haemonchus contortus (HcESPs) intermingle comprehensively with host immune cells and modulate host immune responses. In this study, H contortus ES antigen named as elongation factor 1 alpha (HcEF‐1α) was cloned and expressed. The influences of recombinant HcEF‐1α on multiple functions of goat peripheral blood mononuclear cells (PBMCs) were observed in vitro. Immunoblot analysis revealed that rHcEF‐1α was recognized by the serum of goat infected with H contortus. Immunofluorescence analysis indicated that rHcEF‐1α was bound on surface of PBMCs. Moreover, the productions of IL‐4, TGF‐β1, IFN‐γ and IL‐17 of cells were significantly modulated by the incubation with rHcEF‐1α. The production of interleukin IL‐10 was decreased. Cell migration, cell proliferation and cell apoptosis were significantly increased; however, nitric oxide production (NO) was significantly decreased. The MHC II molecule expression of cells incubated with rHcEF‐1α was increased significantly, whereas MHC‐I was not changed as compared to the control groups (PBS control and pET32a). These findings indicated that rHcEF‐1α protein might play essential roles in functional regulations of HcESPs on goat PBMC and mediate the immune responses of the host during host‐parasite relationship.

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

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

Extreme-QTL mapping of monepantel resistance in Haemonchus contortus

Background

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

Methods

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

Results

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

Conclusions

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

Electronic supplementary material

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

Microsatellite analysis reveals extensive gene flow, and lack of population structure in the farm populations of Haemonchus contortus in northern China

The parasitic nematode Haemonchus contortus is economically an important parasite of small ruminants across the globe. China is the world's largest producer, consumer, and importer of mutton. With ubiquitous distribution across the country H. contortus is one of the potential candidates to cause huge economic losses to small ruminant farming industry in China. We herein investigated genetic diversity and population structure of six farm populations of H. contortus in northern China, and also compared them to H. contortus isolates from UK and Australia. We first prepared individual DNA samples from 240 adult worms, and generated genotyping data using eight microsatellite markers. Obtained data was then subjected to allelic frequency and population genetic analyses. The overall allelic richness (mean/locus/pop = 7.375 ± 0.844-10.125 ± 1.109), and expected heterozygosity (mean/locus/pop = 0.646 ± 0.040-0.735 ± 0.025) indicated high degree of population genetic variation across the Chinese isolates. Low level of genetic differentiation (Fst = 0.010-0.066) was observed across all the populations. AMOVA results showed high level of variation (93%) within the populations. PCA analysis revealed mixed clustering of all the populations with no visible geographical sub-structuring. Finally the population admixture analysis resulted in extensive admixing of genotypes across all the populations. With these findings we conclude that there is no obvious population genetic structure with extensive gene flow across all the farm populations of H. contortus in northern China.

Mating barriers between genetically divergent strains of the parasitic nematode Haemonchus contortus suggest incipient speciation

Haemonchus contortus, in common with many nematode species, has extremely high levels of genetic variation within and between field populations derived from distant geographical locations. MHco10(CAVR), MHco3(ISE) and MHco4(WRS) are genetically divergent H. contortus strains, originally derived from Australia, Kenya and South Africa, respectively, that have been maintained by numerous rounds of in vivo experimental infection of sheep. In order to explore potential pre-zygotic competition or post-zygotic incompatibility between the strains, we have investigated the ability of MHco10(CAVR) to interbreed with either MHco3(ISE) or MHco4(WRS) during dual strain co-infections. Sheep were experimentally co-infected with 4000 infective larvae (L₃) per os of the MHco10(CAVR) strain and an equal number of either the MHco3(ISE) or the MHco4(WRS) strain L₃. The adult worm establishement rates and the proportions of F₁ progeny resulting from intra- and inter-strain mating events were determined by admixture analysis of microsatellite multi-locus genotypes. Although there was no difference in adult worm establishment rates, the proportions of F₁ progeny of both the MHco10(CAVR) × MHco3(ISE) and MHco10(CAVR) × MHco4(WRS) dual strain co-infections departed from Mendelian expectations. The proportions of inter-strain hybrid F₁ progeny were lower than the expected 50%, suggesting either pre-zygotic competition or post-zygotic incompatibility between the co-infecting strains. To investigate this further, both eggs and hatched L₁ of broods from single adult female worms recovered from each dual co-infection were genotyped. Unhatched eggs from the broods revealed no inter-strain hybrid genotype deficit, suggesting there is no pre-zygotic competition between the strains. In contrast, there was a deficit in L₁ inter-strain hybrid genotypes in the broods derived from MHco3(ISE) or MHco4(WRS) maternal parents, but not from MHco10(CAVR) maternal parents. This suggests that hybrid progeny of MHco10(CAVR) paternal parents have reduced post-zygotic development and/or viability consistent with incipient speciation of the MHco10(CAVR) strain. The presence of mating barriers between allopatric H. contortus strains has important implications for parasite ecology, including the ability of newly introduced anthelmintic-resistant parasite populations to compete and interbreed with populations already established in a region.

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.

The transcription factor NHR-8: A new target to increase ivermectin efficacy in nematodes

Resistance to the anthelmintic macrocyclic lactone ivermectin (IVM) has a great impact on the control of parasitic nematodes. The mechanisms by which nematodes adapt to IVM remain to be deciphered. We have identified NHR-8, a nuclear hormone receptor involved in the xenobiotic response in Caenorhabditis elegans, as a new regulator of tolerance to IVM. Loss-of-function nhr-8(ok186) C. elegans mutants subjected to larval development assays and electropharyngeogram measurements, displayed hypersensitivity to IVM, and silencing of nhr-8 in IVM-resistant worms increased IVM efficacy. In addition, compared to wild-type worms, nhr-8 mutants under IVM selection pressure failed to acquire tolerance to the drug. In addition, IVM-hypersensitive nhr-8(ok186) worms displayed low transcript levels of several genes from the xenobiotic detoxification network and a concomitant low Pgp-mediated drug efflux activity. Interestingly, some pgp and cyp genes known to impact IVM tolerance in many nematode species, were down regulated in nhr-8 mutants and inversely upregulated in IVM-resistant worms. Moreover, pgp-6 overexpression in nhr-8(ok186) C. elegans increased tolerance to IVM. Importantly, NHR-8 function was rescued in nhr-8(ok186) C. elegans with the homolog of the parasitic nematode Haemonchus contortus, and silencing of Hco-nhr-8 by RNAi on L2 H. contortus larvae increased IVM susceptibility in both susceptible and resistant H. contortus isolates. Thus, our data show that NHR-8 controls the tolerance and development of resistance to IVM in C. elegans and the molecular basis for this relates to the NHR-8-mediated upregulation of IVM detoxification genes. Since our results show that Hco-nhr-8 functions similarly to Cel-nhr-8, this study helps to better understand mechanisms underlying failure in drug efficacy and open perspectives in finding new compounds with NHR-8 antagonist activity to potentiate IVM efficacy.

IVM is the most important broad-spectrum deworming drug used today but resistance to this drug has appeared in parasites of both animals and humans. This seriously jeopardizes the success of current parasite control. Preserving IVM efficacy is a public health issue, whose outcome depends on the understanding of the molecular basis of selection for resistance to these drugs. We unambiguously show that the nuclear hormone receptor NHR-8, is crucial for protection of the nematode model Caenorhabditis elegans against IVM toxicity. Worms deficient in NHR-8 are hypersensitive to IVM and fail to become resistant to IVM under drug pressure. NHR-8 functions in the parasitic nematode of ruminants Haemonchus contortus and similar mechanisms could occur in other target pathogens. By controlling the xenobiotic detoxification network, NHR-8 may contribute to the biotransformation and elimination of IVM and help to desensitize the worm to the drug. This provides novel molecular targets involved in IVM drug tolerance in parasitic nematodes. Such findings could be exploited for targeted therapeutic intervention to treat parasitic nematode infections and delay the process of resistance development to anthelmintic drugs.

Extreme allelic heterogeneity at a Caenorhabditis elegans beta-tubulin locus explains natural resistance to benzimidazoles

Benzimidazoles (BZ) are essential components of the limited chemotherapeutic arsenal available to control the global burden of parasitic nematodes. The emerging threat of BZ resistance among multiple nematode species necessitates the development of novel strategies to identify genetic and molecular mechanisms underlying this resistance. All detection of parasitic helminth resistance to BZ is focused on the genotyping of three variant sites in the orthologs of the β-tubulin gene found to confer resistance in the free-living nematode Caenorhabditis elegans. Because of the limitations of laboratory and field experiments in parasitic nematodes, it is difficult to look beyond these three sites to identify additional mechanisms that might contribute to BZ resistance in the field. Here, we took an unbiased genome-wide mapping approach in the free-living nematode species C. elegans to identify the genetic underpinnings of natural resistance to the commonly used BZ, albendazole (ABZ). We found a wide range of natural variation in ABZ resistance in natural C. elegans populations. In agreement with known mechanisms of BZ resistance in parasites, we found that a majority of the variation in ABZ resistance among wild C. elegans strains is caused by variation in the β-tubulin gene ben-1. This result shows empirically that resistance to ABZ naturally exists and segregates within the C. elegans population, suggesting that selection in natural niches could enrich for resistant alleles. We identified 25 distinct ben-1 alleles that are segregating at low frequencies within the C. elegans population, including many novel molecular variants. Population genetic analyses indicate that ben-1 variation arose multiple times during the evolutionary history of C. elegans and provide evidence that these alleles likely occurred recently because of local selective pressures. Additionally, we find purifying selection at all five β-tubulin genes, despite predicted loss-of-function variants in ben-1, indicating that BZ resistance in natural niches is a stronger selective pressure than loss of one β-tubulin gene. Furthermore, we used genome-editing to show that the most common parasitic nematode β-tubulin allele that confers BZ resistance, F200Y, confers resistance in C. elegans. Importantly, we identified a novel genomic region that is correlated with ABZ resistance in the C. elegans population but independent of ben-1 and the other β-tubulin loci, suggesting that there are multiple mechanisms underlying BZ resistance. Taken together, our results establish a population-level resource of nematode natural diversity as an important model for the study of mechanisms that give rise to BZ resistance.

Identification and preliminary characterization of Hc-clec-160, a novel C-type lectin domain-containing gene of the strongylid nematode Haemonchus contortus

Background

The strongylid parasite Haemonchus contortus causes severe anemia in domestic animals worldwide. Effective preventive and therapeutical agents are lacking, because of drug resistance and that little is known about the molecular mechanism of the interaction between H. contortus and host cells.

Methods

A new gene, Hc-clec-160, was discovered with RT-PCR. Transcriptional levels of Hc-clec-160 and Ce-clec-160 throughout different growth phases of corresponding nematodes were assayed by qPCR. Immunofluorescence staining of paraffin section were performed to determine the protein localization in adult worms of H. contortus. To monitor the promoter capacity of the 5'-flanking region of Ce-clec-160, micro-injection was used. Overexpression and RNAi constructs was carried out in the N2 strain of Caenorhabditis elegans to find out the gene function of Hc-clec-160.

Results

The full-length cDNA of 1224 bp of Hc-clec-160 was cloned by RT-PCR. The corresponding gene contained twelve exons. Its transcripts peaked in male adult worms. Hc-CLEC-160 was predicted to have a Willebrand factor type A (vWA) domain and a C-type lectin domain. The proteins were not detected by expression in C. elegans or paraffin section experiments in adult of H. contortus. Knockdown of Ce-clec-160 expression in C. elegans by RNAi resulted in shortened body length and decreased brood size.

Conclusions

In this experiment, a new gene Hc-clec-160 was obtained in H. contortus and its function was addressed using a model organism: C. elegans. Our study showed that Hc-clec-160 possesses characteristics similar to those of Ce-clec-160 and plays an important role in the growth and reproduction of this parasitic nematode.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3005-3) contains supplementary material, which is available to authorized users.

Probing the opportunities for designing anthelmintic leads by sub-structural topology-based QSAR modelling

A quantitative structure-activity (QSAR) model has been developed for enriched tubulin inhibitors, which were retrieved from sequence similarity searches and applicability domain analysis. Using partial least square (PLS) method and leave-one-out (LOO) validation approach, the model was generated with the correlation statistics of [Formula: see text] and [Formula: see text] of 0.68 and 0.69, respectively. The present study indicates that topological descriptors, viz. BIC, CH_3_C, IC, JX and Kappa_2 correlate well with biological activity. ADME and toxicity (or ADME/T) assessment showed that out of 260 molecules, 255 molecules successfully passed the ADME/T assessment test, wherein the drug-likeness attributes were exhibited. These results showed that topological indices and the colchicine binding domain directly influence the aetiology of helminthic infections. Further, we anticipate that our model can be applied for guiding and designing potential anthelmintic inhibitors.

Discovery of genomic intervals that underlie nematode responses to benzimidazoles

Parasitic nematodes impose a debilitating health and economic burden across much of the world. Nematode resistance to anthelmintic drugs threatens parasite control efforts in both human and veterinary medicine. Despite this threat, the genetic landscape of potential resistance mechanisms to these critical drugs remains largely unexplored. Here, we exploit natural variation in the model nematodes Caenorhabditis elegans and Caenorhabditis briggsae to discover quantitative trait loci (QTL) that control sensitivity to benzimidazoles widely used in human and animal medicine. High-throughput phenotyping of albendazole, fenbendazole, mebendazole, and thiabendazole responses in panels of recombinant lines led to the discovery of over 15 QTL in C. elegans and four QTL in C. briggsae associated with divergent responses to these anthelmintics. Many of these QTL are conserved across benzimidazole derivatives, but others show drug and dose specificity. We used near-isogenic lines to recapitulate and narrow the C. elegans albendazole QTL of largest effect and identified candidate variants correlated with the resistance phenotype. These QTL do not overlap with known benzimidazole target resistance genes from parasitic nematodes and present specific new leads for the discovery of novel mechanisms of nematode benzimidazole resistance. Analyses of orthologous genes reveal conservation of candidate benzimidazole resistance genes in medically important parasitic nematodes. These data provide a basis for extending these approaches to other anthelmintic drug classes and a pathway towards validating new markers for anthelmintic resistance that can be deployed to improve parasite disease control.

The treatment of roundworm (nematode) infections in both humans and animals relies on a small number of anti-parasitic drugs. Resistance to these drugs has appeared in veterinary parasite populations and is a growing concern in human medicine. A better understanding of the genetic basis for parasite drug resistance can be used to help maintain the effectiveness of anti-parasitic drugs and to slow or to prevent the spread of drug resistance in parasite populations. This goal is hampered by the experimental intractability of nematode parasites. Here, we use non-parasitic model nematodes to systematically explore responses to the critical benzimidazole class of anti-parasitic compounds. Using a quantitative genetics approach, we discovered unique genomic intervals that control drug effects, and we identified differences in the effects of these intervals across compounds and doses. We were able to narrow a major-effect genomic region associated with albendazole resistance and to establish that candidate genes discovered in our genetic mappings are largely conserved in important human and animal parasites. This work provides new leads for understanding parasite drug resistance and contributes a powerful template that can be extended to other anti-parasitic drug classes.

A Genome Resequencing-Based Genetic Map Reveals the Recombination Landscape of an Outbred Parasitic Nematode in the Presence of Polyploidy and Polyandry

The parasitic nematode Haemonchus contortus is an economically and clinically important pathogen of small ruminants, and a model system for understanding the mechanisms and evolution of traits such as anthelmintic resistance. Anthelmintic resistance is widespread and is a major threat to the sustainability of livestock agriculture globally; however, little is known about the genome architecture and parameters such as recombination that will ultimately influence the rate at which resistance may evolve and spread. Here, we performed a genetic cross between two divergent strains of H. contortus, and subsequently used whole-genome resequencing of a female worm and her brood to identify the distribution of genome-wide variation that characterizes these strains. Using a novel bioinformatic approach to identify variants that segregate as expected in a pseudotestcross, we characterized linkage groups and estimated genetic distances between markers to generate a chromosome-scale F1 genetic map. We exploited this map to reveal the recombination landscape, the first for any helminth species, demonstrating extensive variation in recombination rate within and between chromosomes. Analyses of these data also revealed the extent of polyandry, whereby at least eight males were found to have contributed to the genetic variation of the progeny analyzed. Triploid offspring were also identified, which we hypothesize are the result of nondisjunction during female meiosis or polyspermy. These results expand our knowledge of the genetics of parasitic helminths and the unusual life-history of H. contortus, and enhance ongoing efforts to understand the genetic basis of resistance to the drugs used to control these worms and for related species that infect livestock and humans throughout the world. This study also demonstrates the feasibility of using whole-genome resequencing data to directly construct a genetic map in a single generation cross from a noninbred nonmodel organism with a complex lifecycle.

Metabolism of albendazole, ricobendazole and flubendazole in Haemonchus contortus adults: Sex differences, resistance-related differences and the identification of new metabolites

Haemonchus contortus (family Trichostrongylidae, Nematoda), a hematophagous gastrointestinal parasite found in small ruminants, has a great ability to develop resistance to anthelmintic drugs. We studied the biotransformation of the three benzimidazole anthelmintics: albendazole (ABZ), ricobendazole (albendazole S-oxide; RCB) and flubendazole (FLU) in females and males of H. contortus in both a susceptible ISE strain and resistant IRE strain. The ex vivo cultivation of living nematodes in culture medium with or without the anthelmintics was used. Ultrasensitive UHPLC/MS/MS analysis revealed 9, 7 and 12 metabolites of ABZ, RCB and FLU, respectively, with most of these metabolites now described in the present study for the first time in H. contortus. The structure of certain metabolites shows the presence of biotransformation reactions not previously reported in nematodes. There were significant qualitative and semi-quantitative differences in the metabolites formed by male and female worms. In most cases, females metabolized drugs more extensively than males. Adults of the IRE strain were able to form many more metabolites of all the drugs than adults of the ISE strain. Some metabolites were even found only in adults of the IRE strain. These findings suggest that increased drug metabolism may play a role in resistance to benzimidazole drugs in H. contortus.

Recombinant protein of Haemonchus contortus small GTPase ADP-ribosylation factor 1 (HcARF1) modulate the cell mediated immune response in vitro

ADP-ribosylation factors (ARFs) are members of the Ras-related small GTPase family involved in the vesicular trafficking regulation. Immunomodulatory effects of these proteinson host cell arenot being addressed yet. H. contortus small GTPase ADP-ribosylation 1 gene (HcARF1) was cloned and recombinant protein of HcARF1 (rHcARF1) was successfully expressed in Escherichia coli. Binding activity of rHcARF1 to goat PBMCs was confirmed by immunofluorescence assay (IFA) and its immunomudulatory effects on cytokine secretion, cell proliferation, cell migration and nitric oxide production (NO) were observed by co-incubation of rHcARF1. IFA results revealed that rHcARF1 could bind to the PBMCs. The interaction of rHcARF1 modulated the cytokine production, the production of IL-4, IL-10 and IL-17 was increased in a dose dependent manner, however, the IFN-γ production was significantly decreased. Cell migration and NO production were significantly increased by rHcARF1, whereas, rHcARF1 treatment significantly suppressed the proliferation of the PBMC in a dose dependent manner. Our findings showed that the rHcARF1 play important roles on the goat PBMCs.

Increased Expression of a MicroRNA Correlates with Anthelmintic Resistance in Parasitic Nematodes

Resistance to anthelmintic drugs is a major problem in the global fight against parasitic nematodes infecting humans and animals. While previous studies have identified mutations in drug target genes in resistant parasites, changes in the expression levels of both targets and transporters have also been reported. The mechanisms underlying these changes in gene expression are unresolved. Here, we take a novel approach to this problem by investigating the role of small regulatory RNAs in drug resistant strains of the important parasite Haemonchus contortus. microRNAs (miRNAs) are small (22 nt) non-coding RNAs that regulate gene expression by binding predominantly to the 3′ UTR of mRNAs. Changes in miRNA expression have been implicated in drug resistance in a variety of tumor cells. In this study, we focused on two geographically distinct ivermectin resistant strains of H. contortus and two lines generated by multiple rounds of backcrossing between susceptible and resistant parents, with ivermectin selection. All four resistant strains showed significantly increased expression of a single miRNA, hco-miR-9551, compared to the susceptible strain. This same miRNA is also upregulated in a multi-drug-resistant strain of the related nematode Teladorsagia circumcincta. hco-miR-9551 is enriched in female worms, is likely to be located on the X chromosome and is restricted to clade V parasitic nematodes. Genes containing predicted binding sites for hco-miR-9551 were identified computationally and refined based on differential expression in a transcriptomic dataset prepared from the same drug resistant and susceptible strains. This analysis identified three putative target mRNAs, one of which, a CHAC domain containing protein, is located in a region of the H. contortus genome introgressed from the resistant parent. hco-miR-9551 was shown to interact with the 3′ UTR of this gene by dual luciferase assay. This study is the first to suggest a role for miRNAs and the genes they regulate in drug resistant parasitic nematodes. miR-9551 also has potential as a biomarker of resistance in different nematode species.

A method for single pair mating in an obligate parasitic nematode

Parasitic nematode species have extremely high levels of genetic diversity, presenting a number of experimental challenges for genomic and genetic work. Consequently, there is a need to develop inbred laboratory strains with reduced levels of polymorphism. The most efficient approach to inbred line development is single pair mating, but this is challenging for obligate parasites where the adult sexual reproductive stages are inside the host, and thus difficult to experimentally manipulate. This paper describes a successful approach to single pair mating of a parasitic nematode, Haemonchus contortus. The method allows for polyandrous mating behaviour and involves the surgical transplantation of a single adult male worm with multiple immature adult females directly into the sheep abomasum. We used a panel of microsatellite markers to monitor and validate the single pair mating crosses and to ensure that the genotypes of progeny and subsequent filial generations were consistent with those expected from a mating between a single female parent of known genotype and a single male parent of unknown genotype. We have established two inbred lines that both show a significant overall reduction in genetic diversity based on microsatellite genotyping and genome-wide single nucleotide polymorphism. There was an approximately 50% reduction in heterozygous SNP sites across the genome in the MHco3.N1 line compared with the MoHco3(ISE) parental strain. The MHco3.N1 inbred line has subsequently been used to provide DNA template for whole genome sequencing of H. contortus. This work provides proof of concept and methodologies for forward genetic analysis of obligate parasitic nematodes.

Recent Research Progress in China on Haemonchus contortus

Haemonchus contortus is one of the most important parasites of ruminants with worldwide distribution that can bring huge economic losses to the breeding industry of cattle, sheep, and goats. In recent 20 years, studies on H. contortus in China mainly focused on the epidemiology, population genetics, anthelmintic resistance, structural and functional studies of important genes regulating the development of this parasite, interaction between parasite molecules and host cells and vaccine development against haemonchosis, and achieved good progress. However, there is no systematic review about the studies by Chinese researchers on H. contortus in China. The purpose of this review is to bring together the findings from the studies on H. contortus in China in order to obtain the knowledge gained from the recent studies in China and provide foundation for identifying future research directions to establish novel diagnostic methods, discover new drug targets and vaccine candidates for use in preventing and controlling H. contortus in China.

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.

Ivermectin - Old Drug, New Tricks?

Ivermectin is one of the most important drugs in veterinary and human medicine for the control of parasitic infection and was the joint focus of the 2015 Nobel Prize in Physiology or Medicine, some 35 years after its remarkable discovery. Although best described for its activity on glutamate-gated chloride channels in parasitic nematodes, understanding of its mode of action remains incomplete. In the field of veterinary medicine, resistance to ivermectin is now widespread, but the mechanisms underlying resistance are unresolved. Here we discuss the history of this versatile drug and its use in global health. Based on recent studies in a variety of systems, we question whether ivermectin could have additional modes of action on parasitic nematodes.

Ligand-gated ion channels, particularly glutamate-gated chloride channels, are well characterised as the targets of IVM in nematodes and insects.

Nematode genomes are helping to cast light on the diversity of ion-channel subunits in different parasite species of human and veterinary importance.

Resistance to IVM is an increasing problem in the control of parasitic nematodes, and resolving the mechanisms is an important research priority.

Recent studies in other biological systems suggest that IVM can affect a number of additional pathways.

IVM may have novel applications in the treatment and control of important human diseases.

Genome-wide SNP analysis using 2b-RAD sequencing identifies the candidate genes putatively associated with resistance to ivermectin in Haemonchus contortus

Background

The excessive and uncontrolled use of anthelmintics, e.g. ivermectin (IVM) for the treatment of livestock parasites has led to widespread resistance in gastrointestinal nematodes, such as Haemonchus contortus. There is an urgent need for better management of drug-use in nematode control and development of novel anthelmintics. Discovery and identification of anthelmintic resistance-associate molecules/markers can provide a basis for rational anthelmintics-use and development of novel drugs. Recent studies have shown that ivermectin resistance in H. contortus is likely to be multi-genic in nature except for several genes coding for IVM target and efflux pump. However, no other IVM resistance-associated genes were characterized by conventional methods or strategies. In the present study we adopted a new strategy, i.e. using genome-wide single nucleotide polymorphism (SNP) analysis based on 2b-RAD sequencing, for discovering SNPs markers across the genomes in both IVM susceptible and resistant isolates of H. contortus and identifying potential IVM resistance-associated genes.

Results

We discovered 2962 and 2667 SNPs within both susceptible and resistant strains of H. contortus, respectively. A relative lower and similar genetic variations were observed within both resistant and susceptible strains (average π values were equal to 0.1883 and 0.1953, respectively); whereas a high genetic variation was found across both strains (average π value was equal to 0.3899). A significant differentiation across 2b-RAD tags nucleotide sites was also observed between the two strains (average F_ST value was equal to 0.3076); the larger differences in average F_ST were observed at SNPs loci between coding and noncoding (including intronic) regions. Comparison between resistant and susceptible strains revealed that 208 SNPs loci exhibited significantly elevated F_ST values, 24 SNPs of those loci were located in the CDS regions of the nine genes and were likely to have signature of IVM directional selection. Seven of the nine candidate genes were predicted to code for some functional proteins such as potential IVM target and/or efflux pump proteins, component proteins of receptor complex in membrane on neuromuscular cells, and transcriptional regulation proteins. Those genes might be involved in resistance to IVM.

Conclusions

Our data suggest that candidate genes putatively associated with resistance to IVM in H. contortus may be identified by genome-wide SNP analysis using 2b-RAD sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1959-6) contains supplementary material, which is available to authorized users.

Recombinant protein of Haemonchus contortus 14-3-3 isoform 2 (rHcftt-2) decreased the production of IL-4 and suppressed the proliferation of goat PBMCs in vitro

14-3-3 proteins have been found to be an excreted/secreted antigen and assumed to be released into the host-parasite interface and described in several unicellular and multicellular parasites. However, little is known about the immunomodulatory effects of H. controtus 14-3-3 protein on host cell. In present study, 14-3-3 isoform 2 gene, designated as Hcftt-2, was amplified by reverse transcription-polymerase chain reaction (RT-PCR) from the adult H. contortus cDNA and cloned into expression plasmid pET32a (+) and expression of the recombinant protein (rHcftt-2) was induced by IPTG. Binding activity of rHcftt-2 to goat peripheral blood mononuclear cells (PBMCs) was confirmed by immunofluorescence assay (IFA) and modulatory effects on cytokine production, cell proliferation, cell migration and nitric oxide (NO) production were observed by co-incubation of rHcftt-2 with goat PBMCs. Sequence analysis showed that it had significant homology with the known 14-3-3 protein isoform 2. Results of IFA revealed that, the rHcftt-2 was bound to the cell surface. We found that, the productions of IL10, IL-17, IFN-γ and cell migration of PBMCs were increased after the cells were incubated with rHCftt-2. However, the productions of IL-4, NO and cell proliferation of the PBMCs were significantly decreased in dose depended manner. Our results showed that the Hcftt-2 played important suppressive regulatory effects on the goat PBMCs.

Proteomic Analysis of the Excretory and Secretory Proteins of Haemonchus contortus (HcESP) Binding to Goat PBMCs In Vivo Revealed Stage-Specific Binding Profiles

Haemonchus contortus is a parasitic gastrointestinal nematode, and its excretory and secretory products (HcESPs) interact extensively with the host cells. In this study, we report the interaction of proteins from HcESPs at different developmental stages to goat peripheral blood mononuclear cells (PBMCs) in vivo using liquid chromatography-tandem mass spectrometry. A total of 407 HcESPs that interacted with goat PBMCs at different time points were identified from a H. contortus protein database using SEQUEST searches. The L4 and L5 stages of H. contortus represented a higher proportion of the identified proteins compared with the early and late adult stages. Both stage-specific interacting proteins and proteins that were common to multiple stages were identified. Forty-seven interacting proteins were shared among all stages. The gene ontology (GO) distributions of the identified goat PBMC-interacting proteins were nearly identical among all developmental stages, with high representation of binding and catalytic activity. Cellular, metabolic and single-organism processes were also annotated as major biological processes, but interestingly, more proteins were annotated as localization processes at the L5 stage than at the L4 and adult stages. Based on the clustering of homologous proteins, we improved the functional annotations of un-annotated proteins identified at different developmental stages. Some unnamed H. contortus ATP-binding cassette proteins, including ADP-ribosylation factor and P-glycoprotein-9, were identified by STRING protein clustering analysis.

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.