Life History and Demographic Drivers of Reservoir Competence for Three Tick-Borne Zoonotic Pathogens

Abstract

Animal and plant species differ dramatically in their quality as hosts for multi-host pathogens, but the causes of this variation are poorly understood. A group of small mammals, including small rodents and shrews, are among the most competent natural reservoirs for three tick-borne zoonotic pathogens, Borrelia burgdorferi, Babesia microti, and Anaplasma phagocytophilum, in eastern North America. For a group of nine commonly-infected mammals spanning >2 orders of magnitude in body mass, we asked whether life history features or surrogates for (unknown) encounter rates with ticks, predicted reservoir competence for each pathogen. Life history features associated with a fast pace of life generally were positively correlated with reservoir competence. However, a model comparison approach revealed that host population density, as a proxy for encounter rates between hosts and pathogens, generally received more support than did life history features. The specific life history features and the importance of host population density differed somewhat between the different pathogens. We interpret these results as supporting two alternative but non-exclusive hypotheses for why ecologically widespread, synanthropic species are often the most competent reservoirs for multi-host pathogens. First, multi-host pathogens might adapt to those hosts they are most likely to experience, which are likely to be the most abundant and/or frequently bitten by tick vectors. Second, species with fast life histories might allocate less to certain immune defenses, which could increase their reservoir competence. Results suggest that of the host species that might potentially be exposed, those with comparatively high population densities, small bodies, and fast pace of life will often be keystone reservoirs that should be targeted for surveillance or management.

About the Authors:

Richard S. Ostfeld

* E-mail: [email protected]

Affiliation: Cary Institute of Ecosystem Studies, Millbrook, New York, United States of America

Taal Levi

Affiliation: Department of Fisheries and Wildlife, Oregon State University, Corvallis, Oregon, United States of America

Anna E. Jolles

Affiliation: Departments of Biomedical Sciences and Zoology, Oregon State University, Corvallis, Oregon, United States of America

Lynn B. Martin

Affiliation: Department of Integrative Biology, University of South Florida, Tampa, Florida, United States of America

Parviez R. Hosseini

Affiliation: EcoHealth Alliance, New York, New York, United States of America

Felicia Keesing

Affiliation: Biology Program, Bard College, Annandale-on-Hudson, New York, United States of America

Introduction

Blacklegged ticks (Ixodes scapularis) are the primary vectors responsible for transmitting a suite of zoonotic pathogens in eastern and central North America. Lyme disease, caused by the spirochete bacterium, Borrelia burgdorferi, is the most prevalent of these tick-borne diseases, with approximately 300,000 cases occurring annually in the United States [1]. Two other emerging tick-borne diseases are human babesiosis, caused by the piroplasm, Babesia microti, and human granulocytic anaplasmosis, caused by the rickettsial bacterium, Anaplasma phagocytophilum. Despite substantial phylogenetic separation, these three pathogens co-circulate in host and vector communities and show similar transmission dynamics [2], the causes of which we address here.

All three of these pathogens have a broad host range, routinely infecting most species in the community of terrestrial mammals and birds that have been tested [3]–[6]. For all three pathogens, ticks acquire infection when they take a blood meal from an infected host; transovarial transmission has not been demonstrated [7]–[10]. Therefore, considerable interest has focused on the capacity of different vertebrate host species to transmit infection to feeding ticks. The probability that a particular host species will transmit an infection to uninfected ticks under natural conditions has been termed “realized reservoir competence”. Realized reservoir competence is the product of the probability that a host is infected (prevalence) and the probability that it transmits given that it is infected (infectivity) [11]. When prevalence in a given host species is high, which is the case in many endemic zones, its realized reservoir competence and infectivity converge [11]. In regards to the above pathogens, realized reservoir competence is estimated by collecting fed larval blacklegged ticks after they have detached from a known host, allowing them to molt to the nymphal stage, and then testing nymphs individually for infection with each of the zoonotic agents (e.g., [3]). Any nymphs that test positive must have acquired the infection from the host from which they were collected.

For each of the three above pathogens, realized reservoir competence varies dramatically among vertebrate species [3]–[6]. But a striking pattern has emerged from recent studies in the hyperendemic zone for tick-borne zoonoses in the northeastern United States. Two rodents, the white-footed mouse (Peromyscus leucopus) and the eastern chipmunk (Tamias striatus), and two soricimorphs, the masked shrew (Sorex cinereus) and the northern short-tailed shrew (Blarina brevicauda), have among the highest realized reservoir competence values of all hosts for all three pathogens [3]–[6]. Here, we ask whether these hosts share attributes that facilitate the success of this group of very different zoonotic agents.

