Uncovering drivers of dose-dependence and individual variation in malaria infection outcomes

Abstract

To understand why some hosts get sicker than others from the same type of infection, it is essential to explain how key processes, such as host responses to infection and parasite growth, are influenced by various biotic and abiotic factors. In many disease systems, the initial infection dose impacts host morbidity and mortality. To explore drivers of dose-dependence and individual variation in infection outcomes, we devised a mathematical model of malaria infection that allowed host and parasite traits to be linear functions (reaction norms) of the initial dose. We fitted the model, using a hierarchical Bayesian approach, to experimental time-series data of acute Plasmodium chabaudi infection across doses spanning seven orders of magnitude. We found evidence for both dose-dependent facilitation and debilitation of host responses. Most importantly, increasing dose reduced the strength of activation of indiscriminate host clearance of red blood cells while increasing the half-life of that response, leading to the maximal response at an intermediate dose. We also explored the causes of diverse infection outcomes across replicate mice receiving the same dose. Besides random noise in the injected dose, we found variation in peak parasite load was due to unobserved individual variation in host responses to clear infected cells. Individual variation in anaemia was likely driven by random variation in parasite burst size, which is linked to the rate of host cells lost to malaria infection. General host vigour in the absence of infection was also correlated with host health during malaria infection. Our work demonstrates that the reaction norm approach provides a useful quantitative framework for examining the impact of a continuous external factor on within-host infection processes.

About the Authors:

Tsukushi Kamiya

Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Software, Validation, Visualization, Writing – original draft

* E-mail: [email protected]

Current address: Department of Ecology & Evolutionary Biology, University of Toronto, Toronto, ON M5S 3B2, Canada

Affiliation: Department of Ecology & Evolutionary Biology, University of Toronto, Toronto, ON M5S 3B2, Canada

ORCID logo http://orcid.org/0000-0002-9020-6699

Megan A. Greischar

Roles Conceptualization, Investigation, Methodology, Supervision, Writing – review & editing

Affiliation: Department of Ecology Evolutionary Biology, Cornell University, United States of America

ORCID logo http://orcid.org/0000-0002-7521-9344

David S. Schneider

Roles Conceptualization, Investigation, Resources, Supervision, Writing – review & editing

Affiliations Program in Immunology, Stanford University, Stanford, California, United States of America, Department of Microbiology and Immunology, Stanford University, Stanford, California, United States of America

ORCID logo http://orcid.org/0000-0002-2391-9963

Nicole Mideo

Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing

Affiliation: Department of Ecology & Evolutionary Biology, University of Toronto, Toronto, ON M5S 3B2, Canada

ORCID logo http://orcid.org/0000-0002-8719-3620

Introduction

Infections produce divergent outcomes. In human malaria, for example, outcomes of infection with the same parasite, Plasmodium falciparum, range from sub-clinical to fatal [1]. Understanding drivers of variation in infection outcomes is central to explaining why some hosts get sicker than others. Some host and parasite factors underlying this variation have well-understood mechanisms. For example, heterozygosity in the haemoglobin coding gene (i.e., sickle-cell trait, or HbAS) confers partial protection against falciparum malaria: sickle-cell trait individuals experience lower parasite load and reduced likelihood of life-threatening cerebral malaria and severe anaemia [2]. The resistance mechanism of this single locus trait has been corroborated by four decades of research demonstrating that sickling enhances clearance of infected red blood cells (iRBCs) by host immune effectors like macrophages [2]. However, unlike the sickle-cell trait, there are numerous sources of heterogeneity—including in the initial infection dose, nutrition, coinfection, and other genetic factors—for which quantitative impacts on outcomes have been observed, but causal mechanisms have remained elusive [3–6]. To establish a causal link between complex factors and infection outcomes, a key challenge is to quantify how a factor of interest mediates key host and parasite processes, such as host responses to infection and parasite growth [7].

In many disease systems, the initial infection dose is a key biotic factor that varies widely across infection events [8, 9]. Experimental infections in diverse systems show that increasing infection dose negatively impacts host fitness through reduced host vigour, survival and fecundity [10–17], which is likely linked to variation in the within-host parasite dynamics due to dose-dependence in parasite growth and host immune responses [4, 12, 18–21]. Complex interactions between immune responses and the initial infection dose have been revealed by molecular immunology studies of viral and bacterial systems. For example, the expression of pro-inflammatory cytokines (i.e., signalling molecules) and immune cells, can decrease with infection dose [22, 23], presumably due to enhanced evasion and escalated damage of host immune machinery through an increased abundance of reactive oxygen species [24]. Conversely, it has been shown that higher doses trigger distinct, and sometimes more robust, activation of innate and adaptive immune pathways [25–28]. From an applied perspective, a better understanding of the immunogenic dose-response relationship is pertinent for optimising vaccine dosing to ensure improved safety and efficacy [21, 29–31]. However, because the functional output of immune activities (e.g., the rate of immune-mediated iRBC clearance) is difficult to measure directly, it remains an open question how the initial infection dose influences host responses overall.

In malaria infections, the initial density of iRBCs at the start of blood-stage infection likely ranges in the order of hundreds to over millions [32–36], with greater numbers generally shown to increase mortality and worsen morbidity [13, 15]. Experimental manipulations of the initial infection dose, ranging from 100 to 100 million iRBCs, have demonstrated that larger doses increase the pace of infection with each order of magnitude reducing the time until peak infection by roughly 24 hours [15] (Fig 1a). Dose also impacts the peak parasite load: mice initially infected with 100 million parasites harbour roughly 60% more iRBCs at peak compared to those infected with 100 parasites [15] (Fig 1b). Furthermore, high doses induced more severe anaemia measured by the minimum red blood cell (RBC) count in mice on average [15] (Fig 1c). Previous mathematical modelling studies interpret these dose-dependent infection outcomes as a reflection of the underlying dose-dependence in host immune responses [4, 37].

[Figure omitted. See PDF.]

Fig 1.

Higher initial infection doses increase a) the pace of infection (i.e., time until peak iRBC density), b) peak parasite load measured as the maximum iRBC count, and c) severity of anaemia during malaria infection measured as the minimum RBC count. There is also considerable variation in quantitative infection outcomes (i.e., parasite load and anaemia) within infection dose treatments. Data from Timms et al. [15] with 5 to 9 mice infected with the CW strain of P. chabaudi, at each dose. The crosses indicate data and the white lines and grey bands correspond to the means and 95% confidence intervals.

https://doi.org/10.1371/journal.pcbi.1008211.g001

While dose clearly influences malaria infection dynamics, these experimental data also reveal striking variability within dose treatments [15] (Fig 1), meaning that quantitatively diverse infection outcomes were observed across individuals receiving the same infection dose. This is despite the fact that hosts were inbred to homozygosity and parasites were also of single strain origin in the experiment [15]. Identifying the sources of such individual variation—usually considered experimental “noise”—may reveal biologically interesting, subtle trait variation among hosts and/or parasites, and thus new therapeutic targets (e.g., host responses to boost).

