# Meta-analysis of the parasitic phase traits of Haemonchus contortus infection in sheep

Mathilde Saccareau, Guillaume Salle, Christèle Robert-Granié, Tom

Duchemin, Philippe Jacquiet, Alexandra Blanchard, Jacques Cabaret, Carole

Moreno-Romieux

► To cite this version:

Mathilde Saccareau, Guillaume Salle, Christèle Robert-Granié, Tom Duchemin, Philippe Jacquiet, et al.. Meta-analysis of the parasitic phase traits of Haemonchus contortus infection in sheep. Parasites & Vectors, 2017, 10, 14 p. 10.1186/s13071-017-2131-7. hal-01607543

HAL Id: hal-01607543

https://hal.science/hal-01607543

Submitted on 26 May 2020

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
RESEARCH

Open Access

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

CrossMark

# Meta-analysis of the parasitic phase traits of *Haemonchus contortus* infection in sheep

Mathilde Saccareau¹*, Guillaume Sallé², Christèle Robert-Granié¹, Tom Duchemin¹, Philippe Jacquiet³,⁴, Alexandra Blanchard², Jacques Cabaret² and Carole R. Moreno¹

## Abstract

**Background:** The parasitic nematode *Haemonchus contortus* shows highly variable life history traits. This highlights the need to have an average estimate and a quantification of the variation around it to calibrate epidemiological models.

**Methods:** This paper aimed to quantify the main life history traits of *H. contortus* and to identify explanatory factors affecting these traits using a powerful method based on a systematic review and meta-analysis of current literature. The life history traits considered are: (i) the establishment rate of ingested larvae; (ii) the adult mortality rate; (iii) the fertility (i.e. the number of eggs laid/female/day); and (iv) fecundity of female worms (i.e. the number of eggs per gram of faeces).

**Results:** A total of 37 papers that report single experimental infection with *H. contortus* in sheep and published from 1960 to 2015, were reviewed and collated in this meta-analysis. This encompassed 115 experiments on 982 animals. Each trait was analysed using a linear model weighted by its inverse variance. The average (± SE) larval establishment rate was 0.24 ± 0.02, which decreased as a function of the infection dose and host age. An average adult mortality rate of 0.021 ± 0.002) was estimated from the literature. This trait varied as a function of animal age, breed and protective response due to prior exposure to the parasite. Average female fertility was 1295.9 ± 280.4 eggs/female/day and decreased in resistant breeds and previously infected hosts. Average faecal egg count at necropsy was 908.5 ± 487.1 eggs per gram of faeces and varied as a function of infection duration and host resistance. The average sex ratio of *H. contortus* was 0.51 ± 0.006.

**Conclusion:** This work is the first systematic review to summarise the available information on the parasitic phase of *H. contortus* in sheep. The results of the meta-analysis provide robust estimates of life history traits for parametrization of epidemiological models, their expected variation according to experimental factors, and provides correlations between these.

**Keywords:** *Haemonchus contortus*, Gastro-intestinal nematode, Sheep, Meta-analysis, Resistance, Female fertility, Larval establishment, Adult mortality, Sex ratio

## Background

The blood-feeding nematode *Haemonchus contortus* is one of the most pathogenic parasites in sheep [1]. The wide geographic distribution and increasing resistance against anthelmintic control measures has made this species a leading threat to the sustainability of sheep industries in Oceania, warm-humid regions [2] and Europe [3].

The critical impact on livestock production combined with the relative ease of producing parasite material in the laboratory have made this species a primary model of interest. For instance, specific genomic resources have been built to elucidate the genetic architecture of drug resistance mechanisms [4] and this species has also been used to evaluate alternative control strategies such as selective breeding for resistance to infection [5]. The success of these research avenues relies heavily upon the population dynamics and biology of *H. contortus* that both influence the diffusion of drug resistance alleles [6].

<footnote>* Correspondence: mathilde.saccareau@gmail.com
¹INRA - INPT-ENSAT - INPT-ENVT - Université de Toulouse, UMR 1388 GenPhySE, Castanet Tolosan, France
Full list of author information is available at the end of the article</footnote>
Mathematical models of parasite life-cycles provide an effective tool to advance epidemiological understanding [7, 8] and aid decision-making in parasite control strategies [9–12]. Indeed, the in silico approach allows for multiple scenarios to be considered without having to resort to experimentation once the model was validated [13–15]. These mechanistic models rely on the input of specific parameters, such as life history traits of the parasitic phase. Robust and average values for these traits are lacking, despite an extensive body of literature available on the subject. Recent systematic reviews and meta-analyses have been published to build knowledge of these parameters in other important gastrointestinal parasites such as *Cooperia oncophora* and *Ostertagia ostertagi* [16, 17] and *Teladorsagia circumcincta* [18]. Yet, despite the comprehensive literature available tackling the life-cycle of *H. contortus* [1], including both the free-living (eggs to infective larvae) and parasitic stages (infective larvae to reproductive adult worms), a quantitative assessment of *H. contortus* population dynamics and ecology is still missing.

The aim of this paper was to conduct a systematic review and meta-analysis of data collected on experimentally infected sheep with *H. contortus*, published before 2016, using the same methodology as Verschave et al. [16]. The main life history traits selected for study in the parasitic phase of *H. contortus* included: (i) larval establishment; (ii) adult mortality; (iii) female fertility and population level fecundity. To obtain robust estimates of the life history traits, several experimental variation factors were considered i.e. infection dose, duration of the infection, host age, genetic background of host resistance and previous exposure to *H. contortus*.

## Methods

Verschave et al. [16] performed the systematic review and meta-analysis of *Ostertagia ostertagi* life traits in bovines. In this paper, the same approach was used for *Haeemonchus contortus* life history traits. However papers reporting trickle infection results were removed as overlapping worm generations obscure worm age and thus prevent reliable estimation of larval establishment and adult mortality rates. Additional explanatory factors were studied, i.e. the influence of host resistance and of a previous exposure to *H. contortus* before the studied experiment. Some variance equations used in Verschave et al. [16] were corrected in this paper (Additional file 1: Text).

## Searching strategy and eligibility criteria

An exhaustive data collection was performed using the Web of Sciences database (https://www.webofknowledge.com) using the following keywords: ((gastrointestinal AND nematode) or (gastrointestinal AND parasite*) OR

haemonchus OR contortus)) AND (infect* OR pathogen* OR mortality OR establishment OR fecundity) AND (sheep OR lamb* OR ewe* OR wether* OR ram* OR dam*). Selection of the papers used in the present study was completed on January 7th, 2016.

In the first review phase, papers were selected based on title compatibility. Papers without full text accessibility, unavailable through library services or published in languages other than French or English were discarded. Finally, the eligibility of the remaining papers was decided following a full reading of the text to meet the following criteria: (i) natural infections were excluded due to unknown parasite exposure; (ii) experiments which used anthelmintic treatments prior to sheep slaughter were not considered; and (iii) papers reporting their results graphically i.e. without precise values were removed.

In the second review phase, additional selection criteria were applied to the papers based on experimental design conditions. Insufficiently detailed experiments (i.e. the time between infection and slaughter) or nutritional supplementation studies (e.g. protein, zinc, selenium, papaya latex, condensed tannins) were withdrawn. Studies based on trickle infection protocols were also discarded, as larval stages from consequent infections will overlap with adult stages from the initial challenge making it impossible to correctly estimate larval establishment and adult mortality rates, particularly for infections of long duration. To ensure a robust meta-analysis, only experiments reporting the arithmetic means and associated measures of variance (standard deviation or standard error of the mean) for worm burden were kept for further analysis. Studies lacking experimental power (less than 3 animals) were ignored. Finally, the data were weighed according to their associated variances during estimation procedures. As atypical behaviour of variances were obtained for underpowered experiments (less than 10 animals), experiments showing too large or too small worm burden variance estimate (i.e. two standard deviations from the mean associated variance) were discarded.

