Skip to main content

Functional investigation of a QTL affecting resistance to Haemonchus contortus in sheep


This study reports a functional characterization of a limited segment (QTL) of sheep chromosome 12 associated with resistance to the abomasal nematode Haemonchus contortus. The first objective was to validate the identified QTL through the comparison of genetically susceptible (N) and resistant (R) sheep produced from Martinik × Romane back-cross sheep. The R and N genotype groups were then experimentally infected with 10 000 H. contortus larvae and measured for FEC (every three days from 18 to 30 days post-challenge), haematocrit, worm burden and fertility. Significant differences in FEC and haematocrit drop were found between R and N sheep. In addition, the female worms recovered from R sheep were less fecund. The second step of the characterization was to investigate functional mechanisms associated with the QTL, thanks to a gene expression analysis performed on the abomasal mucosa and the abomasal lymph node. The gene expression level of a candidate gene lying within the QTL region (PAPP-A2) was measured. In addition, putative interactions between the chromosome segment under study and the top ten differentially expressed genes between resistant MBB and susceptible RMN sheep highlighted in a previous microarray experiment were investigated. We found an induction of Th-2 related cytokine genes expression in the abomasal mucosa of R sheep. Down-regulation of the PAPP-A2 gene expression was observed between naïve and challenged sheep although no differential expression was recorded between challenged R and N sheep. The genotyping of this limited region should contribute to the ability to predict the intrinsic resistance level of sheep.


The failure of anthelmintic drugs is an issue of major concern throughout the world, especially for the control of small ruminants nematodes such as Haemonchus contortus[1]. Breeding animals with a better ability to resist infection by gastro-intestinal nematodes (GIN) has been proposed as an alternative strategy to drug treatment, and has already been implemented in Australia and New-Zealand [2, 3].

This selection relies on the existing between-animal variation in the acquired immune response against GIN [4] which is mostly related, in murine models or in the sheep, to mounting an efficient Th-2 biased immune response driven by the IL4, IL5 and IL13 cytokines [57]. The humoral profile associated with this Th2 response involves IgA, IgG and IgE antibodies that control larval colonization, worm development and fecundity [5, 8, 9]. However, an innate component is also involved in the anti-nematode response. For instance, recent findings suggest that lectins contribute to entrapping worms in a mucus sheath, and would subsequently facilitate their elimination [10, 11]. A previous microarray experiment comparing gene expression levels in resistant Martinik (MBB) and susceptible Romane (RMN) sheep infected by H. contortus found a stronger induction of Th2-related cytokines and also of lectin genes in MBB [12].

Numerous genetic mapping studies have identified regions of the genome explaining a non-negligible part of the inter-individual variation (known as Quantitative Trait Loci, QTL) in resistance to nematode infection [1319]. The use of the ovine-specific DNA SNP chip showed that resistance to nematodes was determined by many genes with weak effect and some limited regions explaining a higher proportion of the genetic variation [17, 19]. Candidate gene approaches have been carried out for the interferon gamma [2022] and the major histocompatibility complex loci [2326], although none of the other regions identified by genetic mapping strategies have been mined further. Identifying the mutations controlling ovine resistance to H. contortus should improve the ability to perform genetic selection by directly targeting the genes of interest through marker-assisted selection.

In a previous QTL mapping study for resistance to Haemonchus contortus, Sallé et al. found five QTL of greater interest on OAR5, 7, 12, 13 and 21 that affected Faecal Egg Count (FEC) and other parameters measured in a 1000 Martinik Black-Belly × Romane (MBB × RMN) back-cross (BC) lamb population [19]. A QTL region on OAR12, between 47 and 56 Mbp, explained 4% of observed variation in FEC at first and second infection and was detected in two different subsets of back-cross sheep and three independent sheep populations (Sarda*Lacaune [27], Merino [13] and in a free living Soay sheep population [14]). The purpose of this study was to perform a functional validation of this QTL region. Based on the within-family QTL detection results, some BC sheep were selected to produce BCxBC progeny that would carry two alleles associated with resistance or two alleles associated with susceptibility. BCxBC progenies were subsequently selected based on their genotypes for the investigated QTL and submitted to an exhaustive parasitological and haematological data collection. A gene expression study on the top ten differentially expressed genes between resistant MBB and susceptible RMN sheep highlighted in a previous microarray experiment [12] and on a candidate gene lying within the QTL region was also performed.

Materials and methods

Association between FEC at first infection and 4-SNP haplotype

Previous association analysis performed in a MBB × RMN BC flock (reported in [19]) identified a significant association between a 4-SNP haplotype (namely s39968, OAR12_62301297, OAR12_62347621, OAR12_62371899) located at 56.06 Mbp on OAR12 and FEC at first infection. Briefly, this analysis consisted in testing the effect of each 4-SNP haplotype window on the trait of interest, at every 0.05 Mbp. Haplotypes whose frequency was below 1% were discarded to limit standard error of the estimation. In addition, breed origin of the haplotype was taken into account, so that two identical-by-state (IBS) haplotypes were considered different if their breed origin was different.

From these results, two clusters of alleles with significant contrasted effects were defined (Figure 1, Additional files 1 and 2). A first group of two rare RMN alleles (GGAGRMN and AAAGRMN), subsequently denoted S, showed the most unfavourable effects, i.e. associated to the highest FEC during H. contortus infection (Figure 1). In contrast to these S alleles, a cluster of three alleles (subsequently denoted R) was significantly more favourable toward limiting H. contortus infection. A difference of 0.58 σp was estimated between the S and the R alleles for FEC at first infection. The most favourable allele was segregating in the MBB breed (AGCAMBB) but one RMN allele (GGCARMN) also belonged to this cluster supporting the resistance potential of this breed. Remaining alleles were considered as being neutral with respect to resistance for H. contortus infection (denoted N).

Figure 1
figure 1