Life history theory suggests that species of animals and plants that invest heavily in growth and reproduction are subject to trade-offs that reduce allocation to self-maintenance, including immune defense [12]. Consequently, species with fast life histories (early maturity, rapid breeding, short longevity) might generally be more likely to become infected and be infectious than are their slower-lived counterparts, all else equal [13]–[20]. If investments in fast-paced lives indeed reduce allocation to immune defense, then relatively fast-lived species might be more competent reservoirs for multi-host pathogens [20]. Support for this expectation recently was found for the Lyme disease system by Previtali et al. [21], who examined some immune functions in white-footed mice, eastern chipmunks, and eastern gray squirrels (Sciurus carolinensis); these hosts show, respectively, high, moderate, and low reservoir competence for B. burgdorferi. Antibody responses to a novel antigen challenge were highest in squirrels, which have the slowest life history, and lowest in mice, which have the fastest life history [21]. Huang et al. [22] tested the hypothesis that the most competent reservoirs for the agent of Lyme disease are the species with the fastest life histories; for the nine host species they examined, reservoir competence was inversely related to body size. They found similar patterns for two mosquito-transmitted zoonotic pathogens. However, whether body size is somehow causally related to reservoir competence, or simply a correlate of other, more direct influences, is unknown.

In addition to the potential for host life history traits to influence reservoir competence via immune functioning, host traits might also influence the ability of pathogens to adapt to and exploit hosts. Because an extreme host generalist vectors them, all three of the tick-borne pathogens are distributed each generation to many species of hosts. Adaptation to evade immunity and proliferate in one host might impede the ability of the pathogen to infect and proliferate in other hosts [23]. Ostfeld and Keesing [24] proposed that pathogens transmitted by highly non-specific vectors might adapt to the hosts they are most likely to experience. This process would be akin to adaptation by passively-dispersed, free-living organisms to the habitat type they are most likely to encounter. They hypothesized that hosts that are widespread and abundant might be encountered by multi-host pathogens more frequently, leading to pathogen adaptation to exploit most common hosts. Such adaptation would produce high reservoir competence. Clearly though, the potential exists for both host life history traits and pathogen adaptation to affect variation in reservoir competence among host species. In this study we used a model selection approach to ask whether reservoir competence for three emerging tick-borne pathogens is more strongly correlated to: (1) life-history traits (i.e. fast pace of life predicts high competency), or (2) ecological traits (i.e., probability of encounter with infected ticks predicts high competency).

Materials and Methods

Life history traits

