Free

# Natural Selection on Antihelminth Antibodies in a Wild Mammal Population

## Abstract

An effective immune response is expected to confer fitness benefits through improved resistance to parasites but also incur energetic costs that negatively impact fitness-related traits, such as reproduction. The fitness costs and benefits of an immune response are likely to depend on host age, sex, and levels of parasite exposure. Few studies have examined the full extent to which patterns of natural selection on immune phenotypes vary across demographic groups and environments in the wild. Here, we assessed natural selection on plasma levels of three functionally distinct isotypes (IgA, IgE, and IgG) of antibodies against a prevalent nematode parasite measured in a wild Soay sheep population over 26 years. We found little support for environment-dependent selection or reproductive costs. However, antibody levels were negatively associated with parasite egg counts and positively associated with subsequent survival, albeit in a highly age- and isotype-dependent manner. Raised levels of antiparasite IgA best predicted reduced egg counts, but this did not predict survival in lambs. In adults increased antiparasite IgG predicted reduced egg counts, and in adult females IgG levels also positively predicted overwinter survival. Our results highlight the potential importance of age- and sex-dependent selection on immune phenotypes in nature and show that patterns of selection can vary even among functionally related immune markers.

## Introduction

Parasites have a major impact on host condition and fitness and thereby represent a strong selective force for individuals in natural populations (Schmid-Hempel 2011). The main defense against parasites is the host immune system (Murphy 2012), and consequently it is expected for there to be a strong directional selection in favor of robust and effective immune responses in natural populations (Nunn et al. 2000; Lindström et al. 2004; Scharsack et al. 2007). Directional selection on a phenotypic trait is predicted to erode genetic variation in that trait (Falconer and Mackay 1996), but considerable genetic variation underlying parasite burden, infection risk, and immune phenotypes has been observed in humans and laboratory studies and in the wild (Lazzaro and Little 2009; Maizels and Nussey 2013). A key explanation for the maintenance of genetic variation in immune responses in the face of selection lies in the fact that mounting an immune response can incur considerable energetic costs (Lochmiller and Deerenberg 2000). Although investing in immunity may promote survival in the face of parasite infection, it draws resources away from other key fitness functions, such as growth and reproduction (Sheldon and Verhulst 1996). Alternatively, mounting strong immune responses can lead to damage to host tissue and autoimmunity, which would impose a fitness cost (Graham et al. 2005; Viney et al. 2005). These processes could result in a pattern of stabilizing—rather than directional selection—on immune phenotypes (Seppälä 2015).

Both the costs of investing in immune responses and the pattern of selection on immune traits are expected to vary with environment, sex, and age (Seppälä 2015). Parasite exposure risk is expected to vary over space and time, and the cost-to-benefit balance of investing in immunity to that parasite will therefore also vary (Altizer et al. 2006). Similarly, heightened intrasexual competition and immunosuppressive effects of testosterone experienced by males (particularly in polygynous species) leads to the expectation of differences in the fitness cost and benefits of investing in immunity between the sexes (Zuk 1990). Both immune phenotype and the ability to mount an effective response to infection also vary profoundly with host age (Palacios et al. 2011; Simon et al. 2015). In particular, young immunologically naïve hosts tend to have the highest infection rates and parasite burdens, and the fitness costs and benefits of investing resources in immunity during early life may differ profoundly from those in adulthood (Jackson et al. 2014). Given the complexity of the vertebrate immune system, it is also important to consider that the pattern of natural selection may also vary considerably across different components of the immune response and depend critically on the immune assays chosen in a particular study (e.g., Gonzalez et al. 1999; Råberg and Stjernman 2003; Parejo and Silva 2009; Nussey et al. 2014; Watson et al. 2016). For instance, the magnitude of response to a nonspecific challenge—such as the swelling response to phytohemagglutinin challenge widely used in ecoimmunology—may be informative about selection on the host’s general capacity to respond to immune challenge but may not necessarily reflect how selection operates on specific immune responses to ecologically relevant infectious agents (Demas et al. 2011; Graham et al. 2011).

In natural populations, variation in resource acquisition will influence the total availability of resources for investment in immunity and other life-history traits (Siva-Jothy and Thompson 2002; Seppälä and Jokela 2010). This may generate positive associations between immune traits and fitness-related traits because of differences in host condition and resource acquisition and obscure trade-offs between immunity and life-history traits (van Noordwijk and de Jong 1986; Kraaijeveld and Godfray 1997; Reznick et al. 2000). Despite this, studies in natural vertebrate populations do provide support for the hypothesis that increases in measures associated with immune responsiveness predict improved survival (Saino et al. 1997; Christe et al. 1998, 2001; Merino et al. 2000; Nussey et al. 2014; Watson et al. 2016) but reduced reproductive performance (Ilmonen et al. 2000; Råberg et al. 2000; Bonneaud et al. 2003; Uller et al. 2006; Marzal et al. 2007; Gasparini et al. 2009; Graham et al. 2010; Hayward et al. 2014). There is also emerging evidence that selection on immune measurements may be environment dependent (Svensson et al. 2001; Calsbeek et al. 2008; Graham et al. 2010) and that patterns of selection may vary among immune measures (Gonzalez et al. 1999; Parejo and Silva 2009; Nussey et al. 2014; Watson et al. 2016). For example, antibody responses of wild blue tits to diphtheria and tetanus antigens showed positive linear and curvilinear relationships, respectively, with overwinter survival (Råberg and Stjernman 2003). To date, few studies have comprehensively tested for sex-, environment-, and age-dependent associations between immune measurements and subsequent survival and fecundity, while correcting for host body weight, in wild vertebrate populations.

Here, we examine the relationships among three functionally distinct isotypes (IgA, IgE, and IgG) of antinematode antibodies and nematode parasite burden, host body weight, survival, and fecundity in the free-living St. Kilda Soay sheep population. The Soay sheep are infected with a number of parasitic helminths, including a variety of gastrointestinal strongyle nematodes, largely comprised of the species Teladorsagia circumcincta, Trichostrongylus axei, and Trichostrongylus vitrinus (Wilson et al. 2004; Craig et al. 2006). Strongyle nematode burden, in combination with harsh winter weather and low food availability, are a strong selective force on the sheep (Gulland 1992; Wilson et al. 2004; Craig et al. 2006; Hayward et al. 2011). An exhaustive examination of both directional and stabilizing selection on strongyle fecal egg count (FEC) found evidence of age and environment dependencies: high FEC predicted reduced survival but only in lambs when population size was low (Hayward et al. 2011). Immunological research on domestic sheep and laboratory mice demonstrates that antibodies—particularly those of the functionally distinct isotypes IgA, IgE, and IgG—play an important role in both the development and maintenance of immunity to these parasites (Stear et al. 2009; Harris and Gause 2011). IgG is the most common antibody isotype in blood and extracellular fluid and has the longest half-life of all antibody isotypes. IgG is effective at activating the complement system, neutralizing toxins, and acting as an opsonin (Murphy 2012). In mice, IgG is crucially linked to protective immunity to helminth infections (Appleton and McGregor 1987; Blackwell and Else 2001; Harris et al. 2006; McCoy et al. 2008), which may be mediated by complement (Esser-von Bieren et al. 2013). IgA is the predominant antibody isotype at mucosal secretions and is associated with neutralizing pathogens in the gut. IgE has the shortest half-life of all antibody isotypes and is typically found bound to mast cells below the skin and mucosa. IgE is the antibody isotype most frequently linked to defense against helminths and mediates allergic inflammation via mast cell degranulation and the release of cytotoxic molecules from basophils and eosinophils (Murphy 2012). In domestic lambs, reduced FEC are associated with increased parasite-specific IgA levels directed at worm growth and subsequent female fecundity (Stear et al. 1995; Strain et al. 2002) as well as the arrestment of development of L4 larvae (Stear et al. 2004). In addition, IgE-dependent mast cell degranulation is associated with reducing adult worm burden (Stear et al. 1995; Huntley et al. 2001). In previous cross-sectional studies using samples collected in 1 year or a few years in wild Soay sheep, we have assayed antibodies against antigens from larval stages of T. circumcincta. We detected highly isotype-specific relationships with body weight and FEC and positive associations between anti–T. circumcincta IgG antibody levels and subsequent survival (Nussey et al. 2014; Watson et al. 2016). A larger-scale study spanning 11 years found that anti–T. circumcincta pan-isotype antibody levels were associated with reduced breeding success in the subsequent year for both adult males and females of high body weight (Hayward et al. 2014). This suggests that variation in immune defenses in Soay sheep may be maintained by opposing selection on survival and fecundity and that selection may be sex and environment dependent.

