Estimation of heterogeneity in malaria transmission by stochastic modelling of apparent deviations from mass action kinetics
- Thomas A Smith^{1}Email author
https://doi.org/10.1186/1475-2875-7-12
© Smith; licensee BioMed Central Ltd. 2008
Received: 25 July 2007
Accepted: 11 January 2008
Published: 11 January 2008
Abstract
Background
Quantifying heterogeneity in malaria transmission is a prerequisite for accurate predictive mathematical models, but the variance in field measurements of exposure overestimates true micro-heterogeneity because it is inflated to an uncertain extent by sampling variation. Descriptions of field data also suggest that the rate of Plasmodium falciparum infection is not proportional to the intensity of challenge by infectious vectors. This appears to violate the principle of mass action that is implied by malaria biology. Micro-heterogeneity may be the reason for this anomaly. It is proposed that the level of micro-heterogeneity can be estimated from statistical models that estimate the amount of variation in transmission most compatible with a mass-action model for the relationship of infection to exposure.
Methods
The relationship between the entomological inoculation rate (EIR) for falciparum malaria and infection risk was reanalysed using published data for cohorts of children in Saradidi (western Kenya). Infection risk was treated as binomially distributed, and measurement-error (Poisson and negative binomial) models were considered for the EIR. Models were fitted using Bayesian Markov chain Monte Carlo algorithms and model fit compared for models that assume either mass-action kinetics, facilitation, competition or saturation of the infection process with increasing EIR.
Results
The proportion of inocula that resulted in infection in Saradidi was inversely related to the measured intensity of challenge. Models of facilitation showed, therefore, a poor fit to the data. When sampling error in the EIR was neglected, either competition or saturation needed to be incorporated in the model in order to give a good fit. Negative binomial models for the error in exposure could achieve a comparable fit while incorporating the more parsimonious and biologically plausible mass action assumption. Models that assume negative binomial micro-heterogeneity predict lower incidence of infection at a given average exposure than do those assuming exposure to be uniform. The negative binomial model moreover provides an estimate of the variance of the within-cohort distribution of the EIR and hence of within cohort heterogeneity in exposure.
Conclusion
Apparent deviations from mass action kinetics in parasite transmission can arise from spatial and temporal heterogeneity in the inoculation rate, and from imprecision in its measurement. For parasites like P. falciparum, where there is no plausible biological rationale for deviations from mass action, this provides a strategy for estimating true levels of heterogeneity, since if mass-action is assumed, the within-population variance in exposure becomes identifiable in cohort studies relating infection to transmission intensity. Statistical analyses relating infection to exposure thus provide a valid general approach for estimating heterogeneity in transmission but only when they incorporate mass action kinetics and shrinkage estimates of exposure. Such analyses make it possible to include realistic levels of heterogeneity in dynamic models that predict the impact of control measures on transmission intensity.
Background
Heterogeneity in transmission rates of parasitic infections is a major source of uncertainty in predictions of the impact of control measures. In general, heterogeneity decreases the effectiveness of untargeted interventions, and increases the effectiveness of those that are targeted to high transmission groups [1]. In order to make quantitative predictions of the impact of control measures we need models with appropriate allowance for heterogeneity.
For most infectious diseases it is difficult, if not impossible, to directly measure heterogeneity. For instance, Plasmodium falciparum malaria transmission is measured by the Entomological Inoculation Rate (EIR), which is the product of the human biting rate (ma) (mosquito bites per person per night) and the sporozoite rate (the proportion of mosquitoes carrying the infectious sporozoite stage in their salivary glands, s) [2]. Both these quantities can be measured only with difficulty and with considerable sampling error. In particular, ma is assessed either by human landing collections, which rely on the assumption that young men deliberately expose themselves to biting have similar exposure to the general population nearby, or more indirectly, by using mosquito traps.
Many studies have recorded human landing rates in different malaria endemic communities, and it is possible to partition the variation observed into spatial, temporal, and observer specific components, and to estimate the effects of housing construction, or other interventions. However there is always a substantial residual variance that cannot be accounted for by any measured factor and because any given individual's exposure can only be assessed using one method at any one time, it is unclear how much of the variation between measurements reflects real heterogeneity in exposure and how much reflects the variability inherent in the sampling methods. A frequently cited estimate is that 20% of the population receive 80% of the exposure [1], but this figure is based on an empirical summary of light trap data from which the true heterogeneity cannot be estimated.
It might be the case that there is less real heterogeneity than the heterogeneity in the measurements, because the variation in measurements is inflated by differences in the way the measuring tools are applied on different occasions (random measurement error) and by sampling variation. Conversely, if there are many replications of the measurements the variance of the true exposures might be much greater than that of the averaged measurements.
One approach to estimating the degree of heterogeneity that has been attempted is to fit models to age-prevalence data as a function of EIR, incorporating a parameter for heterogeneity[3], however this entails fitting rather speculative models for immunity (the fact that such a model fits well does not mean that it is correct[4]). Moreover, it is hard to see how this approach can incorporate the effects of unknown imprecision in EIR measurements.
Analysis of the shape of the relationship between infection incidence and measured exposure provides a strategy for estimating the true level of heterogeneity. Most models of infectious diseases assume that transmission occurs at a rate directly proportional to the number or density of both susceptible and of infectious hosts. This mass action assumption corresponds to the established biology of infection in homogeneous populations, but where transmission is heterogeneous the population average force of infection is no longer proportional to the average density of infectious hosts, or (in the case of vector borne diseases) to the vectorial capacity [5] and a lower proportion of individuals is infected at any given exposure than is predicted by mass action. If the average exposure and the susceptibility of the hosts are known, then the reduction in the force of infection below the mass action expectation can be used to estimate the amount of micro-heterogeneity. A complication in the case of P. falciparum is that susceptibility in the field is uncertain, and the imprecision of exposure estimates also needs to be taken into consideration.
This article applies a statistical approach to estimate the level of heterogeneity in P. falciparum transmission from exposure-infection relationships in field data from Kenya. The method uses negative binomial models to simultaneously estimate the level of heterogeneity and the susceptibility, incorporating exposure estimates adjusted for effects of sampling variation and assuming underlying mass action kinetics. The estimate of the degree of micro-heterogeneity in transmission provided by this model gives an empirical basis for the level of heterogeneity to incorporate in predictive models. This model is contrasted with other models fitted to the same data that assume either facilitation or competition between different inocula, implying dose-dependence in the infection process [6].
Methods
Field data
The study of Beier and others [7] in Saradidi, Kenya, represents the most informative dataset available on the relationship between P. falciparum exposure and infection.
In this study, P. falciparum parasites were cleared using sulphadoxine-pyrimethamine in 21 cohorts children aged between six months and six years. Follow-up of each of the 21 cohorts started at four week intervals, and infection incidence in each cohort was recorded subsequently in four fortnightly parasitological surveys. Each of these four inter-survey intervals, nested within the cohort, defines one sampling unit, i = 1,2....84, characterized by a number of children at risk, a proportion, h_{ i }, who became infected, and a true average EIR, e_{t(i)}, where t(i) indicates the fortnight (out of the 44 possible) for which the sampling unit was studied.
The entomological data for estimating each value of e_{t(i)}were obtained using weekly all-night human biting collections made inside six houses by six teams [8]. Each of the 44 estimates, z_{t(i)}, corresponding to different values of t(i), thus used a total of 12 human bait collections from which both ma and s were calculated and the published EIR values [7] represent the sum of the sporozoite positive mosquitoes captured, scaled to give values in infectious bites per fortnight.
The dataset as published by Beier and others[7], includes the numbers of children infected, and those at risk, from which the proportion of children infected (h_{ i }) can be estimated for each sampling unit, together with point estimates of e_{t(i)}for each fortnight. The present paper reports a series of statistical models for the relationship between h_{ i }and estimates of e_{t(i)}made using different models for the distribution of inoculations in the population and mprecision in the EIR measurements.
Patterns of seasonality in the data were reported by Beier and others [8] and the original publication was followed in relating recorded infection incidence to concurrently measured exposure. No details of the sampled houses, of day-to-day variation within each two week period, or of inter observer variation in biting rates are available.
Models for the relationship of infection incidence to EIR
Children attract fewer mosquitoes than adults [9, 10], and hence e_{t(i)}is higher than the average number of infectious bites received by children. The Saradidi dataset does not provide any disaggregation of either infection or exposure data by age, and so the extent of this effect must be inferred from other data. In previous analyses[11], average weights of children were used in an East African population to obtain a scale factor, w = 0.302, to be applied to the children in the Saradidi study to allow for this effect, so that the average age adjusted EIR value, e_{at(i)}, is:
e_{at(i)}= we_{t(i)} (1)
Let the random variable x_{ ij }be the actual number of infectious bites received by individual j in unit i any one time period, and let S_{ i }be the probability that an inoculation in sampling unit i survives to give a patent blood-stage infection. Computationally, it is more straightforward to treat the effect of host size as equivalent to a reduction in this survival probability, so that the number of inoculations experienced by children is nominally also x_{ ij }, sampled from the same distribution (with mean e_{t(i)}) as for adults, while the effective survival of the inoculum in children is::
S_{ ai }= wS_{ i } (2)
Various models were consider that differ in terms of the functional form of S_{ i }as follows:
a. Mass Action
The simplest stochastic model for the infection process assumes that S_{ i }is a constant, S, and hence S_{ ai }= wS.
b. Facilitation
The model for human infection thus has two parameters, ${E}_{f}^{*}$ and γ.
c. Saturation
Descriptive analyses of the Saradidi data indicated that the survival probability of the inoculum appears to decrease with increasing transmission intensity [7], suggesting that there is density dependence in the infection probability. In principle, this might arise because of saturation at high inoculation rates, as assumed in the malaria model of the Garki project [12]. [13] adapted the function used in the Garki model and fitted it to the Saradidi data assuming:
E(h_{ ij }) = g(1 - exp(-βe_{at(i)})) (5)
where the two parameters g and β are both measures of susceptibility. This model was generalized by considering uncertainty in e_{at(i)}as well as in the proportion of children infected. This leads to different expressions for S_{ ai }depending on the distribution of e_{at(i)}or equivalently e_{t(i)}(see below).
d. Competition
The parameter S_{∞} is then the lower limit attained by S_{ ai }as the inoculation rate becomes large, and ${E}_{c}^{*}$ is the value of e_{at(i)}at which half the reduction in S_{ ai }is achieved.
Models for imprecision in the estimation of the EIR
The assumptions that we make about the distribution of imprecision in the measurements of the EIR are reflected in the distribution used for the number of bites on an individual, x_{ ij }, with consequences for the probability that a child is infected.
a. No shrinkage of estimates of e_{t(i)}
In the first set (models 1–4) following previous analyses of the Saradidi data, we treat the true average EIR for the sampling unit, e_{t(i)}, as given by the data, ignoring both variation between individuals in exposure and imprecision in the determination of the average, so that both the individual exposures x_{ ij }, and the average EIR, e_{t(i)}, are assumed equal to the measured value, z_{t(i)}.
Correspondingly the force of infection in children, equal to the intensity of the infection process, is assumed known, as S_{ ai }z_{t(i)}, and treating the number of these infections per individual as Poisson distributed the probability that a child receives at least one infection is:
E(h_{ ij }) = 1 - exp (-S_{ ai }z_{t(i)}) (7)
b. Poisson variation in estimates of e_{t(i)}
In the second set (models 5–8), recognising that z_{t(i)}is obtained by counting sporozoite positive mosquitoes, the count of infections recorded in two weeks per mosquito collector is assumed to be a Poisson distributed random variable with expectation equal to the true average for the sampling unit:
z_{t(i)}~ Poisson(e_{t(i)}) (8)
The expected value of z_{t(i)}remains equal to e_{t(i)}, but there is random Poisson variation around this expectation.
We again assume that the number of infections per individual, x_{ ij }, is Poisson distributed, but the mean of the Poisson is now the unknown quantity e_{t(i)}, rather than the observed number of inoculations z_{t(i)}The probability that a child receives at least one infection is therefore:
E(h_{ ij }) = 1 - exp (-S_{ ai }e_{t(i)}) (9)
e_{t(i)}is now a vector of parameters to be estimated conditional on the observed values of z_{ i }
c. Negative binomial variation in estimates of e_{t(i)}
Finally (models 9–12), we consider z_{t(i)}as negative binomially distributed with mean e_{t(i)}and overdispersion parameter κ i.e.:
For the model with saturation (model 11), E(h_{ ij }), is given explicitly by equation (5).
Model fitting
Treating e_{t(i)}as known (models 1–4) leads to non-linear models with binomial errors in the infection probability, for which maximum likelihood methods could easily be applied. The remaining models represent hierarchical random effects models, which can be fitted most easily using a Bayesian Markov chain Monte Carlo (MCMC) algorithm. For consistency MCMC was used to fit all the models, using WinBugs 1.4 software [14] and treating the proportion of children infected in each sampling unit (i.e. each 14 day period within each cohort) as binomial with probability E(h_{ ij }). Uniform(0,1) prior distributions were assigned for the proportions S, g, S_{∞}, log-normal priors with large variance to ${E}_{c}^{*}$, ${E}_{f}^{*}$ and β, and a gamma prior for γ. For the models with negative binomial variation in the inoculation rate we assign inverse gamma priors to κ. Point and interval estimates of all the quantities in equations (4) to (10) were obtained by sampling their posterior distributions. The deviance information criterion (DIC) [15] was used to compare the statistical fit of the models.
Results
Initial analyses considered whether these 83 sampling units can be considered as statistically independent. At the time of the study, parasitological resistance to sulphadoxine-pyrimethamine was extremely infrequent in Kenya, so it is assumed that this treatment cleared all pre-existing infections, and hence there should be no carry-over effects from previous infections. The children with parasites at the start of each period of follow-up were were not included in the risk set for subsequent time points, so no more than one infection event was recorded for any one individual, and the successive observations are statistically independent. However, it remained a priori possible that the relative risk of infection in a cohort might decrease over time, because the more susceptible children will be infected (and censored from the cohort) sooner. This possiblity was evaluated by using logistic regression to compare the risk of infection (h_{ i }) during the 3^{rd} and 4^{th} fortnights of follow-up, with that during the 1^{st} and 2^{nd} fortnights respectively in the parallel cohorts that started follow-up 4 weeks later. This analysis gave an odds ratio of 1.25 (95% confidence intervals 0.98, 1.59), suggesting that if anything, infection rates tend to be higher in the latter part of the follow-up, thus justifying treating the different cohorts and sampling periods as statistically independent in the subsequent analyses.
Summary of model fit
Model for imprecision in inoculation rate | Model for relationship of infection to inoculation rate | Deviance Information Criterion | |||
---|---|---|---|---|---|
Inoculation rate measurement | Infection of children | Total | |||
1 | none | Mass Action | - | 772.7 | - |
2 | none | Facilitation | - | 1840.3 | - |
3 | none | Saturation | - | 528.1 | - |
4 | none | Competition | - | 486.9 | - |
5 | Poisson | Mass Action | 270.2 | 392.1 | 662.3 |
6 | Poisson | Facilitation | 356.7 | 411.9 | 768.6 |
7 | Poisson | Saturation | 252.8 | 383.6 | 636.4 |
8 | Poisson | Competition | 245.8 | 392.5 | 638.3 |
9 | negative binomial | Mass Action | 267.4 | 373.5 | 640.9 |
10 | negative binomial | Facilitation | 277.5 | 376.3 | 653.8 |
11 | negative binomial | Saturation | 268.8 | 371.8 | 640.8 |
12 | negative binomial | Competition | 266.8 | 373.9 | 640.7 |
At high EIR e_{ i }is determined more precisely, so there is less effect of the shrinkage in the EIR estimates at high e_{t(i)}than at low e_{t(i)}. Corresponding to this, the fitted curves relating infection rates to these EIR values are much more similar in Figure 1b (models 5–8) than in Figure 1a. Mass action (model 5) now estimates that 23% of inoculations result in infections. The deviance information criteria (DIC) for these models comprise two terms, one measuring the fit to the EIR data, and one to the data for infection. All four DIC for the fit to the infection data are lower than the value for the best fitting model without the Poisson error term, with saturation (model 7) fitting best among those including Poisson variation in EIR (Table 1), estimating an upper limit of 65% of children becoming infected. The best fit to the EIR data is achieved with competition (model 8), but overall this model has a slightly inferior fit to that of saturation (model 7).
When negative binomial variation in the EIR is assumed, all of the models 9–12 estimate similar degrees of overdispersion, and the estimated underlying EIR values are even more tightly distributed around the average (compare Figure 1c with Figure 1b). A consequence is that the fitted curves for different variations from mass action are also closer to each other than the curves shown in either Figures 1a or 1b. In the range of 1–10 infectious bites per fortnight, into which most of the estimates of EIR fall, the curves for mass action, saturation and competition models are almost superimposed. Over this range of EIR only the facilitation model gives a markedly different estimate of the relationship between infection and EIR, predicting lower infection rates that those observed.
The DIC values for the different negative binomial models are also similar to each other, though facilitation (model 10) again fits markedly worse than mass action (model 9), saturation (model 11), or competition (model 12). Saturation (model 11) now has the best fit to the infection rates (as measured by the DIC), though the level at which saturation occurs is almost the same whether Poisson or negative binomial errors for the EIR are assumed. However, the overall best fit among negative binomial models is achieved assuming competition (model 12), with the mass action only a little inferior. When all twelve models are compared, saturation with Poisson errors (model 7) gives the best fit, as assessed by the DIC (which penalizes the negative binomial models for the extra variance parameter estimated). From a statistical point of view, the negative binomial model is more complicated and the trade off between improvement in fit and complexity of the model tends to speak in favour of a Poisson model.
Parameter values measuring overdispersion of the negative binomial models.
Model | κ | 95% interval | r | (95% interval) |
---|---|---|---|---|
Mass Action | 0.60 | (0.34, 0.93) | 1.78 | (1.08, 2.95) |
Facilitation | 0.84 | (0.56, 1.16) | 1.23 | (0.86, 1.77) |
Saturation | 0.62 | (0.31, 1.10) | 1.78 | (1.68, 3.22) |
Competition | 0.58 | (0.32, 0.90) | 1.85 | (1.11, 3.09) |
Estimated parameter values
Model for inoculum survival | Units | No imprecision in inoculation rate | Poisson imprecision in inoculation rate | Negative binomial imprecision in inoculation rate | ||||
---|---|---|---|---|---|---|---|---|
point estimate | 95% interval | point estimate | 95% interval | point estimate | 95% interval | |||
mass action | S | proportion | 0.18 | 0.15, 0.20 | 0.18 | 0.15, 0.20 | 0.23 | 0.17, 0.31 |
facilitation | ${E}_{f}^{*}$ | inoculations/fortnight | 13.4 | 11.8, 15.1 | 15.1 | 11.6, 19.2 | 27.2 | 10.7, 52.8 |
facilitation | γ | dimension-less | 1.17 | 1.00, 1.60 | 1.02 | 1.00, 1.08 | 1.00 | 1.00, 1.01 |
saturation | g | proportion | 0.43 | 0.39, 0.47 | 0.65 | 0.56, 0.77 | 0.76 | 0.61, 0.96 |
saturation | β | fortnights/inoculation | 1.89 | 1.40, 2.51 | 0.42 | 0.28, 0.62 | 0.29 | 0.17, 0.45 |
competition | ${E}_{c}^{*}$ | inoculations/fortnight | 0.40 | 0.30, 0.52 | 0.22 | 0.10, 0.39 | 0.01 | 0.00, 0.03 |
competition | S _{ inf } | proportion | 0.06 | 0.03, 0.08 | 0.11 | 0.07, 0.15 | 0.06 | 0.04, 0.09 |
Discussion
Micro-heterogeneity in exposure to infectious mosquito bites contributes to sampling variation in estimates of the average exposure to P. falciparum and hence inflates the variation between exposure estimates for different time periods. It also modifies observed relationships between infection rates and average exposure, because heterogeneous exposures reduce the proportion of individuals that get infected. This effect, the result of superinfections occuring disproportionately in some individuals, becomes greater as exposure increases.
The analyses of this paper compare different distributional assumptions for the sampling variation that contributed to the observed average time-period specific exposures and make shrinkage estimates of the true average exposures that allow for this sampling variation.
The same distributions of entomological inoculations around the population average are used to infer the dependence of the probability of infection on the average EIR. By assuming mass action kinetics and taking advantage of the availability of both infection and exposure data from the Saradidi site, the models provides estimates of the true degree of micro-heterogeneity that best accounts for observed infection-exposure relationship.
Mass action kinetics must apply at the level of individual host parasite interactions but field data on the kinetics of P. falciparum transmission from the mosquito vector to the human host seem at first sight to be incompatible with this. This is not a consequence of facilitation, rather the risk for a human of becoming infected for each sporozoite positive mosquito bite appears to decrease steeply with increased transmission intensity. This result is confirmed by both studies of highly exposed cohorts[7, 16], and comparisons of the force of infection between laboratory challenge, low transmission areas, and high transmission areas.
Although epidemiological evidence suggests that acquired immunity offers little or no protection against pre-erythrocytic stages of the parasite in children[11] apparent competition between inocula could result from immunological cross reactions. In particular, highly exposed individuals are likely to have been more highly exposed in the past and so may have reduced levels of detection of parasites, which could lead to an overall bias in estimates of S. But in the present dataset, the exposure variation arises as a result of short-term temporal variation, and the cohorts are sampled from a common age group in the same population, so acquired immunity cannot account for differences between sampling units.
Innate immunity does seem to play a role in limiting super-infections of P. falciparum at the hepatic stage. In particular, asexual parasites stimulate production of interferon-γ, which leads to killing of infected hepatocytes[17, 18]. Some malaria models, such as those of [12] and [11] have used deterministic functions to capture such apparent density dependent regulation of the inoculation rate. However only if the host is already infected does it mount any of these defensive responses. Feedback from pre-existing infections cannot explain the apparent density dependence in the Saradidi dataset, where all the asexual parasites were cleared at the start of the follow-up.
The known biology of malaria thus gives us no reason to suppose that there is any dependence in the infection probability for successive parasites challenging a given uninfected host, and hence all attempts to model deviations from mass action by treating the success probability of the parasite as an explicit function of the intensity of the parasitological challenge introduce a divorce between the mathematical model and the biology of the infection process. Heterogeneities in transmission therefore provide a more parsimonious explanation for apparent competition than do density-dependent host responses. In sub-populations with higher contact rates the probability of becoming infected is higher, but the increase is not proportional to the contact rate because there is a higher chance of an infection being 'wasted' on a host that has already been colonized, so that heterogeneity mimics the effect of density dependence on infection risk.
The present analyses show that models of heterogeneity can fit field data as well as models of density dependence. Humans are known to be heterogeneous in their attractiveness to mosquitoes[19, 20] even in standardized environments. In the field setting, heterogeneities are to be expected both in the exposure and susceptibility of the children within any one cohort and in the sampling efficiency of the different mosquito collectors.
Quantifying focality in transmission via the negative binomial provides a parsimonious model for this micro-heterogeneity that readily accounts for the deviations from linearity in Figure 1A, and gives an estimate of κ that is otherwise unobtainable. Although the dataset is the best available with which to estimate of κ by this approach, the study was not designed with this purpose in mind and some approximations were required. The most straightforward field design to estimate the level of micro-heterogeneity would have involved estimating the EIR for each fortnight by sampling a single random individual over all 14 nights, corresponding directly to a sample from z_{t(i)}. The data were treated as though they had been generated this way by applying the same estimate of κ = 0.6 (for the mass action model in Table 2) both to the distribution of fortnightly exposure in children, translated into the formula for infection risk via equation (15), and to the distribution of fortnightly EIR estimates compiled by summing 12 distinct nightly human landing collections (and rescaling them). In the actual design, the employment of 6 different mosquito collectors rather than one only, presumably brought the overall EIR estimates closer to the population average. Conversely the shorter sampling periods, involving sampling only 12 of the 14 nights during each fortnight, leads to inflation of the variance of the EIR estimate relative to that in the ideal design. These two effects appear to have approximately cancelled each other out: when the negative binomial mass action model was extended to allow distinct values of κ for the sampling of the EIR and for the formula for infection risk, both κ estimates were close to the original value of 0.6, with the difference between them only 0.06 (95% CL for the difference: -0.57, 0.65).
The Saradidi environment can be assumed to be reasonably typical of rural East Africa, and the estimation of κ represents a critical step in parameterising mathematical models for heterogeneity in malaria transmission in such environments. The next step, of partitioning the variance into variation in space, variation between different human hosts, and between different nights during the fortnight, requires data from other sites.
A complete model for heterogeneity in malaria transmission also needs to consider transmission from human to vector. This also seems to defy mass action. Detailed records of both gametocyte densities and transmissibility to vectors are available from neurosyphilis patients who received therapeutic malaria infections prior to the availability of antibiotics. Higher asexual parasite densities are associated with higher gametocyte densities, and these in turn are associated with greater transmission to the vector [21], but the force of infection of the vector does not increase in proportion to host parasite densities, even allowing for the need of oogenesis for gametocytes of both sexes. A stochastic model was recently fitted to these data that accounts for the observed lack of proportionality as a consequence of random variation in gametocytogenesis[22]. In conjunction with the estimate of κ = 0.6 for vector-human transmission, this makes it possible to use field based estimates of heterogeneity in models of the whole P. falciparum transmission cycle.
Conclusion
The present study presents a method for estimating the extent of micro-heterogeneity in exposure to infected malaria vectors, at the same time demonstrating that heterogeneity can explain apparent deviations from mass action kinetics. The estimates obtained can be used in formulating dynamic models of the whole malaria transmission cycle, in which stochastic variations in transmission rates are used to mimic the actual levels of heterogeneity found in the field. This will make it possible to model the impact of malaria control interventions while allowing for realistic levels of heterogeneity in transmission.
Declarations
Acknowledgements
The author would like to acknowledge helpful comments from John Beier, Dieter Ebert and Amanda Ross and from two anonymous reviewers. This work forms part of Swiss National Science Foundation project 3300C0-105994/1.
Authors’ Affiliations
References
- Woolhouse ME, Dye C, Etard JF, Smith T, Charlwood JD, Garnett GP, Hagan P, Hii J, Ndhlovu PD, Quinnell RJ, Watts CH, Chandiwana SK, Anderson RM: Heterogeneities in the transmission of infectious agents: implications for the design of control programs. Proc Natl Acad Sci USA. 1997, 94: 338-342. 10.1073/pnas.94.1.338.PubMed CentralView ArticlePubMedGoogle Scholar
- Macdonald G: The epidemiology and control of malaria. 1957, London, Oxford University PressGoogle Scholar
- Smith DL, Dushoff J, Snow RW, Hay SI: The entomological inoculation rate and Plasmodium falciparum infection in African children. Nature. 2005, 438: 492-495. 10.1038/nature04024.PubMed CentralView ArticlePubMedGoogle Scholar
- Filion GJ, Paul RE, Robert V: Transmission and immunity: the importance of heterogeneity in the fight against malaria. Trends Parasitol. 2006, 22: 345-348. 10.1016/j.pt.2006.06.005.View ArticlePubMedGoogle Scholar
- Wilson EB, Worcester J: The Law of Mass Action in Epidemiology. Proc Natl Acad Sci USA. 1945, 31: 24-34. 10.1073/pnas.31.1.24.PubMed CentralView ArticlePubMedGoogle Scholar
- Regoes RR, Ebert D, Bonhoeffer S: Dose-dependent infection rates of parasites produce the Allee effect in epidemiology. Proc Biol Sci. 2002, 269: 271-279. 10.1098/rspb.2001.1816.PubMed CentralView ArticlePubMedGoogle Scholar
- Beier JC, Oster CN, Onyango FK, Bales JD, Sherwood JA, Perkins PV, Chumo DK, Koech DV, Whitmire RE, Roberts CR: Plasmodium falciparum incidence relative to entomologic inoculation rates at a site proposed for testing malaria vaccines in western Kenya. Am J Trop Med Hyg. 1994, 50: 529-536.PubMedGoogle Scholar
- Beier JC, Perkins PV, Onyango FK, Gargan TP, Oster CN, Whitmire RE, Koech DK, Roberts CR: Characterization of malaria transmission by Anopheles (Diptera: Culicidae) in western Kenya in preparation for malaria vaccine trials. J Med Entomol. 1990, 27: 570-577.View ArticleGoogle Scholar
- Port GR, Boreham PFL, Bryan JH: The relationship of host size to feeding by mosquitos of the Anopheles-gambiae giles complex (Diptera, Culicidae). Bull Entomol Res. 1980, 70: 133-144.View ArticleGoogle Scholar
- Carnevale P, Frezil JL, Bosseno MF, Lepont F, Lancien J: Study of Agressivity of Anopheles-Gambiae A in Relation to Age and Sex of Human Subjects. Bull WHO. 1978, 56: 147-154.PubMed CentralPubMedGoogle Scholar
- Smith T, Maire N, Dietz K, Killeen GF, Vounatsou P, Molineaux L, Tanner M: Relationship between the entomologic inoculation rate and the force of infection for Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006, 75 (2 Suppl): 11-18.PubMedGoogle Scholar
- Dietz K, Molineaux L, Thomas A: A malaria model tested in the African savannah. Bull World Health Org. 1974, 50: 347-357.PubMed CentralPubMedGoogle Scholar
- Smith T, Killeen G, Lengeler C, Tanner M: Relationships between the outcome of Plasmodium falciparum infection and the intensity of transmission in Africa. Am J Trop Med Hyg. 2004, 71 (2 Suppl): 80-86.PubMedGoogle Scholar
- Spiegelhalter DJ, Thomas A, Best N, Lunn D: Winbugs Version 1.4. 2003, Cambridge, England, MRC-BSUGoogle Scholar
- Spiegelhalter DJ, Best N, Carlin BP, van der Linde A: Bayesian measures of model complexity and fit (with discussion). J Roy Statist Soc B. 2002, 64: 583-640. 10.1111/1467-9868.00353.View ArticleGoogle Scholar
- Charlwood JD, Smith T, Lyimo E, Kitua A, Masanja H, Booth M, Alonso PL, Tanner M: Incidence of Plasmodium falciparum infection in infants in relation to exposure to sporozoite-infected anophelines. Am J Trop Med Hyg. 1998, 59: 243-251.PubMedGoogle Scholar
- Schofield L, Ferreira A, Altszuler R, Nussenzweig V, Nussenzweig RS: Interferon-gamma inhibits the intrahepatocytic development of malaria parasites in vitro. J Immunol. 1987, 139: 2020-2025.PubMedGoogle Scholar
- Gonzalez-Aseguinolaza G, de Oliveira C, Tomaska M, Hong S, Bruna-Romero O, Nakayama T, Taniguchi M, Bendelac A, Van Kaer L, Koezuka Y, Tsuji M: alpha -galactosylceramide-activated Valpha 14 natural killer T cells mediate protection against murine malaria. Proc Natl Acad Sci U S A. 2000, 97: 8461-8466. 10.1073/pnas.97.15.8461.PubMed CentralView ArticlePubMedGoogle Scholar
- Knols BG, De Jong R, Takken W: Differential attractiveness of isolated humans to mosquitoes in Tanzania. Trans R Soc Trop Med Hyg. 1995, 89: 604-606. 10.1016/0035-9203(95)90406-9.View ArticlePubMedGoogle Scholar
- Mukabana WR, Takken W, Killeen GF, Knols BG: Allomonal effect of breath contributes to differential attractiveness of humans to the African malaria vector Anopheles gambiae. Malar J. 2004, 3: 1-10.1186/1475-2875-3-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Jeffery GM, Eyles DE: Infectivity to mosquitoes of Plasmodium falciparum as related to gametocyte density and duration of infection. Am J Trop Med Hyg. 1955, 781-789.Google Scholar
- Ross A, Killeen G, Smith T: Relationships between host infectivity to mosquitoes and asexual parasite density in Plasmodium falciparum. Am J Trop Med Hyg. 2006, 75 (2 Suppl): 32-37.PubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.