In the study of acute malaria infection, mathematical models have been fitted to the time course of experimental infections in mice to provide a quantitative understanding of parasite growth, pathogenesis and host responses to infection (e.g., [38–41]; see [42] for a review of earlier work). Infection triggers a variety of host responses, for example, indiscriminate clearance of RBCs, targeted clearance of iRBCs [39] and production of new RBCs to compensate for those lost to infection [39, 41]. It is well documented that these responses involve a complex cascade of interactions across multiple organisational scales from molecules, cells, and tissues to organ systems [43, 44]. However, it remains a challenge to scale up the details of finer level processes to an understanding of the net effect of host responses on parasite load and host health [45]. Data-driven mathematical modelling allows for the inference of these net effects, without necessarily requiring a detailed understanding of the underlying mechanisms.

Here, we fitted a dynamical model of within-host malaria infection to experimental data spanning seven orders of magnitude of initial doses, using a hierarchical Bayesian approach. By modelling the influence of dose on model parameters as a reaction norm, which describes the pattern of phenotypic expression of an organism across an environmental gradient, we identified drivers of the observed dose-dependent malaria parasite load and severity of malaria-induced anaemia. By explicitly modelling individual variation as model parameters, we also examined the origin of quantitatively diverse infection outcomes observed within single initial infection dose treatments.

Methods

System and experimental set-up

The rodent malaria system offers unique opportunities to investigate infection ecology, pathogenesis and host responses [46]. We examined previously published experimental data of C57BL/6 female mice infected with the CW strain of Plasmodium chabaudi [15]. In this experiment [15], infection was initiated with an intraperitoneal injection of iRBCs at seven different doses: 102, 103, 104, 105, 106, 107, 108—and considerable variation in quantitative infection outcomes was observed both among and within dose treatments (Fig 1b & 1c). Details of the experiment are provided by Timms et al. [15] and data are available on Dryad [47].

Model

Innate host responses to malaria infection.

Hosts trigger a variety of responses to resist, tolerate and/or recover from infections. Here, we focused on two forms of rapid immune responses (on the order of minutes [48, 49]) that have been identified as the most pertinent to describing the acute blood-stage malaria infection [38–40].

The first response we modelled was general clearance of RBCs which may involve mechanisms such as retention of RBCs by the spleen and destruction of RBCs by immune effector cells [50, 51]. Clearing RBCs indiscriminately has been proposed as a host adaptation in the presence of malaria parasites to directly clear the parasites (i.e., top-down effect) as well as to limit the resource for the parasite (i.e., bottom-up effect) [52]. The second response we considered was the induction of innate immunity targeting iRBCs only, which is considered predominantly responsible for controlling the acute phase of malaria infection [43].

We modelled regulation of host responses without delving into fine mechanistic details (i.e., we avoided mathematical descriptions of cytokine storms and subsequent cascades of effector responses). In part, this modelling choice was out of necessity because there is no complete map of innate immune responses against malaria [46]. Yet, it was also by design so that we would gain a functional understanding of host responses with minimal complexity. Biologically, responses modelled here may reflect the output of an entire module of proteins and signal transduction pathways that lead to the production of effector cells. Specifically, we used a single ordinary differential equation (ODE) to describe the change in the functional output of each response (Fig 2a; yellow and green block, respectively). We assumed that the maximum possible activity of each response is fixed, at one, and we tracked the dynamics of its proportional activity, Ni, where i indicates the response identity (general RBC clearance, i = 1; targeted iRBC clearance i = 2)(1)

[Figure omitted. See PDF.]

Fig 2.

a) A dynamical regulation model of host responses against blood-stage malaria. We condensed the complexity of the vertebrate innate response against malaria into two independent pathways responsible for general RBC clearance and targeted iRBC clearance (represented by the yellow and green block, respectively). We modelled each pathway using a single differential equation, the activity of which is denoted Ni where the subscript i indicates the identity of each response: general RBC clearance, i = 1; and targeted iRBC clearance, i = 2. (The schematic shows the example of i = 1.) For each response type, the host detects a within-host cue, Ci. The product of the cue and the strength constant, ψi, activates the response. The activity of a response decays spontaneously with a half-life, ϕi. The output of each host response feeds back to influence the within-host infection dynamics (indicated by the coloured arrows in panel b). b) Dynamics of RBCs and blood-stage malaria parasites within the host. Recruitment into and transitions among components of the asexual cycle are indicated with black arrows. Background mortality for different components is indicated by grey arrows. General clearance of RBCs and targeted clearance of iRBCs are marked with yellow and green arrows, respectively. Replenishment of RBCs (erythropoiesis) is indicated in blue.

https://doi.org/10.1371/journal.pcbi.1008211.g002

We defined the activity of N1 and N2 as the proportion of RBCs and iRBCs cleared by indiscriminate and targeted mechanisms per day, respectively. We assumed that there is no response output in the absence of infection, i.e., N1(t = 0) = N2(t = 0) = 0, consequently assuming a stable RBC population and that there is no constitutive immune activity.

We modelled the signalling input that activates each response as a function of a within-host cue, Ci(t) and a constant determining the strength of activation, ψi (Eq 1). Host innate immune responses against malaria are thought to be triggered by pathogen-associated molecular patterns (PAMPs) such as GPI anchors, haemozoin, parasite DNA and RNA [44, 46]. Assuming that the abundance of PAMPs reflects that of iRBCs, our model considers the relative density of iRBCs compared to its observed maximum in any infection across all treatment groups, as the within-host cue for general RBC clearance and targeted iRBC clearance, i.e., where I(t) and max I are the iRBC density at time t and the maximum observed iRBC density, respectively. The latter was reported at 2.18 × 106 from this dataset [15]. We assumed that each response activity decays spontaneously with a half-life of ϕi.

Our two-parameter approach—involving only an activation constant, ψi, and activity half-life, ϕi—to modelling each host response was inspired by Kochin et al. [38] who used a single ODE to model innate immunity against malaria parasites governed by density-dependent response activation and constant decay. However, we interpreted host response activity differently from their study: i.e., we modelled the proportion of RBCs and iRBCs cleared per day whereas they modelled the number of immune cells. Our study also extends the approach to the dynamics of general RBC clearance.

Within-host infection dynamics of blood-stage rodent malaria.

We used a system of ODEs to model the blood-stage asexual cycle of P. chabaudi, tracking the density of uninfected RBCs (hereafter, uRBCs), R, iRBCs, I, and extracellular parasites called merozoites, M, in a microlitre of blood (Fig 2b). In this model, we assumed that RBCs are constantly replenished to maintain a homeostatic equilibrium, thus the daily rate of erythropoiesis in the absence of infection is defined as Rc μR, where μR is the baseline RBC mortality rate. We estimated Rc at 8.89 × 106 from data [15] as the average RBC density of 10 uninfected mice between Day 7 and 14 during which time the RBC density appears stable. During the acute phase of malaria infection, the host upregulates erythropoiesis to restore RBCs lost to malaria-induced anaemia [53, 54]. Following a previous study [39], we modelled this upregulation as a product of the deviation from the homeostatic equilibrium, Rc − R(t) where R(t) is the RBC density at time t, and the proportion of the deviation from the homeostatic equilibrium restored by the host per day, ρ.