Allelic effect of the 4-SNP haplotype estimated with the association analysis performed in the BC population. The thirteen alleles of the 4SNP haplotype associated with Faecal Egg Count at first infection that were segregating in the back-cross population are plotted against their estimated effect (given in phenotypic standard deviation). Alleles inherited from the resistant Martinik breed are plotted in black and corresponding sequences are annotated with “MBB”, while alleles from the Romane breed are shown in red and annotated with ‘RMN”. R, N and S alleles stand for alleles being favourable, neutral and unfavourable toward H. contortus infection respectively.

Production of the R and N sheep

To investigate the biological properties of the identified QTL region, a marker-assisted mating of BC sheep was performed to produce lambs carrying particular combination of QTL alleles, i.e. RR, RN or NN. BC sheep were selected according to the QTL allele they carried. Chromosomes of every BC sheep were reconstructed using their 50 K SNP genotypes (described in [19]) and the LinkPHASE software [28], so that the QTL region could be traced from pure breed grand-parents to BC lambs. Two BC sires with RN genotype and one BC sire with NN genotype were selected for mating with 73 BC ewes (45 NN, 26 RN and two RR ewes). Because of their low frequency, the S alleles were not segregating in the remaining BC population. To randomize as much as possible the distribution of the other QTL in the BCxBC progenies, the three sires were mated to RN and NN ewes. The two RR ewes were mated to the RN sires to increase the number of RR genotypes. In the end 130 BCxBC sheep were born at the La Sapinière experimental farm (Osmoy, France).

Sorting R and N sheep according to the association analysis

BCxBC sheep were genotyped with the 50 K ovine SNP chip (Illumina Inc, San Diego, CA, USA) and the same workflow as applied for their parents SNP data (described in [19]) was followed to select genotypes of interest. After data processing, 85 NN, 32 RN and 13 RR BCxBC progenies were counted. We retained the 54 sheep (22 NN, 20 RN and 12 RR) for which the chance of having inherited the QTL fragment from their ancestors was the highest. We finally compared the NN sheep (N group) against carriers of the R allele (R group, composed of the RN and RR sheep).

Infection procedure, pathophysiological measurements and tissue sampling

Lambs were kept indoors from birth to the end of the experiment, thus remaining totally worm-free before their infection. Lambs were transferred from the experimental farm (La Sapinière farm, Osmoy France) where they were born to the experimental facilities (Langlade farm, Pompertuzat, France). Upon arrival in the experimental facilities, the selected BCxBC sheep were given a Vecoxan ND treatment (diclazuril, 1 mg/kg bodyweight, Janssen) at the recommended dose to prevent any coccidiosis outbreak. They were subsequently left indoors for a one-month acclimation period. After checking that no strongyle eggs were excreted, 44 sheep were infected orally with 10 000 infective L3 larvae of the H. contortus strain used in the previous QTL mapping study [19]. Ten additional uninfected lambs, i.e. five of each susceptibility group, were not challenged and reserved for the gene expression analysis to determine basal gene expression level within each susceptibility group. For practical purpose, control lambs were sacrificed two days after the challenge took place, whereas the 44 other infected animals were euthanized at 30 and 31 dpi. Euthanasia was performed by a veterinary surgeon with a lethal intra-venous injection of embutramide (T61, 6 mL/50 kg bodyweight, Intervet). In agreement with the current French regulations at the time of the experiment (2011), INRA procedures for the care of experiment animals were applied.

Intra-rectal collection of faeces was performed every three days from 18 days post infection (dpi) until 30 dpi for FEC counting following the McMaster method modified by Raynaud [29]. These traits were denoted FEC18, FEC21, FEC24, FEC27 and FEC30. Blood samples were collected just before infection, at 14 dpi and 27 dpi. Samples were processed by the Sysmex XT-2000iV haematology analyser (calibrated for sheep) hence providing a complete screening of haematological parameters. Reticulocytes (denoted RET0, RET14 and RET27 for samples taken before, at 14 and 27 dpi respectively) and white blood cell counts were obtained, i.e. lymphocytes (LYMPH0, LYMPH14 and LYMPH27), monocytes (MONO0, MONO14, and MONO27), neutrophils (NEUT0, NEUT14, NEUT27), eosinophils (EO0, EO14, EO27) and basophils (BASO0, BASO14 and BASO27). Haematocrit was also determined (denoted HCT0, HCT14 and HCT27).

Following euthanasia, abomasal (gastric) lymph nodes (ALN) and a patch of the abomasal fundic mucosa (AFM) were sampled and stored at -20°C in RNAlater (Ambion, USA). Abomasal contents and washings were collected and put into absolute alcohol. Worm burden (WB) was determined using 10% of the total abomasal content. The lengths of 35 intact adult female worms were determined and averaged (denoted FL) for each lamb. To determine the average number of eggs in utero, 20 female worms were digested individually into a bleaching mixture (40 mL of Milton agent diluted into 160 mL of distilled water) and 10% of the resulting mixture was sampled for eggs counting using an optical microscope (denoted FF). The average number of eggs in utero per female H. contortus was subsequently obtained by averaging over the twenty individual females. In addition, sheep were weighed before and at the end of the experimental challenge and the difference between the two measurements divided by 30 was considered as average daily gain (ADG).

Gene expression measure

A gene expression comparison was performed between carriers of the GAAGMBB allele (Ri group, n = 9) and carriers of the GAAGRMN allele (Ni group, n = 8). These two particular alleles were chosen to be the most distinct in their effects and the most frequently distributed among the R and N sheep. Two groups of five uninfected control sheep (denoted Ru and Nu) were used to provide a basal gene expression level for the infected animals from the corresponding genotype.

To identify what genes underly the QTL, the annotated genes located within a 2-Mbp region centred at 56.06 Mbp on OAR12 were retrieved using the ovine genome website [30]. These genes are provided in Table 1. Given their known biological functions [31], we selected the PAPP-A2 gene for gene expression analysis as it was the only gene with obvious link to the immune response (Table 1).

