Impact Factor 4.151

Frontiers journals are at the top of citation and impact metrics

This article is part of the Research Topic

Advances in farm animal genomic resources

Original Research ARTICLE

Front. Genet., 24 March 2015 |

Fecal egg counts for gastrointestinal nematodes are associated with a polymorphism in the MHC-DRB1 gene in the Iranian Ghezel sheep breed

Rahman Hajializadeh Valilou1, Seyed A. Rafat1*, David R. Notter2, Djalil Shojda1, Gholamali Moghaddam1 and Ahmad Nematollahi3
  • 1Department of Animal Science, Faculty of Agriculture, University of Tabriz, Tabriz, Iran
  • 2Department of Animal and Poultry Sciences, Virginia Tech, Blacksburg, VA, USA
  • 3Department of Pathobiology, College of Veterinary Medicine, University of Tabriz, Tabriz, Iran

Genetic variation among sheep breeds in resistance to gastrointestinal nematodes (GIN) has been demonstrated in several production environments. Relationships between the ovine major histocompatibility complex and resistance to GIN have been studied, but few studies have systematically examined this issue in less-developed and semi-arid regions. The aim of the current study was to explore associations between fecal worm egg counts (FEC) for several GIN and polymorphisms in the DRB1 gene. One hundred male lambs were selected at 4–6 months of age from weaned animals in five flocks (n = 20 per flock). Body weights were determined, FAMACHA scores based on color of the ocular mucous membranes were assigned as an indicator of anemia, and blood and fecal samples were collected twice to evaluate FEC and blood packed cell volume (PCV) and for DNA isolation. A repeated-measures analysis of variance was used to test effects of genotype on FEC. The model included fixed effects of flock, genotype, time of measurement (1 or 2), and flock × time and genoype × time interactions, and a random (repeated) effect of lamb. Two genotypes (A1A1 and A1A2) were observed following digestion of Region 1 of Ovar-DRB1 with PstI. Genotypic frequencies were 0.73 for A1A1 and 0.27 for A1A2. FEC differed between Ovar_DRB1 genotypes A1A1 and A1A2 for Marshallagia marshalli, Strongyle, and total nematode FEC. Observed FEC were 30–41% lower for genotype A1A1. Differences among genotypes were consistent across measurement times, with no effect of genotype × measurement time interaction for any parasite class (P ≥ 0.34). A significant association was observed between FAMACHA scores and lamb PCV, and the residual correlation between these two variables was -0.51 (P < 0.001). FAMACHA scores can thus be used to detect differences among lambs in PCV, and polymorphic markers of Ovar-DRB1 have potential value as an indicator of parasite resistance in applied animal breeding programs on sheep farms in this region.


The Ghezel sheep is one of the 27 mainly fat-tailed native breeds of Northwestern Iran (Tavakolian, 2000). Animals of this breed graze for much of the year and are therefore continuously exposed to natural nematode infections. Gastrointestinal nematodes (GIN) of sheep and goats are widespread, diverse, and highly pathogenic, and can also infect other ruminant species such as cattle and reindeer (Jacquiet et al., 1998; Achi et al., 2003; Hrabok et al., 2006). Effects of GIN are most extreme in young animals, and therefore represent a real threat to the sheep industry (Waller, 2006).

The development of multi-drug resistance by GIN has driven research into alternative control measures, including selection of sheep that are genetically resistant to GIN infection. Genetic variation among breeds in resistance to GIN has been demonstrated in a variety of production environments. Genetically resistant sheep (either representing resistant local breeds or developed by selection within commercial breeds) are increasingly being used to improve animal production and well-being (Amarante et al., 2009). Genetically resistance sheep types also provide an opportunity to study novel mechanisms of resistance that may not be present in susceptible commercial breeds (Piedrafita et al., 2010). However, to date the mechanisms underlying genetic resistance of sheep to GIN infections are not precisely known.

Evidence for host genetic variation in aspects of disease resistance has now been documented for many diseases and in all major domestic livestock species (Bishop, 2005). In particular, small ruminants are notable for the large number of diseases where host genetic variation has been documented. Because parasite resistance in sheep has a moderate heritability (0.2–0.6; Baker, 1998; Stear et al., 2007), selective breeding has been used successfully with several breeds of sheep in different climates (Gray, 1987). Most sheep breeding programs for GIN resistance are based on recording of faecal egg counts (FEC), but this type of phenotype measurements is costly and difficult to collect on a large scale. In these situations, use of molecular genetic information is an interesting option. Use of molecular markers of resistance to GIN in sheep breeding programs has shown some promise, but difficulties remain, mainly because effects or previously identified quantitative trait loci (QTL) have not been consistent across breeds (Matika et al., 2011).

In sheep, class II genes of the major histocompatibility complex (MHC) are located on chromosome 20 and encode polymorphic glycoproteins composed of nine covalently linked subunits. Gruszczynska et al. (2000) found significant effect of OLA–DRB1 (MHC class II) on body weight at birth of Polish Heath sheep. There is also a body of scientific literature linking genes in the sheep MHC with the ability of sheep to resist infection by GIN as measured by FEC (Schwaiger et al., 1995; Buitkamp et al., 1996; Stear et al., 1996; Paterson et al., 1998). These findings are perhaps not surprising given the role of these genes in controlling specific immune responses. This “MHC effect” is never-the-less thought to be relatively small, accounting for an estimated 11% of total phenotypic variation in traits associated with GIN resistance (Buitkamp et al., 1996), although it accounts for a somewhat larger proportion of the additive genetic variation (Stear et al., 1997). These results have led to speculation that the MHC contains genes that could be used as markers for breeding to reduce FEC but that do not fully explain genetic resistance to GIN.

