Prevalence of anthelmintic resistance of gastrointestinal nematodes in Polish goat herds assessed by the larval development test
BMC Veterinary Research volume 17, Article number: 19 (2021)
Helminthic infections, in particular those caused by gastrointestinal nematodes (GIN), are found worldwide and are among the most economically important diseases of goats. Anthelmintic resistance (AR) in GIN of goats is currently present worldwide, and single- or multidrug resistant species are widespread. The aim of this study was to determine the prevalence of AR to benzimidazoles (BZ), macrocyclic lactones (ML) and imidazothiazoles represented by levamisole (LEV) in the Polish goat herds by using an in vitro larval development test, which is useful especially in large-scale epidemiological surveys.
This cross-sectional study was conducted from September 2018 to June 2019 and enrolled 42 dairy goat herds scattered over the entire country. The most commonly used anthelmintic class in goat herds in Poland were BZ (92%), followed by ML (85%) and LEV (13%). BZ-resistant GIN populations were found in 37 herds (88%, CI 95%: 75 to 95%), ML-resistant GIN populations in 40 herds (95%, CI 95, 84 to 99%), and LEV-resistant GIN populations in 5 herds (12%, CI 95%: 5 to 25%). Multidrug resistance involving all three anthelmintic classes was found in 5 herds (12%, CI 95, 5 to 25%). Based on the morphological features of stage 3 larvae the main resistant GIN turned out to be Haemonchus contortus and Trichostrongylus spp. The use of BZ and frequency of anthelmintic treatments were significantly related to the presence of AR to BZ in Polish goat herds.
This cross-sectional study demonstrates the existence of AR to BZ, ML and LEV on Polish goat farms. Resistance to BZ and ML is widespread, while AR to LEV is currently at a low level. A considerable proportion of herds harbours multidrug resistant GIN, which requires further consideration. An effective anthelmintic treatment strategy, reasonable preventive measures and better understanding of the resistance-related management practices by farmers and veterinarians may delay further development of AR.
Parasitic infections, especially those caused by gastrointestinal nematodes (GIN), are one of the main factors responsible for economic losses in goat farming around the world . Their control is mainly based on the use of three chemical classes of anthelmintics: benzimidazoles (BZ), macrocyclic lactones (ML), and imidazothiazoles including levamisole (LEV). The widespread and uncontrolled use of anthelmintics has resulted in the emergence of anthelmintic resistance (AR). In some countries the proportion of resistant GIN strains is currently so high that precludes effective control of parasitic diseases .
Several factors are responsible for the development of AR in parasites, of which most important are a high treatment frequency , underdosing of the anthelmintics, and continuous use of the same anthelmintic class over several years . These factors, together with certain types of farm management, can promote the development of AR, especially in goats which need higher doses to ensure anthelmintic efficacy since they metabolise and eliminate various medicines quicker than sheep and cattle [5, 6].
Several in vivo and in vitro tests have been developed for detection of AR. An in vivo fecal egg count reduction test (FECRT) is recommended by the World Association for the Advancement of Veterinary Parasitology (WAAVP). However, this test requires animals in a herd to be either tested twice which makes it expensive, time-consuming and laborious, or a randomly selected group of animals in a herd be left untreated which, in turn, is impractical and hardly acceptable to farmers. Moreover, high inter-animal variation in the pharmacokinetics of anthelmintics in goats may lower the quality of FECRT results [5, 6]. The use of some in vitro methods like egg hatch test (EHT) or molecular tests (RT-PCR or pyrosequencing) is currently limited to BZ . Moreover, EHT can be performed only on fresh fecal samples containing only eggs in early stage of development, while molecular tests are expensive and require specialized equipment in the laboratory, which makes them unsuited for routine AR diagnostics. Therefore, the most efficient in vitro test is the larval development test (LDT) which currently exists in several modifications allowing detection of AR to all three anthelmintic classes [7,8,9,10]. The LDT offers an alternative to the laborious in vivo FECRT and allows investigation of AR to all anthelmintic classes in a single test regardless of the herd size . Moreover, the LDT is the only in vitro AR diagnostic test that has been commercialized and registered on the market as DrenchRite® .
Thus far, a number of reports have been published on the occurrence of AR in goat herds in Europe, and the prevalence of AR especially to BZ appears to be very high in some European countries (Rose et al. 2015). In France studies have shown a prevalence of AR to BZ to vary between 70 to 100% [13, 14]. Very close are results of a recent Slovakian study . Many studies have demonstrated the existence of AR to the main anthelmintic classes in goat herds outside of Europe, namely in the United States [11, 16], Cuba , Kenya , South Africa , Uganda , Ethiopia , Malaysia [22, 23], India [24,25,26], and Pakistan .
Data on the AR situation in Polish small ruminant population are only fragmentary. In the last two decades, sporadic cases of AR to BZ in sheep, cattle, horses, pigs and goats have been reported [28,29,30,31]. Recently, first cases of AR to ML and LEV, as well as multidrug resistance (MDR) in goat herds have been described [32, 33].
On the other hand, our long cooperation with a number of Polish goat farmers has shown not only that GIN infections constitute a considerable clinical problem but also that the factors predisposing to the development of AR are present in the vast majority of goat herds. Therefore, we hypothesized that those cases of AR reported so far were only the tip of the iceberg, and we decided to carry out a large-scale epidemiological observational study to determine the prevalence of AR to the three basic anthelmintic classes in Polish goat herds and identify the genus or species of resistant larvae.
Goat herd characteristics
Forty-two dairy goat herds were enrolled in this cross-sectional study (49% of 85 herds invited). They were evenly scattered over the entire country (located in 13 of 16 provinces of Poland). They counted from 4 to 155 adult goats with the median (IQR) of 15 (10 to 23) goats. Only 5 herds (12%) consisted of more than 50 adult goats and only 2 of them counted more than 100 heads. In most of them two Polish local breeds, Polish White Improved and Polish Fawn Improved, were kept. Nine herds (21%) kept only Anglo-Nubian goats and three herds (7%) kept a traditional local goat breed – Carpathian. In the vast majority of herds (39 of 42; 93%) goats were grazed from April to October. In 10 herds (24%) goats were grazed together on the same pasture with sheep and in 8 herds (19%) with cattle. Only 13 herds (31%) relied solely on their own replacement, the remaining 29 herds (69%) purchased goats from other herds.
Only 3 herds (7%) did not practice routine deworming. In the remaining 39 herds (93%) at least one anthelmintic was routinely used. Most of the herds which practiced routine deworming (n = 26; 67%) dewormed twice a year, 7 herds (18%) once a year, and 6 herds (15%) three or four times a year. BZ (albendazole and fenbendazole) were used in 36 herds (92% of herds practicing routine deworming), ML (eprinomectin, rarely ivermectin) in 33 herds (85%) and LEV in only 5 herds (13%).
Prevalence of anthelmintic resistance
AR to each of the anthelmintic classes was detected using LDT in Polish goat herds. The most prevalent was AR to ML detected in 40 herds (95%, CI 95%: 84 to 99%) and to BZ detected in 37 herds (88%, CI 95%: 75 to 95%). AR to LEV was detected in 5 herds (12%, CI 95%: 5 to 25%) (Table 1).
Only one herd was free from AR to any of anthelmintics. AR to one anthelmintic was found in 5 herds (AR to ML in 4 herds, and to BZ in one herd), AR to two anthelmintics in 31 herds (BZ and ML in all) and AR to all three anthelmintics (synonymous to MDR) in 5 herds, which yielded the prevalence of MDR to be 12% (CI 95%: 5 to 25%). In most of herds AR to BZ and ML was severe (> 50% of larvae developed at the discriminating concentration, DC), while AR to LEV was mostly mild (< 10% of larvae developed at the DC) (Fig. 1).
Corrected percentage of larvae developing in tested wells (cPD) at the DC was significantly positively correlated with EC50 for all three anthelmintic classes and this correlation was high for BZ (Rs = 0.94, p < 0.001) and ML (Rs = 0.90, p < 0.001), and moderate for LEV (Rs = 0.48, p = 0.001) (Fig. S1). EC50 and EC99 for susceptible and resistant GIN populations are summarized in Table 2.
cPD at the DC in BZ-resistant populations ranged from 6 to 99% with the median (IQR) of 63% (32 to 91%), and was not significantly different from cPD at the DC in ML-resistant populations which ranged from 5 to 100% with the median (IQR) of 71% (52 to 92%) (p = 0.883). cPD at the DC in LEV-resistant populations ranged from 4 to 54% with the median of 9%, and was significantly lower than both the former (p = 0.047) and the latter GIN populations (p = 0.010).
Gastrointestinal nematode populations
H. contortus was present in 41 herds (98%, CI 95%: 88 to 100%), Trichostrongylus spp. in 37 herds (88%, CI 95%: 75 to 95%), Oesophagostomum spp. in 22 herds (52%, CI 95%: 38 to 67%), and Teladorsagia spp. in 14 herds (33%, CI 95%: 21 to 48%). Median fecal egg count (FEC) in the herds ranged from 75 to 2450 epg.
At the DC, no Oesophagostomum larvae developed irrespective of the herd’s AR status. In the case of BZ-resistance H. contortus larvae developed at the DC in 95% of the populations which initially harboured this GIN, Trichostrongylus larvae in 81% of the populations, and Teladorsagia larvae in 36% of the populations. In the case of ML-resistance Trichostrongylus larvae developed at the DC in 91% of the populations, Teladorsagia larvae in 79% of the populations and H. contortus larvae in 74% of the populations. In the case of LEV-resistance Trichostrongylus larvae developed at the DC in 100% of the populations, while H. contortus larvae in 20% of the populations, and Teladorsagia larvae in no population. Of 5 herds resistant to LEV, 3 had Trichostrongylus spp. as the only GIN left in wells with LEV at the DC, and 2 had Trichostrongylus spp. and H. contortus. These results indicated that H. contortus spp. and Trichostrongylus spp. were mainly responsible for AR to BZ, all three GIN were responsible for AR to ML, while Trichostrongylus spp. was the main GIN resistant to LEV (Fig. 2).
Factors linked to AR
Presence of AR to BZ and ML in virtually all investigated herds reduced the reliability of most statistical analyses. The number of times goats had been dewormed was significantly positively linked with the intensity of AR to BZ (p = 0.007), but not to ML (p = 0.643) or LEV (p = 0.849). The use of BZ was significantly related to the presence of AR to BZ (p = 0.002), however there was no significant relationship between the use of ML and LEV and the presence of relevant AR (p = 0.449 and p = 0.099, respectively). Interestingly, AR to LEV was observed in 3 herds whose owners declared not to have used LEV in previous 3 years. The same applied to 9 herds with AR to ML, and 3 herds with AR to BZ. None of the following: access to pasture, introducing goats from other herds, and grazing with sheep or cattle on common pasture was significantly linked with AR to any of the anthelmintics.
Our cross-sectional study shows that AR to BZ and ML is widespread in Polish goat herds, whereas AR to LEV is uncommon. This is the first study systematically investigating the prevalence of AR to three basic classes of anthelmintics in small ruminants in Poland. Similar surveys based on in vitro tests have so far been conducted in Europe on sheep flocks in Lithuania [9, 34] and Ireland , and sheep and goat flocks in Slovakia [15, 36]. Moreover, the LDT is a commonly applied technique in the AR prevalence studies in the United States [11, 16, 37, 38], Philippines [39, 40], Canada  and the New Zealand [42, 43] in both goat and sheep flocks. The LDT allows to estimate efficacy of anthelmintics in nematode populations from small farms where performance of a FECRT for even one anthelmintic would be impractical [11, 40]. This is the case in Poland, where most of goat herds consist of only a few animals.
Anthelmintics most often used in goats in Poland are albendazole, fenbendazole and eprinomectin, while levamisole is used rarely and this fact appears to directly correspond to the prevalence of AR. BZ are relatively inexpensive and have short withdrawal periods, which matters a lot to Polish farmers as their herds mainly produce milk and cheese. Moreover, BZ are believed to eliminate tapeworms and liver flukes which are commonly considered by farmers widespread and dangerous for goats, even though our observations seem to contradict those statements . Anyway, several decades of relying virtually solely on BZ must have led to the point at which we are now – the vast majority of herds harbour resistant GIN. Cold comfort is the fact that Polish situation does not seem to be much distinct from other European countries where prevalence of AR to BZ in goats is equally high [3, 15, 45].
Surprisingly, AR to ML is equally widespread in Poland, yet the first case was described only a year ago . The prevalence of AR to ML in goats varies between European countries from 0% in Norway , 20% in Italy  to 100% in some regions of Switzerland . Situation in Poland is closer to this in Switzerland and AR to ML in goats has already exceeded 90%. The most important cause seems to be the presence on the Polish veterinary market of eprinomectin in pour-on formulation for cattle with milk withdrawal period of 0 days (routinely extrapolated to goats). Pour-on formulation makes its use very simple and eprinomectin is also commonly used in Poland to control external parasite invasions without taking into account its simultaneous effect on GIN. These facts make it one of the most popular anthelmintic agents in Poland. However, eprinomectin is very likely to be commonly underdosed in Polish goat population since the appropriate dose of eprinomectin against GIN in goats is at least two times higher than that for cattle [49, 50], Unfortunately, our experience shows that this fact is unknown or routinely neglected by farmers and veterinarians, also for financial reasons. As all avermectins have the same mechanism of action the cross-resistance with ivermectin is complete . Moreover, it has been suggested that ML may select single nucleotide polymorphism at BZ-resistance-associated codons of β-tubulin [52, 53]. As a consequence, AR to ML and BZ appear to be mutually related phenomena and frequent use of ML may stimulate the development of BZ resistance.
AR to LEV in Polish goat population is still uncommon and mild, which appears to result from its rare use. LEV is viewed by veterinarians as an old-fashioned drug with narrow margin of safety and low effectiveness so they prefer to avoid it. Data on the prevalence of AR to LEV in goats in Europe are limited, mainly case reports have been published in Denmark , France [13, 14], and the UK . Recent, large-scale Dutch study did not reveal AR to LEV in sheep flocks .
Interestingly, AR to certain anthelmintic classes was observed on some farms where their owners declared not to have used this particular anthelmintic class at all. We suppose there are three possible explanations to this phenomenon, none of which excludes the others. First, the farmers may have not been fully aware of anthelmintics used since some of them are administered by veterinarians during consultation on the farm, especially in the case of injectable anthelmintics. Secondly, our questionnaire applied to quite a long period of time (3 years) and it is likely that farmers could have simply forgotten the facts of using of anthelmintics on their farms. And thirdly, there could be some goats purchased from other herds in Poland or abroad which had been dewormed with this particular anthelmintic and thus carried resistant GIN population. It is mainly the case in herds with AR to LEV and without previous use of this anthelmintic as in all these herds Anglo-Nubian goats were kept and goats of this breed are usually imported to Poland from abroad, mostly from Germany where AR to all three anthelmintic classes in small ruminants has so far been detected . Moreover, owners of these three herds declared in the questionnaire to have purchased goats in the past.
Multidrug resistance of GIN in goats in Europe is not widespread , and only a few cases have so far been reported in Denmark , France , Switzerland , UK  and Poland . Given that on most of farms BZ and ML are used, either separately, simultaneously or interchangeably, multidrug AR will most likely develop in those herds which, seeking the solution for dropping effectiveness of routine deworming, will switch to levamisole and use it intensively for some time. As AR to LEV seems to develop quite quickly  a year of routine LEV use suffices to produce AR.
In our study we decided to use the DC rather than the EC50 criterion for differentiation between resistant and susceptible GIN as an application of the EC50 may sometimes underestimate the resistance and thus some herds might have been incorrectly classed by this criterion as susceptible . It has been suggested that using the EC99 values or the DC in the in vitro tests such as the LDT or EHT can substantially increase test sensitivity of the test and identify resistance when only a small proportion of the GIN population is resistant [8, 61, 62]. To maintain high accuracy of our classification we decided to use two criteria, one based on the corrected percentage of development at the DC and the another based on EC99, simultaneously.
AR to BZ and ML seems to be conveyed by three main GIN which we identified in our study – H. contortus, Trichostrongylus spp., and Teladorsagia spp. Only Oesophagostomum spp. showed no potential to become resistant to any anthelmintic which is at odds with some recent studies in which Oesophagostomum spp. was together with H. contortus and Trichostrongylus spp. resistant to BZ and ML [63, 64]. Moreover, our observation has to be treated with caution as the only DC which has been scientifically proven to indicate AR is the DC of IVM-AG for H. contortus. Initially, it had been set at 5.4 ng/ml, however then its elevation to 21.6 ng/ml was suggested , and as a consequence we decided to use the higher DC as a safer option. However, nothing is known about the DC indicating AR in Trichostrongylus spp. or Teladorsagia sp. so we cannot be sure that the development of these GIN at this DC corresponds to their resistance to ML. Even though it is very likely to be so, it undoubtedly warrants further investigation. Interestingly, it seems that Trichostrongylus spp., not H. contortus as we initially expected, is the main culprit responsible for the early stage of development of AR to LEV. We made this observation first in our recent study  and this survey appears to confirm it, however we do not know any studies regarding this topic in goats. Studies so far conducted in sheep and goats have indicated that both Trichostrongylus spp. and H. contortus are responsible for AR to LEV [11, 13, 14, 16, 25, 37, 38, 54, 57, 63,64,65,66,67,68,69].
Our work presents the results of an in vitro test (LDT) which may in some circumstances differ from those obtained in an in vivo method such as FECRT. Such discrepancies have multifactorial background and most of these factors remain beyond researcher’s control [15, 70, 71]. Nevertheless, LDT and FECRT have been proven to show moderate to good agreement, especially regarding BZ and LEV resistance detection [10, 15, 72, 73]. Therefore, our results are likely to present the true AR status of the examined goat herds.
The main drawback to our study is the fact that herds enrolled in the survey were not randomly selected. Their owners had sent fecal samples to our laboratory at least once before so they were probably somehow aware of the importance of parasite control in goats. Given that most of them used to deworm all their goats routinely, usually using the same drug for a long time, they were more likely to promote the development of AR on their farms. Therefore, it is possible that our results may to some extent overestimate AR prevalence in the country, and we hope they do so. Anyway, they allow to realize that this problem already exists in Poland and unless actions are undertaken to stop it, soon no anthelmintic will remain effective in some goat herds. The another drawback to our study is that when we performed LDT with LEV we wrongly used 1.3% DMSO in control wells while LEV was diluted in water (contrary to TBZ and IVM-AG which were diluted in 1.3% DMSO). As a result if 1.3% DMSO had a potential to reduce larval development, this mistake could have falsely increased the corrected percentage of developing larvae in tested wells at the DC. Our previous observations have not indicated that DMSO at such a low concentration impacts on larval development to any noteworthy extent (unpublished data). However, even if it did we hope that the double criterion we used for classification GIN populations as resistant would reduce the risk of falsely positive result to an acceptable minimum.
This study provides the first comprehensive data on the prevalence of AR in goat herds in Poland. AR to BZ and ML is widespread, while AR to LEV is currently at a low level. Therefore, an appropriate strategy of GIN control should be applied by farmers and veterinarians as soon as possible not to let the situation slip out of control.
This was a cross-sectional study and lasted virtually one year, from September 2018 to June 2019. During this time the owners of 85 private dairy goat herds who had cooperated with our laboratory in previous 4 years (2015–2018), were invited to voluntarily enrol in the study. Those who declared participation were asked to collect fecal samples after at least 8 weeks since the last deworming had passed, and detailed instructions on proper collection and preparation of anaerobically stored fecal samples were sent to them. In herds which consisted of less than 15 adult goats fecal samples were collected from all animals, otherwise samples were collected from 15 to 20 goats selected by the owners. The owners were asked to pick only these goats which had been present on their farm for at least 8 previous weeks to avoid including goats dewormed within this period of time by the former owner. Fecal samples were collected by the farmers (owners of the goats) directly from the rectum. Pooled fecal samples weighing 50–100 g were placed in plastic containers filled with tap water to ensure anaerobic conditions. Samples were delivered to the laboratory within 24 h at room temperature and processed within next 24 h so that the total time which had elapsed between sample collection and examination did not exceed 48 h.
In each herd the following information was collected using the questionnaire: the number of adult goats (> 6 month-old) and kids, access to the pasture, co-grazing or rotational grazing with sheep and cattle, purchase of goats from other herds, and the history of deworming including the number of times the herd had been dewormed in previous 3 years and drugs used.
According to Polish legal regulations (The Act of the Polish Parliament of 15 January 2015 on the Protection of Animals Used for Scientific or Educational Purposes, Journal of Laws 2015, item 266) no formal ethics consent was required for this study except for the informed consent of participants, which we obtained in written from each participating farmer.
Larval development test
Based on the technique defined by Hubert and Kerboeuf , the larval development test (LDT) was performed with the modifications of Várady et al. [75, 76]. Nematode eggs were collected by sequentially sieving the anaerobically stored fecal samples through stacked sieves with apertures of 250, 100 and 25 μm. The material collected on the 25 μm sieve was washed with tap water, sedimented, and then flotation method with saturated sodium chloride was applied in order to obtain the suspension of the nematodes eggs . After extraction, eggs were inspected microscopically to ensure that embryonation had not yet begun and suspended in deionized water at a concentration of 70–100 eggs in 10 μl per well. The pure thiabendazole (Sigma-Aldrich, Merck, Germany; TBZ), ivermectin aglycone (Santa Cruz Biotechnology, USA; IVM-AG), and levamisole (Sigma-Aldrich, Merck, Germany; LEV) were firstly dissolved in ≥99.5% dimethyl sulfoxide (Sigma-Aldrich, Merck, Germany; DMSO) in case of TBZ and IVM-AG or deionized water in case of LEV and serially diluted 1:2 in DMSO (TBZ, IVM-AG) or in deionized water (LEV). The overall DMSO concentration in the pre-dilution plate was 20%. Then, diluted anthelmintics and DMSO (control) were moved to the test plate to yield final concentration of 1.3% in the tested wells. The 12 concentrations of TBZ, IVM-AG and LEV finally used in the LDT ranged from 0.0006 to 1.28 μg/ml, from 0.084 to 173.6 ng/ml, and from 0.020 to 32 μg/ml, respectively. Test was performed on 96-well cell culture plates (Sarsted, Germany) with culture medium (150 μl) which consisted of 10 μl of (all in one test plate) TBZ, IVM-AG, LEV or DMSO (1.3%; control wells) solution, 110 μl of deionised water, 20 μl of culture medium as described by Hubert and Kerboeuf  and 10 μl of a suspension (approximately 70–100 eggs) containing Amphotericin B (Sigma-Aldrich, Merck, Germany) at a concentration of 5 μg/ml. Each anthelmintic concentration was tested in duplicate. The LDT plates (sealed – to prevent drying) were incubated for 7 days at 25 °C (Cooled incubator, INCU-Line® Standard, VWR). After the incubation period, 10 μl of Lugol’s solution was added to each well to terminate development of larvae. The unhatched eggs and L1-L3 larvae were counted under an inverted microscope in each well (Olympus, CKX53, Poland) and the L3 larvae were classified at species/genus level in the tested and control wells following the procedure detailed elsewhere . The arithmetic mean of the percentage of larvae developing in two tested wells at each anthelmintic concentration (percentage development in tested wells, PDT) was corrected by the percentage of developing larvae in the control wells (PDC) according to the following formula:
cPD = PDT / PDC,
where cPD stood for corrected percentage of larvae developing in tested wells.
The concentration of each anthelmintic inhibiting development of 50% (median effective concentration, EC50) and 99% (EC99) of larvae was estimated using the 4 parameter logistic curve .
The results of LDT were interpreted with respect to the discriminating concentration (DC) of anthelmintic agent which was defined as the concentration of anthelmintic at which the development of at least 99% of susceptible larvae (corrected by the larvae developing in control wells as mentioned above) would have been inhibited [7, 62]. The following DC were used in the study: 0.08 μg/ml for TBZ , 21.6 ng/ml [34, 36] for IVM-AG, and 2.0 μg/ml for LEV .
Two criteria had to be simultaneously fulfilled to classify a GIN population from a given goat herd as resistant to the particular anthelmintic: cPD at the DC significantly higher than 1% (meaning that the entire 95% confidence interval (CI 95%) was above 1%) and EC99 > DC.
AR was arbitrarily classed in terms of severity as follows: cPD at the DC from > 1 to 10% – mild AR, from > 10 to 50% – moderate AR, above 50% – severe AR. Multidrug AR (MDR) was defined as the presence of AR to all three different classes of anthelmintics.
Numerical variables were presented as the median, interquartile range (IQR) and range, and compared between groups using the Mann-Whitney U test (unpaired groups) or Wilcoxon signed rank test (paired groups). Categorical variables were given as count and percentage in a group, and compared between groups using the Pearson’s chi-square test or the Fisher exact test. The 95% confidence intervals (CI 95%) for percentages were calculated using the Wilson score method . Correlations between cPD at the DC and EC50 were assessed using the Spearman’s rank correlation coefficient (Rs). The significance level (α) was set at 0.05 and the Bonferroni correction was applied in the case of multiple comparisons. Statistical analysis was performed in TIBCO Statistica 13.3.0 (TIBCO Software Inc., Palo Alto, CA).
Availability of data and materials
The data sets used and/or analyzed are available from the corresponding author on reasonable request.
- CI 95%:
Confidence interval for the level of confidence of 95%
Corrected percentage of developing larvae
- EC50 :
Median effective concentration
- EC99 :
The 99th percentile effective concentration
Egg hatch test
Eggs per gram
Fecal egg count
Fecal egg count reduction test
Larval development test
- PDC :
Percentage development in control wells
- PDT :
Percentage development in tested wells
- Rs :
Spearman’s rank correlation coefficient
Hoste H, Sotiraki S, de Jesús Torres-Acosta JF. Control of endoparasitic nematode infections in goats. Vet Clin North Am Food Anim Pract. 2011;27:163–73.
Várady M, Papadopoulos E, Dolinská M, Königová A. Anthelmintic resistance in parasites of small ruminants: sheep versus goats. Helmintologia. 2011;48:137–44.
Chartier C, Pors I, Hubert J, Rocheteau D, Benoit C, Bernard N. Prevalence of anthelmintic resistant nematodes in sheep and goats in Western France. Small Rum Res. 1998;29:33–41.
Sargison ND, Jackson F, Bartley DJ, Wilson DJ, Stenhouse LJ, Penny CD. Observations on the emergence of multiple anthelmintic resistance in sheep flocks in the south-east of Scotland. Vet Parasitol. 2007;145:65–76.
Aksit D, Yalinkilinc HS, Sekkin S, Boyacioğlu M, Cirak VY, Ayaz E, Gokbulut C. Comparative pharmacokinetics and bioavailability of albendazole sulfoxide in sheep and goats, and dose dependent plasma disposition in goats. BMC Vet Res. 2015;11:124.
González Canga A, Sahagún Prieto AM, José Diez Liébana M, Martínez NF, Vega MS, Vieitez JJ. The pharmacokinetics and metabolism of ivermectin in domestic animal species. Vet J. 2009;179:25–37.
Coles GC, Jackson F, Pomroy WE, Prichard RK, von Samson-Himmelstjerna G, Silvestre A, Taylor MA, Vercruysse. The detection of anthelmintic resistance in nematodes of veterinary importance. Vet Parasitol. 2006;136:167–85.
Dolinská M, Königová A, Letková V, Molnár L, Várady M. Detection of ivermectin resistance by a larval development test – Back to the past or step forward? Vet Parasitol. 2013;198:154–8.
Kupčinskas T, Stadalienė I, Šarkūnas M, Petkevičius S. Anthelmintic resistance in sheep farms in Lithuania detected by in vitro micro-agar larval development test. Vet Med Zoot. 2015;72:21–3.
Ruffell A, Raza A, Elliott TP, Kotze AC. The use of the larval development assay for predicting the in vivo efficacy of levamisole against Haemonchus contortus and Trichostrongylus colubriformis. Vet Parasitol. 2018;30:6–11.
Crook EK, O'Brien DJ, Howell SB, Storey BE, Whitley NC, Burke JM, Kaplan RM. Prevalence of anthelmintic resistance on sheep and goat farms in the mid-Atlantic region and comparison of in vivo and in vitro detection methods. Small Rumin Res. 2016;143:89–96.
Anonymous. Drenchrite®©, Larval Development Assay Standard Operating Procedures, Horizon Technology Pty Ltd, Roseville, NSW, Australia, 1996.
Chartier C, Soubirac F, Pors I, Silvestre A, Hubert J, Couquet C, Cabaret J. Prevalence of anthelmintic resistance in gastrointestinal nematodes of dairy goats under extensive management conditions in southwestern France. J Helminthol. 2001;75:325–30.
Paraud C, Kulo A, Pors I, Chartier C. Resistance of goat nematodes to multiple anthelmintics on a farm in France. Vet Rec. 2009;164:563–4.
Babják M, Konigöva A, Urda-Dolinská M, Vadlejch J, Varády M. Anthelmintic resistance in goat herds – in vivo versus in vitro detection methods. Vet Parasitol. 2018;254:10–4.
Howell SB, Burke JM, Miller JE, Terrill TH, Valencia E, Williams MJ, Williamson LH, Zajac AM, Kaplan RM. Prevalence of anthelmintic resistance on sheep and goat farms in the south-eastern United States. J Am Vet Med Assoc. 2008;233:1913–9.
Arece-García J, López-Leyva Y, Olmedo-Juárez A, Ramírez-Vargas G, Reyes-Guerrero D, López Arellano ME, De Gives PM, Várady M, Rojo-Rubio R, González-Garduño R. First report of multiple anthelmintic resistance in goat farm in Cuba. Helminthologia. 2017;54:358–62.
Wanyangu SW, Bain RK, Rugutt MK, Nginyi JM, Mugambi JM. Anthelmintic resistance amongst sheep and goats in Kenya. Prev Vet Med. 1996;25:285–90.
Bakunzi FR, Nkomo LK, Motsei LE, Ndou RV, Nyirenda M. A survey on anthelmintic resistance in nematode parasites of communally grazed sheep and goats in a rural area of north West Province, Republic of South Africa. Life Sci J. 2013;10:391–3.
Nabukenya I, Rubaire-Akiiki C, Olila D, Muhangi D, Höglund J. Anthelmintic resistance in gastrointestinal nematodes in goats and evaluation of FAMACHA diagnostic marker in Uganda. Vet Parasitol. 2014;205:666–75.
Wakayo BU, Dewo TF. Anthelmintic resistance of gastrointestinal parasites in small ruminants: a review of the case of Ethiopia. J Vet Sci Technol. 2015;S10:001.
Dorny P, Claerebout E, Vercruysse J, Sani R, Jalila A. Anthelmintic resistance in goats in peninsular Malaysia. Vet Parasitol. 1994;55:327–42.
Basripuzi HB, Sani RA, Ariff OM. Anthelmintic resistance in selected goat farms in Kelantan. Mal J Anim Sci. 2012;15:47–56.
Das M, Singh S. Anthelmintic resistance to nematodes in sheep and goat farms in Hisar. J Vet Parasitol. 2005;19:103–6.
Manikkavasagan I, Binosundar ST, Raman M. Survey on anthelmintic resistance to gastrointestinal nematodes in unorganized goat farms of Tamil Nadu. J Parasit Dis. 2015;39:258–61.
Singh R, Bal MS, Singla LD, Kaur P. Detection of anthelmintic resistance in sheep and goat against fenbendazole by faecal egg count reduction test. J Parasit Dis. 2016;41:463–6.
Jabbar A, Iqbal Z, Saddiqi HA, Babar W, Saeed M. Prevalence of multiple anthelmintic resistant gastrointestinal nematodes in dairy goats in a desolated tract (Pakistan). Parasitol Res. 2008;103:29–35.
Balicka-Ramisz A, Ramisz A, Malecki J, Supera K. Case of the resistance of the gastrointestinal nematodes to fenbendazole in sheep (in polish with summary in English). Magazyn Weterynaryjny. 1997;6:442–3.
Balicka-Ramisz AK, Ramisz AZ. Benzimidazoles resistance in nematode parasites in domesticated animals in north-west part of Poland. Electron J Pol Agric Univ. 1999;2:2 http://www.ejpau.media.pl/volume2/issue2/animal/art-02.html.
Kowal J, Wyrobisz A, Nosal P, Kucharski M, Kaczor U, Skalska M, Sendor P. Benzimidazole resistance in the ovine Haemonchus contortus from southern Poland –coproscopical and molecular findings. Ann Parasitol. 2016;62:119–23.
Mickiewicz M, Czopowicz M, Górski P, Kaba J. The first reported case of resistance of gastrointestinal nematodes to benzimidazole anthelmintic in goats in Poland. Ann Parasitol. 2017;63:317–22.
Mickiewicz M, Czopowicz M, Moroz A, Szaluś-Joardanow O, Górski P, Várady M, Königová A, Spinu M, Lefkaditis M, Kaba J. Development of resistance to eprinomectin in gastrointestinal nematodes in a goat herd with pre-existing resistance to benzimidazoles. Pol J Vet Sci. 2019;22:753–60.
Mickiewicz M, Czopowicz M, Kawecka-Grochocka E, Moroz A, Szaluś-Jordanow O, Várady M, Königová A, Spinu M, Górski P, Bagnicka E, Kaba J. The first report of multidrug resistance in gastrointestinal nematodes in goat population in Poland. BMC Vet Res. 2020;16:270.
Kupčinskas T, Stadalienė I, Šarkūnas M, Riškevičienė V, Várady M, Höglund J, Petkevičius S. Prevalence of anthelmintic resistance on Lithuanian sheep farms assessed by in vitro methods. Acta Vet Scand. 2015;57:88.
Good B, Hanrahan JP, De Waal DT, Patten T, Kinsella A, Lynch CO. Anthelmintic-resistant nematodes in Irish commercial sheep flocks – the state of play. Ir Vet J. 2012;65:21.
Dolinská M, Ivanišinová O, Königová A, Várady M. Anthelmintic resistance in sheep gastrointestinal nematodes in Slovakia detected by in-vitro methods. BMC Vet Res. 2014;10:233.
Terrill TH, Kaplan RM, Larsen M, Samples OM, Miller JE, Gelaye S. Anthelmintic resistance on goat farms in Georgia: efficacy of anthelmintics against gastrointestinal nematodes in two selected goat herds. Vet Parasitol. 2001;97:261–8.
Mortensen LL, Williamson LH, Terrill TH, Kircher RA, Larsen M, Kaplan RM. Evaluation of prevalence and clinical implications of anthelmintic resistance in gastrointestinal nematodes in goats. J Am Vet Med Assoc. 2003;223:495–500.
Venturina VM, Ancheta PB, Dobson RJ, Eduardo SL, Gray GD. Use of a larval development assayto estimate anthelmintic efficacy in goats in smallholder farming systems. Philippine Agricultural Scientist. 2003;86:134–9.
Ancheta PB, Dumilon RA, Venturina VM, Cerbito WA, Dobson RJ, LeJambre LF, Villar EC, Gray GD. Efficacy of benzimidazole anthelmintics in goats and sheep in the Philippines using a larval development assay. Vet Parasitol. 2004;120:107–21.
Falzon LC, Menzies PI, Shakya KP, Jones-Bitton A, Vanleeuwen J, Avula J, Stewart H, Jansen JT, Taylor MA, Learmount J, Peregrine AS. Anthelmintic resistance in sheep flocks in Ontario, Canada. Vet Parasitol. 2013;193:150–62.
Shyam S. A survey of anthelmintic resistance and parasite management practices on sheep farms in New Zealand: a dissertation presented in partial fulfilment of the requirements for the degree of Master of Veterinary Studies (Veterinary Parasitology) at Massey University. Palmerston North: Master of Veterinary Studies (M.V.S.); 2004. p. 37–54.
Kamaludeen J. Studies on the occurrence of anthelmintic resistance in goat parasites in New Zealand: a thesis presented in partial fulfilment of the requirements for the degree of Master of Veterinary Studies in Parasitology at Massey University, Palmerston North, New Zealand. Palmerston North: Master of Veterinary Studies (M.V.S.); 2010. p. 52–73.
Mickiewicz M, Czopowicz M, Moroz A, Witkowski L, Szaluś-Jordanow O, Nalbert T, Markowska-Daniel I, Górski P, Kaba J. Common endoparasitic infections in goats in Poland – diagnostics and treatment (in polish with summary in English). Życie weterynaryjne. 2017;92:665–8.
Rose H, Rinaldi L, Bosco A, Mavrot F, de Waal T, Skuce P, Charlier J, Torgerson PR, Hertzberg H, Hendrickx G, Vercruysse J, Morgan ER. Widespread anthelmintic resistance in European farmed ruminants: a systematic review. Vet Rec. 2015;176:546.
Domke AVM, Chartier C, Gjerde B, Höglund J, Leine N, Vatn S, Stuen S. Prevalence of anthelmintic resistance in gastrointestinal nematodes of sheep and goats in Norway. Parasitol Res. 2012;111:185–93.
Zanzani S, Gazzonis A, Di Cerbo A, Várady M, Manfredi M. Gastrointestinal nematodes of dairy goats, anthelmintic resistance and practices of parasite control in northern Italy. BMC Vet Res. 2014;10:114.
Murri S, Knubben-Schweizer G, Torgerson P, Hertzberg H. Frequency of eprinomectin resistance in gastrointestinal nematodes of goats in canton Berne, Switzerland. Vet Parasitol. 2014;203:114–9.
Gawor J, Borecka A, Malczewski A. Use of eprinomectin (Eprinex pour-on) to control natural infection by gastro-intestinal nematodes in goats. Med Weter. 2000;56:398–400.
Rostang A, Devos J, Chartier C. Review of the Eprinomectin effective doses required for dairy goats: where do we go from here? Vet Parasitol. 2020;277:1089–92.
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.
Mottier ML, Prichard RK. Genetic analysis of a relationship between macrocyclic lactone and benzimidazole anthelmintic selection on Haemonchus contortus. Pharmacogen Genomics. 2008;18:129–40.
Williamson SM, Storey B, Howell S, Harper KM, Kaplan RM, Wolstenholme AJ. Candidate anthelmintic resistance-associated gene expression and sequence polymorphisms in a triple-resistant field isolate of Haemonchus contortus. Mol Biochem Parasitol. 2011;180:99–105.
Maingi N, Bjørn H, Thamsborg SM, Bøgh HO, Nansen P. A survey of anthelmintic resistance in nematode parasites of goats in Denmark. Vet Parasitol. 1996;66:53–66.
Hong C, Hunt KR, Coles GC. Occurrence of anthelmintic resistant nematodes on sheep farms in England and goat farms in England and Wales. Vet Rec. 1996;139:83–6.
Ploeger HW, Everts RR. Alarming levels of anthelmintic resistance against gastrointestinal nematodes in sheep in the Netherlands. Vet Parasitol. 2018;262:11–5.
Voigt K, Scheuerle M, Hamel D, Pfister K. High perinatal mortality associated with triple anthelmintic resistance in a German sheep flock. Tierarztl Prax Ausg G Grosstiere Nutztiere. 2012;40:107–11.
Schnyder M, Torgerson PR, Schönmann M, Kohler L, Hertzberg H. Multiple anthelmintic resistance in Haemonchus contortus isolated from south African Boer goats in Switzerland. Vet Parasitol. 2005;128:285–90.
Jackson F, Jackson E, Coop RL. Evidence of multiple anthelmintic resistance in a strain of Teladorsagia circumcincta (Ostertagia circumcincta) isolated from goats in Scotland. Res Vet Sci. 1992;53:371–4.
Zajac AM, Gipson TA. Multiple anthelmintic resistance in a goat herd. Vet Parasitol. 2000;87:163–72.
Várady M, Cudeková P, Corba J. In vitro detection of benzimidazole resistance in Haemonchus contortus: egg hatch test versus larval development test. Vet Parasitol. 2007;149:104–10.
Čudeková P, Várady M, Dolinská M, Königová A. Phenotypic and genotypic characterisation of benzimidazole susceptible and resistant isolates of Haemonchus contortus. Vet Parasitol. 2010;172:155–9.
Dey AR, Begum N, Anisuzzaman Alim MA, Alam MZ. Multiple anthelmintic resistance in gastrointestinal nematodes of small ruminants in Bangladesh. Parasitol Int. 2020;77:102–5.
Herrera-Manzanilla FA, Ojeda-Robertos NF, González-Garduño R, Cámara-Sarmiento R, Torres-Acosta JFJ. Gastrointestinal nematode populations with multiple anthelmintic resistance in sheep farms from the hot humid tropics of Mexico. Vet Parasitol Reg Stud Reports. 2017;9:29–33.
Chandrawathani P, Yusoff N, Wan LC, Ham A, Waller PJ. Total anthelmintic failure to control nematode parasites of small ruminants on government breeding farms in Sabah. East Malaysia Vet Res Commun. 2004;28:479–89.
Waghorn T, Knight J, Leathwick D. The distribution and anthelmintic resistance status of Trichostrongylus colubriformis, T. vitrinus and T. axeiin lambs in New Zealand. N Z Vet J. 2014;62:152–9.
Chaparro JJ, Villar D, Zapata JD, López S, Howell SB, López A, Storey BE. Multi-drug resistant Haemonchus contortus in a sheep flock in Antioquia, Colombia. Vet Parasitol Reg Stud Reports. 2017;10:29–34.
Hodgson B, Mulvaney CJ. Resistance to a triple-combination anthelmintic in Trichostrongylus spp. on a commercial sheep farm in New Zealand. N Z Vet J. 2017;65:1–5.
Oliveira PAD, Riet-Correa B, Estima-Silva P, Coelho ACB, Santos BLD, Costa MAP, Ruas JL, Schild AL. Multiple anthelmintic resistance in southern Brazil sheep flocks. Rev Bras Parasitol Vet. 2017;26:427–32.
Vidyaszankar AN, Kaplan RM, Chan S. Statistical approach to measure the efficacy of anthelmintic treatment on horse farms. Parasitology. 2007;134:2027–39.
Vidyaszankar AN, Hanlon BM, Kaplan RM. Statistical and biological considerations in evaluating drug efficacy in equine strongyle parasites using fecal egg count data. Vet Parasitol. 2012;185:45–56.
Taylor MA, Learmount J, Lunn E, Morgan C, Craig BH. Multiple resistance to anthelmintics in sheep nematodes and comparison of methods used for their detection. Small Rumin Res. 2009;86:67–70.
Maingi N, Bjorn H, Dangolla A. The relationship between faecal egg count reduction and the lethal dose 50% in the egg hatch assay and larval development assay. Vet Parasitol. 1998;77:133–45.
Hubert J, Kerboeuf D. A micro-larval development assay for the detection of anthelmintic resistance in sheep nematodes. Vet Rec. 1992;130:442–6.
Várady M, Bjørn H, Nansen P. In vitro characterization of anthelmintic susceptibility of field isolates of the pig nodular worm Oesophagostomum spp., susceptible or resistant to various anthelmintics. Int J Parasitol. 1996;26:733–40.
Várady M, Corba J, Letková V, Kovác G. Comparison of two versions of larval development test to detect anthelmintic resistance in Haemonchus contortus. Vet Parasitol. 2009;160:267–71.
Hubert J, Kerboeuf D. A new method for culture of larvae used in diagnosis of ruminant gastrointestinal strongylosis: comparison with fecal cultures. Can J Comp Med. 1984;48:63–71.
van Wyk JA, Mayhew E. Morphological identification of parasitic nematode infective larvae of small ruminants and cattle: a practical lab guide. Onderstepoort J Vet Res. 2013;80:1–14.
Dobson RJ, Griffiths DA, Donald AD, Waller PJ. A genetic model describing the evolution of levamisole resistance in Trichostrongylus colubriformis, a nematode parasite of sheep. IMA J Math Appl Med Biol. 1987;4:279–93.
Taylor MA. A larval development test for the detection of anthelmintic resistance in nematodes of sheep. Res Vet Sci. 1990;49:198–202.
Altman D, Machin D, Bryant T, Gardner M. Statistics with Confidence: Confidence Intervals and Statistical Guidelines. 2nd ed. Bristol: BMJ Books; 2000.
This article is based upon work from COST Action COMBAR CA16230, supported by COST (European Cooperation in Science and Technology).
The study was partly financed by the grant from the Ministry of Science and Higher Education of the Republic of Poland, decision no. 9506/E-385/R/2018. The funding body provided financial support for modernizing and equipping the parasitological laboratory in which parasitological tests were performed. The funding body had no direct role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
The study was carried out in accordance with the standards recommended by the EU Directive 2010/63/EU for animal experiments and Good Laboratory Practice and The Act of the Polish Parliament of 15 January 2015 on the Protection of Animals Used for Scientific or Educational Purposes (Journal of Laws 2015, item 266). According to Polish legal regulations (The Act of the Polish Parliament of 15 January 2015 on the Protection of Animals Used for Scientific or Educational Purposes, Journal of Laws 2015, item 266) no formal ethics consent was required for this study except for the informed consent of participants. Therefore, the written informed consent for participation in the study, including processing fecal samples and data collected in the questionnaire, was obtained by us from all farmers who decided to participate in the study.
Consent for publication
The authors declare that they have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Correlations (presented as the Spearman’s rank correlation coefficients, Rs) between the corrected percentage of larvae developing in tested wells (cPD) at the discriminating concentration (DC) of each anthelmintic agent and the median effective concentration (ED50) of each anthelmintic agent. Scatter plots presenting correlations between cPD at DC and ED50 for each of anthelmintic agents.
About this article
Cite this article
Mickiewicz, M., Czopowicz, M., Moroz, A. et al. Prevalence of anthelmintic resistance of gastrointestinal nematodes in Polish goat herds assessed by the larval development test. BMC Vet Res 17, 19 (2021). https://doi.org/10.1186/s12917-020-02721-9