Hostname: page-component-7c8c6479df-7qhmt Total loading time: 0 Render date: 2024-03-26T14:18:55.827Z Has data issue: false hasContentIssue false

Seasonality of helminth infection in wild red deer varies between individuals and between parasite taxa

Published online by Cambridge University Press:  09 March 2018

Gregory F. Albery*
Affiliation:
Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JT, UK
Fiona Kenyon
Affiliation:
Moredun Research Institute, Pentlands Science Park, Bush Loan, Midlothian, EH26 0PZ, UK
Alison Morris
Affiliation:
Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JT, UK
Sean Morris
Affiliation:
Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JT, UK
Daniel H. Nussey
Affiliation:
Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JT, UK
Josephine M. Pemberton
Affiliation:
Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JT, UK
*
Author for correspondence: Gregory F. Albery, E-mail: gfalbery@gmail.com

Abstract

Parasitism in wild mammals can vary according to myriad intrinsic and extrinsic factors, many of which vary seasonally. However, seasonal variation in parasitism is rarely studied using repeated samples from known individuals. Here we used a wild population of individually recognized red deer (Cervus elaphus) on the Isle of Rum to quantify seasonality and intrinsic factors affecting gastrointestinal helminth parasitism over the course of a year. We collected 1020 non-invasive faecal samples from 328 known individuals which we then analysed for propagules of three helminth taxa: strongyle nematodes, the common liver fluke Fasciola hepatica and the tissue nematode Elaphostrongylus cervi. Zero-inflated Poisson models were used to investigate how season, age and sex were associated with parasite prevalence and count intensity, while Poisson models were used to quantify individual repeatability within and between sampling seasons. Parasite intensity and prevalence varied according to all investigated factors, with opposing seasonality, age profiles and sex biases between parasite taxa. Repeatability was moderate, decreased between seasons and varied between parasites; both F. hepatica and E. cervi showed significant between-season repeatability, while strongyle nematode counts were only repeatable within-season and showed no repeatability within individuals across the year.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © Cambridge University Press 2018

Introduction

Gastrointestinal helminths include a range of nematode, trematode and cestode species and are an important selective force in wild vertebrate populations (Poulin, Reference Poulin2007). Although they can cause severe pathology in some cases, adult helminths can survive for years within a host; infections are commonly chronic and associated with immunosuppression and minimal overt pathology (Maizels et al. Reference Maizels, Hewitson and Smith2012). The life cycles of helminth parasites vary, but broadly involve the parasite living and feeding within the host and producing propagules (eggs or larvae) into the gut lumen which are shed into the environment in the faeces, from which the parasite spreads to other hosts – sometimes via an intermediate host (Bohm et al. Reference Bohm, White, Chambers, Smith and Hutchings2007). Helminth infection is commonly quantified by counting these propagules in non-invasively collected faecal samples. Faecal egg counts (FECs) are often found to correlate well with burden (McKenna, Reference McKenna1981; Budischak et al. Reference Budischak, Jolles, Ezenwa, Hoberg, Abrams, Jolles and Ezenwa2015), despite egg output being a complex product of both host and parasite biology (Sargison, Reference Sargison2013) which can fluctuate over time (Turner et al. Reference Turner, Cizauskas and Getz2010). In both livestock and wild mammals, one striking feature of FEC is their distribution, which is typically strongly right-skewed with a small number of individuals with relatively high counts (Wilson et al. Reference Wilson, Grenfell, Pilkington, Boyd, Gulland, Clutton-Brock and Pemberton2004). A major challenge for helminth epidemiology is to determine the degree to which variation in FEC is driven by factors intrinsic to the host (e.g. age, sex and genotype) vs extrinsic factors (e.g. season, annual climate and host density). Although there is good cross-sectional evidence that FEC varies with season, age and sex in the wild, longitudinal data are essential to properly separate within- and among-individual processes, and to determine the repeatability of helminth infection within the same host. Longitudinal FEC data remain relatively scarce in wild or free-living mammals, although prior studies suggest that counts are to some degree repeatable within individuals over years and seasons (Wilson et al. Reference Wilson, Grenfell, Pilkington, Boyd, Gulland, Clutton-Brock and Pemberton2004; Wood et al. Reference Wood, Matthews, Stephenson, Slote and Nussey2013; Debeffe et al. Reference Debeffe, McLoughlin, Medill, Stewart, Andres, Shury, Wagner, Jenkins, Gilleard and Poissant2016). Importantly, as seasonal changes can affect hosts through changes in nutrition, immunity and investment in tradeoffs (Martin et al. Reference Martin, Weil and Nelson2008), parasite intensities may differ seasonally across the host population, resulting in different individuals peaking in burden or transmission at different times of the year. It is therefore of interest to compare and contrast individual repeatability of FEC across seasons and parasite groups. We aim to use repeated counts at multiple timescales to separate real between-season differences in parasitism from baseline variation arising from experimental error or small-scale temporal fluctuations.

Parasites are well known to be influenced by seasonal extrinsic factors: weather patterns can impact parasites by affecting survival and movement of transmission stages, and regular changes in host biology such as parturition or mating seasons can affect resource availability (Altizer et al. Reference Altizer, Dobson, Hosseini, Hudson, Pascual and Rohani2006). In temperate regions, helminths usually show a peak in egg output during the spring and summer, with lower or arrested transmission in the winter as a result of: (i) reduced survival and mobility of parasites (Stromberg, Reference Stromberg1997) and (ii) adaptive coincidence between infective parasites and immunologically naïve young in the spring (Wilson et al. Reference Wilson, Grenfell, Pilkington, Boyd, Gulland, Clutton-Brock and Pemberton2004). However, different species do show different adaptive seasonality: for example, in Svalbard reindeer the nematode Marshallagia marshalli transmits throughout the winter months despite extreme cold and different host feeding patterns (Carlsson et al. Reference Carlsson, Irvine, Wilson, Piertney, Halvorsen, Coulson, Stien and Albon2012). Few comparisons of seasonality in egg counts across helminth groups have been made in wild vertebrate systems to date. Studies that rely on post-mortem sampling are often seasonally restricted to windows of mortality (e.g. Craig et al. Reference Craig, Pilkington and Pemberton2006) or culling seasons (e.g. Irvine et al. Reference Irvine, Corbishley, Pilkington and Albon2006), and are therefore unable to investigate seasonal trends. Two important host factors associated with parasitism are age and sex. Age-dependent parasitism occurs across the animal kingdom, with higher intensities typically observed in the young and/or the elderly (Hayward et al. Reference Hayward, Wilson, Pilkington, Pemberton and Kruuk2009). Considering gastrointestinal helminths in wild vertebrates, there is strong evidence for declining infection intensity from birth to adulthood, which is typically ascribed to two processes: first, animals gain adaptive immunity as they age (Stear and Murray, Reference Stear and Murray1994; Turner and Getz, Reference Turner and Getz2010); second, animals that are most susceptible to parasites and therefore show the highest egg counts are more likely to die, meaning the animals that survive to old age are those with lower parasite infections (Wilson et al. Reference Wilson, Grenfell, Pilkington, Boyd, Gulland, Clutton-Brock and Pemberton2004). Increasing helminth parasitism in elderly animals has also been ascribed to immunosenescence, although the selective loss of highly infected individuals can also confound estimates of within-individual change in later life (van de Pol and Verhulst, Reference van de Pol and Verhulst2006). A particular advantage of longitudinal studies is that they allow differentiation between the within- and among-individual processes responsible for age-related variation (Clutton-Brock and Sheldon, Reference Clutton-Brock and Sheldon2010). Sex differences in parasitism are another common phenomenon in wild mammal species (Poulin, Reference Poulin1996), with higher burdens in males than females commonly observed in polygynous mammals in particular (Moore and Wilson, Reference Moore and Wilson2002). This is typically linked to sexual dimorphism and increased investment in short-term reproduction rather than immunity, rendering males the ‘sicker sex’ (Zuk, Reference Zuk2009).