Given that N1 and N2 were defined as the proportion of cells cleared by indiscriminate and targeted mechanisms per day, respectively, it was convenient to convert the proportions into daily rates at which cells are cleared in the dynamical within-host model. To do this, we equate , where X is the daily rate of clearance. Solving for X, we obtained the rate of general RBC clearance as −ln(1 − N1) and the rate of targeted iRBC clearance as −ln(1 − N2). Therefore, the sum of the baseline rate, μR and −ln(1 − N1) constitutes the daily mortality rate of uRBCs. uRBCs then become infected at a rate proportional to the density, M, and invasion rate, p, of merozoites. Together the dynamics of uRBCs is expressed as: (2)

Following merozoite invasion, iRBCs remain subjected to background mortality, μR and general RBC clearance, −ln(1 − N1). In addition, infected cells are cleared by targeted immunity at a rate −ln(1 − N2): here, we note that estimates of −ln(1 − N2) may be slightly inflated by iRBCs that commit to sexual reproduction (usually less than 2% of iRBCs [55]) because our model does not consider Plasmodium sexual reproduction. We modelled the development of iRBCs using a gamma-chain trick (also known as linear chain trick) [56, 57], which consists of a series of ODEs:(3)(4)where , and α is the average cell cycle duration, which is 24 hours for P. chabaudi [58]. The number of compartments in the series, n, reflects the assumption about the variance in the developmental time, which is inversely proportional to n [56]. At one compartment per cycle (i.e., n = 1, e.g., [59, 60]), the assumption is that the developmental time is exponentially distributed, with large variance, i.e., α2 [61]. The variance decreases with the number of compartments, and the variance tends to disappear as the number of compartments approaches infinity (i.e., n → ∞)—a scenario equivalent to discrete-time models (e.g., [62, 63]) and delay-differential equation models [64] that assume there is no variation in the developmental period [56]. It has been shown that models with few compartments tend to estimate greater asexual multiplication, compared to discrete-time models [61, 65]. With enough compartments (e.g., one compartment per hour of parasite development [56]), however, the outcome of a gamma-chain ODE model converges to the discrete-time model [65]. We arrived at the choice of n = 12 (one compartment every two hours) for computational efficiency, and because our preliminary analysis showed that the infection dynamics were quantitatively comparable to that of a 24 compartment model (one compartment per hour).

Finally, the production of merozoites is determined by parasite burst size (i.e., the number of progeny parasites emerging from an iRBC), β, and the number of bursting cells, τ In(t). Merozoites are lost as they invade new RBCs and through background mortality, μM. We ignored immune-mediated clearance of merozoites because its effect on the parasite dynamics is functionally similar to clearing iRBCs and comparatively less important for describing the acute phase of malaria infection than general RBC clearance or targeted iRBC clearance [4, 39],(5)

Initial conditions.

We set the initial RBC density, R(t = 0), to the values reported per mouse by Timms et al. [15]. In the experiment [15], malaria infection was initiated with an intraperitoneal injection to mimic the initial cohort of blood-stage malaria parasites following release from the liver. Assuming that the initial parasite growth is near-exponential (for the first three records of iRBCs per mouse), we estimated the initial infection dose per μl of blood in each mouse as the intercept of a linear regression model with the natural logarithm of iRBCs as the response and the time since infection as a predictor. At a preliminary phase, we also estimated the initial infection dose simultaneously with the rest of the model parameters. Because these two methods for estimating initial dose yielded analogous results (see S1a Appendix) and the regression method allowed us to estimate one fewer parameter in the main Bayesian parameter inference procedure (described below), we present results based on the initial infection dose estimated by the regression method. Finally, we defined the age structure of the inoculated iRBCs to schedule the bursting of iRBCs in the initial cohort. Assuming that all inoculated iRBCs commit to producing merozoites (i.e., ignoring the possibility that some of them produce transmission stages instead), the timing of bursting can be described using a symmetrical beta distribution [64]. We assumed a moderately synchronised blood-stage cycle (i.e., with the shape parameter, s = 10). Our preliminary exploration indicated that the daily RBC and iRBC measurements were insensitive to the s parameter. We discretised the beta distribution into 12 compartments (for Eqs 3 and 4) by dividing the cumulative density function into 12 intervals.

Hierarchical Bayesian inference

Bayesian causal inference is an effective tool to paint a picture of processes that generated data [67]. We fitted our ODE model describing the dynamics of RBCs (Eq 2) and iRBCs (Eqs 3 and 4) to the corresponding time-series data from 51 mice [15] using a Bayesian statistical approach. Statistical inference in a Bayesian framework incorporates prior knowledge and uncertainty about model parameters and updates the belief about them (by computing the posterior probability of parameters) based on observations (through a likelihood function) [68]. Besides its conceptual merits [68, 69], Bayesian inference has a pragmatic appeal for its ability to estimate parameters in high dimensional spaces, for example, in hierarchical models, which are used when observations are organised in multiple levels of sampling units [70, 71]. In this study, we considered two levels of sampling units: treatments (i.e., infection doses) and subjects (i.e., individual mice).

Dose and individual-specific effects.

For the dynamical model (Eqs 1 to 5), we estimated dose- and individual-specific effects in a set of seven parameters including the response activation strength of host responses N1 and N2 (ψ1 and ψ2, respectively), half-life of those responses (ϕ1 and ϕ2, respectively), erythropoiesis upregulation (ρ) and parasite burst size (β). Below, we collectively refer to the parameter set as θ (θ ∋ ψ1, ψ2, ϕ1, ϕ2, ρ, β) and refer to a parameter in the set using an index, k. The prior distributions for these parameters are provided in Table 1 and further detailed in S1b Appendix.

[Figure omitted. See PDF.]

Table 1. Descriptions of model parameters and their fixed values, or prior distributions used in Bayesian statistical inference.

Estimated parameters are indicated by an asterisk on the description. We assigned a generic, weakly informative prior, except for erythropoiesis upregulation, ρ, burst size, β, and standard deviations of log10 RBC and iRBC density, σRBC and σiRBC, for which there exist specific prior information from previous studies. Further details of the priors and comparisons with estimated posterior probability densities, and prior sensitivity analyses are provided in S1b Appendix.

https://doi.org/10.1371/journal.pcbi.1008211.t001

Instead of modelling each dose treatment as a discrete group (or a character state), we consider the entire range of the initial infection dose as an environmental gradient. In other words, we modelled each parameter of the dynamical model, θk, as a reaction norm, which describes the pattern of phenotypic expression of an organism across an environmental gradient [72, 73]. The simplest and most commonly used reaction norm assumes a linear relationship between the environment and phenotype expression and consists of two components: the mean intercept, i.e., the phenotype expressed against the “average” environment, and the slope, i.e., the degree of phenotypic change along with the environment. We estimated the mean intercept, (dose was centred so that the intercept is at the middle dose of 105) and the mean slope of the reaction norms, .