In this study, we use a much larger 26-year anti–T. circumcincta antibody data set to comprehensively assess how selection via survival and fecundity operates on these immune markers, including tests for stabilizing, sex-, age-, and environment-dependent patterns of selection. We also tested the associations between antibody levels at the time of measurement and strongyle FEC (as a proxy for parasite burden) and body weight (as a proxy for body condition) and determined whether associations with survival and fecundity were statistically independent of these two measures. We predict that (1) there will be a negative association between antibody measures and strongyle FEC, indicative of their involvement in resistance to infection (Stear et al. 2009); (2) there will be a positive association between body weight and antibody levels, indicative of animals in better condition being able to mount stronger immune responses (Siva-Jothy and Thompson 2002; Seppälä and Jokela 2010); (3) selection on immune responses will be strongest in lambs and males because they have the highest parasite infection intensities and lower annual survival rates (Clutton-Brock and Pemberton 2004; Craig et al. 2008); and (4) selection will be stronger in years of high population density when infection risk is higher (Seppälä 2015).

## Methods

### Study Population

The Soay sheep is a primitive breed of domestic sheep that was isolated on the island of Soay in the remote St. Kilda archipelago several millennia ago and has been living under unmanaged conditions and evolving under natural selection since then (Clutton-Brock and Pemberton 2004). In 1932, 107 Soay sheep were moved to the larger island of Hirta after the evacuation of all human residents. Approximately one-third of the population of these sheep live in the Village Bay area of Hirta, and these individuals have been the subject of a long-term study since 1985 (Clutton-Brock and Pemberton 2004). In April of each year, around 95% of all lambs born in this area are caught and individually tagged. Each August, as many sheep as possible from the study population are recaptured using temporary traps (Clutton-Brock and Pemberton 2004). At capture, animals are weighed and blood and fecal samples are collected. Whole blood samples are collected into heparin tubes and centrifuged at 1,008 g for 10 min, and plasma are removed and stored at −20°C. The majority of sheep deaths occur over winter, and the population dynamics of the Soay sheep are characterized by periods of low but rising population sizes followed by high mortality crash winters in which more than half the population can die (Clutton-Brock and Pemberton 2004). Regular censuses and mortality searches during the winter months result in the majority of carcasses being discovered, providing accurate death date information.

Breeding success for study individuals is known with high precision but is calculated differently depending on sex. For female annual breeding success, we used the number of offspring born to the female in the subsequent spring (derived from observational data), while for males we calculated the number of offspring sired by the ram and born in the subsequent spring (derived using genetic pedigree construction). The pedigree was constructed using 315 unlinked single-nucleotide polymorphisms (linkage disequilibrium $r2<0.05$) with a minor allele frequency >0.4 using the R library sequoia (Huisman 2017). Total Village Bay population (referred to as population density) was used as a measure of environmental conditions because it captures annual variation in both food availability and nematode parasite exposure (Wilson et al. 2004; Hayward et al. 2014).

In this study, we included all animals that were caught and had plasma samples taken in August between 1990 and 2015. At the same time point, fecal samples are collected for the estimation of strongyle FEC. FEC is estimated from fecal samples as the number of eggs per gram using a modified McMaster technique (for full details, see Gulland and Fox 1992; Hayward et al. 2014).

### Measurement of Antibodies

All results were recorded as OD values. In order to minimize confounding of capture year and age effects with plate to plate variation, each plate included samples from 2 years paired at random with different age groups on each plate. All plates were run in duplicate, and duplicate sample ODs were removed if the coefficient of variation between duplicates was >0.2 or the absolute difference between the duplicates was >0.2 OD units. We also checked the correlation of ODs across duplicate plates and reran both plates if the correlation coefficient was <0.8. No duplicate plates failed this check, and the number of samples that failed quality control per assay was 13 for IgA (seven lambs and six adults), eight for IgE (six lambs and two adults), and 27 for IgG (five lambs and 22 adults).

To reduce error due to within-plate variation, we included two sample free wells per plate (50 µL TBST added instead of sheep plasma) as blanks and two wells of positive controls. Positive controls for the IgE assay were serum from ewes trickle infected with T. circumcincta and, for the IgA and IgG assay, plasma from normal healthy nonimmunized domestic sheep. For subsequent analyses, the mean optical density ratio of each sample was taken according to the formula

where the numerator was set to zero if the blank OD was greater than the sample OD (in order to avoid negative values).

### Statistical Analyses

The final data set used in the analyses encompassed 6,398 samples from 3,114 individuals whose age was known. Five samples from very late-born lambs caught in August that were within 50 days of birth were excluded from the data set because of the potential presence of maternal antibodies and differences in development stage to other lambs. The final data set is deposited in the Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.9qr0401 (Sparks et al. 2018).

All analyses were conducted using linear and generalized linear mixed models in R v3.4.4 (R Development Core Team 2016). Because of the immaturity of the immune system at 4 months (Smith et al. 1985) and the large increase in antibody levels between 4 and 16 months, we ran three models per response variable for lambs (both sexes), adult females, and adult males. Adults were separated into females and males because of the differences in range and distribution of phenotypic traits, such as weight and FEC in adults, as well as differences in survival rates and breeding success (). We built models with August weight, August strongyle FEC, overwinter survival, and annual breeding success as response variables with anti–T. circumcincta IgA, IgE, and IgG levels as explanatory variables. Continuous explanatory variables were rescaled to mean 0 and 1 SD before inclusion in each model subset. We first fitted a model containing the fixed effects previously identified as important predictors of these variables. We then simplified the model by stepwise deletion, sequentially removing fixed effects with the lowest t values and determining statistical significance using likelihood ratio tests until a base model containing only significant ($P<.05$) fixed effects was left. All dropped terms were then retested against this base model using the same criteria. Population size and sex (where fitted) were kept in all models in order to test for interactions with antibody levels. To the minimal model containing only significant terms, we separately added linear and quadratic functions of anti-Tc IgA, IgE, and IgG and assessed them for significance using likelihood ratio tests. Next we tested whether there were sex- or density-dependent effects of antibody levels in each of the models by looking at interactions with sex and population density against the linear and, if significant, quadratic antibody terms.

August strongyle FEC was heavily right skewed and zero inflated. Since the accuracy of the method used to estimate FEC is limited to 100 eggs per gram (Gulland and Fox 1992), we binned August FEC into multiples of 100, with zeros counted as 0 and values >6,000 eggs grouped into a single final bin (fig. A1A, A1B). Associations between antibody measures and August strongyle FEC were modeled using generalized linear mixed models via the glmmADMB package v0.8.3.3 (Fournier et al. 2012; Skaug et al. 2016). To assess which distribution best fitted the data set, we tested Poisson distributions and two negative binomial distributions (NB2 parameterization where $variance=μ(1+μ/k)$ and NB1 fit where $variance=φμ$) with and without zero inflation and compared the Akaike information criterion (AIC) between all six models. We found that the negative binomial NB2 distribution without zero inflation had the lowest AIC for all models of FEC: lambs, adult females, and adult males. We included fixed effects for sex (lamb model only), weight, age, and population density. Because previous studies of FEC in Soay sheep have identified curvilinear relationships between FEC and both age and weight (Hayward et al. 2009), we tested both linear and quadratic effects of these terms in our models. We included capture year as a random effect in both models and individual identity as a random effect in the adult models.

We investigated potential associations between August weight and anti-Tc antibody levels using linear mixed effects models via the lme4 package v1.1–17 (Bates et al. 2015), assuming a Gaussian distribution (fig. A1C, A1D). For the lamb models, we included sex, twin status, age in days, maternal age (quadratic), and population size as fixed effects and capture year and maternal identity as random effects. For adult weight models, we included age (linear and quadratic terms) and population size as fixed effects and individual identity and capture year as random effects.