Wild ungulates are commonly infected with a range of helminth parasite species and have formed the basis of many important individual-based studies in disease ecology (Jolles and Ezenwa, Reference Jolles and Ezenwa2015). Red deer have been the subject of multiple parasitological studies, due partly to their abundance in the wild where population management often requires culling (Bohm et al. Reference Bohm, White, Chambers, Smith and Hutchings2007), and as farm animals, principally in New Zealand (Mason, Reference Mason1994). They are of particular interest due to their ability to act as reservoirs and vectors of parasites that commonly infect livestock (Alberti et al. Reference Alberti, Gioia, Sironi, Zanzani, Riccaboni, Magrini and Manfredi2011; Chintoan-Uta et al. Reference Chintoan-Uta, Morgan, Skuce and Coles2014; Davidson et al. Reference Davidson, Kutz, Madslien, Hoberg and Handeland2014), and cross-sectional culling studies have formed a useful knowledge base for red deer parasitology. Previous studies of helminth parasites in wild red deer have demonstrated through cross-sectional post-mortem sampling that parasite burdens vary with age (Vicente et al. Reference Vicente, Fernández De Mera and Gortazar2006), sex (Irvine et al. Reference Irvine, Corbishley, Pilkington and Albon2006; Vicente et al. Reference Vicente, Pérez-Rodríguez and Gortazar2007a; French et al. Reference French, Zadoks, Skuce, Mitchell, Gordon-Gibbs, Craine, Shaw, Gibb and Taggart2016), host density (Vicente et al. Reference Vicente, Höfle, Fernández-De-Mera and Gortazar2007b), supplementary feeding (Hines et al. Reference Hines, Ezenwa, Cross and Rogerson2007), testosterone level (Malo et al. Reference Malo, Roldan, Garde, Soler, Vicente, Gortazar and Gomendio2009) and several measures of condition (Irvine et al. Reference Irvine, Corbishley, Pilkington and Albon2006; Vicente et al. Reference Vicente, Pérez-Rodríguez and Gortazar2007a,Reference Vicente, Höfle, Fernández-De-Mera and Gortazarb). Although such cross-sectional approaches provide detailed data on the parasite community and adult parasite burden, they cannot separate within- and among-host processes involved in helminth epidemiology and host–parasite interactions. Furthermore, male deer are culled earlier in the year than females, therefore confounding sex and seasonal effects (e.g. Irvine et al. Reference Irvine, Corbishley, Pilkington and Albon2006; French et al. Reference French, Zadoks, Skuce, Mitchell, Gordon-Gibbs, Craine, Shaw, Gibb and Taggart2016). Longitudinal studies monitoring and comparing the within-host repeatability and seasonality of counts are currently lacking for helminth parasites of wild red deer, although seasonality of E. cervi and strongyles has been demonstrated using the collection of fresh pellets though without the full benefits of individual-level data (Vicente et al. Reference Vicente, Fierro and Gortazar2005; Hines et al. Reference Hines, Ezenwa, Cross and Rogerson2007), and strongyles in red deer are known to undergo a season of arrested development (Connan, Reference Connan1997). An earlier study examining gastrointestinal parasites of red deer living across the Isle on Rum in Scotland (Irvine et al. Reference Irvine, Corbishley, Pilkington and Albon2006) used animals culled for management purposes in the late summer–autumn (males) or autumn–winter (females). It found a high prevalence, but low burden, of strongyle nematodes Ostertagia spp. and Oesophagostomum venulosum as well as the nematodes Nematodirus sp., Capillaria sp., Trichuris ovis, Elaphostrongylus cervi, Dictyocaulus sp. and coccidian Eimeria sp. The tissue nematode E. cervi and the generalist liver fluke Fasciola hepatica, both of which have life cycles that involve an intermediate snail host, have been documented using egg counts in wild red deer (Vicente et al. Reference Vicente, Pérez-Rodríguez and Gortazar2007a; French et al. Reference French, Zadoks, Skuce, Mitchell, Gordon-Gibbs, Craine, Shaw, Gibb and Taggart2016). However, strongyle nematodes in wild red deer have rarely been studied using non-invasive methods, which is surprising given their detailed study in related livestock hosts including cattle and sheep (Hoberg et al. Reference Hoberg, Kocan, Rickard, Samuel, Pybus and Kocan2001). In the present study, we repeatedly collected faeces from known study individuals of different ages and sexes on Rum within and across seasons in 2016. We analysed propagule counts of the three most prevalent helminth species groups: strongyle nematodes, E. cervi and F. hepatica. Our aims were to: (i) examine the parasite fauna of the Rum red deer and identify taxa of high prevalence for statistical analysis; (ii) assess the repeatability of non-invasive parasitological measures in the deer at multiple timescales, particularly between seasons; and (iii) investigate how prevalence and intensity of infection with the abundant taxa are associated with season, host age and host sex.

Methods

Study area and sample collection

The study was conducted in the North block of the Isle of Rum National Nature Reserve in the Inner Hebrides, Scotland (57°N 6°20′W). The island has a mild, wet climate and the vegetation consists of a mosaic of high-quality grassland and low-quality dry and wet heath and blanket bog. The study population comprises ~350 animals at any one time. Neonates are caught during the calving period May–July and individually marked with collars, ear tags, coloured flashes and ear punches, enabling life-long individual identification. Censuses are carried out five times a month for 9 months of the year with more frequent informal monitoring between censuses allowing compilation of individual life histories. The study area population has not been culled since 1973 and runs at the carrying capacity determined by the ground and prevailing weather conditions (Clutton-Brock et al. Reference Clutton-Brock, Guinness and Albon1982). Faecal sampling was conducted on a seasonal basis, with 2-week trips carried out in winter (January), spring (April), summer (August) and autumn (November); each trip was considered to be representative of the 3 month season in which it occurred. Data for this study were all collected in 2016; as red deer are born in May–June this resulted in the study sampling two different cohorts of calves, born 2015 (sampled in winter and spring before their first birthday) and 2016 (sampled in summer and autumn). Groups of individually recognized deer were observed for defecation events from 15 to 250 m using binoculars and telescopes, with the samples recovered as quickly as possible without unduly disturbing the deer. Samples were stored in ziplock bags until processing. Efforts were made to sample as many different individuals on a trip as possible, with a subset of individuals deliberately sampled more than once in each season in order to examine the within-season repeatability of parasitological measures. Faecal analysis can be affected by the hatching, development and death of parasite propagules, influenced by temperature and oxygen availability (Nielsen et al. Reference Nielsen, Vidyashankar, Andersen, DeLisi, Pilegaard and Kaplan2010). For this reason, time of defecation, time of collection and date of the count were all recorded. Following a return to the field station, samples were weighed and homogenized by hand in their ziplock bags to minimize oxygen exposure. These bags were then put inside date-specific larger bags to keep them as anaerobic as possible. All samples were kept refrigerated at 4 °C until parasitological analysis. Upon return to the laboratory, the analysis was carried out within the next 8–10 weeks (see details below). Over the course of four seasons in 2016, Ns = 1020 faecal samples were collected from Ni = 328 individuals, equating to 783 different individual–season combinations with 237 within-season repeats. The sampled individuals were a mixture of calves, yearlings, 2-year olds and adults of both sexes – although adult males were sampled much less frequently (Ns = 43, Ni = 20) than adult females (Ns = 522, Ni = 137) because relatively few adult males live in the study area. The age range of animals sampled was 0–21 years old (median 3 years old).

Parasitology

