Joint associations of blood plasma proteins with overwinter survival of a large mammal.

In many wild animal populations, hosts are at risk of parasites and malnutrition and resource costs of defence may be difficult to afford. We postulate that proteins, important in homeostasis and immunity, play a complex but central role in condition dependence and resource costs of mammalian immune defence. To test this, we measured plasma concentrations of albumin, total proteins. Self-reactive antibodies and parasite-specific IgG in female Soay sheep. Using a principal component analysis, we found a new metric of condition reflecting individual variation in acquisition, assimilation and/or recycling of plasma proteins that predicted overwinter survival. Controlling for this metric, an age-dependent trade-off between antibody titres and protein reserves emerged, indicating costs of mounting an antibody response: younger individuals survived best when prioritising immunity while older individuals fared better when maintaining high-protein nutritional plane. These findings suggest fascinating roles for protein acquisition and allocation in influencing survival in wild animal populations.

Natural selection on individual variation in tolerance of gastrointestinal nematode infection.

Hosts may mitigate the impact of parasites by two broad strategies: resistance, which limits parasite burden, and tolerance, which limits the fitness or health cost of increasing parasite burden. The degree and causes of variation in both resistance and tolerance are expected to influence host-parasite evolutionary and epidemiological dynamics and inform disease management, yet very little empirical work has addressed tolerance in wild vertebrates. Here, we applied random regression models to longitudinal data from an unmanaged population of Soay sheep to estimate individual tolerance, defined as the rate of decline in body weight with increasing burden of highly prevalent gastrointestinal nematode parasites. On average, individuals lost weight as parasite burden increased, but whereas some lost weight slowly as burden increased (exhibiting high tolerance), other individuals lost weight significantly more rapidly (exhibiting low tolerance). We then investigated associations between tolerance and fitness using selection gradients that accounted for selection on correlated traits, including body weight. We found evidence for positive phenotypic selection on tolerance: on average, individuals who lost weight more slowly with increasing parasite burden had higher lifetime breeding success. This variation did not have an additive genetic basis. These results reveal that selection on tolerance operates under natural conditions. They also support theoretical predictions for the erosion of additive genetic variance of traits under strong directional selection and fixation of genes conferring tolerance. Our findings provide the first evidence of selection on individual tolerance of infection in animals and suggest practical applications in animal and human disease management in the face of highly prevalent parasites.

Marker-dependent associations among oxidative stress, growth and survival during early life in a wild mammal.

Oxidative stress (OS) is hypothesized to be a key physiological mechanism mediating life-history trade-offs, but evidence from wild populations experiencing natural environmental variation is limited. We tested the hypotheses that increased early life growth rate increases OS, and that increased OS reduces first-winter survival, in wild Soay sheep (Ovis aries) lambs. We measured growth rate and first-winter survival for four consecutive cohorts, and measured two markers of oxidative damage (malondialdehyde (MDA), protein carbonyls (PC)) and two markers of antioxidant (AOX) protection (total AOX capacity (TAC), superoxide dismutase (SOD)) from blood samples. Faster lamb growth was weakly associated with increased MDA, but not associated with variation in the other three markers. Lambs with higher SOD activity were more likely to survive their first winter, as were male but not female lambs with lower PC concentrations. Survival did not vary with MDA or total TAC. Key predictions relating OS to growth and survival were therefore supported in some OS markers, but not others. This suggests that different markers capture different aspects of the complex relationships between individual oxidative state, physiology and fitness, and that overarching hypotheses relating OS to life-history variation cannot be supported or refuted by studying individual markers.

Heritable, heterogeneous, and costly resistance of sheep against nematodes and potential feedbacks to epidemiological dynamics.

Infected hosts may preserve fitness by resisting parasites (reducing parasite burden) and/or tolerating them (preventing or repairing infection-induced damage). Theory predicts that these individual-level defense strategies generate divergent population-level feedbacks that would maintain genetic heterogeneity for resistance but purge heterogeneity for tolerance. Because resistance reduces parasite abundance, selection for costly resistance traits will weaken as resistance becomes common. Such negative frequency-dependent selection contrasts with predictions for tolerance, which maintains parasite abundance and so is expected to generate positive frequency-dependent selection, unless, for example, tolerance trades off with resistance. Thus far, there have been few tests of this theory in natural systems. Here, we begin testing the predictions in a mammalian field system, using data on individual gastrointestinal nematode burdens, nematode-specific antibody titers (as a resistance metric), the slope of body weight on parasite burden (as a tolerance metric), and fitness from an unmanaged population of Soay sheep. We find that nematode resistance is costly to fitness and underpinned by genetic heterogeneity, and that resistance is independent of tolerance. Drawing upon empirical metrics such as developed here, future work will elucidate how resistance and tolerance feedbacks interact to generate population-scale patterns in the Soay sheep and other field systems.

