Figures
Abstract
Sex hormones, reproductive status, and pathogen load all affect stress. Together with stress, these factors can modulate the immune system and affect disease incidence. Thus, it is important to concurrently measure these factors, along with their seasonal fluctuations, to better understand their complex interactions. Using steroid hormone metabolites from fecal samples, we examined seasonal correlations among zebra and springbok stress, reproduction, gastrointestinal (GI) parasite infections, and anthrax infection signatures in zebra and springbok in Etosha National Park (ENP), Namibia, and found strong seasonal effects. Infection intensities of all three GI macroparasites examined (strongyle helminths, Strongyloides helminths, and Eimeria coccidia) were highest in the wet season, concurrent with the timing of anthrax outbreaks. Parasites also declined with increased acquired immune responses. We found hormonal evidence that both mares and ewes are overwhelmingly seasonal breeders in ENP, and that reproductive hormones are correlated with immunosuppression and higher susceptibility to GI parasite infections. Stress hormones largely peak in the dry season, particularly in zebra, when parasite infection intensities are lowest, and are most strongly correlated with host mid-gestation rather than with parasite infection intensity. Given the evidence that GI parasites can cause host pathology, immunomodulation, and immunosuppression, their persistence in ENP hosts without inducing chronic stress responses supports the hypothesis that hosts are tolerant of their parasites. Such tolerance would help to explain the ubiquity of these organisms in ENP herbivores, even in the face of their potential immunomodulatory trade-offs with anti-anthrax immunity.
Citation: Cizauskas CA, Turner WC, Pitts N, Getz WM (2015) Seasonal Patterns of Hormones, Macroparasites, and Microparasites in Wild African Ungulates: The Interplay among Stress, Reproduction, and Disease. PLoS ONE 10(4): e0120800. https://doi.org/10.1371/journal.pone.0120800
Academic Editor: Srinand Sreevatsan, University of Minnesota, UNITED STATES
Received: September 24, 2014; Accepted: February 6, 2015; Published: April 15, 2015
Copyright: © 2015 Cizauskas et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: Data is now available online from Dryad, doi:10.5061/dryad.cr24d.
Funding: Funding was provided by National Institutes of Health Ecology and Evolution of Infectious Diseases grant #GM83863 to WMG (http://www.fic.nih.gov/programs/Pages/ecology-infectious-diseases.aspx); Andrew and Mary Thompson Rocca Scholarships to CAC (http://africa.berkeley.edu/fellowships/Rocca.php); University of California, Berkeley, Department of Environmental Science, Policy, and Management block grants and continuing fellowships to CAC (http://ourenvironment.berkeley.edu/current-graduate-students/funding/); and University of California, Berkeley, Graduate Division Summer Grants to CAC (http://grad.berkeley.edu/news/egrad/february-2013/summer-grant-2013/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Stressors can take the form of acute, unpleasant, and usually unpredictable stimuli that precipitate a stress response. The resulting hormonal, physiological, and behavioral changes help animals cope with the stressor and enhance survival [1,2]. Glucocorticoids (GCs), adrenal hormones of the hypothalamic-pituitary-adrenal (HPA) axis that are the primary mediators of the stress response, are helpful during acute stressors but can be problematic when elevated for several hours a day over weeks to months [1]. This chronic increase in GCs has several negative physiological effects, including metabolic disruption, reproductive suppression, immunosuppression and immunomodulation, and increased susceptibility to disease [3–5].
Pathogens themselves can also cause chronic stress in their hosts. While evidence for the ability of macroparasites to cause host stress is equivocal, several studies have found significant, positive relationships between parasite infection intensity or prevalence and GC concentration (e.g. [6,7]). Common gastrointestinal (GI) nematode parasites, such as strongyle and Strongyloides helminths, can have significant negative fitness impacts on hosts, such as impairing growth and causing weight loss, undernutrition, and reduced fecundity [8,9]. Coccidian GI parasites can negatively affect total body mass, impair growth and food utilization, cause epithelial cell damage and inflammation, and result in diarrhea, weight loss, and even death [9,10]. This interplay between stress and disease, however—particularly between stress and macroparasites—has been understudied in natural systems and is poorly understood.
Etosha National Park (ENP), Namibia, is a highly seasonal environment; it experiences one rainy season for approximately three to four months per year, followed by a cold dry and then a hot dry season [11,12]. Two of the predominant plains ungulates in ENP, plains zebra (Equus quagga) and springbok (Antidorcas marsupialis) are affected by these seasonal changes. While there is nearly 100% prevalence of strongyle nematodes in zebra in this system year-round, zebra experience a significant increase in strongyle infection intensity in the wet season compared to in the dry [13]. Springbok in ENP also experience significantly greater infection intensity in three gastrointestinal parasites (strongyle nematodes, Strongyloides nematodes, and Eimeria coccidia) during the wet season [13]. Zebra and springbok also experience annual anthrax outbreaks in the wet season, caused by the bacterium Bacillus anthracis [14]; during this season, it has been estimated that up to 400 zebras die of anthrax, out of a population of approximately 13,000 animals (estimations of springbok anthrax mortality have not been conducted using these methods, but springbok and zebra are the primary anthrax hosts in ENP) [15]. ENP anthrax outbreaks occur shortly after parasite infection intensities peak [16,17], and we have found that GI parasites may contribute to increased zebra susceptibility to this pathogen at this time due to immunomodulatory effects [18].
Hosts can experience other seasonally-driven potential stressors as well. Dry seasons often bring nutritional depletion [17], potentially leading to increased GC secretion [19,20]. While both zebra and springbok are not strictly seasonal breeders, both typically give birth during the wet season; zebra often mate again soon thereafter, while springbok often mate in the cold dry season (July). Thus, most mares and ewes experience mid-gestation during the hot dry season [11,13,21–23]. Reproductive behaviors such as rutting [24], and gestation and lactation [25] are particularly energy demanding, and most last for an extended period, potentially leading to chronic stress [26]. Reproductive hormones themselves can also be immunosuppressive, potentially influencing susceptibility to disease [27].
While studies have examined the impact of chronic stress in wildlife and the interaction between stress and seasonality, or reproduction, or disease individually (e.g. [26,28–30], we are unaware of any that have examined the seasonality of stress responses in combination with environmental changes, reproductive efforts, pathogen pressures and interactions, and immune function. Because all of these factors are closely interlinked, it is important that studies begin examining this complex array of factors in concert with each other. This study examines the seasonal fluctuations in stress hormones, in conjunction with environmental parameters, reproductive and gestational status, and pathogen interactions in zebra and springbok in ENP. As there is such strong seasonality in this system for most parameters examined, we hypothesized that stress responses would be significantly higher in the wet season compared to in the dry season: zebra appear to engage in most of their intense reproductive behaviors in the wet season (birth, lactation, estrus), springbok largely birth early in and lactate through the wet season, and both zebra and springbok experience significantly higher infection intensities of macroparasite and microparasite infections during this time. In addition, because we have previously found that zebra GI parasites likely exert more selection pressure than do the other pathogens in this system [31], and that they influence host immunity during peak infection intensities [18], we hypothesized that high levels of parasitism would correspond most strongly with stress hormone responses out of all the factors examined. Understanding the primary drivers of host stress is important for determining the full extent of pathogen effects, as well as for determining potential avenues for intervention in conservation programs.
Methods
Study Area and Species
Etosha National Park (ENP) is a 22,915 km2 fenced conservation area in northern Namibia, located between 18°30’S-19°30’S and 14°15’E-17°10’E (Fig. A in S1 File). Rainfall in ENP is highly seasonal: the rainy season lasts from December through April, with the greatest rainfall occurring during January and February, and 80% of all rain falling between December and March [11,32] (Fig. B in S1 File). The only perennial water available to the park’s wildlife is found in man-made boreholes, or in natural artesian or contact springs [33]. Permits for this research were issued by the Namibian Ministry of Environment and Tourism (MET) and the Etosha Ecological Institute (EEI).
Plains zebra (Equus quagga) and springbok (Antidorcas marsupialis; members of the family Bovidae) are the two most abundant plains ungulate species in ENP, with populations fluctuating around 13,000 (95% CI rounded to nearest 100: 10,900–15,000) and 15,600 (13,200–17,900) (EEI unpublished aerial survey data 2005), though these numbers are likely an underestimate for springbok. Neither of the species involved in this research is endangered nor protected.
Zebra Capture and Sampling
We obtained whole blood, serum, fecal, and ectoparasite samples from collared zebra over five seasons between 2008 and 2010 (Table A in S1 File). We immobilized, sampled, and released all animals safely under animal handling protocol AUP R217-0509B (University of California, Berkeley). A Namibian Ministry of Environment and Tourism veterinarian was present for all captures, as was author Cizauskas, who is also a licensed veterinarian. Animals were anesthetized for a short time (5–10 minutes), using a combination of etorphine, azaperone, and detomidine delivered via Pneu-Darts (Williamsport, PA). Anesthesia was monitored in the field and animals were given reversal medication at the end of sample collection and observed until they awoke and moved off normally. The work in Etosha National Park, Namibia, was approved by the Namibian Ministry of Environment and Tourism (MET) and the Etosha Ecological Institute (EEI). We originally sampled 38 animals in the wet season, intending to recapture them as many times as possible in five subsequent seasons roughly six months apart to both control for individual variation (Table A in S1 File). If we could not re-capture all initially collared individuals, we captured and collared new individuals in order to have a similar sample size for each season. We conducted a total of 144 zebra capture events; all captured animals for this study were adult females to control for sex differences. We further determined age to half a year by combining tooth eruption observations, caliper measurements of upper incisors, and patterns of wear, and adding on the time between subsequent captures and the time of first dental evaluation for each individual [34]. We recorded, when possible, whether a mare was actively caring for a foal (0 for no foal; 1 for foal present); whether a mare was pregnant (0 for not pregnant; 1 for possibly pregnant based on palpation and visual determination; and 2 for definitely pregnant); and whether a mare was lactating (0 for no lactation; 1 for watery, straw-colored liquid able to be expressed from teats; 2 for milk able to be expressed from teats).
We sampled 69 individuals overall, with 20 sampled twice, 11 sampled three times, 12 sampled four times, and two sampled five times (Table B in S1 File). We collected blood from peripheral veins for blood smears and hematocrit (HCT) analysis, and removed serum for antibody analyses. From blood smears we obtained neutrophil and lymphocyte counts (S1 File). Neutrophils are white blood cells instrumental in the innate immune and inflammatory responses against primarily extracellular pathogens and possibly against certain intracellular bacteria [35]. Lymphocytes include the different types of T and B cells that are activated in an antigen-specific manner as components of the adaptive immune system [36].
From serum we measured anti-anthrax antibody titers using enzyme-linked immunosorbent assays (ELISAs) (S1 File). Antibodies against the protective antigen (PA) toxin component of anthrax have been shown to be essential for adaptive protection against anthrax [37,38]. As a bacterial infection, anthrax provokes primarily a Th1-type immune response, in contrast to the Th2-type response involved in macroparasite infections [39]. We used an ELISA procedure with wildtype Bacillus anthracis PA as coating antigen to measure host anti-PA titers ([40]; S1 File).
When possible, we collected feces by observing an individual prior to capture and collecting a homogenized sub-sample within ten minutes of defecation. For animals that were not observed defecating, we collected fecal samples when possible by inserting a gloved hand into the rectum. As capture events took place between 9:00 and 14:00, fecal samples were collected within this same time window, thereby controlling for potential differences in timing of fecal egg shedding [41].
We collected and counted all visible ectoparasites, regardless of life stage, on zebras during capture events. The vast majority of ticks observed were Rhipicephalus evertsi mimeticus, a tick species found throughout Namibia in wild equids and greater kudu [42]. These ticks parasitize hosts year-round, with more adults present from November to May and immature stages peaking from February to March and May to September [43].
Non-Captured Zebra and Springbok Sampling
Fecal samples were collected noninvasively from zebra and springbok in the Okaukuejo area approximately monthly between August 2005 and May 2007 [12,13]. Fecal collection took place in the first week of each month, with a mean of 37 and 35 samples collected per month for zebra and springbok, respectively. Samples were collected within 100m of roads and between 7:00 and 13:00. For each fecal sample, date, time, species, sex, and age category were recorded. Age and sex were determined within one of three age classes: juveniles for those <1 year old; yearlings for those 1–2 years of age; and adults for those 2+ years old using physical characteristics as previously described [13]. Samples were refrigerated immediately and evaluated for GI parasites within 48hrs, after which they were frozen until hormone extraction (within 18 months of sampling).
We subsampled this large sample set for hormone analyses. We selected only yearlings and adults to control for effects of sexual maturity or near sexual maturity on hormone profiles. We selected male and female samples, and zebra and springbok samples in approximately equal numbers. We chose samples representing similar seasons to those of captured zebra samplings: February-May for the "wet" season, and August-October for the "dry" season. Ultimately, 318 zebra (22 male yearlings, 11 female yearlings, 136 male adults, and 149 female adults) and 272 springbok extracted samples (22 male yearlings, 35 female yearlings, 106 male adults, and 107 female adults) were used for hormone analysis.
Rainfall Quantification
Previous studies in this system determined that helminth infection intensity in zebras and springbok is significantly related to rainfall one and two months prior, with both strongly correlated with helminth eggs and coccidian oocysts shed per gram of feces [12]. We determined cumulative rainfall over the two months prior to sampling by adding up daily rainfall amounts (in mm) over the 60 days prior to each individual sampling event. We linked each fecal sample to daily rainfall records from the closes management station to the sampling location. For some analyses, we also grouped individuals into wet or dry season categories, based on rainfall quantification (see Supplementary Methods, “Rain Group Determination,” and Fig. B in S1 File).
Gastrointestinal Parasite Species and Quantification
The gastrointestinal (GI) nematodes examined in these host species were strongyle helminth parasites (abbreviated in models and figures as GIP), Strongyloides helminths (abbreviated as GIS), and Eimeria coccidia (abbreviated as GIC) (S1 File). Zebra experience infections with strongyles, while springbok experience infections with all three parasite types [13]. All of these parasites exhibit direct life cycles with free-living infectious stages that require 1–2 weeks of development in the external environment prior to becoming infectious to new hosts. These environmental stages are highly susceptible to desiccation [44,45]; thus, it is unsurprising that previous studies in ENP found a strongly seasonal pattern in GI parasite infection intensities, with hosts exhibiting greater new infections during the wet season than in the dry [13].
We evaluated all fecal samples for strongyle eggs, and springbok samples also for Strongyloides eggs and Eimeria oocysts using a modified McMaster technique for fecal egg counts [46], a commonly used non-invasive method for quantifying parasitism [47] (S1 File).
Hormone Metabolite Extraction and Measurement
Steroid hormones can be detected in blood and their metabolites can be detected in feces; however, measurements of hormone metabolites in feces are often best in wildlife because fecal collection can often be done noninvasively and is free from stress hormone feedback mechanisms, and fecal samples reflect longer term patterns of hormone trends, dampening potential acute stressor responses or diurnal variations in hormone levels [48,49]. Fecal assays for concentrations of stress and reproductive hormone metabolites have been developed and used extensively in wildlife (e.g. [50,51]). Fecal hormone metabolite concentrations reflect the state of circulating hormones an individual experienced approximately 10–12 hours prior for most ruminants [52], or approximately one day prior for equids [53]. Elevations in hormone levels were interpreted as reflecting host physiological status over at least the previous week; the probability that we collected a sample exactly 12–24 hours after an animal had experienced an acute, transient hormone fluctuation was extremely low, and thus we made the assumption that elevations in hormone metabolites reflected a longer physiological state.
We extracted fecal hormone metabolites from samples using methods modified from a previous similar study [54] (S1 File). We conducted radioimmunoassays for fecal glucocorticoid metabolites (FGM), fecal progestin metabolites (FPM), fecal estrogen metabolites (FEM), and fecal testosterone metabolites (FTM). Hereafter, these will be referred to simply as stress hormones, progesterone, estrogen, and testosterone or by the acronyms listed above. Commercially available kits were used to conduct these assays, and assays were validated for each species and each hormone (S1 File).
Statistical Analyses
Multiple Imputation of Missing Data—Captured Zebras.
Imputation is a method of replacing missing observations with plausible estimates based on available data. Multiple imputation methods are particularly useful for imputing missing multivariate data, and rely on available data from multiple predictors and covariates to create a set of datasets for each missing value [55,56]. This method, unlike single imputation methods, provides a variance of an estimate and accounts for uncertainty, given that the value in question was imputed rather than observed [57,58]. While multiple imputation is common in public health research, it has rarely thus far been used in ecological studies.
We imputed missing values from an individual capture event (e.g. we imputed a missing value for estrogen concentration for an individual if the majority of other values were present for that individual at that time); however, if an animal was not captured during one season, we did not impute missing values for that individual in that season. We imputed two missing values for HCT (1.4% of total HCT data), 15 (9.7%) for ectoparasite count, 33 (22.9%) for GIP, 47 (32.6%) for FGM, 59 (40.9%) for FPM, 53 (36.8%) for FEM, 16 (11.1%) for presence or absence of a foal, 21 (14.6%) for pregnancy status, and 16 (11.1%) for lactation status. Few of these missing variables were overlapping for the same individual-capture; thus, discarding cases with imputed data would have resulted in disregarding upwards of 90% of the capture events in our analyses, thereby discarding all the valuable information associated with individuals that had some missing data.
For imputation, we used the Multiple Imputation by Chained Equations (MICE) method with the 'mice' package [56] in R v2.15.2 [59]. We performed 100 imputation cycles and generated five imputations. We averaged the five estimates for each data point to produce single mean estimates for use in t tests, while we used the multiply imputed datasets directly in our generalized estimating equations and then determined mean estimates for the model parameters [60]. We adjusted standard errors, Wald statistics, and p values according to Rubin's rules [61]. See S1 File for a further discussion of multiple imputation and methods used.
We did not impute missing data for non-captured (NC) zebra or for springbok analyses, due to both the large number of missing hormone measurements and the fact that most missing variables were missing for the same individuals and/or same seasons. For NC zebra, there were 5 missing FGM (1.6% of total FGM data), 92 missing FPM (28.8%), 99 missing FEM (30.9%), and 55 missing FTM measurements (17.2%), nearly all missing from the same 50–100 individuals; for springbok, there were 2 missing FGM (0.73%), 78 missing FPM (28.6%), 102 missing FEM (37.4%), and 51 missing FTM measurements (18.8%), nearly all missing from the same 50–80 individuals. Instead of imputing missing variables, we partitioned these data in subsets based on the parameters of interest and performed several nested analyses of models with decreasing overall sample size (see below).
Seasonal Comparisons.
For each group of animals examined, we wanted to determine if immune factors, pathogen intensity or exposure, stress hormones, and reproductive hormone levels changed seasonally. We first examined the effects of rainfall for various pathogen and immune factors by comparing data between wet and dry seasons from the first sampling event for each captured zebra. This ensured that we were only comparing unique individuals, and therefore population-level effects, between seasons, without the possibility of autocorrelation and the issue of repeated measures (Table Ba in S1 File). We transformed variables by taking logarithmic, square-root, and fourth-root values to improve assumptions of normality when necessary and possible (Table C in S1 File). We compared normalized GI parasite counts, WBC counts, neutrophil and lymphocyte counts, HCT, and FGM, FEM, and FPM concentrations between seasons using Welch's two-tailed t tests. We compared ectoparasite counts and log2 anti-PA titers between rain groups using two-sided Wilcoxon rank sum tests. We adjusted p values to control for the familywise error rate by using the Holm-Bonferroni method [62]. We then did similar pairwise comparisons across the same animals, between seasons, for first and second captures to examine individual level effects between seasons given individual-level variation that otherwise might have obscured patterns at the population level (Table Bb in S1 File).
Similar to the seasonal comparisons for captured zebra, we performed bulk group comparisons for NC zebra and springbok between rain groups, sexes and age classes. For those response variables that could be successfully transformed to normality (GIP and FTM for NC zebra; FPM, FEM, and FTM for springbok), we performed Type III ANOVAs for unbalanced data to explore seasonal differences and potential interaction effects. We used R v2.15.2 and the 'car' package for these analyses [63]. We used Tukey pairwise post-hoc tests to determine the significance of multiple comparisons. For other response variables that were incapable of being transformed to complete normality, we used Type II ANOVAs and Tukey post-hoc tests to explore potentially significant interaction effects. We then used Wilcoxon rank sum tests to examine rain group, sex, and age differences, and performed these same comparison tests for any interactions significant in the ANOVAs at a p≤0.1 level. We adjusted p values from Wilcoxon rank sum tests to control for the familywise error rate by using the Holm-Bonferroni method discussed above.
For all GI parasite and FGM comparisons, we included interactions between sex and age, rainfall and sex, and rainfall and age. For FPM, FEM, and FTM, we did not include interactions between sex and age; numbers of yearlings were quite low for this dataset for NC zebras (N = 9), and these reproductive hormones were expected to be higher in adults of reproductive age in both species than in yearlings. We did, however, explore interactions between rainfall and sex for FPM and FEM. For FTM, only males were examined. The yearling age group in the NC zebra dataset was also quite low (N = 15), so we explored interactions between rainfall and age with caution. The yearling males in the springbok dataset only included 10 individuals, and thus we did not explore rainfall-age interactions for springbok testosterone.
Coinfections, Immunity, and Hormone Relationships
Captured Zebra.
We wanted to determine how seasonal factors, hormones, and coinfections affected the prevalence and intensity of each of our pathogens. Thus, we allowed each parasite to play the role of response variable in separate models to avoid biasing the directionality of parasite interactions and to allow each parasite outcome to be examined from the standpoint of potential predisposing immune factors toward that infection alone (Table D in S1 File). For these same reasons, we then used each of our immune parameters and the FGM measurements as the response variable in separate models to allow us to more directly examine the cross-relationships between them and the other physiological parameters and parasites (Table D in S1 File). We did not fit individual models with FPM, FEM, or FTM as response variables because these hormones are more clearly directionally determined by maturity, mating, and pregnancy status. These hormones were measured as means of controlling for underlying seasonal changes in mating and pregnancy, rather than as individual response variables of interest.
We developed generalized estimating equation (GEE) models using R v2.15.2 and the 'geepack' [64] to examine the correlations between pathogen types, immune parameters, and hormone concentrations in captured zebras. In all models we used a working correlation matrix with a first-order autoregressive relationship (AR-1) because, while individual immune, disease, and hormonal factors are likely correlated through time, these correlations should decrease between later time points and earlier samplings [65,66].
We developed each GEE by using a backwards stepwise refinement method based on comparing the quasi-likelihood under the independence model criterion (QIC) values between maximal models and models with variables removed [67]. The QIC is equivalent to Akaike's information criterion (AIC) for GEEs, which are not strictly likelihood based. After stepwise selection of the main terms, we added interaction terms between the remaining explanatory variables and further refined the models through backwards, stepwise selection. We subsequently validated the models using established methods [66].
NC Zebra and Springbok.
Similar to the GEE models built for captured zebra, we used each parasite type and FGM as response variables in their own models (Tables E and F in S1 File) to determine how seasonal factors, host immunity, hormone concentrations, and coinfections affected the intensity and prevalence of each pathogen, or how they affected stress hormone concentrations. In addition, because there were several missing variables for many NC zebra and springbok samples (see discussion above), we built a series of nested models for each response parameter. This allowed us to take advantage of the largest possible dataset for each iteration. We began with models of the response variable in question explained solely by simple host and environmental factors (rainfall, age, sex) and FGM (for parasite models) or by those factors and the GI parasite counts (for FGM models). This allowed us to use a sample size of 302 for NC zebra and greater than 260 for springbok (Tables E and F in S1 File). We then added the FPM and FEM reproductive hormones to each model (N = 163 complete cases for NC zebra; N = 113–114 for springbok). Finally, we built males-only models to account for the effects of FTM on parasite counts and FGM concentrations (N = 103 complete cases for NC zebra; N = 77 for springbok) (Tables E and F in S1 File).
We developed generalized linear models (GLM) using R v2.15.2 and the 'MASS' package [68] to examine the correlations between parasite types (Strongyles, Strongyloides and Coccidia for springbok; Strongyles for zebra), seasonal factors and hormone metabolite concentrations. We used negative binomial (NB) GLMs with a log link function for NC zebra GIP models and springbok GIP and GIC models to deal with large amounts of overdispersion in these counts variables [66] (Tables E and F in S1 File). We used a gamma distribution with an inverse link function for FGM GLMs for both NC zebra and springbok (Tables E and F in S1 File). For Strongyloides parasites (GIS) models, zero inflated negative binomial (ZINB) mixture models were used to deal with the large numbers of zero GIS counts (80 samples out of 273, 29.3%) and the large overdispersion present in these data [66] (Table F in S1 File); We used the 'pscl' package in R to develop these models [69]. We verified that ZINB models were more useful for modeling the GIS data than were regular negative binomial GLMs by applying Vuong tests to compare both types of models; a high positive Vuong statistic favors the ZINB over the NB GLM [70].
We developed each GLM or ZINB model by using a backwards, stepwise refinement method based on comparing the AIC values between maximal models and models with variables removed. After stepwise selection of the main terms, we added interaction terms between the remaining explanatory variables and further refined the models through backwards, stepwise selection. We validated the models using previously established methods [66]. For the final, best-fit ZINB models, we determined the fitted values for the logistic regression model components and determined the odds ratios of observing excess false zeros in conjunction with each predictive parameter [71].
Results
Seasonality of Pathogen Prevalence, Immune Measures, and Hormone Concentrations
Captured Zebras.
For both unique animals and paired recaptures, we found statistically significantly higher GI strongyle counts, total white blood cell counts, and neutrophil counts in the wet season (Tables G and H in S1 File; Fig. 1). Zebra experienced significantly higher estrogen concentrations in the dry season, indicating that the majority of mares were in mid-gestation in the middle of the hot dry season (Fig. 1).
Center lines in boxplots represent medians, with notches extending to ±1.58 IQR/ where IQR is the interquartile range and N is the sample size. Notches that do not overlap provide strong evidence that the medians differ (Chambers et al. 1983) (not shown for log2PA because notches grossly overlapped the ends of the boxes). The box hinges represent nearly the first and third quartiles, with whiskers showing the largest and smallest observations 1.5 times the box size from the nearest hinge. Points outside this range are shown as open circles. Boxplots in white are for variables that are significantly different from each other by t tests or Wilcoxon rank sum tests, whereas gray boxplots are not significantly different. sqrtHCT is light gray because the rain group for this variable was nearly significant (p = 0.080).
We found higher lymphocyte counts in the wet season, but only for paired individuals (Tables G and H in S1 File; Fig. 1); this indicates that controlling for individual variation allowed these subtler, yet present seasonal differences to be revealed. HCT was significantly and nearly significantly higher in the wet season for unique and paired animals, respectively. Stress hormone concentrations were also significantly higher in drier seasons compared to in the wetter ones for paired animals, though only prior to the Holm's p correction; this corroborates possibly higher stress in the dry season observed for unique animals.
Thus, overall, GI strongyle infection intensities, total white blood cell counts, lymphocyte and neutrophil cell counts, and HCT were higher in the wet season for these animals, whereas stress hormone concentrations were higher in the dry season while most mares were in mid-gestation.
Non-Captured Zebras.
Similar to the results for comparisons in captured zebra, we found that GI strongyle counts were significantly higher in wetter seasons that in drier ones. The highest overall GIP infection intensities were males and yearlings (Table I in S1 File).
Estrogen concentrations were significantly higher during the dry seasons (Fig. 2), and progesterone was significantly higher in wetter seasons (Fig. 3). Testosterone concentrations were not significantly different between seasons. Females had significantly higher stress hormone concentrations than did males, overall and for two age group comparisons (Table I in S1 File; Fig. 2). Adults in general had higher stress hormone levels.
For clarity, comparisons between interaction effects are not shown. Boxplot statistics are interpreted as in Fig. 1. Boxplots in white are for variables that are significantly different from each other by Tukey's t tests or Wilcoxon rank sum tests, whereas gray boxplots are not significantly different. Though non-transformed parasite counts (GIP) were compared with Wilcoxon rank sum tests, GIP is square root transformed here (as close to normality as possible) for clarity.
For clarity, comparisons between interaction effects are not shown. Boxplot statistics are interpreted as in Fig. 1. Boxplots in white are for variables that are significantly different from each other by Tukey's t tests or Wilcoxon rank sum tests, whereas gray boxplots are not significantly different. Though non-transformed FPM and FEM data were compared with Wilcoxon rank sum tests, they are fourth root and log transformed here (as close to normality as possible) for clarity.
Thus, overall, GI strongyle infection intensities were highest in the wet season and in yearling males, with males overall carrying more parasites than did females. Most mares were in mid-gestation in the dry season, and adult females had the highest stress hormone concentrations.
Springbok.
Similar to the results for comparisons in captured and NC zebra, we found that GI strongyle counts were significantly higher in wetter seasons. Strongyloides and Eimeria counts were also both significantly higher in springbok in the wet season compared to in the dry (Fig. 4). Yearlings harbored significantly more Strongyloides and Eimeria parasites than did adults. Male yearlings had significantly higher strongyle counts than did male and female adults, and possibly more than did female yearlings.
For clarity, comparisons between interaction effects are not shown. Boxplot statistics are interpreted as in Fig. 1. Boxplots in white are for variables that are significantly different from each other by Tukey's t tests or Wilcoxon rank sum tests, whereas gray boxplots are not significantly different. Though non-transformed parasite counts for all parasites were compared with Wilcoxon rank sum tests, they are fourth root and square root transformed here (as close to normality as possible) for clarity.
Estrogen concentrations were significantly higher during the dry seasons (Table J in S1 File; Figs. 5 and 6). Similar to the trends measured in NC zebra, springbok progesterone concentrations were significantly higher in the wet season (Fig. 5).
For clarity, comparisons between interaction effects are not shown. Boxplot statistics are interpreted as in Fig. 1. Boxplots in white are for variables that are significantly different from each other by Tukey's t tests or Wilcoxon rank sum tests, whereas gray boxplots are not significantly different. Though non-transformed FGM data were compared with Wilcoxon rank sum tests, they are fourth-root transformed here (to help meet assumptions of normality) for clarity.
FGM (a measure of host chronic stress responses) is used here as a proxy of host health, and is measured against the range of various GI parasite types found in these hosts. A) Zebra and strongyle parasites; B) Springbok and strongyle parasites; C) Springbok and Strongyloides parasites; D) Springbok and coccidian parasites. With a lack of host tolerance to these parasites, one would expect to see an increase in FGM (decrease in host health) with increased parasite infection intensities (as indicated by increasing epg—eggs per gram of feces—counts) [101] Instead, these hosts experience either no increase in FGM, or a decrease in FGM with increased parasite infection intensities in all instances.
Thus, overall, GI parasite infection intensities of all three kinds of parasites were higher in the wet season, with male yearlings being the highest parasite aggregators. The majority of adult ewes were in mid-gestation during the dry season, similar to zebra mares.
Interrelationships Between Pathogens, Physiology, and Season
Captured Zebras.
Pathogen Interactions: Higher cumulative rainfall two months prior to sampling (hereafter simply "rainfall") predicted higher GI parasite loads (GIP), as expected from rain group analyses (Fig. 1). Mares in full milk were more likely to harbor more strongyles than were those lactating water or not lactating at all. As FEM increased, PA titer prevalence decreased. Strongyle counts also negatively predicted anti-PA titers, as we have previously demonstrated [18] (Table 1). Thus, overall, GI strongyles were highest in the wet season and more intense infections were associated with lactation, and estrogen and strongyle infections were negatively associated with host ability to mount an anti-anthrax immune response.
Immune Interactions: Higher rainfall significantly predicted higher neutrophil counts (Table 2), corroborating similar trends we observed with our group comparisons (Fig. 1). Lactating mares were more likely to have higher neutrophil counts, whereas those with higher estrogen levels were negatively associated with neutrophil counts (Table 2). Contrary to what we observed in our group comparisons, increased rainfall in the GEE model predicted that hosts would carry fewer lymphocytes (Table 2). Neutrophils and strongyles positively predicted lymphocyte counts, similar to the relationships we observed in our neutrophil and GIP models. Pregnant mares had significantly fewer lymphocytes. Thus, overall, zebra in the wet season had higher neutrophil counts but lower lymphocyte counts, estrogen and pregnancy were associated with lower specific immune cell counts, and animals with higher strongyle infection intensities had higher lymphocyte counts.
Hormone Interactions: Rainfall negatively predicted stress hormone levels in our GEE model (Table 2). Both higher rainfall and higher progesterone levels together, however, predicted higher FGM. These opposite effects between rainfall and its interaction with FPM may reflect gestational timing factors. In accordance with this idea, increasing FEM positively predicted FGM levels. Thus, overall, the wet season was predominantly associated with less stress in zebra hosts, though zebra in the periparturient period in the wet season experienced higher stress, as did those with the highest estrogen levels.
NC Zebras.
Pathogen Interactions: In accordance with all other zebra strongyle models and comparisons, more rainfall predicted higher GI parasite loads in all GIP models (Table 3). However, higher rainfall in conjunction with higher stress hormone concentrations predicted less intense strongyle infections in the GIP-hormones model. This relationship was perhaps driven primarily by the stress effects; higher stress levels also significantly, negatively predicted strongyle counts in the simplest GIP-FGM model, corroborating the interaction effects in the more complex GIP-hormones model. This highly significant, negative relationship between FGM and strongyle count also held for the GIP-males only model (Table 3). Males were more likely to have higher strongyle counts in the simplest GIP-FGM model, similar to what we observed in our group comparisons (Fig. 2). Thus, overall, these models indicate that strongyles in zebras were predominantly positively associated with rainfall and males, and were negatively associated with stress hormone concentrations and females.
Hormone Interactions: More rainfall predicted lower stress hormone concentrations in all models, in concordance with the patterns seen in captured zebra (Table 2). In all three FGM models, higher strongyle counts predicted higher stress. However, higher stress levels predicted lower strongyle counts in GIP models; this negative correlation was of greater magnitude than was the positive correlation between strongyles and stress in the FGM models. (Table 3). In addition, rainfall-GIP interaction terms were nearly and highly significant in negatively predicting stress hormone levels in both the FGM-GIP and the FGM-males only models, respectively. Older females significantly predicted higher stress hormone concentrations in the FGM-GIP model, and males with higher testosterone levels were associated with lower FGM in the FGM-males only model. These results, taken together with the fact that females had significantly higher FGM in group comparisons (Fig. 2), suggest that female zebras likely do experience higher stress hormone levels than do males. Thus, overall, stress hormone concentrations were higher in the dry season (as they were in the captured zebra group), and stress hormones were typically negatively associated with GI strongyle infection intensity and maleness.
Springbok.
Pathogen Interactions: In all GIP models, higher rainfall predicted more strongyles (Table 4). Strongyloides counts were inversely associated with rainfall in the final GIS-FGM model; however, increased rainfall also predicted fewer false zeros in Strongyloides counts, indicating a complex relationship between these parasites and rainfall (Table 4). A simple linear model of rainfall on transformed GIS was highly significant and indicated a positive relationship between rainfall and Strongyloides infection intensities (N = 266, F1,264 = 45.9, R2 = 0.15, p = 0.000, coeffiecient = 0.06). Eimeria counts were associated with high rainfall and younger animals. Thus, overall, macroparasite infection intensities in springbok were highest in the wet season.
Hormone Interactions: Higher rainfall significantly, positively predicted higher stress hormones in the FGM-parasites model (Table 5). Higher stress hormones nearly significantly (p = 0.052) predicted higher strongyle counts in the GIP-FGM model; however, strongyle loads significantly, negatively predicted FGM in FGM models, similar to the inverse relationship between GIP and FGM in most NC zebra GIP models (Table 3). Increased stress hormone concentrations significantly predicted fewer Strongyloides in the GIS-FGM model (Table 4), and higher rainfall in conjunction with higher GIS significantly, negatively predicted FGM (Table 5), corroborating trends observed in GIS models. Eimeria counts alone did not significantly influence FGM in any models, and stress and other hormones did not play a role in predicting Eimeria counts in Eimeria models (Table 5). While higher testosterone concentrations in males predicted lower stress, higher Eimeria counts in conjunction with higher testosterone levels predicted higher stress hormone concentrations. Finally, females with high estrogen levels significantly predicted higher stress hormone concentrations in the FGM-hormones model. Thus, overall, higher stress hormone concentrations were associated with the wet season and estrogen, and were negatively associated with strongyles and Strongyloides. Stress hormones were not associated with Eimeria infections except in males with both high testosterone levels and Eimeria infections.
Discussion
This study evaluated the seasonality of stress responses in conjunction with environmental, reproductive, parasitic, and immunologic factors. We aimed to determine how the strong seasonality of environmentally transmitted parasites such as helminths and coccidian parasites correlated with changes in host reproduction, immune function, and measures of chronic stress, as well as with the occurrence of coinfections. As we have previously found that wet season-driven GI strongyles are likely immunomodulatory in zebra hosts and therefore may influence susceptibility to anthrax in Etosha National Park [18], we wanted to determine if these pathogens also influence host susceptibility by invoking stress responses and stress-related immunosuppression. In addition, such correlations would provide evidence that these GI parasites likely affect host fecundity and survival through mechanisms other than direct pathology and resource use. Instead, we found that, while stress responses do exhibit seasonality in zebra, stress hormone levels largely peak in the dry season when parasite infection intensities are lowest. Stress responses in springbok were more complex, with stress alone positively related to rainfall in some models, but significantly, negatively related to rainfall in conjunction with helminth infection intensities, and unrelated to coccidian infections. Rather than by parasite prevalence or infection intensity, peak stress hormone levels are most likely driven by dry season environmental challenges and host reproductive status in both species examined.
Seasonal Patterns
As expected, we found strong seasonality in nearly all parameters examined. Previous studies in this system have found that GI parasite infection intensities in zebra and springbok peak in the wet season, and that prevalence of strongyles, Strongyloides, and Eimeria in springbok reaches nearly 90% during times of highest rainfall [13]. While the seasonality of GI macroparasites in ENP is largely determined by environmental conditions [45,72,73], transmission rates can change seasonally due to changes in host behaviors and the birth of large numbers of immune naïve hosts, and host susceptibility can change due to physiological factors such as reproductive hormonal fluxes and changing nutrition [74,75]. Thus, our examination of the seasonality of these other factors was also important for determining what may influence parasite infections, and how parasite infections may interact temporally with physiological changes.
Though technically nonseasonal breeders [23], both zebra and springbok in ENP have been observed to breed on a largely seasonal basis, with birth peaks in the wet season [11,13]. Estrogen in pregnant zebras peaks in mid-gestation before declining sharply in the last trimester [21,51]. Estrogen also reaches a sustained peak in mid-pregnancy in springbok, with progesterone peaking just prior to parturition [22,76]. Our findings that FEM levels were significantly, negatively associated with rainfall for both species therefore indicates that the majority of the zebra and springbok examined were in mid-gestation in the hot dry season (late September-early November), whereas the significantly increased FPM in the wet season indicates that the majority of females of both species were giving birth and/or coming into estrus at this time. Thus, this study provides definitive hormonal evidence of strong reproductive seasonality for both of these species in this system. Interestingly, both male zebra and springbok in the dry season had significantly higher FEM than did females in the wet season, indicating, unlike in other similar studies [51], that there is likely an environmental, seasonal component to reproductive hormone secretion as well.
We also found seasonal patterns in WBC counts, HCT, lymphocytes, and neutrophils, in corroboration with the seasonality of several other immune factors previously measured in these zebra [18]. Higher HCT in the wet season likely reflects a higher nutritional state, especially in light of the fact that dehydration often results in an elevated HCT [77]. Higher WBC in the wet season may also reflect this higher plane of nutrition, as higher WBCs in general often represent a higher investment in immunity [78].
Despite individuals experiencing significantly more intense macroparasite infections in the wet season, higher stress hormones were overall correlated with decreased rainfall, particularly for zebras. Some of this seasonality was likely driven by environmental factors such as decreased nutrition and water availability in the dry season. Stress hormones may have been more seasonally linked in zebra than in springbok due to the fact that zebra are much more water-dependent than the more arid-adapted springbok [79,80]; zebras may be able to persist at their relatively high population level in ENP due, in part, to the fact that man-made water holes provide accessible drinking water year-round [33]. Elevations of GCs in the face of altered nutrition could be adaptive to help mobilize energy stores, rather than indications of distress: there is evidence that predictable aversive conditions such as seasonal changes may cause animals to preemptively raise GC levels [2]. However, the interaction between FGM and other factors during the dry season indicates that animals also likely do experience adverse chronic stress responses at this time, more so than they do during the wet season. Our more complex models examining the interrelationships between variables helped to shed light on the potential sources of these patterns (discussed below).
Stress and Reproduction
Stress hormone levels were not only higher during the same season most females were in mid-gestation, but were also significantly correlated with reproductive hormone concentrations. Increased FEM significantly predicted higher FGM for zebra (Table 2), and increased FEM predicted higher FGM for female springbok (Table 5). Thus, it appears that pregnancy and its attendant hormone increases correlate with a significant stress response in these animals. The lack of strong correlations in this study between researcher-determined pregnancy status and FGM are likely due to the subjective nature of pregnancy assessment by visual inspection and light palpation. While there are many studies illustrating the effects of stress on reproductive hormone concentrations and reproductive success, this is one of the first studies (also see [81]) examining how pregnancy itself may correlate with chronic stress in wildlife. Progesterone in zebras peaks in the last 50 days prior to foaling, and declines sharply in the short periparturient period [51,53].
In models examining only female zebra (captured animals), FPM was positively associated with FGM, and high rainfall interacting with high FPM significantly predicted increased stress hormone levels (Table 2); this likely reflects females in the periparturient period and just prior to estrus, and the distinct possibility that these reproductive efforts cause stress in mares. Progesterone in springbok was not related to stress, perhaps because these animals tend to give birth earlier in the wet season and mate in the early dry season, while samples examined in this study were predominantly taken from the mid wet season. In accordance with these findings, female zebras overall had higher stress hormone levels than did males. Other wildlife studies have found similar sex-stress patterns and have suggested that, in addition to the stress of reproduction, higher GC secretion in females may be adaptive to promote longer-term survival and the energy demands of reproduction [26].
Reproduction, Pathogens, and Immunity
Gestation in this study correlated with significantly higher stress hormone concentrations in female animals. While such chronic stress often leads to immunosuppression, and while immunosuppression can lead to increased disease susceptibility, we did not find evidence for this kind of cascade in our study. Instead, we found that reproductive hormones themselves influenced immunity more directly, and that reproductive hormones and sex were important factors in predicting pathogen infection intensity. Non-pregnant zebras had higher lymphocyte counts, and lower FEM was significantly correlated with higher neutrophil counts. Breeding and immunity are both expensive and somewhat incompatible, as elevations of one often compromise the other [82]. However, despite the fact that lactation and estrus behaviors are also very energetically costly [83], both lactation and FPM significantly predicted higher neutrophil counts in zebras; this may reflect an inflammatory response to the potentially injurious activities against mares of suckling young and mating stallions.
Pregnant and lactating mares, and mares with higher FEM in the dry season also had higher tick infestations, despite the fact that ticks are significantly more abundant in the dry season and mares lactate most intensively in the wet. Zebras with higher FEM also had lower anti-PA titers, perhaps indicating reproductive hormone immunomodulatory effects. Lactating zebra and those with higher FEM levels in the wet season also had more intense strongyle infections, likely indicating a periparturient rise in parasitism. Such increases in parasite infections are common around the time of parturition due to energy trade-offs between reproduction and immunity, and the immunosuppressive actions of fluctuating reproductive hormones [84,85].
Despite the potential immunosuppressive effects of pregnancy, parturition, and lactation, however, male zebra had significantly higher strongyle counts, and strongyles were most aggregated in younger animals of both species, the age classes that had the lowest stress hormone levels. While springbok group comparisons for males and females were not significantly different for any parasite, male yearlings again shed significantly more strongyle eggs and Eimeria oocysts than did other sex-age groups. These sex aggregations, however, were not related to testosterone levels, as FTM was not included in any final male parasite models. Male-biased parasitism is relatively common in wildlife [86], which has often been attributed to the immunosuppressive effects of testosterone [87]. Our results, however, suggest that Rolff’s application of Bateman’s rule [88,89] (i.e. that females gain fitness through increasing longevity, whereas males gain fitness through increasing mating rates) may be more applicable in this system: he suggested that females are overall more immunocompetent and resistant to infections while males trade off immunity for quicker reproductive fitness, and that this can manifest itself as higher parasite infection intensities in male hosts. The increased parasite infection intensities in younger animals may be due to seasonal reproduction resulting in a pulse of new, immunologically naïve hosts in the wet season, and parasites that have synchronized their life cycles with those of their hosts to ensure that pastures are contaminated with infective stages when transmission is most likely [90]. This decrease in parasite loads with host age also provides evidence for acquired immunity against GI parasites in this system. We also found that more intense strongyle infections predicted higher lymphocyte counts in zebra, despite the fact that lymphocytes were negatively correlated with rainfall; this perhaps correlates with new strongyle infections triggering a memory lymphocyte response in adult hosts, leading to some anti-parasite protection [10,91,92].
Stress, Pathogens, and Immunity
While stress hormone levels were significantly, positively correlated with several reproductive parameters, increased FGM largely did not correlate with increased pathogen infection intensity either alone or in conjunction with reproductive hormones. In fact, decreased FGM significantly predicted increased strongyles in zebra and increased Strongyloides in springbok, and higher strongyle counts in springbok significantly predicted lower FGM levels (note there is no evidence that parasites can directly or indirectly downregulate FGM production; [30]). Strongyle infection intensity in zebras and Strongyloides infection intensity in springbok did predict higher FGM levels; this could indicate that being able to acutely mount a small stress response is actually beneficial for animals fighting against new parasite infections. However, the overall magnitude of the negative relationships between FGM and strongyles in zebras was much higher. Thus, while there is some indication here that helminth parasites may correlate with some host stress, there are clearly more complex interactions between parasites and the HPA axis.
This negative or complete lack of association between stress hormones and macroparasites in this study is likely not due to mechanisms of acclimation or decoupling of physical stress and the HPA axis. While acclimation to a stressor can occur when animals experience the same stressor repeatedly, causing the stressor to be perceived to be less noxious and thereby decreasing the GC response [93], acclimation to parasites as predictable stressors is not likely in ENP; pregnancy and mating are also clearly seasonally predictable events in these populations, and stress hormones correlated significantly with these life events. Alternatively, chronic stress can actually cause GCs to decrease with intrinsic control of adrenal function [94]. This is likely not the case here, however, as lowered GCs often represent a highly distressed situation, and zebra and springbok were still capable of secreting significant GCs in response to reproductive and environmental stressors in both seasons.
Overall, then, parasites in this system do not appear to affect host stress levels. This is even more apparent as male yearlings of both species harbored the most parasites of any group, but females of nearly all age-sex comparisons had higher stress levels in zebras. Male adult springbok had nearly significantly higher stress hormone concentrations than did male yearlings, despite yearlings harboring significantly higher parasite intensities. In addition, Turner et al. 2012 found that body condition of adult female springbok was significantly, negatively related to strongyle intensity, whereas that of yearlings and adult males was significantly, negatively related to rainfall. Therefore, the springbok age-sex class with the highest condition impact from parasites did not have the highest FGM levels, even in the face of concurrent parturition and lactation. In fact, the only relationship we found between stress, reproduction, and parasites all together was seen in male springbok experiencing both high Eimeria intensities and high testosterone concentrations. Thus, only through the interaction with significant amounts of testosterone were parasites correlated with a stress response. While there were significant macroparasite coinfection interactions, none of these interactions were significantly correlated with host stress levels. As increasing GI parasite species richness has been associated with poorer host body condition, the lack of association between coinfection richness in springbok and stress hormones supports the idea that parasites do not cause these hosts significant stress [9].
Given the evidence that GI parasites are capable of causing host pathology [12], immunomodulation [18], and immunosuppression, and that different types of parasites in numerous systems have been found to be strongly correlated with host stress responses in conjunction with parasite effects on host health [95,96], the fact that these parasites persist in ENP hosts with no relationships to chronic stress responses suggests that hosts are tolerant of their parasites (Fig. 6).
Hosts that are good at reducing or clearing parasite burdens are not necessarily the healthiest; eliminating parasites is immunologically costly, and hosts risk immunopathology with vigorous immune responses [97]. In addition, as immunity and reproductive efforts fight for the same resources, increased immune responses, particularly in the face of external resource limitations, can reduce host fecundity. Thus, hosts that can tolerate parasites may have advantages over those that cannot. Tolerance reflects the rate of decline in fitness as parasite burden increases: less fitness decline with increased parasite infection intensity indicates higher parasite tolerance [98]. Tolerance mechanisms are selected for when there are high rates of infection but relatively low virulence [99]. While GI parasites in ENP do affect host condition [13], and likely affect immune function and susceptibility to other pathogens [18], the fact that they are nearly ubiquitous in herbivore hosts suggests that their overall virulence is low (due to transmission-virulence tradeoffs [100]), particularly compared to other pathogens in the system such as B. anthracis. In addition, while the evolution of resistance should reduce parasite prevalence in host populations, the evolution of tolerance to parasites should have a neutral or positive effect on parasite prevalence [98], as is observed in herbivores in ENP. Thus, it appears that, while GI parasites may influence susceptibility to anthrax and ticks in ENP hosts through immunomodulation [18,31], these parasites do not appear to increase host coinfection susceptibility through stress immunosuppression.
Conclusions
In summary, we found evidence that stress, reproduction, immunity, and pathogens exhibit strong seasonality in ENP. Most zebra and springbok in ENP are indeed seasonal breeders, with both species experiencing parturition and lactation primarily in the wet season, zebra experiencing most estrus and breeding behaviors in the wet season, and both species experiencing mid-gestation in the hot dry season. Overall immune allocation is highest during the wet season, likely driven, at least in part, by nutritional factors [17]. Stress hormone levels are highest overall in the dry season, though these seasonal signals are stronger in zebras. While peaks in reproductive hormone levels are correlated with immunosuppression, increased parasite infection intensity, and decreased anti-anthrax antibody titers, stress hormone levels are largely unaffected by macroparasite prevalence or infection intensity. Instead, stress hormone levels are most strongly correlated with environmental and reproductive factors, particularly with the demands of mid-gestation. This likely indicates that these hosts are evolutionarily tolerant of their GI parasites, despite the fact that they are not neutral actors in these hosts; parasites do affect body condition and use host resources, hosts do allocate immune resources toward keeping parasite infections in check, and parasites likely cause host immunomodulation and may even increase susceptibility to B. anthracis in this system. However, parasite tolerance strategies indeed make sense in a system in which macroparasites are nearly ubiquitous but almost never lethal, while pathogens such as B. anthracis are often lethal but considerably more rare.
Supporting Information
S1 File. Supplemental Information.
This file contains supplementary methods for laboratory protocols and statistical analyses, a map of the study region, and tables with information about animal captures, variables and abbreviations, maximal GEE and GLM models, and statistical results for seasonal and sex and age group comparisons.
https://doi.org/10.1371/journal.pone.0120800.s001
(ZIP)
Acknowledgments
We thank the Namibian Ministry of Environment and Tourism for permission to do this research, the Directorate of Parks, Wildlife and Management for permission to work throughout Etosha, and the staff in the Directorate of Scientific Services at the Etosha Ecological Institute for logistical support and assistance. We would like to give special thanks to veterinarians Mark Jago, Conrad Brain, Peter Morkel, and Ortwin Aschenborn for their assistance with animal captures, as well as to Martina Küsters, Shayne Kötting, Gabriel Shatumbu, Wilferd Versfeld, Marthin Kasaona, Royi Zidon, and Werner Kilian, among others, for their tremendous help in the field. We thank Steve Bellan for assistance with the anti-PA ELISA endpoint titer determinations, and Kevin Burgio for editing suggestions.
Author Contributions
Conceived and designed the experiments: CAC NP. Performed the experiments: CAC WCT NP. Analyzed the data: CAC. Contributed reagents/materials/analysis tools: NP WMG. Wrote the paper: CAC WCT WMG.
References
- 1. Dhabhar FS, McEwen BS (1997) Acute stress enhances while chronic stress suppresses cell-mediated immunity in vivo: a potential role for leukocyte trafficking. Brain Behav Immun 11: 286–306. pmid:9512816
- 2. Martin LB (2009) Stress and immunity in wild vertebrates: timing is everything. Gen Comp Endocrinol 163: 70–76. pmid:19318107
- 3. Cohen S, Janicki-Deverts D, Doyle WJ, Miller GE, Frank E, Rabin BS, et al. (2012) Chronic stress, glucocorticoid receptor resistance, inflammation, and disease risk. Proc Natl Acad Sci 109: 5995–5999. pmid:22474371
- 4. Dobson H, Smith RF (2000) What is stress, and how does it affect reproduction? Anim Reprod Sci 60–61: 743–752.
- 5. Khansari D, Murgo AJ, Faith RE (1999) Effects of stress on the immune system. Immunol Today 11: 170–175.
- 6. Muehlenbein MP (2006) Intestinal parasite infections and fecal steroid levels in wild chimpanzees. Am J Phys Anthropol 130: 546–550. pmid:16444733
- 7. Chapman CA, Saj TL, Snaith TV (2007) Temporal dynamics of nutrition, parasitism, and stress in colobus monkeys: implications for population regulation and conservation. Am J Phys Anthropol 134: 240–250. pmid:17596853
- 8. Gulland FM (1992) The role of nematode parasites in Soay sheep (Ovis aries L.) mortality during a population crash. Parasitology 105: 493–503. pmid:1461688
- 9. Lello J, Boag B, Hudson PJ (2005) The effect of single and concomitant pathogen infections on condition and fecundity of the wild rabbit (Oryctolagus cuniculus). Int J Parasitol 35: 1509–1515. pmid:16051247
- 10. Yun CH, Lillehoj HS, Lillehoj EP (2000) Intestinal immune responses to coccidiosis. Dev Comp Immunol 24: 303–324. pmid:10717295
- 11. Gasaway WC, Gasaway KT, Berry HH (1996) Persistent low densities of plains ungulates in Etosha National Park, Namibia: testing the food-regulating hypothesis. Can J Zool 74: 1556–1572.
- 12. Turner WC, Versfeld WD, Kilian JW, Getz WM (2012) Synergistic effects of seasonal rainfall, parasites and demography on fluctuations in springbok body condition. J Anim Ecol 81: 58–69. pmid:21831195
- 13. Turner WC, Getz WM (2010) Seasonal and demographic factors influencing gastrointestinal parasitism in ungulates of Etosha National Park. J Wildl Dis 46: 1108–1119. pmid:20966262
- 14. Lindeque PM, Turnbull PC (1994) Ecology and epidemiology of anthrax in the Etosha National Park, Namibia. Onderstepoort J Vet Res 61: 71–83. pmid:7898901
- 15. Bellan SE, Gimenez O, Choquet R, Getz WM (2013) A hierarchical distance sampling approach to estimating mortality rates from opportunistic carcass surveillance data. Methods Ecol Evol 4: 361–369.
- 16.
Turner WC (2009) The ecology of orally ingested parasites in ungulates of Etosha National Park University of California, Berkeley.
- 17. Turner WC, Imologhome P, Havarua Z, Kaaya GP, Mfune JKE, Mpofu IDT, et al. (2013) Soil ingestion, nutrition and the seasonality of anthrax in herbivores of Etosha National Park. Ecosphere 4: 1–19.
- 18. Cizauskas CA, Turner WC, Wagner B, Küsters M, Vance RE, Getz WM (2014) Gastrointestinal helminths affect host susceptibility to anthrax through seasonal immune trade-offs. BMC Ecol 14: 1–15. pmid:24438134
- 19. Saltz D, White GC (1991) Urinary cortisol and urea nitrogen responses to winter stress in mule deer. J Wildl Manage 55: 1–16.
- 20. DelGiudice GD, Mech LD, Kunkel KE, Gese EM, Seal ES (1992) Seasonal patterns of weight, hematology, and serum characteristics of free-ranging female white-tailed deer in Minnesota. Can J Zool 70: 974–983.
- 21. Barnier F, Grange S, Ganswindt A, Ncube H, Duncan P (2012) Inter-birth interval in zebras is longer following the birth of male foals than after female foals. Acta Oecologica 42: 11–15.
- 22. Skinner DC, Richter TA, Malpaux B, Skinner JD (2001) Annual ovarian cycles in an aseasonal breeder, the springbok (Antidorcas marsupialis). Biol Reprod 64: 1176–1182. pmid:11259265
- 23.
Estes RD (1991) The behavior guide to African mammals; including hoofed mammals, carnivore, primates. Berkeley, California: University of California Press.
- 24. Nagy KA, Knight MH (1994) Energy, water, and food use by springbok antelope (Antidorcas marsupialis) in the Kalahari Desert. J Mammal 75: 860–872.
- 25.
Robbins CT (1993) Wildlife feeding and nutrition. 2nd ed. New York, NY: Academic Press.
- 26.
Vera F, Zenuto RR, Antenucci CD (2013) Seasonal variations in plasma cortisol, testosterone, progesterone and leukocyte profiles in a wild population of tuco-tucos. J Zool 289: 111–118. https://doi.org/10.1111/j.1469-7998.2012.00967.x
- 27. Martin LB, Weil ZM, Nelson RJ (2008) Seasonal changes in vertebrate immune activity: mediation by physiological trade-offs. Philos Trans R Soc B Biol Sci 363: 321–339. pmid:17638690
- 28. Archie EA, Altmann J, Alberts SC (2012) Social status predicts wound healing in wild baboons. Proc Natl Acad Sci 109: 9017–9022. pmid:22615389
- 29. Romero LM (2002) Seasonal changes in plasma glucocorticoid concentrations in free-living vertebrates. Gen Comp Endocrinol 128: 1–24. pmid:12270784
- 30. Monello RJ, Millspaugh JJ, Woods RJ, Gompper ME (2010) The influence of parasites on faecal glucocorticoid metabolite levels in raccoons: an experimental assessment in a natural setting. J Zool 282: 100–108.
- 31. Kamath PL, Turner WC, Küsters M, Getz WM (2014) Parasite-mediated selection drives an immunogenetic trade-off in plains zebras (Equus quagga). Proc R Soc B—Biol Sci 281: 1–10.
- 32. Engert S (1997) Spatial variability and temporal periodicity of rainfall in the Etosha National Park and surrounding areas in northern Namibia. Madoqua 20: 115–120.
- 33. Auer C (1997) Chemical quality of water at waterholes in the Etosha National Park. Madoqua 20: 121–128.
- 34. Smuts GL (1974) Age determination in Burchell’s zebra (Equus burchelli antiquorum) from the Kruger National Park. J South African Wildl Manag Assoc 4: 103–115.
- 35. Mantovani A, Cassatella MA., Costantini C, Jaillon S (2011) Neutrophils in the activation and regulation of innate and adaptive immunity. Nat Rev Immunol 11: 519–531. pmid:21785456
- 36. Viney ME, Riley EM, Buchanan KL (2005) Optimal immune responses: immunocompetence revisited. Trends Ecol Evol 20: 665–669. pmid:16701455
- 37. Little SF, Ivins BE, Fellows PF, Friedlander AM (1997) Passive protection by polyclonal antibodies against Bacillus anthracis infection in guinea pigs. Infect Immun 65: 5171–5175. pmid:9393812
- 38. Marcus H, Danieli R, Epstein E, Velan B, Shafferman A, Reuveny S (2004) Contribution of immunological memory to protective immunity conferred by Bacillus anthracis protective antigen-based vaccine. Infect Immun 72: 3471–3477. pmid:15155654
- 39. Pickering AK, Merkel TJ (2004) Macrophages release tumor necrosis factor alpha and interleukin-12 in response to intracellular Bacillus anthracis spores. Infect Immun 72: 3069–3072. pmid:15102824
- 40. Cizauskas CA, Bellan SE, Turner WC, Vance RE, Getz WM (2014) Frequent and seasonally variable sublethal anthrax infections are accompanied by short-lived immunity in an endemic system. J Anim Ecol 83: 1078–1090.
- 41. Villanúa D, Pérez-Rodríguez L, Gortázar C, Höfle U, Viñuela J (2006) Avoiding bias in parasite excretion estimates: the effect of sampling time and type of faeces. Parasitology 133: 251–259. pmid:16707034
- 42. Horak IG, Anthonissen M, Krecek RC, Boomker J (1992) Arthropod parasites of springbok, gemsbok, kudus, giraffes and Burchell’s and Hartmann's zebras in the Etosha and Hardap Nature Reserves, Namibia. Onderstepoort J Vet Res 59: 253–257. pmid:1297955
- 43. Biggs HC, Langenhoven JW (1984) Seasonal prevalence of Ixodid ticks on cattle in the Windhoek district of South West Africa/Namibia. Onderstepoort J Vet Res 51: 175–182. pmid:6533509
- 44. Nielsen MK, Kaplan RM, Thamsborg SM, Monrad J, Olsen SN (2007) Climatic influences on development and survival of free-living stages of equine strongyles: implications for worm control strategies and managing anthelmintic resistance. Vet J 174: 23–32. pmid:16815051
- 45. Martinaud G, Billaudelle M, Moreau J (2009) Circadian variation in shedding of the oocysts of Isospora turdi (Apicomplexa) in blackbirds (Turdusmerula): an adaptative trait against desiccation and ultraviolet radiation. Int J Parasitol 39: 735–739. pmid:19100268
- 46.
Gibbons LM, Jacobs DE, Fox MT, Hansen J (2005) McMaster egg counting technique. R Vet Coll Agric Organ United Nations Guid to Vet Diagnostic Parasitol. Available: http://www.rvc.ac.uk/Review/Parasitology/EggCount/Purpose.htm.
- 47.
Bowman DD (2003) Georgis’ parasitology for veterinarians, 8th edition. 8th editio. Philadelphia, PA: W.B. Saunders.
- 48. Keay JM, Singh J, Gaunt MC, Kaur T (2006) Fecal glucocorticoids and their metabolites as indicators of stress in various mammalian species: a literature review. J Zoo Wildl Med 37: 234–244. pmid:17319120
- 49. Millspaugh JJ, Washburn BE (2004) Use of fecal glucocorticoid metabolite measures in conservation biology research: considerations for application and interpretation. Gen Comp Endocrinol 138: 189–199. pmid:15364201
- 50. Monfort SL, Wasser SK, Mashburn KL, Burke M, Brewer BA, Creel SR (1997) Steroid metabolism and validation of noninvasive endocrine monitoring in the African wild dog (Lycaon pictus). Zoo Biol 16: 533–548.
- 51. Ncube H, Duncan P, Grange S, Cameron EZ, Barnier F, Ganswindt A (2011) Pattern of faecal 20-oxopregnane and oestrogen concentrations during pregnancy in wild plains zebra mares. Gen Comp Endocrinol 172: 358–362. pmid:21463629
- 52. Möstl E, Maggs JL, Schrötter G, Besenfelder U, Palme R (2002) Measurement of cortisol metabolites in faeces of ruminants. Vet Res Commun 26: 127–139. pmid:11922482
- 53. Asa CS, Bauman JE, Houston EW, Fischer MT, Read B, Brownfield CM, et al. (2001) Patterns of excretion of fecal estradiol and progesterone and urinary chorionic gonadotropin in Grevy’s zebras (Equus grevyi): ovulatory cycles and pregnancy. Zoo Biol 20: 185–195.
- 54. Palme R, Möstl E (1997) Measurement of cortisol metabolites in the faeces of sheep as a parameter of cortisol concentration in blood. Int J Mamm Biol 62: 192–197.
- 55. Little RJ, Rubin DB (2000) Causal effects in clinical and epidemiological studies via potential outcomes: Concepts and analytical approaches. Annu Rev Public Health 21: 121–145. pmid:10884949
- 56. Van Buuren S, Groothuis-Oudshoorn K (2011) mice: Multivariate imputation by chained equations in R. J Stat Softw 45: 1–67.
- 57. Rubin DB (1996) Multiple imputation after 18+ years. J Am Stat Assoc 91: 473–489.
- 58. Van der Heijden GJMG, Donders a RT, Stijnen T, Moons KGM (2006) Imputation of missing values is superior to complete case analysis and the missing-indicator method in multivariable diagnostic research: a clinical example. J Clin Epidemiol 59: 1102–1109. pmid:16980151
- 59.
Team R (2013) R Development Core Team. R A Lang Environ Stat Comput.
- 60. Vergouw D, Heymans MW, van der Windt DAWM, Foster NE, Dunn KM, van der Horst HE, et al. (2012) Missing data and imputation: a practical illustration in a prognostic study on low back pain. J Manipulative Physiol Ther 35: 464–471. pmid:22964020
- 61.
Rubin DB (1987) Multiple imputation for nonresponse in surveys. New York, NY: John Wiley & Sons, Ltd.
- 62. Holm S (1979) A simple sequentially rejective multiple test procedure. Scand J Stat 6: 65–70.
- 63.
Fox J, Weisberg S (2011) An R companion to applied regression. Second edi. Thousand Oaks, California: Sage.
- 64. Højsgaard S, Halekoh U, Yan J (2006) The R package geepack for generalized estimating equations. J Stat Softw 15: 1–11.
- 65. Liang K-Y, Zeger SL (1986) Longitudinal data analysis using generalized linear models. Biometrika 73: 13–22.
- 66.
Zuur A, Ieno E, Walker N, Saveliev A, Smith G (2009) Mixed Effects Models and Extensions in Ecology with R. New York, NY: Springer Science+Business Media, LLC.
- 67. Pan W (2001) Akaike’s information criterion in generalized estimating equations. Biometrics 57: 120–125. pmid:11252586
- 68.
Venables WN, Ripley BD (2002) Modern applied statistics with S. Fourth edi. New York, NY: Springer.
- 69. Zeileis A, Kleiber C, Jackman S (2008) Regression models for count data in R. J Stat Softw 27.
- 70. Vuong QH (1989) Likelihood ratio tests for model selection and non-nested hypotheses. Econometrica 57: 307–333.
- 71. Pilosof S, Dick CW, Korine C, Patterson BD, Krasnov BR (2012) Effects of anthropogenic disturbance and climate on patterns of bat fly parasitism. PLoS One 7: e41487. pmid:22829953
- 72. Cattadori IM, Boag B, Bjørnstad ON, Cornell SJ, Hudson PJ (2005) Peak shift and epidemiology in a seasonal host-nematode system. Proc R Soc B—Biol Sci 272: 1163–1169.
- 73. Banks DJ, Singh R, Barger IA, Pratap B, le Jambre LF (1990) Development and survival of infective larvae of Haemonchus contortus and Trichostrongylus colubriformis on pasture in a tropical environment. Int J Parasitol 20: 155–160. pmid:2332275
- 74. Grassly NC, Fraser C (2006) Seasonal infectious disease epidemiology. Proc R Soc B 273: 2541–2550. pmid:16959647
- 75. Cornell SJ, Bjornstad ON, Cattadori IM, Boag B, Hudson PJ (2008) Seasonality, cohort-dependence and the development of immunity in a natural host-nematode system. Proc R Soc B—Biol Sci 275: 511–518.
- 76. Skinner DC, Cilliers SD, Skinner JD (2002) Effect of ram introduction on the oestrous cycle of springbok ewes (Antidorcas marsupialis). Reproduction 124: 509–513. pmid:12361469
- 77. Voyles J, Vredenburg VT, Tunstall TS, Parker JM, Briggs CJ, Rosenblum EB (2012) Pathophysiology in mountain yellow-legged frogs (Rana muscosa) during a Chytridiomycosis outbreak. PLoS One 7: e35374. pmid:22558145
- 78. Nunn CL, Gittleman JL, Antonovics J (2003) A comparative study of white blood cell counts and disease risk in carnivores. Proc R Soc B—Biol Sci 270: 347–356.
- 79. Nagy K (1994) Seasonal water, energy and food use by free-living, arid-habitat mammals. Aust J Zool 42: 55.
- 80. Cain JW, Owen-Smith N, Macandza VA (2011) The costs of drinking: comparative water dependency of sable antelope and zebra. J Zool 286: 58–67. pmid:21619910
- 81. Foerster S, Cords M, Monfort SL (2012) Seasonal energetic stress in a tropical forest primate: proximate causes and evolutionary implications. PLoS One 7: e50108. pmid:23209651
- 82. Bonneaud C, Mazuc J, Gonzalez G, Haussy C, Hastel O, Faivre B, et al. (2003) Assessing the cost of mounting an immune response. Am Nat 161: 367–379. pmid:12703483
- 83. Olléová M, Pluháček J, King S (2012) Effect of social system on allosuckling and adoption in zebras. J Zool 288: 127–134.
- 84. Coop RL, Kyriazakis I (1999) Nutrition–parasite interaction. Vet Parasitol 84: 187–204. pmid:10456415
- 85. Beasley AM, Kahn LP, Windon RG (2010) The periparturient relaxation of immunity in Merino ewes infected with Trichostrongylus colubriformis. Vet Parasitol 168: 60–70. pmid:19775817
- 86. Grear DA, Perkins SE, Hudson PJ (2009) Does elevated testosterone result in increased exposure and transmission of parasites? Ecol Lett 12: 528–537. pmid:19392718
- 87. Ezenwa VO, Stefan Ekernas L, Creel S (2012) Unravelling complex associations between testosterone and parasite infection in the wild. Funct Ecol 26: 123–133.
- 88. Bateman AJ (1948) Intra-sexual selection in Drosophila. Heredity (Edinb) 2: 349–368. pmid:18103134
- 89. Rolff J (2002) Bateman’s principle and immunity. Proc R Soc B—Biol Sci 269: 867–872.
- 90. Altizer S, Dobson A, Hosseini P, Hudson P, Pascual M, Rohani P (2006) Seasonality and the dynamics of infectious diseases. Ecol Lett 9: 467–484. pmid:16623732
- 91. Rowe A, Mcmaster K, Emery D, Sangster N (2008) Haemonchus contortus infection in sheep: parasite fecundity correlates with worm size and host lymphocyte counts. Vet Parasitol 153: 285–293. pmid:18337014
- 92. Daugschies A, Najdrowski M (2005) Eimeriosis in cattle: current understanding. J Vet Med B 52: 417–427. pmid:16364016
- 93. Romero LM (2004) Physiological stress in ecology: lessons from biomedical research. Trends Ecol Evol 19: 249–255. pmid:16701264
- 94. Linklater WL, MacDonald EA, Flamand JRB, Czekala NM (2010) Declining and low fecal corticoids are associated with distress, not acclimation to stress, during the translocation of African rhinoceros. Anim Conserv 13: 104–111.
- 95. Mougeot F, Martínez-Padilla J, Bortolotti GR, Webster LMI, Piertney SB (2010) Physiological stress links parasites to carotenoid-based colour signals. J Evol Biol 23: 643–650. pmid:20074170
- 96. Martínez-de la Puente J, Merino S, Tomás G, Moreno J, Morales J, Lobato E, et al. (2011) Nest ectoparasites increase physiological stress in breeding birds: an experiment. Naturwissenschaften 98: 99–106. pmid:21120448
- 97. Allen JE, Maizels RM (2011) Diversity and dialogue in immunity to helminths. Nat Rev Immunol 11: 375–388. pmid:21610741
- 98. Råberg L, Graham AL, Read AF (2009) Decomposing health: tolerance and resistance to parasites in animals. Philos Trans R Soc B—Biol Sci 364: 37–49.
- 99. Restif O, Koella JC (2004) Concurrent evolution of resistance and tolerance to pathogens. Am Nat 164: E90–E102. pmid:15459887
- 100. Leggett HC, Buckling A, Long GH, Boots M (2013) Generalism and the evolution of parasite virulence. Trends Ecol Evol 28: 592–596. pmid:23968968
- 101. Schneider DS, Ayres JS (2008) Two ways to survive infection: what resistance and tolerance can teach us about infectious diseases. Nat Rev Immunol 8: 889–895. pmid:18927577