Parasitological terms will be used as defined in Margolis et al. (Reference Margolis, Esch, Holmes, Kuris and Schad1982), with intensity based on propagule counts. ‘Burden’ refers to the number of worms of a species infecting an individual, which could not be measured directly. Strongyle FECs were carried out within 3 weeks of collection using a modified sedimentation-salt flotation method (Taylor et al. Reference Taylor, Coop and Wall2016) with an accuracy of 1 egg per gram (EPG). The method was modified from Kenyon et al. (Reference Kenyon, McBean, Greer, Burgess, Morrison, Bartley, Bartley, Devin, Nath and Jackson2013). Briefly, 2–15 g of faecal matter was mixed with 10 mL water per gram of faeces and the mixture thoroughly homogenized to suspend the eggs. About 10 mL of this suspension was filtered through a tea strainer and washed through with 5 mL water. The resulting liquid was decanted into a 15 mL polyacrylate test tube, which was centrifuged at 1500 rpm for 2 min and the supernatant removed. The resulting pellet was mixed with the saturated salt solution and resuspended, then centrifuged again, leaving the eggs and light debris at the surface of the liquid. Using medical forceps to clamp below the meniscus, this surface layer was poured off into a cuvette which was topped up with saturated salt and then a lid was added to seal the contents. The entire surface area of the cuvette was counted at 4× magnification to give a count of eggs in 1 g faeces, revealing strongyle nematode eggs and a selection of other species the eggs of which are less dense than the salt solution. This included Nematodirus sp. and Capillaria sp. The assay also revealed oocysts of the coccidian parasite Eimeria sp., which came in two varieties (‘large’ and ‘small’), and segments of the cestode Moniezia expansa. 730 strongyle FECs were repeated to estimate the technical repeatability of this method in our hands. These counts were averaged to give an EPG value for the sample.

Fasciola hepatica eggs were detected through a sedimentation method (French et al. Reference French, Zadoks, Skuce, Mitchell, Gordon-Gibbs, Craine, Shaw, Gibb and Taggart2016; Taylor et al. Reference Taylor, Coop and Wall2016) conducted on 0.5–2 g of faecal matter, from the same homogenate as the strongyle FEC, within ten weeks of collection. In the interim, all samples were kept refrigerated and oxygen-deprived, both of which prevent hatching (Hurtrez-Boussès et al. Reference Hurtrez-Boussès, Meunier, Durand and Renaud2001). After filtering the faecal matter through a tea strainer, the sample was left to sediment in a conical beaker for three minutes. During this time the heavier debris (including the eggs) settles to the bottom of the beaker, and the lighter material can then be removed to leave the eggs with as little debris as possible. The remaining filtrate was pipetted onto a Petri dish and stained with methylene blue (1% w/v). The Petri dish was then examined microscopically at 4× magnification and the eggs, which are yellow against blue debris background, were counted, with the counts divided by the weight in grams of the sample used.

Elaphostrongylus cervi and Dictyocaulus sp. larvae were isolated via a modified Baermannisation assay (Gajadhar et al. Reference Gajadhar, Tessaro and Yates1994) within four weeks of collection. 1–14 g of faecal matter was wrapped in muslin cloth and submerged fully in a 50 mL falcon tube filled with water. This was left at room temperature for 20–24 h for the L1 larvae to emerge from the faeces and fall to the bottom of the tube. The supernatant was then carefully removed to leave <2 mL containing the larvae. This fraction was preserved with Lugol's iodine and kept refrigerated at 4 °C until counting. Counts were performed on a subsample under 40× magnification and divided by the weight of faeces used to give a measure of larvae per gram.

Some samples were not large enough to be analysed for all parasite types – hence final sample sizes were 1014 (strongyles), 991 (F. hepatica) and 1003 (E. cervi). Fluctuations in faecal water content can lead to variation in per-gram FEC, particularly across host sex and age classes and across seasons (Turner et al. Reference Turner, Cizauskas and Getz2010). For this reason, the proportion faecal dry matter (FDM) per gram of collected faeces was calculated for each sample by drying a known weight of faecal matter in an oven at 60 °C for 48 h and then weighing the resulting solid. FDM was mean-centred around 1 to prevent changing the distribution of non-zero counts relative to the zero counts. Counts were divided by FDM to give a measure of eggs or larvae per gram of dry matter, which was rounded to the nearest whole number to allow the use of integer count-based models.

Statistical analysis

Correlations of repeated counts

Since many of the parasite groups identified were present at low prevalence (<30%, see Table 1), we restricted further analyses to the three most prevalent parasite groups: strongyle nematodes, the liver fluke F. hepatica, and the tissue worm Elaphostronylus cervi. Initially, we estimated the repeatability of propagule counts at different temporal scales by calculating Spearman's rank correlation coefficients (r) among: (1) repeated counts from the same sample (‘technical repeatability’); (2) repeated counts from the same individual within a season (‘within-season repeatability’); and (3) averaged within-season measures from the same individual in different seasons (‘between-season repeatability’). We also investigated whether factors related to the collection and processing of the samples influenced egg counts, with simple linear models for each parasite taxon. We investigated the influence of: (a) time of collection, (b) time to processing and (c) time to counting. However, these factors had no effect and so were not included in further analysis.

Table 1. Seasonal prevalence and mean and maximum intensity in propagules/g faecal dry matter of each parasite found.

Minimum count of every season–parasite combination was zero.

Intrinsic and extrinsic factors influencing parasite prevalence and intensity

To test how season, age and sex were associated with parasite egg counts and to decompose the within- and among-individual variation in these counts, we used Generalized Linear Mixed Models (GLMMs). Parasite egg count distributions are often non-normal and strongly overdispersed (Alexander, Reference Alexander2012) and can feature significantly more zero counts than would be expected given their distribution (Chipeta et al. Reference Chipeta, Ngwira, Simoonga and Kazembe2014). The analysis was carried out in R version 3.4.0 (R Core Team, 2017), with the Bayesian statistical package MCMCglmm (Hadfield, Reference Hadfield2010) which is flexible with respect to error structures. In this set of models, we used averages of any repeat measures from the same individual within a season. We fitted zero-inflated Poisson (ZIP) GLMMs, with which we generated two estimates for each fixed or random effect within the model. The first estimate (Poisson) calculated the effect a factor has on the data assuming an overdispersed Poisson distribution, including an expected number of zero counts, while the second (zero-inflation) estimated the effect that a factor has on the number of zeros in the data. We consider these two estimates as reflecting variation in parasite intensity and prevalence, respectively. We follow convention by presenting results for prevalence before intensity for each parasite taxon. This method was chosen as factors can have contrasting effects on parasite prevalence and intensity (Chipeta et al. Reference Chipeta, Ngwira and Kazembe2013). The explanatory variables fitted in the models include age category (with four levels: Calf, Yearling, 2-year-old and adult), sex (female and male) and season (winter, spring, summer and autumn), with individual identity as a random effect to control for variation between individuals (Paterson and Lello, Reference Paterson and Lello2003). Models were run for 2.6 million iterations (thinning interval 2000, burnin 600 000). The significance of differences among factor level means was calculated by comparing the proportion overlap of the posterior distributions of the MCMC estimates for each level and then doubled to give P MCMC following Palmer et al. (Reference Palmer, Hadfield and Obbard2017).

Model-derived repeatability

The variance component associated with the individual random effect within mixed effects models is often used to calculate within-individual repeatabilities (Falconer and Mackay, Reference Falconer and Mackay1996). This method accounts for variation between individuals that occurs as a result of the model's fixed effects (i.e. originating from differences between age, sex and seasonal categories) to estimate the proportion of variation which is explained by differences between individuals. However, there is currently no accepted method of extracting repeatability from the random effects structure of zero-inflated models. We therefore re-ran our models of each parasite count on a non-zero-inflated subset of the data using a standard Poisson model featuring additive overdispersion and applied the method described by Nakagawa and Schielzeth (Reference Nakagawa and Schielzeth2010) to calculate repeatability. We included an individual identity random effect to estimate among-season variation within hosts and an individual-by-season interaction as a second random term to estimate within-season variation within hosts. In the absence of between-season repeatability, the significance of this latter term would demonstrate consistency within individuals of those repeat samples collected within seasons. Repeatabilities were calculated on the count scale rather than the latent scale. The data analysed with this model included all samples including 237 within-season repeat samples. Prior to analysis we removed the winter season data (when counts were very low for all parasite taxa) and a subset of mainly prepatent individuals (calves in the summer for F. hepatica and calves in the summer and autumn for E. cervi), as the repeatability of these counts would be of little biological interest and they were the major source of zero-inflation.