We calculated overwinter survival as survival from capture in August to May 1 in the subsequent year using death date, capture, and census information. Survival probability varied with age and sex (fig. A1E) as well as among years, as previously documented (Clutton-Brock and Pemberton 2004). Analyses of survival were performed using generalized linear mixed models (GLMMs) in the lme4 package with a binomial error distribution and a logit link function. Binomial GLMMs were run with the built-in bobyqa optimizer, and the maximum number of iterations increased to 2E−6 to improve model convergence. For the lamb model, we included sex, twin status, weight (linear and quadratic), and population size as fixed effects and capture year as a random effect. For the adult models, we included age (linear and quadratic), weight (linear and quadratic), and population size as fixed effects, with individual identity and year as random effects.

For females that survived the winter, we calculated annual breeding success as the number of offspring born to a female in the subsequent spring from observational data. For male annual breeding success, we calculated the number of offspring sired by the ram and born in the subsequent spring using the pedigree. The male annual breeding success measure included males that died over the winter but excluded all males that were not seen to participate in the rut in the census records taken between October and December. The probability of producing a lamb was highly age and sex dependent (fig. A1F), with female and male lambs having particularly low probabilities of breeding (44% and 9%, respectively). Since the distributions of fecundity measures are very different for males and females (fig. A1F), sex-specific models were run for breeding success in both lambs and adults in this case. Female annual breeding success measures were treated as binary in both lambs and adults (lambed/did not lamb), while male fecundity measures were treated as binary only in lambs because of the low number of individuals siring lambs in their first year. All three models were run as a GLMM with a binomial error structure as for survival. For adult males, we analyzed annual breeding success in the glmmADMB package, where the best-fitting model was the negative binomial with the NB2 parameterization without zero inflation. For male and female lamb models, we included weight (linear and quadratic) and population size as fixed effects and year as a random effect. For adult male and female breeding success, we included age (linear and quadratic), weight (linear and quadratic), and population size as fixed effects and individual identity and year as random effects.

## Results

The three anti-Tc antibody isotypes were positively but weakly to moderately associated with one another, with higher correlations among isotypes observed in lambs than in adults ($r=0.28–0.42$ in lambs and 0.07–0.11 in adults; see fig. A2). We detail the associations found between antibody isotypes and strongyle FEC, weight, overwinter survival, and annual breeding success; a summary of these results can be found in figure A3.

### Associations between Antibody Levels and Strongyle FEC

In lambs, August strongyle FEC was weakly but significantly negatively associated with all three anti-Tc antibody isotype (IgA, IgG, and IgE) levels (Table 1; fig. 1A–1C). For all three isotypes, August FEC decreased linearly with increasing antibody levels. In a model containing all three antibody isotypes, IgA levels were independently negatively associated with FEC ($b=−0.074±0.021$ SE, $χ2(1)=11.860$, $P=.001$), while IgE and IgG were no longer significantly associated with FEC (IgE: $b=−0.030±0.020$ SE, $χ2(1)=2.180$, $P=.140$; IgG: $b=−0.034±0.022$ SE, $χ2(1)=2.420$, $P=.120$). We found no evidence for sex- or density-dependent effects of antibody levels on FEC in lambs (Table 1).

FemalesMales
EstimateSELRTPEstimateSELRTPEstimateSELRTP
Intercept1.556.069  .096.052  1.309.069
Sex (male).423.039113.220<.001
Twin.183.04914.240<.001
Age    −.633.111.579.172
Age2    .600.10234.440<.001−.411.1398.740.003
Weight−.912.127−.399.039101.340<.001−1.055.219
Weight2.700.13030.080<.001    .698.2289.320.002
Population size.234.06410.980.001−.039.048.640.424.054.064.700.403
Dropped:
Weight2      .900.343
IgA−.096.01925.460<.001−.348.120  .460.498
IgA2  .660.417.309.1216.560.010  .060.807
IgE−.059.0199.340.002  .200.655  1.220.269
IgE2  1.160.282  .840.359  .200.655
IgG−.075.01915.220<.001−.120.02916.880<.001−.089.0327.760.005
IgG2  2.140.144  .060.807  .360.549
Sex × IgA  .080.777
Sex × IgE  .200.655
Sex × IgG  .020.888
Population × IgA  2.800.094  1.760.185  3.600.058
Population × IgE  1.020.313  .040.842  2.840.092
Population × IgG  .040.842  3.520.061  .580.446
Population × IgA2      .420.517

Note. Included are the estimated effects (estimate), standard errors (SEs), and significance of fixed effects based on a likelihood ratio test (LRT; P) where df = 1. Dropped fixed effects show the significance of adding dropped terms back into the minimal model. For the added fixed effects, each antibody measure was added separately to the minimal model, including their linear and quadratic terms. Where the quadratic terms were significant, estimated effects are stated for both the linear and the quadratic terms combined, and interactions were tested in addition to both terms. Lambs: $n=1,850$ individuals; adult females: $n=2,694$ observations, $n=875$ individuals; adult males: $n=845$ observations, $n=451$ individuals.

View Table Image

In adults, there was a negative but weak linear association between anti-Tc IgG levels and August FEC for both males and females (Table 1; fig. 1E, 1F). In females, there was also a weak curvilinear relationship of IgA with FEC (Table 1; fig. 1D). A model with both significant antibody terms for females showed these two associations to be independent (IgG: $b=−0.116±0.029$ SE, $χ2(1)=15.780$, $P<.001$; IgA: $b(IgA2)=0.312±0.120$ SE, $b(IgA)=−0.339±0.119$ SE, $χ2(1)=6.720$, $P=.010$). Model predictions suggested that FEC was lowest at intermediate anti-Tc IgA levels in adult females (fig. 1D), and without the quadratic term in the model there was not a significant linear association between anti-Tc IgA and FEC ($b=−0.050±0.029$ SE, $χ2(1)=2.980$, $P=.084$). We found no evidence for density-dependent effects of antibody levels on FEC in adults (Table 1).

### Associations between Antibody Levels and Weight

There was a significant quadratic association between August weight and anti-Tc IgE levels in lambs (Table 2). Model predictions suggested that weight increased and then plateaued with increasing IgE levels in lambs (fig. 2A). Without the quadratic in the model, there was a significant positive linear association between IgE and weight ($b=0.128±0.042$ SE, $χ2(1)=9.399$, $P=.002$). This association was not driven by a few individuals with high anti-Tc IgE levels because a quadratic association remained after removing the 40 individuals with scaled IgE levels >3 ($b(IgE2)=−1.019±0.294$ SE, $b(IgE)=0.704±0.154$ SE, $χ2(1)=11.996$, $P=.001$). We found no association of the other antibody isotypes and weight in lambs or antibody × sex or density interactions (Table 2).

FemalesMales
EstimateSELRTPEstimateSELRTPEstimateSELRTP
Intercept12.869.221  21.290.130  27.194.232
Sex (males)1.448.079306.950<.001
Twin−3.299.112711.190<.001
Age (days).383.05154.670<.001
Age (years)    5.840.10714.191.328
Age (years)2    −4.337.1061259.800<.001−9.361.316555.040<.001
Maternal age4.236.166
Maternal age2−3.936.166482.220<.001
Population size−.409.2014.153.042−.549.11019.051<.001−.444.1875.543.019
IgA  .306.580.753.189.361.1337.275.007
IgA2  .161.688−.671.18812.820<.001  2.765.096
IgE.281.086.185.05212.654<.001  1.744.187
IgE2−.169.0834.191.041  .003.957  .001.971
IgG  .531.466.754.153  .292.589
IgG2  3.155.076−.660.14520.796<.001  1.905.168
Sex × IgA  .036.851
Sex × IgE  .794.373
Sex × IgG  1.788.181
Sex × IgE2  .441.507
Population × IgA  .023.881.077.0316.196.013  .192.662
Population × IgE  .152.697  .571.450  3.062.080
Population × IgG  1.100.294  .785.376  .088.766
Population × IgA2      .448.503
Population × IgE2  1.402.237
Population × IgG2      .708.400

