Effects of in vitro exposure to ivermectin and levamisole on the expression patterns of ABC transporters in Haemonchus contortus larvae

The mrp-3 gene is involved in haem efflux and detoxification in a blood-feeding nematode

BACKGROUND

Haem is essential but toxic for metazoan organisms. Auxotrophic nematodes can acquire sufficient haem from the environment or their hosts in the meanwhile eliminate or detoxify excessive haem through tightly controlled machinery. In previous work, we reported a role of the unique transporter protein HRG-1 in the haem acquisition and homeostasis of parasitic nematodes. However, little is known about the haem efflux and detoxification via ABC transporters, particularly the multiple drug resistance proteins (MRPs).

RESULTS

Here, we further elucidate that a member of the mrp family (mrp-3) is involved in haem efflux and detoxification in a blood-feeding model gastrointestinal parasite, Haemonchus contortus. This gene is haem-responsive and dominantly expressed in the intestine and inner membrane of the hypodermis of this parasite. RNA interference of mrp-3 resulted in a disturbance of genes (e.g. hrg-1, hrg-2 and gst-1) that are known to be involved in haem homeostasis and an increased formation of haemozoin in the treated larvae and lethality in vitro, particularly when exposed to exogenous haem. Notably, the nuclear hormone receptor NHR-14 appears to be associated the regulation of mrp-3 expression for haem homeostasis and detoxification. Gene knockdown of nhr-14 and/or mrp-3 increases the sensitivity of treated larvae to exogenous haem and consequently a high death rate (> 80%).

CONCLUSIONS

These findings demonstrate that MRP-3 and the associated molecules are essential for haematophagous nematodes, suggesting novel intervention targets for these pathogens in humans and animals.

Comparison of P-glycoprotein gene expression of two Haemonchus contortus isolates from Yucatan, Mexico, with resistant or susceptible phenotype to ivermectin in relation to a susceptible reference strain

The variability in the expression of different P-glycoprotein (P-gp) genes in parasitic nematodes of ruminants such as Haemonchus contortus (Hco-pgp) may be caused by different factors including nematode biology, geographical region and anthelmintic pressure. This study analysed the relative expression level of 10 P-gp genes in two H. contortus (Hco-pgp) field isolates from Yucatan, Mexico: 1) PARAISO (IVM-resistant) and 2) FMVZ-UADY (IVM-susceptible). These isolates were compared with a susceptible reference isolate from Puebla, Mexico, namely "CENID-SAI". In all cases H. contortus adult males were used. The Hco-pgp genes (1, 2, 3, 4, 9, 10, 11, 12, 14 and 16) were analysed for each isolate using the RT-qPCR technique. The Hco-pgp expressions were pairwise compared using the 2-ΔΔCt method and a t-test. The PARAISO isolate showed upregulation compared to the CENID-SAI isolate for Hco-pgp 1, 3, 9, 10 and 16 (P < 0.05), and the PARAISO isolate showed upregulation vs. FMVZ-UADY isolate for Hco-pgp 2 and 9 (P < 0.05), displaying 6.58- and 5.93-fold differences (P < 0.05), respectively. In contrast, similar Hco-pgp gene expression levels were recorded for FMVZ-UADY and CENID-SAI isolates except for Hco-pgp1 (P <0.1), which presented a significant upregulation (6.08-fold). The relative expression of Hco-pgp allowed confirming the IVM-resistant status of the PARAISO isolate and the IVM-susceptible status of the FMVZ-UADY isolate when compared to the CENID-SAI reference isolate. Therefore, understanding the association between the Hco-pgp genes expression of H. contortus and its IVM resistance status could help identifying the genes that could be used as molecular markers in the diagnosis of IVM resistance. However, it is important to consider the geographic origin of the nematode isolate and the deworming history at the farm of origin.

Curcumin-PVP improves the in vitro efficacy of ivermectin against resistant and susceptible Haemonchus contortus

Ivermectin (IVM) resistance in parasitic nematodes such as Haemonchus contortus has spurred a search for substances that help to recover its efficacy. One potential agent is the natural product curcumin (CUR). In this study, CUR was combined with polyvinylpyrrolidone (PVP) (CUR/PVP) to improve its solubility and biological applicability. This study determined the effect of CUR preincubation on the effective concentration 50% (EC50) of IVM in three H. contortus isolates with different susceptibilities to IVM. The IVM EC50 was determined for three H. contortus isolates with different IVM susceptibilities using the larval migration inhibition (LMI) test. The three isolates were (i) PARAISO (IVM resistant), (ii) FMVZ-UADY (IVM susceptible), and (iii) CENID-SAI INIFAP (reference IVM susceptible). The L3 of each isolate were preincubated for 3 h with one of three concentrations of CUR (μg curcumin/mL): CONC-1 (3.67), CONC-2 (5.67), or CONC-3 (8.48). Corresponding controls were performed without CUR. The EC50 of IVM was determined for each isolate after they were exposed to the different CUR concentrations. The EC50 of IVM differed between the isolates PARAISO > FMVZ-UADY > CENID-SAI INIFAP (P < 0.05). The CUR preincubation at CONC-1 did not decrease the EC50 of IVM for any of the three isolates, suggesting a hormetic effect. By contrast, CUR preincubation at CONC-2 or CONC-3 decreased the IVM EC50 for the PARAISO isolate (P < 0.05) compared with the reference isolate and reduced the EC50 of IVM for the FMVZ-UADY and CENID-SAI INIFAP isolates below the EC50 for the CENID-SAI INIFAP isolate without CUR preincubation. In conclusion, preincubation of H. contortus L3 with CUR reduced the EC50 of IVM for field isolates classified as resistant and susceptible to IVM. The CUR preincubation reduced the IVM resistance factor in the different isolates tested.