Within each dose treatment, Timms et al’s dataset [15] contains repeated measures from replicate mice that showed marked individual variability (Fig 1b & 1c). We sought to identify the source of this variability by explicitly modelling individual variation in θk among mice through partial pooling. This means that a given parameter was considered a sample from a common population distribution with a mean—in this case the intercept, , and the slope, —and the deviation of the parameter from the mean for each individual, i, which we express as uk,i and vk,i for the intercept and slope variation, respectively. We assumed that the individual deviations, uk,i and vk,i, are samples of a normal distribution with standard deviations, σu,k and σv,k that we estimated from data [74].

Together, dose- and individual-specific parameter, θk,i, is expressed as:(6)where Dosei indicates the dose treatment applied to individual i. Dose was coded as {−3, −2, −1, 0, 1, 2, 3} such that the intercept was centred at the initial infection dose of 105. As is customary in quantitative genetics [72] and hierarchical modelling [75], we modelled a covariance structure describing the association among parameters in individual deviation following [75].

Likelihood.

A Bayesian approach requires a likelihood function to assess the probability of observing the data given model parameters and associated predictions. Our log-likelihood function assumed that the measurement error for the total density of RBCs (i.e., sum of uRBCs and iRBCs), and iRBCs is distributed normally and log10-normally, respectively [40, 76]:(7)where and are the observed count of total RBCs and iRBCs, and are the model predictions of total RBCs and iRBCs for individual i at time t. We estimated standard deviations, σRBC and σiRBC for the total RBC and iRBC count, respectively, with specific informative priors based on [76] (Table 1).

Our modelling focused on the first wave of infection before iRBCs recrudesce and adaptive immunity starts to take effect [39]. Thus, ntime was defined {16, 16, 13, 11, 11, 10, 8} for the seven dose treatments, respectively, noting that higher doses lead to shorter time series because of the faster pace of infection (Fig 1a). We further subsetted the dataset by removing instances of atypical dynamics (S1c Appendix for more details). In total, we fitted data from 51 individual mice (nmice = 51).

MCMC sampling.

Estimating the posterior probability density of parameters of a complex model requires a Markov Chain Monte Carlo (MCMC) sampling algorithm, which can be computationally intensive for large hierarchical models. Our model was written in Stan 2.18.2 and fitted through the RStan interface [77, 78], which provides an efficient general-purpose MCMC sampler (No-U-Turn Hamiltonian Monte Carlo) and a Bayesian inference environment. The model was fitted in parallel in four independent chains, each with 4000 sampled iterations and 1000 warmup iterations. For diagnostics, we confirmed over 400 effective samples and ensured convergence of independent chains using the metric (values below 1.1 are considered an indication of multi-chain convergence) for all parameters [68, 79]. The computer programmes used in the present study are available in the S1d Appendix.

Results and discussion

Our fitted model accurately describes the daily time course of RBCs and iRBCs during the acute phase of malaria infection in mice, initiated at doses spanning seven orders of magnitude, from 102 to 108 iRBCs (Fig 3, S1e Appendix). We found evidence for dose-dependence in key parameters of host responses underlying the dynamics of RBCs and iRBCs. Additionally, even under the highly controlled condition of Timms et al.’s experiment [15]—with the single strain combination of hosts (C57/BL6) and P. chabaudi parasites (CW)—we identified individual variation in host and parasite traits that influences the variation in infection outcomes independent of the dose treatment. Below, we closely examine the sources and impacts of dose-dependence and individual variation in different parameters of infection and initial conditions.

[Figure omitted. See PDF.]

Fig 3. The fit of the full model (with parameters defined by Eq 6) to the density of RBCs (orange) and iRBCs (blue) for individual mice inoculated at 102, 103, 104, 105, 106, 107 and 108 infected cells.

Each column corresponds to the initial infection dose given at the top. The crosses indicate data and grey bands correspond to 95% predictive intervals of the full model, incorporating uncertainty in parameter estimation and sampling.

https://doi.org/10.1371/journal.pcbi.1008211.g003

Dependence on the initial infection dose

We found evidence that increasing dose had two opposing effects on general RBC clearance. First, increasing dose reduces the activation strength of this response (Fig 4a). The lower peak responses estimated for high doses (Fig 5a) were attributable to a negative influence of the initial dose on the response activation strength (Fig 4a), supporting the notion that higher infection doses enhance the parasites’ ability to evade/suppress host responses and/or damage host machinery [24]. In contrast, we found that dose facilitates the same response by inducing a longer action (Fig 4b). The positive relationship between dose and activity half-life of these responses (Fig 4b) explained the relatively low peak response at low doses (Fig 5a). As a result of these two opposing effects of dose, the strongest host response in general RBC clearance was predicted at an intermediate dose, 105 (white squares in Fig 5a). This finding draws a comparison to another data-driving modelling work demonstrating that the maximum immune protection against influenza is generally achieved at an intermediate dose due to an interplay between innate and adaptive responses [21]. Further studies in other disease systems are desired to probe the generality of the intermediate peak and understand diverse mechanisms that may generate this pattern.

[Figure omitted. See PDF.]

Fig 4. Dose dependence in some parameters of host responses, but not parasite traits.

Plotted are the estimated relationship between the initial infection dose and a) the activation strength (ψ1 and ψ2), b) activity half-life (ϕ1 and ϕ2) of general RBC clearance (yellow) and targeted iRBC clearance (green), c) erythropoiesis upregulation, ρ, and d) parasite burst size, β. The line and band indicates the median prediction and 95% predictive intervals, respectively. The solid line indicates a statistically significant sign of dose-dependence (See S1b Appendix for the prior and posterior distributions).

https://doi.org/10.1371/journal.pcbi.1008211.g004

[Figure omitted. See PDF.]

Fig 5. Host responses against malaria depend on the initial infection dose.

Shown are the median predicted host responses (colour) as a function of time centred around the dose-dependent day of peak infection (x-axis) and the initial infection dose (y-axis) for a) general RBC clearance (i.e., proportion of RBCs, indiscriminately of their infection status, cleared by the host per day), b) targeted iRBC clearance (i.e., proportion of iRBCs cleared by the host per day), c) the total iRBC clearance (the sum of a and b). The grey region indicates that the model was not fitted for the day because either data do not exist or the day is beyond the first wave of infection.

https://doi.org/10.1371/journal.pcbi.1008211.g005

