Skip to main content

Phenotypic and genotypic analysis of benzimidazole resistance in the ovine parasite Nematodirus battus

Abstract

Benzimidazole resistance is common amongst many ovine trichostrongylid nematodes species globally. Although anthelmintics have been used for over half a century in some areas of the world for the control of Nematodirus battus, resistance has never been detected. Veterinary investigations conducted in 2010 demonstrated reduced efficacy in a flock that had been treated previously with fenbendazole (FBZ), suggesting probable resistance in N. battus. Infective larvae (L3; designated MNba2) were generated from the original material to conduct a controlled efficacy test (CET). Faecal egg counts showed an average of 37% reduction in the FBZ treated group 7 days post treatment compared to the untreated lambs. Average worm burden results showed no reduction after FBZ treatment compared to the untreated group (3850 and 3850 worms respectively). A molecular assay to assess the frequency of the commonly associated single nucleotide polymorphisms (SNP) in the β-tubulin isotype 1 gene, F200Y and E198A, was developed. Larval genotypes were predominantly homozygous resistant at codon 200 SNP, ranging from 56%-83% and remained stable at 70% for adult worm populations taken from treated and control lambs in the CET. Only susceptible genotypes were found at codon 198. The allele frequency for F200Y ranged between 80-83% in adult worms taken from the CET from treated and control lambs. The results confirmed initial findings and demonstrated the first report of FBZ resistance in N. battus whilst providing evidence that the P200 point mutation in the β-tubulin isotype 1 gene is a potential mechanism of resistance in the species.

Introduction

There are several species of Nematodirus that cause disease and production losses in young lambs, although most of the losses in the UK are caused by Nematodirus battus[1]-[3]. Nematodirosis is a disease that is normally seen in young lambs, the signs of infection are often acute and, if left unchecked, can result in a high (typically 5-20%) mortality rate [3]. Acute disease is the consequence of very heavy larval challenge and the effects of the developing larvae. As a result, lambs can present with severe scouring and other clinical signs before there are well established levels of adult worms and, hence, eggs appearing in faeces. Currently five classes of broad spectrum anthelmintics are available for use in sheep in the UK. All classes have licensed efficacy against N. battus, however, activity against certain life stages (particularly immatures) can vary by formulation and mode of application [4], particularly within the macrocyclic lactone class. Many farmers in the UK opt to use benzimidazoles due to its wide safety margin and good efficacy [5]. Reports of resistance to this class of anthelmintic are commonplace in ovine trichostrongylid nematode species such as Teladorsagia circumcincta throughout sheep flocks in the UK [6]-[9] and globally [10]-[14]; however it has never previously been reported in N. battus. The only reports of benzimidazole resistance involving Nematodirus spp. (Nematodirus spathiger) are cases of oxfendazole resistance in New Zealand, Australia and Tasmania [15]-[17] and one case of thiabendazole resistance in a mixed species isolate of N. spathiger and Nematodirus filicollis in Australia [18].