RESULTS

Correlations of repeated counts

The seasonal prevalence, intensity and maximum count of each of the parasite taxa are displayed in Table 1. All parasites were found throughout the year except Dictyocaulus sp. and M. expansa, which were not found in the winter. Repeated strongyle counts of the same sample were strongly correlated (Spearman's R = 0.95, Fig. 1). Within-season resampling of the same individuals revealed strong correlations for all parasites (Spearman's R > 0.6, Fig. 2). Between-season averaged counts from the same individuals correlated less, and varied more between parasites (see Fig. 3 for coefficients). Strongyle nematodes showed the lowest between-season repeatability, with E. cervi and F. hepatica higher. These results were qualitatively similar to the model variance-derived repeatability estimates; see below.

Fig. 1. Correlation between first and second strongyle faecal egg count of the same sample (Ns = 730). R is the Spearman’s rank correlation.

Fig. 2. Within-season correlations of individuals’ counts from different faecal samples (A: Strongylates, B: F. hepatica, C: E. cervi). Axes have been log(x+1) transformed for display purposes. The dashed line represents equal counts, y = x. R is the Spearman’s rank correlation.

Fig. 3. Between-season correlations of parasite counts (A: Strongyles, B: F. hepatica, C: E. cervi). Values are Spearman’s rank correlations; values above the diagonal represent the number of pairs of samples the comparisons are based on.

Intrinsic and seasonal correlates of counts

Model results revealed contrasting trends for all three parasites according to all investigated factors (Table 2; Fig. 4). Strongyle prevalence and intensity peaked in the spring and summer (Fig. 5, P MCMC < 0.001) and decreased in the autumn, with intensity remaining higher than winter levels (P MCMC < 0.001) but not prevalence (P MCMC = 0.114). Spring and summer were not significantly different in prevalence (P MCMC = 0.988) or intensity (P MCMC = 0.086). There was a persistent age effect in that older individuals tended to be less often infected and at lower intensities (Fig. 5). Calves showed a higher prevalence of infection than 2-year-olds (P MCMC = 0.038) and adults (P MCMC < 0.001) and had a higher intensity than all age categories (P MCMC < 0.001). Yearlings also had a higher prevalence than adults (P MCMC < 0.001), and higher intensity than both 2-year-olds and adults (P MCMC = 0.01; P MCMC < 0.001, respectively). Sex also had an effect in strongyles (Fig. 6), with males showing higher intensity infections (P MCMC = 0.024) but no difference in prevalence (P MCMC = 0.436).

Fig. 4. MCMCglmm zero-inflated Poisson model output for each of the three analysed parasite taxa. Points represent posterior estimates for mean effect sizes; error bars represent the 95% credibility intervals of the mean. Symbol corresponds to the statistic being estimated – zero-inflation (prevalence) or Poisson (intensity). Zero-inflation coefficients have been multiplied by −1 to aid interpretation; that is, a positive value represents a decrease in zero-inflation and therefore an increase in prevalence.

Fig. 5. Seasonal mean intensities (±s.e.) for each parasite (A: Strongyles, B: F. hepatica, C: E. cervi) in each age category. Numbers correspond to per cent prevalence. Groups on the x-axis are calves, yearlings, 2-year-olds and adults in order. Figures were created using raw faecal dry matter-transformed data. The calf category represents two different cohorts: those born in 2015 (winter and spring) and those born in 2016 (summer and autumn).

Fig. 6. Mean infection intensity (±s.e.) for each parasite (A: Strongyles, B: F. hepatica, C: E. cervi) according to sex and age category, calculated as in Fig. 5. Numbers above the bars correspond to per cent prevalence. Groups on the x-axis are calves, yearlings, 2-year-olds and adults in order. Figures were created using raw faecal dry matter-transformed data.

Table 2. MCMCglmm model estimates for zero-inflated Poisson GLMMs for each analysed parasite taxon.

The estimates represent the posterior mean, with 95% credibility intervals in brackets. Prevalence (zero-inflation) values have been multiplied by −1 to aid interpretation: a positive value in the model corresponds to an increase in zero-inflation, i.e. a decrease in infection probability. Asterisks represent significance intervals: *, ** and *** refer to P < 0.05, P < 0.01 and P < 0.001, respectively.

Fasciola hepatica increased in intensity after the winter (P MCMC < 0.001) as did strongyles, but decreased in prevalence in the summer (P MCMC < 0.001) and prevalence was highest in the autumn compared with the winter (P MCMC = 0.036), spring (P MCMC = 0.042) and summer (P MCMC < 0.001). Unlike strongyles, F. hepatica prevalence peaked in yearlings (P MCMC = 0.01) rather than in calves, and decreased in prevalence and intensity thereafter. The intensity of infection was lower in 2-year-olds than in calves (P MCMC = 0.024), and adults had lower prevalence and intensity of infection compared to all other age classes (P MCMC = 0.002 for prevalence; P MCMC < 0.001 for intensity). There was no evidence of a sex bias in either prevalence (P MCMC = 0.496) or intensity (P MCMC = 0.386) of F. hepatica infection, in contrast to the results of the strongyle model.

Prevalence of E. cervi was lower in the summer and autumn than in the winter and spring (P MCMC < 0.001, Fig. 5), differing again from the patterns shown by either strongyles or F. hepatica. However, like the other parasites, E. cervi intensity was highest in the spring (P MCMC < 0.001). The E. cervi age trend differed from the other parasites in that older age classes had higher prevalence but the lower intensity of infection. Calves showed a lower E. cervi prevalence than all age classes (P MCMC < 0.001) and lower intensity compared to yearlings (P MCMC < 0.001) and 2-year-olds (P MCMC = 0.022). Adults had a higher prevalence than yearlings (P MCMC = 0.044) but a lower intensity than yearlings and 2-year-olds (Both P MCMC < 0.001). Unlike both strongyles and F. hepatica, E. cervi showed a weak female sex bias, with increased intensity in females (P MCMC = 0.04). Prevalence was also higher in females, although this was not significant (P MCMC = 0.15).

Model-derived repeatability estimates

After accounting for differences in age, sex and season using overdispersed Poisson models, the individual component of variation in count differed between all parasites (P MCMC = 0.014 for F. hepatica and E. cervi; P MCMC < 0.001 otherwise). The direction and significance of the fixed effects from these models closely followed the Poisson component of the ZIP models. Strongyle counts had very low repeatability between-seasons with posterior mode R = 0 (95% credibility intervals 0–0.07). F. hepatica and E. cervi repeatabilities were higher: R = 0.19 (0.12–0.27) and 0.30 (0.25–0.38) respectively. However, strongyle counts were repeatable within-season, as demonstrated by variance accounted for by the ID:Season interaction term, R = 0.38 (0.24–0.49). This term was low for F. hepatica and E. cervi (R = 0 for both), showing that while both were repeatable at all levels this was expressed in the between-season individual identity term, with no additional repeatability within seasons. These results differ markedly from the raw correlations (Figs 2 and 3) as they take into account variation that is attributable to age, sex and season, thereby estimating the within-individual correlations given these factors.

Discussion