Note. Included are the estimated effects (estimate), standard errors (SEs), and significance of fixed effects based on a likelihood ratio test (LRT; P) where df = 1. For the added fixed effects, each antibody measure was added separately to the minimal model, including their linear and quadratic terms. Where the quadratic terms were significant, estimated effects are stated for both the linear and the quadratic terms combined, and interactions were tested in addition to both terms. Lambs: $n=1,992$ individuals; adult females: $n=2,975$ observations, $n=918$ individuals; adult males: $n=927$ observations, $n=480$ individuals.

View Table Image

In adult females, there was a positive but weak linear association between August weight and anti-Tc IgE levels, while a significant curvilinear association was found between IgA and IgG levels and weight (Table 2; fig. 2B–2D). While higher IgE levels were associated with heaviest August weights, it was the intermediate anti-Tc IgA and IgG levels that were associated with heaviest weights (fig. 2B–2D). These associations were not driven by a few individuals with high anti-Tc antibody levels because the association remained after removal of individuals with high antibody levels of each isotype (after removal of IgA levels >2.5 [$n=15$]: $χ2(1)=10.235$, $P=.001$; after removal of IgE levels >2.5 [$n=61$]: $χ2(1)=6.981$, $P=.008$; after removal of IgG levels >3.75 [$n=5$]: $χ2(1)=26.621$, $P<.001$). Linear associations between weight and IgA and IgG were weaker and less significant than quadratic effects (IgA: $b=0.102±0.052$ SE, $χ2(1)=3.873$, $P=.049$; IgG: $b=0.085±0.044$ SE, $χ2(1)=3.628$, $P=.057$). When all the significant antibody terms were in the same model, all remained significant, suggesting that each antibody isotype had independent associations with August weight (IgA: $b(IgA2)=−0.621±0.187$ SE, $b(IgA)=0.667±0.189$ SE, $χ2(1)=11.075$, $P=.001$; IgE: $b=0.156±0.053$ SE, $χ2(1)=8.934$, $P=.003$; IgG: $b(IgG2)=−0.605±0.145$ SE, $b(IgG)=0.674±0.154$ SE, $χ2(1)=17.413$, $P<.001$). A marginal density-dependent association with IgA was found, indicating a slightly stronger increase with weight at high sheep densities (Table 2). The interaction remained significant when the other antibody terms were included in the model ($b=0.077±0.031$ SE, $χ2(1)=6.262$, $P=.012$). In adult males we found a positive but weak linear association of IgA with weight but no association with the other antibody isotypes or density interactions (Table 2; fig. 2E).

### Associations between Antibody Levels and Overwinter Survival

We found little evidence for associations between any of the three anti-Tc antibody levels and survival in lambs or adult males (Table 3). Although there were no significant interactions between antibody levels and sex for lambs, we found a marginally significant interaction between anti-Tc IgA level and population density. In adult females, we found that anti-Tc IgG and IgE levels were significant predictors of overwinter survival (Table 3; fig. 3). Females with higher IgG and IgE levels were more likely to survive the winter, independent of age and weight (Table 3; fig. 3). In a model with both IgG and IgE fitted, IgE became nonsignificant, but IgG remained significant (IgE: $b=0.149±0.083$ SE, $χ2(1)=3.301$, $P=.069$; IgG: $b=0.217±0.084$ SE, $χ2(1)=7.151$, $P=.007$). In addition, IgG remained significant if strongyle FEC was included in the model ($b=0.278±0.095$ SE, $χ2(1)=9.544$, $P=.002$) and if the one individual with unusually high IgG (>5) was dropped ($b=0.235±0.085$ SE, $χ2(1)=8.240$, $P=.004$). There was no evidence for environment-dependent selection on adult females (Table 3).

FemalesMales
EstimateSELRTPEstimateSELRTPEstimateSELRTP
Intercept.133.242  3.103.278  1.558.344
Sex (male)−1.033.11978.850<.001
Age (years)    −1.463.151261.080<.001−1.082.21439.792<.001
Weight.797.066166.440<.0011.415.137210.990<.0011.205.19251.924<.001
Population size−1.241.23020.141<.001−.764.2508.471.004−.929.3357.652.006
Dropped:
Twin  .920.337
Weight2  .441.506  3.793.051  .496.481
Age2      .320.571  .599.439
IgA−.093.0573.208.073  .386.535  2.444.118
IgA2  .004.952  .005.942  .014.906
IgE  .223.637.178.0824.893.027  <.001.996
IgE2  .115.734  .990.320  .118.731
IgG  .017.895.239.0848.743.003  <.001.990
IgG2  <.001.987  1.693.193  .159.690
Sex × IgA  1.625.202
Sex × IgE  .248.618
Sex × IgG  .242.623
Population × IgA−.129.0663.862.049  .013.908  .137.712
Population × IgE  .336.562  .039.843  1.073.300
Population × IgG  .594.441  1.790.181  .148.701

Note. Included are the estimated effects (estimate), standard errors (SEs), and significance of fixed effects based on a likelihood ratio test (LRT; P value) where df = 1. Dropped fixed effects show the significance of adding dropped terms back into the minimal model. For the added fixed effects, each antibody measure was added separately to the minimal model, including their linear and quadratic terms. Where the quadratic terms were significant, effects are stated for both the linear and the quadratic terms combined, and interactions were tested in addition to both terms. Lambs: $n=2,224$ individuals; adult females: $n=2,935$ observations, $n=903$ individuals; adult males: $n=856$ observations, $n=441$ individuals.

View Table Image

### Associations between Antibody Levels and Annual Breeding Success

In females, we found a negative association between anti-Tc IgE levels as lambs and the probability of females breeding in their first year (Table 4). Although our initial model detected a significant quadratic association between IgE and female first year breeding probability (Table 4), visual inspection of this relationship and the model predictions suggested that it could be driven by a handful of lambs with very high antibody levels. When the four individuals with high anti-Tc IgE values were dropped, the quadratic was no longer significant (after removal of IgE levels >4: $χ2(1)=0.906$, $P=.341$), but there was a significant negative linear relationship between IgE levels and first-year breeding probability ($b=−0.294±0.118$ SE, $χ2(1)=6.465$, $P=.011$; fig. A4A). We found no associations with any other antibody isotype or density interactions in female lambs. In males, we found a positive association between IgE levels as lambs and breeding success in their first year (Table 4). Male lambs with higher IgE antibody levels were more likely to sire offspring in their first year (fig. A4B), and this association remained after dropping the six individuals with high IgE levels ($b=0.353±0.113$ SE, $χ2(1)=9.037$, $P=.003$). We found no associations with any other antibody isotype or environment-dependent interactions in lambs.

Fixed effectsFemalesMales
EstimateSELRTPEstimateSELRTPEstimateSELRTPEstimateSELRTP
Intercept−.434.199  6.468.615  −3.100.213  −.241.123
Age    1.857.486    .203.06110.820.001
Age2    −2.349.43832.231<.001
Weight.614.10935.592<.0012.1481.0594.5401.560.785.076102.700<.001
Weight2    −2.1741.0554.062.044−3.5401.3997.954.005
Population size−.666.1949.553.002−.134.1251.069.301−.926.14526.737<.001−.240.1024.900.027
Dropped:
Age2              1.820.177
Weight2  2.477.116          .500.480
IgA  .766.382−1.176.648  2.619.106  1.760.185
IgA2  .410.5221.244.6494.075.044  .149.700  <.0011.000
IgE−.645.230  .600.439.262.0907.634.006  1.880.170
IgE2.561.2546.183.013  .916.339  1.492.222  .020.888
IgG  .087.768−1.025.557  .891.345.100.0504.040.044
IgG2  .481.488.991.5313.871.049  .329.566  .860.354
Population × IgA  .089.765  .981.322  .118.731  .220.639
Population × IgE  .839.360  .895.344  .054.816  1.760.185
Population × IgG  .003.955  .603.437  .665.415  3.400.065
Population × IgA2      1.130.288
Population × IgE2  .023.878
Population × IgG2      .787.375