Multivariate immune defences and fitness in the wild: complex but ecologically important associations among plasma antibodies, health and survival.

Despite our rapidly advancing mechanistic understanding of vertebrate immunity under controlled laboratory conditions, the links between immunity, infection and fitness under natural conditions remain poorly understood. Antibodies are central to acquired immune responses, and antibody levels circulating in vivo reflect a composite of constitutive and induced functional variants of diverse specificities (e.g. binding antigens from prevalent parasites, self tissues or novel non-self sources). Here, we measured plasma concentrations of 11 different antibody types in adult females from an unmanaged population of Soay sheep on St Kilda. Correlations among antibody measures were generally positive but weak, and eight of the measures independently predicted body mass, strongyle parasite egg count or survival over the subsequent winter. These independent and, in some cases, antagonistic relationships point to important multivariate immunological heterogeneities affecting organismal health and fitness in natural systems. Notably, we identified a strong positive association between anti-nematode immunoglobulin (Ig) G antibodies in summer and subsequent over-winter survival, providing rare evidence for a fitness benefit of helminth-specific immunity under natural conditions. Our results highlight both the evolutionary and ecological importance and the complex nature of the immune phenotype in the wild.

A mixed-model approach for estimating drivers of microbiota community composition and differential taxonomic abundance.

Next-generation sequencing (NGS) and metabarcoding approaches are increasingly applied to wild animal populations, but there is a disconnect between the widely applied generalized linear mixed model (GLMM) approaches commonly used to study phenotypic variation and the statistical toolkit from community ecology typically applied to metabarcoding data. Here, we describe the suitability of a novel GLMM-based approach for analyzing the taxon-specific sequence read counts derived from standard metabarcoding data. This approach allows decomposition of the contribution of different drivers to variation in community composition (e.g., age, season, individual) via interaction terms in the model random-effects structure. We provide guidance to implementing this approach and show how these models can identify how responsible specific taxonomic groups are for the effects attributed to different drivers. We applied this approach to two cross-sectional data sets from the Soay sheep population of St. Kilda. GLMMs showed agreement with dissimilarity-based approaches highlighting the substantial contribution of age and minimal contribution of season to microbiota community compositions, and simultaneously estimated the contribution of other technical and biological factors. We further used model predictions to show that age effects were principally due to increases in taxa of the phylum Bacteroidetes and declines in taxa of the phylum Firmicutes. This approach offers a powerful means for understanding the influence of drivers of community structure derived from metabarcoding data. We discuss how our approach could be readily adapted to allow researchers to estimate contributions of additional factors such as host or microbe phylogeny to answer emerging questions surrounding the ecological and evolutionary roles of within-host communities. IMPORTANCE NGS and fecal metabarcoding methods have provided powerful opportunities to study the wild gut microbiome. A wealth of data is, therefore, amassing across wild systems, generating the need for analytical approaches that can appropriately investigate simultaneous factors at the host and environmental scale that determine the composition of these communities. Here, we describe a generalized linear mixed-effects model (GLMM) approach to analyze read count data from metabarcoding of the gut microbiota, allowing us to quantify the contributions of multiple host and environmental factors to within-host community structure. Our approach provides outputs that are familiar to a majority of field ecologists and can be run using any standard mixed-effects modeling packages. We illustrate this approach using two metabarcoding data sets from the Soay sheep population of St. Kilda investigating age and season effects as worked examples.

Lactoferrin quantification in cattle faeces by ELISA.