Following the peak infection day, malaria parasites find themselves in a hostile within-host environment as over 70% of the iRBC population per day is cleared by host immunity targeting iRBCs (Fig 5b). Combining targeted and indiscriminate responses, up to 80% of iRBCs were cleared per day by the host one to three days after the peak of infection (Fig 5a). The activation of targeted iRBC clearance was estimated to be independent of the initial infection dose (Fig 4a) while we found that the half-life of this response decreased with infection dose (Fig 4b). This result is consistent with the faster waning of parasite clearance with dose observed by Metcalf et al. [37] who speculate mechanisms including enhanced antigen escape, reduced immune memory due to low RBC availability, depleted immune effectors and downregulation by the host or parasites. Without this dose-dependent effect, our sensitivity analysis demonstrated that the peak targeted iRBC clearance would occur at high doses, where higher iRBC density triggers an elevated response (S1f Appendix).

To assess the relative importance of dose-dependence in the host responses, we examined the sensitivity of the model fit to whether host response parameters were dose-dependent (by setting , where k = {ψ1, ϕ1}, k = {ψ2, ϕ2}, for the two responses, respectively). We found that the dynamics of RBCs were most sensitive to dose-dependence in general RBC clearance: goodness of model fit declined by 2.6-fold and 3.4-fold, respectively, when dose-dependence in activation strength or activity half-life was ignored (Fig 6). Interestingly, even though the targeted response clears more iRBCs than the general response (Fig 5a & 5b), we found that iRBC dynamics were also overwhelmingly more sensitive to dose-dependence in general RBC clearance (Fig 6): goodness of model fit declined by 3.6-fold and 7-fold, respectively, when dose-dependence in either activation strength or activity half-life (Fig 6) was ignored. While the mechanism through which the host clears RBCs indiscriminately remains an open question, the functional importance of general clearance of RBCs is apparent from our work here and the work of others [41]. Furthermore, the distinct patterns of dose-dependence we found (Fig 4a & 4b) suggest that there is separate host machinery for specifically clearing iRBCs and clearing RBCs indiscriminately.

[Figure omitted. See PDF.]

Fig 6. RBC and iRBC dynamics are most sensitive to dose-dependence in general RBC clearance.

Plotted is the sensitivity of RBC and iRBC dynamics to dose-dependence of a) activation strength and b) activity half-life in general RBC clearance and targeted iRBC clearance (green), calculated as the sum of absolute standardised residuals (see S1e Appendix for details) relative to the full model (with parameters defined by Eq 6).

https://doi.org/10.1371/journal.pcbi.1008211.g006

Previously, dose-dependence in infection dynamics has been attributed to ineffective clearance of iRBCs at high doses due to handling time in innate immune effectors targeting iRBCs [37]. Our findings offer an alternative explanation that dose promotes some aspects of some response while debilitating others. Of largest effects, we found that dose increases the half-life, but reduces activation strength, of general RBC clearance (Figs 4 & 6). This finding likely reflects complex mechanisms through which the initial infection dose impacts different aspects of host response machinery [24–28, 80]. Regardless of the mode of defence involved (i.e., indiscriminate or targeted iRBC clearance) and mechanisms (i.e., handling time, damage to host machinery, or immune evasion), the two distinct modelling frameworks demonstrate that malaria parasites are at an advantage at high initial infection dose due to less efficient host responses to clear iRBCs. We did not find significant evidence that upregulation in erythropoiesis or parasite burst size depends on the initial iRBC density (Fig 4c & 4d). Therefore, our findings suggest that dose-dependent variation in infection dynamics observed across initial infection dose treatments in Timms et al. [15] was driven by the interaction between the initial infection dose and host responses, but not plasticity in malaria parasites injected at different doses as predicted by Mideo et al. [76].

Individual variation

We explicitly modelled individual variation in each fitted parameter among the inbred mice using hyperparameters, σu,k and σv,k. These sources of variation—analogous to subject-level random effects in regression analyses—capture unobserved heterogeneity among individuals in a sample, independent of the experimental manipulation in the initial infection dose. We found evidence of individual variation in every host and parasite parameter of the model describing malaria infection ecology (Fig 7a; S1b Appendix). There was no evidence of moderate or stronger correlation (r > 0.3) among individual variation in parameter values (S1g Appendix), suggesting that there are no clear trade-offs nor facilitation among different arms of host responses.

[Figure omitted. See PDF.]

Fig 7. Individual variation in parameters of host responses and parasite growth impacts quantitative infection outcomes.

a) Coefficient of variation among individuals in specific model parameters and initial conditions. Variation in the dose-independent (intercept) term, σu,k and dose-dependent (slope) term σv,k. is indicated in black and grey, respectively. b) Correlation coefficients between individual-level estimates of parameters or initial conditions, and the infection outcome measured as the maximum and minimum densities of iRBCs (blue) and RBCs (red), respectively. The error bars indicate the 95% confidence interval. Statistical significance is indicated when the confidence interval does not intersect 0 (dashed line).

https://doi.org/10.1371/journal.pcbi.1008211.g007

To understand the functional importance of individual variation in estimated model parameters, as well as random variation in the reported initial RBC density, host weight, and initial infection dose, we computed the correlation coefficient, r, between individual variation and two key infection outcomes: parasite load (maximum iRBC density) and anaemia severity (minimum RBC density). We identified two individual-level correlates of the maximum iRBC density (Fig 7b). First, we found that better suppression of the peak parasite load was associated with mice that activated the response to clear iRBCs more strongly than average through a targeted mechanism (Fig 7b). While it is not clear what underlies this within-strain variation, our finding, nonetheless, demonstrates the functional importance of non-genetic heterogeneity in immune responses, drawing parallels to invertebrate systems in which intrinsic within-clone variation (e.g., differences in sizes at birth and molecular mechanisms of immune responses) are thought to impact susceptibility to infection [81–83]. With the rodent malaria system, it may be possible to identify the mechanistic causal agents, for example, by characterising within-strain variation in immune effector expression prior to and throughout the course of infection. We also found that variation in the initial infection dose within dose treatments correlated positively with the peak parasite load (Fig 7b). In the context of Timms et al.’s [15] experiments, variation in the initial infection dose may be attributable, for example, to random experimental variability in sampling by a syringe and in the parasites’ ability to reach blood vessels following intraperitoneal injection. While neither of these sources of random variation are relevant outside the lab, the observation that higher doses—both among and within dose treatments—increase parasite burden [15] (Figs 1 & 7b) highlights the need to better understand any causes of variation in pre-blood-stage parasite densities in natural malaria infections, both within and among Plasmodium species [58]. One such source of variation is pre-blood-stage host immunity that develops in response to exposure to the parasite stage injected by mosquitoes (sporozoites) and has been shown effective in reducing the number of liver-stage parasites from which blood-stage merozoites originate [84].

There were also three significant individual-level correlates of anaemia severity (Fig 7b). First, we found a significant negative effect of burst size meaning that infections initiated with a parasite population that happens to proliferate and exhaust RBCs faster than average caused more severe anaemia (Fig 7b). In addition, mice that were heavier or had more RBCs before infection suffered less severe anaemia (Fig 7), suggesting that general host vigour in the absence of infection is an indicator of host health under malaria infection. Because these measurements of host vigour did not correlate with the peak parasite load, host vigour can be interpreted as an indicator of tolerance (i.e., a host’s ability to minimise anaemia), rather than a resistance (i.e., a host’s ability to minimise parasite burden) mechanism [85].