Note. For lambs (both sexes) and adult females, the response is binary (lambed/did not lamb), while for adult males the response is a count of all lambs sired in the following spring. Included are the estimated effects (estimate), standard errors (SEs), and significance of fixed effects based on a likelihood ratio test (LRT; P) where df = 1. Dropped fixed effects show the significance of adding dropped terms back into the minimal model. For the added fixed effects, each antibody measure was added separately to the minimal model, including their linear and quadratic terms. Where the quadratic terms were significant, estimated effects are stated for both the linear and the quadratic terms combined, and interactions were tested in addition to both terms. Female lambs: $n=592$ individuals; female adults: $n=2,579$ observations, $n=792$ individuals; male lambs: $n=1,108$ individuals; male adults: $n=887$ observations, $n=457$ individuals.

View Table Image

In adult females, we found a curvilinear association of both anti-Tc IgA and anti-Tc IgG on breeding probability, but these associations were extremely weak and marginally significant (Table 4). There was no underlying linear association between either isotypes and annual breeding probability (IgA: $b=0.030±0.153$ SE, $χ2(1)=0.039$, $P=.843$; IgG: $b=−0.015±0.132$ SE, $χ2(1)=0.013$, $P=.908$). There was a weak and marginally significant positive association between IgG levels and the breeding success of adult males the following spring (Table 4). We found no significant associations between annual breeding success and any other antibody isotypes or density interactions in adults.

## Discussion

In this study, we have provided an unusually comprehensive assessment of how circulating antibody levels against a prevalent nematode parasite relate to parasite egg counts, host body weight, and host fitness in a wild mammal population. Our data offer a rare opportunity to explore the relationships between markers of immunity against an ecologically important parasite group and fitness in both sexes under a highly variable natural environment. Unsurprisingly, given the considerable environmental heterogeneity in natural systems such as ours, the patterns of associations we observe are complex and vary even among closely functionally related antibody isotypes (for a summary of the results, see fig. A3). Our results provide important evidence for age-, sex- and isotype-dependent selection on antihelminth antibody levels, which may help explain how genetic variation in immune phenotypes is maintained in wild populations. Furthermore, several interesting and important patterns emerge from our analyses, most saliently (1) a clear age-dependent switch in the antibody isotype most strongly associated with FEC and fitness, which may offer important insights into how selection is operating on antihelminth immunity in natural vertebrate populations; (2) little consistent evidence for the expected consistent positive association between antibody measures and host body weight; and (3) limited support for environment-dependent selection, strong stabilizing selection, or costs associated with reproduction with any of the antibody measures we examined. We discuss each of these results and their implications in more detail.

### Age-Dependent Associations among Antibodies, Strongyle Egg Counts, and Fitness

Generally, we found significant negative associations between anti-Tc antibody levels and strongyle FEC and positive associations with survival and reproduction, which is consistent with levels of these antibodies reflecting, in part, host resistance to nematode parasites rather than simply variation in parasite exposure (Hayward et al. 2014). However, the associations between antibody levels and strongyle FEC were typically weak, highlighting that numerous other factors beyond antibodies influence strongyle FEC, some of which have previously been identified in this study population (Wilson et al. 2004; Craig et al. 2008; Hayward et al. 2009, 2014). We also found evidence of a striking shift in the relationships between different antibody isotypes, worm egg counts, and overwinter survival probabilities with host age. IgA antibodies were the best predictor of FEC in lambs, but none of the isotypes predicted lamb overwinter survival, while in adults IgG antibodies were the only isotype to predict FEC and also were positively associated with female overwinter survival and weakly with male annual breeding success. This pattern may reflect an important shift in selection on the immune response to parasitic helminths from initial exposure and the development of an immune response during early life to the mature adult immune response, which is thought to be more geared toward tolerating infection, wound healing, and limiting damage (Allen and Maizels 2011). Our results suggest that while antiparasite IgA levels may represent the best marker of the development of immunity to these parasites during early life, IgG represents the best marker of an adult’s long-term ability to cope with persistent exposure and infection with nematode parasites. To our knowledge, our results provide the first evidence of such an age-related change in selection on immunity from the wild and provide an important basis for further studies to explore age-dependent changes in the immune response and its consequences for health and fitness under natural conditions.

Although all three antibody isotypes were moderately positively correlated in lambs (fig. A2) and were all significantly negatively associated with FEC when fitted separately in models, only anti-Tc IgA levels remained significant when all three isotypes were included in the model. This suggests shared covariance between isotypes and FEC, with IgA providing the strongest independent association. Veterinary studies of domestic sheep have demonstrated an important role of the mucosal IgA and IgE responses in the development of immunity to Teladorsagia circumcincta in lambs (Stear et al. 1995; Murphy et al. 2010). The hypersensitivity response (allergic inflammation caused by IgE-mediated mast cell degranulation) associated with production of antiparasite IgE is generally thought to be involved in the expulsion of larval parasites (McNeilly et al. 2009), while mucosal antiparasite IgA, which is positively correlated with circulating levels of antiparasite IgA (Henderson and Stear 2006), is thought to be involved in reducing adult worm growth and fecundity and inhibiting the development of worm larvae (Stear et al. 1995, 2004; Strain et al. 2002). Negative associations between FEC and circulating anti-Tc IgA levels have been previously reported in both domestic and Soay sheep lambs (Stear et al. 1995, 2004; Coltman et al. 2001; Strain et al. 2002). Despite its negative association with FEC, we found little evidence that anti-Tc IgA levels in summer predicted first-winter survival or the probability of breeding the following year. This may be due to the immature, dynamic nature of the immune response at 4 months old, which may not reflect their anti-Tc antibody levels going into winter (Smith et al. 1985). Overall, the absence of relationships between anti-Tc antibody levels and first-winter survival in lambs was surprising, especially given that a previous study of Soay sheep found a negative association between strongyle FEC and annual fitness in lambs but not in adults (Hayward et al. 2011). This seems likely to reflect the greater importance of variation in exposure to helminths for health and fitness in immunologically immature lambs and contrasts with the patterns observed in immunologically mature adult sheep.

In adult Soay sheep there was a significant negative relationship between anti-Tc IgG and FEC in both sexes and significant positive associations between anti-Tc IgG and overwinter survival in females and annual breeding success in males. Anti-Tc IgA did independently predict FEC in adult females, but the relationship was not linear and implied slightly raised FEC at the extremes of the antibody distribution. Two previous smaller-scale, cross-sectional studies of this population investigating samples from different years have found that anti-Tc IgG levels predict overwinter survival in adult females, independent of other immune measures and body weight (Nussey et al. 2014; Watson et al. 2016). Observing the relationship again in our much larger-scale longitudinal study suggests that this association is observed only in adult females but is robust and independent of environmental conditions and individual body mass. Although studies on domestic ruminants have focused mainly on young animals (Stear et al. 2009), negative associations between anti-Tc IgG levels and FEC or worm number have been reported for adults (Williams et al. 2010; McBean et al. 2016), and there is a clear role for IgG in mediating protective immunity against helminths in laboratory mice (Blackwell and Else 2001; McCoy et al. 2008). It has been hypothesized that in early life, when animals are most vulnerable to helminth-associated morbidity, an aggressive immune response to limit worm establishment may develop, while in later life, once immune responses have matured and helminths have established infections, a more tolerant immune phenotype that limits worm fecundity, promotes healing of damaged tissues, and prevents immunopathology would be favored by natural selection (Allen and Maizels 2011; Medzhitov et al. 2012). Age-dependent selection on a tolerance-associated immune phenotype has been documented in wild voles, where increased GATA3 gene expression (a master regulator of the Th2 response) was positively associated with body condition and survival but only in older and larger males (Jackson et al. 2014). Our results provide support for age-dependent changes in the role of different antibody isotypes in resistance and tolerance to nematode parasites under natural infection. Further detailed immunological study is now required to understand why IgG antibodies provide the best marker of an individual’s ability to cope with helminth infection in adulthood but not during early life.

### Sex- and Environment-Dependent Selection on Antibody Levels