We discovered a significant individual component of parasitism despite extensive variation in parasite counts between individuals and between seasons. Both F. hepatica and E. cervi were repeatable between seasons once age and sex category were accounted for using Poisson models, while strongyle counts were only repeatable within-season. Seasonal trends varied, with strongyles showing the most extensive seasonality, though counts of all parasites were lowest in winter. As expected younger individuals tended to have higher intensity infections than adults, but this was not mirrored in increased prevalence apart from in strongyles. Strikingly, sex biases differed between parasites, with higher strongyle intensities in males and marginally higher E. cervi prevalence in females. An overriding feature of this study is the contrasting and asynchronous effects shown by the different parasite taxa: despite all being helminths and having relatively similar life cycles, they exhibited substantially different seasonality and intrinsic trends as well as showing different levels of repeatability which are not attributable to their different detection assays. Studies in wild mammals which investigate multiple pathogens (e.g. Vicente et al. Reference Vicente, Höfle, Fernández-De-Mera and Gortazar2007b) often use distantly related microparasites and macroparasites with very different life cycles. Here we demonstrate the value of investigating multiple high-prevalence parasite taxa even where the chosen parasites are ostensibly similar. In accordance with previous studies (e.g. Chipeta et al. Reference Chipeta, Ngwira and Kazembe2013), zero-inflated Poisson models successfully revealed that some factors affected prevalence and intensity in different directions. While factors that increased prevalence of infections also tended to increase intensity whether their effects were significant or not (Fig. 4), there were notable exceptions in the season categories: the summer (and autumn for E. cervi) featured increased intensity but reduced prevalence of F. hepatica and E. cervi compared with winter as a result of very low prevalence in calves. Similarly, the models distinguished between increased prevalence but reduced intensities of E. cervi in adults compared with younger age classes. We advocate the use of these models where sample sizes are sufficient, particularly in situations where differences between seasons and classes will result in divergent processes affecting helminth prevalence and intensity (e.g. age-related infection prepatency).

Propagule output is a function of both host and parasite biology (Sargison, Reference Sargison2013), representing a combination of adult worm burden and host health as well as fluctuations in worm reproduction; hence faecal propagule counts are subject to fluctuations through time. We have shown that despite this, helminth counts can be repeatable amongst wild individuals within and across seasonal time frames. In the repeated FECs strongyle egg count had a Spearman's rank correlation of R = 0.95, demonstrating high technical repeatability of this counting method. Within-season correlations were high for all parasites, showing that each sample taken was largely representative of an individual's parasite count regardless of when the sample was taken in the day or within the sampling trip; this also demonstrates a high reliability of the assays used and low importance for potential nuisance factors such as time to processing or time of sampling. This is in accordance with helminth studies in other ruminants (e.g. Rinaldi et al. Reference Rinaldi, Veneziano, Morgoglione, Pennacchio, Santaniello, Schioppi, Musella, Fedele and Cringoli2009) – however, there may still be effects of time to collection and analysis that were not detected in this dataset and may reduce repeatability. In addition, it is noteworthy that a number individuals switched between zero and non-zero counts for all parasites (Fig. 2) and therefore a negative propagule count is not necessarily indicative of an uninfected individual, demonstrating the potential value of repeated sampling and the difficulty of diagnosing helminth infection using non-invasive faecal sampling. Similar high repeatability for different parasite taxa is nevertheless surprising, given that propagule counts were performed using three different assays which may differ in their reliability, the low intensity and possible low burden of infection and the fact that propagule shedding of all three taxa is intermittent (Gajadhar et al. Reference Gajadhar, Tessaro and Yates1994; Vercruysse and Claerebout, Reference Vercruysse and Claerebout2001; Schär et al. Reference Schär, Hattendorf, Khieu, Muth, Char, Marti and Odermatt2014).

Between-season correlations were lower as a result of the differences that emerged between individuals between seasons, and these correlations decreased further when repeatability was derived from Poisson models, thereby accounting for variance due to the season, age category and sex. This shows that while our counts were accurate and repeatable, much of the variation was due to certain classes (e.g. calves with strongyles) showing higher prevalence and intensity than others, and when looking within-category individuals’ counts were less consistent. Nevertheless, model-derived R 2 estimates associated with individual identity did not overlap with zero for either F. hepatica or E. cervi, demonstrating consistent differences between individuals throughout the year. Strongyle counts were not as repeatable between seasons but did show individual consistency within seasons (R 2 = 0.38), demonstrating that individual repeatability decreased between seasons rather than being absent at all levels. This R 2 value is intermediate compared to other studies investigating strongyle FEC of wild horses within and between seasons (Wood et al. Reference Wood, Matthews, Stephenson, Slote and Nussey2013; Debeffe et al. Reference Debeffe, McLoughlin, Medill, Stewart, Andres, Shury, Wagner, Jenkins, Gilleard and Poissant2016) but low compared to between-season repeatability in farmed horses (Scheuerle et al. Reference Scheuerle, Stear, Honeder, Becher and Pfister2016) and goats (Hoste et al. Reference Hoste, Le Frileux, Goudeau, Chartier, Pors, Broqua and Bergeaud2002). Low repeatability compared with farmed animals is unsurprising given the large range of different conditions that wild individuals experience. High variability of parasite counts within individuals has implications for their quantification in wildlife. For example, an individual that shows high strongyle FEC in the spring may not do so in the summer. This demonstrates the value of multiple sampling seasons as well as reflecting the evolutionary ecology of the deer and their helminths. For example, different seasonal peaks may be associated with between-individual variation in seasonal trade-offs with immunity through e.g. reproduction (Martin et al. Reference Martin, Weil and Nelson2008) or with varying levels of tolerance to infection. Similarly, asynchronous peaks of egg output across the host population may be adaptive for the helminths in encouraging year-round transmission and bet-hedging to buffer for unfavourable climatic conditions.

All parasites showed some transmission in each sampling trip, although as the expected intensity of infection was lowest in the winter for all three parasites, with peaks in the spring and summer which continued into the autumn for F. hepatica and E. cervi. This low transmission in colder seasons likely reflects a reduction in egg production rather than solely burden, as strongyle burden in Spanish red deer stays constant or increases in the winter (Santín-Durán et al. Reference Santín-Durán, Alunda, Hoberg and de la Fuente2008) and F. hepatica is found regularly in necropsies of the Rum deer throughout the mortality period (personal observation). Freezing, which regularly occurs on the ground in the winter on Rum, is known to damage strongyle eggs and larvae (Foreyt, Reference Foreyt1986; Wharton and Allan, Reference Wharton and Allan1989) and F. hepatica eggs (French et al. Reference French, Zadoks, Skuce, Mitchell, Gordon-Gibbs, Craine, Shaw, Gibb and Taggart2016). However, overwinter transmission can occur in some parasite species (e.g. Carlsson et al. Reference Carlsson, Irvine, Wilson, Piertney, Halvorsen, Coulson, Stien and Albon2012); future work identifying the strongyle species present may be able to identify whether winter and autumn strongyle transmission involves a few frost-resistant species. As well as increasing the survival of environmental stages, high helminth transmission in the spring is a possible adaptive strategy resulting in coincidence between vulnerable calves and maximum infective parasites in the environment, similarly to the periparturient rise in sheep strongyle FEC (Armour, Reference Armour1980). Fasciola hepatica and E. cervi showed less variation between seasons than did strongyles, which may be linked to their reduced reliance on the influx of young naïve individuals, as well as reflecting a longer lifespan of individual parasites and/or hardier propagules more capable of year-round transmission. Low seasonality of these parasites is surprising given that both go through intermediate snail hosts, their infection of which would be expected to rely on weather conditions (Olsen et al. Reference Olsen, Frankena, Bødker, Toft, Thamsborg, Enemark and Halasa2015; Kim et al. Reference Kim, Echaubard, Suwannatrai, Kaewkes, Wilcox and Sripa2016), and temperatures above 10 °C are required for F. hepatica egg development (Ollerenshaw and Smith, Reference Ollerenshaw and Smith1969). Elaphostrongylus cervi output varies according to monthly rainfall patterns in Spanish red deer (Vicente et al. Reference Vicente, Fierro and Gortazar2005), so larval output may be dependent on environmental cues on a shorter timescale rather than fluctuating annually as do strongyles. It is important to clarify that verifying seasonal dynamics such as this would require continuous sampling rather than employing discrete seasons as we do here for practical reasons. We therefore cannot confirm exactly during which period each parasite peaks at an individual- or population-level; however, the high repeatability of F. hepatica and E. cervi between seasons supports higher population-level synchrony, with more variable strongyle transmission. Asynchrony in seasonal peaks between parasites will affect host-parasite interactions by necessitating different immune responses at different times of the year. In this case, immunity to strongyles is likely to rise in the warmer months; lower seasonality of F. hepatica and E. cervi transmission will necessitate year-round immunity to these parasites, possibly interacting with seasonal costs of variation in nutrition, mating, reproduction and maternal care.