Table 1 List of the genes lying within the QTL region under study

Additionally, a previous microarray study [12] found that, some genes (Table 2) were differentially expressed between infected MBB and RMN sheep, either in AFM (LGALS15, ITLN2, TFF3) or in ALN (TNFRSF4, CCL26, CXCL14) or both (IL4, IL5, IL13, TNFα, IFNγ). We thus performed a gene expression analysis on these candidates to assess whether they would be correlated to the QTL region under study.

Table 2 List of the genes considered for gene expression analysis

Total mRNA from abomasal fundic mucosa (AFM) and draining lymph nodes (ALN) of the 27 R and N sheep was extracted following the commercial RNeasy Mini Kit (Qiagen). The quality of the recovered RNA was monitored by A260/A280 spectrophotometry. RNAs were subsequently reverse-transcribed to cDNA with the Reverse Transcriptase commercial kit (Invitrogen).

Primers were designed for expression analysis using the “primer-BLAST” NCBI website [32]. Secondary structures were identified with the Mfold website [33] and selected primer sequences were blasted against the third version of the ovine genome to ensure specificity of their target. The qPCR was performed with three replicates per sample. Gene expression levels of five reference genes, namely β-ACTIN, TYQ, SDH, S26Q and HPRT, were measured. Their respective gene-wise stability values were estimated as reported in [34] and most stable genes were kept for subsequent analyses.

Statistical analyses

Strong departures from normality were found for FEC data (using the UNIVARIATE procedure implemented in the SAS software, SAS 2001, Cary NC) that were corrected by a 4th root transformation.

Resistance to nematodes is known to be polygenic and other QTL did segregate in this population [19]. Therefore, a genomic value (gEBV) evaluating the effect of the genomic background on FEC at first infection was estimated for each BCxBC sheep. FEC was considered as a proxy for resistance to nematodes in order to maximize the data available for estimating the gEBV by using the whole set of BC and BCxBC data (n = 1200 records). The computation was performed with the Bayes C genomic selection method [35] implemented in the GS3 software [36]. FEC was modelled as the sum of a mean, the already identified environmental fixed effects [19], i.e. sex, age at sampling and litter size, the marker effect and a residual term. None of the markers located on OAR12, i.e. the chromosome harbouring the QTL region, were included in the analysis as it was expected to partially sweep the genetic variance explained by the QTL into the gEBV as reported elsewhere [37].

The PROC MIXED procedure implemented in the SAS software [38] was used to test for significant differences between R and N genotypic groups for each of the measured parasitological and haematological data. Usually encountered environmental effects were considered as fixed effects and the computed gEBV was fitted to the model as a covariate to account for the effect of the rest of the genome. For haematological parameters, basal value of the considered parameter (indexed by 0) was considered as covariate to account for potential inter-individual variation before the beginning of the experiment.

Normality of the cycle time (CT) values distribution was checked using the Shapiro-Wilk test implemented in the R software [39]. One outlier from the Ri group (sheep 12493, see Additional file 3), with CT values for house-keeping genes exceeding a ± three standard deviations range, was discarded from the AFM dataset. Differential expression was tested following the 2-ΔΔCT method [40]. For each individual the mean CT value from the three replicates was expressed as a relative abundance to the average expression level from the reference genes, denoted ΔCT.

To assess any modifications before or after infection related to the QTL genotype, pair-wise comparisons between Ru and Nu and Ri and Ni groups respectively were tested. To assess whether any differences between Ri and Ni groups was due to a lack of change between the infected sheep and their naïve counterparts, the Ri to Ru and Ni to Nu comparisons were tested.

Pair-wise comparisons were computed as the difference between ΔCT of the two groups as ΔΔCT = ΔCT1 – ΔCT2. For comparing Ri and Ni gene expression levels, ΔCT of the infected sheep was expressed as a relative abundance to the average expression of the corresponding susceptibility control group and the subsequent ΔΔCT value computed as:

Δ Δ C T = Δ C TRi mean Δ C TRu - Δ C TNi mean Δ C TNu

This correction was applied to correct for putative differences in gene expression levels between RU and NU sheep before infection took place.

Fold change in gene expression between considered groups was computed as 2-ΔΔCT[40]. Subsequently, a Wilcoxon test was applied to determine any significant difference between the compared groups. To account for multiple testing, a nominal p-value of 1% was considered for significance leading to less than one significant difference occurring by chance. An additional suggestive threshold was considered for nominal p-value below 5%. The complete data processing was performed using an in-house R script [39].


Phenotypic comparison of the R and N sheep

The average gEBV of the R and N groups were equivalent (p = 0.34) hence allowing a standardized comparison of their QTL genotype.

A selection of parasitological and haematological data for the two groups is provided in Table 3 (the complete list of recorded traits and associated statistics are provided in Additional file 4). The 4-SNP-based clustering of the BCxBC sheep allowed prediction of true “high-” and “low-FEC” sheep, as illustrated by the 601 eggs/g and 11546 eggs/g difference obtained between the R and N groups at 18 and 30 dpi (p = 0.02 and 0.01 respectively in both cases, Table 3). Further, the 4-SNP genotype was associated with strong differences in the length and fertility of female worms (Table 3). Female H. contortus collected from R sheep were 1.6 mm shorter on average (p = 0.0013) and showed 1.5 times fewer eggs in utero (p = 5.10-4) than those recovered in N sheep (525 and 358 eggs in utero/female for the N and R groups respectively, Table 3). These parasitological findings also correlated with the reduced blood loss at 14 and 27 dpi in R lambs (p = 0.03), and a higher production of reticulocytes in the N lambs (Table 3). However no significant differences were observed for WB between the two groups (p = 0.73). From a production perspective, no differences of growth rates could be found between the two groups (p = 0.23).

Table 3 Phenotypic comparison of R and N sheep