## Parameter definition and specific eligibility criteria

The life history traits of *H. contortus* considered were: larval establishment which was estimated as the proportion of ingested larvae retrieved at necropsy, irrespective of the developmental stage (i.e. *L*₄, immature stages, adult worms), the adult mortality rate (i.e. the proportion of adult worms that die per day), the female fertility and the population level fecundity. As defined by Southwood [19], fertility corresponds to the individual capacity of a female to lay eggs (measured by the total daily egg output in sheep faeces divided by the total number of female worms), whereas fecundity corresponds to the capacity of
a population to lay eggs (measured by the total daily egg output in faeces). For each considered life history trait, the equations used to calculate mean estimates and their associated variances are summarised in Table 1 and detailed in Additional file 1: Text.

To disentangle larval establishment rate from adult mortality rate, we assumed according to expert knowledge that no adult worms die before 30 days post-infection (dpi), hence estimating larval establishment rate from experiments lasting less than 30 days and adult mortality rate from longer experiments. The worm's sex ratio was estimated from experiments reporting separately female and male worm numbers. The number of female worms was assumed by this sex ratio estimation when unknown. Daily faeces production (DFP, g/day) was estimated as a direct function of host body weight (BW, kg) [20]:

$$DFP (g/day) = 0.041 * 1000 * BW^{0.75}$$

Body weight that was missing in one paper [21] using Texel sheep was inferred from the paper of Freetly et al. [22] who derived a function linking body weight and age in Texel sheep.

The general study description (country, year, number of sheep), sheep-host details (breed, gender, age, body weight, previous exposure to *H. contortus*), *H. contortus* details (isolate, drug resistance status, infection dose) and experimental conditions (infection duration, concomitant infection with other nematode species) were included in our database. Host resistance to gastrointestinal nematodes was inferred from the literature, where Barbados Black Belly, Saint Croix, Florida Native [23], Santa Ines [24, 25], Canaria hair breed [26] and divergent selection lines were considered resistant breeds. Two *H. contortus* isolates were included: anthelmintic resistant (thiabendazole-resistant strain) or *H. contortus* isolates (ISE isolates) obtained after serial passages in resistant hosts.

## Meta-analysis

Statistical analyses were performed using R version 3.1.3. [27]. The life history trait estimates and their associated

variances were calculated for each experiment of our final database (Table 1). Each life history trait estimate was analysed using a linear fixed model (stats package) weighted by the inverse of their associated variance. The considered explanatory variables were the host age (in months), the host resistance (i.e. susceptible, resistant, unknown status), the infection dose (i.e. the number of L₃ inoculated per animal), the infection duration (i.e. number of days post-infection), previous exposure to *H. contortus* (i.e. immunologically naïve, natural pre-infection assumed, experimental pre-infection). Due to the skewed distribution of the life history traits, they were studied with, and without, a log transformation in further analyses. Quantitative explanatory variables (i.e. host age, infection dose and duration of infection) were clustered into balanced classes as the hypothesis of linearity between life history traits and explanatory variables was rejected.

First, a backward variable selection was implemented based on AIC criterion [28]. Then, only factors with a *P*-value < 0.05 were kept in the final models. *F*-tests were calculated to select significant factors using the *car* R package [29]. The estimation of means for each level of factors was performed using the *lsmeans* R package [30].

As each life history trait had a different inverse variance weight, it was impossible to use the same criteria to weight the correlation between life history traits. Consequently, correlations between life history traits were weighted by the number of sheep involved in the experiment using the 'weight' R package [31].

## Results

Figure 1 provides a description of the paper selection steps and the associated selection criteria. A total of 9480 references were retrieved from Web of Sciences that were narrowed down to 383 following title-based selection. The first review phase of the papers led to a selection of 81 eligible publications, reduced with the second review phase to 37 papers that estimate at least one *H. contortus* life history trait [21, 32–67]. This final database (see Additional file 2) comprised a total of 982 animals used in 115 different experiments (using different

<table><caption>**Table 1** Mean estimate and associated variance of the main life history traits of *H. contortus*</caption><thead><tr><th></th><th>Larval establishment (E)</th><th>Adult mortality (μ)</th><th>Female fertility (F)</th><th>Population level fecundity (f)</th></tr></thead><tbody><tr><td>Definition</td><td>Proportion of ingested larvae that develop into immature or adult worms</td><td>Proportion of adult worms that die per day</td><td>Mean number of eggs laid in sheep faeces by an adult female per day</td><td>Mean number of eggs counted per gram of faeces</td></tr><tr><td>Mean estimate</td><td>WB/ID</td><td>-ln(WB/t)</td><td>FEC<sub>n</sub>*DFP/WB+F<sub>p</sub></td><td>FEC<sub>n</sub></td></tr><tr><td>Variance</td><td>SE(WB)<sup>2</sup>/ID<sup>2</sup></td><td>SE(WB)<sup>2</sup>/(t*WB)<sup>2</sup></td><td>DFP<sup>2</sup>/F<sub>p</sub><sup>2</sup> * FEC<sub>n</sub><sup>2</sup>/WB<sup>2</sup> * (SE(FEC<sub>n</sub>)<sup>2</sup>/FEC<sub>n</sub><sup>2</sup> + SE(WB)<sup>2</sup>/WB<sup>2</sup>)</td><td>SE(FEC<sub>n</sub>)</td></tr></tbody></table>

Abbreviations: WB Worm burden, ID Infection dose, t days after infection, FECn Faecal egg count at necropsy, DFP Daily faeces production, Fp, proportion of female worms, SE Standard error
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Diagram of the paper selection from the database Web of Science and associated eligibility criteria</figcaption></figure>

experimental conditions as the duration of infection and/or the infection dose).

Table 2 summarizes the details and available life history traits for each paper kept in the final database and the distribution of the experiment number over the different explanatory factors is shown in Additional file 3: Table S1.

In addition, the proportion of female worms was estimated from 28 experiments that indicated female and male worm numbers separately. The mean (± standard error, SE) proportion of female worms was 0.51 ± 0.006, which corresponds to a female: male sex ratio of 1.04:1. Estimates for other life history traits of the parasitic phase of *H. contortus* are provided in Table 3.

## Larval establishment rate

Average larval establishment rate was 0.24 ± 0.02. The final model, including all the significant explanatory factors, explained 49.5% of total variance (Table 3).

The *H. contortus* establishment rate peaked between 10 and 19 dpi (Fig. 2a). Larval establishment rate significantly decreased with the size of the infection dose from 10,000 to 20,000 L₃ (Fig. 2b), but was highest for the greatest infection dose (result not significant). A reduced proportion of larvae reached adult stages in older hosts (Fig. 2c). This result was consistently obtained, with or without experiments using older sheep than usual age, i.e. 36-month-old vs 7 months on average for remaining studies (results not shown).

In our meta-analysis, the proportion of immature stages significantly decreased after 20 dpi and almost disappeared after 40 dpi (Additional file 4: Figure S1), except within a few underpowered experiments from the same paper [47] performed at 49 and 70 dpi.

## Adult mortality rate