We found little evidence for sex differences in the associations between antibody levels and FEC or body weight, but marginal significant trends in our analysis suggest that relationships with survival and fecundity might differ between males and females. Perhaps most strikingly, in lambs there was a positive association between anti-Tc IgE and the probability of males producing offspring in the first year of life but a negative association in females (fig. A4). IgE antibodies are associated with a hypersensitivity response and mast cell degranulation in the gut, which is an effective mediator of nematode expulsion. It is also associated with inappetence, diarrhea, and reduced weight gain in domestic lambs (Stear et al. 2009). Given that males grow more rapidly than females in this species and body size is under strong sexual selection in male Soay sheep (Clutton-Brock and Pemberton 2004), it is surprising that an immune measure potentially associated with a cost to weight gain should be associated with improved reproductive performance in males but not females. We also found evidence for a weak and marginally significant positive association between anti-Tc IgG levels and adult male annual breeding success, while this isotype predicted overwinter survival in adult females. In wild ruminants generally, males are thought to invest a large amount of energetic resources in the autumn rut at a clear cost to their overwinter survival prospects, while females are more conservative and balance investment between their developing fetus and their own maintenance over winter (Clutton-Brock and Pemberton 2004). If antibody levels in summer reflect an individual’s ability to resist nematode infections and limit the damage they cause, this may reduce the cost of infection in autumn and winter, freeing up resources that adult males would be expected to invest in the rut, while females invest in their own maintenance over winter. Overall, these results may reflect important and previously unknown differences in the costs and benefits between the sexes of mounting an aggressive hypersensitivity response to infection during early life and mounting a more tolerance-oriented response in adulthood. These trends warrant further investigation in both wild and domestic ruminant systems, involving more regular seasonal monitoring of reproductive hormones, immunity, and worm burdens across seasons in both sexes.

Given the considerable annual variation in environmental conditions and nematode parasite exposure experienced by Soay sheep on St. Kilda that is associated with often dramatic changes in population size (Clutton-Brock and Pemberton 2004), the general lack of evidence for interactions between antibody measures and sheep density in our models was surprising. We would predict the benefits, in terms of parasite burden and fitness, of improved immune responses to nematodes to be most strongly felt in years of high sheep density, when animals are in the poorest conditions and exposure to nematode larvae is highest (Hayward et al. 2014). However, there was very limited evidence for such relationships. This lack of environment-dependent associations may reflect the fact that anti-Tc antibody levels in Soay sheep appear highly repeatable in adulthood (Hayward et al. 2014; A. M. Sparks, unpublished data). Although there has been a far greater theoretical and empirical focus on environment-dependent selection on immune phenotypes as a force maintaining genetic variation in the immune response in natural systems (Svensson et al. 2001; Calsbeek et al. 2008; Graham et al. 2010; Seppälä 2015), our analyses provide far stronger support for age-dependent rather than environment-dependent selection. This argues for more theoretical and empirical effort being directed at investigating the evolutionary and ecological causes and consequences of such age-dependent variation in selection on immunity in the future.

### Reproductive Costs and Condition Dependence

Unlike a number of experimental studies in wild birds, we found no strong evidence for the predicted negative association between parasite-specific antibodies and subsequent reproductive success, indicative of a reproductive cost of immunity (Ilmonen et al. 2000; Bonneaud et al. 2003; Marzal et al. 2007; Gasparini et al. 2009). There was a negative association between IgE levels and the probability of female sheep breeding in their first year; however, this result was only marginally significant. A previous study in Soay sheep found that high levels of self-reactive antinuclear antibodies were associated with reduced breeding success in both adult male and female Soay sheep (Graham et al. 2010). Furthermore, previous work using an anti-Tc antibody ELISA that simultaneously captures all isotypes (pan-isotype) found that adult males with high antibody levels were less likely to sire offspring the following year, and females of high body weight and high antibody levels were also less likely to have a lamb the subsequent year (Hayward et al. 2014). This study also found a marginal association between anti-Tc antibodies and survival in female lambs only (Hayward et al. 2014). This pan-isotype assay is expected to largely capture variation in IgG, which is the predominant isotype in circulation. However, these results contrast with the results of the current and previous isotype-specific studies (Nussey et al. 2014; Watson et al. 2016), suggesting that it is actually providing a different immunological signal to IgG when measured alone. One interesting possibility is that the pan-isotype measure better captures overall protein costs associated with antiparasite immune responses in general because it captures total immunoglobulin levels directed at these parasites. On the other hand, specific isotypes offer a clearer signal of specific immune responses directed at resisting or tolerating worm infection and thus are more strongly associated with survival. This difference among the current and previous study further highlights the complex and variable patterns of natural selection that can be observed even on apparently functionally closely related immune measures within the same study population.

There is a strong expectation for immune responses to be tightly linked with the condition of the host, with individuals in a better nutritional state more able to invest resources in immunity and therefore to resist infection (van Noordwijk and de Jong 1986; Kraaijeveld and Godfray 1997; Reznick et al. 2000). Interestingly, we found evidence for both positive linear and curvilinear relationships between body weight, depending on the age or sex group, and antibody isotype. However, these associations were typically weak, suggesting that anti-Tc antibody levels are not consistently or strongly positively associated with body weight of the host. A number of other studies have shown that body mass and food intake are closely related to immune responses (Siva-Jothy and Thompson 2002; Seppälä and Jokela 2010). However, in a previous study of a smaller sample of adult female Soay sheep, we did not find associations between anti-Tc antibody levels and weight (Nussey et al. 2014). Similarly, a recent study found that protein markers expected to be associated with nutritional state, total protein, and albumin were only weakly associated with body mass (Garnier et al. 2017). While body weight is often used as a proxy for condition in field studies, both our results and previous studies suggest that the relationship between body weight and levels of proteins in the blood that we expect to be intimately linked to nutritional state and physiological condition may be variable and weak (Hayward et al. 2014; Garnier et al. 2017). While our results highlight the point that using a few catch-all measures of immunocompetence is very unlikely to capture the complexity of the immune response and selection acting on it in wild systems (Demas et al. 2011), the complex relationships observed with body weight emphasize the need for biomarkers of nutritional state and condition within ecology (Garnier et al. 2017). Our study shows that closely related immune traits may have very different associations with health and fitness measures and suggests that multivariate immune phenotypes must be measured in order to fully understand how genetic variation is maintained in these traits in the wild.

We thank National Trust for Scotland, Scottish Natural Heritage, QinetiQ, Eurest, and Kilda Cruises for logistical support on St. Kilda. We further thank the many field volunteers and all those who have contributed to keeping the project going over many years, including T. Clutton-Brock, B. Grenfell, T. Coulson, L. Kruuk, and M. Crawley as well as Ian Stevenson for data management support. We also thank Jisca Huisman for constructing the genetic pedigree used in this study. We are very grateful to Andrea Graham for comments on an early version of the manuscript. The Soay sheep project is supported by grants from the Natural Environment Research Council; D.H.N. was supported by a Biotechnology and Biological Sciences Research Council David Phillips fellowship, and A.M.S. was supported by a Medical Research Council PhD studentship.

## Literature Cited

• Allen, J. E., and R. M. Maizels. 2011. Diversity and dialogue in immunity to helminths. Nature Reviews Immunology 11:375–388.

• Altizer, S., A. Dobson, P. Hosseini, P. Hudson, M. Pascual, and P. Rohani. 2006. Seasonality and the dynamics of infectious diseases. Ecology Letters 9:467–484.

• Appleton, J. A., and D. D. McGregor. 1987. Characterization of the immune mediator of rapid expulsion of Trichinella spiralis in suckling rats. Immunology 62:477–484.

• Bates, D., M. Mächler, B. M. Bolker, and S. C. Walker. 2015. Fitting linear mixed-effects models using lme4. Journal of Statistical Software 67:1–48.

• Blackwell, N. M., and K. J. Else. 2001. B cells and antibodies are required for resistance to the parasitic gastrointestinal nematode Trichuris muris. Infection and Immunity 69:3860–3868.

• Bonneaud, C., J. Mazuc, G. Gonzalez, C. Haussy, O. Chastel, B. Faivre, and G. Sorci. 2003. Assessing the cost of mounting an immune response. American Naturalist 161:367–379.