Understanding anthelmintic resistance in livestock using "omics" approaches

Widespread and improper use of various anthelmintics, genetic, and epidemiological factors has resulted in anthelmintic-resistant (AR) helminth populations in livestock. This is currently quite common globally in different livestock animals including sheep, goats, and cattle to gastrointestinal nematode (GIN) infections. Therefore, the mechanisms underlying AR in parasitic worm species have been the subject of ample research to tackle this challenge. Current and emerging technologies in the disciplines of genomics, transcriptomics, metabolomics, and proteomics in livestock species have advanced the understanding of the intricate molecular AR mechanisms in many major parasites. The technologies have improved the identification of possible biomarkers of resistant parasites, the ability to find actual causative genes, regulatory networks, and pathways of parasites governing the AR development including the dynamics of helminth infection and host-parasite infections. In this review, various "omics"-driven technologies including genome scan, candidate gene, quantitative trait loci, transcriptomic, proteomic, and metabolomic approaches have been described to understand AR of parasites of veterinary importance. Also, challenges and future prospects of these "omics" approaches are also discussed.

Effect of Ivermectin on the Expression of P-Glycoprotein in Third-Stage Larvae of Haemonchus contortus Isolated from China

Haemonchus contortus poses a severe hazard to the healthy development of the sheep industry and threatens the welfare of sheep. Ivermectin is the primary anthelmintic used for the prevention and treatment of H. contortus parasitism. However, the widespread and uncontrolled application of ivermectin has resulted in the development and spread of resistant strains of H. contortus. P-glycoprotein (P-gp) plays important roles in the pharmacology and toxicology of ivermectin, and changes in P-gp expression levels can be used to analyze the resistance of H. contortus to ivermectin. This study aimed to analyze the effects of ivermectin on P-gp expression in H. contortus L3 larvae isolated from China and to evaluate whether changes in P-gp expression levels can be used to analyze resistant H. contortus strains. In the absence of drug treatment, the ivermectin-resistant strains isolated in China showed increased expression of P-gp11 (p < 0.01) compared with sensitive strains from elsewhere, whereas the expressions of P-gp2 and P-gp9.1 were downregulated (p < 0.01). When the same strain was compared before and after drug treatment, obvious differences in expression were observed between the different strains. Ivermectin-induced P-gp expression was found to be very complex among the L3 larvae of different strains. In addition, it was confirmed that using P-gp to determine ivermectin resistance in H. contortus strains from different geographic environments can yield different results.

Importance of ABC Transporters in the Survival of Parasitic Nematodes and the Prospect for the Development of Novel Control Strategies

ABC transporters, a family of ATP-dependent transmembrane proteins, are responsible for the active transport of a wide range of molecules across cell membranes, including drugs, toxins, and nutrients. Nematodes possess a great diversity of ABC transporters; however, only P-glycoproteins have been well-characterized compared to other classes. The ABC transport proteins have been implicated in developing resistance to various classes of anthelmintic drugs in parasitic nematodes; their role in plant and human parasitic nematodes still needs further investigation. Therefore, ABC transport proteins offer a potential opportunity to develop nematode control strategies. Multidrug resistance inhibitors are becoming more attractive for controlling nematodes due to their potential to increase drug efficacy in two ways: (i) by limiting drug efflux from nematodes, thereby increasing the amount of drug that reaches its target site, and (ii) by reducing drug excretion by host animals, thereby enhancing drug bioavailability. This article reviews the role of ABC transporters in the survival of parasitic nematodes, including the genes involved, their regulation and physiological roles, as well as recent developments in their characterization. It also discusses the association of ABC transporters with anthelmintic resistance and the possibility of targeting them with next-generation inhibitors or nutraceuticals (e.g., polyphenols) to control parasitic infections.

Repertoire of P-glycoprotein drug transporters in the zoonotic nematode Toxocara canis