The average (± SE) daily per capita adult mortality rate was 0.021 ± 0.002, resulting in a mean life expectancy (1/mortality rate) of 50 days. The proportion of worms surviving (worm burden/infection dose) with
<table><caption>**Table 2** Characteristics of the studies included in the final database to estimate at least one of the key life history traits of *Haemonchus contortus*</caption><thead><tr><th>Reference</th><th>Country</th><th>No. of experiments</th><th>No. of animals</th><th>Breed</th><th>Age (months)</th><th>Immune past</th><th>Duration of infection (days)</th><th>Infection dose</th><th>No. of experiments</th></tr><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>E</th><th>μ</th><th>F</th><th>f</th></tr></thead><tbody><tr><td>[32]</td><td>Australia</td><td>6</td><td>22</td><td>Merino</td><td></td><td>Naïve: 20,000 L<sub>3</sub></td><td>7-21-24-84</td><td>5000-10,000</td><td>5</td><td>1</td><td></td><td></td></tr><tr><td>[33]</td><td>Australia</td><td>5</td><td>21</td><td>Merino</td><td>5</td><td>Naïve-Natural</td><td>7-21-24</td><td>5000-20,000</td><td>5</td><td></td><td></td><td></td></tr><tr><td>[34]</td><td>Australia</td><td>5</td><td>22</td><td>Merino</td><td>2.5-6.5</td><td>Natural: 5000 L<sub>3</sub>-6000 L<sub>3</sub></td><td>9-12-30-58</td><td>5000-6000-10,000</td><td>4</td><td>1</td><td></td><td></td></tr><tr><td>[35]</td><td>Australia</td><td>4</td><td>15</td><td>Merino</td><td></td><td>Naïve</td><td>21</td><td>10,000</td><td>4</td><td></td><td></td><td></td></tr><tr><td>[36]</td><td>Australia</td><td>2</td><td>12</td><td>Merino</td><td>8.5-36</td><td>Naïve: 20,000 L<sub>3</sub></td><td>21</td><td>10,000</td><td>2</td><td></td><td></td><td></td></tr><tr><td>[37]</td><td>Australia</td><td>2</td><td>10</td><td>Merino</td><td>9.25-14</td><td>Naïve</td><td>24-28</td><td>5000-10,000</td><td>2</td><td></td><td></td><td></td></tr><tr><td>[38]</td><td>Australia</td><td>2</td><td>9</td><td>Merino</td><td></td><td>Naïve-Natural</td><td>21-26</td><td>5000-10,000</td><td>2</td><td></td><td></td><td></td></tr><tr><td>[39]</td><td>Iraq</td><td>2</td><td>10</td><td>Awassi-Merino</td><td>5.5</td><td>Natural</td><td>39</td><td>10,600-13,300</td><td>2</td><td></td><td></td><td></td></tr><tr><td>[40]</td><td>Brazil</td><td>1</td><td>6</td><td>Santa Ines</td><td>3</td><td>Natural</td><td>40</td><td>4000</td><td>1</td><td></td><td></td><td></td></tr><tr><td>[41]</td><td>Mexico</td><td>1</td><td>4</td><td>Suffolk x Dorset</td><td>6.5</td><td>Naïve</td><td>42</td><td>10,000</td><td>1</td><td></td><td></td><td></td></tr><tr><td>[42]</td><td>Brazil</td><td>2</td><td>20</td><td>Corriedale-Crioula Lanada</td><td>3</td><td>Naïve</td><td>84</td><td>3560-4900</td><td>2</td><td></td><td></td><td></td></tr><tr><td>[43]</td><td>UK</td><td>1</td><td>5</td><td>Blueface Leicester x Scottish blackface</td><td>4.5</td><td>Naïve</td><td>35</td><td>5000</td><td>1</td><td></td><td></td><td></td></tr><tr><td>[44]</td><td>Brazil</td><td>1</td><td>7</td><td>Suffolk</td><td>3</td><td>Natural</td><td>28</td><td>4000</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>[45]</td><td>Scotland</td><td>1</td><td>3</td><td>Scottish Blackface (3/4)</td><td>8</td><td>Natural</td><td>27</td><td>5000</td><td>1</td><td></td><td></td><td></td></tr><tr><td>[46]</td><td>France</td><td>4</td><td>20</td><td>Barbados Black Belly-Romane</td><td>12</td><td>20,000 L<sub>3</sub></td><td>60</td><td>10,000</td><td>4</td><td></td><td></td><td></td></tr><tr><td>[47]</td><td>USA</td><td>23</td><td>69</td><td>Dorset</td><td>6.5</td><td>Naïve: 30,000</td><td>28-49-70</td><td>30,000</td><td>7</td><td>16</td><td>23</td><td></td></tr><tr><td>[48]</td><td>Mexico</td><td>1</td><td>6</td><td>Colombian breed</td><td>7</td><td>Naïve</td><td>51</td><td>5000</td><td>1</td><td></td><td></td><td></td></tr><tr><td>[49]</td><td>USA</td><td>1</td><td>9</td><td>Dorset</td><td>8</td><td>10,000 L<sub>3</sub></td><td>35</td><td>10,000</td><td>1</td><td></td><td></td><td></td></tr><tr><td>[50]</td><td>Mexico</td><td>1</td><td>5</td><td>Pelibuey x Dorper</td><td>3</td><td>Naïve</td><td>41</td><td>3000</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>[51]</td><td>Australia</td><td>4</td><td>22</td><td>Merino (genetically resistant) Merino</td><td>5.5</td><td>Naïve: 20,000</td><td>28-42</td><td>20,000</td><td>2</td><td>2</td><td></td><td></td></tr><tr><td>[52]</td><td>Australia</td><td>1</td><td>3</td><td>Merino</td><td>12</td><td>Natural</td><td>28</td><td>20,000</td><td>1</td><td></td><td></td><td></td></tr><tr><td>[53]</td><td>Spain</td><td>2</td><td>17</td><td>Canaria-Canaria Hair Breed</td><td>8</td><td>Natural</td><td>28</td><td>20,000</td><td>2</td><td></td><td></td><td>2</td></tr><tr><td>[54]</td><td>UK</td><td>1</td><td>4</td><td>Dorset</td><td>3</td><td>Naïve</td><td>22</td><td>5000</td><td>1</td><td></td><td></td><td></td></tr><tr><td>[55]</td><td>Germany</td><td>1</td><td>100</td><td>German Merino and 4 crossbreds</td><td>3</td><td>Naïve</td><td>49</td><td>5000</td><td>1</td><td>1</td><td></td><td></td></tr><tr><td>[56]</td><td>Hungary</td><td>3</td><td>70</td><td>Hungarian Merino</td><td>3</td><td>Naïve: 7000 L<sub>3</sub></td><td>48-85</td><td>7000</td><td>3</td><td></td><td></td><td></td></tr><tr><td>[57]</td><td>Ethiopia</td><td>2</td><td>12</td><td>Ethiopian highland sheep</td><td>11</td><td>Natural</td><td>91</td><td>4000-6000</td><td>2</td><td></td><td></td><td></td></tr><tr><td>[58]</td><td>France</td><td>2</td><td>10</td><td>Préalpes Ile de France</td><td>2.5</td><td>Naïve</td><td>73-87</td><td>10,000-50,000</td><td>2</td><td></td><td></td><td></td></tr><tr><td>[59]</td><td>USA</td><td>2</td><td>24</td><td>St Croix-Dorset x Rambouillet x Finnysheep</td><td>4.5</td><td>12,000</td><td>27</td><td>10,000</td><td>2</td><td></td><td></td><td>2</td></tr><tr><td>[60]</td><td>South Africa</td><td>1</td><td>3</td><td>Mutton Merino</td><td>6</td><td>Naïve</td><td>42</td><td>50,000</td><td>1</td><td></td><td></td><td></td></tr><tr><td>[61]</td><td>France</td><td>4</td><td>268</td><td>Romane x Barbados BlackBelly BC</td><td>3</td><td>Natural</td><td>42</td><td>5,000</td><td>4</td><td></td><td></td><td>4</td></tr></tbody></table>
<table><caption>**Table 2** Characteristics of the studies included in the final database to estimate at least one of the key life history traits of *Haemonchus contortus* (Continued)</caption><tr><td>[62]</td><td>France</td><td>2</td><td>44</td><td>Romane x Barbados BlackBelly BC</td><td>5</td><td>Náive</td><td>30</td><td>10,000</td><td>2</td><td>2</td></tr><tr><td>[21]</td><td>The Netherlands</td><td>4</td><td>18</td><td>Texel</td><td>8</td><td>Náive: 5000-10,000-40,000</td><td>28</td><td>5000</td><td>4</td><td>4</td></tr><tr><td>[63]</td><td>France</td><td>8</td><td>38</td><td>Barbados Black Belly-INRA 401</td><td>6</td><td>Náive: 10,000 L<sub>3</sub></td><td>4-30</td><td>10,000</td><td>8</td><td>4</td></tr><tr><td>[64]</td><td>France</td><td>2</td><td>8</td><td>Barbados Black Belly-INRA 401</td><td>8</td><td>Náive</td><td>16</td><td>10,000</td><td>2</td><td></td></tr><tr><td>[65]</td><td>New Zealand</td><td>2</td><td>26</td><td>Romney-Romney selected for greasy fleece</td><td>14.5</td><td>Natural</td><td>28</td><td>4000</td><td>2</td><td></td></tr><tr><td>[66]</td><td>Ethiopia</td><td>1</td><td>6</td><td>Indigenous breed</td><td>1</td><td>Náive</td><td>84</td><td>5000</td><td>1</td><td></td></tr><tr><td>[67]</td><td>USA</td><td>8</td><td>34</td><td>St Croix-Florida Native-Dorset x Rambouillet</td><td>9.5</td><td>Náive: 20,000</td><td>7-14</td><td>16,000-20,000</td><td>8</td><td></td></tr></table>
<table><caption>**Table 3** The results of the final linear models to estimate each life history trait of the parasitic phase of *H. contortus*</caption><thead><tr><th>Life history traits</th><th>Explanatory factors</th><th>Estimation of means</th><th>SE</th><th>P</th></tr></thead><tbody><tr><td colspan="5">Larval establishment rate (proportion of ingested larvae that develop into immature or adult worms)</td></tr><tr><td rowspan="3">Infection duration</td><td>4-9 dpi</td><td>0.09</td><td>0.08</td><td rowspan="3">0.04</td></tr><tr><td>12-16 dpi</td><td>0.27</td><td>0.05</td></tr><tr><td>21-30 dpi</td><td>0.23</td><td>0.03</td></tr><tr><td rowspan="4">Infection dose</td><td>4000-5000 L<sub>3</sub></td><td>0.25</td><td>0.05</td><td rowspan="4">0.004</td></tr><tr><td>10,000 L<sub>3</sub></td><td>0.22</td><td>0.03</td></tr><tr><td>16,000-20,000 L<sub>3</sub></td><td>0.05</td><td>0.07</td></tr><tr><td>30,000 L<sub>3</sub></td><td>0.26</td><td>0.08</td></tr><tr><td rowspan="3">Host age</td><td>2.5-5 months</td><td>0.29</td><td>0.06</td><td rowspan="3">0.007</td></tr><tr><td>5.5-10 months</td><td>0.21</td><td>0.02</td></tr><tr><td>12-36 months</td><td>0.08</td><td>0.07</td></tr><tr><td colspan="5">Adult mortality rate (Proportion of adult worms that die per day)</td></tr><tr><td rowspan="3">Host GIN resistance status</td><td>Susceptible</td><td>0.039</td><td>0.009</td><td rowspan="3">0.003</td></tr><tr><td>Resistant</td><td>0.065</td><td>0.011</td></tr><tr><td>Unknown</td><td>0.048</td><td>0.009</td></tr><tr><td rowspan="4">Isolate status</td><td>Susceptible to anthelmintics or non-adapted to sheep resistance isolate of H. contortus</td><td>0.044</td><td>0.009</td><td rowspan="4">0.005</td></tr><tr><td>Isolate adapted to sheep resistance<sup>a</sup></td><td>0.069</td><td>0.033</td></tr><tr><td>Isolate resistant to anthelmintics</td><td>0.062</td><td>0.012</td></tr><tr><td>Unknown</td><td>0.026</td><td>0.003</td></tr><tr><td rowspan="3">Host age</td><td>1.2-3 months</td><td>0.045</td><td>0.009</td><td rowspan="3">0.005</td></tr><tr><td>4.5-8 months</td><td>0.059</td><td>0.010</td></tr><tr><td>11-12 months</td><td>0.047</td><td>0.009</td></tr><tr><td rowspan="3">Previous host exposure</td><td>Immunologically naive</td><td>0.049</td><td>0.009</td><td rowspan="3">&lt; 0.001</td></tr><tr><td>Suspected natural pre-infection</td><td>0.039</td><td>0.010</td></tr><tr><td>Experimental abbreviated pre-infection</td><td>0.063</td><td>0.009</td></tr><tr><td colspan="5">Female fertility (number of eggs/female/day)</td></tr><tr><td rowspan="3">Host GIN resistance status</td><td>Susceptible</td><td>4545.2</td><td>310.0</td><td rowspan="3">&lt; 0.001</td></tr><tr><td>Resistant</td><td>2740.9</td><td>361.5</td></tr><tr><td>Unknown</td><td>3465.6</td><td>513.1</td></tr><tr><td rowspan="4">Infection duration</td><td>27-30 dpi</td><td>2136.5</td><td>308.9</td><td rowspan="4">&lt; 0.001</td></tr><tr><td>41-50 dpi</td><td>2270.6</td><td>204.3</td></tr><tr><td>70 dpi</td><td>6344.7</td><td>821.2</td></tr></tbody></table>
<table><caption>**Table 3** The results of the final linear models to estimate each life history trait of the parasitic phase of *H. contortus* (Continued)</caption><tbody><tr><td>Previous exposure of host</td><td></td><td></td><td>&lt; 0.001</td></tr><tr><td>Immunologically naive host</td><td>4626.8</td><td>372.0</td><td></td></tr><tr><td>Experimental pre-infection</td><td>2541.0</td><td>377.8</td><td></td></tr><tr><td>Population level fecundity (number of eggs/gram of faeces)<sup>b</sup></td><td></td><td></td><td></td></tr><tr><td>Host GIN resistance status</td><td></td><td></td><td>&lt; 0.001</td></tr><tr><td>Susceptible</td><td>13536.8</td><td>1429.1</td><td></td></tr><tr><td>Resistant</td><td>4837.4</td><td>918.8</td><td></td></tr><tr><td>Unknown</td><td>9563.8</td><td>992.7</td><td></td></tr><tr><td>Infection duration</td><td></td><td></td><td>0.002</td></tr><tr><td>27–28 dpi</td><td>5032.6</td><td>1034.9</td><td></td></tr><tr><td>30 dpi</td><td>18032.4</td><td>2575.9</td><td></td></tr><tr><td>41–49 dpi</td><td>4873.0</td><td>415.6</td><td></td></tr></tbody></table>