Testing for differential gene expression

The gene expression profiles of each Ri and Ni as well as their respective control, Ru and Nu, are shown in Table 4, Figures 2 and 3. CT values for each measured gene are provided in Additional file 3.

Table 4 Differential expression of the selected gene set in abomasal mucosa (AFM) and draining lymph node (ALN)
Figure 2
figure 2

Gene expression level for a subset of candidate genes in the abomasal mucosa of each allelic group 30 days post-infection by H. contortus. A boxplot of the gene expression levels is provided for each of the considered allelic group. The horizontal bar within the box is the median value while open circles indicate values falling above the upper (or below the lower) quartile plus (minus) 1.5 times the interquartile distance. Ri and Ni are the infected resistant and susceptible lambs while Ru and Nu are the non-infected control resistant and susceptible lambs. Asterisks indicate suggestive (p < 0.05, shown as “*”) or significant (p < 0.01, shown as “**”) differences between groups.

Figure 3
figure 3

Gene expression level for a subset of candidate genes in the abomasal lymph node of each allelic group 30 days post-infection by H. contortus. A boxplot of the gene expression levels is provided for each of the considered allelic group. The horizontal bar within the box is the median value while open circles indicate values falling above the upper (or below the lower) quartile plus (minus) 1.5 times the interquartile distance. Ri and Ni are the infected resistant and susceptible lambs while Ru and Nu are the non-infected control resistant and susceptible lambs. Asterisks indicate suggestive (p < 0.05, shown as “*”) or significant (p < 0.01, shown as “**”) differences between groups.

Among the annotated genes retrieved within a 2 Mbp region centred at 56.06 Mbp (Table 1), the pappalysin gene (PAPP-A2) was identified as a putative candidate due to its regulatory effect on IGF-1. This gene was found down-regulated in the Ni in comparison to their uninfected counterparts in both AFM and ALN (p = 0.02 and 0.05 respectively). However the similar down-regulation trend observed between Ri and Ru sheep was not significant (Table 4, Figures 2 and 3). Neither was the PAPP-A2 expression differential between Ri and Ni groups (Table 4, Figures 2 and 3).

Ru and Nu sheep gene expression levels were compared to test for differences of the basal expression pattern linked to the QTL genotype. However no significant nor suggestive differences in gene expression levels was recorded between the Ru and Nu sheep.

Similarly, Ri and Ni sheep exhibited similar expression levels of AFM-specific genes involved in the innate response, i.e. lectin genes (LGALS15 and ITLN2) and the TFF3 gene (Table 4). A significant induction of LGALS15 mediated by H. contortus infection was observed between infected groups and their respective controls (p = 0.04 between the Ri and Ru groups and p = 0.02 between the Ni and Nu sheep, Table 4 and Figure 2).

Both suggestive and significant variations were observed between Ri and Ni sheep for some components of the acquired response (Table 4). Indeed, Ri sheep demonstrated a 4-fold increase in the IL-4 and IL-13 gene expression level in AFM in comparison to the Ni sheep (p = 0.04 and p = 0.02 for IL4 and IL13 respectively, Table 4 and Figure 2). The difference of IL13 expression between Ri and Ni seemed to be related to the lack of IL13 induction in Ni sheep as no significant difference was observed between Ni and Nu sheep (p = 0.72). Also, a slight suggestive down-regulation of TNFα was observed between Ri and Ni sheep (Table 4, Figure 3). The apparent stronger Th2 cytokine environment in the Ri sheep was also reinforced by a significant reduction by a 1/3 factor of the IFNγ expression (p = 0.01) in ALN between Ri and Ru sheep which could not be found between Ni and Nu sheep (Table 4, Figure 3). In addition, a non-significant down-regulation of TNFRSF4 was also observed between the same groups (Table 4, Figure 3).


The reported study aimed at mining the functional properties of a QTL region associated with resistance to H. contortus. The first objective was to validate the identified QTL through the comparison of individuals selected on their particular QTL alleles. The second related goal was to investigate functional mechanisms associated with the QTL, thanks to a wider range of phenotypes, including gene expression analysis. Abundant literature has been produced on the role of particular loci like the MHC [2426, 41, 42] or IFNγ [20, 22]. Other research teams mined the functional differences between divergent lines of sheep selected for low or high FEC [4346], but to our knowledge, this is the first attempt to functionally investigate the properties of a positional candidate without prior evidence of a functional candidate affecting resistance to GIN in sheep.

Our results demonstrated that a 4-SNP haplotype of OAR12 could discriminate between resistant and susceptible lambs. The genotypic groups based on this 4-SNP genotype showed significant reduction in FEC (ranging between 600 and 11 000 eggs/g difference between R and N groups) haematocrit drop and worm fecundity (1.5 times less eggs in utero in H. contortus females in the R sheep). Given that the two genotypic groups had identical genomic background, our findings provide strong functional support for the QTL signal detected by the mapping approach.

As a wider RNAseq experiment is to be undertaken, only a limited set of targeted genes were measured for gene expression as a primary screening. These genes were either directly lying within the QTL region under study, or had been previously involved as major candidates in the resistance difference between MBB and RMN breeds.

The nearest annotated gene underlying the 4SNP region was the PAPP-A2 gene. It encodes a protease cleaving Insulin Growth Factor Binding Protein-4 (IGFBP4), therefore increasing the bioavailability of the Insulin Growth Factor (IGF) [47]. The IGF gene plays a role in the immune response [4749] and was recently highlighted as a key player during wound healing associated with a nematode-induced Th2 response in a murine model [50]. In addition, PAPP-A2 expression is highly dependent of pro-inflammatory cytokines such as IFNγ and TNFα [47, 48] and it is induced during wound healing [47]. Interestingly, PAPP-A2 gene expression was significantly reduced in infected R and N sheep, between Nu and Ni (p = 0.02 and p = 0.05 in AFM and ALN respectively) although the reduction was not significant between Ri and Ru sheep (p = 0.09 and p = 0.07 in AFM and ALN respectively). To date, this is the first report of a relationship between PAPP-A2 and H. contortus infection. However no differences were observed between resistant and susceptible sheep after infection, and expression levels of PAPP-A2 did not follow the trends of other cytokines known to have regulatory effects, such as IL4 (up-regulation), IFNγ (down-regulation) or TNFα (up-regulation). Therefore the putative role of this gene still remains to be confirmed. Further, the observed gene expression levels might reflect some process that occurred earlier during the infection process preventing identification of differences between the different susceptibility groups. In particular, it would be interesting to examine IGF-1 and associated PAPP-A2 expression levels early in the infection process and as early as four dpi as reported for another murine model [50].