Toxocara canis has a complex lifecycle including larval stages in the somatic tissue of dogs that tolerate macrocyclic lactones. In this study, we investigated T. canis permeability glycoproteins (P-gps, ABCB1) with a putative role in drug tolerance. Motility experiments demonstrated that while ivermectin failed to abrogate larval movement, the combination of ivermectin and the P-gp inhibitor verapamil induced larval paralysis. Whole organism assays revealed functional P-gp activity in larvae which were capable of effluxing the P-gp substrate Hoechst 33342 (H33342). Further investigation of H33342 efflux demonstrated a unique rank order of potency for known mammalian P-gp inhibitors, suggesting that one or more of the T. canis transporters has nematode-specific pharmacological properties. Analysis of the T. canis draft genome resulted in the identification of 13 annotated P-gp genes, enabling revision of predicted gene names and identification of putative paralogs. Quantitative PCR was used to measure P-gp mRNA expression in adult worms, hatched larvae, and somatic larvae. At least 10 of the predicted genes were expressed in adults and hatched larvae, and at least 8 were expressed in somatic larvae. However, treatment of larvae with macrocyclic lactones failed to significantly increase P-gp expression as measured by qPCR. Further studies are needed to understand the role of individual P-gps with possible contributions to macrocyclic lactone tolerance in T. canis.

Pharmacokinetic and pharmacodynamic considerations for treating sarcoptic mange with cross-relevance to Australian wildlife

Sarcoptes scabiei is the microscopic burrowing mite responsible for sarcoptic mange, which is reported in approximately 150 mammalian species. In Australia, sarcoptic mange affects a number of native and introduced wildlife species, is particularly severe in bare-nosed wombats (Vombatus ursinus) and an emerging issue in koala and quenda. There are a variety of acaricides available for the treatment of sarcoptic mange which are generally effective in eliminating mites from humans and animals in captivity. In wild populations, effective treatment is challenging, and concerns exist regarding safety, efficacy and the potential emergence of acaricide resistance. There are risks where acaricides are used intensively or inadequately, which could adversely affect treatment success rates as well as animal welfare. While reviews on epidemiology, treatment strategies, and pathogenesis of sarcoptic mange in wildlife are available, there is currently no review evaluating the use of specific acaricides in the context of their pharmacokinetic and pharmacodynamic properties, and subsequent likelihood of emerging drug resistance, particularly for Australian wildlife. This review critically evaluates acaricides that have been utilised to treat sarcoptic mange in wildlife, including dosage forms and routes, pharmacokinetics, mode of action and efficacy. We also highlight the reports of resistance of S. scabiei to acaricides, including clinical and in vitro observations.

Transcriptomics of ivermectin response in Caenorhabditis elegans: Integrating abamectin quantitative trait loci and comparison to the Ivermectin-exposed DA1316 strain

Parasitic nematodes pose a significant threat to human and animal health, as well as cause economic losses in the agricultural sector. The use of anthelmintic drugs, such as Ivermectin (IVM), to control these parasites has led to widespread drug resistance. Identifying genetic markers of resistance in parasitic nematodes can be challenging, but the free-living nematode Caenorhabditis elegans provides a suitable model. In this study, we aimed to analyze the transcriptomes of adult C. elegans worms of the N2 strain exposed to the anthelmintic drug Ivermectin (IVM), and compare them to those of the resistant strain DA1316 and the recently identified Abamectin Quantitative Trait Loci (QTL) on chromosome V. We exposed pools of 300 adult N2 worms to IVM (10-7 and 10-8 M) for 4 hours at 20°C, extracted total RNA and sequenced it on the Illumina NovaSeq6000 platform. Differentially expressed genes (DEGs) were determined using an in-house pipeline. The DEGs were compared to genes from a previous microarray study on IVM-resistant C. elegans and Abamectin-QTL. Our results revealed 615 DEGs (183 up-regulated and 432 down-regulated genes) from diverse gene families in the N2 C. elegans strain. Of these DEGs, 31 overlapped with genes from IVM-exposed adult worms of the DA1316 strain. We identified 19 genes, including the folate transporter (folt-2) and the transmembrane transporter (T22F3.11), which exhibited an opposite expression in N2 and the DA1316 strain and were deemed potential candidates. Additionally, we compiled a list of potential candidates for further research including T-type calcium channel (cca-1), potassium chloride cotransporter (kcc-2), as well as other genes such as glutamate-gated channel (glc-1) that mapped to the Abamectin-QTL.

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.

Proteomic Comparison of Ivermectin Sensitive and Resistant Staphylococcus aureus Clinical Isolates Reveals Key Efflux Pumps as Possible Resistance Determinants

Ivermectin (IVM) is a versatile drug used against many microorganisms. Staphylococcus aureus is one of the most devastating microorganisms. IVM sensitive and resistant S. aureus strains were recently reported. However, the underlying molecular mechanisms of resistance are unknown. Clinical isolates of S. aureus were used for determination of the sensitivities against IVM by growth curve analysis and time-kill kinetics. Then, proteomic, and biochemical approaches were applied to investigate the possible mechanisms of resistance. Proteomic results showed a total of 1849 proteins in the dataset for both strains, 425 unique proteins in strain O9 (IVM sensitive), and 354 unique proteins in strain O20 (IVM resistant). Eight proteins with transport functions were differentially expressed in the IVM resistant strain. Among them, three efflux pumps (mepA, emrB, and swrC) were confirmed by qPCR. The IVM resistant S. aureus may overexpress these proteins as a key resistance determinant. Further experiments are required to confirm the exact mechanistic relationship. Nevertheless, the possibility of blocking these transporters to reverse or delay the onset of resistance and reduce selection pressure is potentially appealing.

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.