Young individuals tended to experience higher intensity infections than adults, likely playing an important role in maintaining and transmitting helminth infections in the population. Strongyles were particularly age-biased, with calves showing a higher prevalence and intensity that decreased with each successive life stage in a similar pattern to that seen in other studies of ungulate strongyles. F. hepatica and E. cervi also closely followed previously-seen age profiles, increasing in prevalence at the yearling stage (Vicente et al. Reference Vicente, Pérez-Rodríguez and Gortazar2007a; French et al. Reference French, Zadoks, Skuce, Mitchell, Gordon-Gibbs, Craine, Shaw, Gibb and Taggart2016). However, F. hepatica prevalence and intensity decreased in adults while E. cervi increased in prevalence but decreased in intensity past the yearling stage. This age-biased infection implies that studies based on selective culling regimes that focus on, for example, adults, may indeed be missing relevant season-group categories which are important in determining the extent of parasitism within a population. For example, the high numbers of calves and yearlings sampled in our study contributed to the high prevalence of F. hepatica seen here (>50%) compared to that in a study of Scottish deer which largely used culled adults (French et al. Reference French, Zadoks, Skuce, Mitchell, Gordon-Gibbs, Craine, Shaw, Gibb and Taggart2016, mean 26% prevalence). The stronger age bias in strongyles may result from strongyle infection causing more mortality or more effective adaptive immunity than F. hepatica and especially E. cervi. This concurs with the view of some strongyle species as highly pathogenic (Hoberg et al. Reference Hoberg, Kocan, Rickard, Samuel, Pybus and Kocan2001) while E. cervi is often asymptomatic in red deer (Alberti et al. Reference Alberti, Gioia, Sironi, Zanzani, Riccaboni, Magrini and Manfredi2011). Despite having high-intensity infections in spring, and in contrast to strongyles, calves exhibited a very low prevalence of F. hepatica and E. cervi in the summer and autumn, as revealed by the zero-inflation of the models (Fig. 4). This was a likely result of infection prepatency: F. hepatica has an 8-week prepatent period in cattle (de León et al. Reference de León, Quiñones and Hillyer1981), while E. cervi can take 80–200 days to develop depending on dose (Gajadhar et al. Reference Gajadhar, Tessaro and Yates1994). This influence of prepatency on prevalence patterns reinforces the need to understand seasonality when investigating helminths using non-invasive methods: it is likely that many calves with zero or low counts in the summer were in fact heavily infected with both parasites, but this was not yet detectable using faecal examination. Thus, quantifying the repeatability of infection with these parasites in individual calves was not possible until the autumn (for F. hepatica) or the spring (for E. cervi).

We expected to see a male bias in all three parasites (Moore and Wilson, Reference Moore and Wilson2002; Zuk, Reference Zuk2009); however, all three parasites differed in their distribution between the sexes, with a male bias in strongyle intensity and a female bias in E. cervi intensity, both evident in the first year, and with no effect evident in F. hepatica. The male bias in strongyle infection arose late in life compared to that in Soay sheep, in which male counts greatly increase relative to female counts within months of birth (Wilson et al. Reference Wilson, Grenfell, Pilkington, Boyd, Gulland, Clutton-Brock and Pemberton2004). The difference is also small (Fig. 4), and F. hepatica showed no sex difference, despite a previous (seasonally confounded) study showing a higher prevalence in males than females (French et al. Reference French, Zadoks, Skuce, Mitchell, Gordon-Gibbs, Craine, Shaw, Gibb and Taggart2016). On Rum most adult males live outside the study area, in areas which are at a lower density as a result of culling (Clutton-Brock et al. Reference Clutton-Brock, Coulson, Milner-Gulland, Thomson and Armstrong2002) and with different geography and grazing. This may reduce exposure levels while influencing susceptibility through differences in diet. The weak female bias in E. cervi infection is the opposite of the expected pattern, and is particularly surprising as other red deer studies have shown male biases in E. cervi (Vicente et al. Reference Vicente, Fernández De Mera and Gortazar2006, Reference Vicente, Pérez-Rodríguez and Gortazar2007a); opposing sex effects between parasites have been reported, e.g. in ectoparasites and endoparasites of African ground squirrels (Hillegass et al. Reference Hillegass, Waterman and Roth2008) and in stickleback helminths (Reimchen and Nosil, Reference Reimchen and Nosil2001), but due to a dearth of studies investigating multiple similar parasites we are not aware of any in ungulates. A possible explanation for this female bias is lower maternal care: it has been shown that hinds invest slightly less in female than male calves (Froy et al. Reference Froy, Walling, Pemberton, Clutton-brock and Kruuk2016), and earlier weaning may result in female calves suffering higher exposure to E. cervi, leading to earlier patent infections, and/or they may be weaker and therefore more susceptible so show higher-intensity infections. This disparity in sex bias disagrees with the expected male bias in helminth infection and implies that male deer may in fact feature a difference in the community rather than the extent of parasitism.

Acknowledgements

We thank Scottish Natural Heritage for permission to work on the Isle of Rum NNR and for the support of the reserve management team on the island. Thanks to Dave McBean and Gillian Mitchell at the Moredun Research Institute for their help with parasitological methods. Thanks also to Jarrod Hadfield for statistical advice and to Olly Gibb, William Mason, Jonty Lord and Alice Thomson for help in the field.

Financial Support

The long-term red deer study is funded by the Natural Environment Research Council (grant number NE/L00688X/1), as is GFA's PhD studentship through the E3 Doctoral Training Partnership (grant number NE/L002558/1).

Author contributions

G.F.A. collected samples, performed counts, analysed the data and drafted the manuscript; A.M. and S.M. helped collect samples; D.N., J.P. and F.K. offered statistical advice and helped draft the manuscript.

References

