Skip to content


  • Research
  • Open Access

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

  • 1, 2Email author,
  • 3,
  • 4,
  • 3,
  • 3,
  • 5,
  • 6,
  • 6,
  • 7,
  • 7,
  • 8,
  • 7,
  • 7, 9 and
  • 7, 9
Veterinary Research201445:68

  • Received: 11 December 2013
  • Accepted: 4 June 2014
  • Published:


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.


  • Worm Burden
  • Female Worm
  • Pure Breed
  • Haemonchus Contortus
  • Abomasal Mucosa


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

Gene name


Related GO functions


55 199 014

DNA damage response


55 509 260

Cell cycle


55 757 260

Regulation of Insulin-like Growth Factor (IGF) activity


56 095 325

Nervous system development


56 545 767

Cell cycle, nervous system development


56 864 943

Initiation of transcription by RNA polymerase II

Annotated genes comprised within the 55.1 and 57.1 Mbp interval of OAR12 are provided. Their related gene ontology function has been retrieved from the NCBI’s website [31]. Position is given in bp and corresponds to the starting position as indicated on the ovine sheep genome website [30]. PAPP-A2*: as this gene was the only one within the QTL region with an obvious role in immunity, only this gene was examined in gene expression analysis.

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


Gene name



Related GO functions

Tissue in which the gene was tested

Within QTL




55 757 260

Regulation of Insulin-like Growth Factor (IGF) activity



Innate response




Carbohydrate binding





110 480 137

Response to nematode





260 633 891

Molecular function





19 237 699

Cytokine activity





19 436 357

Cytokine activity





19 261 902

Cytokine activity


Outside QTL

Acquired response



26 854 381

Cytokine activity





151 528 005

Cytokine activity





49 359 903

Tumor necrosis factor-activated receptor activity





33 952 859

Chemokine (C-C motif) ligand 26





44 432 648

Chemokine (C-X-C motif) ligand 14


The genes whose expression level has been measured are listed with their respective gene name, chromosome (OAR) and starting position (given in bp) and gene ontology function as provided by the NCBI’s website [31]. Expression levels were measured in the abomasal fundic mucosa (AFM)or in the abomasal lymph node (ALN) or both.

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





Phenotypic difference (σp)





FEC18 (eggs/g)







FEC21 (eggs/g)







FEC24 (eggs/g)







FEC27 (eggs/g)







FEC30 (eggs/g)







WB (no. worms)







FL (in mm)







FF (no. eggs in utero)







HCT14 (%)







HCT27 (%)







RET27 (%)







ADG (kg/day)







Raw parasitological data and haematological profiles of the R and N groups are reported in this Table. P-values were computed after correction for fixed effect and statistical transformation when appropriate. Phenotypic differences are the difference between the R and the N groups expressed in phenotypic standard deviation (σp).

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)



FC Ri vs Ni

FC Ri vs Ru

FC Ni vs Nu

FC Ru vs Nu






































Not tested in this tissue


Not tested in this tissue


Not tested in this tissue




























Not tested in this tissue


Not tested in this tissue


Not tested in this tissue





















Fold changes (FC) between the four considered groups (resistant allele carriers infected and uninfected (Ri and Ru respectively) and susceptible allele carriers infected and uninfected (Ni and Nu respectively) are provided. Significant differences (nominal p-value < 0.01) are indicated with **, and suggestive differences (nominal p-value < 0.05) are marked as *.

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.



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.

Authors’ Affiliations

INRA, UMR1282, Infectiologie et Santé Publique, F-37380 Nouzilly, France
Université François Rabelais de Tours, UMR1282 Infectiologie et Santé Publique, F-37000 Tours, France
INRA, UR631, Station d'Amélioration Génétique des Animaux, BP 27, F-31326 Castanet-Tolosan, France
INRA, UMR444, Laboratoire de Génétique Cellulaire, BP 27, F-31326 Castanet-Tolosan, France
INRA, UE332, Domaine de la Sapinière, F-18390 Osmoy, France
INRA, UE65, Domaine de Langlade, F-31450 Pompertuzat, France
INRA, UMR1225, Interactions Hôtes – Agents Pathogènes (IHAP), F-31076, BP 87614, Toulouse, France
Department of Clinical Sciences, INP-ENVT, Toulouse, France
Université de Toulouse, Institut National Polytechnique (INP), Ecole Nationale Vétérinaire de Toulouse (ENVT), UMR1225, IHAP, Toulouse, France


  1. Kaplan RM, Vidyashankar AN: An inconvenient truth: global worming and anthelmintic resistance. Vet Parasitol. 2012, 186: 70-78.View ArticlePubMedGoogle Scholar
  2. Bishop SC, Morris CA: Genetics of disease resistance in sheep and goats. Small Ruminant Res. 2007, 70: 48-59.View ArticleGoogle Scholar
  3. Karlsson LJE, Greeff JC: Genetic aspects of sheep parasitic diseases. Vet Parasitol. 2012, 189: 104-112.View ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle Scholar
  5. Allen JE, Maizels RM: Diversity and dialogue in immunity to helminths. Nat Rev Immunol. 2011, 11: 375-388.View ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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. 2011Google 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticleGoogle 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.View ArticleGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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-PubMed CentralView ArticlePubMedGoogle 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-PubMed CentralView ArticlePubMedGoogle 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–05Google 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.PubMed CentralView ArticlePubMedGoogle 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)PubMedGoogle 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-PubMed CentralView ArticlePubMedGoogle Scholar
  33. Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415.PubMed CentralView ArticlePubMedGoogle 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-PubMed CentralView ArticlePubMedGoogle Scholar
  35. Habier D, Fernando RL, Kizilkaya K, Garrick DJ: Extension of the bayesian alphabet for genomic selection. BMC Bioinformatics. 2011, 12: 186-PubMed CentralView ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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-PubMed CentralView ArticlePubMedGoogle Scholar
  47. Conover CA: Key questions and answers about pregnancy-associated plasma protein-A. Trends Endocrinol Metab. 2012, 23: 242-249.PubMed CentralView ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle Scholar
  55. Stear MJ, Bishop SC: The curvilinear relationship between worm length and fecundity of Teladorsagia circumcincta. Int J Parasitol. 1999, 29: 777-780.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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 pressGoogle Scholar


© Sallé et al.; licensee BioMed Central Ltd. 2014

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.