In general, we found that the magnitude of individual variation, estimated here as the coefficient of variation, did not coincide with whether or not a particular trait impacted an infection outcome. For example, we estimated large individual-level variation in the half-life of host response for targeted clearance of iRBCs (Fig 7a), yet infection outcomes were insensitive to this variation (Fig 7b). Conversely, we detected a small coefficient of variation in burst size (Fig 7a), but this small variation showed a marked impact on the severity of anaemia (Fig 7b). The misalignment between the magnitude and impact of individual variation poses a challenge from a clinical perspective because it is not clear whether those parameters can be estimated with sufficient precision from patient data. Finally, despite the fact that the dynamics of RBCs and iRBCs are ecologically coupled (i.e., malaria parasites are consumers of RBCs), we found that the impact of individual variation on anaemia severity was not coupled to that of parasite load. These findings indicate that resistance and tolerance are likely uncorrelated at the individual level in the rodent malaria system.

Conclusion

We examined drivers of dose-dependent malaria parasite load and the severity of malaria-induced anaemia. We also shed light on the role of unobserved heterogeneity in producing diverse infection outcomes, identifying sources of subtle variation beneath an experimental treatment. More often than not, infection experiments are structured in multiple levels of sampling units (e.g., host, parasite genotypes, presence of coinfection, and drug treatments) with many replicates within treatment groups. Observations from such experiments contain multiple sources of variation whose effects on infection outcomes are difficult to disentangle directly. Our study demonstrates that the combination of dynamical within-host model and a Bayesian approach is a powerful tool for causal inference of infection outcome variation.

Supporting information

S1 Appendix.

a) Estimation of the initial infection dose, b) Graphical summary of prior and posterior distributions, c) Mice excluded from model fitting, d) Computer programmes, e) Assessment of model fit: standardised residuals, f) Sensitivity of targeted iRBC clearance to dose-dependent half-life, and g) Correlations of individual variation.

https://doi.org/10.1371/journal.pcbi.1008211.s001

(PDF)

Acknowledgments

The project was made possible thanks to the experimental data collected and published by Timms and coworkers in the Read Group. We are also grateful to the community of Stan developers and users, who were extremely generous with their time on The Stan Forums https://discourse.mc-stan.org/. We thank Troy Day, Devin Kirk, Katie O’Dwyer and members of the Mideo group at the University of Toronto and the Schneider Group at Stanford University for helpful discussion. This work also benefited from discussions at an RCN-IDEAS-funded workshop on “Immunity Across Scales”, organised by Clay Cressler in 2018.

Citation: Kamiya T, Greischar MA, Schneider DS, Mideo N (2020) Uncovering drivers of dose-dependence and individual variation in malaria infection outcomes. PLoS Comput Biol 16(10): e1008211. https://doi.org/10.1371/journal.pcbi.1008211

References

1. Galatas B, Bassat Q, Mayor A. Malaria parasites in the asymptomatic: looking for the hay in the haystack. Trends in Parasitology. 2015;32(4):296–308. pmid:26708404

2. Luzzatto L. Sickle cell anaemia and malaria. Mediterranean Journal of Hematology and Infectious Diseases. 2012;4(1). pmid:23170194

3. Shankar AH. Nutritional modulation of malaria morbidity and mortality. Journal of Infectious Diseases. 2000;182(Supplement_1):S37–S53. pmid:10944483

4. Haydon DT, Matthews L, Timms R, Colegrave N. Top–down or bottom–up regulation of intra–host blood–stage malaria: do malaria parasites most resemble the dynamics of prey or predator? Proceedings of the Royal Society B: Biological Sciences. 2003;270(1512):289–298. pmid:12614579

5. Weatherall D. Genetic variation and susceptibility to infection: the red cell and malaria. British Journal of Haematology. 2008;141(3):276–286. pmid:18410566

6. Hochman S, Kim K. The impact of HIV and malaria coinfection: what is known and suggested venues for further study. Interdisciplinary Perspectives on Infectious Diseases. 2009;2009. pmid:19680452

7. Lazzaro BP, Little TJ. Immunity in a variable world. Philosophical Transactions of the Royal Society B: Biological Sciences. 2008;364(1513):15–26. pmid:18926975

8. Schmid-Hempel P, Frank SA. Pathogenesis, virulence, and infective dose. PLOS Pathogens. 2007;3(10). pmid:17967057

9. Leggett HC, Cornwallis CK, West SA. Mechanisms of pathogenesis, infective dose and virulence in human parasites. PLOS Pathogens. 2012;8(2). pmid:22359500

10. Nie P, Kennedy C. Infection dynamics of larval Bothriocephalus claviceps in Cyclops vicinus. Parasitology. 1993;106(5):503–509.

11. Ashworth S, Kennedy C, Blanc G. Density-dependent effects of Anguillicola crassus (Nematoda) within and on its copepod intermediate hosts. Parasitology. 1996;113(3):303–309. pmid:8811853

12. Michael E, Bundy D. Density dependence in establishment, growth and worm fecundity in intestinal helminthiasis: the population biology of Trichuris muris (Nematoda) infection in CBA/Ca mice. Parasitology. 1989;98(3):451–458. pmid:2771451

13. Mackinnon MJ, Read AF. Genetic relationships between parasite virulence and transmission in the rodent malaria Plasmodium chabaudi. Evolution. 1999;53(3):689–703. pmid:28565637

14. Ebert D, Zschokke-Rohringer CD, Carius HJ. Dose effects and density-dependent regulation of two microparasites of Daphnia magna. Oecologia. 2000;122(2):200–209. pmid:28308373

15. Timms R, Colegrave N, Chan B, Read A. The effect of parasite dose on disease severity in the rodent malaria Plasmodium chabaudi. Parasitology. 2001;123(1):1–11. pmid:11467777

16. Tan GK, Ng JK, Trasti SL, Schul W, Yip G, Alonso S. A non mouse-adapted dengue virus strain as a new model of severe dengue infection in AG129 mice. PLOS Neglected Tropical Diseases. 2010;4(4). pmid:20436920

17. Brunner JL, Richards K, Collins JP. Dose and host characteristics influence virulence of ranavirus infections. Oecologia. 2005;144(3):399–406. pmid:15891818

18. McLean AR, Bostock CJ. Scrapie infections initiated at varying doses: an analysis of 117 titration experiments. Philosophical Transactions of the Royal Society of London Series B: Biological Sciences. 2000;355(1400):1043–1050. pmid:11186305

19. Hughes G, Kitching R, Woolhouse M. Dose-dependent responses of sheep inoculated intranasally with a type O foot-and-mouth disease virus. Journal of Comparative Pathology. 2002;127(1):22–29. pmid:12354542