Ivermectin-induced gene expression changes in adult Parascaris univalens and Caenorhabditis elegans: a comparative approach to study anthelminthic metabolism and resistance in vitro

Background

The nematode Parascaris univalens is one of the most prevalent parasitic pathogens infecting horses but anthelmintic resistance undermines treatment approaches. The molecular mechanisms underlying drug activity and resistance remain poorly understood in this parasite since experimental in vitro models are lacking. The aim of this study was to evaluate the use of Caenorhabditis elegans as a model for P. univalens drug metabolism/resistance studies by a comparative gene expression approach after in vitro exposure to the anthelmintic drug ivermectin (IVM).

Methods

Twelve adult P. univalens worms in groups of three were exposed to ivermectin (IVM, 10^–13 M, 10^–11 M, 10^–9 M) or left unexposed for 24 h at 37 °C, and total RNA, extracted from the anterior end of the worms, was sequenced using Illumina NovaSeq. Differentially expressed genes (DEGs) involved in metabolism, transportation, or gene expression with annotated Caernorhabditis elegans orthologues were identified as candidate genes to be involved in IVM metabolism/resistance. Similarly, groups of 300 adult C. elegans worms were exposed to IVM (10^–9 M, 10^–8 M and 10^–7 M) or left unexposed for 4 h at 20 °C. Quantitative RT-PCR of RNA extracted from the C. elegans worm pools was used to compare against the expression of selected P. univalens candidate genes after drug treatment.

Results

After IVM exposure, 1085 DEGs were found in adult P. univalens worms but the relative gene expression changes were small and large variabilities were found between different worms. Fifteen of the DEGs were chosen for further characterization in C. elegans after comparative bioinformatics analyses. Candidate genes, including the putative drug target lgc-37, responded to IVM in P. univalens, but marginal to no responses were observed in C. elegans despite dose-dependent behavioral effects observed in C. elegans after IVM exposure. Thus, the overlap in IVM-induced gene expression in this small set of genes was minor in adult worms of the two nematode species.

Conclusion

This is the first time to our knowledge that a comparative gene expression approach has evaluated C. elegans as a model to understand IVM metabolism/resistance in P. univalens. Genes in P. univalens adults that responded to IVM treatment were identified. However, identifying conserved genes in P. univalens and C. elegans involved in IVM metabolism/resistance by comparing gene expression of candidate genes proved challenging. The approach appears promising but was limited by the number of genes studied (n = 15). Future studies comparing a larger number of genes between the two species may result in identification of additional candidate genes involved in drug metabolism and/or resistance.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05260-4.

Sensitivity of Haemonchus contortus to anthelmintics using different in vitro screening assays: a comparative study

BACKGROUND

Helminthiasis and resistance to commercial anthelmintic compounds are major causes of economic losses for livestock producers, resulting in an urgent need for new drugs and reliable in vitro screening tests capable of detecting potentially active products. Considering this, a series of novel benzimidazole derivatives (5-methylbenzimidazole 1,2-disubstituted, 5-carboxybenzimidazole, 5-methylbenzimidazole 2-one) was screened on exsheathed L3 (xL3) and on the adult stage of Haemonchus contortus (Kirby anthelmintic-susceptible McMaster isolate).

METHODS

This work presents the set-up of an automated motility assay on the xL3 stage of H. contortus using an infrared tracking device (WMicrotracker One) together with a larval development test (xL3 to L4) and a motility assay on the adult stage of H. contortus. A comparative study of the sensitivity of these in vitro assays using commercial anthelmintics with different mechanisms of action was carried out, also evaluating anthelmintic activity of a series of novel benzimidazole derivatives.

RESULTS

The automated xL3 assay had the great advantage of being able to analyze many compounds simultaneously, but it showed the limitation of having lower sensitivity, requiring higher concentrations of the commercial anthelmintics tested compared to those needed for the adult motility or development assays. Although none of the novel 1,2,5-tri-substituted benzimidazole derivatives could significantly decrease the motility of xL3s, one of them (1e) significantly affected the development of xL3s to L4, and five new compounds (1b, 1d, 1e, 2a and 2c) reduced the motility of H. contortus adult stage.

CONCLUSIONS

The analysis of the results strongly suggests that the in vitro xL3 to L4 development test, particularly for the L4 stage, could be closer to the pharmacological sensitivity of the adult stage of H. contortus (target of interest) for commercial anthelmintic selected, with different mechanisms of action, and for the series of benzimidazole derivatives assayed. Therefore, an automated motility assay on L4 using the infrared tracking device is being set up. Further studies will be conducted to evaluate the in vivo anthelmintic activity of the most active novel benzimidazole derivatives.

Transporter gene expression and Wolbachia quantification in adults of Dirofilaria immitis treated in vitro with ivermectin or moxidectin alone or in combination with doxycycline for 12 hours

Due to their marked larvicidal activity, macrocyclic lactones (MLs) are used for the prevention of heartworm disease ( Dirofilaria immitis) in dogs. They have also been shown to eliminate adult parasites after long-term administration, with a so-called "slow-kill" effect. In addition, recent studies have established that a combination of doxycycline, which eliminates the endosymbiont Wolbachia, and MLs has superior adulticide effects when compared to MLs alone. It has been hypothesized that the apparent synergism between doxycycline/MLs may be due to interaction with drug efflux transport proteins. The aim of the present study was to evaluate gene expression of several transport proteins in D. immitis adults treated in vitro either with doxycycline alone, ivermectin alone, moxidectin alone, or a combination of ivermectin or moxidectin with doxycycline for 12h. Quantitative PCR analysis showed a sex-dependent response to treatments. In female worms, Dim-pgp-10, Dim-haf-1 and Dim-haf-5 were upregulated compared to controls with doxycycline alone and when combined with ivermectin. Moxidectin did not induce any changes in gene expression. In males, moxidectin administered alone induced a slight increase in Dim-pgp-10, Dim-pgp-11and Di-avr-14, while ivermectin in combination with doxycycline produced significant upregulation of the ML receptor Di-avr-14. These results suggest possible synergism between the two drug classes and different susceptibility of males vs. females to adulticide effects.

Constitutive and differential expression of transport protein genes in Parascaris univalens larvae and adult tissues after in vitro exposure to anthelmintic drugs

The equine roundworm Parascaris univalens has developed resistance to the three anthelmintic substances most commonly used in horses. The mechanisms responsible for resistance are believed to be multi-genic, and transport proteins such as the P-glycoprotein (Pgp) family have been suggested to be involved in resistance in several parasites including P. univlaens. To facilitate further research into the mechanisms behind drug metabolism and resistance development in P. univalens we aimed to develop an in vitro model based on larvae. We developed a fast and easy protocol for hatching P. univalens larvae for in vitro studies, resulting in a hatching rate of 92 %. The expression of transport protein genes pgp-2, pgp-9, pgp-11.1, pgp-16.1 and major facilitator superfamily (MFS) genes PgR006_g137 and PgR015_g078 were studied in hatched larvae exposed to the anthelmintic drugs ivermecin (IVM) 10^-9 M, pyrantel citrate (PYR) 10^-6 M and thiabendazole (TBZ) 10^-5 M for 24 h. In comparison, the expression of these transport protein genes was studied in the anterior end and intestinal tissues of adult worms in vitro exposed to IVM, TBZ and PYR, at the same concentrations as larvae, for 3 h, 10 h and 24 h. Larval exposure to sub-lethal doses of IVM for 24 h did not affect the expression levels of any of the investigated genes, however larvae exposed to PYR and TBZ for 24 h showed significantly increased expression of pgp-9. In vitro drug exposure of adult worms did not result in any significant increases in expression of transport protein genes. Comparisons of constitutive expression between larvae and adult worm tissues showed that pgp-9, pgp-11.1, pgp-16.1 and MFS gene PgR015_g078 were expressed at lower levels in larvae than in adult tissues, while pgp-2 and MFS gene PgR006_g137 had similar expression levels in larvae and adult worms. All investigated transport protein genes were expressed at higher rates in the intestine than in the anterior end of adult worms, except pgp-11.1 where the expression was similar between the two tissues. This high constitutive expression in the intestine suggests that this is an important site for xenobiotic efflux in P. univalens. Despite the fact that the results of this study show differences in expression of transport protein genes between larvae and adult tissues, we believe that the larval assay system described here will be an important tool for further research into the molecular mechanisms behind anthelmintic resistance development and for other in vitro studies.

The diversity of ABC transporter genes across the Phylum Nematoda

It is estimated that one billion people globally are infected by parasitic nematodes, with children, pregnant women, and the elderly particularly susceptible to morbidity from infection. Control methods are limited to de-worming, which is hampered by rapid re-infection and the inevitable development of anthelmintic resistance. One family of proteins that has been implicated in nematode anthelmintic resistance are the ATP binding cassette (ABC) transporters. ABC transporters are characterized by a highly conserved ATP-binding domain and variable transmembrane regions. A growing number of studies have associated ABC transporters in anthelmintic resistance through a protective mechanism of drug efflux. Genetic deletion of P glycoprotein type ABC transporters in Caenorhabditis elegans demonstrated increased sensitivity to anthelmintics, while in the livestock parasite, Haemonchus contortus, anthelmintic use has been shown to increase the expression of ATP transporter genes. These studies as well as others, provide evidence for a potential role of ABC transporters in drug resistance in nematodes. In order to understand more about the family of ABC transporters, we used hidden Markov models to predict ABC transporter proteins from 108 species across the phylum Nematoda and use these data to analyze patterns of diversification and loss in diverse nematode species. We also examined temporal patterns of expression for the ABC transporter family within the filarial nematode Brugia malayi and identify cases of differential expression across diverse life-cycle stages. Taken together, our data provide a comprehensive overview of ABC transporters in diverse nematode species and identify examples of gene loss and diversification in nematodes based on lifestyle and taxonomy.

Sub-lethal doses of albendazole induce drug metabolizing enzymes and increase albendazole deactivation in Haemonchus contortus adults

The efficacy of anthelmintic therapy of farm animals rapidly decreases due to drug resistance development in helminths. In resistant isolates, the increased expression and activity of drug-metabolizing enzymes (DMEs), e.g. cytochromes P450 (CYPs), UDP-glycosyltransferases (UGTs) and P-glycoprotein transporters (P-gps), in comparison to sensitive isolates have been described. However, the mechanisms and circumstances of DMEs induction are not well known. Therefore, the present study was designed to find the changes in expression of CYPs, UGTs and P-gps in adult parasitic nematodes Haemonchus contortus exposed to sub-lethal doses of the benzimidazole anthelmintic drug albendazole (ABZ) and its active metabolite ABZ-sulfoxide (ABZSO). In addition, the effect of ABZ at sub-lethal doses on the ability to deactivate ABZ during consequent treatment was studied. The results showed that contact of H. contortus adults with sub-lethal doses of ABZ and ABZSO led to a significant induction of several DMEs, particularly cyp-2, cyp-3, cyp-6, cyp-7, cyp-8, UGT10B1, UGT24C1, UGT26A2, UGT365A1, UGT366C1, UGT368B2, UGT367A1, UGT371A1, UGT372A1 and pgp-3, pgp-9.1, pgp-9.2, pgp-10. This induction led to increased formation of ABZ metabolites (especially glycosides) and their increased export from the helminths' body into the medium. The present study demonstrates for the first time that contact of H. contortus with sub-lethal doses of ABZ (e.g. during underdose treatment) improves the ability of H. contortus adults to deactivate ABZ in consequent therapy.

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

Background

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

Methods

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

Results

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

Conclusions

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

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

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

Differential ABC transporter gene expression in adult Dirofilaria immitis males and females following in vitro treatment with ivermectin, doxycycline or a combination of both

Background

Combination doxycycline/macrocyclic lactone (ML) protocols have been shown to provide a more rapid adulticidal and microfilaricidal effect than either MLs or doxycycline alone, although female worms were reported to have a higher tolerance to treatments compared to male worms. The present study aimed to evaluate how ABC transporters may be involved in the synergic effect of the combination treatment. Adult worms of D. immitis were treated in vitro for 24 hours with doxycycline (DOXY), ivermectin (IVM) and a combination of both, and changes in the modulation of ABC transporter genes were measured. Levels of doxycycline inside different treatment media, post-treatment, were determined through HPLC analysis.

Results

Quantitative RT-PCR analysis showed the presence of changes in the modulation of ABC transporter genes evaluated in this study. In particular, in female worms, the combination treatment induced a substantial increase in gene expressions, especially of Dim-pgp-10 and Dim-haf-4; whereas in male worms, the greatest increase in gene expression was observed for Dim-pgp-10 and Dim-pgp-11 when treated with DMSO + IVM and DMSO + DOXY/IVM. HPLC analysis of the DOXY concentrations in the media after in vitro treatments of male worms showed a slight difference between the DMSO + DOXY samples and the combination (DMSO + DOXY + IVM), while no difference was observed among females.

Conclusions

Further studies are required to explain whether the modulation of cellular efflux plays a role, even partially, in the adulticide effect of doxycycline/macrocyclic lactone combinations in heartworm-infected dogs. To the authors’ knowledge, this is the first study to evaluate P-gp expression in adult D. immitis.

Y75B8A.8 (HC8) protein of Haemonchus contortus: A functional inhibitor of host IL-2

Interleukin 2 (IL-2) is an important immune regulatory factor in the immune response of the host. However, little is known about the inhibitor of host IL-2 in Haemonchus contortus infection. In this study, we found that globin domain-containing protein (HCGB) and Protein Y75B8A.8 (HC8) from H contortus excretory and secretory products are two binding proteins of IL-2 in goats. The yeast two-hybrid screening further validated the positive interactions of IL-2 with HCGB and HC8. Meanwhile, we found that HC8 had inhibitory effects on IL-2-induced peripheral blood mononuclear cell (PBMC) proliferation, while HCGB did not. Furthermore, transcriptional analysis revealed that HC8 could block the IL-2-activated signalling pathway. Our results showed that HC8 was a functional inhibitor of goat IL-2.

Minimal modulation of macrocyclic lactone susceptibility in Caenorhabditis elegans following inhibition of cytochrome P450 monooxygenase activity

