# First report of eprinomectin-resistant isolates of Haemonchus contortus in 5 dairy sheep farms from the Pyrénées Atlantiques département in France

S Jouffroy, L Bordes, C Grisez, J F Sutra, T Cazajous, J Lafon, N Dumont, M Chastel, C Vial-Novella, D Achard, et al.

► To cite this version:

S Jouffroy, L Bordes, C Grisez, J F Sutra, T Cazajous, et al. First report of eprinomectin-resistant isolates of Haemonchus contortus in 5 dairy sheep farms from the Pyrénées Atlantiques département in France. Parasitology, 2023, 150 (4), pp.365-373. 10.1017/s0031182023000069. hal-04031790

HAL Id: hal-04031790

https://hal.inrae.fr/hal-04031790

Submitted on 16 Mar 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Research Article

**Cite this article:** Jouffroy S *et al* (2023). First report of eprinomectin-resistant isolates of *Haemonchus contortus* in 5 dairy sheep farms from the *Pyrénées Atlantiques* département in France. *Parasitology* **150**, 365–373. https://doi.org/10.1017/S0031182023000069

Received: 23 November 2022

Revised: 5 January 2023

Accepted: 9 January 2023

First published online: 20 January 2023

**Keywords:**

Dairy sheep; eprinomectin; *Haemonchus contortus*; resistance

**Author for correspondence:**

S. Jouffroy, E-mail: sophie.jouffroy@envt.fr

# First report of eprinomectin-resistant isolates of *Haemonchus contortus* in 5 dairy sheep farms from the *Pyrénées Atlantiques* département in France

S. Jouffroy<sup>1,2,3</sup>, L. Bordes<sup>2</sup>, C. Grisez<sup>2</sup>, J. F. Sutra<sup>1</sup>, T. Cazajous<sup>4</sup>, J. Lafon<sup>5</sup>, N. Dumont<sup>6</sup>, M. Chastel<sup>7</sup>, C. Vial-Novella<sup>8</sup>, D. Achard<sup>3</sup>, H. Karembe<sup>3</sup>, M. Devaux<sup>2</sup>, M. Abbadie<sup>2</sup>, C. Delmas<sup>2</sup>, A. Lespine<sup>1</sup> and P. Jacquiet<sup>2</sup>

**Abstract**

Infection of sheep by gastrointestinal nematodes (GIN) in pastoral systems such as those found in the South Western area of France, the *Pyrénées Atlantiques*, is one of the main reasons for economic loss and degradation of their welfare. In the present study, the efficacy of eprinomectin (EPN) was monitored on farms from this area following suspicion of lack of anthelmintic efficacy. Suspicions were raised by veterinarians, based on clinical signs ranging from milk and body condition loss, to anaemia, and mortality. Resistance was evaluated according to the World Association for the Advancement for Veterinary Parasitology (WAAVP) guidelines using fecal egg count reduction tests reinforced by individual analysis of drug concentration in the serum of all treated ewes by high-performance liquid chromatography (HPLC). EPN was administered by subcutaneous (SC) and topical (T) route according to manufacturer's requirements, as well as by the oral route (O) with the topical solution according to off-labelled practices in the field. For the first time in France, the presence of resistant isolates of *Haemonchus contortus* to EPN was observed in 5 dairy sheep farms. The HPLC dosages showed exposure of worms to concentrations compatible with anthelmintic activity for animals treated by the SC and O routes. By contrast, they showed under exposure to the drug of most individuals treated by the T route. EPN is the only null milk withdrawal anthelmintic molecule currently available. The presence of resistant isolates of the pathogenic *H. contortus* to EPN in this important dairy region requires an urgent change in grazing, and sometimes production, systems.

**Introduction**

The most south-western *département* of mainland France, the *Pyrénées Atlantiques*, is the country's second largest dairy sheep production area after the Roquefort perimeter. It is where the Protected Designation of Origin (PDO) European labelled cheese Ossau Iraty is made, using milk from 3 local breads: *Manech Tête Rousse, Manech Tête Noire* and *Basco Bearnais* sheep. Sheep breeding and grazing for cheese making is a vital part of the local culture and economy, so much so that it has been included in the requirements for the production of PDO Ossau Iraty: 'ewes should graze for at least 240 days per lactation period' (INAO, 2015). Locally, the climate is oceanic, with mild temperatures year long and is one of the most humid parts of mainland France, receiving about 1300-1600 mm of rain per year (Meteo France). In this setting, sheep are frequently infected with heavy loads of gastrointestinal nematodes (GIN) when grazing, and farmers have to deal with the challenge of trying to control the parasite load almost all year long. Parasite control has relied on benzimidazoles for several decades, but the use of this family of molecules diminished with the increasing appearance of resistant GIN strains (Geurden *et al.*, 2014; Rose Vineer *et al.*, 2020) and when the milk withdrawal period changed from zero to at least 4 days in 2014 (Zoetis France, 2009), using these molecules during the lactation period was no longer a financially sound option.

Consequently, the macrocyclic lactone (ML) eprinomectin (EPN) has become the main treatment option during lactation: it has a very low blood to milk partition (Imperiale and Lanusse, 2021), making it the only available molecule in France with a zero milk withdrawal period. First commercialized for cattle as a topical formulation in 1996, and later on as an injectable formulation in 2015 in France, it was not until 2016 and 2020 that EPN was approved for small ruminants, for the topical (HPRA, 2016) and injectable formula, respectively (HPRA, 2020). Before 2016, the topical formulation was administered to dairy sheep and goats off label. Ineffectiveness of the topical formula was rapidly reported in goats on the basis of fecal egg count reduction tests (FECRT) (Murri *et al.*, 2014), which prompted veterinarians

© The Author(s), 2023. Published by Cambridge University Press. This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.

<figure><img src="image_3.png" /></figure>

CAMBRIDGE
UNIVERSITY PRESS

<footnote n="1"><sup>1</sup>INTHERES, Université de Toulouse, INRAE, ENVT, 31027 Toulouse Cedex 3, France; <sup>2</sup>IHAP, Université de Toulouse, INRAE, ENVT, 31027 Toulouse Cedex 3, France; <sup>3</sup>CEVA Santé Animale, 33500 Libourne, France; <sup>4</sup>Selarl Vétérinaire du Piémont, 64800 Mirepeix, France; <sup>5</sup>Clinique Vétérinaire du Haut Béarn, 64400 Oloron Ste Marie, France;</footnote>

<footnote n="6"><sup>6</sup>Vétérinaires Garazi, 64220 St Jean le Vieux, France; <sup>7</sup>Clinique du Saison, 64470 Tardets-Sorholus, France and</footnote>