BACKGROUND: Promoting and maintaining health is critical to ruminant welfare and productivity. Within human medicine, faecal lactoferrin is quantified for routine assessment of various gastrointestinal illnesses avoiding the need for blood sampling. This approach might also be adapted and applied for non-invasive health assessments in animals. METHODS: In this proof-of-concept study, a bovine lactoferrin enzyme-linked immunosorbent assays (ELISA), designed for serum and milk, was applied to a faecal supernatant to assess its potential for quantifying lactoferrin in the faeces of cattle. Faecal lactoferrin concentrations were compared to background levels to assess the viability of the technique. A comparison was then made against serum lactoferrin levels to determine if they were or were not reflective of one another. RESULTS: The optical densities of faecal samples were significantly greater than background readings, supporting the hypothesis that the assay was effective in quantifying faecal lactoferrin (T 13, 115 = 11.99, p < 0.0005). The mean faecal lactoferrin concentration was 0.269 µg mL-1 (S.E. 0.031) and the mean serum concentration 0.074 µg mL-1 (S.E. 0.005). Lactoferrin concentrations of faecal and serum samples, taken from the same animals on the same day, were significantly different (T 21 = 2.20, p = 0.039) and did not correlate (r = 0.2699, p = 0.238). CONCLUSION: Results support the hypothesis that lactoferrin can be quantified in cattle faeces by ELISA. Whilst further research is required to determine the physiological source of the lactoferrin, this highlights the potential of the method for non-invasive assessment of cattle immunology and pathology.

A candidate gene approach to study nematode resistance traits in naturally infected sheep.

Sheep naturally acquire a degree of resistant immunity to parasitic worm infection through repeated exposure. However, the immune response and clinical outcome vary greatly between animals. Genetic polymorphisms in genes integral to differential T helper cell polarization may contribute to variation in host response and disease outcome. A total of twelve single nucleotide polymorphisms (SNPs) were sequenced in IL23R, RORC2 and TBX21 from genomic DNA of Scottish Blackface lambs. Of the twelve SNPs, six were non-synonymous (missense), four were within the 3' UTRs and two were intronic. The association between nine of these SNPs and the traits of body weight, faecal egg count (FEC) and relative T. circumcincta L3-specific IgA antibody levels was assessed in a population of domestic Scottish Blackface ewe lambs and a population of free-living Soay ewe lambs both naturally infected with a mixture of nematodes. There were no significant associations identified between any of the SNPs and phenotypes recorded in either of the populations after adjustment for multiple testing (Bonferroni corrected P value≤0.002). In the Blackface lambs, there was a nominally significant association (P=0.007) between IL23R p.V324M and weight at 20 weeks. This association may be worthy of further investigation in a larger sample of sheep.

Fecal antibody levels as a noninvasive method for measuring immunity to gastrointestinal nematodes in ecological studies.

Among-individual variation in antibody-associated immunity to gastrointestinal nematode parasites (GIN) is known be associated with life-history traits and vital rates in wild vertebrate systems. To date, measurement of levels of antibodies against GIN antigens in natural populations has exclusively been based on invasive blood sampling techniques. Previous work in laboratory rodents and ruminant livestock suggests that antibody measures from feces may provide a viable noninvasive approach. We measured total and anti-GIN antibodies of different isotypes (immunoglobulin (Ig) G, IgA and IgE) from paired samples of plasma and feces from free-living Soay sheep of different ages and sexes. We tested the correlations among these measures as well as their associations with body mass and Strongyle nematode fecal egg counts (FEC). Significant positive correlations were present among plasma and fecal anti-GIN antibody levels for IgG and IgA. Generally, correlations between total antibody levels in plasma and feces were weaker and not significant. No significant relationships were found between any antibody measures and body mass; however, fecal anti-GIN antibody levels were significantly negatively correlated with FEC. Our data clearly demonstrate the feasibility of measuring anti-GIN antibodies from fecal samples collected in natural populations. Although associations of fecal antibody levels with their plasma counterparts and FEC were relatively weak, the presence of significant correlations in the predicted direction in a relatively small and heterogeneous sample suggests fecal antibody measures could be a useful, noninvasive addition to current eco-immunological studies.

Cellular and humoral immunity in a wild mammal: Variation with age & sex and association with overwinter survival.