Alberti, EG, Gioia, G, Sironi, G, Zanzani, S, Riccaboni, P, Magrini, M and Manfredi, MT (2011) Elaphostrongylus cervi in a population of red deer (Cervus elaphus) and evidence of cerebrospinal nematodiasis in small ruminants in the province of Varese, Italy. Journal of Helminthology 85, 313318.Google Scholar
Alexander, N (2012) Analysis of parasite and other skewed counts. Tropical Medicine & International Health 17, 684693.Google Scholar
Altizer, S, Dobson, A, Hosseini, P, Hudson, P, Pascual, M and Rohani, P (2006) Seasonality and the dynamics of infectious diseases. Ecology Letters 9, 467484.Google Scholar
Armour, J (1980) The epidemiology of helminth disease in farm animals. Veterinary Parasitology 6, 746.Google Scholar
Bohm, M, White, PCL, Chambers, J, Smith, L and Hutchings, MR (2007) Wild deer as a source of infection for livestock and humans in the UK. Veterinary Journal 174, 260276.Google Scholar
Budischak, SA, Jolles, AE, Ezenwa, VO, Hoberg, EP, Abrams, A, Jolles, AE and Ezenwa, VO (2015) A combined parasitological-molecular approach for non-invasive characterization of parasitic nematode communities in wild hosts. Molecular Ecology Resources 15, 11121119.Google Scholar
Carlsson, AM, Irvine, RJ, Wilson, K, Piertney, SB, Halvorsen, O, Coulson, SJ, Stien, A and Albon, SD (2012) Disease transmission in an extreme environment: nematode parasites infect reindeer during the Arctic winter. International Journal for Parasitology 42, 789795.Google Scholar
Chintoan-Uta, C, Morgan, ER, Skuce, PJ and Coles, GC (2014) Wild deer as potential vectors of anthelmintic-resistant abomasal nematodes between cattle and sheep farms. Proceedings of the Royal Society B, Biological Sciences 281, 20132985.Google Scholar
Chipeta, MG, Ngwira, B and Kazembe, LN (2013) Analysis of Schistosomiasis haematobium infection prevalence and intensity in Chikhwawa, Malawi: an application of a two part model. PLoS Neglected Tropical Diseases 7, e2131.Google Scholar
Chipeta, MG, Ngwira, BM, Simoonga, C and Kazembe, LN (2014) Zero adjusted models with applications to analysing helminths count data. BMC Research Notes 7, 856.Google Scholar
Clutton-Brock, T and Sheldon, BC (2010) Individuals and populations: the role of long-term, individual-based studies of animals in ecology and evolutionary biology. Trends in Ecology and Evolution 25, 562573.Google Scholar
Clutton-Brock, TH, Guinness, FE and Albon, SD (1982) Red Deer: Behavior and Ecology of Two Sexes. Chicago, IL: University of Chicago Press.Google Scholar
Clutton-Brock, TH, Coulson, TN, Milner-Gulland, EJ, Thomson, D and Armstrong, HM (2002) Sex differences in emigration and mortality affect optimal management of deer populations. Nature 415, 633637.Google Scholar
Connan, RM (1997) Hypobiosis in the ostertagids of red deer and the efficacy of ivermectin and fenbendazole against them. Veterinary Record 140, 203205.Google Scholar
Craig, BH, Pilkington, JG and Pemberton, JM (2006) Gastrointestinal nematode species burdens and host mortality in a feral sheep population. Parasitology 133, 485496.Google Scholar
Davidson, RK, Kutz, SJ, Madslien, K, Hoberg, E and Handeland, K (2014) Gastrointestinal parasites in an isolated Norwegian population of wild red deer (Cervus elaphus). Acta Veterinaria Scandinavica 56, 18.Google Scholar
Debeffe, L, McLoughlin, PD, Medill, SA, Stewart, K, Andres, D, Shury, T, Wagner, B, Jenkins, E, Gilleard, JS and Poissant, J (2016) Negative covariance between parasite load and body condition in a population of feral horses. Parasitology 143, 983997.Google Scholar
de León, D, Quiñones, R and Hillyer, GV (1981) The prepatent and patent periods of Fasciola hepatica in cattle in Puerto Rico. Journal of Parasitology 67, 734735.Google Scholar
Falconer, DS and Mackay, TFC (1996) Introduction to Quantitative Genetics. 4th edn. Harlow, Essex, UK: Longman.Google Scholar
Foreyt, WJ (1986) Recovery of nematode eggs and larvae in deer: evaluation of fecal preservation methods. Journal of the American Veterinary Medical Association 189, 10651067.Google Scholar
French, AS, Zadoks, RN, Skuce, PJ, Mitchell, G, Gordon-Gibbs, DK, Craine, A, Shaw, D, Gibb, SW and Taggart, MA (2016) Prevalence of liver fluke (Fasciola hepatica) in wild red deer (Cervus elaphus): coproantigen ELISA is a practicable alternative to faecal egg counting for surveillance in remote populations. PLoS ONE 11, 118.Google Scholar
Froy, H, Walling, CA, Pemberton, JM, Clutton-brock, TH and Kruuk, LEB (2016) Relative costs of offspring sex and offspring survival in a polygynous mammal. Biology Letters 12, 20160417.Google Scholar
Gajadhar, AA, Tessaro, SV and Yates, WDG (1994) Diagnosis of Elaphostrongylus cervi infection in New Zealand red deer (Cervus elaphus) quarantined in Canada, and experimental determination of a new extended prepatent period. Canadian Veterinary Journal 35, 433437.Google Scholar
Hadfield, JD (2010) MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. Journal of Statistical Software 33, 122.Google Scholar
Hayward, AD, Wilson, AJ, Pilkington, JG, Pemberton, JM and Kruuk, LEB (2009) Ageing in a variable habitat: environmental stress affects senescence in parasite resistance in St Kilda Soay sheep. Proceedings of the Royal Society B, Biological Sciences 276, 34773485.Google Scholar
Hillegass, MA, Waterman, JM and Roth, JD (2008) The influence of sex and sociality on parasite loads in an African ground squirrel. Behavioral Ecology 19, 10061011.Google Scholar
Hines, AM, Ezenwa, VO, Cross, P and Rogerson, JD (2007) Effects of supplemental feeding on gastrointestinal parasite infection in elk (Cervus elaphus): preliminary observations. Veterinary Parasitology 148, 350355.Google Scholar
Hoberg, EP, Kocan, AA and Rickard, LG (2001) Gastrointestinal strongyles in wild ruminants. In Samuel, MD, Pybus, MJ and Kocan, AA (eds) Parasitic Diseases of Wild Mammals. Iowa City, IA: Iowa State University Press, pp. 193227.Google Scholar
Hoste, H, Le Frileux, Y, Goudeau, C, Chartier, C, Pors, I, Broqua, C and Bergeaud, JP (2002) Distribution and repeatability of nematode faecal egg counts in dairy goats: a farm survey and implications for worm control. Research in Veterinary Science 72, 211215.Google Scholar
Hurtrez-Boussès, S, Meunier, C, Durand, P and Renaud, F (2001) Dynamics of host – parasite interactions: the example of population biology of the liver fluke (Fasciola hepatica). Microbes and Infection 3, 841849.Google Scholar
Irvine, RJ, Corbishley, H, Pilkington, JG and Albon, SD (2006) Low-level parasitic worm burdens may reduce body condition in free-ranging red deer (Cervus elaphus). Parasitology 133, 465475.Google Scholar
Jolles, AE and Ezenwa, VO (2015) Ungulates as model systems for the study of disease processes in natural populations. Journal of Mammalogy 96, 415.Google Scholar
Kenyon, F, McBean, D, Greer, AW, Burgess, CGS, Morrison, AA, Bartley, DJ, Bartley, Y, Devin, L, Nath, M and Jackson, F (2013) A comparative study of the effects of four treatment regimes on ivermectin efficacy, body weight and pasture contamination in lambs naturally infected with gastrointestinal nematodes in Scotland. International Journal for Parasitology: Drugs and Drug Resistance 3, 7784.Google Scholar
Kim, CS, Echaubard, P, Suwannatrai, A, Kaewkes, S, Wilcox, BA and Sripa, B (2016) Seasonal and spatial environmental influence on Opisthorchis viverrini intermediate hosts, abundance, and distribution: insights on transmission dynamics and sustainable control. PLoS Neglected Tropical Diseases 10, e0005121.Google Scholar
Maizels, RM, Hewitson, JP and Smith, KA (2012) Susceptibility and immunity to helminth parasites. Current Opinion in Immunology 24, 459466.Google Scholar
Malo, AF, Roldan, ERS, Garde, JJ, Soler, AJ, Vicente, J, Gortazar, C and Gomendio, M (2009) What does testosterone do for red deer males? Proceedings of the Royal Society B, Biological Sciences 276, 971980.Google Scholar
Margolis, L, Esch, GW, Holmes, JC, Kuris, AM and Schad, GA (1982) The use of ecological terms in parasitology (report of an ad hoc committee of the American society of parasitologists). Journal of Parasitology 68, 131133.Google Scholar
Martin, LB, Weil, ZM and Nelson, RJ (2008) Seasonal changes in vertebrate immune activity: mediation by physiological trade-offs. Philosophical Transactions of the Royal Society B: Biological Sciences 363, 321339.Google Scholar
Mason, P (1994) Parasites of deer in New Zealand. New Zealand Journal of Zoology 21, 3947.Google Scholar
McKenna, PB (1981) The diagnosis value and interpretation of faecal egg counts in sheep. New Zealand Veterinary Journal 29, 129132.Google Scholar
Moore, SL and Wilson, K (2002) Parasites as a viability cost of sexual selection in natural populations of mammals. Science 297, 20152018.Google Scholar
Nakagawa, S and Schielzeth, H (2010) Repeatability for Gaussian and non-Gaussian data: a practical guide for biologists. Biological Reviews 85, 935956.Google Scholar
Nielsen, MK, Vidyashankar, AN, Andersen, UV, DeLisi, K, Pilegaard, K and Kaplan, RM (2010) Effects of fecal collection and storage factors on strongylid egg counts in horses. Veterinary Parasitology 167, 5561.Google Scholar
Ollerenshaw, CB and Smith, LP (1969) Meteorological factors and forecasts of helminthic disease. Advances in Parasitology 7, 283323.Google Scholar
Olsen, A, Frankena, K, Bødker, R, Toft, N, Thamsborg, SM, Enemark, HL and Halasa, T (2015) Prevalence, risk factors and spatial analysis of liver fluke infections in Danish cattle herds. Parasites & Vectors 8, 160.Google Scholar
Palmer, WH, Hadfield, J and Obbard, DJ (2017) RNA interference pathways display high rates of adaptive protein evolution across multiple invertebrates. Genetics. Available at https://doi.org/10.1534/genetics.117.300567Google Scholar
Paterson, S and Lello, J (2003) Mixed models: getting the best use of parasitological data. Trends in Parasitology 19, 370375.Google Scholar
Poulin, R (1996) Helminth growth in vertebrate hosts: does host sex matter? International Journal for Parasitology 26, 13111315.Google Scholar
Poulin, R (2007) Evolutionary Ecology of Parasites. Princeton, NJ: Princeton University Press.Google Scholar
R Core Team (2017) R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.Google Scholar
Reimchen, TE and Nosil, P (2001) Ecological causes of sex-biased parasitism in threespine stickleback. Biological Journal of the Linnean Society 73, 5163.Google Scholar
Rinaldi, L, Veneziano, V, Morgoglione, ME, Pennacchio, S, Santaniello, M, Schioppi, M, Musella, V, Fedele, V and Cringoli, G (2009) Is gastrointestinal strongyle faecal egg count influenced by hour of sample collection and worm burden in goats? Veterinary Parasitology 163, 8186.Google Scholar
Santín-Durán, M, Alunda, JM, Hoberg, EP and de la Fuente, C (2008) Age distribution and seasonal dynamics of abomasal helminths in wild red deer from central Spain. Journal of Parasitology 94, 10311037.Google Scholar
Sargison, ND (2013) Understanding the epidemiology of gastrointestinal parasitic infections in sheep: what does a faecal helminth egg count tell us? Small Ruminant Research 110, 7881.Google Scholar
Schär, F, Hattendorf, J, Khieu, V, Muth, S, Char, MC, Marti, HP and Odermatt, P (2014) Strongyloides stercoralis larvae excretion patterns before and after treatment. Parasitology 141, 892897.Google Scholar
Scheuerle, MC, Stear, MJ, Honeder, A, Becher, AM and Pfister, K (2016) Repeatability of strongyle egg counts in naturally infected horses. Veterinary Parasitology 228, 103107.Google Scholar
Stear, MJ and Murray, M (1994) Genetic resistance to parasitic disease: particularly of resistance in ruminants to gastrointestinal nematodes. Veterinary Parasitology 54, 161176.Google Scholar
Stromberg, BE (1997) Environmental factors influencing transmission. Veterinary Parasitology 72, 247256–264.Google Scholar
Taylor, MA, Coop, RL and Wall, RL (2016) Laboratory diagnosis of parasitism. In Veterinary Parasitology. Hoboken, NJ: John Wiley & Sons, pp. 259312.Google Scholar
Turner, WC and Getz, WM (2010) Seasonal and demographic factors influencing gastrointestinal parasitism in ungulates of Etosha National Park. Journal of Wildlife Diseases 46, 11081119.Google Scholar
Turner, WC, Cizauskas, CA and Getz, WM (2010) Variation in faecal water content may confound estimates of gastro-intestinal parasite intensity in wild African herbivores. Journal of Helminthology 84, 99105.Google Scholar
van de Pol, M and Verhulst, S (2006) Age-dependent traits: a new statistical model to separate within- and between-individual effects. The American Naturalist 167, 766773.Google Scholar
Vercruysse, J and Claerebout, E (2001) Treatment vs non-treatment of helminth infections in cattle: defining the threshold. Veterinary Parasitology 98, 195214.Google Scholar
Vicente, JJ, Fierro, Y and Gortazar, C (2005) Seasonal dynamics of the fecal excretion of Elaphostrongylus cervi (Nematoda, Metastrongyloidea) first-stage larvae in Iberian red deer (Cervus elaphus hispanicus) from southern Spain. Parasitology Research 95, 6064.Google Scholar
Vicente, J, Fernández De Mera, IG and Gortazar, C (2006) Epidemiology and risk factors analysis of elaphostrongylosis in red deer (Cervus elaphus) from Spain. Parasitology Research 98, 7785.Google Scholar
Vicente, J, Pérez-Rodríguez, L and Gortazar, C (2007a) Sex, age, spleen size, and kidney fat of red deer relative to infection intensities of the lungworm Elaphostrongylus cervi. Naturwissenschaften 94, 581587.Google Scholar
Vicente, J, Höfle, U, Fernández-De-Mera, IG and Gortazar, C (2007b) The importance of parasite life history and host density in predicting the impact of infections in red deer. Oecologia 152, 655664.Google Scholar
Wharton, DA and Allan, GS (1989) Cold tolerance mechanisms of the free-living stages of Trichostrongylus colubriformis (Nematoda). Journal of Experimental Biology 145, 353369.Google Scholar
Wilson, K, Grenfell, BT, Pilkington, JG, Boyd, HEG and Gulland, FMD (2004). Parasites and their impact. In Clutton-Brock, T and Pemberton, J (eds). Soay Sheep: Dynamics and Selection in an Island Population. Cambridge, UK: Cambridge University Press, pp. 113165.Google Scholar
Wood, ELD, Matthews, JB, Stephenson, S, Slote, M and Nussey, DH (2013) Variation in fecal egg counts in horses managed for conservation purposes: individual egg shedding consistency, age effects and seasonal variation. Parasitology 140, 115128.Google Scholar
Zuk, M (2009) The sicker sex. PLoS Pathogens 5, 13.Google Scholar
Figure 0