20. Li Y, Handel A. Modeling inoculum dose dependent patterns of acute virus infections. Journal of Theoretical Biology. 2014;347:63–73. pmid:24440713

21. Handel A, Li Y, McKay B, Pawelek KA, Zarnitsyna V, Antia R. Exploring the impact of inoculum dose on host immunity and morbidity to inform model-based vaccine design. PLOS Computational Biology. 2018;14(10):e1006505. pmid:30273336

22. Schmidt C, Schneble N, Müller J, Bauer R, Perino A, Marone R, et al. Phosphoinositide 3-kinase γ mediates microglial phagocytosis via lipid kinase-independent control of cAMP. Neuroscience. 2013;233:44–53. pmid:23276671

23. Morris MC, Gilliam EA, Button J, Li L. Dynamic modulation of innate immune response by varying dosages of lipopolysaccharide (LPS) in human monocytic cells. Journal of Biological Chemistry. 2014;289(31):21584–21590. pmid:24970893

24. Schmidt C, Schneble N, Wetzker R. The fifth dimension of innate immunity. Journal of Cell Communication and Signaling. 2014;8(4):363–367. pmid:25278167

25. Hatta Y, Hershberger K, Shinya K, Proll SC, Dubielzig RR, Hatta M, et al. Viral replication rate regulates clinical outcome and CD8 T cell responses during highly pathogenic H5N1 influenza virus infection in mice. PLOS Pathogens. 2010;6(10):e1001139. pmid:20949022

26. Marois I, Cloutier A, Garneau É, Richter MV. Initial infectious dose dictates the innate, adaptive, and memory responses to influenza in the respiratory tract. Journal of Leukocyte Biology. 2012;92(1):107–121. pmid:22504848

27. Segueni N, Tritto E, Bourigault ML, Rose S, Erard F, Le Bert M, et al. Controlled Mycobacterium tuberculosis infection in mice under treatment with anti-IL-17A or IL-17F antibodies, in contrast to TNFα neutralization. Scientific Reports. 2016;6:36923. pmid:27853279

28. Pagan AJ, Levitte S, Berg RD, Hernandez L, Zimmerman J, Tobin DM, et al. mTOR deficiency reveals an immunological trade-off in innate resistance to mycobacterial infection in vivo. Journal of Immunology. 2016.

29. Rhodes SJ, Guedj J, Fletcher HA, Lindenstrøm T, Scriba TJ, Evans TG, et al. Using vaccine Immunostimulation/Immunodynamic modelling methods to inform vaccine dose decision-making. NPJ Vaccines. 2018;3(1):1–7. pmid:30245860

30. Rhodes SJ, Knight GM, Kirschner DE, White RG, Evans TG. Dose finding for new vaccines: The role for immunostimulation/immunodynamic modelling. Journal of theoretical biology. 2019;465:51–55. pmid:30639297

31. Afrough S, Rhodes S, Evans T, White R, Benest J. Immunologic dose-response to adenovirus-vectored vaccines in animals and humans: a systematic review of dose-response studies of replication incompetent adenoviral vaccine vectors when given via an intramuscular or subcutaneous route. Vaccines. 2020;8(1):131. pmid:32192058

32. Garnham PCC. Malaria parasites and other haemosporidia. Blackwell Scientific; 1966.

33. Rosenberg R, Wirtz RA, Schneider I, Burge R. An estimation of the number of malaria sporozoites ejected by a feeding mosquito. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1990;84(2):209–212. pmid:2202101

34. Beier J, Onyango F, Koros J, Ramadhan M, Ogwang R, Wirtz R, et al. Quantitation of malaria sporozoites transmitted in vitro during salivation by wild Afrotropical Anopheles. Medical and Veterinary Entomology. 1991;5(1):71–79. pmid:1768903

35. Lines J, Armstrong J. For a few parasites more: inoculum size, vector control and strain-specific immunity to malaria. Parasitology Today. 1992;8(11):381–383. pmid:15463547

36. Beier J. Malaria sporozoites: survival, transmission and disease control. Parasitology Today. 1993;9(6):210–215. pmid:15463756

37. Metcalf C, Graham A, Huijben S, Barclay V, Long G, Grenfell B, et al. Partitioning regulatory mechanisms of within-host malaria dynamics using the effective propagation number. Science. 2011;333(6045):984–988. pmid:21852493

38. Kochin BF, Yates AJ, De Roode JC, Antia R. On the control of acute rodent malaria infections by innate immunity. PLOS One. 2010;5(5):e10444. pmid:20463903

39. Miller MR, Råberg L, Read AF, Savill NJ. Quantitative analysis of immune response and erythropoiesis during rodent malarial infection. PLOS Computational Biology. 2010;6(9):e1000946. pmid:20941388

40. Mideo N, Savill NJ, Chadwick W, Schneider P, Read AF, Day T, et al. Causes of variation in malaria infection dynamics: insights from theory and data. American Naturalist. 2011;178(6):174–188. pmid:22089879

41. Wale N, Jones MJ, Sim DG, Read AF, King AA. The contribution of host cell-directed vs. parasite-directed immunity to the disease and dynamics of malaria infections. Proceedings of the National Academy of Sciences. 2019;116(44):22386–22392. pmid:31615885

42. Mideo N, Day T, Read AF. Modelling malaria pathogenesis. Cellular Microbiology. 2008;10(10):1947–1955. pmid:18647174

43. Stevenson MM, Riley EM. Innate immunity to malaria. Nature Reviews Immunology. 2004;4(3):169. pmid:15039754

44. Gazzinelli RT, Kalantari P, Fitzgerald KA, Golenbock DT. Innate sensing of malaria parasites. Nature Reviews Immunology. 2014;14(11):744. pmid:25324127

45. Schneider DS. Tracing personalized health curves during infections. PLOS Biology. 2011;9(9):e1001158. pmid:21957398

46. Stephens R, Culleton RL, Lamb TJ. The contribution of Plasmodium chabaudi to our understanding of malaria. Trends in Parasitology. 2012;28(2):73–82. pmid:22100995

47. Timms R, Colegrave N, Chan B, Read A. The effect of parasite dose on disease severity in the rodent malaria Plasmodium chabaudi; 2001. http://dx.doi.org/10.5061/dryad.stqjq2c1k.

48. Safeukui I, Correas JM, Brousse V, Hirt D, Deplaine G, Mulé S, et al. Retention of Plasmodium falciparum ring-infected erythrocytes in the slow, open microcirculation of the human spleen. Blood. 2008;112(6):2520–2528. pmid:18579796

49. White NJ. Malaria parasite clearance. Malaria Journal. 2017;16(1):88. pmid:28231817

50. Price RN, Simpson JA, Nosten F, Luxemburger C, Hkirjaroen L, ter Kuile F, et al. Factors contributing to anemia after uncomplicated falciparum malaria. The American Journal of Tropical Medicine and Hygiene. 2001;65(5):614–622. pmid:11716124