A previous microarray experiment identified a few genes differentially expressed between pure breed MBB and RMN sheep. Expression of these genes was measured to identify any relationship with the QTL genotype. From the expression data, no particular relationship could be drawn between this QTL and components of the innate response, such as lectins or TFF3. Such a relationship, if any, might occur at a different time point from the anti-H. contortus response. The only similar expression variation to the array study was the significant induction of LGALSL15 between naïve and infected sheep [12].

On the contrary, the observed cytokine gene expression ratios suggest an effect on worm fecundity through the mounting of a stronger Th-2 type environment as illustrated by the increase in IL4 and IL13 expression [5, 46]. Further, the IFNγ expression, known to be associated with susceptibility to GIN in murine models [5] was also down-regulated in Ri sheep in comparison to the Ru sheep. The additional slight down-regulation of TNFRSF4 (p = 0.07) between the same two groups constitutes another factor contributing to the mounting of a Th-2 environment against GIN while repressing the Th-1 response. Indeed the inhibition of the TNFRSF4 cytokine is known to induce a more efficient expulsion of helminths in mice models of nematode infection [51, 52]. Same findings of a stronger Th-2 environment were already reported while comparing pure breed Martinik and Romane sheep either at 4 and 30 dpi [53] or at 8 dpi [12]. IL5 is an additional Th-2 associated cytokine that was also found to be over-expressed in pure breed Martinik infected by H. contortus[12, 53] and in resistant Blackface sheep facing T. circumcincta infection [46]. In our experiment however, IL5 expression level was undetected after 40 qPCR cycles in 88% of the replicates, either using the primers published in [53] or after designing new primers (data not shown).

Proposing a detailed mode of action of the investigated QTL region would require more detailed investigation. However, extreme QTL alleles showed great contrasts in their effect on FEC, sheep blood loss and female worm fecundity. Since every lamb was inoculated with a similar infecting dose and no differences in worm burden could be found, the strong differences in female fertility reported here cannot be related to reduction in fecundity associated with density dependent effects [9, 54]. Hence, observed differences in worm fecundity seem to be directly mediated by the QTL under study. However, the QTL allelic group explained 26% of the variation for this trait at most, suggesting other factors are also involved as reported elsewhere [9, 54]. Further, significant differences in haematocrit were consistently observed between R and N groups at 14 and 27 dpi. Provided reticulocyte production was significantly higher in susceptible lambs after infection only, observed differences reflected a difference in true blood loss and not a higher regeneration ability of the resistant lambs. In addition, haematocrit at 27 dpi was negatively correlated to worm size (-0.43) and worm fecundity (-0.53). Both findings lead to the hypothesis that the investigated QTL region could limit worm feeding hence reducing their growth (shorter females) and fecundity (lower eggs recovered in utero). This could be mediated by the stronger Th-2 response that seems to be mounted by R sheep after infection. Indeed, local IgA response limits worm fecundity in T. circumcincta infection [9, 55] and in H. contortus infection [7]. More recent findings also support the relationship between T-cells number and worm female length [56]. Additional investigations on the characterization of subpopulations of T-cells in allelic carriers of each type could bring additional insights. As well, histological examination of abomasal mucosa of extreme animals could also confirm the stronger Th-2 response by measuring the eosinophilic infiltration and the number of mast cells. Additional detrimental factors like an increase of lectins in mucus could not be demonstrated in this study.