• Calsbeek, R., C. Bonneaud, and T. B. Smith. 2008. Differential fitness effects of immunocompetence and neighbourhood density in alternative female lizard morphs. Journal of Animal Ecology 77:103–109.

• Christe, P., F. De Lope, G. González, N. Saino, and A. P. Møller. 2001. The influence of environmental conditions on immune responses, morphology and recapture probability of nestling house martins (Delichon urbica). Oecologia 126:333–338.

• Christe, P., A. P. Møller, and F. de Lope. 1998. Immunocompetence and nestling survival in the house martin: the tasty chick hypothesis. Oikos 83:175–179.

• Clutton-Brock, T., and J. Pemberton. 2004. Soay sheep: dynamics and selection in an island population. Cambridge University Press, Cambridge.

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

• Craig, B. H., J. G. Pilkington, and J. M. Pemberton. 2006. Gastrointestinal nematode species burdens and host mortality in a feral sheep population. Parasitology 133:485–496.

• Craig, B. H., L. J. Tempest, J. G. Pilkington, and J. M. Pemberton. 2008. Metazoan-protozoan parasite co-infections and host body weight in St Kilda Soay sheep. Parasitology 135:433–441.

• Demas, G. E., D. A. Zysling, B. R. Beechler, M. P. Muehlenbein, and S. S. French. 2011. Beyond phytohaemagglutinin: assessing vertebrate immune function across ecological contexts. Journal of Animal Ecology 80:710–730.

• Esser-von Bieren, J., I. Mosconi, R. Guiet, A. Piersgilli, B. Volpe, F. Chen, W. C. Gause, A. Seitz, J. S. Verbeek, and N. L. Harris. 2013. Antibodies trap tissue migrating helminth larvae and prevent tissue damage by driving IL-4Rα-independent alternative differentiation of macrophages. PLoS Pathogens 9:e1003771.

• Falconer, D. S., and T. F. C. Mackay. 1996. Introduction to quantitative genetics. 4th ed. Pearson Prentice Hall, Essex.

• Fournier, D. A., H. J. Skaug, J. Ancheta, J. Ianelli, A. Magnusson, M. N. Maunder, A. Nielsen, and J. Sibert. 2012. AD Model Builder: using automatic differentiation for statistical inference of highly parameterized complex nonlinear models. Optimization Methods and Software 27:233–249.

• Garnier, R., C. K. Cheung, K. A. Watt, J. G. Pilkington, J. M. Pemberton, and A. L. Graham. 2017. Joint associations of blood plasma proteins with overwinter survival of a large mammal. Ecology Letters 20:175–183.

• Gasparini, J., P. Bize, R. Piault, K. Wakamatsu, J. D. Blount, A.-L. Ducrest, and A. Roulin. 2009. Strength and cost of an induced immune response are associated with a heritable melanin-based colour trait in female tawny owls. Journal of Animal Ecology 78:608–616.

• Gonzalez, G., G. Sorci, A. P. Møller, P. Ninni, C. Haussy, and F. De Lope. 1999. Immunocompetence and condition-dependent sexual advertisement in male house sparrows (Passer domesticus). Journal of Animal Ecology 68:1225–1234.

• Graham, A. L., J. E. Allen, and A. F. Read. 2005. Evolutionary causes and consequences of immunopathology. Annual Review of Ecology, Evolution, and Systematics 36:373–397.

• Graham, A. L., A. D. Hayward, K. A. Watt, J. G. Pilkington, J. M. Pemberton, and D. H. Nussey. 2010. Fitness correlates of heritable variation in antibody responsiveness in a wild mammal. Science 330:662–665.

• Graham, A. L., D. M. Shuker, L. C. Pollitt, S. K. J. R. Auld, A. J. Wilson, and T. J. Little. 2011. Fitness consequences of immune responses: strengthening the empirical framework for ecoimmunology. Functional Ecology 25:5–17.

• Gulland, F. M. 1992. The role of nematode parasites in Soay sheep (Ovis aries L.) mortality during a population crash. Parasitology 105:493–503.

• Gulland, F. M., and M. Fox. 1992. Epidemiology of nematode infections of Soay sheep (Ovis aries L.) on St Kilda. Parasitology 105:481–492.

• Harris, N., and W. Gause. 2011. To B or not to B: B cells and the Th2-type immune response to helminths. Trends in Immunology 32:80–88.

• Harris, N. L., I. Spoerri, J. F. Schopfer, C. Nembrini, P. Merky, J. Massacand, J. F. Urban, et al. 2006. Mechanisms of neonatal mucosal antibody protection. Journal of Immunology 177:6256–6262.

• Hayward, A. D., R. Garnier, K. A. Watt, J. G. Pilkington, B. T. Grenfell, J. B. Matthews, J. M. Pemberton, D. H. Nussey, and A. L. Graham. 2014. Heritable, heterogeneous, and costly resistance of sheep against nematodes and potential feedbacks to epidemiological dynamics. American Naturalist 184(suppl.):S58–S76.

• Hayward, A. D., A. J. Wilson, J. G. Pilkington, T. H. Clutton-Brock, J. M. Pemberton, and L. E. B. Kruuk. 2011. Natural selection on a measure of parasite resistance varies across ages and environmental conditions in a wild mammal. Journal of Evolutionary Biology 24:1664–1676.

• Hayward, A. D., A. J. Wilson, J. G. Pilkington, J. M. Pemberton, and L. E. B. Kruuk. 2009. Ageing in a variable habitat: environmental stress affects senescence in parasite resistance in St Kilda Soay sheep. Proceedings of the Royal Society B 276:3477–3485.

• Henderson, N. G., and M. J. Stear. 2006. Eosinophil and IgA responses in sheep infected with Teladorsagia circumcincta. Veterinary Immunology and Immunopathology 112:62–66.

• Huisman, J. 2017. Pedigree reconstruction from SNP data: parentage assignment, sibship clustering and beyond. Molecular Ecology Resources 17:1009–1024.

• Huntley, J. F., J. Redmond, W. Welfare, G. Brennan, F. Jackson, F. Kooyman, and L. Vervelde. 2001. Studies on the immunoglobulin E responses to Teladorsagia circumcincta in sheep: purification of a major high molecular weight allergen. Parasite Immunology 23:227–235.

• Ilmonen, P., T. Taarna, and D. Hasselquist. 2000. Experimentally activated immune defence in female pied flycatchers results in reduced breeding success. Proceedings of the Royal Society B 267:665–670.

• Jackson, J. A., A. J. Hall, I. M. Friberg, C. Ralli, A. Lowe, M. Zawadzka, A. K. Turner, et al. 2014. An immunological marker of tolerance to infection in wild rodents. PLoS Biology 12:e1001901.

• Kraaijeveld, A. R., and H. C. J. Godfray. 1997. Trade-off between parasitoid resistance and larval competitive ability in Drosophila melanogaster. Nature 389:278–280.

• Lazzaro, B. P., and T. J. Little. 2009. Immunity in a variable world. Philosophical Transactions of the Royal Society B 364:12–26.

• Lindström, K. M., J. Foufopoulos, H. Pärn, and M. Wikelski. 2004. Immunological investments reflect parasite abundance in island populations of Darwin’s finches. Proceedings of the Royal Society B 271:1513–1519.

• Lochmiller, R. L., and C. Deerenberg. 2000. Trade-offs in evolutionary immunology: just what is the cost of immunity? Oikos 88:87–98.

• Maizels, R. M., and D. H. Nussey. 2013. Into the wild: digging at immunology’s evolutionary roots. Nature Immunology 14:879–883.

• Marzal, A., M. Reviriego, F. de Lope, and A. P. Møller. 2007. Fitness costs of an immune response in the house martin (Delichon urbica). Behavioral Ecology and Sociobiology 61:1573–1580.

• McBean, D., M. Nath, F. Kenyon, K. Zile, D. J. Bartley, and F. Jackson. 2016. Faecal egg counts and immune markers in a line of Scottish cashmere goats selected for resistance to gastrointestinal nematode parasite infection. Veterinary Parasitology 229:1–8.