51. Castro-Gomes T, Mourão LC, Melo GC, Monteiro WM, Lacerda MV, Braga ÉM. Potential immune mechanisms associated with anemia in Plasmodium vivax malaria: a puzzling question. Infection and Immunity. 2014;p. IAI–01972. pmid:25092911

52. Metcalf C, Long G, Mideo N, Forester J, Bjørnstad O, Graham A. Revealing mechanisms underlying variation in malaria virulence: effective propagation and host control of uninfected red blood cell supply. Journal of The Royal Society Interface. 2012;p. rsif20120340. pmid:22718989

53. Chang KH, Tam M, Stevenson MM. Modulation of the course and outcome of blood-stage malaria by erythropoietin-induced reticulocytosis. Journal of Infectious Diseases. 2004;189(4):735–743. pmid:14767829

54. Bunn HF. Erythropoietin. Cold Spring Harbor Perspectives in Medicine. 2013;3(3):a011619. pmid:23457296

55. Greischar MA, Mideo N, Read AF, Bjørnstad ON. Quantifying transmission investment in malaria parasites. PLOS Computational Biology. 2016;12(2). pmid:26890485

56. Gravenor MB, Lloyd AL, Kremsner PG, Missinou MA, English M, Marsh K, et al. A model for estimating total parasite load in falciparum malaria patients. Journal of Theoretical Biology. 2002;217(2):137–148. pmid:12202108

57. Smith H. Distributed delay equations and the linear chain trick. In: An Introduction to Delay Differential Equations with Applications to the Life Sciences. Springer; 2011. p. 119–130.

58. Paul R, Ariey F, Robert V. The evolutionary ecology of Plasmodium. Ecology Letters. 2003;6(9):866–880.

59. Anderson R, May R, Gupta S. Non-linear phenomena in host-parasite interactions. Parasitology. 1989;99(S1):S59–S79. pmid:2682486

60. Hellriegel B. Modelling the immune response to malaria with ecological concepts: short-term behaviour against long-term equilibrium. Proceedings of the Royal Society of London Series B: Biological Sciences. 1992;250(1329):249–256. pmid:1362993

61. Saul A. Models for the in-host dynamics of malaria revisited: errors in some basic models lead to large over-estimates of growth rates. Parasitology. 1998;117(5):405–407. pmid:9836304

62. White N, Chapman D, Watt G. The effects of multiplication and synchronicity on the vascular distribution of parasites in falciparum malaria. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1992;86(6):590–597. pmid:1287908

63. Molineaux L, Diebner H, Eichner M, Collins W, Jeffery G, Dietz K. Plasmodium falciparum parasitaemia described by a new mathematical model. Parasitology. 2001;122(4):379–391. pmid:11315171

64. Greischar MA, Read AF, Bjørnstad ON. Synchrony in malaria infections: how intensifying within-host competition can be adaptive. American Naturalist. 2013;183(2):E36–E49. pmid:24464205

65. Crooks L. Problems with continuous-time malaria models in describing gametocytogenesis. Parasitology. 2008;135(8):881–896. pmid:18598574

66. McAlister RO. Time-dependent loss of invasive ability of Plasmodium berghei merozoites in vitro. Journal of Parasitology. 1977;p. 455–463. pmid:325188

67. McElreath R. Statistical rethinking: A Bayesian course with examples in R and Stan. Chapman and Hall/CRC; 2018.

68. Gelman A, Stern HS, Carlin JB, Dunson DB, Vehtari A, Rubin DB. Bayesian data analysis. Chapman and Hall/CRC; 2013.

69. Gelman A. Bayesian Statistical Pragmatism. Statistical Science. 2011;26(1):10–11.

70. Mugglin AS, Cressie N, Gemmell I. Hierarchical statistical modelling of influenza epidemic dynamics in space and time. Statistics in Medicine. 2002;21(18):2703–2721. pmid:12228886

71. Cressie N, Calder CA, Clark JS, Hoef JMV, Wikle CK. Accounting for uncertainty in ecological analysis: the strengths and limitations of hierarchical statistical modeling. Ecological Applications. 2009;19(3):553–570. pmid:19425416

72. Martin JG, Nussey DH, Wilson AJ, Réale D. Measuring individual differences in reaction norms in field and experimental studies: a power analysis of random regression models. Methods in Ecology and Evolution. 2011;2(4):362–374.

73. Dingemanse NJ, Dochtermann NA. Quantifying individual variation in behaviour: mixed-effect modelling approaches. Journal of Animal Ecology. 2013;82(1):39–54. pmid:23171297

74. Gelman A. Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper). Bayesian Analysis. 2006;1(3):515–534.

75. Sorensen T, Hohenstein S, Vasishth S. Bayesian linear mixed models using Stan: A tutorial for psychologists, linguists, and cognitive scientists. Quantitative Methods for Psychology. 2016;12(3):175–200.

76. Mideo N, Barclay VC, Chan BH, Savill NJ, Read AF, Day T. Understanding and predicting strain-specific patterns of pathogenesis in the rodent malaria Plasmodium chabaudi. American Naturalist. 2008;172(5):E214–E238. pmid:18834302

77. Carpenter B, Gelman A, Hoffman MD, Lee D, Goodrich B, Betancourt M, et al. Stan: A probabilistic programming language. Journal of Statistical Software. 2017;76(1).

78. Stan Development Team. RStan: the R interface to Stan, Version 2.18.2.; 2019. http://mc-stan.org.

79. Stan Development Team. Stan Modeling Language User’s Guide and Reference Manual, Version 2.18.0; 2018. Available from: http://mc-stan.org/.

80. Chong KKL, Tay WH, Janela B, Yong AMH, Liew TH, Madden L, et al. Enterococcus faecalis modulates immune activation and slows healing during wound infection. Journal of Infectious Diseases. 2017;216(12):1644–1654. pmid:29045678

81. Dwyer G, Elkinton JS, Buonaccorsi JP. Host heterogeneity in susceptibility and disease dynamics: tests of a mathematical model. American Naturalist. 1997;150(6):685–707. pmid:18811331

82. Watson FL, Püttmann-Holgado R, Thomas F, Lamar DL, Hughes M, Kondo M, et al. Extensive diversity of Ig-superfamily proteins in the immune system of insects. Science. 2005;309(5742):1874–1878. pmid:16109846

83. Ben-Ami F, Regoes RR, Ebert D. A quantitative test of the relationship between parasite dose and infection probability across different host–parasite combinations. Proceedings of the Royal Society B: Biological Sciences. 2008;275(1636):853–859. pmid:18198145

84. Nahrendorf W, Spence PJ, Tumwine I, Levy P, Jarra W, Sauerwein RW, et al. Blood-stage immunity to Plasmodium chabaudi malaria following chemoprophylaxis and sporozoite immunization. eLife. 2015;4:e05165. pmid:25714922

85. Schneider DS, Ayres JS. Two ways to survive infection: what resistance and tolerance can teach us about treating infectious diseases. Nature Reviews Immunology. 2008;8(11):889. pmid:18927577