An additional meta-analysis of other parasite infection datasets gathered in the European 3SR project [57] is currently in process to confirm this region as a key player across multiple European breed, and eventually refine its location [58]. However, the back-cross design will not permit further refinement of the causative mutation explaining these variations. Therefore, pure breed association analyses are currently in progress to help refine the QTL position and a RNAseq experiment will be undertaken to widen the scope of functional candidates.


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

    Article  PubMed  Google Scholar 

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

    Article  Google Scholar 

  3. Karlsson LJE, Greeff JC: Genetic aspects of sheep parasitic diseases. Vet Parasitol. 2012, 189: 104-112.

    Article  CAS  PubMed  Google Scholar 

  4. Stear MJ, Fitton L, Innocent GT, Murphy L, Rennie K, Matthews L: The dynamic influence of genetic variation on the susceptibility of sheep to gastrointestinal nematode infection. J R Soc Interface. 2007, 4: 767-776.

    Article  PubMed Central  PubMed  Google Scholar 

  5. Allen JE, Maizels RM: Diversity and dialogue in immunity to helminths. Nat Rev Immunol. 2011, 11: 375-388.

    Article  CAS  PubMed  Google Scholar 

  6. Anthony RM, Rutitzky LI, Urban JF, Stadecker MJ, Gause WC: Protective immune mechanisms in helminth infection. Nat Rev Immunol. 2007, 7: 975-987.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Lacroux C, Nguyen TH, Andreoletti O, Prevot F, Grisez C, Bergeaud JP, Gruner L, Brunel JC, Francois D, Dorchies P, Jacquiet P: Haemonchus contortus (Nematoda: Trichostrongylidae) infection in lambs elicits an unequivocal Th2 immune response. Vet Res. 2006, 37: 607-622.

    Article  CAS  PubMed  Google Scholar 

  8. Kooyman FN, Schallig HD, Van Leeuwen MA, MacKellar A, Huntley JF, Cornelissen AW, Vervelde L: Protection in lambs vaccinated with Haemonchus contortus antigens is age related, and correlates with IgE rather than IgG1 antibody. Parasite Immunol. 2000, 22: 13-20.

    Article  PubMed  Google Scholar 

  9. Stear MJ, Bishop SC, Doligalska M, Duncan JL, Holmes PH, Irvine J, McCririe L, McKellar QA, Sinski E, Murray M: 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-652.

    Article  CAS  PubMed  Google Scholar 

  10. French AT, Knight PA, Smith WD, Brown JK, Craig NM, Pate JA, Miller HRP, Pemberton AD: Up-regulation of intelectin in sheep after infection with Teladorsagia circumcincta. Int J Parasitol. 2008, 38: 467-475.

    Article  CAS  PubMed  Google Scholar 

  11. French AT, Bethune JA, Knight PA, McNeilly TN, Wattegedera S, Rhind S, Miller HRP, Pemberton AD: The expression of intelectin in sheep goblet cells and upregulation by interleukin-4. Vet Immunol Immunopathol. 2007, 120: 41-46.

    Article  CAS  PubMed  Google Scholar 

  12. Liénard E, Foucras G, Prévot F, Grisez C, Bergeaud JP, François D, Bouvier F, Jacquiet P: Comparison of gene expression profiles between resistant and susceptible ovine breeds to Haemonchus contortus. 23rd World Association for the Advancement of Veterinary Parasitology, Buenos Aires, Argentina. 2011

    Google Scholar 

  13. Beh KJ, Hulme DJ, Callaghan MJ, Leish Z, Lenane I, Windon RG, Maddox JF: A genome scan for quantitative trait loci affecting resistance to Trichostrongylus colubriformis in sheep. Anim Genet. 2002, 33: 97-106.

    Article  CAS  PubMed  Google Scholar 

  14. Beraldi D, McRae AF, Gratten J, Pilkington JG, Slate J, Visscher PM, Pemberton JM: Quantitative trait loci (QTL) mapping of resistance to strongyles and coccidia in the free-living Soay sheep (Ovis aries). Int J Parasitol. 2007, 37: 121-129.

    Article  CAS  PubMed  Google Scholar 

  15. Davies G, Stear MJ, Benothman M, Abuagob O, Kerr A, Mitchell S, Bishop SC: Quantitative trait loci associated with parasitic infection in Scottish blackface sheep. Heredity. 2006, 96: 252-258.

    Article  CAS  PubMed  Google Scholar 

  16. Dominik S, Hunt PW, McNally J, Murrell A, Hall A, Purvis IW: Detection of quantitative trait loci for internal parasite resistance in sheep. I. Linkage analysis in a Romney x Merino sheep backcross population. Parasitol. 2010, 137: 1275-1282.

    Article  CAS  Google Scholar 

  17. Kemper KE, Emery DL, Bishop SC, Oddy H, Hayes BJ, Dominik S, Henshall JM, Goddard ME: The distribution of SNP marker effects for faecal worm egg count in sheep, and the feasibility of using these markers to predict genetic merit for resistance to worm infections. Genet Res. 2011, 93: 203-219.

    Article  CAS  Google Scholar 

  18. Matika O, Pong-Wong R, Wooliams JA, Bishop SC: Confirmation of two quantitative trait loci regions for nematode resistance in commercial British terminal sire breeds. Animal. 2011, 5: 1149-1156.

    Article  CAS  PubMed  Google Scholar 

  19. Sallé G, Jacquiet P, Gruner L, Cortet J, Sauvé C, Prévot F, Grisez C, Bergeaud JP, Schibler L, Tircazes A, François D, Pery C, Bouvier F, Thouly JC, Brunel JC, Legarra A, Elsen JM, Bouix J, Rupp R, Moreno CR: A genome scan for QTL affecting resistance to Haemonchus contortus in sheep. J Anim Sci. 2012, 90: 4690-4705.

    Article  PubMed  Google Scholar 

  20. Dervishi E, Uriarte J, Valderrábano J, Calvo JH: Structural and functional characterisation of the ovine interferon gamma (IFNG) gene: its role in nematode resistance in Rasa Aragonesa ewes. Vet Immunol Immunopathol. 2011, 141: 100-108.

    Article  CAS  PubMed  Google Scholar 

  21. Sayers G, Good B, Hanrahan JP, Ryan M, Sweeney T: Intron 1 of the interferon gamma gene: Its role in nematode resistance in Suffolk and Texel sheep breeds. Res Vet Sci. 2005, 79: 191-196.

    Article  CAS  PubMed  Google Scholar 

  22. Coltman DW, Wilson K, Pilkington JG, Stear MJ, Pemberton JM: A microsatellite polymorphism in the gamma interferon gene is associated with resistance to gastrointestinal nematodes in a naturally-parasitized population of Soay sheep. Parasitology. 2001, 122: 571-582.

    Article  CAS  PubMed  Google Scholar 

  23. Buítkamp J, Filmether P, Stear MJ, Epplen JT: Class I and class II major histocompatibility complex alleles are associated with faecal egg counts following natural, predominantly Ostertagia circumcincta infection. Parasitol Res. 1996, 82: 693-696.

    Article  PubMed  Google Scholar 

  24. Paterson S, Wilson K, Pemberton JM: Major histocompatibility complex variation associated with juvenile survival and parasite resistance in a large unmanaged ungulate population. Proc Natl Acad Sci U S A. 1998, 95: 3714-3719.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Hassan M, Good B, Hanrahan JP, Campion D, Sayers G, Mulcahy G, Sweeney T: The dynamic influence of the DRB1*1101 allele on the resistance of sheep to experimental Teladorsagia circumcincta infection. Vet Res. 2011, 42: 46-

    Article  PubMed Central  PubMed  Google Scholar 

  26. Hassan M, Hanrahan JP, Good B, Mulcahy G, Sweeney T: A differential interplay between the expression of Th1/Th2/Treg related cytokine genes in Teladorsagia circumcincta infected DRB1*1101 carrier lambs. Vet Res. 2011, 42: 45-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Moreno CR, Gruner L, Scala A, Mura L, Schibler L, Amigues Y, Sechi T, Jacquiet P, François D, Sechi S, Roig A, Casu S, Barillet F, Brunel JC, Bouix J, Carta A, Rupp R: QTL for resistance to internal parasites in two designs based on natural and experimental conditions of infection. 8th World Congress on Genetics Applied to Livestock Productions, Belo Horizonte, Brazil. 2006, 15–05

    Google Scholar 

  28. Druet T, Fritz S, Boussaha M, Ben-Jemaa S, Guillaume F, Derbala D, Zelenika D, Lechner D, Charon C, Boichard D, Gut IG, Eggen A, Gautier M: Fine mapping of quantitative trait loci affecting female fertility in dairy cattle on BTA03 using a dense single-nucleotide polymorphism map. Genetics. 2008, 178: 2227-2235.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Raynaud JP: Study of the efficiency of a quantitative coproscopic technic for the routine diagnosis and control of parasitic infestations of cattle, sheep, horses and swine. Ann Parasitol Hum Comp. 1970, 45: 321-342. (in French)

    CAS  PubMed  Google Scholar 

  30. Ovine (Texel) version 3.1 Genome Assembly. []

  31. Home - Gene - NCBI. []

  32. Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL: Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012, 13: 134-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034-

    Article  PubMed Central  PubMed  Google Scholar 

  35. Habier D, Fernando RL, Kizilkaya K, Garrick DJ: Extension of the bayesian alphabet for genomic selection. BMC Bioinformatics. 2011, 12: 186-

    Article  PubMed Central  PubMed  Google Scholar 

  36. Duchemin SI, Colombani C, Legarra A, Baloche G, Larroque H, Astruc J-M, Barillet F, Robert-Granié C, Manfredi E: Genomic selection in the French Lacaune dairy sheep breed. J Dairy Sci. 2012, 95: 2723-2733.

    Article  CAS  PubMed  Google Scholar 

  37. Riggio V, Matika O, Pong-Wong R, Stear MJ, Bishop SC: Genome-wide association and regional heritability mapping to identify loci underlying variation in nematode resistance and body weight in Scottish Blackface lambs. Heredity. 2013, 110: 420-429.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. AS 9.1 Documentation. []

  39. The Comprehensive R Archive Network. []

  40. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408.

    Article  CAS  PubMed  Google Scholar 

  41. Sayers G, Good B, Hanrahan JP, Ryan M, Angles JM, Sweeney T: Major histocompatibility complex DRB1 gene: its role in nematode resistance in Suffolk and Texel sheep breeds. Parasitology. 2005, 131: 403-409.

    Article  CAS  PubMed  Google Scholar 

  42. Stear MJ, Innocent GT, Buitkamp J: The evolution and maintenance of polymorphism in the major histocompatibility complex. Vet Immunol Immunopathol. 2005, 108: 53-57.

    Article  CAS  PubMed  Google Scholar 

  43. Ingham A, Reverter A, Windon R, Hunt P, Menzies M: Gastrointestinal nematode challenge induces some conserved gene expression changes in the gut mucosa of genetically resistant sheep. Int J Parasitol. 2008, 38: 431-442.

    Article  CAS  PubMed  Google Scholar 

  44. Ingham A, Menzies M, Hunt P, Reverter A, Windon R, Andronicos N: Divergent ghrelin expression patterns in sheep genetically resistant or susceptible to gastrointestinal nematodes. Vet Parasitol. 2011, 181: 194-202.

    Article  CAS  PubMed  Google Scholar 

  45. Gruner L, Bouix J, Vu Tien Khang J, Mandonnet N, Eychenne F, Cortet J, Sauvé C, Limouzin C: A short-term divergent selection for resistance to Teladorsagia circumcincta in Romanov sheep using natural or artificial challenge. Genet Sel Evol. 2004, 36: 217-242.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Gossner AG, Venturina VM, Shaw DJ, Pemberton JM, Hopkins J: Relationship between susceptibility of Blackface sheep to Teladorsagia circumcincta infection and an inflammatory mucosal T cell response. Vet Res. 2012, 43: 26-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Conover CA: Key questions and answers about pregnancy-associated plasma protein-A. Trends Endocrinol Metab. 2012, 23: 242-249.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Resch ZT, Chen B-K, Bale LK, Oxvig C, Overgaard MT, Conover CA: Pregnancy-associated plasma protein a gene expression as a target of inflammatory cytokines. Endocrinology. 2004, 145: 1124-1129.

    Article  CAS  PubMed  Google Scholar 

  49. Borghetti P, Saleri R, Mocchegiani E, Corradi A, Martelli P: Infection, immunity and the neuroendocrine response. Vet Immunol Immunopathol. 2009, 130: 141-162.

    Article  CAS  PubMed  Google Scholar 

  50. Chen F, Liu Z, Wu W, Rozo C, Bowdridge S, Millman A, Van Rooijen N, Urban JF, Wynn TA, Gause WC: An essential role for TH2-type responses in limiting acute tissue damage during experimental helminth infection. Nat Med. 2012, 18: 260-266.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Ierna MX, Scales HE, Schwarz H, Bunce C, McIlgorm A, Garside P, Lawrence CE: OX40 interactions in gastrointestinal nematode infection. Immunology. 2006, 117: 108-116.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Ekkens MJ, Liu Z, Liu Q, Whitmire J, Xiao S, Foster A, Pesce J, VanNoy J, Sharpe AH, Urban JF, Gause WC: The role of OX40 ligand interactions in the development of the Th2 response to the gastrointestinal nematode parasite Heligmosomoides polygyrus. J Immunol. 2003, 170: 384-393.

    Article  CAS  PubMed  Google Scholar 

  53. Terefe G, Lacroux C, Andreoletti O, Grisez C, Prevot F, Bergeaud JP, Penicaud J, Rouillon V, Gruner L, Brunel JC, Francois D, Bouix J, Dorchies P, Jacquiet P: Immune response to Haemonchus contortus infection in susceptible (INRA 401) and resistant (Barbados Black Belly) breeds of lambs. Parasite Immunol. 2007, 29: 415-424.

    Article  CAS  PubMed  Google Scholar 

  54. Fleming MW: Size of inoculum dose regulates in part worm burdens, fecundity, and lengths in ovine Haemonchus contortus infections. J Parasitol. 1988, 74: 975-978.

    Article  CAS  PubMed  Google Scholar 

  55. Stear MJ, Bishop SC: The curvilinear relationship between worm length and fecundity of Teladorsagia circumcincta. Int J Parasitol. 1999, 29: 777-780.

    Article  CAS  PubMed  Google Scholar 

  56. Rowe A, McMaster K, Emery D, Sangster N: Haemonchus contortus infection in sheep: parasite fecundity correlates with worm size and host lymphocyte counts. Vet Parasitol. 2008, 153: 285-293.

    Article  PubMed  Google Scholar 

  57. SR Project - Sustainable Solutions for Small Ruminants.,

  58. Riggio V, Pong-Wong R, Sallé G, Usai MG, Casu S, Moreno CR, Matika O, Bishop SC: A joint-analysis to identify loci underlying variation in nematode resistance in three European sheep populations. J Anim Breed Genet. in press