time was plotted for every experiment (Fig. 3). As expected, an exponential decrease of the survival probability was observed. The explanatory factors selected as significant in the final model (Table 3) explained 64.6% of the heterogeneity in adult mortality. Strong discrepancies in study's contribution were observed (Fig. 4).

Significant variations in adult mortality rates were observed according to the *H. contortus* isolate used in the experiment. Susceptible isolates had a significantly higher adult mortality rate than isolates with an unknown

resistance. Although not significant, anthelmintic resistant isolates (thiabendazole-resistant strain) demonstrated higher estimated adult mortality rates, and isolates resistant to the immune response had the highest adult mortality rate estimate.

Other explanatory factors were host resistance ($F_{(2,37)} = 6.82, P = 0.003$), previous exposure of the host to *H. contortus* ($F_{(2,37)} = 11.385, P < 0.001$) and host age ($F_{(2,37)} = 6.039, P = 0.005$). Resistant sheep and sheep previously infected with *H. contortus* experimentally eliminated worm burdens faster than their susceptible and immunologically naive

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

<footnote><sup>a</sup>By serial passages in resistant host</footnote>

<footnote><sup>b</sup>The host age was not included in the models for this trait due to a limited dataset</footnote>

<footnote>Abbreviation: SE standard error</footnote>
<figure><img src="image_5.png" /></figure>