Immune defenses are expected to be crucial for survival under the considerable parasite pressures experienced by wild animals. However, our understanding of the association between immunity and fitness in nature remains limited due to both the complexity of the vertebrate immune system and the often-limited availability of immune reagents in nonmodel organisms. Here, we use methods and reagents developed by veterinary researchers for domestic ungulates on blood samples collected from a wild Soay sheep population, to evaluate an unusually broad panel of immune parameters. Our evaluation included different innate and acquired immune cell types as well as nematode parasite-specific antibodies of different isotypes. We test how these markers correlate with one another, how they vary with age-group and sex, and, crucially, whether they predict overwinter survival either within or among demographic groups. We found anticipated patterns of variation in markers with age, associated with immune development, and once these age trends were accounted for, correlations among our 11 immune markers were generally weak. We found that females had higher proportions of naïve T cells and gamma-delta T cells than males, independent of age, while our other markers did not differ between sexes. Only one of our 11 markers predicted overwinter survival: sheep with higher plasma levels of anti-nematode IgG antibodies were significantly more likely to survive the subsequent high mortality winter, independent of age, sex, or weight. This supports a previous finding from this study system using a different set of samples and shows that circulating antibody levels against ecologically relevant parasites in natural systems represent an important parameter of immune function and may be under strong natural selection. Our data provide rare insights into patterns of variation among age- and sex groups in different T-cell subsets and antibody levels in the wild, and suggest that certain types of immune response-notably those likely to be repeatable within individuals and linked to resistance to ecologically relevant parasites-may be most informative for research into the links between immunity and fitness under natural conditions.

Vitamin D status predicts reproductive fitness in a wild sheep population.

Vitamin D deficiency has been associated with the development of many human diseases, and with poor reproductive performance in laboratory rodents. We currently have no idea how natural selection directly acts on variation in vitamin D metabolism due to a total lack of studies in wild animals. Here, we measured serum 25 hydroxyvitamin D (25(OH)D) concentrations in female Soay sheep that were part of a long-term field study on St Kilda. We found that total 25(OH)D was strongly influenced by age, and that light coloured sheep had higher 25(OH)D3 (but not 25(OH)D2) concentrations than dark sheep. The coat colour polymorphism in Soay sheep is controlled by a single locus, suggesting vitamin D status is heritable in this population. We also observed a very strong relationship between total 25(OH)D concentrations in summer and a ewe's fecundity the following spring. This resulted in a positive association between total 25(OH)D and the number of lambs produced that survived their first year of life, an important component of female reproductive fitness. Our study provides the first insight into naturally-occurring variation in vitamin D metabolites, and offers the first evidence that vitamin D status is both heritable and under natural selection in the wild.

Plasma markers of oxidative stress are uncorrelated in a wild mammal.

Oxidative stress, which results from an imbalance between the production of potentially damaging reactive oxygen species versus antioxidant defenses and repair mechanisms, has been proposed as an important mediator of life-history trade-offs. A plethora of biomarkers associated with oxidative stress exist, but few ecological studies have examined the relationships among different markers in organisms experiencing natural conditions or tested whether those relationships are stable across different environments and demographic groups. It is therefore not clear to what extent studies of different markers can be compared, or whether studies that focus on a single marker can draw general conclusions regarding oxidative stress. We measured widely used markers of oxidative damage (protein carbonyls and malondialdehyde) and antioxidant defense (superoxide dismutase and total antioxidant capacity) from 706 plasma samples collected over a 4-year period in a wild population of Soay sheep on St Kilda. We quantified the correlation structure among these four markers across the entire sample set and also within separate years, age groups (lambs and adults), and sexes. We found some moderately strong correlations between some pairs of markers when data from all 4 years were pooled. However, these correlations were caused by considerable among-year variation in mean marker values; correlation coefficients were small and not significantly different from zero after accounting for among-year variation. Furthermore, within each year, age, and sex subgroup, the pairwise correlation coefficients among the four markers were weak, nonsignificant, and distributed around zero. In addition, principal component analysis confirmed that the four markers represented four independent axes of variation. Our results suggest that plasma markers of oxidative stress may vary dramatically among years, presumably due to environmental conditions, and that this variation can induce population-level correlations among markers even in the absence of any correlations within contemporaneous subgroups. The absence of any consistent correlations within years or demographic subgroups implies that care must be taken when generalizing from observed relationships with oxidative stress markers, as each marker may reflect different and potentially uncoupled biochemical processes.

Fitness correlates of heritable variation in antibody responsiveness in a wild mammal.

A functional immune system is important for survival in natural environments, where individuals are frequently exposed to parasites. Yet strong immune responses may have fitness costs if they deplete limited energetic resources or cause autoimmune disease. We have found associations between fitness and heritable self-reactive antibody responsiveness in a wild population of Soay sheep. The occurrence of self-reactive antibodies correlated with overall antibody responsiveness and was associated with reduced reproduction in adults of both sexes. However, in females, the presence of self-reactive antibodies was positively associated with adult survival during harsh winters. Our results highlight the complex effects of natural selection on immune responsiveness and suggest that fitness trade-offs may maintain immunoheterogeneity, including genetic variation in autoimmune susceptibility.