Download references


Staff of the Langlade and La Sapinière experimental units are acknowledged for the care and management of sheep. The authors are greatly in-debt to S.C. Bishop and A. Legarra for their fruitful comments and advice while analysing the dataset. Authors would also like to thank L. Bodin, M. Lasalle (Cryopic), S. Coppin and V. Jacquiet for their support during the experimental infection. Critical reviewing and constructive comments of this manuscript by anonymous reviewers, A. Blanchard-Letort and R. Beech were particularly appreciated. This project was funded by the European 3SR project and G. Sallé was granted by the INRA Animal Health and Animal Genetics departments.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Guillaume Sallé.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

GS designed the experiment, sampled animals, performed parasitological and haematological measurements and molecular biology, analysed data and drafted the manuscript. CM and PJ designed the experiment, sampled animals and were the grant holders. SB performed the hard selection signatures analyses. JR supervised the experiment and participated in animal samplings. MA and JL took care of the animals and FB carried out the sheep mating. AT, FP and CG sampled animals and performed molecular biology. CT performed haematological analyses. JP sampled animals and performed parasitological measures. EL sampled animals and analysed the data. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Frequency and estimated effects of the 4-SNP haplotype associated to Faecal Egg Count at first infection in the back-cross population. The estimated effects of the 4-SNP haplotypes identified in the back-cross population are provided in this file and given in phenotypic standard deviation. Frequencies of each allele are reported for both the back-cross population and the BCxBC progeny. a: AGCAMBB, AGCA allele from the Martinik Black-Belly breed; GGCARMN, GGCA allele from the Romane breed; b: allelic effect is given in phenotypic standard deviation; standard errors of the estimates are indicated in brackets. (DOC 38 KB)