counterparts. Adult mortality rates increased with host age (except for individuals aged from 8 to 12 months old, result not significant). In this class, two experiments were high-weighted and decreased the adult mortality rate estimation by erasing the trend of the other experiments.

## Female fertility and fecundity

The average (± SE) female fertility was 1295.9 ± 280.4 eggs per female worm per day. It was strongly associated

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

with the duration of infection ($F_{(2,31)} = 12.343, P < 0.001$): female fertility increased with the duration of infection (Fig. 5). Female fertility was less in previously challenged hosts and in resistant breeds compared to their immunologically naive counterparts and susceptible breeds. The explanatory factors considered in the model explained 85.11% of the heterogeneity.

The average faecal egg count at necropsy was 908.5 ± 487.1 eggs per gram of faeces. This fecundity was also affected by the host resistance status ($F_{(2,9)} = 72.341, P < 0.001$) and by the duration of infection ($F_{(2,9)} = 13.171, P = 0.002$). Fecundity strongly decreased from 30 dpi (Additional file 5: Figure S2) whereas the female fertility continued to increase until 70 dpi (Fig. 5).

## Phenotypic correlations between life history traits

As shown on Table 4, the adult mortality rate was negatively correlated with female fertility ($r_{(22)} = -0.49, P < 0.05$). Sex ratio was significantly correlated with all other life history traits: negatively correlated with both larval establishment rate ($r_{(12)} = -0.77, P < 0.01$) and female worm fertility ($r_{(5)} = -0.9, P < 0.05$) and positively correlated with adult mortality rate ($r_{(12)} = 0.69, P < 0.05$).

## Discussion

To our knowledge, this meta-analysis is the first to collect and summarize the available information on the main life history traits of the parasitic phase of *H. con-tortus*. The quantification of each parameter, associated with a variation, provides a good foundation to calibrate

<figure><img src="image_7.png" /></figure>
<table><caption>**Table 4** Pearson's correlation coefficients between the life history traits of the parasitic phase of *H. contortus* (weighted by the number of animals in the experiment)</caption><thead><tr><th></th><th>Larval establishment rate (E)</th><th>Adult mortality rate (μ)</th><th>Female fertility (F)</th><th>Population level fecundity (f)</th><th>Sex ratio</th></tr></thead><tbody><tr><td>E</td><td>ne</td><td></td><td>-0.07</td><td>0.39</td><td>-0.77**</td></tr><tr><td>μ</td><td></td><td></td><td>-0.49*</td><td>-0.88*</td><td>0.69*</td></tr><tr><td>F</td><td></td><td></td><td></td><td>0.58*</td><td>-0.9*</td></tr><tr><td>f</td><td></td><td></td><td></td><td></td><td>-0.86</td></tr><tr><td>Sex ratio</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

future epidemiological models describing the parasitic phase of *H. contortus*.

To ensure the quality of our estimates, several criteria were integrated. First, we selected papers giving arithmetic means and variance for two traits: worm burden and faecal egg counts. For the few papers, which gave individual measurements for these two traits, aggregated values were calculated (arithmetic mean and variance). This conservative strategy allowed for increased information from a greater number of studies to be included. This is in contrast to a similar meta-analysis carried out for *T. circumcincta*, which only considered individual worm burden measurements [18]. Like many parametric tests, the analysis of variance assumed the data were normally distributed. However, the traits studied here were non-normal and presented a skewed distribution. Fortunately, the ANOVA test is not very sensitive to moderate deviations from normality; simulation studies, using a variety of non-normal distributions have shown that the false positive rate is not much affected by not adhering to this assumption [68]. In our case, a log-transformation was applied to normalize the life history traits; the results were the same with and without this transformation confirming the robustness of ANOVA. Secondly, using the inverse variance of each life history trait to weight the mean estimate is the most robust way to analyse data reported as a mean with measurement standard deviation and animal number [69]. Studies which only reported means (without the associated measure of variance) or graphical values had to be excluded from further analysis.

Larval establishment and adult mortality are two interwoven traits. Indeed, while most larval stages will evolve into mature adults within 15 days [70], delays can occur that result in larval recovery for up to 30 days after infection [63]. To provide the best larval establishment rate estimates, i.e. the total proportion of infective larvae reaching the adult stage, worm count data up to 30 dpi were considered. Usually, *H. contortus* egg output increases between this stage of infection and 30 dpi as seen in longitudinal faecal egg count sampling (Additional file 5: Figure S2), suggesting that either

female worm fertility increases, or that the total number of mature adult stages increases, or both. Therefore, we considered that adult mortality between 15 and 30 dpi was negligible and not to be taken into account. To calculate the number of female worms, a sex ratio of 0.5 was assumed in the material and methods and the present meta-analysis corroborates this sex ratio. This value is a slightly lower than estimations from another meta-analysis which considered several nematode species simultaneously [71].

The success of the nematode in terms of fitness, based on the ability to establish and to produce offspring, relies on both parasite density-dependent constraints and host immune response [72]. Density-dependent stunting of female fertility has been reported for various parasitic species like *H. contortus* in sheep [73], *O. ostertagi* in cattle [16], *Trichostrongylus retortaeformis* in rabbits [72], *Strongyloides robustus* in squirrels [74] or *Syngamus trachea* in crows [75]. From our data, the infection dose did not significantly impact *H. contortus* female fertility, nor its adult mortality rate, but it affected the larval establishment rate. This unexpected result may be explained by the homogenous and relatively small infection doses used in our database to estimate female fertility and adult mortality which were mainly between 5000 L₃ and 10,000 L₃ and thus lower than the 10,000 L₃ threshold where a density dependence effect is generally expected [73]. The only experiments with an infection dose of 30,000 L₃ used a small sample size (three animals) leading to highly variable estimations. In addition, all these experiments came from the same paper ([47]) with specific environmental conditions (i.e. parasite strain, sheep breed, body weight, sheep sex). In natural infections, the cumulative level of infection is probably lower than 10,000 L₃ larvae (J. Cabaret, personal communication). The results gathered in our study have confirmed the negative impact of the host immune status on both *H. contortus* fertility and mortality, as reported for other trichostrongylid species [76]. These findings tend to favor a negative effect of steric crowding on larval establishment. It remains unresolved whether the observed effects on *H. contortus* mortality and

<footnote>*P < 0.05; **P < 0.01</footnote>

<footnote>Abbreviation: ne No estimate (less than two experiments)</footnote>
fertility were due to direct host responses, or from density-dependent effects mediated by the host immune response as demonstrated in *Strongyloides ratti* infections [77]. Further, the infection dose is only an indirect approximation of the actual established worm burden so that density dependence may be overlooked.

The negative impact of host genetic resistance on *H. contortus* mortality rate and female fertility corroborates previous results comparing resistant and susceptible breeds [78, 79] or lines of sheep with divergent resistance status [62, 80]. Such an observation also links with the proposed framework for the immune response against *T. circumcincta* [81].

Surprisingly, naturally pre-challenged sheep seemed to exert a significantly lighter response on adult mortality rate than their experimentally pre-challenged or naïve counterparts. This may be a result of the misclassification of studies due to a lack of information in the material and methods section. For example, some experimental animals were explicitly mentioned as naturally infected (even if pasture infectivity was not stated) and treated to remove presence of parasites, but some others were mentioned as receiving preventative treatments but not really faced parasite challenge.