Table 1. Seasonal prevalence and mean and maximum intensity in propagules/g faecal dry matter of each parasite found.

Figure 1

Fig. 1. Correlation between first and second strongyle faecal egg count of the same sample (Ns = 730). R is the Spearman’s rank correlation.

Figure 2

Fig. 2. Within-season correlations of individuals’ counts from different faecal samples (A: Strongylates, B: F. hepatica, C: E. cervi). Axes have been log(x+1) transformed for display purposes. The dashed line represents equal counts, y = x. R is the Spearman’s rank correlation.

Figure 3

Fig. 3. Between-season correlations of parasite counts (A: Strongyles, B: F. hepatica, C: E. cervi). Values are Spearman’s rank correlations; values above the diagonal represent the number of pairs of samples the comparisons are based on.

Figure 4

Fig. 4. MCMCglmm zero-inflated Poisson model output for each of the three analysed parasite taxa. Points represent posterior estimates for mean effect sizes; error bars represent the 95% credibility intervals of the mean. Symbol corresponds to the statistic being estimated – zero-inflation (prevalence) or Poisson (intensity). Zero-inflation coefficients have been multiplied by −1 to aid interpretation; that is, a positive value represents a decrease in zero-inflation and therefore an increase in prevalence.

Figure 5

Fig. 5. Seasonal mean intensities (±s.e.) for each parasite (A: Strongyles, B: F. hepatica, C: E. cervi) in each age category. Numbers correspond to per cent prevalence. Groups on the x-axis are calves, yearlings, 2-year-olds and adults in order. Figures were created using raw faecal dry matter-transformed data. The calf category represents two different cohorts: those born in 2015 (winter and spring) and those born in 2016 (summer and autumn).

Figure 6

Fig. 6. Mean infection intensity (±s.e.) for each parasite (A: Strongyles, B: F. hepatica, C: E. cervi) according to sex and age category, calculated as in Fig. 5. Numbers above the bars correspond to per cent prevalence. Groups on the x-axis are calves, yearlings, 2-year-olds and adults in order. Figures were created using raw faecal dry matter-transformed data.

Figure 7

Table 2. MCMCglmm model estimates for zero-inflated Poisson GLMMs for each analysed parasite taxon.