• McCoy, K. D., M. Stoel, R. Stettler, P. Merky, K. Fink, B. M. Senn, C. Schaer, et al. 2008. Polyclonal and specific antibodies mediate protective immunity against enteric helminth infection. Cell Host and Microbe 4:362–373.

• McNeilly, T. N., E. Devaney, and J. B. Matthews. 2009. Teladorsagia circumcincta in the sheep abomasum: defining the role of dendritic cells in T cell regulation and protective immunity. Parasite Immunology 31:347–356.

• Medzhitov, R., D. S. Schneider, and M. P. Soares. 2012. Disease tolerance as a defense strategy. Science 335:936–941.

• Merino, S., A. P. Møller, and F. de Lope. 2000. Seasonal changes in cell-mediated immunocompetence and mass gain in nestling barn swallows: a parasite-mediated effect? Oikos 90:327–332.

• Murphy, K. M. 2012. Janeway’s immunobiology. 8th ed. Garland Science, Abingdon.

• Murphy, L., P. D. Eckersall, S. C. Bishop, J. J. Pettit, J. F. Huntley, R. Burchmore, and M. J. Stear. 2010. Genetic variation among lambs in peripheral IgE activity against the larval stages of Teladorsagia circumcincta. Parasitology 137:1249–1260.

• Nunn, C. L., J. L. Gittleman, and J. Antonovics. 2000. Promiscuity and the primate immune system. Science 290:1168–1170.

• Nussey, D. H., K. A. Watt, A. Clark, J. G. Pilkington, J. M. Pemberton, A. L. Graham, and T. N. McNeilly. 2014. Multivariate immune defences and fitness in the wild: complex but ecologically important associations among plasma antibodies, health and survival. Proceedings of the Royal Society B 281:20132931.

• Palacios, M. G., D. W. Winkler, K. C. Klasing, D. Hasselquist, and C. M. Vleck. 2011. Consequences of immune system aging in nature: a study of immunosenescence costs in free-living tree swallows. Ecology 92:952–966.

• Parejo, D., and N. Silva. 2009. Immunity and fitness in a wild population of Eurasian kestrels Falco tinnunculus. Naturwissenschaften 96:1193–1202.

• Råberg, L., J.-A. Nilsson, P. Ilmonen, M. Stjernman, and D. Hasselquist. 2000. The cost of an immune response: vaccination reduces parental effort. Ecology Letters 3:382–386.

• Råberg, L., and M. Stjernman. 2003. Natural selection on immune responsiveness in blue tits Parus caeruleus. Evolution 57:1670–1678.

• R Development Core Team. 2016. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna.

• Reznick, D., L. Nunney, and A. Tessier. 2000. Big houses, big cars, superfleas and the costs of reproduction. Trends in Ecology and Evolution 15:421–425.

• Saino, N., A. M. Bolzern, and A. P. Møller. 1997. Immunocompetence, ornamentation, and viability of male barn swallows (Hirundo rustica). Proceedings of the National Academy of Sciences of the USA 94:549–552.

• Scharsack, J. P., M. Kalbe, C. Harrod, and G. Rauch. 2007. Habitat-specific adaptation of immune responses of stickleback (Gasterosteus aculeatus) lake and river ecotypes. Proceedings of the Royal Society B 274:1523–1532.

• Schmid-Hempel, P. 2011. Evolutionary parasitology: the integrated study of infections, immunology, ecology, and genetics. Oxford University Press, New York.

• Seppälä, O. 2015. Natural selection on quantitative immune defence traits: a comparison between theory and data. Journal of Evolutionary Biology 28:1–9.

• Seppälä, O., and J. Jokela. 2010. Maintenance of genetic variation in immune defense of a freshwater snail: role of environmental heterogeneity. Evolution 64:2397–2407.

• Sheldon, B. C., and S. Verhulst. 1996. Ecological immunology: costly parasite defences and trade-offs in evolutionary ecology. Trends in Ecology and Evolution 11:317–321.

• Simon, A. K., G. A. Hollander, and A. McMichael. 2015. Evolution of the immune system in humans from infancy to old age. Proceedings of the Royal Society B 282:20143085.

• Siva-Jothy, M. T., and J. J. W. Thompson. 2002. Short-term nutrient deprivation affects immune function. Physiological Entomology 27:206–212.

• Skaug, H., D. Fournier, B. Bolker, A. Magnusson, and A. Nielsen. 2016. Generalized linear mixed models using AD Model Builder. R package version 0.8.3.3.

• Smith, W. D., F. Jackson, E. Jackson, and J. Williams. 1985. Age immunity to Ostertagia circumcincta: comparison of the local immune responses of 4 1/2- and 10-month-old lambs. Journal of Comparative Pathology 95:235–245.

• Sparks, A. M., K. Watt, R. Sinclair, J. G. Pilkington, J. M. Pemberton, S. E. Johnston, T. N. McNeilly, and D. H. Nussey. 2018. Data from: Natural selection on antihelminth antibodies in a wild mammal population. American Naturalist, Dryad Digital Repository, http://dx.doi.org/10.5061/dryad.9qr0401.

• Stear, M. J., K. Bairden, G. T. Innocent, S. Mitchell, S. Strain, and S. C. Bishop. 2004. The relationship between IgA activity against 4th-stage larvae and density-dependent effects on the number of 4th-stage larvae of Teladorsagia circumcincta in naturally infected sheep. Parasitology 129:363–369.

• Stear, M. J., S. C. Bishop, M. Doligalska, J. L. Duncan, P. H. Holmes, J. Irvine, L. McCririe, Q. A. McKellar, E. Sinski, and M. Murray. 1995. Regulation of egg production, worm burden, worm length and worm fecundity by host responses in sheep infected with Ostertagia circumcincta. Parasite Immunology 17:643–652.

• Stear, M. J., B. Boag, I. Cattadori, and L. Murphy. 2009. Genetic variation in resistance to mixed, predominantly Teladorsagia circumcincta nematode infections of sheep: from heritabilities to gene identification. Parasite Immunology 31:274–282.

• Strain, S. A. J., S. C. Bishop, N. G. Henderson, A. Kerr, Q. A. McKellar, S. Mitchell, and M. J. Stear. 2002. The genetic control of IgA activity against Teladorsagia circumcincta and its association with parasite resistance in naturally infected sheep. Parasitology 124:545–552.

• Svensson, E., B. Sinervo, and T. Comendant. 2001. Density-dependent competition and selection on immune function in genetic lizard morphs. Proceedings of the National Academy of Sciences of the USA 98:12561–12565.

• Uller, T., C. Isaksson, and M. Olsson. 2006. Immune challenge reduces reproductive output and growth in a lizard. Functional Ecology 20:873–879.

• van Noordwijk, A. J., and G. de Jong. 1986. Acquisition and allocation of resources: their influence on variation in life history tactics. American Naturalist 128:137–142.

• Viney, M. E., E. M. Riley, and K. L. Buchanan. 2005. Optimal immune responses: immunocompetence revisited. Trends in Ecology and Evolution 20:665–669.

• Watson, R. L., T. N. McNeilly, K. A. Watt, J. M. Pemberton, J. G. Pilkington, M. Waterfall, P. R. T. Hopper, D. Cooney, R. Zamoyska, and D. H. Nussey. 2016. Cellular and humoral immunity in a wild mammal: variation with age and sex and association with overwinter survival. Ecology and Evolution 6:8695–8705.

• Williams, A. R., D. G. Palmer, I. H. Williams, P. E. Vercoe, D. L. Emery, and L. J. E. Karlsson. 2010. Relationships between immune indicators of parasitic gastroenteritis, nematode burdens and faecal dry matter in sheep. Animal Production Science 50:219–227.

• Wilson, K., B. T. Grenfell, J. G. Pilkington, H. E. G. Boyd, and F. M. Gulland. 2004. Parasites and their impact. Pages 17–51 in T. H. Clutton-Brock and J. M. Pemberton, eds. Soay sheep: dynamics and selection in an island population. Cambridge University Press, Cambridge.

• Zuk, M. 1990. Reproductive strategies and disease susceptibility: an evolutionary viewpoint. Parasitology Today 6:231–233.

Associate Editor: Gregory E. Demas

Editor: Daniel I. Bolnick