*Haemonchus contortus* is well known as one of the most prolific sheep parasite [82]. Our estimates indicate a high female fertility of 1295.9 ± 280.4 eggs/female/day on average, which outperforms other trichostrongylids [82]. It is usually postulated that this high female fertility will result in a huge population size showing a high degree of genetic variability [83]. However, our estimates also suggest that larval establishment constitutes a strong bottleneck as only a quarter of the ingested larvae will have an opportunity to mate and pass on their genetic material to subsequent generations. This may constrain the expected population size of *H. contortus*.

The mean estimated life expectancy of 50 days by this meta-analysis underestimated the predicted half-life of an established worm burden of 69 days as described by Barger & Le Jambre [84]. Indeed, the equation used to estimate the adult mortality rate [47] did not take into account the establishment of the ingested larvae (assumed as 100%) leading to an overestimation of the adult mortality rate. The analyses were also made with the equation of Barger & Le Jambre [84] using the average larval establishment rate estimated by our meta-analysis. Unfortunately, for a third of the experiments this estimate is higher than the effective larval establishment rate of the experiment, leading to a negative adult mortality rate (replaced by 0 in this case) and thus biasing the adult mortality estimates.

Interestingly, the more female-biased *H. contortus* population was, the lower the establishment and fertility rates were, yet with higher adult mortality rates. These

correlations are difficult to interpret as the higher proportion of females could be a cause, or a consequence, of the two other traits. Female worms are thought to live longer than males, which could explain why populations exhibiting high adult mortality rate also show more females than males [71]. The higher proportion of females in low larval establishment rates may be a result of a better female survival capacity. The decrease in female fertility observed for female-biased populations may be due to either a competition for males, or competition for resources as their nutrient requirements are usually higher than for males [71].

## Conclusions

To our knowledge, this comprehensive systematic review of the parasitic phase of *H. contortus* is the first to summarize these extensive data and to provide average overall life history trait estimates. These informative estimations are very helpful to parameterize epidemiological transmission models more accurately. This meta-analysis also gives evidence for density dependence of larval establishment, for host age affecting larval establishment and adult mortality and shows that an improved immune response, induced by the resistance status of the breed or by a previous exposure to the parasites, affects adult mortality and female fertility. An improved understanding of the parasitic life-cycle will allow us to evaluate the impact of different control strategies on parasitic infection with an increased confidence in the output predictions of dynamic models. For example, the sustainability of treatments and genetic selection strategies could be evaluated both on the infection level of the flock and on the increase of resistant alleles in parasite strains. The impact of environmental factors on free living stages also lacks precise estimates and a meta-analysis of these would certainly help to conclude on this point.

## Additional files

**Additional file 1:** Text: Calculation method to estimate the variance of each life history trait. (DOCX 30 kb)

**Additional file 2:** Data: Final database of all the papers included in the meta-analysis after all the selection criteria. (CSV 13 kb)

**Additional file 3: Table S1.** Distribution of the experiment number, of the animal number and of the experiment weight over different levels of the explanatory factors. (DOCX 17 kb)

**Additional file 4: Figure S1.** Proportion of immature stages on the total worm burden according to the infection duration post infection. Point size represents the weight associated to each experiment. The solid black line represents the fitting by polynomial model (degree 2) and the dashed lines represent its confidence interval. (TIFF 2774 kb)

**Additional file 5: Figure S2.** Faecal egg count according to the infection duration post infection. Black cross sizes represent the weight associated to each faecal egg count reported. Grey crosses represent faecal egg count reported without the associated measure of variance. (TIFF 2600 kb)
**Abbreviations**

DFP: Daily faeces production; FEC: Faecal egg count; FEC~n~: Faecal egg count at necropsy; F~p~: Proportion of female worms; ID: Infection dose; SE: Standard error; WB: Worm Burden.

**Acknowledgments**

We are grateful to Dr Caroline Chylinski for helpful comments on the manuscript.

**Funding**

This project was funded by the INRA program GISA. M. Saccareau receives a grant from the INRA Animal Health and Animal Genetics departments.

**Availability of data and materials**

The data supporting the conclusions of this article are included within the article and its additional files.

**Authors' contributions**

MS performed literature review, statistical analysis and the main redaction of the paper. CM and CR-G participated to the choice of statistical methods. All other co-authors participated to the interpretations of results, discussions and corrections of the redaction. All authors read and approved the final manuscript.

**Competing interests**

The authors declare that they have no competing interests.

**Consent for publication**

Not applicable.

**Ethics approval and consent to participate**

Not applicable.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Author details**

¹INRA - INPT-ENSAT - INPT-ENVT - Université de Toulouse, UMR 1388 GenPhySE, Castanet Tolosan, France. ²INRA - Université de Tours, UMR 1282 ISP, 37380, Nouzilly, France. ³INRA - INPT-ENVT - Université de Toulouse, UMR 1225 IHAP, Toulouse, France. ⁴INRA - INPT-ENVT, UMT Santé des petits ruminants, Centre de recherche de Toulouse, Toulouse, France.

Received: 23 September 2016 Accepted: 5 April 2017

Published online: 24 April 2017

**References**

1. Besier RB, Kahn LP, Sargison ND, Van Wyk JA. Chapter Four-The pathophysiology, ecology and epidemiology of *Haemonchus contortus* infection in small ruminants. Adv Parasitol. 2016;93:95-143.

2. Kaplan RM, Vidyashankar AN. An inconvenient truth: global worming and anthelmintic resistance. Vet Parasitol. 2012;186:70-8.

3. Blake N, Coles G. Flock cull due to anthelmintic-resistant nematodes. Vet Rec. 2007;161:36-36.

4. Laing R, Kikuchi T, Martinelli A, Tsai IJ, Beech RN, Redman E, et al. The genome and transcriptome of *Haemonchus contortus*, a key model parasite for drug and vaccine discovery. Genome Biol. 2013;14:1-16.

5. Bishop SC, Morris CA. Genetics of disease resistance in sheep and goats. Small Rumin Res. 2007;70:48-59.

6. Cornelius MP, Jacobson C, Dobson R, Besier RB. Computer modelling of anthelmintic resistance and worm control outcomes for refugia-based nematode control strategies in Merino ewes in Western Australia. Vet Parasitol. 2016;220:59-66.

7. Callinan APL, Morley FHW, Arundel JH, White DH. A model of the life cycle of sheep nematodes and the epidemiology of nematodiasis in sheep. Agric Syst. 1982;9:199-225.

8. Barger IA. Genetic resistance of hosts and its influence on epidemiology. Vet Parasitol. 1989;32:21-35. Spec. Issue 13th WAAVP Conf.

9. Saccareau M, Moreno CR, Kyriazakis I, Faivre R, Bishop SC. Modelling gastrointestinal parasitism infection in a sheep flock over two reproductive seasons: *in silico* exploration and sensitivity analysis. Parasitology. 2016;143: 1509-31.

10. Learmount J, Taylor MA, Smith G, Morgan C. A computer model to simulate control of parasitic gastroenteritis in sheep on UK farms. Vet Parasitol. 2006; 142:312-29.

11. Barnes EH, Dobson RJ, Barger IA. Worm control and anthelmintic resistance: adventures with a model. Parasitol Today. 1995;11:56-63.

12. Leathwick DM, Vlassoff A, Barlow ND. A model for nematodiasis in New Zealand lambs: The effect of drenching regime and grazing management on the development of anthelmintic resistance. Int J Parasitol. 1995;25: 1479-90.

13. Laurenson YCSM, Kyriazakis I, Forbes AB, Bishop SC. Exploration of the epidemiological consequences of resistance to gastro-intestinal parasitism and grazing management of sheep through a mathematical model. Vet Parasitol. 2012;189:238-49.

14. Laurenson YCSM, Kyriazakis I, Bishop SC. *In silico* exploration of the impact of pasture larvae contamination and anthelmintic treatment on genetic parameter estimates for parasite resistance in grazing sheep. J Anim Sci. 2012;90:2167-80.