Additional file 2: Pair-wise comparison of the effects for every allele identified in the back-cross population. For each couple of allele, a t-test has been applied considering respective estimated allelic effects and associated standard error from the QTL detection analysis and inferring the number of observations from the allelic frequency. The associated p-values are reported in the Table. A p-value below 0.05 (in bold) was considered as significant. The two clusters of alleles with significant contrasted effects are in italic. The MBB subscript indicates alleles inherited from the resistant Martinik breed (all other alleles segregated in the susceptible Romane breed). (DOC 42 KB)


Additional file 3: Geometric mean and associated standard deviation of measured Ct values in abomasal fundic mucosa (A) and abomasal lymph node (B) for each sheep × gene combination. The geometric means and associated standard deviations of the measured CT values for the three replicates are given for every sheep × gene combination. The five housekeeping genes were also provided. Any Ct value above 40 cycles were not considered for analyses; NA values for s indicates that either no or one Ct value was retained for analysis hence making it impossible to compute the mean or the standard deviation or both. Ri and Ru are respectively experimentally infected and uninfected resistant sheep while Ni and Nu are the genetically susceptible counterparts. A. Geometric mean (μ) and associated standard deviation (s) of measured Ct values for each sheep × gene combination in abomasal fundic mucosa (AFM). Any Ct value above 40 was discarded; NA values for s indicates that either no or one Ct value was retained for analysis. Ri and Ru are respectively experimentally infected and uninfected resistant sheep while Ni and Nu are the genetically susceptible counterparts. B. Geometric mean (μ) and associated standard deviation (s) of measured Ct values for each sheep × gene combination in abomasal lymph node (ALN). Any Ct value above 40 was discarded; NA values for s indicates that either no or one Ct value was retained for analysis. Ri and Ru are respectively experimentally infected and uninfected resistant sheep while Ni and Nu are the genetically susceptible counterparts. (DOC 248 KB)


Additional file 4: List of parasitological and haematological traits in the R and N sheep. The average performance of the R and N sheep is provided for every single trait monitored during experimental challenge by H. contortus and after necropsy. Average raw data are provided for the ease of reading, while statistical tests were performed after correction for fixed effect (e.g. sheep sex) and statistical transformation when appropriate. (DOC 66 KB)

Authors’ original submitted files for images

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

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sallé, G., Moreno, C., Boitard, S. et al. Functional investigation of a QTL affecting resistance to Haemonchus contortus in sheep. Vet Res 45, 68 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: