Increased expression of ATP binding cassette transporter genes following exposure of Haemonchus contortus larvae to a high concentration of monepantel in vitro
© The Author(s). 2016
Received: 4 July 2016
Accepted: 19 September 2016
Published: 29 September 2016
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.
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.
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.
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.
KeywordsMonepantel ABC transporters Gene induction Haemonchus contortus Ivermectin tolerance
Monepantel (MPL) (trade name Zolvix®) was the first new anthelmintic for livestock use for over 25 years when it appeared in 2009. It was first introduced in New Zealand for the control of gastrointestinal nematodes (GINs) and associated diseases in small ruminants, and later in Australia and the United Kingdom in 2010. It shows broad-spectrum anthelmintic activity, targeting larval and adult stages of the most important species of sheep GINs [1, 2]. The compound is an amino-acetonitrile derivative (AAD) that acts as a positive allosteric modulator of a nematode-specific clade of nicotinic acetylcholine receptor (nAChR) subunits. Genetic screening of Caenorhabditis elegans revealed nAChR subunit ACR-23 as a target for MPL action , while Rufener et al.  suggested that Hco-MPTL-1 and other nAChR subunits of the DEG-3 subfamily are the targets for MPL in H. contortus.
Resistance to MPL was first reported by Scott et al.  in New Zealand. The authors reported that MPL showed no efficacy in terms of egg count reduction in goats or reduction of worm burden in sheep, and was ineffective against at least two GINs (Teladorsagia circumcincta and Trichostrongylus colubriformis). There is also evidence of MPL resistance in Teladorsagia spp. and T. colubriformis on a goat farm in New South Wales, Australia . Resistance to MPL has also been reported in H. contortus in Uruguay, where the authors reported poor efficacy against this nematode with 42 and 82 % faecal egg count reduction at two different sheep farms . A more recent study has also reported H. contortus population resistant to MPL at a sheep farm in Netherlands . Early work on MPL resistance in laboratory-selected isolates of H. contortus showed that the resistance was associated with a number of changes to the target site. These mutations resulted in the generation of various truncated versions of the target receptor that would be expected to be non-functional . However, there is no information at present on whether such target site changes are involved in any of the cases of resistance reported from the field. We were interested in whether ABC transporters may interact with MPL in nematodes, as a number of anthelmintics are known to be substrates of these transporters, and increases in efflux activity are associated with anthelmintic resistance in some cases, particularly for macrocyclic lactones [8–10]. As AAD compounds are hydrophobic, they may also be substrates for the drug transport proteins .
There is a considerable evidence that exposure to anthelmintics results in increased transcription of some ABC transporter genes, both in vitro [11–14] and in vivo [15–17]. These studies have generally concluded that such a transcriptional response suggests a role for the transporters in the efflux of the transcription-inducing anthelmintic, that is, the anthelmintic is likely a substrate for the transporters. Therefore, as a first step in understanding the interaction of efflux pathways with MPL in nematodes, the present study aimed to explore the effects of MPL exposure on the transcription patterns of ABC transporters in third-stage larvae (L3) of two MPL-susceptible isolates of H. contortus. We further investigated the phenotypic effects of this drug exposure by (i) measuring the efflux of the fluorescent dye rhodamine-123 (R-123) from the L3 pre-exposed to MPL compared to controls, and (ii) examining whether MPL pre-exposure enabled the L3 to tolerate subsequent exposure to ivermectin (IVM) and levamisole (LEV) in larval migration assays.
Two isolates of H. contortus were used in the present study: (i) Kirby, a drug susceptible isolate recovered from the field at the University of New England, Kirby Research Farm in 1986 ; and (ii) Wallangra (WAL), a multi-drug resistant isolate collected from the New England region of New South Wales (NSW) in 2003; at the time of isolation, it was resistant to benzimidazole, closantel, LEV, IVM and moxidectin . The isolate has been further selected over at least five generations with moxidectin (Cydectin®) and is now unaffected by the recommended dose of moxidectin.
Both isolates were recovered from the field before the introduction of MPL, and show equivalent sensitivity to this drug in larval development assays . Infected animals were housed at the Commonwealth Scientific and Industrial Research Organisation (CSIRO) Agriculture and Food FD McMaster laboratory at Armidale, NSW. Faeces were collected from infected animals and sent to the CSIRO Agriculture laboratories at the Queensland Bioscience Precinct, Brisbane, QLD, in zip-lock bags. Third-stage larvae (L3) were harvested from faecal cultures held at 27 °C for 7 days. The L3 were stored at 15 °C, and used for experiments within 4 weeks.
In vitro monepantel exposure
The commercially available drench product Zolvix® (Novartis Animal Health, Australia) (25 mg MPL/ml) was used as a source of MPL in this study. Multiple separate anthelmintic solutions were prepared by two-fold serial dilutions in dimethyl-sulfoxide (DMSO).
250 μg/ml: the highest concentration that L3 could be exposed to while maintaining an ability to migrate at equivalent levels to controls in subsequent migration assays, that is, the highest concentration that could be used without compromising the fitness of the larvae as measured by migration assays (migration at 250 μg/ml approximately 90 % of controls) (unpublished data); and
2.5 μg/ml: as an approximation of the range of MPL concentrations of 2–4 μg/g measured by Lifschitz et al.  in the abomasal contents of sheep within the first 48 h following administration of Zolvix®.
For use in RNA extraction, larvae were centrifuged (3000× g, 1 min), and the pellets were washed with 1 ml water, centrifuged again (3000× g, 1 min), and the pelleted L3 were stored at -80 °C.
For use in R-123 and migration assays, larvae were centrifuged (3000× g, 1 min), and the pellets were washed with 1 ml water, centrifuged again (3000× g, 1 min), and the L3 were resuspended in water and used immediately.
RNA extraction, cDNA synthesis and qPCR were performed as described by Raza et al. . Briefly, total RNA was extracted from 30,000 L3 using the RNeasy mini kit (Qiagen, Hilden, Germany) as recommended by the manufacturer. Turbo-DNase (Ambion, Carlsbad, USA) was used to remove genomic DNA. cDNA was synthesized using DNase-treated RNA with superscript III reverse transcriptase (Invitrogen, Carlsbad, USA) according to the manufacturer’s instructions. The cDNA samples were diluted to a concentration of 4 ng/ml for downstream applications.
The primers used in this study for quantitative PCR of 11 pgp genes, two multidrug-resistance protein genes (mrp-1 and mrp-5), two genes from the ABCF family(abcf-1 and abcf-2) and one half transporter gene (haf-6) were as reported by Raza et al. . Three housekeeping genes (GAPDH, actin and β-tubulin) were used as reference genes for the qPCR analysis (see Additional file 1: Table S1). The SYBR Green dye system (Applied Biosystems, Warrington, UK) was used in a Vii A7 thermocycler (Applied Biosystems, USA) under the following PCR cycling conditions: 50 °C for 2 min, 95 °C for 10 min (stage I), followed by 40 cycles of 95 °C for 15 s, 60 °C for 1 min (stage II) and a melt curve stage at 95 °C for 15 s, 60 °C for 1 min and 95 °C for 15 s (stage III). Three separate extractions for each treatment were examined, with each PCR run in quadruplicate. Reaction efficiencies, determined by standard curves, were in the range between 80 and 99 %. The homogeneity of the PCR products was ensured by (i) analysing the melt-curves in each run, (ii) visualizing a single band for each PCR product on 2 % electrophoresis gels, and (iii) cloning of the selected primer products into the vector PCR2.1 (Invitrogen, USA) followed by sequencing (Big-Dye terminator, V3.1; Applied Biosystems, USA) using M13 forward/reverse primers. Expression values for all genes in each sample were normalised to the housekeeping genes using REST 2009 (version v2.0.13) to determine the effects of MPL pre-exposure on the transcription of ABC transporters (drug exposure regimen described above) using DMSO treated samples as control for each isolate. The triplicate expression values were log10 transformed and analysed using repeated measures ANOVA with Fischer’s Least Significant Difference (LSD) post-hoc test (P < 0.05) to compare the transcription profiles of control and MPL-exposed samples in GraphPad Prism(version 6.01).
Rhodamine (R-123) efflux assay
R-123 efflux was measured in order to determine the effects of MPL pre-exposure (2.5 μg/ml and 250 μg/ml) on the efflux activity of ABC transporters in Kirby and WAL L3stage larvae, using the protocols described earlier . Briefly, approximately 20,000 L3 were exposed to MPL or DMSO for 3 and 6 h (as described above) and then placed into a solution of R-123 (2 ml of 1.5 μM) for 15 min in the dark at room temperature on a roller-mixer. The worms were centrifuged (3000× g, 1 min) and washed with 1 ml of distilled water. The pellet was resuspended in distilled water (1 ml) and placed on a roller for 60 min in the dark. Finally, the worms were centrifuged (3000× g, 1 min) and the supernatant was collected and stored in dark for 60 min before use. R-123 concentration in the supernatant was detected by measuring the specific fluorescence (λ for excitation = 495 nm and λ for emission = 525 nm) using a spectrophotometer (Spectra Max M3, Molecular Devices®, Sunnyvale, CA, USA). The concentration of R-123 in each experimental sample was calculated from a standard curve, determined using a range of R-123 concentrations. Each experiment consisted of duplicate incubations for each treatment. Three separate experiments were performed and data (n = 3 for each treatment) were log10-transformed and analysed using repeated measures ANOVA followed by Fischer’s LSD test (P < 0.05) (GraphPad Prism, version 6.01).
Larval migration assay (LMA)
LMAs were used to measure the changes in tolerance of Kirby and WAL L3 to IVM and LEV following pre-exposure to MPL for 3 or 6 h. We were unable to assess whether MPL pre-exposure resulted in tolerance to MPL itself as this drug does not inhibit larval migration, even at high concentrations, and is therefore not suitable for use in LMAs . The assay was performed following the procedure as reported by Raza et al. . This method was based on that described earlier by Kotze et al. , except with the drug exposure and migration periods reduced from 24–48 h to 30 min each in order to allow for assessment of migration inhibition without providing any time for drug-induced gene expression changes to occur during the time course of the migration assay itself.
Briefly, MPL pre-treated (2.5 μg/ml and 250 μg/ml) or control (pre-exposure to 1 % DMSO only) larvae (3800 L3/ ml of water) were exposed to a range of IVM and LEV concentrations in 96-well microtiter plates for 30 min at 27 °C. Final concentrations for IVM and LEV ranged from 100-0.195 μg/ml and 12.5–0.024 μg/ml, respectively. After 30 min, the larvae were collected using a multi-tip pipette and placed into migration plates with 20 μm filters above receiver plate wells (Millipore, Bayswater, Australia). The larvae were allowed to migrate through the filters for a period of 30 min, before the filter plates were removed and the larvae in the receiver wells killed using Lugol’s iodine, and counted.
Migration assays consisted of triplicate assay wells at each IVM or LEV concentration. Three separate experiments were performed for each isolate. Data were analysed using non-linear regression in GraphPad Prism (version 6.01). IC50 values and 95 % confidence intervals were calculated based on the pooled data from each set of nine assays, and significant differences were determined by the overlap of 95 % confidence intervals.
The dose response to IVM in the short-term migration assays used in the present study showed the presence of a plateau in the response at the highest drug concentrations, as described previously by Raza et al. . The percent migration remained at a constant level (above 0 %) over the highest 2 to 3 drug concentrations. Hence, for the analysis of the IVM dose-response data, we used a non-normalised model in GraphPad (‘top to bottom’), with a variable slope. The output of this analysis provided us with two parameters with which to compare populations: first, the percent migration at the dose-response plateau (that is, the percentage of the population unaffected by the highest concentrations of IVM in the assay), and secondly, the IC50 in the remaining proportion of the population that had shown a dose-response to the drug (for instance, if the plateau existed at a level of 30 % migration, then the IC50 value defined the response to the drug in the remaining 70 % of the worm population alone). In addition, to further determine the significant differences between the percentage of larval migration following exposure to DMSO and MPL, the percent migration in MPL-pre-exposed compared to DMSO-pre-exposed worms at each IVM concentration were compared using a t-test with Welch’s correction (P < 0.05) (GraphPad Prism).
Transcriptional response of ABC transporters to MPL exposure in Kirby larvae
Relative transcription levels of ABC transporter genes in the 250 μg/ml MPL-treated Kirby and WAL larvae compared to DMSO-treated controls. Significant (P < 0.05) increases or decreases in transcription in the MPL-treated larvae compared to the DMSO-treated larvae are indicated by * and #, respectively
Following exposure to MPL250 μg/ml for 3 and 6 h, pgp-12 transcription was increased by 5.4-fold (P = 0.007) and 9.1-fold (P = 0.002), respectively, while pgp-14 showed 8.7-fold (P = 0.002) and 4.6-fold (P = 0.03) increases, respectively (Fig. 1d, e) (Table 1). Exposure to MPL250 μg/ml for 3 h resulted in significant downregulation of pgp-2 and haf-6 by 3.4-fold (P = 0.002) and 1.7-fold (P = 0.02), respectively. These decreases in transcription were short-lived, as, by 6 h, transcription levels for these two genes had returned to levels equivalent to DMSO-controls (Fig. 1d, e). Four P-gp genes and mrp-1 were significantly upregulated (P < 0.0001) following 24 h exposure to MPL250 μg/ml; the increases were as follows: 8.8-fold for pgp-11, 2.4-fold for pgp-14, 2.3-fold for pgp-12, 2.2-fold for pgp-2 and 1.8-fold for mrp-1. In addition, there were some instances of downregulation, with pgp-10 (1.6-fold, P < 0.0001), pgp-16 (2-fold, P < 0.0001) and mrp-5 (1.4-fold, P = 0.009) showing significantly decreased transcriptions relative to DMSO controls (Fig. 1f) (Table 1). The most consistent changes in larval gene expression following exposure to MPL250 μg/ml were the increased transcription of pgp-12 and pgp-14 at each of the three time points. There was a single instance of inconsistent changes in gene expression values between the two MPL concentrations: after 24 h exposure, pgp-11 showed an increased expression level in response to MPL250 μg/ml (8.8-fold), compared to a decrease in expression in response toMPL2.5 μg/ml (1.6-fold) (compare Fig. 1c, f).
Transcriptional response of ABC transporters to MPL exposure in WAL larvae
In contrast, exposure of WAL larvae to MPL250 μg/ml for 3 h resulted in significantly increased transcription of five P-gp genes and the haf-6 (Fig. 2d) (Table 1). Of note here was the greater number of upregulated genes in WAL compared to Kirby isolate following 3 h exposure to MPL (compare Figs. 1d and 2d). The magnitudes of the increases for WAL were as follows: 7.9-fold for pgp-2 (P = 0.002), 8.4-fold for pgp-9.1 (P = 0.002), 11.7-fold for pgp-11 (P = 0.0004), 9.4-fold for pgp-12 (P = 0.001), 12.2-fold for pgp-14 (P = 0.0002) and 20.5-fold for haf-6 (P = 0.0007). There were no instances of significant downregulation at the 3 h and 6 h time points. The increased transcription of pgp-2, pgp-9.1 and haf-6 genes observed at 3 h was short lived as it had returned to DMSO control levels by the 6 h time point for these three genes (Fig. 2d, e) (Table 1). In contrast, the upregulation of pgp-11, pgp-12 and pgp-14 was maintained at the 6 h time point; pgp-11 7.4-fold (P = 0.0007), pgp-12 8.2-fold (P = 0.001), and pgp-14 5.5-fold (P = 0.002). These increases also occurred at the 24 h time point; pgp-11 6.9-fold (P < 0.0001), pgp-12 2.9-fold (P < 0.0001) and pgp-14 2.6-fold (P < 0.0001). For each of these genes, the magnitude of fold change in expression relative to controls decreased over the sequential time points from 3 to 24 h (for example from 12.2-fold to 2.6-fold for pgp-14). Following 24 h exposure to MPL, transcription of pgp-9.3, pgp-10, pgp-16 and mrp-5 were significantly downregulated at levels of 1.4-fold (P = 0.003), 1.7-fold (P < 0.0001), 1.3-fold (P = 0.02) and 1.3-fold (P = 0.03), respectively (Fig. 2f).
There were several instances of inconsistent changes in gene expression values between the two MPL concentrations, with expression reduced at the lower MPL concentration and increased at the higher concentration for haf-6 at 3 h, as well as pgp-2 and pgp-11 at 24 h. As described above, this difference in response to the two MPL concentrations was also observed for Kirby with pgp-11 at 24 h.
Overall, for WAL larvae, the most marked effects of MPL exposure were the significant upregulation of a number of transporter genes at the 3 h time point for MPL250 μg/ml, and the sustained nature of these increases for pgp-11, pgp-12 and pgp-14, alongside the absence of any upregulation by MPL2.5 μg/ml. This pattern was very similar to that seen with Kirby larvae, with the exception of the unchanged expression for pgp-11 at 3 h for the Kirby larvae (Table 1).
Over this period, the fold increase in this gene in MPL-treated larvae compared to controls decreased from 11.7-fold (P = 0.0004) to 3-fold (P = 0.004) (compare Figs. 2d and 3). There was some variability in the expression patterns of pgp-9.1 (as indicated by the large SEM bars), however, statistical analysis across the three replicate experiments showed that the expression level of this gene was not significantly different to controls.
Phenotypic characterization of the MPL-treated larvae
Rhodamine-123 efflux assay
The levels of R-123 efflux were equivalent in Kirby and WAL L3 in the absence of any drug exposure (DMSO-treated) at both 3 and 6 h (P = 0.24 and P = 0.12) (data not shown).
Larval migration assay
The functional consequences of changes in transcription of transporter genes following exposure to the two concentrations of MPL were further examined using migration assays. Following 3 h exposure to MPL2.5 μg/ml, there were no significant changes to the tolerance of both Kirby and WAL L3to IVM and LEV, both in terms of the IVM IC50 values and the percent migration at the dose-response plateau at the high IVM concentrations (data not shown). Given, this lack of effect of 3 h exposure to MPL2.5 μg/ml, alongside the absence of gene expression changes in 6 h MPL2.5 μg/ml-treated Kirby and WAL larvae (Figs. 1b, 2b), we did not measure the response of 6 h MPL2.5 μg/ml pre-exposed larvae to IVM or LEV.
Response of Haemonchus contortus larvae to ivermectin following pre-exposure to monepantel
Drug conc.a (μg/ml)
95 % CI
95 % CI
In contrast to these changes in response to the highest IVM concentrations, the response of the remaining worm population did not change significantly, as shown by equivalent IC50 values, and a lack of statistical differences in percent migration at individual IVM concentrations (t-test analysis) (Table 2, Fig. 5). That is, the IVM tolerance was only observed in a component of the worm population. The only exception to this, in terms of the t-test analysis, was the decreased migration for MPL treated larvae at an IVM concentration of 0.2 μg/ml for WAL larvae (Fig. 5c). There were no significant changes in LEV IC50 values following exposure of Kirby and WAL larvae to MPL250 μg/ml (see Additional file 2: Figure S1; Additional file 3: Table S2).
In this study, we have demonstrated that exposure of H. contortus L3 to a high concentration of MPL (250 μg/ml) results in increased transcription of a number of ABC-transporter genes. Furthermore, exposure to MPL at this concentration results in greater efflux of R-123 from the L3, as well as increased tolerance to IVM in a proportion of the larval population. In contrast, pre-exposure to a 100-fold lower concentration of the compound (2.5 μg/ml) only resulted in a single instance of low-level upregulation of a transporter gene (pgp-11), an increase in R-123 efflux in only one out of four isolate/ time point instances, and did not result in any IVM tolerance in subsequent migration assays. Hence, the effects of MPL on gene transcription, R-123 efflux and subsequent IVM tolerance were dependent on the concentration of the compound to which the larvae were exposed. This study is the first, to our knowledge, to indicate an interaction between MPL and ABC transporters in nematodes. However, our data describe the effects of MPL on the regulation of ABC transporter gene expression rather than a direct effect of the drug on the transporter proteins themselves. The demonstration of functional consequences of the gene transcription increases, in terms of R-123 efflux and IVM tolerance, indicates that the transcription increases most likely result in increased transporter protein activity towards R-123 and IVM. Hence, the present study provides evidence of an interaction of MPL with the regulatory mechanism for ABC transporter gene expression, and subsequent protein synthesis, rather than evidence as to whether MPL itself is a substrate for the transporters.
As described above, the transporter gene upregulation/R-123 efflux / IVM tolerance relationships were only observed at the high MPL concentration. This concentration was well above the levels that an adult worm would encounter in the abomasum of a sheep, as measured by Lifschitz et al.  to be 2–4 μg/g of abomasal content 48 h after administration of the drug. On the other hand, the lower concentration of MPL used in the present study, which was selected based on it approximating this in vivo concentration range, resulted in only one instance of transporter gene upregulation, and did not result in any IVM tolerance. It is clear that while the present study indicates an interaction between MPL at 250 μg/ml and the regulatory mechanism for ABC transporter transcription in H. contortus, it does not provide evidence for such a role in vivo. Further experiments would be needed to determine if MPL interacts with transporters in adult worms in vivo.
While the patterns of upregulation in Kirby and WAL L3s were quite similar at 6 and 24 h time points, the WAL larvae showed a much greater level of upregulation (in terms of fold increases as well as the number of transporters affected) at the 3 h time point compared to Kirby (compare Figs. 2d and 1d). Our earlier study, reported the upregulation of ABC transporters in the WAL isolate after exposure to IVM and LEV, alongside no changes in Kirby at the 3 and 6 h time points examined in that study . This suggested that the rapid upregulation in response to these two drugs in the WAL isolate may be a component of the resistance shown by this isolate to the two anthelmintics. WAL and Kirby are both susceptible to MPL , having been isolated from the field before the use of this drug. The present study therefore indicates that this increased responsiveness of WAL larvae also occurs with respect to its response to drugs to which it is not resistant, and to which it has never been exposed. Thus, if the rapid upregulation following exposure to IVM and LEV is a component of the resistance shown towards these drugs as suggested by Raza et al. , then the present study suggests that this increased responsiveness may be quite general in nature, extending to xenobiotics beyond the anthelmintics that were involved in the original resistance-selection process in the field. Importantly though, while WAL larvae respond more rapidly and to a greater extent than Kirby to MPL under our experimental conditions, the two isolates show equivalent sensitivities to the drug in larval development assays . Hence, there is no evidence yet that the increased responsiveness in ABC transporters seen in WAL in the present study provides any subsequent protection against MPL; only protection against IVM has been suggested by our data.
The pattern of transcription changes over the time course of the present experiments was quite different to that observed in our earlier experiments with IVM and LEV . In this earlier study, exposure to either of the two anthelmintics resulted in ABC transporter transcription increases in WAL larvae at the 3 h time point, however, these had returned to control levels either completely (for IVM) or almost completely (with the exception of one gene for LEV) by 6 h. In contrast, in the present study, the upregulation response was sustained over the 24 h time period of the experiments. This difference between the two studies is most likely due to the relative concentrations of the drugs used. The concentrations were chosen on the basis of being the highest concentrations that could be examined without effecting the ability of the larvae to migrate in LMAs in the second phase of the drug-exposure experiments, thereby allowing for the consequences of drug exposure to be assessed using LMAs. As MPL is far less potent as an inhibitor of larval migration than IVM and LEV [20, 22], the higher of the two concentrations used in the present study was much higher than those used previously by Raza et al.  for IVM and LEV (250 μg/ml compared to 0.2–0.8 μg/ml).
It has been suggested that nematodes may use gene upregulation as a mechanism to counter the toxic effects of anthelmintics . There is considerable evidence that anthelmintics act as substrates of ABC transporters and have inducing effects on the expression levels of transporter genes in nematodes and mammals [11, 24, 25]. However, nematode transporter genes do not show a consistent pattern of upregulation following exposure to anthelmintics. Significantly increased transcription of several P-gp genes was reported in wild-type and IVM-resistant isolates of Caenorhabditis elegans post-exposure to moxidectin . On the other hand, exposure of the same isolates of C. elegans to IVM resulted in upregulation of a different set of P-gp genes . Some studies on parasitic nematodes have reported increases in transcription of transporter genes following in vitro and in vivo exposure to IVM [12, 16, 26], whereas no transcriptional changes were observed in IVM-resistant H. contortus and Cooperia oncophora worms collected from animals treated with IVM compared to the worms collected from untreated control animals [27, 28].
We also examined the stability of the observed transcriptional changes by exposing WAL L3 to MPL250 μg/ml for 3 h, rinsing the larvae in water, and then maintaining them in the absence of drug for a further 24 h. We found that only the transcription level of pgp-11remained at significantly higher levels compared to control larvae. We also observed increased tolerance to IVM in a proportion of the WAL larval population 24 h after the removal of MPL following a 3 h pre-treatment period (see Fig. 6). There is very little information available on the time course of ABC transporter gene expression patterns after the removal of the inducing agent. Fardel et al.  reported that expression of P-gp (mdr-1) in rat liver epithelial cells, as measured by qPCR and Western blotting, was significantly higher following exposure to 3-methylcholanthene for 24 h, with a return to almost basal levels 72 h after removal of the inducing agent. The sustained upregulation of pgp-11 alone after removal of the MPL in the present study suggests that this specific ABC transporter may play a more important role in interaction with MPL in WAL larvae than the other transporters. On the other hand, the maintenance of IVM tolerance 24 h after removal of the MPL cannot be linked specifically to the sustained upregulation of pgp-11 as the tolerance may be due to increased efflux activity of other ABC transporter proteins that may be present at higher levels as a result of the earlier period of gene upregulation (from Fig. 2d).
The presence of plateaus in IVM dose-response curves at the highest IVM concentrations was most likely due to the short nature of both the drug incubation and migration phases (30 min for each) of the LMA compared to the 48 and 24 h incubation/migration periods usually used for this assay . As described earlier, the assay was modified to avoid any effects of the drugs in possibly inducing transcription of the transporter genes during long (24 and 48 h) incubation periods. These plateaus were likely due to the subsequent inability of IVM to inhibit migration completely in all larvae in the short 30 min time frame, and were also reported in our earlier study . Tolerance to IVM only occurred in the proportion of the population represented by the dose-response plateau, as the IC50 in the remainder of the population was unchanged compared to DMSO pre-treated L3 (see Table 2). A proportion of the larvae were clearly better equipped to tolerate IVM following MPL250 μg/ml exposure (as indicated by percent migration at plateau and t-test analysis), while the remainder of the population showed no change in IVM sensitivity. Although we have no direct evidence for a role of transporters in the observed IVM tolerance, the known role for ABC transporters in IVM tolerance  in parallel with our gene upregulation/R-123 efflux increases/IVM tolerance observations in the present study, suggest that the MPL-induced transporters were able to subsequently provide protection against IVM in a proportion of the larval population in our experiments. Hence, our data demonstrate the potential role of the transporters in responding to exposure to an anthelmintic, and, in so doing, provide protection against structurally-unrelated compounds. The fact that IVM tolerance was only observed in a proportion of the worm population following MPL pre-exposure may indicate a degree of heterogeneity in the effect of MPL exposure on ABC transporter gene expression across the worm population, as was also suggested by the similar IVM dose-response plateau versus IC50 effects following pre-exposure to IVM or LEV in our previous study .
There were no significant changes in LEV dose-response experiments following exposure to MPL at both concentrations, suggesting that the transporters induced in response to MPL exposure are able to subsequently provide protection against IVM but not LEV. Raza et al.  previously showed that pre-exposure to IVM or LEV resulted in tolerance to IVM, alongside no change in the response to LEV. Hence, the earlier and present studies provide evidence to suggest that IVM may be a better substrate than LEV for H. contortus ABC transporters, as has been reported for mammalian P-gps .
In conclusion, the present study has shown that exposure to a high concentration of MPL increases the transcription of multiple ABC transporter genes in two MPL-susceptible isolates of H. contortus. The most significant interactions were for pgp-11, pgp-12 and pgp-14, based on (i) the magnitude of the transcriptional response, (ii) the occurrence of this response across multiple time points, (iii) the consistency of the response between the two isolates, and (iv) the stability of the response upon removal of the MPL (applicable to pgp-11 only). The subsequent ability of a proportion of MPL-exposed L3 to tolerate higher levels of IVM provides further evidence that ABC transporters play an important role in protection of worms against this anthelmintic, and illustrates the ability of ABC transporters to interact with different chemical entities. However, while the study describes interactions of ABC transporters with a high concentration of MPL in vitro, the effects of MPL on efflux pathways in vivo remains to be determined.
AR is a recipient of an International Postgraduate Research Scholarship from the Australian Government, and a University of Queensland Centennial Scholarship.
Availability of data and material
AR and AK designed the study, interpreted the data and wrote the article, with corrections and suggestions from AJ, SK and NB. NB and AR designed and validated the qPCR primers. AR carried out all the experimental work. All authors read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All animal procedures were approved by the FD McMaster Animal Ethics Committee, CSIRO Agriculture (Animal Ethics Approval Number AEC 13/23).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Lecova L, Stuchlikova L, Prchal L, Skalova L. Monepantel: the most studied new anthelmintic drug of recent years. Parasitology. 2014;141(13):1686–98.View ArticlePubMedGoogle Scholar
- Kaminsky R, Ducray P, Jung M, Clover R, Rufener L, Bouvier J, et al. A new class of anthelmintics effective against drug-resistant nematodes. Nature. 2008;452(7184):176–80.View ArticlePubMedGoogle Scholar
- Rufener L, Maser P, Roditi I, Kaminsky R. Haemonchus contortus acetylcholine receptors of the DEG-3 subfamily and their role in sensitivity to monepantel. PLoS Pathog. 2009;5(4):e1000380.View ArticlePubMedPubMed CentralGoogle Scholar
- Scott I, Pomroy WE, Kenyon PR, Smith G, Adlington B, Moss A. Lack of efficacy of monepantel against Teladorsagia circumcincta and Trichostrongylus colubriformis. Vet Parasitol. 2013;198(1–2):166–71.View ArticlePubMedGoogle Scholar
- Love S. First confirmed case of resistance to new drench in Australia—what does it mean?http://www.wormboss.com.au/news/articles/drench-resistance/first-confirmed-case-of-resistance-to-new-drench-in-australiawhat-does-it-mean.php.Accessed 20 Aug 2015.
- Mederos AE, Ramos Z, Banchero GE. First report of monepantel Haemonchus contortus resistance on sheep farms in Uruguay. Parasit Vectors. 2014;7:598.View ArticlePubMedPubMed CentralGoogle Scholar
- Van den Brom R, Moll L, Kappert C, Vellema P. Haemonchus contortus resistance to monepantel in sheep. Vet Parasitol. 2015;209(3-4):278–80.View ArticlePubMedGoogle Scholar
- Ardelli BF, Prichard RK. Inhibition of P-glycoprotein enhances sensitivity of Caenorhabditis elegans to ivermectin. Vet Parasitol. 2013;191(3-4):264–75.View ArticlePubMedGoogle Scholar
- Lespine A, Ménez C, Bourguinat C, Prichard RK. P-glycoproteins and other multidrug resistance transporters in the pharmacology of anthelmintics: Prospects for reversing transport-dependent anthelmintic resistance. Int J Parasitol Drugs Drug Resist. 2012;2:58–75.View ArticlePubMedGoogle Scholar
- Kotze AC, Hunt PW, Skuce P, von Samson-Himmelstjerna G, Martin RJ, Sager H, et al. Recent advances in candidate-gene and whole-genome approaches to the discovery of anthelmintic resistance markers and the description of drug/receptor interactions. Int J Parasitol Drugs Drug Resist. 2014;4(3):164–84.View ArticlePubMedPubMed CentralGoogle Scholar
- James CE, Davey MW. Increased expression of ABC transport proteins is associated with ivermectin resistance in the model nematode Caenorhabditis elegans. Int J Parasitol. 2009;39(2):213–20.View ArticlePubMedGoogle Scholar
- De Graef J, Demeler J, Skuce P, Mitreva M, Von Samson-Himmelstjerna G, Vercruysse J, et al. Gene expression analysis of ABC transporters in a resistant Cooperia oncophora isolate following in vivo and in vitro exposure to macrocyclic lactones. Parasitology. 2013;140(4):499–508.Google Scholar
- Raza A, Kopp SR, Bagnall NH, Jabbar A, Kotze AC. Effects of in vitro exposure to ivermectin and levamisole on the expression patterns of ABC transporters in Haemonchus contortus larvae. Int J Parasitol Drugs Drug Resist. 2016;6(2):103–15.Google Scholar
- Bygarski EE, Prichard RK, Ardelli BF. Resistance to the macrocyclic lactone moxidectin is mediated in part by membrane transporter P-glycoproteins: Implications for control of drug resistant parasitic nematodes. Int J Parasitol Drugs Drug Resist. 2014;4(3):143–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Heckler RP, Almeida GD, Santos LB, Borges DG, Neves JP, Onizuka MK, Borges FA. P-gp modulating drugs greatly potentiate the in vitro effect of ivermectin against resistant larvae of Haemonchus placei. Vet Parasitol. 2014;205(3-4):638–45.View ArticlePubMedGoogle Scholar
- Tydén E, Skarin M, Höglund J. Gene expression of ABC transporters in Cooperia oncophora after field and laboratory selection with macrocyclic lactones. Mol Biochem Parasitol. 2014;198(2):66–70.View ArticlePubMedGoogle Scholar
- Lloberas M, Alvarez L, Entrocasso C, Virkel G, Ballent M, Mate L, Lanusse C, Lifschitz A. Comparative tissue pharmacokinetics and efficacy of moxidectin, abamectin and ivermectin in lambs infected with resistant nematodes: Impact of drug treatments on parasite P-glycoprotein expression. Int J Parasitol Drugs Drug Resist. 2013;3:20–7.View ArticlePubMedGoogle Scholar
- Albers GA, Burgess SK. Serial passage of Haemonchus contortus in resistant and susceptible sheep. Vet Parasitol. 1988;28(4):303–6.View ArticlePubMedGoogle Scholar
- Love SC, Neilson FJ, Biddle AJ, McKinnon R. Moxidectin-resistant Haemonchus contortus in sheep in northern New South Wales. Aust Vet J. 2003;81(6):359–60.View ArticlePubMedGoogle Scholar
- Raza A, Lamb J, Chambers M, Hunt PW, Kotze AC. Larval development assays reveal the presence of sub-populations showing high- and low-level resistance in a monepantel (Zolvix®)-resistant isolate of Haemonchus contortus. Vet Parasitol. 2016;220:77–82.View ArticlePubMedGoogle Scholar
- Lifschitz A, Ballent M, Virkel G, Sallovitz J, Viviani P, Lanusse C. Accumulation of monepantel and its sulphone derivative in tissues of nematode location in sheep: pharmacokinetic support to its excellent nematodicidal activity. Vet Parasitol. 2014;203(1-2):120–6.View ArticlePubMedGoogle Scholar
- Kotze AC, Le Jambre LF, O'Grady J. A modified larval migration assay for detection of resistance to macrocyclic lactones in Haemonchus contortus, and drug screening with Trichostrongylidae parasites. Vet Parasitol. 2006;137(3-4):294–305.View ArticlePubMedGoogle Scholar
- Ardelli BF, Prichard R. Effects of ivermectin and moxidectin on the transcription of genes coding for multidrug resistance associated proteins and behaviour in Caenorhabditis elegans. J Nematol. 2008;40(4):290–8.Google Scholar
- Godoy P, Lian J, Beech RN, Prichard RK. Haemonchus contortus P-glycoprotein-2: in situ localisation and characterisation of macrocyclic lactone transport. Int J Parasitol. 2015;45(1):85–93.View ArticlePubMedGoogle Scholar
- Menez C, Sutra JF, Prichard R, Lespine A. Relative neurotoxicity of ivermectin and moxidectin in Mdr1ab (-/-) mice and effects on mammalian GABA(A) channel activity. PLoS Neg Trop Dis. 2012;6(11):e1883.View ArticleGoogle Scholar
- Dicker AJ, Nisbet AJ, Skuce PJ. Gene expression changes in a P-glycoprotein (Tci-pgp-9) putatively associated with ivermectin resistance in Teladorsagia circumcincta. Int J Parasitol. 2011;41(9):935–42.View ArticlePubMedGoogle Scholar
- Alvarez L, Suarez G, Ceballos L, Moreno L, Canton C, Lifschitz A, et al. Integrated assessment of ivermectin pharmacokinetics, efficacy against resistant Haemonchus contortus and P-glycoprotein expression in lambs treated at three different dosage levels. Vet Parasitol. 2015;210(1-2):53–63.View ArticlePubMedGoogle Scholar
- Areskog M, Engström A, Tallkvist J, Samson-Himmelstjerna G, Höglund J. PGP expression in Cooperia oncophora before and after ivermectin selection. Parasitol Res. 2013;112(8):3005–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Fardel O, Lecureur V, Corlu A, Guillouzo A. P-glycoprotein induction in rat liver epithelial cells in response to acute 3-methylcholanthrene treatment. Biochem Pharmacol. 1996;51(11):1427–36.View ArticlePubMedGoogle Scholar
- Efferth T, Volm M. Reversal of doxorubicin-resistance in sarcoma 180 tumor cells by inhibition of different resistance mechanisms. Cancer Lett. 1993;70(3):197–202.View ArticlePubMedGoogle Scholar