15. Vagenas D, Doeschl-Wilson A, Bishop SC, Kyriazakis I. *In silico* exploration of the effects of host genotype and nutrition on the genetic parameters of lambs challenged with gastrointestinal parasites. Int J Parasitol. 2007;37: 1617-30.

16. Verschave SH, Vercruysse J, Claerebout E, Rose H, Morgan ER, Charlier J. The parasitic phase of *Ostertagia ostertagi*: quantification of the main life history traits through systematic review and meta-analysis. Int J Parasitol. 2014;44:1091-104.

17. Verschave SH, Rose H, Morgan ER, Claerebout E, Vercruysse J, Charlier J. Modelling *Cooperia oncophora*: Quantification of key parameters in the parasitic phase. Vet Parasitol. 2016;223:111-4.

18. Gaba S, Gruner L, Cabaret J. The establishment rate of a sheep nematode: Revisiting classics using a meta-analysis of 87 experiments. Vet Parasitol. 2006;140:302-11.

19. Southwood TRE. Observational and experimental methods for the estimation of natality, mortality and dispersal. In: Southwood TRE, editor. Ecological methods: With Particular Reference to the Study of Insect Populations. Netherlands: Springer; 1978. p. 302-55.

20. Cabaret J, Ouhelli, H. La fertilite des strongles digestifs lors de l'infestation naturelle des brebis. Rev Med Vet. 1984;135:627-33.

21. Schallig HDF, van der Aar W, Boersema J, Cornelissen AWC. The effect of oxfendazole terminated infections with *Haemonchus contortus* on the development of immunity in sheep. Vet Parasitol. 2000;88:61-72.

22. Freely HC, Nienaber JA, Leymaster KA, Jenkins TG. Relationships among heat production, body weight, and age in Suffolk and Texel ewes. J Anim Sci. 1995;73:1030-7.

23. Bishop SC. Possibilities to breed for resistance to nematode parasite infections in small ruminants in tropical production systems. Animal. 2012;6:741-7.

24. Amarante AFT, Susin I, Rocha RA, Silva MB, Mendes CQ, Pires AV. Resistance of Santa Ines and crossbred ewes to naturally acquired gastrointestinal nematode infections. Vet Parasitol. 2009;165:273-80.

25. do Amarante A, Bricarello PA, Rocha RA, Gennari SM. Resistance of Santa Ines, Suffolk and Ile de France sheep to naturally acquired gastrointestinal nematode infections. Vet Parasitol. 2004;120:91-106.

26. Guo Z, González JF, Hernandez JN, McNeilly TN, Corripio-Miyar Y, Frew D, et al. Possible mechanisms of host resistance to *Haemonchus contortus* infection in sheep breeds native to the Canary Islands. Sci Rep. 2016;6:26200.

27. R: A Language and Environment for Statistical Computing. Vienna: R Core Team; 2014. Available from: http://www.R-project.org/. Accessed 10 Apr 2017.

28. Ripley B, Venables B, Bates DM, Hornik K, Gebhardt A, Firth D, et al. Package "MASS." CRAN Repos. Httpcran R-Proj. OrgwebpackagesMASSMASS Pdf. 2013.

29. Fox J, Weisberg S, Adler D, Bates D, Baud-Bovy G, Ellison S, et al. Package "car". 2016. https://cran.r-project.org/package=car. Accessed 10 Apr 2017.

30. Lenth R. Least-Squares means: The R Package Ismeans. J Stat Softw. 2016; doi:10.18637/jss.v069.i01

31. Pasek J. Package "weights". 2016. https://cran.r-project.org/package=weights. Accessed 10 Apr 2017.
32. Adams DB. Time of onset and target of immune reactions in sheep with acquired immunity against *Haemonchus contortus*. Int J Parasitol. 1982;12:439–43.

33. Adams DB. Investigation with dexamethasone of the processes which moderate immunity against the nematode, *Haemonchus contortus*, in sheep. Aust J Exp Biol Med Sci. 1983;61:345–53.

34. Adams DB. Observations on the self-cure reaction and other forms of immunological responsiveness against *Haemonchus contortus* in sheep. Int J Parasitol. 1983;13:571–8.

35. Adams DB. Developmental arrest of *Haemonchus contortus* in sheep treated with a corticosteroid. Int J Parasitol. 1986;16:659–64.

36. Adams DB. The effect of dexamethasone on a single and a superimposed infection with *Haemonchus contortus* in sheep. Int J Parasitol. 1988;18:575–9.

37. Adams DB. A preliminary evaluation of factors affecting an experimental system for vaccination-and-challenge with *Haemonchus contortus* in sheep. Int J Parasitol. 1989;19:169–75.

38. Adams DB, Davies HI. Enhanced resistance to infection with *Haemonchus contortus* in sheep treated with a corticosteroid. Int J Parasitol. 1982;12:523–9.

39. Khshali MN Al-, Altaif KI. The response of Awassi and Merino sheep to primary infection with *Haemonchus contortus*. Trop Anim Health Prod. 1979; 11:164–70.

40. Almeida FA, Garcia KCOD, Torgerson PR, Amarante AFT. Multiple resistance to anthelmintics by *Haemonchus contortus* and *Trichostrongylus colubriformis* in sheep in Brazil. Parasitol Int. 2010;59:622–5.

41. Bautista-Garfias C, Flores-Hernandez O, Quiroz-Romero H. Non-specific resistance of sheep against *Haemonchus* conforms with Freund's complete adjuvant. Parasite Immunol. 1991;13:565–9.

42. Bricarello PA, Gennari SM, Oliveira-Sequeira TCG, Vaz CMSL, de Gonçalves IG, Echevarria FAM. Response of Corriedale and Crioula Lanada sheep to artificial primary infection with *Haemonchus contortus*. Vet Res Commun. 2002;26:447–57.

43. Buttle DJ, Behnke JM, Bartley Y, Elsheikha HM, Bartley DJ, Garnett MC, et al. Oral dosing with papaya latex is an effective anthelmintic treatment for sheep infected with *Haemonchus contortus*. Parasit Vectors. 2011;4:36.

44. Carvalho N, Das Neves JH, Nazato C, Louvandini H, Amarante AFT. The effects of diet and corticosteroid-induced immune suppression during infection by Haemonchus contortus in lambs. Vet Parasitol. 2015;214: 289–94.

45. Christie MG, Angus KW, Hotson IK. Manifestations of resistance to *Haemonchus contortus* in sheep: Worm populations and abomasal changes in sheep superinfected with 1,000,000 larvae of *H. contortus*. Int J Parasitol. 1975;5:193–8.

46. Chylinski C. What makes a gastrointestinal nematode successful in their sheep host? Exploring the role of the nematode, the sheep host and the farmer. PhD thesis. Tours: l'Université François - Rabelais de Tours; 2014. http://www.theses.fr/2014TOUR4028.

47. Coyne MJ, Smith G. The mortality and fecundity of *Haemonchus contortus* in parasite-naive and parasite-exposed sheep following single experimental infections. Int J Parasitol. 1992;22:315–25.

48. Cuenca-Verde C, Buendía-Jiménez JA, Valdivia-Anda G, Cuéllar-Ordaz JA, Muñoz-Guzmán MA, Alba-Hurtado F. Decrease in establishment of *Haemonchus contortus* caused by inoculation of a *Taenia hydatigena* larvae vesicular concentrate. Vet Parasitol. 2011;177:332–8.

49. De Wolf BM, Zajac AM, Hoffer KA, Sartini BL, Bowdridge S, LaRoith T, et al. The effect of vitamin E supplementation on an experimental *Haemonchus contortus* infection in lambs. Vet Parasitol. 2014;205:140–9.