This analysis was conducted to ask whether reservoir competence for the set of tick-borne, zoonotic pathogens is predictable based on life-history features of the host, under the assumption that life history features influence defense against these pathogens (see Introduction). We began with the set of life history traits selected by Huang et al. [22], namely, adult body mass, gestation length, and litter size, for each of the nine mammalian hosts for which we have gathered data on reservoir competence (see below). We note that body size is, strictly speaking, not a life history trait, although it impacts many life history variables that influence the pace of life [25]. Given little knowledge on which to base hypotheses linking specific life history traits to specific immune traits, we take an inclusive approach here. Following Huang et al. [22], we log transformed gestation length and body mass. We conducted a second, more comprehensive analyses of pace of life impacts on competence that included adult body mass, gestation length, offspring produced per year, age at maturity, longevity, and log (basal metabolic rate * g⁻¹) (Table S1). Life history data were taken from the PanTHERIA database of mammal life history [26]. For eastern chipmunks and white-tailed deer, data on basal metabolic rate were not available in PanTHERIA, so we obtained relevant data from AnAge, the animal aging and longevity database (http://genomics.senescence.info/species/). Because many life history traits are correlated, we used principal components analysis to derive a set of new orthogonal variables that maximized the variance explained. We then conducted linear regressions with the first and second principal components as independent variables and realized reservoir competence as the dependent variable. We performed these analyses for each of the three tick-borne pathogens separately.

Because we were using host species as our units of comparison, the potential exists for inflating degrees of freedom if particular species’ traits are not independent of one another by virtue of phylogenetic similarity [27]. For our group of hosts we think this potential is low, because the maximum number of species per order is three (Rodentia), per family is two (two each for Soricidae and Sciuridae), and per genus is one. Similar analyses of life history correlates of reservoir competence for Lyme disease by Huang et al. [22] showed only modest quantitative (not qualitative) effects of using a phylogenetic correction, and that the uncorrected (conventional) analyses were conservative. For these reasons and for transparency of interpretation, we have chosen to present the uncorrected results. Nevertheless, in Supplementary Online materials, we present results of analyses using a phylogenetic correction based on a cladogram (Fig. S1) for the nine species we analyzed (Table S3).

Surrogates for tick/pathogen encounter rates

These variables were included to address whether reservoir competence for the tick-borne pathogens is potentially affected by the ability of the pathogens to adapt to the hosts they are most likely to encounter. We have no direct measurement of encounter rates between ticks or tick-borne pathogens and hosts, so we chose two surrogates. First, we assumed that overall encounter rates between ticks and a host species are a function of the average population density of that host. We base this assumption on the observation that blacklegged ticks are “sit-and-wait” parasites that move very little on their own and depend on host movements in order to gain access to a blood meal [28]. We are aware that host species-specific features other than population density, such as movement patterns, might additionally affect encounter rates with ticks, but we were unable to determine which behaviors might be important or to find relevant data for this community of hosts.

Although host population density might govern that species’ overall encounter rates with ticks, prior research shows that a large proportion of ticks encountering a host fail to feed to repletion, with some hosts grooming off and killing many ticks [29]. The product of tick encounter rates and the proportion of those ticks that survive to repletion (i.e., permissiveness of that host species; [29] determines the natural “body burden”, defined as the average number of ticks that feeds from a member of each host species under natural conditions [29]. Therefore, the average body burden on a host species times the average population density of that host determines the total number of ticks that feeds on that population of hosts in nature. Consequently, for our second surrogate for tick-host encounter rates we used host density * body burden, which probably is correlated with the opportunity for a tick-borne pathogen to invade a host species. Ideally, body burden data would be specific to the nymphal stage, because this is the stage responsible for most transmission events. However, in the absence of data on nymphal body burdens, we relied on our field data on larval body burdens for the nine mammal species [5], [29]. Values and data sources for the demographic variables density and density * burden are given in Table S2.

Comparing models with life history traits and tick/pathogen encounter rates

We then assessed levels of support for models in which reservoir competence was a function of specific life history traits or our surrogates for tick-pathogen encounter rates, or both. We used model comparison procedures to ask which combinations of the life history variables (the more comprehensive set), principal component scores, average population densities, and density * body burden were best supported by the data for each of the three pathogens separately. We used AIC corrected for small sample sizes (AIC_c) and Akaike weights to compare levels of support for each model, including those models with at least 10% weight.

Reservoir competence

Extensive studies have been undertaken to determine realized reservoir competence for a large component of the host community in an endemic zone for Lyme disease in southeastern New York State. Although these studies are specific to one site, the location is representative of the northeastern region of the United States with high Lyme disease incidence. Methods have been described in detail in [3]–[5]. We use previously published data for B. burgdorferi [5], B. microti [3], and the strain of A. phagocytophilum that infects humans (ha strain: [6]). Reservoir competence values for the three tick-borne pathogens are given in Table S2.

Results

Life history correlates

Principle components 1 and 2, respectively, explained 76.8% and 20.6% of the variation for the life history data set consisting of the three variables used by Huang et al. [22], and 74.5% and 15.3% of the variation for the more complete life-history data set (Fig. 1). For the 3-variable data set, high values of PC1 represented small bodies, short gestation, and large litters, whereas low values for PC2 represented small bodies and small litters. Mice, chipmunks, and the two shrew species were clustered in the lower right quadrant of Fig. 1A, representing fast life histories. For the larger data set, high values of PC1 were similar to the smaller data set (small bodies, short gestation, and many offspring per year) but additionally represented young age at maturity, short life span, and high mass-specific basal metabolic rate (BMR). Low values of PC2 were best represented by low BMR and few offspring/year (Fig. 1B). The four smallest mammals were clustered at high values of PC1 but spread more widely across the PC2 axis.

[Figure omitted. See PDF.]

Figure 1. Results of principle components analysis on life history variables used by Huang et al.

[22] and in the more inclusive data set (all variables). Vectors show the mapping of the original variables onto the first two principle component axes (PC1 and PC2). Host species are given by the following notation: OV = Odocoileus virginianus; PrL = Procyon lotor; MM = Mephitis mephitis; DV = Didelphis virginiana; ScC = Sciurus carolinensis; TS = Tamias striatus; BB = Blarina brevicauda; SoC = Sorex cinereus; and PeL = Peromyscus leucopus.

https://doi.org/10.1371/journal.pone.0107387.g001

For the 3-variable life history data set, multiple regression revealed that reservoir competence for B. burgdorferi was positively correlated with PC1 (t = 2.4, P = 0.05) and negatively correlated with PC2 (t = −3.3, P = 0.01). Reservoir competence for A. phagocytophilum was negatively correlated with PC2 (t = −2.6, P = 0.04) and uncorrelated with PC1 (t = 1.2, P = 0.27). For B. microti reservoir competence was not predicted by species scores on either of the principal components axes (PC1: t = 1.4, P = 0.206; PC2: t = −1.599, P = 0.16) (Fig. 2). For the 6-variable life history data set, multiple regression showed that reservoir competence for B. burgdorferi was positively correlated with PC1 (t = 3.0, P = 0.02), but uncorrelated with PC2 (t = −1.2, P = 0.27) (Table 1). B. microti and A. phagocytophilum were not correlated with either PC1 (B. microti: t = 1.7, P = 0.15; A. phagocytophilum: t = 1.4, P = 0.21) or PC2 (B. microti: t = −1.0, P = 0.37; A. phagocytophilum: t = −1.0, P = 0.35).

[Figure omitted. See PDF.]

Figure 2. Relationship between a host species’ score on the ordination of six life history traits, represented by two principal components (PC1 and PC2) and that host’s reservoir competence for Borrelia burgdorferi (Bb), Babesia microti (Bm), and the human-infectious strain of Anaplasma phagocytophilum.

https://doi.org/10.1371/journal.pone.0107387.g002

[Figure omitted. See PDF.]

Table 1. Multiple regression results for analyses of the relationship between the two principal components (PC1 and PC2) describing host life history and reservoir competence for three tick-borne pathogens.

https://doi.org/10.1371/journal.pone.0107387.t001

Demographic correlates

For B. burgdorferi and A. phagocytophilum we observed a strong positive correlation between population density of a host and its reservoir competence (t = 11.4, P<0.001; and t = 3.4, P = 0.01, respectively). A similar, positive relationship was also apparent for B. microti, but the correlation was marginally non-significant (t = 2.2, P = 0.06) (Fig. 3).

[Figure omitted. See PDF.]

Figure 3. Relationship between average host population density [5] and reservoir competence of that host species.

Host population density is used as one of two proxies for encounter rates between ticks and host species.

https://doi.org/10.1371/journal.pone.0107387.g003

We found that reservoir competence for B. burgdorferi was significantly correlated with host density * tick burden (t = 2.5, P = 0.04), but reservoir competence for B. microti and A. phagocytophilum were not, although we observed a marginally non-significant positive correlation in the case of A. phagocytophilum (t = 1.9, P = 0.10) (Fig. 4).

[Figure omitted. See PDF.]

Figure 4. Relationship between the product of average host population density and measured body burden (number of naturally occurring larval ticks per host individual) and reservoir competence of that host species.

The product of population density and body burden, which represents the average number of ticks fed by a particular population of hosts, is used as one of two proxies for encounter rates between ticks and host species.

https://doi.org/10.1371/journal.pone.0107387.g004

Model comparisons

For each of the three pathogens, we computed AIC_c values and Akaike weights for multiple regressions that included all combinations of the first two principal component axes (PC1 and PC2), the six life history predictors, and both demographic predictor variables. For B. burgdorferi, the best-supported models were those including population density (positive correlation) plus density * body burden (positive) (52.3% model weight) and density alone (33.3% model weight) (Table 2). For B. microti, the best-supported models were those including body mass (negative), the null (intercept only) model, and density (positive), with 20.2%, 14.5%, and 14.3% weight, respectively. For A. phagocytophilum, the best supported model was that including density (positive: 44.8% model weight) followed by density plus age at maturity (positive; 8.6% model weight). No other models had more than 10% model weight.

[Figure omitted. See PDF.]

Table 2. Results of model comparisons in which all the two principal component axes, PC1 and PC2, as well as the individual life history traits and the two surrogates of pathogen-host encounter probabilities are incorporated as explanatory variables for reservoir competence.

https://doi.org/10.1371/journal.pone.0107387.t002

Discussion

For a suite of three tick-borne zoonotic pathogens in eastern North America, two small rodents and two small soricimorphs consistently are among the most competent natural reservoirs for tick infection. Although many other species of terrestrial vertebrates are commonly exposed to these pathogens, they are comparatively inefficient reservoirs, infecting only a small proportion of the ticks that feed on them [11], [29]. Following the analysis of Lyme disease reservoirs by Huang et al. [22], we asked whether life history features of the members of this community of mammalian hosts for blacklegged ticks predict their reservoir competence for two other quite distinct pathogens. We also asked whether two measures of host-tick encounter probability influence reservoir competence. For all three pathogens we found that high average population density and “fast” life history features were predictors of high reservoir competence. However, the importance of life history traits and the specific traits linked to high reservoir competence varied among the pathogens. In the case of B. burgdorferi, model comparison revealed that all of the best-supported models (>10% Akaike weights) included high population density as a determinant of high reservoir competence, with some support for models including high basal metabolic rate. In contrast, for B. microti, small body mass, and to a lesser extent high population density were supported as determinants of reservoir competence. For this pathogen, the null model was almost equally as well supported as were models with body mass and population density. High population density was again supported for models of A. phagocytophilum, with early age at maturity also correlated with high reservoir competence. Surprisingly, neither of the composite life-history variables, represented as principal component axes, entered any of the best-supported models for any of the three pathogens.

Host life history and reservoir competence

Our results based on life history traits alone suggested that high reservoir competence for B. burgdorferi, and possibly A. phagocytophilum, was associated with high BMR and many offspring/year, characteristic of species with fast life histories. This conclusion is broadly consistent with research by Huang et al. [22], who compared reservoir competence metrics across mammalian and avian communities and found that hosts with faster life histories were better reservoirs for the agent of Lyme disease, West Nile virus, and Eastern Equine Encephalitis virus. These results are also consistent with similar analyses of an aphid-transmitted multi-host pathogen (barley-cereal yellow dwarf virus) of herbaceous plants [17] and trematode parasites of amphibians [20], showing that fast-lived host species were more prone to infection or high reservoir potential.

Fast life history strategies are thought to represent heavy investment in rapid growth, maturity, and reproduction at the expense of other factors, namely immune defenses [12], [13], [16], [30]–[32]. According to the “ecoimmunological pace of life” hypothesis [12], [33], [34], as a result of trade-offs between pace of life and self-defense, fast-lived species should be less resistant to infection and proliferation by some pathogens, which could increase reservoir competence [16], [32], [35], [36]. Alternatively, fast-lived species might tend to be more tolerant of infections [37], [38], but see [20]. Regardless, identifying life history features that predict reservoir competence for multi-host pathogens will help in predicting the consequences of changing host communities for disease risk [18], [21], [22]. Ideally, the search for predictive life history features would be informed by more thorough knowledge of specific relationships and trade-offs between immune functions and determinants of survival rates and fecundity.

Future research mechanistically linking life history features to reservoir competence for these pathogens should address two key questions: (1) which specific life history traits influence immune functions; and (2) which specific immune traits determine reservoir competence? Because experimentally manipulating multiple life history features, independently and in combination, seems infeasible, comparative approaches based on large samples of hosts are warranted. Similarly, challenges in manipulating multiple immune traits within hosts would appear to dictate a comparative approach using many species or populations of hosts.

Host availability and reservoir competence

We found support for a positive relationship between host population density and reservoir competence for all three pathogens. This result is consistent with the hypothesis that multi-host pathogens evolve to exploit the hosts that they are most likely to encounter [20], [24], [39]. We found only weak support for the hypothesis that reservoir competence correlates with total tick feeding opportunities, represented by the product of host density and tick burdens. It is important to note that population density and body mass are inversely correlated for our group of hosts. In fact, it is often the case that species with fast life histories tend to occur at higher population densities than their slower-lived counterparts [25]. Therefore, it is not possible to distinguish between features intrinsic to hosts (e.g., life history and immune strategies) and extrinsic features that represent pathogen-encounter probabilities (e.g., host density and tick burden). It is entirely plausible that high reservoir competence is caused both by pathogen adaptation to frequently-encountered hosts and fast-paced hosts being more poorly defended against pathogens. In addition, collinearity among independent variables complicates our ability to interpret results. Although we used ordination to replace collinear life history variables with orthogonal ones, none of the principal components explained patterns of reservoir competence as well as some of the original variables. Quantifying the relative importance of these potential drivers of the association between small-bodied, abundant species and reservoir competence will require a concerted effort to determine how pathogen adaptation and host immune functions independently and together determine reservoir competence.

In addition to limited ability to infer causality from correlations, a key caveat of our approach is the use of average values among individuals to estimate reservoir competence, life history traits and demographic features. Heterogeneities among individuals (and populations) certainly exist and should be incorporated in future efforts [32]. For example, intraspecific variation in reservoir competence for B. burgdorferi has been described [11], [40], but individual heterogeneities in life history from these same populations have not. Variance in host traits by itself could be an important determinant of reservoir competence at the population level. For instance, a host species with high mean susceptibility but low variance might strongly foster infection and proliferation by a pathogen, whereas a species with similar mean but high variance might provide fewer opportunities for pathogen transmission [41], [42]. Simultaneous estimates of reservoir competence, life history, and demography pertaining to individual hosts and host populations would facilitate stronger inference.

The patterns linking demographic and life history features with reservoir competence will be less likely to apply to specialist pathogens than to the generalist pathogens we analyzed. Strict specialists will experience only one host; therefore, interspecific comparisons among host traits will not be relevant to transmission dynamics. On the other hand, such comparisons might pertain if specialist pathogens frequently evolve from more generalized ancestors [43]. With ample opportunity to infect multiple species of host, a generalist pathogen will be likely to experience trade-offs that affect its ability to colonize and proliferate within each host [23]. Such trade-offs could contribute to increasing specialization that would then reduce selection pressures from other (non-specialized) hosts.

One potential application of our results is increased ability to predict which hosts might act as “keystone reservoirs” when generalist pathogens emerge. Our results suggest that of the host species that might potentially be exposed, those with comparatively high population densities, small bodies, and fast pace of life will often be keystone reservoirs that should be targeted for surveillance or management.

Conclusions

We examined the traits of mammals that correlate with competence in transmitting zoonotic, tick-borne pathogens to feeding tick vectors. We found support for the hypothesis that the most competent wildlife reservoirs for three common tick-borne pathogens are hosts that tend to occur at high population densities and have fast life histories compared to other mammalian hosts that are poorer reservoirs. However, the specific traits that correlated best with reservoir competence for the agents of Lyme disease, babesiosis, and anaplasmosis differed to some degree. The relative importance of extrinsic (population density) vs intrinsic (life history) features remains to be fully determined, although we found stronger support for the former. Specific causal relationships between demography, life history, immune functions, and reservoir competence, are important frontiers for future research.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Phylogeny of the mammalian hosts for tick-borne zoonotic pathogens analyzed for their life history features, tick-encounter surrogates, and reservoir competence. A phylogenetic correction, based on this cladogram, was applied before conducting the analyses the results of which are provided in Table S3.

https://doi.org/10.1371/journal.pone.0107387.s001

(TIF)

Table S1.

Data from PanTHERIA and AnAge used in the analyses of effects of life history variables on reservoir competence for the three tick-borne pathogens.

https://doi.org/10.1371/journal.pone.0107387.s002

(DOCX)

Table S2.

Data for each of the nine hosts analyzed for impacts of life history variables and tick-encounter surrogates. RC = reservoir competence for Borrelia burgdorferi, Babesia microti, and the human-infectious (ha) strain of Anaplasma phagocytophilum.

https://doi.org/10.1371/journal.pone.0107387.s003

(DOCX)

Table S3.

Results of model comparisons in which all the two principal component axes, PC1 and PC2, as well as the individual life history traits and the two surrogates of pathogen-host encounter probabilities are incorporated as explanatory variables for reservoir competence. This analysis differs from that summarized in the text and Table 2 in that a phylogenetic correction (see Materials and Methods) was conducted prior to analysis. The coefficients for each explanatory variable, AICc values (corrected for small samples), and model weights are provided for at least the top two models and all models with ≥10% AIC weight.

https://doi.org/10.1371/journal.pone.0107387.s004

(DOCX)

Author Contributions

Conceived and designed the experiments: RSO TL AEJ LBM PRH FK. Performed the experiments: RSO FK. Analyzed the data: RSO TL FK. Contributed reagents/materials/analysis tools: TL AEJ LBM FK. Contributed to the writing of the manuscript: RSO TL AEJ LBM PRH FK.

Citation: Ostfeld RS, Levi T, Jolles AE, Martin LB, Hosseini PR, Keesing F (2014) Life History and Demographic Drivers of Reservoir Competence for Three Tick-Borne Zoonotic Pathogens. PLoS ONE 9(9): e107387. https://doi.org/10.1371/journal.pone.0107387

References

1. Kuehn BM (2013) CDC estimates 300,000 US cases of Lyme disease annually. JAMA-J Am Med Assoc 310: 1110

2. Hersh MH, Ostfeld RS, McHenry DJ, Tibbetts M, Brunner JL, et al. (2014) Co-Infection of blacklegged ticks with Babesia microti and Borrelia burgdorferi is higher than expected and acquired from small mammal hosts. PLoS ONE 9(6): e99348

3. Hersh MH, Tibbetts M, Strauss M, Ostfeld RS, Keesing F (2012) Reservoir competence of wildlife host species for Babesia microti. Emerg Infect Dis 18: 1951–1957.

4. Keesing F, Hersh MH, Tibbetts M, McHenry DJ, Duerr S, et al. (2012) Reservoir competence of vertebrate hosts for Anaplasma phagocytophilum. Emerg Infect Dis 18: 2013–2016.

5. LoGiudice K, Ostfeld RS, Schmidt KA, Keesing F (2003) The ecology of infectious disease: effects of host diversity and community composition on Lyme disease risk. P Natl Acad Sci USA 100: 567–571.

6. Keesing F, McHenry DJ, Hersh M, Tibbetts M, Brunner JL, et al. (2014) Prevalence of human-active and variant-1 strains of the tick-borne pathogen Anaplasma phagocytophilum in hosts and forests of eastern North America. Am J Trop Med Hygiene 91: 302–309.

7. Gray J, von Stedingk LV, Gurtelschmid M, Granstrom M (2002) Transmission studies of Babesia microti in Ixodes ricinus ticks and gerbils. J Clin Microbiol 40: 1259–1263.

8. Oliveira MR, Kreier JP (1979) Transmission of Babesia microti using various species of ticks as vectors. J Parasitol 65: 816–817.

9. Hotopp JCD, Lin M, Madupu R, Crabtree J, Angiuoli SV, et al. (2006) Comparative genomics of emerging human ehrlichiosis agents. PLoS Genetics 2: 208–223.

10. Patrican LA (1997) Absence of Lyme disease spirochetes in larval progeny of naturally infected Ixodes scapularis (Acari: Ixodidae) fed on dogs. J. Med. Entomol. 34: 52–55.

11. Brunner JL, LoGiudice K, Ostfeld RS (2008) Estimating reservoir competence of Borrelia burgdorferi hosts: prevalence and infectivity, sensitivity, and specificity. J Med Entomol 45: 139–147.

12. Ricklefs RE, Wikelski M (2002) The physiology/life-history nexus. Trends Ecol Evol 17: 462–468.

13. Martin LB, Weil ZM, Nelson RJ (2007) Immune defense and reproductive pace of life in Peromyscus mice. Ecology 88: 2516–2528.

14. Martin LB, Weil ZM, Nelson RJ (2008) Seasonal changes in vertebrate immune activity: mediation by physiological trade-offs. Philos T Roy Soc B 363: 321–339.

15. Martin LB, Weil ZM, Nelson RJ (2008) Fever and sickness behaviour vary among congeneric rodents. Funct Ecol 22: 68–77.

16. Lee KA, Wikelski M, Robinson WD, Robinson TR, Klasing KC (2008) Constitutive immune defences correlate with life-history variables in tropical birds. J Anim Ecol 77: 356–363.

17. Cronin JP, Welsh ME, Dekkers MG, Abercrombie ST, Mitchell CE (2010) Host physiological phenotype explains pathogen reservoir potential. Ecol Lett 13: 1221–1232.

18. Keesing F, Belden LK, Daszk P, Dobson A, Harwell CD, et al. (2010) Impacts of biodiversity on the emergence and transmission of infectious diseases. Nature 468: 647–652.

19. Miller BLW, Sinervo B (2007) Heritable body size mediates apparent life-history trade-offs in a simultaneous hermaphrodite. J Evolution Biol 20: 1554–1562.

20. Johnson PTJ, Rohr JR, Hoverman JT, Kellermanns E, Bowerman J, et al. (2012) Living fast and dying of infection: host life history drives interspecific variation in infection and disease risk. Ecol Lett 15: 235–242.

21. Previtali MA, Ostfeld RS, Keesing F, Jolles A, Hanselmann R, et al. (2012) Relationship between pace of life and immune responses of wild rodents. Oikos 121: 1483–1492.

22. Huang ZYX, de Boer WF, van Langevelde F, Olson V, Blackburn TM, et al. (2013) Species’ life-history traits explain interspecific variation in reservoir competence: a possible mechanism underlying the dilution effect. PLoS ONE 8: e54341

23. Schmid-Hempel P (2011) Evolutionary Parasitology. Oxford, UK: Oxford University Press. 496 p.

24. Ostfeld R, Keesing F (2000) The function of biodiversity in the ecology of vector-borne zoonotic diseases. Can J Zool 78: 2061–2078.

25. Stearns SC (1992) The Evolution of Life Histories. Oxford, UK: Oxford University Press. 249 p.

26. Jones KE, Bielby J, Cardillo M, Fritz SA, O’Dell J, et al. (2009) PanTHERIA: a species-level database of life history, ecology, and geography of extant and recently extinct mammals. Ecology 90: 2648

27. Felsenstein J (1985) Phylogenies and the comparative method. Am Nat 125: 1–15.

28. Ostfeld RS (2011) Lyme Disease: The Ecology of a Complex System. New York: Oxford University Press. 232 p.

29. Keesing F, Brunner J, Duerr S, Killilea M, LoGiudice K, et al. (2009) Hosts as ecological traps for the vector of Lyme disease. Proc R Soc B 276: 3911–3919.

30. Martin LB, Weil ZM, Kuhlman JR, Nelson RJ (2006) Trade-offs within the immune systems of female White-footed Mice, Peromyscus leucopus. Funct Ecol 20: 630–636.

31. Viney ME, Riley EM, Buchanan KL (2005) Optimal immune responses: immunocompetence revisited. Trends Ecol Evol 20: 665–669.

32. Hawley DM, Altizer SM (2011) Disease ecology meets ecological immunology: understanding the links between organismal immunity and infection dynamics in natural populations. Funct Ecol 25: 48–60.

33. Martin LB, Pless M, Svoboda J, Wikelski M (2004) Immune activity in temperate and tropical house sparrows: a common-garden experiment. Ecology 85: 2323–2331.

34. Martin LB, Gilliam J, Han P, Lee K, Wikelski M (2005) Corticosterone suppresses cutaneous immune function in temperate but not tropical House Sparrows, Passer domesticus. Gen Comp Endocr 140: 126–135.

35. Knowles SCL, Nakagawa S, Sheldon BC (2009) Elevated reproductive effort increases blood parasitaemia and decreases immune function in birds: a meta-regression approach. Funct Ecol 23: 405–415.

36. Lee KA, Martin LB, Hasselquist D, Ricklefs RE, Wikelski M (2006) Contrasting adaptive immune defenses and blood parasite prevalence in closely related Passer sparrows. Oecologia 150: 383–392.

37. Raberg L, Graham AL, Read AF (2009) Decomposing health: tolerance and resistance to parasites in animals. Philos T Roy Soc B 364: 37–49.

38. Sears BF, Rohr JR, Allen JE, Martin LB (2011) The economy of inflammation: when is less more? Trends Parasitol 27: 382–387.

39. Keesing F, Holt RD, Ostfeld RS (2006) Effects of species diversity on disease risk. Ecol. Lett. 9: 485–498.

40. Brunner JL, Ostfeld RS (2008) Multiple causes of variable tick burdens on small-mammal hosts. Ecology 89: 2259–2272.

41. Beldomenico PM, Begon M (2010) Disease spread, susceptibility and infection intensity: vicious circles? Trends Ecol Evol 25: 21–27.

42. Tompkins DM, Dunn AM, Smith MJ, Telfer S (2011) Wildlife diseases: from individuals to ecosystems. J Anim Ecol 80: 19–38.

43. Wolfe ND, Dunavan CP, Diamond J (2007) Origins of major human infectious diseases. Nature 447: 279–283.