INTRODUCTION
The prevention of new infections by stimulating immunity against sporozoites or liver stages of Plasmodium falciparum parasites is a major strategy being pursued for the development of malaria vaccines.^{1} There is substantial immunologic evidence that immune responses to infected hepatocytes occur in nature, and the effectiveness of inoculation with irradiated sporozoites demonstrates the feasibility of stimulating effective protective responses,^{2} as do the effects demonstrated in recent challenge and field trials of the RTS,S/AS02A vaccine against pre-erythrocytic stages of the parasite.^{3}^{–}^{5}
However, there is little evidence from field studies that naturally acquired pre-erythrocytic immunity is of epidemiologic importance. The most obvious effect of acquired immunity in malaria is to reduce asexual parasite densities without preventing infection, and people exposed repeatedly to P. falciparum infections over a lifetime do not become completely refractory to infection. The incidence of new infections in adults in settings of high transmission is surprisingly similar to that in their children. In Dielmo, Senegal, reinfection rates were lower in adults than children;^{6} in Saradidi in western Kenya, daily attack rates in adults^{7} were approximately half those observed in children;^{8} while in Navrongo in northern Ghana, incidence of infection in a cohort of adults^{9} was similar to that of young children.^{10} The proportion of inoculations that result in infection decreases as the intensity of transmission increases.^{8}^{,}^{11}
The most accepted integrated models of malaria transmission dynamics and immunity to date remain those of the Garki project^{12} and variants of it.^{13}^{,}^{14} These are complex mass action models based on calculation of transition rates between compartments within the host population with allowance for latent periods. These models assume that the relationship between the entomologic and epidemiologic inoculation rates is independent of prior exposure. This is equivalent to assuming that there is effectively no immunologic memory of natural immunity against pre-erythrocytic stages of the parasite.
To estimate the effects of such immunity, and how rapidly this is acquired, data are needed from epidemiologic settings where incidence of infection has been compared across a range of age-groups or exposure histories. This requirement is fulfilled by data from the village of Matsari in northern Nigeria, which was studied as part of the Garki project,^{15} where mass treatment was administered every 10 weeks for a total of 80 weeks during 1972–1973, and reinfection was recorded in the whole village population before each round of treatment.
Anopheline mosquitoes bite larger people more often than smaller ones.^{16}^{,}^{17} When this is allowed for, it is apparent that the success rate of inoculations in young children in endemic areas must be much greater than that in adults in the same communities. Otherwise, the high incidence of clinical infections in the youngest children cannot be explained,^{18} and incidence rates in adults in studies such as those in Saradidi and Navrongo would be much greater than those in children. Therefore, naturally acquired immunity that prevents blood stage infection should not be ignored in malaria models, although epidemiologic data alone cannot strictly distinguish whether this is directed against pre-erythrocytic stages or whether it represents an infection-blocking response to some asexual blood stages. These effects of host size in malaria transmission have been considered in mathematical models of malaria epidemiology.^{12}^{–}^{14}^{,}^{19}
With the objective of making realistic predictions of the potential impact of malaria vaccines and other malaria control interventions, we are developing a new dynamic model for the epidemiology of P. falciparum.^{20} This model considers (among other phenomena) the effects of host size on the rate at which the host is bitten by anophelines, the relationship between the entomologic inoculation rate (EIR) and the infection rate, and decreasing survival of the inoculations as the hosts acquire pre-erythrocytic immunity. We now report the statistical fitting of this model to the re-infection data from both Saradidi and Matsari and use this to make inferences about the importance of immunity that prevents the establishment of infections.
MATERIALS AND METHODS
Description of the model
Differential feeding by mosquitoes depending on body surface area
To estimate the age-specific EIR assuming proportionality between mosquito bites received and the body surface area of the host, we estimated the mean body surface area by age for a typical African rural population (from Tanzania)^{21} using the formula of Mosteller.^{22} Defining A(a(i,t)) as the mean body surface area for individuals in the same age group as individual i, at time t, age group a, we obtained E_{a}(i,t), the age-adjusted EIR, from
where, A(a(i,t)) is the average body surface area estimated for an individual of age (a(i,t) and A_{max} is the average surface area of people ≥ 20 years of age in the same population. E_{max}(t) refers to the usual measure of the EIR computed from human bait collections made by fully grown mosquito collectors.
Control of pre-erythrocytic stages
We use equation 1 to estimate for each time and for each individual i, the number of infective bites received per unit time, adjusted for age, E_{a}(i,t). We then define a survival function S_{p}(i,t) denoting the probability that the progeny of each inoculation survive to give rise to a patent blood stage infection, and model the force of infection as
We then assume that the number of infections introduced in unit time in each individual follows a Poisson distribution with mean λ(i,t), i.e.,
The probability S_{p}(i,t) is formed as the product of two terms,
Respectively:
S_{l}(i,t) captures innate density dependent effects that ensure that as E_{a}(i,t) increases, the proportion of inoculations that result in infections decreases. We assume that at any time all individuals are susceptible. At very low E_{a}(i,t) in the naive host this corresponds to the proportion of infectious bites that lead to an infection, i.e., λ(i,t) = E_{a}(i,t) and S_{l}(i,t) = 1. As E_{a}(i,t) increases, the proportion of inoculations that are successful decreases. A functional form that satisfies these requirements for the relationship between entomologic and epidemiologic inoculations in naive individuals is
where S_{∞}represents the lower limit attained byS_{l}(i,t) as theinoculation rate is large, and E* is the value of E_{a}(i,t) at which half the reduction in S_{l}(i,t) is achieved. These parameters are estimated with the constraints 0 < S_{∞}< 1 and E* > 0.
S_{2}(i,t) measures the effects of acquired pre-erythrocytic immunity, which leads to reduced take in hosts who have been previously exposed and who have accrued some level of exposure to pre-erythrocytic challenge, X_{p}(i,t), where
and a(i,t) is the age at time t.
Exploratory analyses of the data indicate that a good fit is obtained by assuming that at high X_{p}(i,t) there is a lower limit to the susceptibility, and that the relationship of susceptibility to X_{p}(i,t) can be captured using a Hill function as follows:
where S_{imm }is the maximal survival of the inoculum in the most immune individuals, and X_{p}* is a critical value of the cumulative exposure. Assuming independence of the effects of transmission intensity and acquired immunity, we combine equations 1–7 to obtain
Immunologic evidence suggests that the liver stages of the parasite are also controlled by innate mechanisms related to the density of erythrocytic stages. In particular, infected hepatocytes are highly sensitive to interferon-γ, the production of which is stimulated by asexual forms of P. falciparum. We propose that a further control is applied to liver stage infections that have survived the effect of acquired pre-erythrocytic immunity five days after inoculation. At this point, only a random sample, S_{h}(i,t), of the infections proceed where
and Y_{h}* is a parameter determining the strength of this innate immune effect and Y(i,t) is the current density of asexual blood stage parasites.
Data sources
Biting rate in relation to human weight, The Gambia
Port and others^{16} reported a series of experiments in The Gambia in which fed mosquitoes were collected from the mosquito nets of 35 groups of people who normally slept together. The blood meals were assigned to hosts using haptoglobin or ABO typing, and the proportion of mosquitoes that had fed on each host was analyzed in relation to the host’s contribution to the total biomass and surface area of the people sleeping in the mosquito net. Similar relationships were found between the biting rate and both weight and the surface area of the host (Figure 1), but the investigators could not determine which of these was the key determinant of how many mosquitoes bite an exposed host.
Infection rates by EIR, Saradidi, western Kenya
Beier and others^{8} studied the reinfection of 21 cohorts of between 43 and 50 children between 6 months and 6 years of age, each followed-up during four successive two-week periods after initial clearance of P. falciparum parasites (Figure 2). They reported entomologically assessed exposures in terms of the numbers of sporozoite-positive mosquitoes biting during each two-week period.
Infection rates by age and EIR, Matsari village, Nigeria
Matsari village was monitored entomologically for four years (November 1970–November 1973) during the Garki project.^{23} The first two years corresponded to the pre-intervention phase, during which eight cross-sectional malariologic surveys of the whole village population and intensive entomological surveillance were carried out. The latter provided estimates of the EIR computed from human bait collections of mosquitoes and dissections of the mosquito salivary glands for sporozoites (Figure 3) (for details of the methods see Molineaux and Gramiccia^{23}).
The intervention phase began in mid 1972, and an additional eight surveys were carried out at 10-week intervals (surveys 9–16). During the intervention phase, indoor residual spraying with Propoxur was carried out comprehensively in the village, as was mass treatment of the population, again at 10-week intervals, immediately after assessment of their parasitologic status. We analyze the data from cross-sectional surveys 9–16, which provide an assessment of infection during the preceding 10-week period.
Model implementation and fitting
Estimation of the model parameters required fitting of γ_{p}, S_{∞} S_{imm}, E*, X_{p}* to both the data from Matsari and those from Saradidi. Because the incidence estimates in these studies were made by following-up hosts in whom the parasites are cleared, they do not provide any data with which to estimate the parameter Y_{h}*. Thus, we report here models that ignored this control on liver stages of the parasite.
Exploratory analysis of the Matsari data indicated that the effect of acquired immunity, 1 −S_{2}(i,t), must be small in children even in highly endemic areas. To estimate the parameters E* and S_{∞}, we therefore fitted the model of equations 1–8 to the data from Saradidi conditional on a value of 1 for S_{2}(i,t). To obtain the correct effect of host size, we assumed the children in Saradidi to have a uniform age distribution, and the proportion that became infected in each two-week period to be binomially distributed about the proportion predicted by the model. We fitted the model using a Bayesian Markov Chain Monte Carlo algorithm in the software Winbugs.^{24}
To estimate and S_{imm}, X_{p}* γ_{p}, we fitted the full model of equations 1–8 to the data from Matsari conditioning on the estimates of E* and S_{∞} obtained from Saradidi. We simulated infections for a human population exposed since birth to the pre-intervention seasonal pattern of EIR in Matsari using a proposal for the parameter vector and equations 1–8 to determine the infection probability and immunologic status of each individual at each time point. We then continued the simulation to include the intervention period, assuming each mass-treatment to clear all the parasites, and using the measured EIR for Matsari during the intervention phase. In a first model (A) we predicted the age-prevalence pattern at each survey during the intervention period by making deterministic estimates for each individual of λ(i,t) and X_{p}(i,t), and thus of the probability that they became infected during the preceding 10-week (70-day) interval, i.e.,
Taking the mean of this for all individuals in each age group as our prediction of the age-specific prevalence, we computed the likelihood for the proposal assuming binomial errors.
In a second model (B), we used a model for parasite densities to allow for sub-patent infection in semi-immune individuals. The model for the infection process was the same, but the parasite densities were generated via a stochastic simulation model fitted jointly to the Matsari data and to data from other epidemiologic settings.^{25} This simulation provided predictions of the parasite density at each time point and each simulated individual that were compared with the limit of detection of the slide-reading actually used in Matsari to make predictions of the age-specific prevalence of patent parasitemia at each survey.
We used a simulated annealing algorithm^{26}^{,}^{27} to maximize the likelihood for Matsari conditional on the estimates of E* and S_{∞} obtained for Saradidi, and thus to obtain estimates of S_{imm}, X_{p}* and γ_{p}.
RESULTS
We explored the fit of a wide range of models to the Matsari and Saradidi datasets, but report detailed results only for those that conformed to what we know of the biology of malaria and simultaneously gave a good fit to the data. As well as considering the fit to the Matsari and Saradidi data, we also examined whether the models were consistent with age-prevalence patterns in untreated individuals from a range of epidemiologic settings. This consideration led us to prefer models assuming the bites received by the host to be proportional to A(a(i,t)), rather than those that we presented previously^{18} that relied on host weight. Whereas surface area increases with the square of the linear dimensions of the host, weight increases with the cube of the linear dimensions, and hence shows a steeper relationship with age (Figure 4a). When body weight was used we could not easily reconcile the very steep age-inoculation relationship with observed age patterns of parasitemia in untreated individuals even by introducing a model for acquired immunity.
One aim was to obtain a parsimonious model, and we were able to obtain a good fit using only the five parameters γ_{p}, S_{∞,}S_{imm}, E*, and X_{p}*. To model the relationship between EIR and the survival of the inoculum we required a sigmoidal function with the properties that at EIR = 0, the force of infection λ(i,t) = 0, and that as the EIR increases, λ(i,t) must increase, although S_{1}(i,t) should decrease monotonically. Although there is considerable scatter in the data from Saradidi, we were able to obtain a good fit that satisfies these constraints (Figures 2 and 4b). There is no evidence from the data that the proportion of hosts becoming infected saturates at high EIR, although this is a feature both of those malaria models based on that of Garki,^{12} and of the statistical model that we fitted previously to the Saradidi data.^{18} We could not improve the fit of our model by allowing for inoculations that do not result in blood stage infection at the lowest inoculation rates, i.e., our best fitting model constrains S_{1}(i,t) = 1 when E(i,t) = 0.
The EIR in Matsari for the 1971 transmission season was approximately 67 inoculations per person per year (Figure 3).^{23} We estimated that vector control reduced this during the 1972–1973 seasons to an average of 5.5 infectious bites per year. This estimate is approximate because only 5 of 613 An. gambiae s.s. analyzed from this period were sporozoite positive, and our estimate of the variations in EIR therefore depends mainly on the variation in the measured human biting rates.
Even at this relatively low inoculation rate, many individuals became infected, and 226 (8%) of 2,663 follow-up blood slides were positive for P. falciparum. The average proportion of slides positive was maximal in the 10–15-year-old age group. Younger children showed a gradual increase in slide positivity with age, while the infection rate in adults averaged only approximately 5% (Figure 5).
The Saradidi database contains much more information about incidence in the 0.5–6-year-old age group than does the Matsari database because of the large number of observations of non-zero prevalence that could be compared with the model predictions. The parameters that gave the best fit to the Saradidi database (Figure 2) were thus very similar to those of the jointly best-fitting model (Table 1). The model giving the best fit to Matsari conditional on the values of S_{∞}, E* fitted to Saradidi (model A) gave an excellent fit to the data for adolescents and adults (Figure 5), with the shift from the non-immune state to the state of pre-erythrocytic immunity appearing to occur extremely abruptly when the host has been exposed to approximately 420 infectious bites (Figure 4c). Since the transition appeared to be so abrupt, this model estimated that young children had no acquired immunity and consequently could not give a good fit to the observed incidence rates in children in Matsari (Figure 5), for which it predicted higher values than those observed.
A major simplification made when fitting the model initially (model A) was to assume that all new infections are patent at the time of sampling. This is not necessarily the case, and biases the estimates especially for adults (in whom many infections are sub-patent) and in children less than six months of age who still have some maternally derived protection leading to very low parasite densities. For this reason we also present the results (model B) of fitting the model for infection jointly with a model for the control of parasite densities to the Matsari data (allowing for sub-patent parasitemia in Matsari) as well as to an additional six datasets comprising repeated cross-sectional survey data. Parameter estimates from model B are given in Table 1, the fitted curves are again shown in Figure 5, and the functions estimated are shown in Figure 4c.
As a validation exercise we compared the relationship between the proportions of infants converting and the EIR (E_{max}(t)) predicted by model B with the data from the study of Charlwood and others^{11} for the village of Idete in south-central Tanzania (Figure 6). We used our model (including functions representing blood stage immunity^{25}) to simulate this study including the survey frequency and parasitologic procedures. There was a remarkably good congruence between the two sets of proportions, although the field data were lower than the predictions made using model B at very high transmission levels. At these high levels of transmission saturation was observed in Idete, but it was not predicted by our model.
DISCUSSION
The challenge of fitting a model to predict the force of infection for malaria from the EIR has led us to several insights about the factors that limit the frequency of infection, including host size and both acquired and innate immunity. Our strategy allowed us to consider these factors separately both from those determining the levels of immunity that control asexual blood stage densities^{25} and from those modifying the clinical response to specific parasite loads.^{28}
Lower exposure of small children to bites from anophelines has substantial epidemiologic implications that remain to be fully explored. In areas of high transmission, prevalence, multiplicity, and in many situations the incidence of clinical malaria, all show non-monotonic relationships with age,^{29}^{,}^{30} implying that some risk factor increases with age at the same time as immunity decreases the risk. We propose that this risk factor is the increase in exposure to anopheline biting as children grow older. The higher exposure of adults can also perhaps explain why they appear to acquire immunity more quickly than children when exposed to endemic malaria transmission for the first time.^{31}
An upper limit to the rate of increase in the infection rate with age is imposed by the assumption that this increase results entirely from the growth in body size. Given this constraint, we found that only models that assumed a rather abrupt onset of control of infection can provide a good fit to the much lower force of infection in adults than in older children in Matsari. With only the Matsari data we are not able to identify whether this control is caused by acquired immunity or to an innate age-dependent mechanism. However, the high susceptibility to P. falciparum infection of immunologically naive malariatherapy patients^{32} and volunteers in challenge trials^{33} suggests that acquired immunity is the more likely explanation.
When we treated all simulated infections as patent, the best fitting model for Matsari (model A) estimated a very abrupt onset of pre-erythrocytic immunity, with X_{p}* = 523 inoculations corresponding to a good fit for the data for adults (Figure 5). The model (model B), which was fitted jointly with a model for asexual blood stage immunity^{25} and allowed for sub-patent infections, explained much of the decrease in prevalence with age as a result of acquired immunity to asexual blood stages. Model B predicted sub-patency in a high proportion of the infections in adults and more gradual acquisition of pre-erythrocytic immunity. Although its fit to age-prevalence curves in adults in Matsari was inferior to that of model A, it gave better predictions of the force of infection in the Matsari children and of age-prevalence patterns in untreated populations in the Garki project.^{23}^{,}^{25} This model thus appears more plausible than model A and we have adopted it in our further simulations.
An important feature of all our models is the reduction of the survival probability of the inoculum with an increase in the inoculation rate. It has long been known that conversion rates in infants are much less than the EIR,^{34}^{,}^{35} and it was apparent during the development of the Garki model^{12}^{,}^{23} that some pre-erythrocytic density-dependent regulation of the inoculation rate is needed for realistic simulations of patterns of transmission. We have constrained the success probability of the inoculations to decrease as E_{a}(i,t) increases (Figure 4b), but based on the pattern evident in Saradidi, and in distinction from the Garki and related models,^{14}^{,}^{36} saturation is never reached. One likely factor contributing to this pattern is host heterogeneity in susceptibility. The only host heterogeneity that our model explicitly includes is the effect of host age, but humans are known to be heterogeneous in their attractiveness to mosquitoes.^{37}^{,}^{38} It is likely that an improved model could be formulated by incorporating effects of additional sources of variation, such as variation between hosts in innate susceptibility to infection, and explicit modeling of heterogeneity in exposure to mosquitoes. The latter would lead to extra-Poisson variation in the infection process. However, in the absence of appropriate field data on which to base such extensions, we consider that model B represents an appropriate representation of the various factors determining the force of infection for P. falciparum for use in predictive models.
An additional factor that contributes to the apparent reduction in the survival of the inoculum with increased EIR is likely to be blocking of super-infection by innate immunity to hepatic stages. This could be stimulated by either hepatic or erythrocytic stages of P. falciparum. Infected hepatocytes are highly sensitive to interferon-γ, the production of which is stimulated by asexual forms of P. falciparum^{39} and which in turn stimulates hepatic anti-malaria activity by natural killer T cells.^{40} Although we were able to propose a model for feedback from asexual parasitemia on pre-erythrocytic stages, all asexual parasites were cleared at the start of the studies we used to assess infection rates; thus, we could not fit this model to these data. We attempted to estimate this feedback as part of our model for blood stage immunity, but incorporation of this effect did not improve the fit. However, equations 5 and 9 are competing for simulating innate density-dependent regulation, and available data may be inadequate to dissociate their respective effects.
At very low transmission levels distinct inoculations of the same host are widely separated in time. Thus, it is unlikely that there is any real competition between them. However, even under controlled conditions for therapeutic or experimental trials with very high sporozoite loads, there is less than 100% take in inoculations of P. falciparum.^{32}^{,}^{41}^{,}^{42} Our models predict that at the very lowest exposure in naive individuals all inoculations will be successful (Figure 4b). However, this is an extrapolation outside the range of EIRs that can be estimated in field studies, and the success probability at the lower limit of measurable EIRs (of the order 0.003 inoculations per night) should be taken as our best estimate of the upper limit of the survival of inocula.
The validation of model B using the data of Charlwood and others^{11} confirms that it can give estimates of conversion proportions similar to those from field studies other than those to which the model was fitted. At very high transmission intensities these data, in contrast to those of Beier and others,^{8} suggest that saturation occurs. Infections are inapparent in some of the youngest children studied in Idete, who included newborns protected by maternal immunity.^{43} There is a need for further studies to determine whether such saturation occurs in other settings.
We conclude that the epidemiologic evidence suggests that there is naturally acquired immunity that prevents establishment of some blood stage infections of P. falciparum. However, this is acquired even more slowly than immunity that controls asexual blood stage densities and plays an important role only in individuals who have already been heavily exposed to the parasite. Pre-erythrocytic vaccines that aim to prevent infections in the most vulnerable group (very young children) will need to stimulate responses that are stronger or qualitatively different from those observed in nature.
Estimated parameter values
Parameters determined from Saradidi data | Units/dimensions | Point estimate (95% confidence limits) | |
---|---|---|---|
S _{∞} | Lower limit of success probability of inoculations at high E_{a}(i,t) | Proportion | 0.049 (0.026, 0.073) |
E* | Critical value of E_{a}(i,t) | Inoculations/person-night | 0.032 (0.025, 0.041) |
Log likelihood (Saradidi) | −1,159.3 |
Parameters determined from Matsari data | Model A (assuming no sub-patent infections) | Model B (allowing for sub-patent infections in Matsari data)* | ||
---|---|---|---|---|
* Confidence intervals for Model B estimated as described by Maire and others.^{25} | ||||
S _{ imm } | Lower limit of success probability of inoculations in immune individuals | Proportion | 0.12 | 0.14 (0.13, 0.15) |
γ_{p} | Steepness of relationship between S_{3}(i,t) and X_{p}(i,t) | Dimensionless constant | 5.1 | 2.04 (1.45, 2.89) |
X_{p}* | Critical value of cumulative number of entomologic inoculations | Inoculations | 523.0 | 1,514.4 (1,441.3, 1,591.3) |
Log likelihood (Matsari) | −778.2 | −779.7 |
Incidence by entomologic inoculation rate in Saradidi, Kenya. The filled circles indicate observed proportions of children 0.5–6 years of age who became infected in a two-week period. The line gives the fit of the model.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Incidence by entomologic inoculation rate in Saradidi, Kenya. The filled circles indicate observed proportions of children 0.5–6 years of age who became infected in a two-week period. The line gives the fit of the model.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Incidence by entomologic inoculation rate in Saradidi, Kenya. The filled circles indicate observed proportions of children 0.5–6 years of age who became infected in a two-week period. The line gives the fit of the model.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Inoculation rates and intervention history in Matsari, Nigeria. The vertical arrows indicate the time points of cross-sectional surveys of the village population. The letter M indicates that mass treatment of the population was carried out at the same time as the survey.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Inoculation rates and intervention history in Matsari, Nigeria. The vertical arrows indicate the time points of cross-sectional surveys of the village population. The letter M indicates that mass treatment of the population was carried out at the same time as the survey.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Inoculation rates and intervention history in Matsari, Nigeria. The vertical arrows indicate the time points of cross-sectional surveys of the village population. The letter M indicates that mass treatment of the population was carried out at the same time as the survey.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Proposed functions relating infection incidence to exposure. a, Effect of age on number of mosquito bites received. Shown is the ratio of bites received to those received by an adult (A(a(i,t))/A_{max}(t)) by age (a(i,t)). The continuous line corresponds to proportionality between bites received and expected body surface area (the model that we implemented); the dotted line assumes proportionality between bites received and body weight. b, Effect of entomologic inoculation rate (EIR) on survival of the inoculum. Shown is the proportion of inocula surviving in malaria-naive individuals, ( S_{1}(i,t)), by age-adjusted EIR (E_{a}(i,t)). c, Effect of cumulative exposure on survival of the inoculum. Shown is the proportion of inocula surviving (in the limit when (S_{1}(i,t)→1), (S_{2}(i,t)), by cumulative EIR, (X_{p}(i,t). Dashed line = model A; continuous line = model B.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Proposed functions relating infection incidence to exposure. a, Effect of age on number of mosquito bites received. Shown is the ratio of bites received to those received by an adult (A(a(i,t))/A_{max}(t)) by age (a(i,t)). The continuous line corresponds to proportionality between bites received and expected body surface area (the model that we implemented); the dotted line assumes proportionality between bites received and body weight. b, Effect of entomologic inoculation rate (EIR) on survival of the inoculum. Shown is the proportion of inocula surviving in malaria-naive individuals, ( S_{1}(i,t)), by age-adjusted EIR (E_{a}(i,t)). c, Effect of cumulative exposure on survival of the inoculum. Shown is the proportion of inocula surviving (in the limit when (S_{1}(i,t)→1), (S_{2}(i,t)), by cumulative EIR, (X_{p}(i,t). Dashed line = model A; continuous line = model B.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Proposed functions relating infection incidence to exposure. a, Effect of age on number of mosquito bites received. Shown is the ratio of bites received to those received by an adult (A(a(i,t))/A_{max}(t)) by age (a(i,t)). The continuous line corresponds to proportionality between bites received and expected body surface area (the model that we implemented); the dotted line assumes proportionality between bites received and body weight. b, Effect of entomologic inoculation rate (EIR) on survival of the inoculum. Shown is the proportion of inocula surviving in malaria-naive individuals, ( S_{1}(i,t)), by age-adjusted EIR (E_{a}(i,t)). c, Effect of cumulative exposure on survival of the inoculum. Shown is the proportion of inocula surviving (in the limit when (S_{1}(i,t)→1), (S_{2}(i,t)), by cumulative EIR, (X_{p}(i,t). Dashed line = model A; continuous line = model B.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Age-prevalence pattern in the Matsari, Nigeria intervention phase. ▪ = mean of observed prevalence; *—* = mean of predictions from model A; ○—○ mean of predictions from model B.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Age-prevalence pattern in the Matsari, Nigeria intervention phase. ▪ = mean of observed prevalence; *—* = mean of predictions from model A; ○—○ mean of predictions from model B.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Age-prevalence pattern in the Matsari, Nigeria intervention phase. ▪ = mean of observed prevalence; *—* = mean of predictions from model A; ○—○ mean of predictions from model B.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Proportion of infants converting in Idete, Tanzania. ○–○ = predictions from model B; ▪ = observed conversion proportion^{11} (from slide negative to slide positive) with 95% confidence intervals.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Proportion of infants converting in Idete, Tanzania. ○–○ = predictions from model B; ▪ = observed conversion proportion^{11} (from slide negative to slide positive) with 95% confidence intervals.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Proportion of infants converting in Idete, Tanzania. ○–○ = predictions from model B; ▪ = observed conversion proportion^{11} (from slide negative to slide positive) with 95% confidence intervals.
Citation: The American Journal of Tropical Medicine and Hygiene Am J Trop Med Hyg 75, 2_suppl; 10.4269/ajtmh.2006.75.2_suppl.0750011
Address correspondence to Thomas Smith, Swiss Tropical Institute, Socinstrasse 57, PO Box, CH-4002 Basel, Switzerland. E-mail: Thom-as-A.Smith@unibas.ch
Authors’ addresses: Thomas Smith, Nicolas Maire, Penelope Vounatsou, and Marcel Tanner, Swiss Tropical Institute, Socinstrasse 57, PO Box, CH-4002 Basel, Switzerland, Telephone: 41-61-284-8273, Fax: 41-61-284-8105, E-mails: Thomas-A.Smith@unibas.ch, nicolas.maire@unibas.ch, penelope.vounatsou@unibas.ch, and marcel.tanner@unibas.ch. Klaus Dietz, Department of Medical Biometry, University of Tübingen, Westbahnhofstrasse 55, 72070 Tübingen, Germany, Telephone: 49-7071-29-78253, Fax: 49-7071-29-5075, E-mail: klaus.dietz@uni-tuebingen.de. Gerry F. Killeen, Ifakara Health Research and Development Centre, Ifakara, Kilombero District, Tanzania, Telephone: 255-748-477-118, Fax: 255-23-262-5312, E-mail: gkilleen@ihrdc.or.tz. Louis Molineaux, Peney-Dessus, CH-1242 Satigny, Geneva, Switzerland.
Acknowledgments: We thank Paulette Rosé for assistance in retrieving the Garki data and Timothy Haley for assistance with statistical analyses. We also thank the members of the Technical Advisory Group (Michael Alpers, Paul Coleman, David Evans, Brian Greenwood, Carol Levin, Kevin Marsh, F. Ellis McKenzie, Mark Miller, and Brian Sharp), the Project Management Team at the Program for Appropriate Technology in Health (PATH) Malaria Vaccine Initiative, and GlaxoSmithKline Biologicals S.A. for their assistance.
Financial support: The mathematical modeling study was supported by the PATH Malaria Vaccine Initiative and GlaxoSmithKline Biologicals S.A.
Disclaimer: Publication of this report and the contents hereof do not necessarily reflect the endorsement, opinion, or viewpoints of the PATH Malaria Vaccine Initiative or GlaxoSmithKline Biologicals S.A.
REFERENCES
- 1↑
Ballou WR, Arevalo-Herrera M, Carucci D, Richie TL, Corradin G, Diggs C, Druilhe P, Giersing BK, Saul A, Heppner DG, Kester KE, Lanar DE, Lyon J, Hill AV, Pan W, Cohen JD, 2004. Update on the clinical development of candidate malaria vaccines. Am J Trop Med Hyg 71 (Suppl 2):239–247.
- 2↑
Hoffman SL, Goh LM, Luke TC, Schneider I, Le TP, Doolan DL, Sacci J, de la Vega P, Dowler M, Paul C, Gordon DM, Stoute JA, Church LW, Sedegah M, Heppner DG, Ballou WR, Richie TL, 2002. Protection of humans against malaria by immunization with radiation-attenuated Plasmodium falciparum sporozoites. J Infect Dis 185 :1155–1164.
- 3↑
Garcon N, Heppner DG, Cohen J, 2003. Development of RTS,S/AS02: a purified subunit-based malaria vaccine candidate formulated with a novel adjuvant. Expert Rev Vaccines 2 :231–238.
- 4
Bojang K, Milligan PJM, Pinder M, Vigneron L, Alloueche A, Kester KE, Ballou WR, Conway D, Reece WHH, Gothard P, Yamuah L, Delchambre M, Voss G, Greenwood BM, Hill A, McAdam KP, Tornieporth N, Cohen JD, Doherty T, 2001. Efficacy of RTS,S/AS02 malaria vaccine against Plasmodium falciparum infection in semi-immune adult men in The Gambia: a randomised trial. Lancet 358 :1927–1934.
- 5↑
Alonso PL, Sacarlal J, Aponte J, Leach A, Macete E, Milman J, Mandomando I, Spiessens B, Guinovart C, Espasa M, Bassat Q, Aide P, Ofori-Anyinam O, Navia MM, Corachan S, Ceuppens M, Dubois MC, Demoitie MA, Dubovsky F, Menendez C, Tornieporth N, Ballou WR, Thompson R, Cohen J, 2004. Efficacy of the RTS,S/AS02A vaccine against Plasmodium falciparum infection and disease in young African children: randomised controlled trial. Lancet 364 :1411–1420.
- 6↑
Rogier C, Trape JF, 1995. Study of premunition development in holo- and meso-endemic malaria areas in Dielmo and Ndiop (Senegal): preliminary results. Med Trop (Mars) 55 :71–76.
- 7↑
Hoffman SL, Oster CN, Plowe CV, Woollett GR, Beier JC, Chulay JD, Wirtz RA, Hollingdale MR, Mugambi M, 1987. Naturally acquired antibodies to sporozoites do not prevent malaria: vaccine development implications. Science 237 :639–642.
- 8↑
Beier JC, Oster CN, Onyango FK, Bales JD, Sherwood JA, Perkins PV, Chumo DK, Koech DV, Whitmire RE, Roberts CR, 1994. Plasmodium falciparum incidence relative to entomologic inoculation rates at a site proposed for testing malaria vaccines in western Kenya. Am J Trop Med Hyg 50 :529–536.
- 9↑
Owusu-Agyei S, Koram KA, Baird JK, Utz G, Binka F, Nkrumah F, Fryauff DJ, Hoffman SL, 2001. Incidence of symptomatic and asymptomatic Plasmodium falciparum infection following curative therapy in adult residents of northern Ghana. Am J Trop Med Hyg 65 :197–203.
- 10↑
Baird JK, Owusu AS, Utz G, Koram K, Barcus MJ, Jones TR, Fryauff DJ, Binka F, Hoffman SL, Nkrumah F, 2002. Seasonal malaria attack rates in infants and young children in northern Ghana. Am J Trop Med Hyg 66 :280–286.
- 11↑
Charlwood JD, Smith T, Lyimo E, Kitua A, Masanja H, Booth M, Alonso PL, Tanner M, 1998. Incidence of Plasmodium falciparum infection in infants in relation to exposure to sporozoite-infected anophelines. Am J Trop Med Hyg 59 :243–251.
- 12↑
Dietz K, Molineaux L, Thomas A, 1974. A malaria model tested in the African savannah. Bull World Health Organ 50 :347–357.
- 13↑
Nedelman J, 1985. Estimation for a model of multiple malaria infections. Biometrics 41 :447–453.
- 14↑
Struchiner CJ, Halloran ME, Spielman A, 1989. Modeling malaria vaccines. I: New uses for old ideas. Math Biosci 94 :87–113.
- 15↑
Molineaux L, Storey J, Cohen JE, Thomas A, 1980. A longitudinal study of human malaria in the west African savanna in the absence of control measures: relationships between different Plasmodium species, in particular P. falciparum and P. malariae. Am J Trop Med Hyg 29 :725–737.
- 16↑
Port GR, Boreham PFL, Bryan JH, 1980. The relationship of host size to feeding by mosquitos of the Anopheles gambiae Giles complex (Diptera, Culicidae). Bull Entomol Res 70 :133–144.
- 17↑
Carnevale P, Frezil JL, Bosseno MF, Lepont F, Lancien J, 1978. Study of agressivity of Anopheles gambiae A in relation to age and sex of human subjects. Bull World Health Organ 56 :147–154.
- 18↑
Smith T, Killeen G, Lengeler C, Tanner M, 2004. Relationships between the outcome of Plasmodium falciparum infection and the intensity of transmission in Africa. Am J Trop Med Hyg 71 (Suppl 2):80–86.
- 20↑
Smith T, Killeen G, Maire N, Ross A, Molineaux L, Tediosi F, Hutton G, Utzinger J, Dietz K, Tanner M, 2005. Mathematical modeling of the impact of malaria vaccines on the clinical epidemiology and natural history of Plasmodium falciparum malaria: overview. Am J Trop Med Hyg 75 (Suppl 2):1–10.
- 21↑
Keiser O, 2001. Einfluss von Morbidität auf Kinderwachstum in Papua Neu Guinea und Tansania. Basel, Swiss Tropical Institute. Diploma Thesis.
- 25↑
Maire N, Smith T, Ross A, Owusu-Agyei S, Dietz K, Molineaux L, 2005. A model for natural immunity to asexual blood stages of Plasmodium falciparum in endemic areas. Am J Trop Med Hyg 75 (Suppl 2):19–31.
- 26↑
Kirkpatrick S, Gelatt CD Jr, Vecchi MP, 1983. Optimization by simulated annealing. Science 220 :671–680.
- 27↑
Press WH, Flannery BP, Teukolsky SA, Vetterling WT, 1988. Numerical Recipes in C: the Art of Scientific Computing. Cambridge, United Kingdom: Cambridge University Press
- 28↑
Smith T, Ross A, Maire N, Rogier C, Trape JF, Molineaux L, 2005. An epidemiologic model of the incidence of acute illness in Plasmodium falciparum malaria. Am J Trop Med Hyg 75 (Suppl 2):56–62.
- 29↑
Baird JK, Masbar S, Basri H, Tirtokusumo S, Subianto B, Hoffman SL, 1998. Age-dependent susceptibility to severe disease with primary exposure to Plasmodium falciparum. J Infect Dis 178 :592–595.
- 30↑
Smith T, Beck H-P, Kitua A, Mwankusye S, Felger I, Fraser-Hurt N, Irion A, Alonso PL, Teuscher T, Tanner M, 1999. Age dependence of the multiplicity of Plasmodium falciparum infections and of other malariological indices in an area of high endemicity. Trans R Soc Trop Med Hyg 93 (Suppl 1):15–20.
- 31↑
Baird JK, 1998. Age-dependent characteristics of protection v. susceptibility to Plasmodium falciparum. Ann Trop Med Parasitol 92 :367–390.
- 32↑
Boyd MF, Kitchen SF, 1937. Observations on induced falciparum malaria. Am J Trop Med 17 :213–235.
- 33↑
Kester KE, McKinney DA, Tornieporth N, Ockenhouse CF, Heppner DG, Hall T, Krzych U, Delchambre M, Voss G, Dowler MG, Palensky J, Wittes J, Cohen J, Ballou WR, 2001. Efficacy of recombinant circumsporozoite protein vaccine regimens against experimental Plasmodium falciparum malaria. J Infect Dis 183 :640–647.
- 34↑
Macdonald G, 1950. The analysis of malaria parasite rates in infants. Trop Dis Bull 47 :915–938.
- 35↑
Pull JH, Grab B, 1974. A simple epidemiological model for evaluating the malaria inoculation rate and the risk of infection in infants. Bull World Health Organ 51 :507–516.
- 36↑
Nedelman J, 1984. Inoculation and recovery rates in the malaria model of Dietz, Molineaux and Thomas. Math Biosci 69 :209–233.
- 37↑
Knols BG, De Jong R, Takken W, 1995. Differential attractiveness of isolated humans to mosquitoes in Tanzania. Trans R Soc Trop Med Hyg 89 :604–606.
- 38↑
Mukabana WR, Takken W, Killeen GF, Knols BG, 2004. Allomonal effect of breath contributes to differential attractiveness of humans to the African malaria vector Anopheles gambiae. Malar J 3 :1.
- 39↑
Schofield L, Ferreira A, Altszuler R, Nussenzweig V, Nussenzweig RS, 1987. Interferon-gamma inhibits the intrahepatocytic development of malaria parasites in vitro. J Immunol 139 :2020–2025.
- 40↑
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, 2000. Alpha-galactosylceramide-activated Valpha 14 natural killer T cells mediate protection against murine malaria. Proc Natl Acad Sci USA 97 :8461–8466.
- 41↑
James SP, 1931. Some general results of a study of induced malaria in England. Trans R Soc Trop Med Hyg 24 :477–538.
- 42↑
Rickman LS, Jones TR, Long GW, Paparello S, Schneider I, Paul CF, Beaudoin RL, Hoffman SL, 1990. Plasmodium falciparum-infected Anopheles stephensi inconsistently transmit malaria to humans. Am J Trop Med Hyg 43 :441–445.
- 43↑
Kitua A, Smith T, Alonso PL, Masanja H, Urassa H, Menendez C, Kimario J, Tanner M, 1996. Plasmodium falciparum malaria in the first year of life in an area of intense and perennial transmission. Trop Med Int Health 1 :475–484.