50. Galicia-Aguilar HH, Rodríguez-González LA, Capetillo-Leal CM, Cámara-Sarmiento R, Aguilar-Caballero AJ, Sandoval-Castro CA, et al. Effects of *Havardia albicans* supplementation on feed consumption and dry matter digestibility of sheep and the biology of *Haemonchus contortus*. Anim Feed Sci Technol. 2012;176:178–84.

51. Gill HS. Genetic control of acquired resistance to haemonchosis in Merino lambs. Parasite Immunol. 1991;13:617–28.

52. Gill HS, Husband AJ, Watson DL. Localization of immunoglobulin-containing cells in the abomasum of sheep following infection with *Haemonchus contortus*. Vet Immunol Immunopathol. 1992;31:179–87.

53. González JF, Hernández Á, Molina JM, Fernández A, Raadsma HW, Meeusen ENT, et al. Comparative experimental *Haemonchus contortus* infection of two sheep breeds native to the Canary Islands. Vet Parasitol. 2008;153:374–8.

54. Hunter AR, Mackenzie G, Ross JG. The effects of transfer factor, parasitic antigen and precipitated immunoglobulin treatments in lambs infected with *Haemonchus contortus*. Vet Parasitol. 1983;13:45–53.

55. Idris A, Moors E, Budnick C, Herrmann A, Erhardt G, Gauly M. Is the establishment rate and fecundity of *Haemonchus contortus* related to body or abomasal measurements in sheep? Animal. 2011;5:1276–82.

56. Kassai T, Fésüs L, Hendrikx WML, Takáts C, Fok É, Redl P, et al. Is there a relationship between haemoglobin genotype and the innate resistance to experimental *Haemonchus contortus* infection in Merino lambs? Vet Parasitol. 1990;37:61–77.

57. Kelkele FA, Tolossa YH, Kassa GM. Experimental infection of Ethiopian highland sheep by different infective doses of *Haemonchus contortus* (L 3): haematological and parasitological parameters, serum protein concentrations and clinical responses. Ethiop Vet J. 2012;16:41–57.

58. Kerboeuf D. Changes in pepsinogen, proteins and lipids in the serum during experimental haemonchosis in sheep. Ann Rech Vét. 1977;8:257–66.

59. MacKinnon KM, Zajac AM, Kooymان FNJ, Notter DR. Differences in immune parameters are associated with resistance to *Haemonchus contortus* in Caribbean hair sheep. Parasite Immunol. 2010;32:484–93.

60. Owen NC. Some nutritional aspects of haemonchosis in experimentally infected lambs. South Afr J Anim Sci. 1973;3:1–4.

61. Salle G, Jacquiet P, Gruner L, Cortet J, Sauve C, Prevot F, et al. A genome scan for QTL affecting resistance to *Haemonchus contortus* in sheep. J Anim Sci. 2012;90:4690–705.

62. Salle G, Moreno C, Boitard S, Ruesche J, Tircazes-Secula A, Bouvier F, et al. Functional investigation of a QTL affecting resistance to *Haemonchus contortus* in sheep. Vet Res. 2014;45:68.

63. Terefe G, Lacroux C, Andreoletti O, Grisez C, Prevot F, Bergeaud JP, et al. Immune response to *Haemonchus contortus* infection in susceptible (INRA 401) and resistant (Barbados Black Belly) breeds of lambs. Parasite Immunol. 2007;29:415–24.

64. Terefe G, Lacroux C, Prévot F, Grisez C, Bergeaud JP, Bleuart C, et al. Eosinophils in *Haemonchus contortus*-infected resistant and susceptible breeds of sheep: Abomasal tissue recruitment and in vitro functional state. Vet Parasitol. 2009;165:161–4.

65. Williamson JF, Blair HT, Garrick DJ, Pomroy WE, Douch PGC. The relationship between internal parasite burden, faecal egg count and mucosal mast cells in fleeceweight-selected and control sheep. Proc NZ Soc Anim Prod. 1994; 54:9–14.

66. Yacob HT, Mistre C, Adem AH, Basu AK. Parasitological and clinical responses of lambs experimentally infected with *Haemonchus contortus* (L3) with and without ivermectin treatment. Vet Parasitol. 2009;166:119–23.

67. Zajac AM, Krakowka S, Herd RP, McClure KE. Experimental *Haemonchus contortus* infection in three breeds of sheep. Vet Parasitol. 1990;36:221–35.

68. Glass GV, Peckham PD, Sanders JR. Consequences of failure to meet assumptions underlying the fixed effects analyses of variance and covariance. Rev Educ Res. 1972;42:237–88.

69. Jones DR. Meta-analysis: weighing the evidence. Stat Med. 1995;14:137–49.

70. Veglia F. The Anatomy and life history of the *Haemonchus contortus* (Rud). In: Third and Fourth Reports of the Director of Veterinary Research. Pretoria: Government Printer and Stationery Office; 1915. p. 347–500.

71. Poulin R. Population abundance and sex ratio in dioecious helminth parasites. Oecologia. 1997;111:375–80.

72. Chylinski C, Boag B, Stear MJ, Cattadori IM. Effects of host characteristics and parasite intensity on growth and fecundity of *Trichostrongylus retortaeformis* infections in rabbits. Parasitology. 2009;136:117–23.

73. Fleming MW. Size of inoculum dose regulates in part worm burdens, fecundity, and lengths in ovine *Haemonchus contortus* infections. J Parasitol. 1988;6:975–8.

74. Romeo C, Wauters LA, Cauchie S, Martinoli A, Matthysen E, Saino N, et al. Faecal egg counts from field experiment reveal density dependence in helminth fecundity: *Strongyloides robustus* infecting grey squirrels (*Sciurus carolinensis*). Parasitol Res. 2014;113:3403–8.

75. Gethings OJ, Sage RB, Leather SR. Density-dependent regulation of fecundity in *Syngamus trachea* infrapopulations in semi-naturally occurring ring-necked pheasants (*Phasianus colchicus*) and wild Carrion Crows (*Corvus corone*). Parasitology. 2016;143:716–22.

76. Stear MJ, Bishop SC, Doligalska M, Duncan JL, Holmes PH, Irvine J, et al. Regulation of egg production, worm burden, worm length and worm fecundity by host responses in sheep infected with *Ostertagia circumcincta*. Parasite Immunol. 1995;17:643–52.
77. Paterson S, Viney ME. Host immune responses are necessary for density dependence in nematode infections. Parasitology. 2002;125:283-92.

78. Aumont G, Gruner L, Hostache G. Comparison of the resistance to sympatric and allopatric isolates of *Haemonchus contortus* of Black Belly sheep in Guadeloupe (FWI) and of INRA 401 sheep in France. Vet Parasitol. 2003;116:139-50.

79. Gruner L, Aumont G, Getachew T, Brunel JC, Pery C, Cognie Y, et al. Experimental infection of Black Belly and INRA 401 straight and crossbred sheep with trichostrongyle nematode parasites. Vet Parasitol. 2003;116:239-49.

80. McRae KM, Good B, Hanrahan JP, Glynn A, O'Connell MJ, Keane OM. Response to *Teladorsagia circumcincta* infection in Scottish Blackface lambs with divergent phenotypes for nematode resistance. Vet Parasitol. 2014;206:200-7.

81. Stear MJ, Bairden K, Duncan JL, Holmes PH, McKellar QA, Park M, et al. How hosts control worms. Nature. 1997;389:27-27.

82. Cabaret J, Gasnier N, Jacquiet P. Faecal egg counts are representative of digestive-tract strongyle worm burdens in sheep and goats. Parasite. 1998;5:137-42.

83. Gilleard JS, Redman E. Genetic diversity and population structure of *Haemonchus contortus*. Adv Parasitol. 2016;93:31-68.

84. Barger IA, Le Jambre LF. Regulation of *Haemonchus contortus* populations in sheep: Mortality of established worms. Int J Parasitol. 1988;18:269-73.

Submit your next manuscript to BioMed Central
and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central