Few studies have addressed the relationships between the ovine MHC and resistance to GIN in less-developed and semi-arid regions. The current study is part of a multi-national collaborative project sponsored by the International Atomic Energy Agency and the Food and Agricultural Organization of the United Nations and designed to study genetic control of resistance to GIN in local sheep breeds. A particular focus of the study is the abomasal GIN Haemonchus contortus and was justified by this parasite’s ability to produce large numbers of eggs, resulting in extensive pasture contamination; the blood-sucking nature of this nematode, which can causes life-threatening levels of anemia; and the associated potential for very significant reductions in lamb performance and survival (Gatongi et al., 1998; Waller et al., 2004). In Iran, Ashrafi et al. (2014) demonstrated polymorphism in exon 2 of MHC gene OLA-DRB1 in the Makui sheep breed of Northwestern Iran. In this study, our aim was thus to explore possible associations between nematode resistance and polymorphism of DRB1. Because H. contortus is only one of several GIN known to be present in the temperate regions of Iran, we likewise focused our study on a variety of different GIN known to be important in the region (Garedaghi and Bahavarnia, 2013; Moradpour et al., 2013; Yagoob et al., 2013).

Material and Methods

Study Area

The study area is shown in Figure 1. The districts of Eastern and Western Azerbaijan provinces are agro-ecological zones and these zones are the site of origin and the habitat of the Ghezel sheep breed in Iran (Tavakolian, 2000). This breed was therefore selected for the present study. The study region is located at latitude 35–38.8°North and longitude 46–48°East and receives annual rainfall of 150–350 mm. The temperature is highest in June, before the onset of the monsoon season. During late spring and early to mid-summer, the daily maximum temperature rarely declines below 22°C. Relative humidity is lowest during April and May and rises during the monsoon season. The year is commonly divided into four seasons: winter (December–February), spring (March–April), summer (May–September), and autumn (October–November). The summer also includes the monsoon season (July–August; This study was conducted in May and June, a time when GIN numbers were anticipated to be elevated.


FIGURE 1. The study region was located in Northwestern Iran between 35 and 39°North latitude and between 46 and 48°East longitude.

Animals and Scheduling of Phenotype Sampling

One hundred male lambs at an age of 4–6 months were selected for this study. Each lamb was randomly selected from weaned animals within five flocks (n = 20 per flock). After deworming to eliminate existing nematode infection and when a parasite-free condition was confirmed (28 days later), the 20 lambs in each flock were allowed to graze together with untreated contemporaries from the same flock for at least 28 days without deworming. From day 31 post-infection, body weights were determined and blood and fecal samples were collected twice, 1 week apart, to evaluate fecal parasite egg counts and blood packed cell volume (PCV) and for DNA isolation. A single observer assigned FAMACHA scores for all flocks and both sampling times using a scale of 1–5. Scores were based on the color of the ocular mucous membranes surrounding the eye using procedures and color charts described by Vatta et al. (2001). All experimental procedures were approved by the University of Tabriz Animal Care and Ethics Committee.

Sample Processing

Individual fecal samples were collected from the rectum, processed to determine FEC using the modified McMaster technique, and reported as eggs per gram of feces. Observed parasite ova in the feces were categorized by parent species as: (1) Strongyles, (2) Nematodirus sp., (3) Trichuris sp., and (4) Marshallagia marshalli. The Strongyle group potentially included a number of common abomasal and intestinal sheep nematodes typical of mixed nematode infections in small ruminants such as H. contortus, Teladorsagia circumcincta, Ostertagia occidentalis, and Trichostrongylus axei, colubriformis, vitrinus, and rugatus. Fecal egg counts (FEC) were also summed across the four parasite classes for each lamb to derive a total nematode egg count. FEC were determined by the Clayton Lane technique.

Blood was obtained from the jugular vein with sterile vacuum tubes with anticoagulant (EDTA). For each sample, the PCV (%) was determined on the day of collection using the micro-hematocrit method. Blood was then mixed with 0.5 M of EDTA (pH = 8), and frozen at -20°C. DNA was isolated from blood using the protocol of Samadi Shams et al. (2011). Sequences of forward and reverse primers for amplification of the Ovar MHC-DRB1 (Region 2) gene are shown in Table 1 (Amills et al., 1995).


TABLE 1. Characteristics of PCR primers, nucleotide substitutions, restriction enzymes, amino acid changes, PCR producte size, and digested sequence sizes (allele) sizes for RFLP polymorphisms in Exon 2 of the Ovar-DRB1 gene.

The PCR was performed in a 25 μl reaction using the master mix kit (Ampliqon Company) in a T-Personal thermo-cycler (BiometeraPersonal Cycler Version 3.26 co., Germany). The PCR mixture contained: 50–100 ng of DNA, 2.5 μl of 10X PCR buffer (200 mM (NH4)2SO4, 0.1 mM Tween 20%, 750 mM Tris-HCl (pH 8.8), 2.5 mM MgCl2, 200 μM dNTPs, and 3 μl mix of oligonucleotids (10 pmol from each primer), 1U Taq DNA polymerase (Dream Taq polymerase, Ampliqon company) and 11 μl ddH2O. A total of 35 cycles was adapted for denaturation at 94°C/1 min, annealing at 61°C /1 min and polymerization at 72°C/2 min (Table 2). The PCR products were electrophoresed at 85 V for 45 min in 2.5% agarose gels, and visualized under UV light. The power supply for electrophoresis was a PAC1000 (Bio-Rad company; USA). The size of the alleles was determined based on a 100 bp DNA size standard (Ampliqon Company) using the computer software BIO 1D++. The PCR product for each sample was digested with 10 units of PstI and TaqI enzymes at 37 and 65°C, respectively. The characterization of each has been given in Table 3. The digested products were separated in a 2 and 3% agaros gel for 1 h at 85 V. The gels were stained with ethidium bromide.


TABLE 2. PCR temperatures program.


TABLE 3. Characteristics and features of TaqI and PstI enzymes and restriction sites.

Statistical Analysis

Fecal egg counts were analyzed separately for each class of parasite and for the total nematode egg count. For each parasite class, FEC values were expressed as residual deviations from flock × time subclass means, and distributions of residuals was tested for skewness (ω) and kurtosis (κ). If FEC were not normally distributed, FEC values were transformed before further analysis using Box–Cox transformations of the form [(FEC – 1)/; Box and Cox, 1964]. Optimum values of for a range of values between -2 and 2 were determined using a maximum-likelihood criterion (Draper and Smith, 1981) in the TransReg Procedure of SAS.

A repeated-measures analysis of variance was conducted using the MIXED Procedure of SAS to test effects of genotype on transformed FEC, body weights, and PCV. The model included fixed effects of flock, genotype, time of measurement (1 or 2), and flock × time and genoype × time interactions and a random (repeated) effect of lamb. This model was also fitted to untransformed FEC using the GLIMMIX Procedure of SAS and assuming a negative binomial distribution of FEC (O’Hara and Kotze, 2010). Associations between FAMACHA scores and measured variables were evaluated by adding effects of FAMACHA score and associated two-way interactions with other fixed effects to this mixed model.


Two genotypes (A1A1 and A1A2) were observed following digestion of Region 1 of Ovar-DRB1 with PstI (Figures 2 and 3). Two genotypes (B1B1 and B2B2) were likewise obtained following digestion with TaqI, but only one lamb had genotype B2B2 and effects of the TaqI polymorphism were therefore not considered further. For the PstI polymorphism, both genotypes were present in each of the five flocks, with overall genotypic frequencies of 0.73 for A1A1 and 0.27 for A1A2 (Table 4). The observed number of alleles at these two loci was similar to that pervious reported by Amills et al. (1995). The PCR product was sequenced (n = 4; Bioneer, Munpyeongseo-ro, Daedeok-gu, Daejeon 306–220, Republic of Korea). Figure 4 illustrates the sequence.


FIGURE 2. PCR products of the Ovar-DRB1 gene.


FIGURE 3. Genotypes and PCR-RFLP pattern of the digested Ovar-DRB1 gene digested (Left: digested with PstI; Right: digested with TaqI).


FIGURE 4. Nucleotide sequence of the DRB1 allele in one of the samples (animal no 10568). Primer complementary regions are indicated in bold type while the PstI sites are underlined.


TABLE 4. Genotype and gene frequencies.

Descriptive statistics for measured variables are shown in Table 5. The incidence of Trichurus infection in these data was very low. Eggs of Trichurus sp. were observed in fewer than 7% of the samples and FEC for Trichurus sp. were therefore excluded from the statistical analysis. Means and medians for FEC for remaining parasite classes are shown by flock and measurement time in Figure 5. Means for Strongyle FEC were very low in flocks 1 and 2, and means for M. marshalli FEC were very low in flocks 2, 3, and 4. The FEC records for these flocks for these parasite classes were therefore excluded from the final analyses.


TABLE 5. Descriptive statistics for fecal egg counts (FEC; eggs per gram of feces) for various classes of gastrointestinal nematodes and performance for 100 Ghezel lambs from five flocks.


FIGURE 5. Means and medians for Strongyle, Nematodirus sp., Marshallagia marshalli, and total nematode fecal egg counts (FEC; eggs/gram of feces) by flock and time of sampling. Medians are generally shown by the height of the blue bars and means are shown by the total height of the blue plus orange bars. In rare cases where the median exceeded the mean, the median is shown by the combined height of blue plus pink bars and the mean is shown by the height of the blue bars.

After removing FEC records for flocks 1 and 2 and expressing FEC as residual deviations from flock × time subclass means, the distribution of M. marshalli FEC was somewhat skewed to the right (ω = 0.49; P < 0.05) but did not exhibit significant kurtosis (κ = -0.04). This result was in agreement with the similarity between means and medians for M. marshalli FEC in Figure 5. A Box–Cox transformation with λ = 0.5 reduced the observed level of skewness (ω = -0.28), and a square-root transformation was used in the final analysis of M. marshalli FEC. As expected from differences between means and medians in Figure 5, Nematodirus, Strongyle, and total nematode FEC were not normally distributed. For these parasite classes, distributions of FEC were strongly skewed to the right (ω ≥ 1.28) and were leptokurtic (κ ≥ 3.51). Distributions of Nematodirus and Strongyle FEC also exhibited clumping at zero. Estimates of λ for Box–Cox transformations were close to zero for these parasite classes, with λ = 0.1 for total nematode and Strongyle FEC and λ = -0.1 for Nematodirus FEC. The Box–Cox transformation is undefined at λ = 0 but asymptotically approaches a logarithmic transformation as λ approaches zero. We therefore used a simple logarithmic transformation [ln(FEC + 1)] for these parasite classes.

Results of the repeated-measures analysis (Table 6) indicated that flock effects on FEC were large (P < 0.001) for all remaining parasite classes. Significant differnces were observed between Ovar_DRB1 genotypes A1A1 and A1A2 for M. marshalli, Strongyle, and total nematode FEC (Figure 6). Means for M. marshalli FEC in Figure 6 were based on untransformed FEC and were 40% lower for lambs of genotype A1A1 compared to lambs of genotype A1A2. Means for Strongyle, Nematodirus, and total nematode FEC were backtransformed from means of transformed variables (m) as (em-1), and SEs for backtransformed means were approximated by assuming that SEMs for log-transformed FEC were approximately equal to coefficients of variation of backtransformed means. Backtransformed means for these parasite classes indicated that Strongyle and total nematode FEC were 41 and 30% lower, respectively, for lambs of genotype A1A1 compared to lambs of genotype A1A2. No effect of genotype was observed for Nematodirus FEC Differences among genotypes were consistent among flocks and measurement times, with no effect of genotype × flock (P ≥ 0.23) or genotype × measurement time interaction (P ≥ 0.34 for any parasite class.


FIGURE 6. Genotype means for FEC (egg/gram of feces) for the PstI restriction-fragment polymorphism in the Ovar-DRB1 gene for several parasite classes. Values for Strongyle, Nematodirus sp., and total nematode FEC were transformed as ln(FEC + 1) for statistical analysis; least-squares means (m) were back-transformed to the original scale as (em- 1) for presentation and are geometric, rather than arithmetic, means. Values for M. marshalli FEC were transformed as the square root of FEC for analysis; least-squares means shown here are based on untransfomed FEC.


TABLE 6. Results of analysis of variance for fixed effects of flock, genotype, time of measurement, and their two-way interactions on FEC for each parasite class.

The impact of logarithmic transformation of Strongyle and total nematode FEC can be seen by comparing differences among genotypes for untransformed FEC. For untransformed FEC, means for lambs of genotype A1A1 were 29 (P = 0.20) and 23% (P = 0.12) lower compared to lambs of genotype A1A2. These differences in effect of genotype reflect the greater impact of occasional large FEC from the skewed right tail of the FEC distribution on mean differences between genotypes and on the residual variance. This issue could be addressed in untransformed data by removing records with very high FEC as outliers. We considered recoding of extreme values to be preferable to removing them, and use of a normalizing transformation provides an objective strategy to account for the presence of extreme values.

Results from analysis of untransformed FEC assuming a negative binomial distribution did not differ greatly from those from the analysis of transformed FEC. Means for untransformed M. marshalli, Strongyle, and total nematode FEC assuming a negative binomial distribution were 40 (P = 0.02), 18 (P = 0.01), and 43% (P = 0.15) lower, respectively, for lambs of genotype A1A1 compared to lambs of genotype A1A2. A significant effect of genotype was thus confirmed for M. marshalli and Strongyle FEC, but not for total nematode FEC, perhaps in association with different contributions of the various parasite classes to total nemadode FEC among flocks and measurement times.

Consistent effects of measurement time on FEC were not observed for any parasite class, but flock × measurement time interaction was significant for Strongyle and total nematode FEC (Table 6). The interaction was explained by a threefold increase in Strongyle FEC between the first and second measurement time in flock 3 (Figure 5). Significant difference in FEC between measurements times were not observed in other flocks.

Effects of genotype were not observed for lamb body weight or PCV. The distribution of FAMACHA scores revealed little evidence of anemia in these lambs, with frequencies across flocks and measurement times for FAMACHA scores of 1 through 5 of 40, 44,15, 1, and 0%, respectively. No association was observed beteween FAMACHA scores and lamb body weights or FEC for any parasite class. However, a significant association was observed between FAMACHA scores and lamb PCV (Figure 7). The PCV declined linearly as FAMACHA scores declined from 1 to 3 and were much lower for the two lambs that received a FAMACHA score of four. After adjusting for effects of herd and measurement time, a residual correlation of -0.51 (P < 0.001) was observed between PCV and the FAMACHA score.


FIGURE 7. Least-squares means and SEs for effects of FAMACHA score on packed cell volume (PCV).


Breeding for resistance to nematode infection can complement the use of anthelmintics in sheep husbandry. Resistant animals can be selected on the basis of low FEC (Eady et al., 2003; Kahn et al., 2003), and measurement of FEC is generally considered to be the standard method for assessment of the level of resistance to GIN. The number of eggs is easy to measure, indicates parasitism per se, and correlates well to the number of adult nematodes present in lambs (Douch et al., 1996; Baker, 1999). However, FEC is affected by several factors, such as parasite fecundity and egg-laying patterns, variations in egg distribution in feces, diet composition, intestinal transit time, and the level of immunity (Douch et al., 1996). A disadvantage of FEC as a marker of resistance is the requirement that animals be infected, either naturally or artificially, to determine the FEC value. The effort and cost in obtaining FEC measurements can also be a disadvantage, especially under extensive production conditions. Under some production conditions, it is therefore difficult to assess the resistance status of animals for breeding programs and there is, consequently, considerable interest in the evaluation of phenotypic and genetic markers associated with parasite resistance. Other phenotypic measures such as degree of anemia, circulating eosinophil counts, antibody levels to larval, or adult stages, and plasma pepsinogen concentrations can be used to predict worm burdens and resistance levels in infected sheep, but phenotypic markers that allow accurate prediction of an individual’s resistance status in the absence of infection are generally not available (Albers et al., 1987; Beh and Maddox, 1996).

A number of studies around the world have attempted to identify relationships between genetic resistance to GIN and various genes and genetic markers. Some of these studies have been summarized in Table 7. Polymorphisms within the ovine MHC complex were associated with resistance to T. colubriformis (Douch and Outteridge, 1989), T. circumcincta (Schwaiger et al., 1995), and H. contortus (Luffau et al., 1990; Outteridge et al., 1996). However, Cooper et al. (1989), Blattman et al. (1993) and Hulme et al. (1993) did not find evidence of an association between polymorphisms in the ovine MHC locus and resistance or susceptibility to H. contortus. A list of QTL for GIN resistance in sheep was provided by Dominik (2005). This review suggests that there is considerable evidence for an important role for the MHC in parasite susceptibility and resistance to H. contortus.


TABLE 7. Comparision of of results of the current study with other investigations of Ovar-DRB1 or nearby genes.

Sallé et al. (2012) reported four QTL regions on sheep chromosomes (OAR) 5, 12, 13, and 21 in Romane X Martinik Black Belly backcross lamb that had an important role in genetic resistance to H. contortus. Riggio et al. (2014) suggested other regions of OAR1, 3, 4, 5, 7, 19, 20, and 24 that were involved in GIN resistance. Davies et al. (2006) likewise found evidences for QTL on chromosomes 2, 3, 14, and 20 that were associated with parasitic infections in Scottish blackface sheep.

Results of the current study provided additional evidence of an association between polymorphism in the DRB1 gene and GIN FEC, and the first indication of an effect of this locus on M. marshalli FEC. Screening of GIN FEC under natural infection is informative because it corresponds to typical conditions in the production environment. However, possible interactions among effects of different parasite classes and variation among flocks and measurement times in parasite loads preclude a deeper understanding of mechanisms driving the observed associations. More intersive studies, involving controlled infections with individual species of GIN are thus required to confirm hypothesized effects of DRB1 genotype on GIN parasite resistance and to confirm the specificity or generality of observed associations among parasite classes.

A significant negative association between PCV and FAMACHA score confirmed that the FAMACHA score can be used to diagnose differences in PCV in lambs and was consistent with previous results (Kaplan et al., 2004). Recommendations for veterinary intervention based on FAMACHA scores (Vatta et al., 2001) suggest that lambs with scores of one or two do not require attention, but that lambs with scores of four or five require immediate attention. Recommendations for animals with a score of three depend upon the age and nutritional status of the lamb and anticipated cause(s) of anemia. Intervention is generally recommended for lambs, but not adults, with a FAMACHA score of three. In the current study, the lack of association between FAMACHA scores and FEC, low overall FEC levels, and limited evidence for parasitism by H. contortus (the only GIN, of those evaluated, that causes blood loss and anemia) suggest that parasitism probably is not the main cause of subclinical anemia in these flocks. Nonetheless, results in Figure 7 indicate that FAMACHA scores can be used to detect mean differences in PCV among lambs.


Our results reinforce previous evidence that some alleles of the ovine MHC are involved in determining levels of susceptibility or resistance to infection with GIN. This result provides the opportunity to use these alleles as genetic markers of resistance to GIN, leading to the development of that are better adapted to parasite infestations in the environment.

Implication of this research are that FAMACHA tests and polymorphic markers of Ovar-DRB1 can be used in applied animal breeding programs on sheep farms of the region, especially in animals infected with GIN and located in the temperate regions of Asia. Assessment of the precision of genetic evaluations based on molecular information has potential to provide a new perspective on the design of sheep breeding schemes and selection programs (Assenza et al., 2014).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


We acknowledge the assistance of the FAO/IAEA program, Contract No 16101. We acknowledge A. Javanmard, A. Barzgari, S. Radmand, and H. Cheraghi for help in this research. We appreciate constructive comments by two anonymous reviewers regarding options for statistical analysis of FEC.


Achi, Y. L., Zinsstag, J., Yao, K., Yeo, N., Dorchies, P., and Jacquiet, P. (2003). Host specificity of Haemonchus spp. for domestic ruminants in the savanna in northern Ivory Coast. Vet. Parasitol. 116, 151–158. doi: 10.1016/S0304-4017(03)00258-9

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Albers, A. J., Gray, G. D., Piper, L. R., Le Jambre, L. F., and Barger, I. A. (1987). The genetics of resistance and resilience to Haemonchus contortus infection in young merino sheep. Int. J. Parasitol. 17, 1355–1363. doi: 10.1016/0020-7519(87)90103-2

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Amarante, A. F. T., Susin, I., Rocha, R. A., Silva, M. B., Mendes, C. Q., and Pires, A. V. (2009). Resistance of santa ines and crossbred ewes to naturally acquired gastrointestinal nematode infections. Vet. Parasitol. 165, 273–280. doi: 10.1016/j.vetpar.2009.07.009

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Amills, M., Francino, O., and Sanchez, A. (1995). Nested PCR allows the characterization of TaqI and PstI RFLPs in the second exon of the caprine MHC class II DRB gene. Vet. Immunol. Immunopathol. 48, 313–321. doi: 10.1016/0165-2427(95)05442-9

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Ashrafi, F., Hashemi, A., Mardani, K., and Darvishzadeh, R. (2014). Study on genetic variability in MHC-DRB1 second exon in makuie sheep breed population. Genetika 46, 269–275. doi: 10.2298/GENSR1401269A

CrossRef Full Text | Google Scholar

Assenza, F., Elsen, J. M., Legarra, A., Carré, C., Sallé, G., Robert-Granié, C.,et al. (2014). Genetic parameters for growth and faecal worm egg count following Haemonchus contortus experimental infestations using pedigree and molecular information. Genet. Sel. Evol. 46, 13. doi: 10.1186/1297-9686-46-13

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Baker, R. L. (1998). A review of genetic resistance to gastrointestinal nematode parasites in sheep and goats in the tropics and evidence for resistance in some sheep and goat breeds in sub-humid coastal Kenya. Anim. Gen. Res. Inform. Bull. 24, 13–30. doi: 10.1017/S1014233900001103

CrossRef Full Text | Google Scholar

Baker, R. L. (1999). Genetics of resistance to endoparasites and ectoparasites. Int. J. Parasitol. 29, 73–75. doi: 10.1016/S0020-7519(98)00181-7

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Beh, K. J., and Maddox, J. F. (1996). Prospects for development of genetic markersf for resistance to gastrointestinal parasite infection in sheep. Int. J. Parasitol. 26, 879–897. doi: 10.1016/S0020-7519(96)80060-9

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Benavides, M. V., Weimer, T. A., Borba, M. F. S., and Sacco, A. M. S. (2009). Genetic analyses of polymorphisms on ovine chromosomes 5 and 20 and their effect on resistance to internal parasites. Small Rumin. Res. 83, 67–73. doi: 10.1016/j.smallrumres.2009.03.006

CrossRef Full Text | Google Scholar

Bishop, S. C. (2005). “Disease resistance: genetics,” in Encyclopedia of Animal Science, eds W. G. Pond and A. W. Bell (New York: Marcel Dekker), 288–290.

Google Scholar

Blattman, A. N., Hulme, D. J., Kinghorn, B. P., Woolaston, R. R., Gray, G. D., and Beh, K. J. (1993). A search for associations between major histocompatibility complex restriction fragment length polymorphism bands and resistance to Haemonchus confortus infection in sheep. Anim. Genet. 24, 277–282. doi: 10.1111/j.1365-2052.1993.tb00311.x

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Box, G. E. P., and Cox, D. R. (1964). An analysis of transoformations. J. R. Stat. Soc. Ser. B 26, 211–252.

Google Scholar

Buitkamp, J., Filmether, P., Stear, M. J., and Epplen, J. T. (1996). Class I and class II major histocompatibility complex alleles are associated with faecal egg counts following natural, predominantly Ostertagia circumcincta infection. Parasitol. Res. 82, 693–696. doi: 10.1007/s004360050187

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Charon, K. M., Moskwa, B., Rutkowski, R., Gruszczyñska, J., and Świderek W. (2002). Microsatellite polymorphism in DRB1 gene (MHC class II) and its relation to nematode faecal egg count in Polish Heath Sheep. J. Anim. Feed Sci. 11, 47–58.

Google Scholar

Coltman, D. W., Wilson, K., Pilkington, J. G., Stear, M. J., and Pemberton, J. M. (2001). A microsatellite polymorphism in the gamma interferon gene is associated with resistance to gastrointestinal nematodes in a naturally-parasitized population of Soay sheep. Parasitology 122, 571–582. doi: 10.1017/S0031182001007570

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Cooper, D. W., Van Oorschot, R. A. H., Piper, L. R., and Le Jambre, L. F. (1989). No association between the ovine leucocyte antigen (OLA) system in the australian merino and susceptibility to Haemonchuscontortus infection. Int. J. Parasitol. 19, 695–697. doi: 10.1016/0020-7519(89)90051-9

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Davies, G., Stear, M. J., Benothman, M., Abuagob, O., Kerr, A., Mitchell, S.,et al. (2006). Quantitative trait loci associated with parasitic infection in Scottish blackface sheep. Heredity (Edinb) 96, 252–258. doi: 10.1038/sj.hdy.6800788

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Dominik, S. (2005). Quantitative trait loci for internal nematode resistance in sheep: a review. Genet. Sel. Evol. 37(Suppl. 1), S83–S96. doi: 10.1051/gse:2004027

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Douch, P. G. C., Green, R. S., Morris, C. A., McEwan, J. C., and Windon, R. G. (1996). Phenotypic markers for selection of nematode-resistant sheep. Int. J. Parasitol. 26, 899–911. doi: 10.1016/S0020-7519(96)80062-2

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Douch, P. G. C., and Outteridge, P. M. (1989). The relationship between ovinelym-phocyte antigens and parasitological and production parameters in Romney sheep. Int. J. Parasitol. 19, 35–41. doi: 10.1016/0020-7519(89)90019-2

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Draper, N. R., and Smith, H. (1981). Applied Regression Analysis, 2nd Edn. New York: John Wiley & Sons.

Google Scholar

Eady, S. J., Woolaston, R. R., and Barger, I. A. (2003). Comparison of genetic and nongenetic strategies for control of gastrointestinal nematodes of sheep. Livest. Prod. Sci. 81, 11–23. doi: 10.1016/S0301-6226(02)00197-5

CrossRef Full Text | Google Scholar

Figueroa Castillo, J. A., Mendez Mendina, R. D., BerruecosVillalbos, J. M., Gayosso-Vazquez, A., Ulloa-Arvizu, R., Acosta Rodriguez, R.,et al. (2011). Association between major histocompatibility complex microsatellites, fecal egg count, blood packed cell volume and blood eosinophilia in Peibuey sheep infected with Haemonchus contortus. Vet. Parasitol. 177, 339–344. doi: 10.1016/j.vetpar.2010.11.056

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Garedaghi, Y., and Bahavarnia, S. R. (2013). Prevalence and species composition of abomasal nematodes in sheep and goats slaughtered at Tabriz town, Iran. J. Anim. Sci. Adv. 3, 37–41. doi: 10.5455/jasa.20130219031609

CrossRef Full Text | Google Scholar

Gatongi, P. M., Prichard, R. K., Ranjan, S., Gathuma, J. M., Munyua, W. K., Cheruyot, H.,et al. (1998). Hypobiosis of Haemonchus contortus in natural infections of sheep and goats in a semi-arid area of Kenya. Vet. Parasitol. 77, 49–61. doi: 10.1016/S0304-4017(97)00207-0

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Gray, G. D. (1987). Genetic resistance to haemonchosis in sheep. Parasitol. Today 3, 253–255. doi: 10.1016/0169-4758(87)90155-4

CrossRef Full Text | Google Scholar

Gruszczynska, J., Charon, K. M., Kitlinska, J., and Szydlowski, M. (2000). The influence of OLA-DRB1 (MHC class II) gene polymorphism on lamb body weight and weight gain in Polish Heath Sheep. J. Appl. Genet. 41, 101–112.

Google Scholar

Hrabok, J. T., Oksanen, A., Nieminen, M., and Waller, P. J. (2006). Prevalence of gastrointestinal nematodes in winter slaughtered reindeer of northern Finland. Rangifer 27, 133–139. doi: 10.7557/

CrossRef Full Text | Google Scholar

Hulme, D. J., Nicholas, F. W., Windon, R. G., Brown, S. C., and Beh, K. J. (1993). The MHC class II region and resistance to an intestinal parasite in sheep. J. Anim. Breed. Gen. 110, 459–472. doi: 10.1111/j.1439-0388.1993.tb00758.x

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Jacquiet, P., Cabaret, J., Thiam, E., and Cheikh, D. (1998). Host range and the maintenance of Haemoncus spp. in an adverse arid climate. Int. J. Parasitol. 28, 253–261. doi: 10.1016/S0020-7519(97)00185-9

CrossRef Full Text | Google Scholar

Kahn, L. P., Knox, M. R., Walkden- Brown, S. W., and Lea, J. M. (2003). Regulation of the resistance to nematode parasites of single- and twin-bearing Merino ewes through nutrition and genetic selection. Vet. Parasitol. 114, 15–31. doi: 10.1016/S0304-4017(03)00099-2

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Kaplan, R. M., Burke, J. M., Terrill, T. H., Miller, J. E., Getz, W. R., Mobini, S.,et al. (2004). Validation of the FAMACHA© eye color chart for detecting clinical anemia on sheep and goat farms in the southern United States. Vet. Parasitol. 123, 105–120. doi: 10.1016/j.vetpar.2004.06.005

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Luffau, G., Vu TienKhang, J., Bouix, J., Nguyen, T. G., Cullen, P., and Ricordeau, G. (1990). Resistance to experimental infections with Haemonchus contortus in Romanov sheep. Genet. Sel. Evol. 22, 205–229. doi: 10.1186/1297-9686-22-2-205

CrossRef Full Text | Google Scholar

Matika, O., Pong-Wong, R., Wooliams, J. A., and Bishop, S. C. (2011). Confirmation of two quantitative trait loci regions for nematode resistance in commercial British terminal sire breeds. Animal 5, 1149–1156. doi: 10.1017/S175173111100022X

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Moradpour, N., Borji, H., Razmi, G., Maleki, M., and Kazimi, H. (2013). Pathophysiology of Marshallagia marshalli in experimentally infected lambs. Parasitology 140, 1762–1767. doi: 10.1017/S0031182013001042

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

O’Hara, R. B., and Kotze, D. J. (2010). Do not log-transform count data. Meth. Ecol. Evol. 1, 118–122. doi: 10.1111/j.2041-210X.2010.00021.x

CrossRef Full Text | Google Scholar

Outteridge, P. M., Anderson, L., Douch, P. G. C., Green, R. S., Gwakisa, P. S., Hohenhaus, M. A.,et al. (1996). The PCR typing of MHC-DRB genes in the sheep using primers for an intronic microsatellite: application to nematode parasite resistance. Immunol. Cell Biol. 74, 330–336. doi: 10.1038/icb.1996.59

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Paterson, S., Wilson, K., and Pemberton, J. M. (1998). Major histocompatibility complex variation associted with juvenile survival and prasite resistance in a large ungulate population (Ovis aries L.). Proc. Natl. Acad. Sci. U.S.A. 95, 3714–3719. doi: 10.1073/pnas.95.7.3714

CrossRef Full Text

Piedrafita, D., Raadsma, H. W., Gonzalez, J., and Meeusen, E. (2010). Increased production through parasite control: can ancient breeds of sheep teach us new lessons? Trends Parasitol. 26, 568–573. doi: 10.1016/

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Riggio, V., Pong-Wong, R., Salle, G., Usai, M. G., Casu, S., Moreno, C. R.,et al. (2014). A joint analysis to identify loci underlying variation in nematode resistance in three European sheep populations. J. Anim. Breed. Genet. 131, 426–436. doi: 10.1111/jbg.12071

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Sallé, G., Jacquiet, P., Gruner, L., Cortet, J., Sauvé, C., Prévot, F.,et al. (2012). A genome scan for QTL affecting resistance to Haemonchus contortus in sheep. J. Anim. Sci. 90, 4690–4705. doi: 10.2527/jas.2012-5121

PubMed Abstract | Full Text | CrossRef Full Text

Samadi Shams, S., Zununi Vahed, S., Soltanzad, F., Kafil, V., Barzegari, A., Atashpaz, S.,et al. (2011). Highly effective DNA extraction method from fresh, frozen, dried and clotted blood samples. BioImpacts 1, 183–187. doi: 10.5681/bi.2011.025

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Sayers, G., Good, B., Hanrahan, J. P., Ryan, M., Angles, J. M., and Sweeney, T. (2005). Major histocompatibility complex DRB1 gene: its role in nematode resistance in Suffolk and Texel sheep breeds. Parasitology 131, 403–409. doi: 10.1017/S0031182005007778

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Schwaiger, F. W., Gostomski, D., Stear, M. J., Duncan, J. L., Mckellar, Q. A., Epplen, J. T.,et al. (1995). An-ovine major histocompatibility complex DRB1 allele is associated with low faecal egg counts following natural, predominantly Ostertagia circumcincta infection. Int. J. Parasitol. 25, 815–822. doi: 10.1016/0020-7519(94)00216-B

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Stear, M. J., Bairden, K., Bishop, S. C., Buitkamp, J., Duncan, J. L., Gettinby, G.,et al. (1997). The genetic basis of resistance to Ostertagia circumcincta in lambs. Vet. J. 154, 111–119. doi: 10.1016/S1090-0233(97)80049-4

CrossRef Full Text | Google Scholar

Stear, M. J., Doligalska, M., and Donskow Schmelter, K. (2007). Alternatives to anthelmintics for the control of nematodes in livestock. Parasitology 134, 139–151. doi: 10.1017/S0031182006001557

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Stear, M. J., Park, M., and Bishop, S. C. (1996). The key components of resistance to Ostertagia circumcincta in lambs. Parasitol. Today 12, 438–441. doi: 10.1016/0169-4758(96)10069-7

CrossRef Full Text | Google Scholar

Tavakolian, J. (2000). An Introduction to Genetic Resources of Native Farm Animals in Iran. Tehran: Animal Science Genetic Research Institute Press.

Google Scholar

Vatta, A. F., Letty, B. A., van der Linde, M. J., van Wijk, E. F., Hansen, J. W., and Krecek, R. C. (2001). Testing for clinical anaemia caused by Haemonchus spp. in goats farmed under resource-poor conditions in South Africa using an eye colour chart developed for sheep. Vet. Parasitol. 99, 1–14. doi: 10.1016/S0304-4017(01)00446-0

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Waller, P. J. (2006). Sustainable nematode parasite control strategies for ruminant livestock by grazing management and biological control. Anim. Feed Sci. Tech. 126, 277–289. doi: 10.1016/j.anifeedsci.2005.08.007

CrossRef Full Text | Google Scholar

Waller, P. J., Schwan, O., Ljungström, B.-L., Rydzik, A., and Yeates, G. W. (2004). Evaluation of biological control of sheep parasites using Duddingtonia flagrans under commercial arming conditions on the island of Gotland, Sweden. Vet. Parasitol. 126, 299–315. doi: 10.1016/j.vetpar.2004.08.008

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Yagoob, G., Hossein, H., and Asso, F. (2013). Prevalence of abomasal nematodes in sheep slaughtered at Baneh town. Amer. J. Anim. Vet. Sci. 8, 142–145. doi: 10.3844/ajavsp.2013.142.145

CrossRef Full Text | Google Scholar

Keywords: nematodes, genetic resistance, FAMACHA, MHC-DRB1, PCR-RFLP, sheep

Citation: Valilou RH, Rafat SA, Notter DR, Shojda D, Moghaddam G and Nematollahi A (2015) Fecal egg counts for gastrointestinal nematodes are associated with a polymorphism in the MHC-DRB1 gene in the Iranian Ghezel sheep breed. Front. Genet. 6:105. doi: 10.3389/fgene.2015.00105

Received: 16 November 2014; Accepted: 28 February 2015;
Published online: 24 March 2015

Edited by:

Stéphane Joost, École Polytechnique Fédérale de Lausanne, Switzerland

Reviewed by:

Francisco Ruiz-Fons, Spanish National Wildlife Research Institute – Spanish Research Council, Spain
Yuri Tani Utsunomiya, Universidade Estadual Paulista, Brazil

Copyright © 2015 Valilou, Rafat, Notter, Shojda, Moghaddam and Nematollahi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Seyed A. Rafat, Department of Animal Science, Faculty of Agriculture, University of Tabriz, Boulevard 29 Bahman, Tabriz 5166616471, Iran