Anthelmintic and in particular macrocyclic lactone (ML) resistance is a widespread problem in trichostrongyloid parasitic nematodes, yet mechanisms of ML resistance are still poorly understood. In the absence of target-site changes in resistant parasite field populations, increased drug extrusion and xenobiotic metabolism have been implicated in modification of susceptibility to MLs. In addition to P-glycoproteins, cytochrome P450 monooxygenases (CYPs) were considered to be involved in ML resistance. CYPs are highly divergent in nematodes with about 80 genes in the model organism Caenorhabditis elegans. Using larval development assays in the C. elegans model, piperonyl butoxide (PBO) and a temperature-sensitive variant of the emb-8 cytochrome reductase were used for chemical and genetic ablation of CYP activity. Additionally, a loss-of-function variant of cyp-14A5 was characterized to determine whether increased expression of this CYP in an ivermectin (IVM)-tolerant C. elegans line might be related to the phenotype. In a preliminary experiment with PBO, susceptibility to 5 nM IVM was synergistically increased by PBO. However, effects of genetic ablation of CYP activity on the EC₅₀ values were small (1.5-fold decrease) for IVM and not significant for moxidectin (MOX). However, due to the steep concentration-response-curves, there were again strong differences between the wild-type and the CYP deficient genotype at individual IVM but not MOX concentrations. Although these results suggest small but significant effects on the susceptibility level of C. elegans to IVM, the cyp14A5 gene proposed by a previous study as candidate was ruled out since it was neither IVM/MOX inducible nor did a strain with a loss-of-function allele show increased susceptibility to either drug. In conclusion, the effect of the CYP system on IVM susceptibility in C. elegans is at best low while effects on MOX susceptibility were not detected. The previously suggested candidate cyp14A5 could be excluded to be involved in ML metabolism.

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.

Detection and quantification of Parascaris P-glycoprotein drug transporter expression with a novel mRNA hybridization technique

Macrocyclic lactone-resistant Parascaris have been reported throughout the world. In part, the drug resistant phenotype is hypothesized to be associated with ATP-binding cassette transporters known as P-glycoproteins. In many systems, P-glycoproteins efflux drugs out of cells thereby precluding drug binding to target receptors. Parascaris may evade macrocyclic lactone-mediated death by effluxing drugs away from target receptors in the nervous system. Alternatively, P-glycoprotein expression in the gut or body wall could prevent penetration of drugs into the body of the parasite altogether. In the present study, we evaluate expression of Peq-pgp-11 and Peq-pgp-16 using a novel multiple nucleic acid hybridization method. This method allowed for visualization of individual mRNA transcripts within fixed tissue sections of Parascaris adults. Our investigation revealed expression of Peq-pgp-11 and Peq-pgp-16 in the intestine, body wall, nerves, lateral cords, and reproductive tissues of male and female parasites. These results suggest that Pgp could efflux drugs locally at the level of parasite neuronal tissue as well as at sites of entry for drugs such as the hypodermis and intestine. The multiple nucleic acid hybridization method could be useful for providing tissue context for gene expression in a variety of nematode parasites.

Differences in protein expression associated with ivermectin resistance in Caenorhabditis elegans

The indiscriminate administration of synthetic anthelmintics such as ivermectin contributes to the selection of subpopulations capable of resisting the drugs' effects. To understand the mechanisms of ivermectin resistance in Caenorhabditis elegans, this study attempted to identify molecular targets. C. elegans lineages that were sensitive and resistant to ivermectin were used. Collected nematodes were added to an extraction buffer and macerated in liquid nitrogen for protein extraction. The extracted proteins were separated according to molecular weight by SDS-PAGE to verify their integrity. Subsequently, proteins from both lineages were separated using two-dimensional electrophoresis. The gels were analyzed and the relevant spots were excised and identified by mass spectrometry (NanoESI-Q-TOF and MASCOT®) and subsequently assessed by GO enrichment and STRING® analyses. The increased expression of proteins associated with high metabolic activity, such as ATP-2 and ENOL-1, which are responsible for ATP synthesis, was observed. Furthermore, proteins with involvement in mediating muscular function (MLC-1, ACT-1, and PDI-2), signaling (FAR-1 and FAR-2), and embryo development (VHA-2) were identified. Protein interaction analysis indicated that the majority of the identified proteins in the resistant lineages participated in the same reaction triggered by ivermectin.

In vitro activity of ivermectin against Staphylococcus aureus clinical isolates

Background

Ivermectin is an endectocide against many parasites. Though being a macrocyclic lactone, its activity against bacteria has been less known, possibly due to the fact that micromolar concentrations at tissue levels are required to achieve a therapeutic effect. Among pathogenic bacteria of major medical significance, Staphylococcus aureus cause a number of diseases in a wide variety of hosts including humans and animals. It has been attributed as one of the most pathogenic organisms. The emergence of methicillin resistance has made the treatment of S. aureus even more difficult as it is now resistant to most of the available antibiotics. Thus, search for alternate anti-staphylococcal agents requires immediate attention.

Methods

