BACKGROUND
infection by gastrointestinal nematodes  is common in ruminants worldwide, causing major economic losses due to decreased growth and milk production  <cit> .
grazing ruminants are infected by a variety of species of gin with different pathogenicities and geographical distributions  <cit> .
the control of gin in ruminants is largely based on the use of anthelmintics, combined with grazing management strategies.
however, anthelmintic resistance has appeared worldwide .
in northwest  spain, a recent survey showed that gin in  <dig>  % of the sampled flocks were resistant to at least one of the most commonly used drugs  <cit> .
the increasing prevalence of anthelmintic resistance has led to the search for alternative control methods, such as selective breeding for resistance to gin.
however, for this purpose, the identification of an appropriate method to measure resistance to infection is necessary, especially in conditions where the worm burden is low.
hence, a sensitive method for detecting infections is needed.
faecal egg counts  have been the traditional indicator trait used to assess the level of infection, based on the number of eggs per gram  of faeces, and it is related to both the worm burden and the fecundity of female adults in the host .
faecal egg counts have been used to measure genetic resistance to gin, although in natural infections they can be quite variable both within and between populations  <cit> .
however, fec are not particularly sensitive and should be interpreted in conjunction with information about the nutritional status, age and management of sheep flocks  <cit> .
as adult sheep are in general more resistant than naïve young animals, their fecs tend to be lower, adding an additional limitation to the sensitivity problem of the technique.
other phenotypes related to gin infections, such as the levels of iga in serum may be taken into account with the goal of defining resistant animals under natural conditions.
iga is a secreted antibody that plays a major role in gut infections.
animals that display high iga activity have been shown to present lower fec and shorter adult female teladorsagia circumcincta among experimentally and naturally infected sheep  <cit> .
the distribution of fec in naturally infected populations is characteristically over-dispersed within domestic and wild animals  <cit> , as well as in human populations  <cit> .
the negative binomial  distribution has been widely used to describe parasite eggs distribution.
however, when there are more zero fec values than expected, zero-inflated negative binomial  models are more appropriate  <cit> .
a zero-inflated distribution is a mixture of two distributions and can arise if some animals with zero egg counts have been exposed and are resistant to the infection while other animals with zero egg counts have not been exposed or recently infected e.g.
no established worms since the last anthelmintic treatment.
resistant animals tend to have few parasite eggs in their faeces.
due to the mcmaster measurement technique, small egg numbers are difficult to detect and will be counted as zero, whether the animal has really zero eggs or just a small number of them.
we hypothesize that by exploiting additional information, such as that provided by parasite-specific iga activity, we could improve the ability to discriminate animals with low level of infection with zero egg counts from unexposed/recently uninfected animals.
therefore, the objective of the study was to determine the prevalence of gin infections in naturally infected adult sheep showing low levels of infection by combining information from the two widely used indicator traits previously mentioned .
for this purpose, we applied a zinb model and extended it to include data from iga responses.
for the subset of animals that were considered as exposed to the infection based on the zinb model, we calculated the correlations among the two indicator traits related to the infection by gin  and the hidden variable of animal status .
the aim was to test whether we could improve the value of mixture and enhance the utility of the zinb model in animals naturally infected with low doses of parasites.
methods
study area and animal sampling
the study was carried out in the region of castilla y león, in the nw of spain, and included  <dig> commercial dairy flocks distributed in seven out of the nine provinces of the region  .
in the study area, the flocks are reared under a semi-extensive system in which sheep graze on natural pasture for six hours per day and are kept indoors for the rest of the day.
the average size of the sampled flocks was  <dig>  ranging from  <dig> to  <dig> animals per flock.fig.
<dig> map of the region of castilla y leon .
the map shows the location of the farms where the flocks were sampled.
map created in r using data from www.gadm.org
the survey was conducted from december  <dig> to june  <dig>  this period was extremely dry .
two conditions had to be met to include a flock in the study: first, the last anthelmintic treatment must have been administered at least 2 months before collecting the samples, and second, the sheep had to be grazing at the time of sampling.
the animals included in this study were ewes obtained by artificial insemination from farms belonging to the selection nucleus of the national association of churra breeders .
moreover, these animals were a subset of those previously genotyped with the illumina ovinesnp <dig> beadchip by  <cit>  which were still alive during the sampling period and for which both phenotypes related to parasite resistance were available.
faecal samples were collected for each ewe directly from the rectum and blood samples were obtained by venipuncture of the jugular vein.
serum samples were stored at -20 °c until processing.
this study is based on  <dig> adult churra sheep with faecal and blood serum samples.
the mean number of sheep sampled per flock was  <dig> .
the age of the sheep included in the study varied between  <dig> and 11 years.
all of the sheep were undergoing milking at the time of sampling and were experiencing at least their third lactation.
parasitological measures
a modified mcmaster technique  <cit>  using zinc sulphate as a flotation solution was used to determine the number of eggs in faeces.
the minimum detection limit of this technique was  <dig> eggs per gram .
faecal egg counts were determined by multiplying the number of eggs observed microscopically  by  <dig>
in each flock, pooled faeces were cultured to recover and identify third-stage larvae  following standard parasitological techniques  <cit> .
a total of 100 l <dig> were identified per flock to estimate the percentage of each species.
titre of iga
an indirect elisa was carried out to determine the activity of iga in the serum, results were scored as optical density .
the preparation of somatic antigen from fourth-stage larvae  of t. circumcincta was conducted as previously described by  <cit> .
microtitre plates  were coated with 100 μl of pbs containing  <dig>  μg/ml of t. circumcincta l <dig> somatic antigen, after which the plates were stored overnight at 4 °c.
after discarding their contents, the plates were blocked with 250 μl of pt-milk  for 30 min at 37 °c.
then, the blocking buffer was discarded, and 100 μl of serum was added, followed by incubation for 30 min at 37 °c.
after washing the plates four times with pbs-tween, 100 μl of a rabbit anti-sheep iga antibody, conjugated to horseradish peroxidase , at a dilution of 1/ <dig> in pt-milk, was added, followed by incubation for 30 min at 37 °c.
the plates were then washed again four times with pbs-tween and subsequently incubated in a peroxidase substrate and tetramethylbenzidine solution to produce a colour reaction, which was stopped by the addition of 50 μl of 2 m h2so <dig>  finally, the absorbance was measured at 450 nm in a microplate reader .
positive and negative controls were included in every plate.
positive controls were obtained from a pool of serum from experimentally infected sheep with t. circumcincta and negative controls from non-infected sheep that were kept indoors.
the results were expressed as the optical density ratio :  <dig> odr=sampleod‐negativeod/positiveod‐negativeod
descriptive statistics
descriptive statistical analysis for the two traits was conducted for the  <dig> sampled animals with the ‘pastecs’ package  <cit>  in r  <cit> .
the shapiro-wilk test was carried out to determine if the data for each trait was normally distributed.
due to the large number of zero counts in the fec data and the fact that the animals graze during short periods of time , we decided to use a zinb model to estimate the zero-inflation parameter and then extended it to discriminate between exposed and unexposed animals.
the zero inflated model with iga data was compared to a simpler negative binomial model using a likelihood ratio test.
moreover, in this particular study, a zero inflated model is a biologically meaningful description of the system; the adverse climatic conditions for larval development of the year studied will reduce pasture contamination, and the short grazing periods due to the semi-extensive rearing system will reduce exposure, which means that some animals would not have been infected at the time of sampling, and may not have been infected since the last anthelmintic treatment.
the zero inflated model also allows for a more natural extension into discriminating between infected and uninfected animals.
estimation of zero-inflation
in the zero inflated model, positive fec are derived from a nb distribution, while a zero count can arise from either the nb distribution or the zero distribution .
the probability of belonging to the zero distribution is called the zero-inflation parameter.
the animals that have zero counts arising from the zero distribution are assumed to have not been infected since the last anthelmintic treatment, so these animals can be excluded from further analysis.
a markov chain monte carlo model similar to the one described in denwood et al.
<cit>  using the ‘runjags’ package  <cit>  was employed to estimate the zero-inflation parameter.
in this model, the negative binomial distribution arises from a gamma-poisson mixture distribution.
uninformative priors were used for the parameters of the gamma distribution.
the posterior distribution of the zero-inflation parameter is shown in fig.
<dig> fig.
<dig> posterior distribution obtained from the extended zinb model.
each colour represents a different chain.
both chains have a mean around  <dig>  and no sample was recovered from either of the chains with a zero-inflation parameter equal to zero
extending the zinb model
a zero-inflation model does not determine which animals are exposed and resistant .
the classical zinb model was therefore extended to accommodate iga data as additional information for the animal status, i.e.
infected or not recently infected.
the animal status is calculated as,  <dig> status={0;notrecentlyinfectedwithprobability1−pexp,1;infectedwithprobabilitypexp
where status =  <dig> means that the animal has not been recently infected and status =  <dig> means that the animal is infected.
p is the probability of being recently exposed and is equivalent to one minus the zero-inflation parameter.
the raw egg counts  were used and it is assumed that for each animal i, the number of eggs counted arises from the following,  <dig> neggsi∼{0ifstatus= <dig> poissonifstatus= <dig>
where λi is the number of eggs arising from the gamma distribution .
<dig> λi∼gammashape,rate
with the shape and the rate parameters of the gamma being calculated by the model.
similarly the iga data can be partitioned in  <dig> gamma distributions  based on the animal status.
<dig> igai~gammash1rt1ifstatus= <dig> gammash2rt2ifstatus= <dig>
with sh <dig>  sh <dig>  rt <dig> and rt <dig> being the two shapes and two rates respectively that parametrize the two gamma distributions.
in the model, samples are drawn for sh <dig> and sh <dig> as well as for mn <dig> and mn <dig>  which are the two means of the two gamma distributions.
the rates are calculated by rate = shape/mean and the mean for the animals not recently infected  is always smaller than the mean of the infected .
the fully annotated r code of the model is given in the additional file  <dig>
the number of iterations sampled was  <dig> , with the first  <dig>  being discarded , and assessed convergence with the gelman-rubin statistic from the ‘coda’ package  <cit>  being under  <dig> .
using the realisations of the animal status across the iterations , it is possible to calculate the probability for each animal to be in one status or the other, piexp; animals without zero fec will always be in the infected status.
the animals that were estimated to be unexposed, i.e.
the animals with status =  <dig>  in each sample of the markov chain were excluded from further analyses, allowing the use of simple statistical tools to analyse the remaining dataset for each sample.
correlations between phenotypes
considering fec, iga and the realisations of animal status, piexp, the kendall’s rank correlation coefficient was used to estimate the relationships among these three parameters.
we used kendall’s rank because it is an appropriate non-parametric hypothesis test.
correlations were calculated in r, using the ‘ltm’ package  <cit> , for each sample of the markov chain and the average across the samples is reported below.
RESULTS
descriptive statistics of the phenotypic data
faecal egg counts and larval identification: faecal egg counts of gin ranged from  <dig> to  <dig>  epg.
in 64 % of the faecal samples no eggs were detected.
the fec mean and total variance were  <dig>   and  <dig> .
<dig> respectively.
the fec distribution was heavily skewed to the right and showed a high level of over-dispersion .
the shapiro-wilk test for the fec data indicated a clear deviation from normality .
most of the eggs detected in positive samples were strongyle-type.fig.
<dig> distribution of  faecal egg counts and  plasma iga across the  <dig> spanish adult churra ewes
apart from the gin eggs, other parasite eggs were detected in faeces:  <dig>  % of the sheep sampled had d. dendriticum eggs, with a range of 0– <dig>  epg;  <dig>  % had trichuris spp.
eggs , two animals , had moniezia spp.
eggs  and one ewe had capillaria spp eggs at a concentration of  <dig> eggs per gram.
after collecting l <dig> from coprocultures, we identified the following genera of gin: trichostrongylus spp.
, t. circumcincta , nematodirus spp.
and cooperia spp.
.
in all flocks, we confirmed the presence of t. circumcincta.
we also observed a number of lungworm larvae, though they were not identified to the species level.
iga activity in the serum samples: for individual animals, the mean odr was  <dig>  , showing a range between  <dig>  and  <dig> ; the odr variance was  <dig> .
the distribution of iga activities was positively skewed  with most of the sheep displaying relatively low iga values, and only a few sheep presenting particularly high levels of iga.
the shapiro-wilk test indicated a clear deviation from the normality .
the kolmogorov-smirnov test indicated that the iga was not gamma distributed , however this is due to the long tail of high iga values.
if the analysis is done with  <dig> animals less , the test indicates that the data is indeed gamma distributed .
zero-inflation parameter and extension of the zinb model for fec data
to verify that the data is zero inflated, a likelihood ratio test was performed comparing the zinb model to a simpler nb model, with a p-value of the likelihood ratio test =  <dig>  x10- <dig>  which indicates that the zero-inflated model provides a better fit to the data.
the mean of the zero-inflation parameter was  <dig> , this indicates that on average, 38 % of all the animals were not exposed and infected since the last anthelmintic treatment , therefore it was estimated that  <dig> ewes were infected at sampling, even though only  <dig> had non-zero fec.
the zero-inflation parameter credible interval was much narrower when using the extended zinb model as opposed to the zinb model using fec data only .
the distribution of the probability of being exposed across all the animals in the data is shown in fig.
<dig> fig.
<dig> probability of being exposed and infected, piexp, for the  <dig> animals sampled, which is calculated from the realisations of animal status  across the samples of the mcmc chain
associations between phenotypes
the associations between phenotypes was calculated for the subset of animals that were considered exposed to the infection based on the implementation of the extended zinb model  in each sample of the markov chain.
the correlations between neggs, iga and the estimated probability of being exposed to infection  are shown in table  <dig>  the phenotypic correlation between plasma iga and number of eggs was close to zero and not statistically significant, while animal status was positively correlated to the number of eggs and iga.table  <dig> estimated correlations in the churra sheep population
piexp
piexp
neggs is the number of eggs counted, iga is the activity of iga in serum  and piexp is the probability of being exposed
**p <  <dig>
discussion
adult female sheep play a key role in the epidemiology of gin infection because eggs deposited during the periparturient period influence the severity of the infection during the grazing season.
however, outside the periparturient period, egg counts in adult sheep are typically low  <cit> .
in general, gin populations in naturally infected sheep are usually over-dispersed, with the majority of sheep showing low epg values and only a few sheep presenting a high level of infection  <cit> .
in addition, some infected sheep will have low egg counts  <cit> .
therefore, supplementary information is needed as well as egg counts to determine which sheep are infected in adult sheep flocks.
in this study, the mean fec per flock was quite low  compared with other studies carried out in the same area .
gutiérrez-gil et al.
<cit>  reported that the mean fec was  <dig> epg between the years  <dig> and  <dig>  similar records were described by martínez-valladares et al.
<cit> , who showed that the prevalence of gin, based solely on the presence or absence of fec, in sheep flocks was 100 %, and the mean epg was  <dig>   between the years  <dig> and  <dig>  in the current study, the low levels of infection are likely a consequence of the exceptional climatic conditions during this study since the longevity of infective trichostrongylid l <dig> nematodes is related to temperature and humidity  <cit> .
the table in additional file  <dig> displays the mean temperature and precipitation for the period between december-june of the last 5 years  in the region of castilla y león, highlighting the fact that the year 2011/ <dig> was extremely dry.
according to martínez-valladares et al.
<cit> , there is a direct relationship between gin infection levels and the humidity of ambient air.
faecal egg count, which has been for many years the traditional diagnostic tool for assessing gin infection, has a low sensitivity  <cit> , especially for very low counts as is the case in this study.
therefore, when the excretion of eggs in faeces is low, it is necessary to use other, more sensitive, diagnostic methods that might provide a more reliable indicator of infection.
iga activity in the current study is moderately high, and this is presumed to be due to the fact that the antibodies persist for some time after gin infection.
the experimental studies of different breeds of sheep infected with gin showed iga activity for prolonged periods of time post infection.
in an experiment carried out by henderson and stear  <cit> , the peak of iga was at 6– <dig> days after a deliberate infection with t. circumcincta in sheep although detectable iga was evident 6 weeks later.
furthermore in an experiment with churra sheep, martínez-valladares et al.
<cit>  also showed that the elevated level of iga in blood and nasal secretions was maintained 4 weeks post infection with this same parasite species.
in the study by mackinnon et al.
<cit>  iga activity was also evident 4 weeks post infection with haemonchus contortus in caribbean hair sheep.
in this study, a zinb model was used to calculate the extent of zero inflation.
this approach has been applied to several parasitic infections  <cit> .
this model was then extended to identify the animals that are likely to be uninfected.
this was done by adding the iga information to the model.
in a zinb model using only fec data, the model would not be able to assign animals with zero fec as infected or uninfected.
.
there is heterogeneity among animals in the intensity of infection.
some infected animals will be exposed to more parasites than others.
both genetic variation in resistance and variation in exposure will contribute to the observed variation in iga activity and fec in exposed animals.
among animals that have not been exposed to parasites, fec will be zero and parasite-specific iga will be very low or zero.
animals with zero fec and zero or low iga activity are therefore more likely to be unexposed but it is possible that some of these animals have been exposed to low intensities of infection.
therefore the extension of the zinb model to include additional data does not guarantee that every animal will be correctly assigned.
it does however significantly improve the discrimination between exposed and unexposed animals and make subsequent analyses based on exposed animals more reliable .
to our knowledge, this is the first description of a zinb model for the analysis of multiple traits with the aim of discerning which animals are infected and which have not been recently exposed or exposed to a very low infection level.
this procedure is relatively straightforward and allows the study of nematode infections in adult animals and in flocks with low prevalence of infection, such as in mediterranean dairy farms where animals are under a semi-extensive management system.
the approach improves our ability to identify animals that have been infected with gin, even at low fec, which is needed for the study of host resistance in naturally infected individuals and the breeding of resistant sheep.
because the over-dispersion pattern of gin  is also observed in other hosts such as cattle, free-range pigs, chickens, humans and wild animals , the approach described here could also be useful in other systems.
the correlations between the number of eggs and iga and animal status were calculated using the non-parametric kendall’s test.
although the number of eggs has been found negatively correlated with iga in young lambs  <cit> , in the case of adult sheep, this correlation is not as clear and both coltman et al.
<cit>  and gutiérrez-gil et al.
<cit>  reported non-significant correlations in naturally infected adult sheep after comparing logfec and iga against somatic antigen from t. circumcinta l <dig>  our results are similar, and suggest that this correlation is indeed close to zero in adult sheep.
in experimentally infected adult sheep, martinez-valladares et al.
<cit>  showed negative correlations between iga in gastric mucus and fec whereas the correlation between fec and the serum iga levels  were not significant.
the absence of a clear correlation between plasma iga and fec may be a consequence of the fact that plasma iga shows a complex relationship with mucosal iga  <cit> .
alternatively, adult sheep may show greater ige activity; reduced numbers of established parasites would decrease iga responses and the relative importance of iga on egg output would be lowered  <cit> .
the extension of the zinb model has allowed us to combine the information from two different traits that can indicate resistance or susceptibility to gin.
the iga response was added to the model to help discriminate between unexposed and infected animals with zero fec.
recent research has produced an index of the intensity of nematode infection in young lambs  <cit>  and the observed correlations among the parasitological variable are necessary for this process.
as mentioned previously, the use of a reliable indicator trait may be of interest not only for the management of parasite infections but also for the design of breeding programs aimed at achieving resistance to parasites.
CONCLUSIONS
in summary, in the current study, two different phenotypes related to gin infection  were analysed.
there was a high percentage of sheep without eggs in faeces  and a zero inflated model was used to detect the amount of zero inflation in the data.
the zinb model suggested that 38 % of sampled sheep had not been exposed to nematode infection in the previous 2 months, since the last anthelmintic treatment.
therefore, in addition to fec data, the evaluation of iga in serum may help to distinguish adult animals with low level of infection from resistant animals assist selective breeding for resistance to gin.
abbreviations
anche, national association of churra breeders ; fec, faecal egg counts; gin, gastrointestinal nematodes; iga, immunoglobulin a; od/odr, optical density ; zinb/nb,  negative binomial.
additional files
additional file 1: mean temperatures  and precipitation  from december to june during the sampling period , and during the four previous years.
additional file 2: annotated r code for the zinb model.
additional file 3: exposed probability in a classic zinb model.
histogram of probabilities of being exposed for the data  and zoom of only the zero fec  using only fec data in the zinb model.
animals with non-zero fec will always have an “infected” status in the model  while animals with zero fec can be exposed or unexposed.
if only the fec data is used, each animal with zero fec will have a probability of being infected similar to one minus the zero-inflation parameter .
additional file 4: raw data used.

 infection by gastrointestinal nematodes is common in ruminants worldwide, causing major economic losses due to decreased growth and milk production. grazing ruminants are infected by a variety of species of gin with different pathogenicities and geographical distributions. the control of gin in ruminants is largely based on the use of anthelmintics, combined with grazing management strategies. however, anthelmintic resistance has appeared worldwide. in northwest spain, a recent survey showed that gin in  % of the sampled flocks were resistant to at least one of the most commonly used drugs. the increasing prevalence of anthelmintic resistance has led to the search for alternative control methods, such as selective breeding for resistance to gin. however, for this purpose, the identification of an appropriate method to measure resistance to infection is necessary, especially in conditions where the worm burden is low. hence, a sensitive method for detecting infections is needed. faecal egg counts have been the traditional indicator trait used to assess the level of infection, based on the number of eggs per gram of faeces, and it is related to both the worm burden and the fecundity of female adults in the host. faecal egg counts have been used to measure genetic resistance to gin, although in natural infections they can be quite variable both within and between populations. however, fec are not particularly sensitive and should be interpreted in conjunction with information about the nutritional status, age and management of sheep flocks. as adult sheep are in general more resistant than naïve young animals, their fecs tend to be lower, adding an additional limitation to the sensitivity problem of the technique. other phenotypes related to gin infections, such as the levels of iga in serum may be taken into account with the goal of defining resistant animals under natural conditions. iga is a secreted antibody that plays a major role in gut infections. animals that display high iga activity have been shown to present lower fec and shorter adult female teladorsagia circumcincta among experimentally and naturally infected sheep. the distribution of fec in naturally infected populations is characteristically over-dispersed within domestic and wild animals, as well as in human populations. the negative binomial distribution has been widely used to describe parasite eggs distribution. however, when there are more zero fec values than expected, zero-inflated negative binomial models are more appropriate. a zero-inflated distribution is a mixture of two distributions and can arise if some animals with zero egg counts have been exposed and are resistant to the infection while other animals with zero egg counts have not been exposed or recently infected e.g. no established worms since the last anthelmintic treatment. resistant animals tend to have few parasite eggs in their faeces. due to the mcmaster measurement technique, small egg numbers are difficult to detect and will be counted as zero, whether the animal has really zero eggs or just a small number of them. we hypothesize that by exploiting additional information, such as that provided by parasite-specific iga activity, we could improve the ability to discriminate animals with low level of infection with zero egg counts from unexposed/recently uninfected animals. therefore, the objective of the study was to determine the prevalence of gin infections in naturally infected adult sheep showing low levels of infection by combining information from the two widely used indicator traits previously mentioned. for this purpose, we applied a zinb model and extended it to include data from iga responses. for the subset of animals that were considered as exposed to the infection based on the zinb model, we calculated the correlations among the two indicator traits related to the infection by gin and the hidden variable of animal status. the aim was to test whether we could improve the value of mixture and enhance the utility of the zinb model in animals naturally infected with low doses of parasites. methods study area and animal sampling the study was carried out in the region of castilla y león, in the nw of spain, and included commercial dairy flocks distributed in seven out of the nine provinces of the region. in the study area, the flocks are reared under a semi-extensive system in which sheep graze on natural pasture for six hours per day and are kept indoors for the rest of the day. the average size of the sampled flocks was ranging from to animals per flock.fig. map of the region of castilla y leon. the map shows the location of the farms where the flocks were sampled. map created in r using data from www.gadm.org the survey was conducted from december to june this period was extremely dry. two conditions had to be met to include a flock in the study: first, the last anthelmintic treatment must have been administered at least 2 months before collecting the samples, and second, the sheep had to be grazing at the time of sampling. the animals included in this study were ewes obtained by artificial insemination from farms belonging to the selection nucleus of the national association of churra breeders. moreover, these animals were a subset of those previously genotyped with the illumina ovinesnp beadchip by which were still alive during the sampling period and for which both phenotypes related to parasite resistance were available. faecal samples were collected for each ewe directly from the rectum and blood samples were obtained by venipuncture of the jugular vein. serum samples were stored at -20 °c until processing. this study is based on adult churra sheep with faecal and blood serum samples. the mean number of sheep sampled per flock was. the age of the sheep included in the study varied between and 11 years. all of the sheep were undergoing milking at the time of sampling and were experiencing at least their third lactation. parasitological measures a modified mcmaster technique using zinc sulphate as a flotation solution was used to determine the number of eggs in faeces. the minimum detection limit of this technique was eggs per gram. faecal egg counts were determined by multiplying the number of eggs observed microscopically by in each flock, pooled faeces were cultured to recover and identify third-stage larvae following standard parasitological techniques. a total of 100 l were identified per flock to estimate the percentage of each species. titre of iga an indirect elisa was carried out to determine the activity of iga in the serum, results were scored as optical density. the preparation of somatic antigen from fourth-stage larvae of t. circumcincta was conducted as previously described by. microtitre plates were coated with 100 μl of pbs containing  μg/ml of t. circumcincta l somatic antigen, after which the plates were stored overnight at 4 °c. after discarding their contents, the plates were blocked with 250 μl of pt-milk for 30 min at 37 °c. then, the blocking buffer was discarded, and 100 μl of serum was added, followed by incubation for 30 min at 37 °c. after washing the plates four times with pbs-tween, 100 μl of a rabbit anti-sheep iga antibody, conjugated to horseradish peroxidase, at a dilution of 1/ in pt-milk, was added, followed by incubation for 30 min at 37 °c. the plates were then washed again four times with pbs-tween and subsequently incubated in a peroxidase substrate and tetramethylbenzidine solution to produce a colour reaction, which was stopped by the addition of 50 μl of 2 m h2so finally, the absorbance was measured at 450 nm in a microplate reader. positive and negative controls were included in every plate. positive controls were obtained from a pool of serum from experimentally infected sheep with t. circumcincta and negative controls from non-infected sheep that were kept indoors. the results were expressed as the optical density ratio: odr=sampleod‐negativeod/positiveod‐negativeod descriptive statistics descriptive statistical analysis for the two traits was conducted for the sampled animals with the ‘pastecs’ package in r. the shapiro-wilk test was carried out to determine if the data for each trait was normally distributed. due to the large number of zero counts in the fec data and the fact that the animals graze during short periods of time, we decided to use a zinb model to estimate the zero-inflation parameter and then extended it to discriminate between exposed and unexposed animals. the zero inflated model with iga data was compared to a simpler negative binomial model using a likelihood ratio test. moreover, in this particular study, a zero inflated model is a biologically meaningful description of the system; the adverse climatic conditions for larval development of the year studied will reduce pasture contamination, and the short grazing periods due to the semi-extensive rearing system will reduce exposure, which means that some animals would not have been infected at the time of sampling, and may not have been infected since the last anthelmintic treatment. the zero inflated model also allows for a more natural extension into discriminating between infected and uninfected animals. estimation of zero-inflation in the zero inflated model, positive fec are derived from a nb distribution, while a zero count can arise from either the nb distribution or the zero distribution. the probability of belonging to the zero distribution is called the zero-inflation parameter. the animals that have zero counts arising from the zero distribution are assumed to have not been infected since the last anthelmintic treatment, so these animals can be excluded from further analysis. a markov chain monte carlo model similar to the one described in denwood et al. using the ‘runjags’ package was employed to estimate the zero-inflation parameter. in this model, the negative binomial distribution arises from a gamma-poisson mixture distribution. uninformative priors were used for the parameters of the gamma distribution. the posterior distribution of the zero-inflation parameter is shown in fig.  fig. posterior distribution obtained from the extended zinb model. each colour represents a different chain. both chains have a mean around and no sample was recovered from either of the chains with a zero-inflation parameter equal to zero extending the zinb model a zero-inflation model does not determine which animals are exposed and resistant. the classical zinb model was therefore extended to accommodate iga data as additional information for the animal status, i.e. infected or not recently infected. the animal status is calculated as, status={0;notrecentlyinfectedwithprobability1−pexp,1;infectedwithprobabilitypexp where status =  means that the animal has not been recently infected and status =  means that the animal is infected. p is the probability of being recently exposed and is equivalent to one minus the zero-inflation parameter. the raw egg counts were used and it is assumed that for each animal i, the number of eggs counted arises from the following, neggsi∼{0ifstatus= poissonifstatus= where λi is the number of eggs arising from the gamma distribution. λi∼gammashape,rate with the shape and the rate parameters of the gamma being calculated by the model. similarly the iga data can be partitioned in gamma distributions based on the animal status. igai~gammash1rt1ifstatus= gammash2rt2ifstatus= with sh sh rt and rt being the two shapes and two rates respectively that parametrize the two gamma distributions. in the model, samples are drawn for sh and sh as well as for mn and mn which are the two means of the two gamma distributions. the rates are calculated by rate = shape/mean and the mean for the animals not recently infected is always smaller than the mean of the infected. the fully annotated r code of the model is given in the additional file the number of iterations sampled was, with the first being discarded, and assessed convergence with the gelman-rubin statistic from the ‘coda’ package being under. using the realisations of the animal status across the iterations, it is possible to calculate the probability for each animal to be in one status or the other, piexp; animals without zero fec will always be in the infected status. the animals that were estimated to be unexposed, i.e. the animals with status =  in each sample of the markov chain were excluded from further analyses, allowing the use of simple statistical tools to analyse the remaining dataset for each sample. correlations between phenotypes considering fec, iga and the realisations of animal status, piexp, the kendall’s rank correlation coefficient was used to estimate the relationships among these three parameters. we used kendall’s rank because it is an appropriate non-parametric hypothesis test. correlations were calculated in r, using the ‘ltm’ package, for each sample of the markov chain and the average across the samples is reported below.  descriptive statistics of the phenotypic data faecal egg counts and larval identification: faecal egg counts of gin ranged from to epg. in 64 % of the faecal samples no eggs were detected. the fec mean and total variance were and. respectively. the fec distribution was heavily skewed to the right and showed a high level of over-dispersion. the shapiro-wilk test for the fec data indicated a clear deviation from normality. most of the eggs detected in positive samples were strongyle-type.fig. distribution of faecal egg counts and plasma iga across the spanish adult churra ewes apart from the gin eggs, other parasite eggs were detected in faeces:  % of the sheep sampled had d. dendriticum eggs, with a range of 0– epg;  % had trichuris spp. eggs, two animals, had moniezia spp. eggs and one ewe had capillaria spp eggs at a concentration of eggs per gram. after collecting l from coprocultures, we identified the following genera of gin: trichostrongylus spp., t. circumcincta, nematodirus spp. and cooperia spp.. in all flocks, we confirmed the presence of t. circumcincta. we also observed a number of lungworm larvae, though they were not identified to the species level. iga activity in the serum samples: for individual animals, the mean odr was, showing a range between and; the odr variance was. the distribution of iga activities was positively skewed with most of the sheep displaying relatively low iga values, and only a few sheep presenting particularly high levels of iga. the shapiro-wilk test indicated a clear deviation from the normality. the kolmogorov-smirnov test indicated that the iga was not gamma distributed, however this is due to the long tail of high iga values. if the analysis is done with animals less, the test indicates that the data is indeed gamma distributed. zero-inflation parameter and extension of the zinb model for fec data to verify that the data is zero inflated, a likelihood ratio test was performed comparing the zinb model to a simpler nb model, with a p-value of the likelihood ratio test =  x10- which indicates that the zero-inflated model provides a better fit to the data. the mean of the zero-inflation parameter was, this indicates that on average, 38 % of all the animals were not exposed and infected since the last anthelmintic treatment, therefore it was estimated that ewes were infected at sampling, even though only had non-zero fec. the zero-inflation parameter credible interval was much narrower when using the extended zinb model as opposed to the zinb model using fec data only. the distribution of the probability of being exposed across all the animals in the data is shown in fig.  fig. probability of being exposed and infected, piexp, for the animals sampled, which is calculated from the realisations of animal status across the samples of the mcmc chain associations between phenotypes the associations between phenotypes was calculated for the subset of animals that were considered exposed to the infection based on the implementation of the extended zinb model in each sample of the markov chain. the correlations between neggs, iga and the estimated probability of being exposed to infection are shown in table  the phenotypic correlation between plasma iga and number of eggs was close to zero and not statistically significant, while animal status was positively correlated to the number of eggs and iga.table estimated correlations in the churra sheep population piexp piexp neggs is the number of eggs counted, iga is the activity of iga in serum and piexp is the probability of being exposed **p <  discussion adult female sheep play a key role in the epidemiology of gin infection because eggs deposited during the periparturient period influence the severity of the infection during the grazing season. however, outside the periparturient period, egg counts in adult sheep are typically low. in general, gin populations in naturally infected sheep are usually over-dispersed, with the majority of sheep showing low epg values and only a few sheep presenting a high level of infection. in addition, some infected sheep will have low egg counts. therefore, supplementary information is needed as well as egg counts to determine which sheep are infected in adult sheep flocks. in this study, the mean fec per flock was quite low compared with other studies carried out in the same area. gutiérrez-gil et al. reported that the mean fec was epg between the years and similar records were described by martínez-valladares et al., who showed that the prevalence of gin, based solely on the presence or absence of fec, in sheep flocks was 100 %, and the mean epg was between the years and in the current study, the low levels of infection are likely a consequence of the exceptional climatic conditions during this study since the longevity of infective trichostrongylid l nematodes is related to temperature and humidity. the table in additional file displays the mean temperature and precipitation for the period between december-june of the last 5 years in the region of castilla y león, highlighting the fact that the year 2011/ was extremely dry. according to martínez-valladares et al., there is a direct relationship between gin infection levels and the humidity of ambient air. faecal egg count, which has been for many years the traditional diagnostic tool for assessing gin infection, has a low sensitivity, especially for very low counts as is the case in this study. therefore, when the excretion of eggs in faeces is low, it is necessary to use other, more sensitive, diagnostic methods that might provide a more reliable indicator of infection. iga activity in the current study is moderately high, and this is presumed to be due to the fact that the antibodies persist for some time after gin infection. the experimental studies of different breeds of sheep infected with gin showed iga activity for prolonged periods of time post infection. in an experiment carried out by henderson and stear, the peak of iga was at 6– days after a deliberate infection with t. circumcincta in sheep although detectable iga was evident 6 weeks later. furthermore in an experiment with churra sheep, martínez-valladares et al. also showed that the elevated level of iga in blood and nasal secretions was maintained 4 weeks post infection with this same parasite species. in the study by mackinnon et al. iga activity was also evident 4 weeks post infection with haemonchus contortus in caribbean hair sheep. in this study, a zinb model was used to calculate the extent of zero inflation. this approach has been applied to several parasitic infections. this model was then extended to identify the animals that are likely to be uninfected. this was done by adding the iga information to the model. in a zinb model using only fec data, the model would not be able to assign animals with zero fec as infected or uninfected.. there is heterogeneity among animals in the intensity of infection. some infected animals will be exposed to more parasites than others. both genetic variation in resistance and variation in exposure will contribute to the observed variation in iga activity and fec in exposed animals. among animals that have not been exposed to parasites, fec will be zero and parasite-specific iga will be very low or zero. animals with zero fec and zero or low iga activity are therefore more likely to be unexposed but it is possible that some of these animals have been exposed to low intensities of infection. therefore the extension of the zinb model to include additional data does not guarantee that every animal will be correctly assigned. it does however significantly improve the discrimination between exposed and unexposed animals and make subsequent analyses based on exposed animals more reliable. to our knowledge, this is the first description of a zinb model for the analysis of multiple traits with the aim of discerning which animals are infected and which have not been recently exposed or exposed to a very low infection level. this procedure is relatively straightforward and allows the study of nematode infections in adult animals and in flocks with low prevalence of infection, such as in mediterranean dairy farms where animals are under a semi-extensive management system. the approach improves our ability to identify animals that have been infected with gin, even at low fec, which is needed for the study of host resistance in naturally infected individuals and the breeding of resistant sheep. because the over-dispersion pattern of gin is also observed in other hosts such as cattle, free-range pigs, chickens, humans and wild animals, the approach described here could also be useful in other systems. the correlations between the number of eggs and iga and animal status were calculated using the non-parametric kendall’s test. although the number of eggs has been found negatively correlated with iga in young lambs, in the case of adult sheep, this correlation is not as clear and both coltman et al. and gutiérrez-gil et al. reported non-significant correlations in naturally infected adult sheep after comparing logfec and iga against somatic antigen from t. circumcinta l our results are similar, and suggest that this correlation is indeed close to zero in adult sheep. in experimentally infected adult sheep, martinez-valladares et al. showed negative correlations between iga in gastric mucus and fec whereas the correlation between fec and the serum iga levels were not significant. the absence of a clear correlation between plasma iga and fec may be a consequence of the fact that plasma iga shows a complex relationship with mucosal iga. alternatively, adult sheep may show greater ige activity; reduced numbers of established parasites would decrease iga responses and the relative importance of iga on egg output would be lowered. the extension of the zinb model has allowed us to combine the information from two different traits that can indicate resistance or susceptibility to gin. the iga response was added to the model to help discriminate between unexposed and infected animals with zero fec. recent research has produced an index of the intensity of nematode infection in young lambs and the observed correlations among the parasitological variable are necessary for this process. as mentioned previously, the use of a reliable indicator trait may be of interest not only for the management of parasite infections but also for the design of breeding programs aimed at achieving resistance to parasites.  in summary, in the current study, two different phenotypes related to gin infection were analysed. there was a high percentage of sheep without eggs in faeces and a zero inflated model was used to detect the amount of zero inflation in the data. the zinb model suggested that 38 % of sampled sheep had not been exposed to nematode infection in the previous 2 months, since the last anthelmintic treatment. therefore, in addition to fec data, the evaluation of iga in serum may help to distinguish adult animals with low level of infection from resistant animals assist selective breeding for resistance to gin. abbreviations anche, national association of churra breeders; fec, faecal egg counts; gin, gastrointestinal nematodes; iga, immunoglobulin a; od/odr, optical density; zinb/nb, negative binomial. additional files additional file 1: mean temperatures and precipitation from december to june during the sampling period, and during the four previous years. additional file 2: annotated r code for the zinb model. additional file 3: exposed probability in a classic zinb model. histogram of probabilities of being exposed for the data and zoom of only the zero fec using only fec data in the zinb model. animals with non-zero fec will always have an “infected” status in the model while animals with zero fec can be exposed or unexposed. if only the fec data is used, each animal with zero fec will have a probability of being infected similar to one minus the zero-inflation parameter. additional file 4: raw data used. 
the correlations calculated here will be useful in the development of a reliable index of gin resistance that could be of assistance for the study of host resistance in studies based on natural infection, especially in adult sheep, and also the design of breeding programs aimed at increasing resistance to parasites. in this dataset, 64 % of animals had zero fec while the zinb model suggested that 38 % of sheep had not been recently infected with gin. animals with zero fec and high iga activity were considered resistant while animals with zero fec and low iga activity were considered as not recently infected. the model presented here improved the detection of infected animals with zero fec. therefore 26 % of sheep were predicted to be infected animals with egg counts that were zero or below the detection limit and likely to be relatively resistant to nematode infection. iga activities of all animals were then used to decide which of the sheep with zero egg counts had been exposed and which sheep had not been recently exposed. a zero inflated negative binomial model model was used to calculate the extent of zero inflation for fec; the model was extended to include information from the iga responses. in this study, two traits related with resistance to gastrointestinal nematodes were measured in adult sheep: faecal egg count and activity of immunoglobulin a in plasma. gastrointestinal nematodessheepprevalenceegg countigazinbjunta de castilla y leonle245a12-2http://dx.doi.org/ /501100000780european commissionp7-kbbe-2009-3-245140http://dx.doi.org/ /501100000268biotechnology and biological sciences research councilbb/l004004070/1http://dx.doi.org/ /501100003329ministerio de economía y competitividadryc-2012-10230http://dx.doi.org/http://dx.doi.org/ /501100006003fundación general csicissue-copyright-statement© the author 2016
 in this study, two traits related with resistance to gastrointestinal nematodes were measured in adult sheep: faecal egg count and activity of immunoglobulin a in plasma. in dry years, fec can be very low in semi-extensive systems, such as the one studied here, which makes identifying animals that are resistant or susceptible to infection a difficult task. a zero inflated negative binomial model model was used to calculate the extent of zero inflation for fec; the model was extended to include information from the iga responses.  in this dataset, 64 % of animals had zero fec while the zinb model suggested that 38 % of sheep had not been recently infected with gin. therefore 26 % of sheep were predicted to be infected animals with egg counts that were zero or below the detection limit and likely to be relatively resistant to nematode infection. iga activities of all animals were then used to decide which of the sheep with zero egg counts had been exposed and which sheep had not been recently exposed. animals with zero fec and high iga activity were considered resistant while animals with zero fec and low iga activity were considered as not recently infected. for the animals considered as exposed to the infection, the correlations among the studied traits were estimated, and the influence of these traits on the discrimination between unexposed and infected animals was assessed.  the model presented here improved the detection of infected animals with zero fec. the correlations calculated here will be useful in the development of a reliable index of gin resistance that could be of assistance for the study of host resistance in studies based on natural infection, especially in adult sheep, and also the design of breeding programs aimed at increasing resistance to parasites. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords gastrointestinal nematodessheepprevalenceegg countigazinbjunta de castilla y leonle245a12-2http://dx.doi.org/ /501100000780european commissionp7-kbbe-2009-3-245140http://dx.doi.org/ /501100000268biotechnology and biological sciences research councilbb/l004004070/1http://dx.doi.org/ /501100003329ministerio de economía y competitividadryc-2012-10230http://dx.doi.org/http://dx.doi.org/ /501100006003fundación general csicissue-copyright-statement© the author 2016 