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

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. Electronic supplementary material The online version of this article (doi:10.1186/s13071-017-2131-7) contains supplementary material, which is available to authorized users.


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].
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][10][11][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][14][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 Haemonchus 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 7 th , 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 4 , 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 postinfection (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 3 inoculated per animal), the infection duration (i.e. number of days post-infection), previous exposure to H. contortus (i.e. immunologically naïve, natural preinfection 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]. 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,. This final database (see Additional file 2) comprised a total of 982 animals used in 115 different experiments (using different

Results
Abbreviations: WB Worm burden, ID Infection dose, t days after infection, FEC n Faecal egg count at necropsy, DFP Daily faeces production, F p , proportion of female worms, SE Standard error 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 3 (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     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  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 highweighted 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 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.

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. contortus. The quantification of each parameter, associated with a variation, provides a good foundation to calibrate 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 logtransformation 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 3 and 10,000 L 3 and thus lower than the 10,000 L 3 threshold where a density dependence effect is generally expected [73]. The only experiments with an infection dose of 30,000 L 3 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 3 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 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 the 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.