Several genetic mechanisms have been associated with BZ resistance in parasitic nematodes of sheep; for example, loss of isotype 2 of the β-tubulin gene and single nucleotide polymorphism (SNP) within isotype 1 of the β-tubulin gene [19],[20]. SNPs in the β-tubulin gene are responsible for an amino acid transversion at each of the sites; phenylalanine to tyrosine at codon 200 (TTC → TAC, F200Y; [21],[22], phenylalanine to tyrosine at 167 (TTT → TAT, F167Y; [23],[24] and glutamic acid to alanine at 198 (GAA → GCA, E198A; [25]). Not all of these mutations are found in all ovine parasitic nematode species that are resistant and the presence of one SNP is not usually accompanied by a second [26]. At present, no information relating to the mechanisms involved in benzimidazole resistance in N. battus is available.

In 2010, a clinical case of nematodirosis was submitted to Animal Health Veterinary Laboratory Agency (AHVLA) Penrith where a farmer was concerned about a group of 60 lambs showing signs of disease (diarrhoea and ill thrift); three lambs had died from the group, with a fourth being submitted for post mortem examination. Nematodirosis due to N. battus was diagnosed as the cause of death due to the high numbers of N. battus eggs detected in the faeces (1850 eggs per gram) and a small number (300) of adults present in the small intestinal contents. No other gastro-intestinal nematodes were identified in the abomasal or small intestinal contents. Treatment of the surviving lambs with FBZ resulted in little clinical improvement; however the farmer only treated the lambs to the average weight of the group. After seeking further advice, a faecal egg count reduction test (FECRT) was conducted by the attending veterinarian using FBZ at the manufacturer’s recommended dose rate (MRDR) and faecal samples were submitted to the AHVLA at the time of treatment (designated hereafter as MNba2VCF) and 10 days later [27].

Here, the authors report the findings from a controlled efficacy test (CET) performed at Moredun Research Institute using the field isolate of N. battus derived from the MNba2VCF generated material. Secondly, they describe a pyrosequencing assay developed to identify polymorphisms at codons 198 and 200 of the β-tubulin isotype 1 gene to investigate the potential mechanisms involved in BZ resistance in N. battus.

Materials and methods

Parasite isolate

Faecal material was collected from lambs following an on-farm FBZ FECRT which indicated the presence of benzimidazole resistance [27]. N. battus eggs were cultured and extracted from faeces using a modification of methodologies described in the Manual of Veterinary Parasitology Techniques reference book [28]. In brief, eggs were extracted from faeces by differential sieving and salt flotation. Recovered eggs were stored in tap water in 75 cm2 vented culture flasks (Corning B.V. Life Sciences, Amsterdam, The Netherlands) at around 20 °C for over 7 weeks to allow full development of larvae in the eggs. To assist hatching, eggs were sandwiched between two glass plates and gently crushed by applying downward pressure until a faint cracking sound was heard. Successful hatching was assessed by viewing the eggs down a stereo microscope (×20) before Baermannising the suspension. The resultant third stage larvae (L3 n = 4300) were subsequently passaged through a parasite-naïve lamb, which was administered with FBZ at day 43 post-infection (pi) at MRDR. Faecal material was collected before and after treatment and processed as above to generate sufficient L3 for the controlled efficacy trial. Eggs recovered post-treatment were artificially hatched and the L3 used in the controlled efficacy trial.

Experimental design

Ten helminth-free, female, six month-old Texel X Greyface lambs were each artificially infected per os with 6000 L3 (day 0 pi). On day 24 pi, faecal egg counts (FEC) were conducted using a modification of the salt flotation technique described by [29] and [30] with a sensitivity of up to one egg per gram. All lambs were weighed. The lambs were then allocated into one of two groups ensuring they were balanced as closely as possible for both FEC and weight. On day 25 pi, one group was orally administered fenbendazole (Panacur 2.5%, Intervet, Milton Keynes, UK; 5 mg/kg body weight (BW)), whilst the second group remained untreated. These groups were designated MNba2FBZ and MNba2CON, respectively. All anthelmintic treatment doses were calculated according to the respective manufacturer’s instructions, with doses rounded up to the nearest 0.5 mL (dosage range 5.0 - 5.3 mg/kg BW). Faecal samples were taken per rectum daily from each lamb from day 14 pi until day 32 pi and FECs conducted. All lambs were slaughtered on day 32 pi and the small intestines removed for saline digest to estimate total worm burdens. Methods used for necropsy and worm recovery were as described previously [31]. Total burdens of each animal were estimated from a 2% subsample (100 mL) of the small intestinal wash and digest. Worms were separated into adult males, adult females and juvenile stages using the criteria described in the Manual of Veterinary Parasitology Techniques reference book [28].

All experimental procedures conducted at Moredun Research Institute were assessed and approved by the Institute’s Experiments and Ethics Committee and were conducted under the legislation of a UK Home Office License (reference PPL 60/03899) in accordance with the Animals (Scientific Procedures) Act of 1986.

Statistical analysis

Nematode burdens and FECs were square-root transformed to normalize for variance. Burdens were compared using one way ANOVA (Minitab version 13, Coventry, UK), followed by Fisher’s pairwise comparisons when found to be significant (p < 0.05). Statistical significance was accepted as p < 0.05. Fenbendazole efficacy was calculated based on group mean FECs using one of a range of standard formulae, (1 − [T2/C2]) × 100 using arithmetic means [32], (1 − [T2/T1] [C1/C2]) × 100 using geometric means [33], (1 − [T2/T1] [C1/C2]) × 100 using arithmetic means [34], (1 − [T2/T1]) × 100 using arithmetic means [35], where C1 and C2 are the FEC of untreated control animals pre- and post treatment respectively and T1 and T2 are the FEC of animals pre- and post treatment respectively. Bootstrap analysis of the data using the “BootStreat” program was also conducted, with a re-sampling number of 2000 to calculate mean treatment efficacies and upper and lower 95% confidence limits [36]. Efficacy based on total worm burdens was calculated using the WAAVP guidelines formula [32].

Molecular analysis

Due to lack of published sequences available for N. battus β-tubulin gene, partial β-tubulin isotype 1 sequences were amplified from cDNA generated from total RNA that had been extracted using TRIzol reagent (Life Technologies, California, USA) from 100, 000 L3 from a known FBZ-sensitive isolate (designated MNba1) and MNba2VCF isolate L3 obtained after FBZ treatment along with a FBZ susceptible isolate of Haemonchus contortus (MHco3; [37]), as a control. Twenty microlitres of cDNA were generated using Invitrogen superscript III reverse transcriptase kit (Invitrogen, California, USA), where 4 μL of RNA extract was used per sample. Generic primers (Table 1: Gen Beta-tub For1 & Gen Beta-tub Rev1) were designed to amplify an area of the isotype 1 β-tubulin gene that covered the three most common SNP’s of interest using published sequences from other trichostrongylid nematode species; Trichostrongylus colubriformis (Accession number: L23506), Teladorsagia circumcincta (Accession number - Z69258), H. contortus (Accession number - M76491/EF198865) and Cooperia oncophora (Accession number - AY259994). PCR products were amplified using Novataq™ Hot start master mix (Novagen, Madison, USA) in a 50 μL reaction, with primers and MgCl2 at a final concentration of 0.3 μM and 3 mM, respectively, with the following cycle conditions: 94 °C for 10 min, (94 °C for 30 s, 52 °C for 30 s, 72 °C for 30 s) for 35 cycles and 72 °C for 10 min. Beta-tubulin isotype 1 sequences spanning the SNPs of interest were also generated from genomic DNA (gDNA) isolated from 10 individual adult N. battus worms (5 MNba1 and 5 MNba2CON). One microlitre of gDNA was added to a 25 μL PCR reaction using HcPy2 PCR For and HcPy2PCR Rev primers [38], see Table 1) at 0.3 μM and 3.0 mM MgCl2 final concentrations in NovaTaq™ Hotstart master mix. PCR amplicons of the appropriate size were cut from a 1.5% agarose gel (Gel Extraction Kit, Qiagen, Hilden, Germany) and the gDNA ligated into the pGEM-T vector (Promega, Madison, USA) following the manufacturer’s protocol. Transformations were carried out following the standard protocol for JM109 Escherichia coli Competent Cells (Promega) and two (for gDNA sequences) and six (for cDNA sequences) transformed colonies were picked and grown up in 10 mL LB broth overnight before DNA extraction was carried out using a Wizard SV mini prep kit (Promega). DNA concentrations for gDNA and cDNA were assessed using a Nanodrop spectrophotometer (Nanodrop Technologies Inc., Delaware, USA) and adjusted if necessary using DNA/RNA-free water according to Eurofins MWG requirements for sequencing. DNAStar (Madison, USA), Lasergene 9 EditSeq, BLAST and ClustalW were used to analyse and align sequences. Consensus N. battus sequence was obtained using the gDNA sequences and used to design primers that spanned codons 200 and 198 for a pyrosequencing assay using “Pyrosequencing assay design” software, Version 1.0 (Table 1).

Table 1 Primers used in this study (*invitrogen oligo perfect designer used; # pyrosequencing assay design software version 1.0 used;$primer sequences obtained from[38]

Pyrosequencing

To provide gDNA for pyrosequencing, 20 ethanol-fixed adult N. battus were used from three populations; MNba1, MNba2CON and MNba2FBZ. Genomic DNA was extracted using DNeasy Kit (Qiagen) as per kit instructions. For PCR, the adult gDNA extracts were diluted 1:5 in DNA/RNA-free water (Sigma, Dorset, UK). For provision of gDNA from four populations of L3 (MNba1, MNba2VCF, MNba2CON and MNba2FBZ), ethanol fixed L3 were bathed in PBS for 15–30 min prior to transferring individuals into 30 μL of lysis buffer (50 mM KCl, 10 mM Tris (pH 8.3), 2.5 mM MgCl2, 0.45% (v/v) Nonidet P-40, 0.45% (v/v) Tween 20, 0.01% (w/v) gelatin and 0.1 mg/mL Proteinase K) in a 96-well plate (Axygen, California, USA); frozen for 30 min at −20 °C, prior to incubating at 56 °C overnight. Ethanol precipitation was carried out to clean up the gDNA and re-supsended in 25 μL of DNA/RNA-free water (as above). A minimum of 80 L3 were lysed for each study population. For pyrosequencing, 4 μL of gDNA from both adult and L3 were amplified in a mix containing 0.185 μM Nb B-t200 Fbio, 0.2 μM Nb B-t200 Rev, 4.5 mM MgCl2, 25 μL 2 × buffer and made up to 50 μL using DNA/RNA-free water (NovaTaq™ Hot start master mix, Novagen). No template controls were included on each plate. Amplification was performed, following a 15 min 95 °C polymerase activation step, for 45 cycles at 94 °C for 30 s, 58 °C for 30 s and 72 °C for 30 s, followed by a final extension step at 72 °C for 10 min. The pyrosequencing P200/198 assay was conducted according to the manufacturer’s (Qiagen) protocols using a PyroMark ID instrument. In brief, following PCR amplification, 40 μL of the reaction was transferred into a 96-well plate and stored at 4 °C until the remaining PCR reaction (10 μL) was assessed by gel electrophoresis on 2% agarose gels stained with gel red (Biotium, California, USA). The plate was run on the Pyromark ID instrument if clear PCR amplification was observed.

Genotyping/allele frequency

Genotypes (TTC/TTC, TTC/TAC and TAC/TAC) and allele frequencies for resistant (R) and susceptible alleles (S) from each population are expressed as percentages [39].

Chi–square tests and Phi coefficient association of P200 genotypes were conducted using Minitab version 15.

Results

Faecal egg count reduction test analysis

All individual lambs were excreting N. battus eggs by day 14 pi. The FEC counts continued to rise in most animals until day 24 pi (Figure 1). Arithmetic mean FECs at day 24 pi were 265 (range 198 – 338 EPG) and 277 (range 216 – 360 EPG) EPG for the MNba2CON and MNba2FBZ groups, respectively. At post mortem the mean FEC of both groups remained similar, with 167 and 172 EPG for the FBZ treated and control groups, respectively, Table 2. No significant difference was observed between the groups.

Figure 1
figure 1

Faecal egg count profile of lambs infected with 6,000 Nematodirus battus infective larvae. Arithmetic mean faecal egg counts (± standard error of the mean) of two groups of lambs infected with 6000 Nematodirus battus infective larvae. On day 25 post infection (↑) one group was orally administered fenbendazole (□; Panacur, Intervet; 5 mg/kg BW), whilst the second group remained untreated (■).

Table 2 Arithmetic mean (± S.E.M) [range] faecal egg counts of untreated control and fenbendazole treated lambs on day of treatment and seven days post treatment

CET

Arithmetic mean percentage establishment of N. battus in control lambs (MNba2CON) was 64% (arithmetic mean = 3850 worms), with the estimated total number of worms ranged between 2350 and 5750 (Figure 2) within the group. Arithmetic mean worm burden of the FBZ treated (MNba2FBZ) group was 3850 worms (range, 3050 – 4500 worms) and hence the treatment efficacy was 0%.

Figure 2
figure 2

Individual estimated total worm burden of the two groups of lambs infected with 6000 Nematodirus battus infective larvae. The stacked bar chart shows the estimated total worm burden (reported as males (light grey), females (dark grey) and juveniles (black)) for each lamb in the fenbendazole (FBZ) treated (MNba2FBZ; 5 mg/kg FBZ body weight) and control (MNba2CON) groups. The control group (MNba2CON) estimated total number of worms ranged between 2350 and 5750, with an arithmetic mean of 3850 worms. The FBZ treated group (MNba2FBZ) estimated total number of worms ranged between 3050 – 4500 worms, with an arithmetic mean of 3850 worms. This showed a 0% efficacy with FBZ treatment.

Molecular analysis

Sequencing results

A 785 bp fragment was amplified from the cDNA from the resistant (MNba2VCF) and susceptible (MNba1) isolates of N. battus. Sequences from MNba1 were subjected to a Blast search, with sequence homologies of 94% with T. circumcincta (Z69258) and 85% with C. oncophora (AY259994) isotype 1 sequences with hits also showing homology with, H. contortus, O. ostertagi and T. colubriformis isotype 1 β-tubulin sequences. Consensus sequence derived from the MNba2VCF isolate showed 87% homology to H. contortus and C. oncophora (AY259994) isotype 1 β-tubulin mRNA sequences. A smaller fragment of approximately 340 bp was amplified from genomic DNA. The consensus sequence from MNba1 gDNA blast search showed 86% and 85% homology with H. contortus (FJ981633) and T. circumcincta (Z69258) isotype 1 β-tubulin sequences respectively. Similar results were produced with MNba2FBZ consensus sequence.

Genotyping

Cross reactivity of the isotype 1 β-tubulin pyrosequencing assay was examined against other species, N. filicollis, T. circumcincta, H. contortus and T. colubriformis. DNA from these species did not provide amplification using the Nb B-t 200 forward and reverse primers. The mean percentage genotypes (homozygous TTC/TTC, heterozygous TTC/TAC and homozygous TAC/TAC) for L3 and adult worms from the MNba1 and MNba2 populations are shown in Tables 3 and 4, respectively. For all samples tested, the MNba2 isolate was predominantly homozygous resistant genotype (TAC/TAC) at codon 200: for example, 56% of L3 (MNba2VCF) were homozygous resistant in the population obtained from the on-farm FECRT and 83% of L3 (MNba2CON) from the CET control population were homozygous resistant. In the case of the adult worms obtained from the CET study, 70% of all worms tested in the (MNba2CON and MNba2FBZ) populations were identified as homozygous resistant at codon 200.

Table 3 Beta-tubulin codon 200 genotyping results of the infective larvae of two Nematodirus battus isolates (MNba1 FBZ-sensitive and MNba2, FBZ-resistant phenotypes)
Table 4 Beta-tubulin codon 200 genotyping results from Nematodirus battus adult worms of two Nematodirus battus isolates (MNba1 FBZ-sensitive and MNba2, FBZ-resistant phenotypes)

Both L3 and adult worms derived from the MNba2FBZ samples showed a slightly higher number of homozygous susceptible genotypes (TTC/TTC) compared to the MNba2CON population; the Chi-square test showed that there was no significant difference (adults - X2 = 0.44, df = 2, p = 0.801 and L3X2 = 1.025, df = 2, P = 0.599) between the genotypes of these populations. The test showed a difference between the P200 genotypes identified, when worms derived from MNba2VCF and MNba2CON (X2 = 14.833, df = 2, p-value = 0.001) were compared. No individuals from isolate MNba1 were homozygous (TAC/TAC) or heterozygous (TTC/TAC) resistant at codon 200. Only homozygous “susceptible” genotypes (GAA/GAA) were identified at codon198 for both adults and larval populations from all the N.battus isolates examined here.

The R allele frequency for 200Y ranged from 85% - 90% within the MNba2 L3 populations (MNba2VCF, MNba2CON and MNba2FBZ; see Table 3). The R allele frequency was 83% and 80% for MNba2CON and MNba2FBZ adults, respectively (see Table 4).

Discussion

The results here provide confirmation, based on residual worm burden and faecal egg count analysis post benzimidazole treatment, of the preliminary results of the first reported case of BZ resistance in N. battus globally [27]. The zero percentage efficacy noted here, poses the question as to whether any members of the benzimidazole class of anthelmintics would be effective against this isolate. In this context, side resistance within the BZ class has been shown to occur in T. colubriformis, H. contortus[40],[41] and undifferentiated Nematodirus species [18] and may be expected to occur with N. battus. Currently, the most practical method of determining anthelmintic sensitivity is the FECRT. There are a number of drawbacks to this test; it is time-consuming, requiring 10–14 days to complete [42], is relatively insensitive at low levels of resistance [43] and is potentially prone to misdiagnosis if the window of opportunity prior to re-infection is missed. With Nematodirus, in particular, the correlation between FEC and worm burden is particularly poor in treated [44] and untreated lambs [45],[46]. Nematodirus battus historically was characterised as being fully developed by 14 days pi, but that maximum egg production did not occur until slightly later [47]. The results based on the isolate here suggest that significant FEC may be detectable in faeces on or before 14 days pi, even at moderate infection rates. The current findings are contrary to those found in New Zealand where in some cases of BZ resistance in N. spathiger, adult worms were recovered in animals that had zero FEC after BZ administration [48]. If the disparity between FEC and worm burden can occur with N. battus then it may explain the lack of detection of resistance in this parasite in the past and would have a bearing on detection of resistance in the field.

To implement sustainable control strategies it is essential to identify the presence of resistance at an early stage of development. In vitro detection of benzimidazole resistance in Nematodirus species using traditional methodologies such as the egg hatch test is complicated by the parasites development and hatching behavior as well as the egg shell morphology [16]. Therefore to better understand the acquisition and development of FBZ resistance in N. battus the authors here strove to identify whether any of the “commonly associated” point mutations of the β-tubulin isotype 1 gene were present in the FBZ resistant population. The results from MNba2 only yielded the F200Y SNP and found that there was no changes observed at E198A. The relationship between F200Y mutations in N. battus and actual anthelmintic efficacy in the field requires further examination, but numerous studies have shown that in other trichostrongylid species there is a high degree of correlation between phenotypic BZ characterisation and resistance allele frequencies of various life stages [49]-[52]. Researchers have found that H. contortus individuals which are homozygous resistant at codon 198 have a higher level of phenotypic resistance than those with SNP mutations at codon 200 [51]. Further work is required in N. battus to investigate other SNPs in the isotype 1 gene and also to characterize isotype 2 of the β-tubulin gene as it is possible that this is present in the populations here as indicated by amplification of some sequences with high identity to this isotype (data not shown).

The comparison of the genotyping results obtained from the MNba2 isolate over the course of the study, starting from the initial field isolate (MNba2VCF) to the CET populations (MNba2CON and MNba2FBZ), showed that the two treatments with FBZ at the MRDR (one in the donor animal that produced the L3 for the CET and one as part of the CET) only marginally led to a difference in the F200Y allele frequency (85, 90 and 87% respectively), which would indicate that this population is stable and highly selected. The observation of a small number of RS and SS genotypes in adult survivors following FBZ treatment, suggests either that the worms are phenotypically sensitive but were able to “hide” from the lethal effects of treatment and then recover/resume activities as the local anthelmintic concentrations fall over time, or that there are other potential mechanisms of resistance involved in the BZ resistant phenotype. These findings are in agreement with other characterisation studies where analysis of the genotypes of first stage T.circumcincta larvae survivors from an in vitro egg hatch test, at concentrations up to 2 μg/mL thiabendazole, still had around 5% homozygous susceptible and 15% heterozygous genotypes [53].

As to why this isolate of N. battus developed FBZ resistance when previously the species has appeared to be susceptible to treatment is unclear and may be the result of a number of factors. It is possible that detection on this farm was aided by greater awareness of anthelmintic resistance and good co-operation between the farmer and his veterinarian. Benzimidazoles had been used on the farm to control N. battus since 2007 until resistance was diagnosed in 2010. Inadvertent under-dosing could have occurred as lamb weight was estimated and they did not normally calibrate their dosing equipment. Reduced bioavailablity due to rapid gut flow may also have contributed to under-dosing [54],[55] either from the clinical effects of N. battus or other concurrent causes of diarrhoea, such as coccidiosis. One of the major pathogenic effects of N. battus is attributable to disruption of the intestinal mucosa and the villous atrophy associated with larval stage development. The result of this atrophy is a reduced capability for fluid exchange, acute diarrhoea and increased gut flow [56]. Under dosing has been shown to allow heterozygote resistant individuals to survive treatment and contribute genes for resistance to the subsequent populations [57],[58].

In some years benzimidazole treatments had been repeated at approximately monthly intervals in the young lambs on this farm dependent on the farmer and veterinary surgeons perception of risk of disease. Prophylactic treatments three weeks apart have been recommended in high risk areas to coincide with the predicted timing of peak hatch [56]. These high treatment frequencies place a considerable selection pressure on the population and could have, over time, resulted in the development of a resistant population.

As common on many sheep farms successive lamb crops grazed the same permanent pastures year on year. This in combination with frequent treatment administrations may have resulted in the slow increase in gene frequency for FBZ resistance within the population. There was no history of dose and movement of lambs onto new pastures. Previous studies have shown that rapid selection of BZ resistance in nematode species has occurred following dose and move treatment strategies [59]-[61], if lambs are grazed on pastures that had been previously been “seeded” with eggs of survivors of FBZ treatment, the possibility exist that resistance could be selected and propagated [60],[62].

Changes in the parasite population brought about by changes in environmental/climatic conditions leading to changes in human behaviour, land use and/or animal husbandry [63]-[65] could be involved. Changes in climate have been shown to lead to conditions favourable for longer grazing periods [66]. The potential for longer parasite seasons will inevitably result in subsequent changes in treatment patterns associated with controlling infections leading to possible parasite adaptation.

However, no significant risk factor(s) can be identified on this farm with respect to anthelmintic resistance development in N. battus that would not be present on a large number of other sheep farms in the UK. It may be that within N. battus, the polymorphism associated with FBZ resistance conferred a fitness cost and until the introgression of resistance gene(s) within a population was sufficiently stable, the population was unable to thrive and propagate itself leading to a sporadic appearance on the farm. Previous works conducted on other nematode species, namely H. contortus and T. colubriformis, have observed poorer development and survival of eggs in anthelminitic resistant populations at a variety of temperatures [67],[68].

An understanding of whether selection for BZ resistance within this population has occurred as a result of pre-adaptation, spontaneous mutation or gene flow will provide a better opportunity to develop effective sustainable strategies for control and now requires further investigation. This study highlights that for sustainable control of N. battus, farmers need to consider monitoring treatment efficacies, to minimize treatment frequency where possible and avoid indiscriminate use of anthelmintic compounds.

References

  1. Thomas RJ: Field studies on the seasonal incidence of Nematodirus battus and N. filicollis in sheep. Parasitology. 1959, 49: 387-410. 10.1017/S0031182000026937.

    Article  CAS  PubMed  Google Scholar 

  2. Kingsbury PA:Nematodirus infestation - a probable cause of losses amongst lambs. Vet Rec. 1953, 65: 167-169.

    Google Scholar 

  3. Thomas RJ, Stevens AJ: Some observations on Nematodirus disease in Northumberland and Durham. Vet Rec. 1956, 68: 471-475.

    Google Scholar 

  4. NOAH Compendium. [http://www.noahcompendium.co.uk/Compendium/Overview/-21789.html].

  5. Abbott KA, Taylor MA, Stubbings LA (2012) A technical manual for veterinary surgeons and advisors. Sustainable Control of Parasites in Sheep. 4th edition. http://www.scops.org.uk/vets-manual.html,

  6. Britt DP: Benzimidazole-resistant nematodes in Britain. Vet Rec. 1982, 110: 343-344. 10.1136/vr.110.14.343.

    Article  CAS  PubMed  Google Scholar 

  7. Cawthorne RJ, Whitehead JD: Isolation of benzimidazole resistant strains of Ostertagia circumcincta from British sheep. Vet Rec. 1983, 112: 274-277. 10.1136/vr.112.12.274.

    Article  CAS  PubMed  Google Scholar 

  8. Scott EW, Duncan JL, McKellar QA, Coop RL, Jackson F, Mitchell GB: Benzimidazole resistance in sheep nematodes. Vet Rec. 1991, 128: 618-619. 10.1136/vr.128.26.618.

    Article  CAS  PubMed  Google Scholar 

  9. Bartley DJ, Jackson E, Johnston K, Coop RL, Mitchell GBB, Sales J: A survey of anthelmintic resistant nematode parasites in Scottish sheep flocks. Vet Parasitol. 2003, 117: 61-71. 10.1016/j.vetpar.2003.07.023.

    Article  PubMed  Google Scholar 

  10. Vlassoff A, Kettle PR: Benzimidazole resistance in Haemonchus contortus. N Z Vet J. 1980, 28: 23-24. 10.1080/00480169.1980.34683.

    Article  Google Scholar 

  11. McKenna PB: The identity of nematode genera involved in cases of ovine anthelmintic resistance in the southern north island of New Zealand. N Z Vet J. 1995, 43: 225-227. 10.1080/00480169.1995.35897.

    Article  CAS  PubMed  Google Scholar 

  12. Waller PJ: Anthelmintic resistance. Vet Parasitol. 1997, 72: 391-405. 10.1016/S0304-4017(97)00107-6.

    Article  CAS  PubMed  Google Scholar 

  13. Howell SB, Burke JM, Miller JE, Terrill TH, Valencia E, Williams MJ: Prevalence of anthelmintic resistance on sheep and goat farms in the southeastern United States. J Am Vet Med Assoc. 2008, 233: 1913-1919. 10.2460/javma.233.12.1913.

    Article  PubMed  Google Scholar 

  14. Wolstenholme AJ, Fairweather I, Prichard R, von Samson-Himmelstjerna G, Sangster NC: Drug resistance in veterinary helminths. Trends Parasitol. 2004, 20: 469-476. 10.1016/j.pt.2004.07.010.

    Article  CAS  PubMed  Google Scholar 

  15. Middleberg A, McKenna PB: Oxfendazole resistance in Nematodirus spathiger. N Z Vet J. 1983, 31: 65-66. 10.1080/00480169.1983.34971.

    Article  Google Scholar 

  16. Obendorf DL, Parsons J, Nicholls J: An egg development test for the evaluation of benzimidazole resistance in Nematodirus spathiger. Aust Vet J. 1986, 63: 382-383. 10.1111/j.1751-0813.1986.tb02911.x.

    Article  CAS  PubMed  Google Scholar 

  17. Obendorf DL, Nicholls J, Koen T, Lacy E: Benzimidazole resistant Nematodirus sp in Tasmania. Aust Vet J. 1991, 68: 72-73. 10.1111/j.1751-0813.1991.tb03142.x.

    Article  CAS  PubMed  Google Scholar 

  18. Martin PJ, Anderson N, Jarrett RG: Resistance to benzimidazole anthelmintics in field strains of Ostertagia and Nematodirus in sheep. Aust Vet J. 1985, 62: 38-43. 10.1111/j.1751-0813.1985.tb14230.x.

    Article  CAS  PubMed  Google Scholar 

  19. Kwa MSG, Veenstra JG, Roos MH: Molecular characterization of beta-tubulin genes present in benzimidazole resistant populations of Haemonchus contortus. Mol Biochem Parasitol. 1993, 60: 133-144. 10.1016/0166-6851(93)90036-W.

    Article  CAS  PubMed  Google Scholar 

  20. Kwa MSG, Kooyman FNJ, Boersema JH, Roos MH: Effect of selection for benzimidazole resistance in Haemonchus contortus on beta-tubulin isotype-1 and isotype-2 genes. Biochem Biophys Res Commun. 1993, 191: 413-419. 10.1006/bbrc.1993.1233.

    Article  CAS  PubMed  Google Scholar 

  21. Kwa MS, Veenstra JG, Roos MH: Benzimidazole resistance in Haemonchus contortus is correlated with a conserved mutation at amino acid 200 in beta-tubulin isotype 1. Mol Biochem Parasitol. 1994, 63: 299-303. 10.1016/0166-6851(94)90066-3.

    Article  CAS  PubMed  Google Scholar 

  22. Elard L, Humbert JF: Importance of the mutation of amino acid 200 of the isotype 1 beta-tubulin gene in the benzimidazole resistance of the small ruminant parasite Teladorsagia circumcincta. Parasitol Res. 1999, 85: 452-456. 10.1007/s004360050577.

    Article  CAS  PubMed  Google Scholar 

  23. Silvestre A, Cabaret J: Mutation in position 167 of isotype 1 beta-tubulin gene of trichostrongylid nematodes: role in benzimidazole resistance?. Mol Biochem Parasitol. 2002, 120: 297-300. 10.1016/S0166-6851(01)00455-8.

    Article  CAS  PubMed  Google Scholar 

  24. Prichard RK: Genetic variability following selection of Haemonchus contortus with anthelmintics. Trends Parasitol. 2001, 17: 445-453. 10.1016/S1471-4922(01)01983-3.

    Article  CAS  PubMed  Google Scholar 

  25. Ghisi M, Kaminsky R, Maser P: Phenotyping and genotyping of Haemonchus contortus isolates reveals a new putative candidate mutation for benzimidazole resistance in nematodes. Vet Parasitol. 2007, 144: 313-320. 10.1016/j.vetpar.2006.10.003.

    Article  CAS  PubMed  Google Scholar 

  26. Silvestre A, Humbert JF: Diversity of benzimidazole resistance alleles in populations of small ruminant parasites. Int J Parasitol. 2002, 32: 921-928. 10.1016/S0020-7519(02)00032-2.

    Article  CAS  PubMed  Google Scholar 

  27. Mitchell S, Mearns R, Richards I, Donnan AA, Bartley DJ: Benzimidazole resistance in Nematodirus battus. Vet Rec. 2011, 168: 623-624. 10.1136/vr.d3584.

    Article  PubMed  Google Scholar 

  28. MAFF: Ministry of Agriculture, Fisheries and Food, Manual of Veterinary Parasitological Laboratory Techniques, Reference Book 418. 1986, Her Majesty’s Stationery Office (HMSO), London

    Google Scholar 

  29. Christie M, Jackson F: Specific identification of strongyle eggs in small samples of sheep faeces. Res Vet Sci. 1982, 32: 113-117.

    CAS  PubMed  Google Scholar 

  30. Jackson F: New technique for obtaining nematode ova from sheep faeces. Lab Pract. 1974, 23: 65-66.

    CAS  PubMed  Google Scholar 

  31. Patterson DM, Jackson F, Huntley JF, Stevenson LM, Jones DG, Jackson E: Studies on caprine responsiveness to nematodiasis: segregation of male goats into responders and non-responders. Int J Parasitol. 1996, 26: 187-194. 10.1016/0020-7519(95)00121-2.

    Article  CAS  PubMed  Google Scholar 

  32. Coles GC, Bauer C, Borgsteede FHM, Geerts S, Klei TR, Taylor MA, Waller PJ: World association for the advancement of veterinary parasitology (WAAVP) methods for the detection of anthelmintic resistance in nematodes of veterinary importance. Vet Parasitol. 1992, 44: 35-44. 10.1016/0304-4017(92)90141-U.

    Article  CAS  PubMed  Google Scholar 

  33. Presidente PJA: Methods for detection of resistance to anthelmintics. Resistance in Nematodes to Anthelmintic Drugs. Edited by: Anderson N, Waller PJ. 1985, CSIRO Division of Animal Health, and Australian Wool Corporation, Glebe, NSW, 13-27.

    Google Scholar 

  34. Dash KM, Hall E, Barger IA: The role of arithmetic and geometric mean worm egg counts in faecal egg count reduction tests and in monitoring strategic drenching programs in sheep. Aust Vet J. 1988, 65: 66-68. 10.1111/j.1751-0813.1988.tb07359.x.

    Article  CAS  PubMed  Google Scholar 

  35. McKenna PB: The detection of anthelmintic resistance by the fecal egg count reduction test - an examination of some of the factors affecting performance and interpretation. N Z Vet J. 1990, 38: 142-147. 10.1080/00480169.1990.35640.

    Article  CAS  PubMed  Google Scholar 

  36. Cabaret J: Reliable Phenotypic Evaluation of Anthelmintic Resistance in Herbivores: how and When Should They be Done?. Anthelmintics, Clinical Pharmacology, Uses in Veterinary Medicine and Efficacy. Edited by: Quick W. 2014, Nova science publisher, Nova, New York, 1-26.

    Google Scholar 

  37. Roos MH, Otsen M, Hoekstra R, Veenstra JG, Lenstra JA: Genetic analysis of inbreeding of two strains of the parasitic nematode Haemonchus contortus. Int J Parasitol. 2004, 34: 109-115. 10.1016/j.ijpara.2003.10.002.

    Article  CAS  PubMed  Google Scholar 

  38. von Samson-Himmelstjerna G, Walsh TK, Donnan AA, Carriere S, Jackson F, Skuce PJ: Molecular detection of benzimidazole resistance in Haemonchus contortus using real-time PCR and pyrosequencing. Parasitology. 2009, 136: 349-358. 10.1017/S003118200800543X.

    Article  CAS  PubMed  Google Scholar 

  39. Falconer DS, MacKay TFC: Introduction to quantitative genetics. 1996, Pearson Education Ltd, Harlow, 4

    Google Scholar 

  40. Hall CA, Kelly JD, Campbell NJ, Whitlock HV, Martin IC: The dose response of several benzimidazole anthelmintics against resistant strains of Haemonchus contortus and Trichostrongylus colubriformis selected with thiabendazole. Res Vet Sci. 1978, 25: 364-367.

    CAS  PubMed  Google Scholar 

  41. Hall CA, Campbell NJ, Richardson NJ: Levels of benzimidazole resistance in Haemonchus contortus and Trichostrongylus colubriformis recorded from an egg hatch test procedure. Res Vet Sci. 1978, 25: 360-363.

    CAS  PubMed  Google Scholar 

  42. Coles GC, Jackson F, Pomroy WE, Prichard RK, von Samson-Himmelstjerna G, Silvestre A: The detection of anthelmintic resistance in nematodes of veterinary importance. Vet Parasitol. 2006, 136: 167-185. 10.1016/j.vetpar.2005.11.019.

    Article  CAS  PubMed  Google Scholar 

  43. Martin PJ, Anderson N, Jarrett RG: Detecting benzimidazole resistance with fecal egg count reduction tests and in vitro assays. Aust Vet J. 1989, 66: 236-240. 10.1111/j.1751-0813.1989.tb13578.x.

    Article  CAS  PubMed  Google Scholar 

  44. Black WJM: The development of a preventative routine against Nematodirus disease of lambs I. Administration of bephenium compounds at 21-day intervals. Br Vet J. 1964, 120: 301-

    Google Scholar 

  45. Kingsbury PA: Relationship between egg counts and worm burdens of young sheep. Vet Rec. 1965, 77: 900-901. 10.1136/vr.77.31.900.

    Article  CAS  PubMed  Google Scholar 

  46. McKenna PB: The diagnosis value and interpretation of faecal egg counts in sheep. N Z Vet J. 1981, 29: 129-132. 10.1080/00480169.1981.34821.

    Article  CAS  PubMed  Google Scholar 

  47. Thomas RJ: A comparative study of the life histories of Nematodirus battus and N. filicollis, nematode parasites of sheep. Parasitology. 1959, 49: 374-386. 10.1017/S0031182000026925.

    Article  CAS  PubMed  Google Scholar 

  48. Chalmers K: Detection of benzimidazole resistant Nematodirus spathiger. N Z Vet J. 1985, 33: 53-10.1080/00480169.1985.35160.

    Article  CAS  PubMed  Google Scholar 

  49. von Samson-Himmelstjerna G, von Witzendorff C, Sievers G, Schnieder T: Comparative use of faecal egg count reduction test, egg hatch assay and beta-tubulin codon 200 genotyping in small strongyles (cyathostominae) before and after benzimidazole treatment. Vet Parasitol. 2002, 108: 227-235. 10.1016/S0304-4017(02)00197-8.

    Article  CAS  PubMed  Google Scholar 

  50. Drogemuller M, Failing K, Schnieder T, von Samson-Himmelstjerna G: Effect of repeated benzimidazole treatments with increasing dosages on the phenotype of resistance and the beta-tubulin codon 200 genotype distribution in a benzimidazole resistant cyathostomin population. Vet Parasitol. 2004, 123: 201-213. 10.1016/j.vetpar.2004.06.019.

    Article  CAS  PubMed  Google Scholar 

  51. Kotze AC, Cowling K, Bagnall NH, Hines BM, Ruffell AP, Hunt PW: Relative level of thiabendazole resistance associated with the E198A and F200Y SNPs in larvae of a multi-drug resistant isolate of Haemonchus contortus. Int J Parasitol Drugs Drug Resist. 2012, 2: 92-97. 10.1016/j.ijpddr.2012.02.003.

    Article  PubMed Central  PubMed  Google Scholar 

  52. Cudekova P, Varady M, Dolinska M, Konigova A: Phenotypic and genotypic characterisation of benzimidazole susceptible and resistant isolates of Haemonchus contortus. Vet Parasitol. 2010, 172: 155-159. 10.1016/j.vetpar.2010.04.022.

    Article  CAS  PubMed  Google Scholar 

  53. Stenhouse LJ: Characterisation of Anthelmintic Resistance in a Multiple Drug Resistant Teladorsagia circumcincta Isolate. 2007, University of Glasgow, Glasgow

    Google Scholar 

  54. Sargison ND: Pharmaceutical treatments of gastrointestinal nematode infections of sheep–future of anthelmintic drugs. Vet Parasitol. 2012, 189: 79-84. 10.1016/j.vetpar.2012.03.035.

    Article  CAS  PubMed  Google Scholar 

  55. Lanusse CE, Prichard RK: Relationship between pharmacological properties and clinical efficacy of ruminant anthelmintics. Vet Parasitol. 1993, 49: 123-158. 10.1016/0304-4017(93)90115-4.

    Article  CAS  PubMed  Google Scholar 

  56. Taylor MA, Coop RL, Wall RL: Veterinary Parasitology. 2007, Blackwell Publishers, Oxford, 3

    Google Scholar 

  57. Coles GC, Roush RT: Slowing the spread of anthelmintic resistant nematodes of sheep and goats in the United Kingdom. Vet Rec. 1992, 130: 505-510. 10.1136/vr.130.23.505.

    Article  CAS  PubMed  Google Scholar 

  58. Roush RT, McKenzie JA: Ecological genetics of insecticide and acaricide resistance. Annu Rev Entomol. 1987, 32: 361-380. 10.1146/annurev.en.32.010187.002045.

    Article  CAS  PubMed  Google Scholar 

  59. Waghorn TS, Miller CM, Oliver AM, Leathwick DM: Drench and shift is a high-risk practice in the absence of refugia. N Z Vet J. 2009, 57: 359-363. 10.1080/00480169.2009.64723.

    Article  CAS  PubMed  Google Scholar 

  60. Martin PJ, Le Jambre LF, Claxton JH: The impact of refugia on the development of thiabendazole resistance in Haemonchus contortus. Int J Parasitol. 1981, 11: 35-41. 10.1016/0020-7519(81)90023-0.

    Article  CAS  PubMed  Google Scholar 

  61. Martin PJ: Selection for thiabendazole resistance in Ostertagia spp. by low efficiency anthelmintic treatment. Int J Parasitol. 1989, 19: 317-325. 10.1016/0020-7519(89)90143-4.

    Article  CAS  PubMed  Google Scholar 

  62. Van Wyk JA: Refugia - overlooked as perhaps the most potent factor concerning the development of anthelmintic resistance. Onderstepoort J Vet Res. 2001, 68: 55-67.

    CAS  PubMed  Google Scholar 

  63. van Dijk J, Morgan ER: The influence of water and humidity on the hatching of Nematodirus battus eggs. J Helminthol. 2012, 86: 287-292. 10.1017/S0022149X1100040X.

    Article  CAS  PubMed  Google Scholar 

  64. van Dijk J, Morgan ER: Variation in the hatching behaviour of Nematodirus battus: polymorphic bet hedging?. Int J Parasitol. 2010, 40: 675-681. 10.1016/j.ijpara.2009.11.002.

    Article  CAS  PubMed  Google Scholar 

  65. van Dijk J, Morgan ER: The influence of temperature on the development, hatching and survival of Nematodirus battus larvae. Parasitology. 2008, 135: 269-283.

    Article  CAS  PubMed  Google Scholar 

  66. Barnett C, Hossell J, Perry M, Procter C, Hughs G: A handbook of Climate Trends Across Scotland. 2006, Scotland & Northern Ireland Forum for Environmental Research, Scotland, Northern Ireland

    Google Scholar 

  67. Scott EW, Baxter P, Armour J: Fecundity of anthelmintic resistant adult Haemonchus contortus after exposure to ivermectin or benzimidazoles in vivo. Res Vet Sci. 1991, 50: 247-249. 10.1016/0034-5288(91)90117-7.

    Article  CAS  PubMed  Google Scholar 

  68. Echevarria FA, Armour J, Bairden K, Duncan JL: Laboratory selection for ivermectin resistance in Haemonchus contortus. Vet Parasitol. 1993, 49: 265-270. 10.1016/0304-4017(93)90125-7.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We gratefully acknowledge The Scottish Government’s Rural and Environment Science and Analytical Services Division (RESAS) for funding this work and our colleagues in the Bioservices Division, Moredun Research Institute, for expert care and assistance with animals. Scanning surveillance as carried out by AHVLA is funded by the Department for Environment, Food & Rural Affairs (Defra).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alison A Morrison.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SM, RM and IR were instrumental in initial parasite detection and collection of material for analysis. DJB and AAM participated in the design and undertaking of the controlled efficacy test, performed the statistical analysis and drafted the manuscript. AAM and JBM conducted the molecular genetic studies and participated in the sequence alignment. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Morrison, A.A., Mitchell, S., Mearns, R. et al. Phenotypic and genotypic analysis of benzimidazole resistance in the ovine parasite Nematodirus battus. Vet Res 45, 116 (2014). https://doi.org/10.1186/s13567-014-0116-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13567-014-0116-5

Keywords