<footnote n="8"><sup>8</sup>Centre Départemental Elevage Ovin, 64130 Ordiarp, France</footnote>
to use EPN via other routes of administration that are known to be associated with higher overall exposure in plasma and tissues and potentially higher efficacy than topical administration (Lespine et al., 2012). Lack of efficacy of the topical formula has also been recently described in dairy sheep (Bouy et al., 2021; Bordes et al., 2022).

Routinely, anthelmintic treatments are administered to the whole lactating flock, usually at a fixed time of the year determined by habit, production stage and/or season of the year. Animals are treated against GIN infections 3-4 times a year on average, using mainly molecules from the ML family (Centre Départemental de l'Elevage Ovin, unpublished data). Together with inaccurate animal weight measurements (underdosing), high frequency of treatment has been proven to be one of the main drivers of anthelmintic resistance (AR) (Wolstenholme et al., 2004; Falzon et al., 2014; Sangster et al., 2018) and loss of efficacy of EPN was expected to happen in the *Pyrénées Atlantiques* sooner or later. From 2018, veterinarians first reported loss of efficacy of avermectins: benzimidazole/ivermectin multi-resistant isolates of *Haemonchus contortus* have been isolated from an ovine meat production farm in the *Hautes Pyrénées* (Cazajous et al., 2018) and benzimidazole/EPN multi-resistant isolates of the parasite were identified in a dairy goat herd in the *Pyrénées Atlantiques* (Bordes et al., 2020). The implementation of this resistance motivated the creation of a 3-year long project, ANTHERIN for ANTHelmintic Resistance in dairy sheep farms: survey and INnovative solutions. The results presented in this study are linked to this project.

Of the 3 main pathogenic species for sheep and goats, *H. contortus* is the most pathogenic and prolific (Arsenopoulos et al., 2021). It has probably spread across the globe thanks to commercial activities, and has been able to adapt to different climates (Sallé et al., 2019), yet its development remains conditioned by external temperatures and humidity (O'Connor et al., 2006; Arsenopoulos et al., 2021). *Haemonchus contortus* has also been capable of adapting to anthelmintic treatment, and to this day resistance to all major anthelmintic families have been described (Kotze and Prichard, 2016). Adult worms being blood-sucking parasites, infection of sheep by *H. contortus* causes a range of symptoms depending on host susceptibility and parasite load, from loss of milk production and body condition, to anaemia and death (Arsenopoulos et al., 2021).

This study reports for the first time EPN resistance of *H. contortus* in 5 dairy sheep farms in France, investigated between June 2020 and April 2021. In addition to anthelmintic efficacy measured by FECRT, concentrations of EPN were determined in sheep sera 2 and 5 days after treatment, to differentiate cases of loss of efficacy due to drug resistance from those linked to underexposure of GIN to EPN.

## Materials and methods

### Farm selection

Five farms were included in the study based on suspicion of lack of efficacy of EPN in lactating dairy ewes by the veterinary practitioner. These suspicions emerged in February (farm 5), April (farms 1 and 4), June (farm 2) 2020 and April 2021 (farm 3), following oral or injectable EPN treatment. All farmers observed clinical signs compatible with strongylosis that did not improve after EPN treatment in lactating animals. Of these, 3 flocks had symptoms suggestive of haemonchosis (anaemia and on 2 farms mortality), and the remaining 2 flocks showed milk and weight losses. In all cases, the attending veterinarian did a fecal egg count (FEC) about 2 weeks after treatment that revealed the presence of strongyle eggs. Further investigation into the lack of efficacy was conducted to determine whether it was due to

underexposure of the strongyles to EPN, or due to the presence of a resistant strain of worms. The 5 investigated farms had an average of 360 lactating ewes [215-500]. Four out of the 5 farms worked with Basco-Bearnaise sheep, and 1 with Manech Tête Rousse (farm 3). All 5 farms sent their lactating ewes in collective middle (1000 m) to high-altitude (≥2000 m) summer pastures (at least 1 other sheep herd grazing in the same area).

### On-farm protocol

Efficacy of EPN was evaluated using FECRT according to the World Association for the Advancement for Veterinary Parasitology (WAAVP) Guidelines (Coles et al., 1992) in lactating dairy ewes. On-farm visits were done rapidly following suspicion, whenever possible. However, for farms 2 and 4, the first Covid-19 lockdown in France caused a 2-month delay between suspicion and visit. On 1 farm (farm 2), the lactating ewes were not available (e.g. they were already grazing in high-altitude summer pastures) and due to the emergency of the situation, FECRT was conducted on ewe lambs. The animals were randomly allocated to 4 groups of 10-11 animals according to their age to compose homogenous groups representative of the herd or the age group. A control group was left untreated, 1 group received injectable EPN (0.2 mg kg⁻¹ of LBW, Eprecis® injectable, CEVA Santé Animale, Libourne, France; further referred to as the 'SC group') and another group received a topical 'Pour-On' form of EPN (1 mg kg⁻¹ of LBW, Eprinex Multi®, Boehringer Ingelheim, Lyon, France; 'T group') according to the manufacturer's indication. The last group received EPN orally, using the topical formula (Eprinex Multi®, Boehringer Ingelheim; 'O group') off label and at the dose of 0.5 mg kg⁻¹ of body weight (Badie et al., 2015). All animals were treated with a dose rate of 80 kg, which is heavier than the heaviest animal of the group, the 2 encountered breed weights being on average 55-60 kg for female individuals. For the topical treatment, wool was carefully parted so as to apply the solution as well as possible on the skin. For the oral drench, a graduated single-use syringe was used and the absence of regurgitation was verified after treatment. Groups were marked according to their treatment regimen. Feces was collected individually from all animals, samples were identified using the animals' tag 5 digit number and treatment was administered to the pre-defined groups. Fourteen days after treatment, feces was collected individually from all groups. Animals with no fecal samples collected on day 14 after treatment were excluded from the study.

### FECRT

Individual FECs were conducted using the modified McMaster method with a sensitivity of 15 eggs per gram (EpG) (Raynaud et al., 1970) within a maximum of 48 h after sampling. Animals for which no strongyle eggs were detected at D0 or animals for which no feces was collected post-treatment were excluded from the study. Once the FEC for the D0 samples were established, the mean FEC was calculated per group, as described by Coles et al. (1992). Only groups for which the mean FEC was greater than 300 EPG were included in the study (Cabaret and Berrag, 2004).

Fecal egg count reduction (FECR) and 95% confidence intervals (CIs) were calculated using 3 different formulas, as follows (Coles et al., 1992; Dash et al., 1988; McKenna, 2006):

$$ 
\begin{aligned}
\text{FECR}_1 &= 100 \times (1 - \text{EpG}_{\text{T2}}/\text{EpG}_{\text{C2}}) \\
\text{FECR}_2 &= 100 \times (1 - (\text{EpG}_{\text{T2}} \times \text{EpG}_{\text{C1}}/\text{EpG}_{\text{T1}} \times \text{EpG}_{\text{C2}})) \\
\text{FECR}_3 &= 100 \times (1 - \text{EpG}_{\text{T2}}/\text{EpG}_{\text{T1}})
\end{aligned}
$$

where EpG<sub>T1</sub> and EpG<sub>T2</sub> are the arithmetic means of FEC in a treated group at D0 and D14, respectively, and EpG<sub>C1</sub> and
<table><caption>**Table 1.** Interpretation guide for FECRT results (COMBAR, 2021)</caption><thead><tr><td>Efficacy</td><td>Results</td></tr></thead><tbody><tr><td>Reduced</td><td>FECR &lt; 95% and lower limit of the 95% CI &lt; 90%</td></tr><tr><td>Doubtful</td><td>Either FECR &lt; 95% or lower limit of the 95% CI &lt; 90%</td></tr><tr><td>Normal</td><td>FECR ≥ 95% and lower limit of the 95% CI ≥ 90%</td></tr></tbody></table>

EpG<sub>C2</sub> are the arithmetic means of FEC in the control group, at D0 and D14, respectively.

CIs for these 3 formulas were calculated according to the methods described by Coles *et al.* (FECR<sub>1</sub>) (Coles *et al.*, 1992) and by Lyndal-Murphy *et al.* (FECR<sub>2</sub> and FECR<sub>3</sub>) (Lyndal-Murphy *et al.*, 2014).

The results were interpreted as described in Table 1.

## Larvae collection

After FEC, at D0 and D14, stools were combined by group (control, SC, O, T) for fecal culture. Mixing was done so that when the remaining amount after FEC was sufficient, 3–5 g of feces from each animal was combined into the culture. The composite fecal cultures were then incubated for at least 12 days at 24 ± 1° C, and humidified every 2–3 days with tap water. For larvae collection, pots were filled to the brim with tap water and turned up-side down into Petri dishes, which were in turn filled with water. Larvae were collected twice at a 24 h interval in a volume of 40–45 mL and stored vertically at 4°C until DNA extraction (MAFF, 1986).

## Larvae quantification and identification

The supernatant of the tubes stored at 4°C was discarded, and 5 mL of the pellet containing the larvae was kept for further analysis. Furthermore, 500 µL of the pellet was used for the DNA extraction, using the DNeasy PowerSoil kit (QIAGEN, Hilden, Germany). Molecular identification was then performed using a real-time polymerase chain reaction (PCR) according to Milhes *et al.* (2017). Experiments were based on real-time PCR reactions, and standard curves for larval DNA quantitation were established for each PCR run and for 3 species *H. contortus*, *Teladorsagia circumcincta* and *Trichostrongylus colubriformis*.

## EPN analysis in sheep serum

Blood samples were collected 2 and 5 days post-EPN treatment, in dry tubes from the jugular vein of all treated animals. Blood samples were centrifuged at 3500 rpm for 10 min. Serum was collected and stored at –20°C until further analysis.

After extraction from serum with acetonitrile, EPN concentration was measured using high-performance liquid chromatography with fluorescent detection, as previously described by Sutra *et al.* (1998). The quantification limit of the method was 0.07 ng mL<sup>-1</sup>, and the inter-assay coefficient of variation was lower than 5%.

## Statistical analysis

Graphs were executed using GraphPad Prism version for Windows, GraphPad Software, San Diego, California, USA. Statistical analyses were conducted using R [version 4.1.1 (2021-08-10)] and RStudio version 1.4.1106 (RStudio Team, 2021). Mean EPN concentrations were compared for different treatment regimens within farms using a non-parametric Wilcoxon test.

# Results

## FECRT

Average FEC on D0 in every group and on each farm were above 300 EpG. On farms 2–5, after withdrawing animals not responding to inclusion criteria detailed in ‘FECRT’ section of Materials and methods, 7–11 animals remained per group. At the moment of testing, farm 1 had started the transfer of some animals to summer pastures. The number of animals remaining on the farm was sufficient for 3 groups of 8 lactating ewes (control, treated by subcutaneous and with the topical route) (Table 2).

Calculated FECR are presented in Fig. 1. Depending on the farm, the group and the formula used, FECR results varied widely. All values of FECR were lower than 95%, except for FECR<sub>3</sub> for group O of farm 4 and all lower level CIs were inferior to 90%. With the exception of FECR<sub>3</sub> of group O for farm 4 (97%), these criteria indicate reduced efficacy for all 5 farms. Regarding group O of farm 4, interpretation would have been that efficacy of EPN was doubtful. However, given results obtained with 2 other formulas, including FECR<sub>1</sub> recommended by the WAAVP guidelines, efficacy for group O of farm 4 is clearly reduced.

Different FECR formulas yield different results, depending mainly on the value of the mean FEC of different groups. FECR<sub>1</sub> and FECR<sub>3</sub> yield similar results when the mean FEC of the control group on D14 is close to the mean FEC of treated groups on D0, as is the case on farm 2. On this farm, the FEC on D14 for the control group (3756 EpG) is not significantly different from FECs for the SC, O and T groups (respectively 4036, 1445 and 5072 EpG) on D0 (*P* value < 0.05). FECR<sub>1</sub> formula that compares the FEC of treated and control groups on D14 yields similar percentages to FECR<sub>3</sub>, which compares FEC of treated group on D14 and D0, in this case (e.g. 87 and 88% reduction for the SC group with FECR<sub>1</sub> and FECR<sub>3</sub>, respectively). On farm 4, percentages of fecal egg reduction differ between FECR<sub>1</sub> and FECR<sub>3</sub>. There is a significant difference between the FEC of the control group on D0 (194 EpG) and of the treated groups on D14 (1813, 1986 and 1338 EpG for SC, O and T groups, respectively, *P* < 0.05). For example, considering the SC group on this farm, post-treatment fecal egg reduction is –25% using FECR<sub>1</sub>, yet it is 83% using FECR<sub>3</sub>.

Mean EPN concentrations in serum 2 days after administration were significantly higher after subcutaneous injection than after oral administration in 3 out of the 4 farms where both these routes were tested (Fig. 2; farms 2, 3 and 5; *P* value < 0.01). On farm 4, there was no significant difference between both routes. Topical administration of EPN resulted in dramatically low drug concentrations in the serum of ewes in all 5 farms (Fig. 2; *P* value < 0.01). Mean values were between 14.35 [s.d.: 4.43] and 27.84 ng mL<sup>-1</sup> [s.d.: 6.37] for the SC group; 5.06 [s.d.: 5.71] and 25.14 ng mL<sup>-1</sup> [s.d.: 9.74] for the O group and between 0.97 [s.d.: 0.54] and 6.92 ng mL<sup>-1</sup> [s.d.: 4.2] for the T group, 2 days after treatment (Fig. 2). Five days after treatment, EPN serological concentrations were significantly higher after subcutaneous injection than after either of the other routes.

For the T groups, the individual concentrations of EPN varied from being all below 2 ng mL<sup>-1</sup> (farm 5) to all above (farm 1). On farms 2, 3 and 4, respectively, 44, 60 and 50% of individual EPN concentration values in the T group were below this concentration threshold (Fig. 2 and Supplementary data).

## Strongyle species present pre- and post treatment

Before treatment, at D0, the strongyle species most present in the cultured feces of every farm was *H. contortus*. On farm 1, it was the only species of the 3 the qPCR could identify. On the 4 other farms, *T. colubriformis* and *T. circumcincta* were present
<table><caption>**Table 2.** Mean FEC results, with minimal and maximal individual value and final number of animals included, per group and per farm</caption><thead><tr><th rowspan="2">Farm</th><th rowspan="2"></th><th colspan="2">Control</th><th colspan="2">SC</th><th colspan="2">O</th><th colspan="2">T</th></tr><tr><th>D0</th><th>D14</th><th>D0</th><th>D14</th><th>D0</th><th>D14</th><th>D0</th><th>D14</th></tr></thead><tbody><tr><td rowspan="3">1</td><td>Mean</td><td><strong>850</strong></td><td>1956</td><td><strong>2019</strong></td><td>2689</td><td></td><td></td><td><strong>1408</strong></td><td>919</td></tr><tr><td>Min-Max</td><td><strong>100-2250</strong></td><td>0-7100</td><td><strong>100-10 000</strong></td><td>0-11 750</td><td></td><td></td><td><strong>150-5300</strong></td><td>50-4800</td></tr><tr><td>N</td><td>8</td><td></td><td>8</td><td></td><td></td><td></td><td>8</td><td></td></tr><tr><td rowspan="3">2<sup>1</sup></td><td>Mean</td><td><strong>3028</strong></td><td>3756</td><td><strong>4036</strong></td><td>486</td><td><strong>1445</strong></td><td>400</td><td><strong>5072</strong></td><td>394</td></tr><tr><td>Min-Max</td><td><strong>75-6550</strong></td><td>300-7550</td><td><strong>800-17 200</strong></td><td>100-2000</td><td>450-2750</td><td>0-1100</td><td><strong>850-12 700</strong></td><td>50-950</td></tr><tr><td>N</td><td>8</td><td></td><td>11</td><td></td><td>10</td><td></td><td>9</td><td></td></tr><tr><td rowspan="3">3</td><td>Mean</td><td><strong>1428</strong></td><td>800</td><td><strong>662</strong></td><td>650</td><td><strong>1089</strong></td><td>1039</td><td><strong>1205</strong></td><td>885</td></tr><tr><td>Min-Max</td><td><strong>100-6550</strong></td><td>50-2250</td><td><strong>15-2050</strong></td><td>0-1950</td><td>50-2550</td><td>0-2550</td><td><strong>50-2350</strong></td><td>0-2200</td></tr><tr><td>N</td><td>9</td><td></td><td>10</td><td></td><td>9</td><td></td><td>10</td><td></td></tr><tr><td rowspan="3">4</td><td>Mean</td><td><strong>1693</strong></td><td>194</td><td><strong>1813</strong></td><td>302</td><td><strong>1986</strong></td><td>50</td><td><strong>1338</strong></td><td>581</td></tr><tr><td>Min-Max</td><td><strong>100-8500</strong></td><td>45-500</td><td><strong>150-5300</strong></td><td>0-2000</td><td>150-8200</td><td>0-200</td><td><strong>50-3250</strong></td><td>50-1300</td></tr><tr><td>N</td><td>7</td><td></td><td>8</td><td></td><td>7</td><td></td><td>8</td><td></td></tr><tr><td rowspan="3">5</td><td>Mean</td><td><strong>469</strong></td><td>525</td><td><strong>376</strong></td><td>413</td><td><strong>838</strong></td><td>980</td><td><strong>735</strong></td><td>1020</td></tr><tr><td>Min-Max</td><td><strong>50-1550</strong></td><td>0-1800</td><td><strong>60-1700</strong></td><td>0-1600</td><td>50-2650</td><td>100-2700</td><td><strong>350-1200</strong></td><td>300-1900</td></tr><tr><td>N</td><td>10</td><td></td><td>8</td><td></td><td>10</td><td></td><td>10</td><td></td></tr></tbody></table>

In bold, values of mean and minimum-maximum FEC on day 0, for an easier quick reading.

in small proportions: *T. colubriformis* larvae composed at most 16.4% (on farm 3) of the larval culture yield (Table 3).

After treatment, *H. contortus* was the main species identified on the 5 farms (Table 3). For 4 farms out of 5, *H. contortus* was the only species remaining for the groups treated by subcutaneous injection. On farm 3, *T. colubriformis* larvae were also present in the cultures post injectable EPN treatment (Table 3). On farm 2, *T. colubriformis* larvae were the only species present post oral treatment, however in small quantities (Supplementary data, Table S8). *Haemonchus contortus* and *T. colubriformis* larvae were present in the post-treatment fecal cultures of the ewes treated with a topical solution, except in farm 1 where the collected larvae were 100% *H. contortus*.

Discussion

In this field study, EPN was found to have a reduced efficacy whatever the formula tested (SC, O, T) in 5 commercial dairy sheep farms where veterinarians and farmers suspected a lack of efficacy based on the persistence of clinical signs following oral or injectable treatment. Low drug efficacy was observed even in the SC administration group, which showed the highest drug levels in the hosts' serum.

We confirmed that EPN serum levels are highly dependent upon the route of administration. EPN concentration measured in the serum of treated lactating and pre-lactating ewes 2 days after treatment was well above the 2 ng mL⁻<fnref n="1" /> minimal efficacy concentration for all the animals treated with a subcutaneous or oral formula, which supports the idea that these ewes received a dose of the molecule that can be considered sufficient to kill strongyle adults. The poor FECRT performance obtained in these animals strongly supports the presence of strongyles resistant to EPN. Although no pharmacokinetic and pharmacodynamic (PK/PD) study has been conducted to specifically identify the minimal EPN therapeutic dose, the minimal active dose for this drug family (e.g. ivermectin) has been shown to be above 2 ng mL⁻<fnref n="1" /> (Bousquet-Mélou *et al.*, 2011). This concentration is

therefore considered a threshold that guarantees efficacy of EPN in small ruminants (Hoste *et al.*, 2004; Rostang *et al.*, 2020). The time at which EPN concentration reaches its highest averages between 1.2(±0.4) and 3.13(±2.99) days depending on dosage rate and physiology of the animals was considered (Imperiale *et al.*, 2006; Hodošček *et al.*, 2008; Hamel *et al.*, 2017). Average serum EPN concentrations in the groups treated with a topical solution of EPN were low in most of the farms on D2. Differences in the breed and physiology of the animals could explain values below those found by Hamel *et al.* (2017) in dry merino crossed sheep. In the present study, on all farms except farm 2, FECR tests were conducted on lactating ewes and in all farms animals bore a substantial worm burden, and both lactation and body condition have been shown to influence ML pharmacokinetic parameters (Lespine *et al.*, 2004, 2012; Rostang *et al.*, 2020).

The purpose of this study was to provide reliable information about AR status in farms. We set up a feasible protocol to monitor drug efficacy through combining FECR and drug concentration monitoring in treated animals. Measuring concentrations at 2 critical times 2 days (close to maximal concentration) and 5 days (elimination phase), these data points allow simulation of the complete drugs' pharmacokinetics.

FECR values indicate a reduced efficacy with all 3 formulas and the low FECR after a topical treatment is due to the presence of resistant strongyle. However, cases of underexposure of GIN to EPN when using a topical formula have previously been reported by veterinary practitioners and confirmed by 2 recent studies in France. Bouy *et al.* and Bordes *et al.* described cases where FECR after EPN treatment were below 95% when using the topical solution and above 95% when animals of the same flock were treated with a subcutaneous solution (Bouy *et al.*, 2021; Bordes *et al.*, 2022). In the study by Bordes *et al.* (2022) serum concentrations of EPN were below 2 ng mL⁻<fnref n="1" /> for all animals treated with a topical formula. These findings are in line with others (Hoste *et al.*, 2004; Hodošček *et al.*, 2008) that underline the highly variable bioavailability of the topical formula of EPN. The use of topical route for EPN is therefore not recommended.

<footnote n="1">¹On farm 2, FECRT were conducted on ewe lambs.</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 1.** FECR results and confidence intervals for the 5 farms, calculated according to 3 different formulas and for all treatment types.</figcaption></figure>

The significant difference between mean concentrations of EPN in SC and T groups is observed although the dose rate for the topical solution is 5 times higher than for the injection solution (1 mg kg⁻¹ of LBW for the topical solution and 0.2 mg kg⁻¹ LBW for the injection solution). Given the impact of ML on non-target species such as dung beetles, elimination of the molecule, through direct contact or by the fecal matters, could have an impact on pasture quality (Sands and Wall, 2018; Verdú *et al.*, 2018; Weaving *et al.*, 2020).

Upon communication of these results back to farmers and veterinarians, one of the challenges was explaining the discrepancy between the reduction calculated for the SC group and the one calculated for the oral drench group on the same farm at the same date, therefore with the same control group. Questions arose for these 2 treatment regimens although not for the topical formula group for the reasons explained above. Much has already

been said about FECRT and their limits, and the debate remains open as summarized very recently by Morgan *et al.* (2022). However, our field experience teaches us that such differences need explanations, in order to be accepted by farmers and veterinarians. In our study, FECRT variations could hardly be explained by pharmacological factors, as EPN concentrations were well above 2 ng mL⁻¹, a concentration at which the molecule is considered to be efficient (Guillot *et al.*, 1986; Guyonnet *et al.*, 2017). Although ewes of similar ages were evenly distributed between the SC and O group, differences in mean EPG partly explain the differences observed between the FECR of these 2 treatment formulas. Host factors, such as consistency of feces could contribute to variations in egg outputs and contribute to variations from one group to another, and host immunity could also be of importance for the inter-individual variations. *Haemonchus contortus* being the main strongyle species before
<figure><img src="image_5.png" /><figcaption>**Fig. 2.** Serological concentrations of EPN by farm, and by administration route. Horizontal lines indicate mean concentrations per treatment type and per day post-treatment. Black triangle: SC group; medium grey diamond: O group; light grey dot: T group. Red dotted line is at 2 ng mL⁻¹. D2: 2 days after treatment; D5: 5 days after treatment. Individual concentrations and mean concentrations per group ± s.d. (7–11 animals per group).</figcaption></figure>

and after treatment in all 5 farms, variations in FECR due to initial diversity of species is limited. Parasite fitness, including its fecundity, could add to the observed variations, although this should be limited by the fact that it is highly probable all sheep of the same flock harbour the same resistant isolate.

Identification of the 3 main pathogenic strongyle species for small ruminants was done using fecal cultures and qPCR. In all 5 farms, *H. contortus* was the main species present after treatment. Furthermore, resistance to EPR of isolates from farms 1 and 4 was confirmed by infestation and EPN challenge in experimental sheep (G. Salle, unpublished data). Interestingly, before treatment, *H. contortus* was predominant in the fecal culture yields of all 5

farms, and on farm 1 it was the only species present. These farms were investigated because EPN treatments were not resolving the observed symptoms, although the farmers had sometimes drenched the animals with EPN several times before calling their veterinarians to alert them to the problem. Our hypothesis is that these repeated treatments have exerted an important selective pressure upon the present worms, and cleared all susceptible populations. Farms 1, 4 and 5 were included in this study because during the spring of 2020 they were facing dire situations. On farms 2 and 3, health issues were sub-acute. Veterinarians knew the history of these 5 farms and that they are prone to facing GIN issues. This first report of EPN resistance in dairy sheep in
<table><caption>**Table 3.** Proportions (%) of the 3 main pathogenic species in bulk fecal cultures before (D0) and after (D14) for every group of the 5 farms</caption><thead><tr><th rowspan="3">Farm</th><th colspan="4">Control</th><th colspan="4">SC</th><th colspan="4">O</th><th colspan="4">T</th></tr><tr><th colspan="2">D0</th><th colspan="2">D14</th><th colspan="2">D0</th><th colspan="2">D14</th><th colspan="2">D0</th><th colspan="2">D14</th><th colspan="2">D0</th><th colspan="2">D14</th></tr><tr><th>Hc</th><th>Te</th><th>Tr</th><th>Hc</th><th>Te</th><th>Hc</th><th>Te</th><th>Hc</th><th>Tr</th><th>Hc</th><th>Te</th><th>Hc</th><th>Tr</th><th>Hc</th><th>Te</th><th>Hc</th><th>Tr</th></tr></thead><tbody><tr><td>1</td><td>100</td><td>0</td><td>0</td><td>100</td><td>0</td><td>0</td><td>100</td><td>0</td><td>0</td><td>100</td><td>0</td><td>0</td><td>100</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>2<sup>a</sup></td><td>97.4</td><td>2.2</td><td>0.4</td><td>99.2</td><td>0</td><td>0.8</td><td>99.3</td><td>0.6</td><td>0.2</td><td>100</td><td>0</td><td>0</td><td>97.8</td><td>0.4</td><td>1.8</td><td>0</td><td>0</td></tr><tr><td>3</td><td>94.5</td><td>1.2</td><td>4.3</td><td>88.8</td><td>4.6</td><td>6.6</td><td>99.4</td><td>0.6</td><td>0</td><td>65.8</td><td>0</td><td>34.2</td><td>78.2</td><td>15.2</td><td>6.7</td><td>100</td><td>0</td></tr><tr><td>4</td><td>75.5</td><td>12.2</td><td>12.4</td><td>100</td><td>0</td><td>0</td><td>88.8</td><td>0.3</td><td>10.9</td><td>100</td><td>0</td><td>0</td><td>78.8</td><td>8.9</td><td>12.3</td><td>100</td><td>0</td></tr><tr><td>5</td><td>99.9</td><td>0.1</td><td>0</td><td>100</td><td>0</td><td>0</td><td>99.8</td><td>0</td><td>0</td><td>100</td><td>0</td><td>0</td><td>100</td><td>0</td><td>0</td><td>0</td><td>0</td></tr></tbody></table>

this dairy region is a description of how resistance can manifest itself and how it has been investigated. This study was conducted to investigate issues raised by some farmers and veterinarians, and the aim was not to determine the prevalence of EPN resistance. Hence, to this date, it is not known if these cases are the tip of the iceberg regarding EPN resistance in the *Pyrénées Atlantiques*, and if they are, how big the iceberg is.

In post topical treatment coprocultures from farms 2 to 5, *T. colubriformis* and *T. circumcincta* (on farm 3) were also identified. Considering the EPR serum concentration for this treatment formula, we hypothesize these worms are still present because they have been under-exposed to the molecule. Furthermore, the intestinal species *T. colubriformis* has already been described as dose-limiting for EPN, its presence after topical treatment is not surprising (Chartier *et al.*, 1999; Hoste *et al.*, 2004).

In this southwestern area of France, resistance of *H. contortus* to EPN had previously been described in lactating goats (Bordes *et al.*, 2020) and in sheep raised for meat (Cazajous *et al.*, 2018), but this is the first report of the presence of resistant strongyle strains in dairy sheep in the *Pyrénées Atlantiques*. Special attention has been brought to EPN resistance in this area, as well as in the Roquefort area, being the only molecule with a zero milk withdrawal period. Hence, EPN resistance means treating ewes against GIN during their lactation will inevitably come at the cost of, at least, throwing away milk during the withdrawal period. Farmers have to find other ways of controlling parasitism than solely relying on drugs, in order to keep their grazing flocks healthy while maintaining a decent level of production. The presence of resistant strains is also of concern here because of the use of collective pastures that could allow for the dissemination of resistant strains.

Resistance of GIN to ML molecules in France appeared relatively late, compared to other countries: Geurden *et al.* (2014) and Paraud *et al.* (2010) reported no resistance to ML in the 2 main dairy sheep regions and in dairy goats, respectively, and the first case of ivermectin resistance was described for *T. circumcincta* in meat sheep production in central France in 2014 (Paraud *et al.*, 2016). These observations are in line with the general trend in Europe of increasing AR throughout anthelmintic families and GIN species (Rose Vineer *et al.*, 2020). Of the 3 main pathogenic GIN species for sheep, *H. contortus*, originally a parasite of warm and humid climates (Sallé *et al.*, 2019; Arsenopoulos *et al.*, 2021), profits from global warming and is reaching farther north in Europe and infection pressure expands during the year (Rose *et al.*, 2016). Although located in southern France, the area of interest for the present study harbours an oceanic climate with normally mild temperatures. In the last few years, it seems to be hosting more frequent cases of severe haemonchosis (P. Jacquiet, T. Cazajous, personal communication), in line with the rise of average temperatures. Four of the 5 cases of haemonchosis presented here happened after a particularly mild winter (MeteoFrance, 2020), which could have allowed an increased survival of *H. contortus* larvae on pastures. As well as adapting to climate, *H. contortus* has also been adapting to ML treatments, and in addition to the present study, resistance of the Barber Pole Worm has also been recently described in sheep in southwest England (Bull *et al.*, 2022) and in goats in Austria (Hinney *et al.*, 2022).

## Conclusion

The present study reports for the first time in France 5 cases of EPN-resistant *H. contortus* in dairy ewes in the south western *département* of the *Pyrénées Atlantiques*, the country's second largest sheep milk and cheese production area. These isolated cases

<footnote>Hc, Haemonchus contortus; Te, Teladorsagia circumcincta; Tr, Trichostrongylus colubriformis; SC, subcutaneous EPN; O, oral drench EPN; T, topical EPN

ªOn farm 2, FECRT were conducted on ewe lambs.</footnote>
are particularly worrying for the dairy sheep production in this area, where an important percentage of farms rely on summer pastures, usually shared, as part of their forage resource. However, they reflect what is happening elsewhere in Europe, in meat but also dairy sheep production. For milk-producing farms, EPN resistance in GIN does not come as a complete surprise, as only a very small pool of molecules is used for treatment. This study further confirms the variability of EPN serum concentration depending on administration route, with the injectable formula yielding the highest concentration post-treatment. The topical solution yields highly variable and sub-therapeutic EPN concentrations. Therefore, the use of topical EPN should be discouraged so as not to exacerbate the problem of resistance to ML and to better maintain good animal health. Dairy production adds an extra challenge to field management, yet some farmers have already found some encouraging, and hopefully durable, improvements to their system. For the farms not yet facing resistance of GIN to EPN, a simple and robust protocol to target and selectively treat lactating ewes to maintain a refuge population is in trial.

**Supplementary material.** The supplementary material for this article can be found at https://doi.org/10.1017/S0031182023000069.

**Data availability.** Data supporting results are provided within the article and in the supplementary materials.

**Acknowledgements.** The authors are very grateful to the farmers who took part in this study, for their cooperation and for granting them access to their flocks.

**Author's contributions.** P. J. and A. L. conceived and designed the study. L. B., P. J., S. J., T. C., J. L., N. D., M. C., C. V.-N., M. D., M. A. and C. D. carried out the fieldwork. L. B., C. G., S. J., J. F. S., M. D., M. A. and C. D. carried out laboratory analysis. S. J., L. B., N. D., M. D., D. A., H. K., A. L. and P. J. wrote the article.

**Financial support.** Research presented in this study has received financial support from the Nouvelle Aquitaine Region 'Paralut' project, from the France Futur Elevage (F2E Institut-Carnot Santé Animale) 'Antherin' project and CEVA Santé Animale.

**Conflict of interest.** S. J., D. A. and H. K. are employees of CEVA Santé Animale. All work herein was conducted in academic settings, and the authors declare no conflict of interest regarding the present results.

**Ethical standards.** Stool and blood collection and anthelmintic treatments are a part of routine veterinary procedures without any traumatic method. Such procedures are not qualified as animal experimentation involving vertebrates according to French laws, so no specific ethical clearance was required.

## References

* Arsenopoulos KV, Fthenakis GC, Katsarou EI and Papadopoulos E (2021) Haemonchosis: a challenging parasitic infection of sheep and goats. *Animals* **11**, 363.

* Badie C, Lespine A, Devos J, Sutra JF and Chartier C (2015) Kinetics and anthelmintic efficacy of topical eprinomectin when given orally to goats. *Veterinary Parasitology* **209**, 56–61.

* Bordes L, Dumont N, Lespine A, Souil E, Sutra J-F, Prévot F, Grisez C, Romanos L, Dailledouze A and Jacquiet P (2020) First report of multiple resistance to eprinomectin and benzimidazole in *Haemonchus contortus* on a dairy goat farm in France. *Parasitology International* **76**, 102063.

* Bordes L, Ticoulet D, Sutra JF, Lespine A and Jacquiet P (2022) Lack of efficacy of topical administration of eprinomectin against gastrointestinal nematode in a French dairy sheep farm: a case of underexposure of worms. *Veterinary Record Case Reports*, e435. doi: 10.1002/vrc2.435

* Bousquet-Mélou A, Jacquiet P, Hoste H, Clément J, Bergeaud J-P, Alvinerie M and Toutain P-L (2011) Licking behaviour induces partial anthelmintic efficacy of ivermectin pour-on formulation in untreated cattle. *International Journal for Parasitology* **41**, 563–569.

* Bouy M, Fito-Boncompte L, Harinck E, Lukkes S and Heckendorn F (2021) Echecs de traitement anthelminthiques à base d'éprinomectine sur des petits ruminants: resistance ou voie d'administration inappropriée? *Le Nouveau Praticien Vétérinaire - Elevage et santé* **13**, 126-130.

* Bull K, Glover MJ, Rose Vineer H and Morgan ER (2022) Increasing resistance to multiple anthelmintic classes in gastrointestinal nematodes on sheep farms in southwest England. *Veterinary Record* **190**, e1531. doi: 10.1002/vetr.1531

* Cabaret J and Berrag B (2004) Faecal egg count reduction test for assessing anthelmintic efficacy: average versus individually based estimations. *Veterinary Parasitology* **121**, 105-113.

* Cazajous T, Prevot F, Kerbiriou A, Milhes M, Grisez C, Tropee A, Godart C, Aragon A and Jacquiet P (2018) Multiple-resistance to ivermectin and benzimidazole of a *Haemonchus contortus* population in a sheep flock from mainland France, first report. *Veterinary Parasitology: Regional Studies and Reports* **14**, 103-105.

* Chartier C, Etter E, Pors I and Alvinerie M (1999) Activity of eprinomectin in goats against experimental infections with *Haemonchus contortus*, *Teladorsagia circumcincta* and *Trichostrongylus colubriformis*. *Veterinary Record* **144**, 99-100.

* Coles GC, Bauer C, Borgsteede FHM, Geerts S, Klei TR, Taylor MA and Waller PJ (1992) World Association for the Advancement of Veterinary Parasitology (W.A.A.V.P.) methods for the detection of anthelmintic resistance in nematodes of veterinary importance. *Veterinary Parasitology* **44**, 35-44.

* COMBAR (2021). Fecal Egg Count Reduction Test (FECRT) protocol. Gastrointestinal nematodes - Sheep and Goats. https://www.combar-ca.eu/sites/default/files/FECRT_PROTOCOL_sheep.goats_March%202021.pdf

* Dash KM, Hall E and Barger IA (1988) The role of arithmetic and geometric mean worm egg counts in faecal egg count reduction tests and in monitoring strategic drenching programs in sheep. *Australian Veterinary Journal* **65**, 66-68.

* Falzon LC, O'Neill TJ, Menzies PI, Peregrine AS, Jones-Bittoon A, vanLeeuwen J and Mederos A (2014) A systematic review and meta-analysis of factors associated with anthelmintic resistance in sheep. *Preventive Veterinary Medicine* **117**, 388-402.

* Geurden T, Hoste H, Jacquiet P, Traversa D, Sotiraki S, Frangipane di Regalbono A, Tzanidakis N, Kostopoulou D, Gaillac C, Privat S, Giangaspero A, Zanardello C, Noé L, Vanimisetti B and Bartram D (2014) Anthelmintic resistance and multidrug resistance in sheep gastrointestinal nematodes in France, Greece and Italy. *Veterinary Parasitology* **201**, 59-66.

* Guillot FS, Wright FC and Oehler D (1986) Concentration of ivermectin in bovine serum and its effect on the fecundity of psoroptic mange mites. *American Journal of Veterinary Research* **47**, 525-527.

* Guyonnet J, Karembe H, Magnier R and Menudier N (2017) *Pharmacocinétique et déplétion des résidus d'éprinomectine dans le lait après adminstraton sous-cutanée d'Eprecis chez la brebis*. Paris, France: Société Nationale des Groupements Techniques Vétérinaires.

* Hamel D, Bosco A, Rinaldi L, Cringoli G, Kaulfuß K-H, Kellermann M, Fischer J, Wang H, Kley K, Mayr S, Rauh R, Visser M, Wiefel T, Fankhauser B and Rehbein S (2017) Eprinomectin pour-on (EPRINEX® Pour-on, Merial): efficacy against gastrointestinal and pulmonary nematodes and pharmacokinetics in sheep. *BMC Veterinary Research* **13**, 148.

* Hinney B, Wiedermann S, Kaiser W, Krücken J and Joachim A (2022) Eprinomectin and moxidectin resistance of trichostrongyloids on a goat farm in Austria. *Pathogens* **11**, 498.

* Hodošček L, Grabnar I, Milčinski L, Süssinger A, Eržen NK, Zadnik T, Pogačnik M and Cerkvenik-Flajs V (2008) Linearity of eprinomectin pharmacokinetics in lactating dairy sheep following pour-on administration: excretion in milk and exposure of suckling lambs. *Veterinary Parasitology* **154**, 129-136.

* Hoste H, Lespine A, Lemercier P, Alvinerie M, Jacquiet P and Dorchies P (2004) Efficacy of eprinomectin pour-on against gastrointestinal nematodes and the nasal bot fly (*Oestrus ovis*) in sheep. *Veterinary Record* **154**, 782-785. https://doi.org/10.1136/vr.154.25.782

* HPRA (2016) *Publicly Available Assessment Report for a Veterinary Medicinal Product: Eprinex Multi 5 mg/ml Pour-on for Beef and Dairy Cattle, Sheep and Goats.*, IE/V/0347/001/DC. Dublin, Ireland: HPRA.

* HPRA (2020) *Publicly Available Assessment Report for a Veterinary Medicinal Product: Eprecis 20 mg/ml Solution for Injection for Cattle, Sheep and Goats*, IE/V/0340/001/DC. Dublin, Ireland: HPRA.
Imperiale F and Lanusse C (2021) The pattern of blood-milk exchange for antiparasitic drugs in dairy ruminants. *Animals* **11**, 2758.

Imperiale F, Pis A, Sallovitz J, Lisfchitz A, Busetti M, Suárez V and Lanusse C (2006) Pattern of eprinomectin milk excretion in dairy sheep unaffected by lactation stage: comparative residual profiles in dairy products. *Journal of Food Protection* **69**, 2424-2429.

INAO (2015) *Cahier des Charges de l'appellation d'origine 'Ossau-iraty', décret 2015-1225*. Montreuil, France: INAO.

Kotze AC and Prichard RK (2016) Anthelmintic resistance in *Haemonchus contortus*. *Advances in Parasitology* **93**, 397-428. doi: 10.1016/bs.apar.2016.02.012

Lespine A, Sutra J-F, Dupuy J, Alvinerie M and Aumont G (2004) The influence of parasitism on the pharmacokinetics of moxidectin in lambs. *Parasitology Research* **93**, 121-126.

Lespine A, Chartier C, Hoste H and Alvinerie M (2012) Endectocides in goats: pharmacology, efficacy and use conditions in the context of anthelmintics resistance. *Small Ruminant Research* **103**, 10-17.

Lyndal-Murphy M, Swain AJ and Pepper PM (2014) Methods to determine resistance to anthelmintics when continuing larval development occurs. *Veterinary Parasitology* **199**, 191-200.

MAFF (1986) *Manual of Veterinary Parasitological Laboratory Techniques*. London: HMSO.

McKenna P (2006) A comparison of faecal egg count reduction test procedures. *New Zealand Veterinary Journal* **54**, 202-203.

Meteo France (1991-2020). meteofrance.com.

MeteoFrance (2020) *France: l'hiver le plus chaud jamais mesuré*. Saint Mandé, France: MeteoFrance.

Milhes M, Guillerm M, Robin M, Eichstadt M, Roy C, Grisez C, Prévot F, Liénard E, Bouhsira E, Franc M and Jacquiet P (2017) A real-time PCR approach to identify anthelmintic-resistant nematodes in sheep farms. *Parasitology Research* **116**, 909-920.

Morgan ER, Lanusse C, Rinaldi L, Charlier J and Vercruysse J (2022) Confounding factors affecting faecal egg count reduction as a measure of anthelmintic efficacy. *Parasite* **29**, 20.

Murri S, Knubben-Schweizer G, Torgerson P and Hertzberg H (2014) Frequency of eprinomectin resistance in gastrointestinal nematodes of goats in canton Berne, Switzerland. *Veterinary Parasitology* **203**, 114-119.

O'Connor LJ, Walkden-Brown SW and Kahn LP (2006) Ecology of the free-living stages of major trichostrongylid parasites of sheep. *Veterinary Parasitology* **142**, 1-15.

Paraud C, Pors I, Rehby L and Chartier C (2010) Absence of ivermectin resistance in a survey on dairy goat nematodes in France. *Parasitology Research* **106**, 1475-1479.

Paraud C, Marcotty T, Lespine A, Sutra JF, Pors I and Devos I (2016) Cross-resistance to moxidectin and ivermectin on a meat sheep farm in France. *Veterinary Parasitology* **226**, 88-92.

Raynaud J-P, William G and Brunault G (1970) Etude de l'efficacité d'une technique de coproscopie quantitative pour le diagnostic de routine et le contrôle des infestations parasitaires des bovins, ovins, équins et porcins. *Annales de Parasitologie Humaine et Comparée* **45**, 321-342.

Rose H, Caminade C, Bolajoko MB, Phelan P, Van Dijk J, Baylis M, Williams D and Morgan ER (2016) Climate-driven changes to the spatio-temporal distribution of the parasitic nematode. *Global Change Biology* **22**, 1271-1285. doi: 10.1111/gcb.13132

Rose Vineer H, Morgan ER, Hertzberg H, Bartley DJ, Bosco A, Charlier J, Chartier C, Claerebout E, de Waal T, Hendrickx G, Hinney B, Höglund J, Ježek J, Kašný M, Keane OM, Martínez-Valladares M, Mateus TL, McIntyre J, Mickiewicz M, Munoz AM, Phythian CJ, Ploeger HW, Rataj AV, Skuce PJ, Simin S, Sotiraki S, Spinu M, Stuen S, Thamsborg SM, Vadlejch J, Varady M, von Samson-Himmelstjerna G and Rinaldi L (2020) Increasing importance of anthelmintic resistance in European livestock: creation and meta-analysis of an open database. *Parasite* **27**, 69.

Rostang A, Devos J and Chartier C (2020) Review of the eprinomectin effective doses required for dairy goats: where do we go from here? *Veterinary Parasitology* **277**, 108992.

RStudio Team (2021) *RStudio: Integrated Development for R*. Boston, MA: RStudio, PBC.

Sallé G, Doyle SR, Cortet J, Cabaret J, Berriman M, Holroyd N and Cotton JA (2019) The global diversity of *Haemonchus contortus* is shaped by human intervention and climate. *Nature Communications* **10**, 4811.

Sands B and Wall R (2018) Sustained parasiticide use in cattle farming affects dung beetle functional assemblages. *Agriculture, Ecosystems & Environment* **265**, 226-235.

Sangster NC, Cowling A and Woodgate RG (2018) Ten events that defined anthelmintic resistance research. *Trends in Parasitology* **34**, 553-563.

Sutra JF, Galtier P, Alvinerie M and Chartier C (1998) Determination of eprinomectin in plasma by high-performance liquid chromatography with automated solid phase extraction and fluorescence detection. *The Analyst* **123**, 1525-1527.

Verdú JR, Lobo JM, Sánchez-Piñero F, Gallego B, Numa C, Lumaret J-P, Cortez V, Ortiz AJ, Tonelli M, García-Teba JP, Rey A, Rodríguez A and Durán J (2018) Ivermectin residues disrupt dung beetle diversity, soil properties and ecosystem functioning: an interdisciplinary field study. *Science of the Total Environment* **618**, 219-228.

Weaving H, Sands B and Wall R (2020) Reproductive sublethal effects of macrocyclic lactones and synthetic pyrethroids on the dung beetle *Onthophagus similis*. *Bulletin of Entomological Research* **110**, 195-200.

Wolstenholme AJ, Fairweather I, Prichard R, von Samson-Himmelstjerna G and Sangster NC (2004) Drug resistance in veterinary helminths. *Trends in Parasitology* **20**, 469-476.

Zoetis France (2009) Résumé des Caractéristiques du Produit - Valbazen Moutons et Chèvres 1.9%. Available at http://www.ircp.anmv.anses.fr.