Twenty-one clinical isolates of S. aureus were isolated from bovine milk collected from Lahore and Faisalabad Pakistan. Different anthelmintics including levamisole, albendazole and ivermectin were tested against S. aureus to determine their minimum inhibitory concentrations. This was followed-up by growth curve analysis, spot assay and time-kill kinetics.

Results

The results showed that ivermectin but not levamisole or albendazole exhibited a potent anti-staphylococcal activity at the concentrations of 6.25 and 12.5 μg/ml against two isolates. Interestingly, one of the isolate was sensitive while the other was resistant to methicillin/cefoxitin.

Conclusions

Our novel findings indicate that ivermectin has an anti-bacterial effect against certain S. aureus isolates. However, to comprehend why ivermectin did not inhibit the growth of all Staphylococci needs further investigation. Nevertheless, we have extended the broad range of known pharmacological effects of ivermectin. As pharmacology and toxicology of ivermectin are well known, its further development as an anti-staphylococcal agent is potentially appealing.

P-glycoprotein-9 and macrocyclic lactone resistance status in selected strains of the ovine gastrointestinal nematode, Teladorsagia circumcincta

The Teladorsagia circumcincta P-glycoprotein-9 (Tci-pgp-9) gene has previously been implicated in multiple-anthelmintic resistance in this parasite. Here we further characterise genetic diversity in Tci-pgp-9 and its possible role in ivermectin (IVM) and multi-drug resistance using two UK field isolates of T. circumcincta, one susceptible to anthelmintics (MTci2) and the other resistant to most available anthelmintics including IVM (MTci5). A comparison of full-length Tci-pgp-9 cDNA transcripts from the MTci2 and MTci5 isolates (∼3.8 kb in both cases) indicated that they shared 95.6% and 99.5% identity at the nucleotide and amino acid levels, respectively. Nine non-synonymous SNPs were found in the MTci5 sequences relative to their MTci2 counterparts. Twelve genomic sequence variants of the first internucleotide binding domain of Tci-pgp-9 were identified and up to 10 of these were present in some individual worms, strongly supporting previous evidence that amplification of this gene has occurred in T. circumcincta. On average, fewer distinct sequence variants of Tci-pgp-9 were present in individual worms of the MTci5 isolate than in those of the MTci2 isolate. A further reduction in the number of sequence variants was observed in individuals derived from an IVM-treated sub-population of MTci5. These findings suggest that Tci-pgp-9 was under purifying selection in the face of IVM treatment in T. circumcincta, with some sequence variants being selected against.

•

Tci-pgp-9 cDNA sequences from the MTci2 and MTci5 isolates shared 95.6% identity.

•

Comparison of Tci-pgp-9 sequences from MTci5 and MTci2 showed 9 non-synonymous SNPs.

•

Multiple Tci-pgp-9-IBDA sequence variants occurred in the majority of individual larvae.

•

Tci-pgp-9 appears to be under purifying selection after anthelmintics are applied.

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.

Increased expression of ATP binding cassette transporter genes following exposure of Haemonchus contortus larvae to a high concentration of monepantel in vitro

Background

There is some evidence that ATP binding cassette (ABC) transporters play a role in resistance to anthelmintics, particularly against macrocyclic lactones. Some anthelmintics, including ivermectin (IVM), have been shown to induce transcription of multiple ABC transporters in nematodes; however, the effects of monepantel (MPL) on transcription of these transporter genes has not been studied.

Methods

Larvae of two MPL-susceptible isolates of Haemonchus contortus were exposed to MPL at two concentrations (2.5 and 250 μg/ml) for periods of 3, 6 and 24 h. Transcription levels of sixteen ABC transporter genes were measured at the end of the incubation periods. The consequences of MPL exposure were examined by measuring rhodamine-123 efflux from the larvae, and their sensitivity to subsequent treatment with IVM or levamisole.

Results

Multiple ABC transporter genes showed significantly higher transcription in both worm isolates following exposure to MPL at 250 μg/ml for 3, 6 or 24 h, particularly the P-glycoprotein (P-gp) genes pgp-11, pgp-12 and pgp-14. Of these, only pgp-11 maintained the elevated levels 24 h after the end of the drug exposure period. In contrast, there was only a single instance of low-level upregulation as a result of exposure to MPL at 2.5 μg/ml. Larvae exposed to MPL at 250 μg/ml showed an increased efflux of rhodamine-123 and a proportion of the larval population showed an ability to subsequently tolerate higher concentrations of IVM in migration assays. There was no increased tolerance to IVM following pre-exposure to MPL at 2.5 μg/ml.

Conclusions

Exposure of H. contortus larvae to 250 μg/ml MPL results in increased transcription of multiple transporter genes and increased R-123 efflux. The subsequent ability of a proportion of the larvae to tolerate IVM suggests a protective role of ABC transporters across different chemical entities. However, these observations were only made at a concentration of MPL well above that experienced by parasitic life stages in vivo, and hence their significance remains unclear.

Electronic supplementary material

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