|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
There is abundant evidence that most clinical episodes are caused by newly inoculated genetically distinct parasites.8,9 One proposed model is that parasite populations are structured into a limited number of strains, each stimulating long-term clinical immunity.10,11 However most analyses of the population biology of P. falciparum have concluded that there is frequent genetic exchange,1215 many malaria antigens are extremely polymorphic,1618 cross-protection is clearly important, and natural immunity to the immunodominant epitopes is not necessarily lifelong.1921
The adequate modeling of all these complex immunologic phenomena represents a major challenge. However epidemiologic analyses of the tolerance of parasites can be used to predict the likelihood of clinical episodes as a function of densities of peripheral parasitemia without explicitly considering how those densities occur.2225 In a study carried out in Dielmo, where parasitemia was assessed twice weekly Rogier and others estimated well-defined pyrogenic thresholds for different ages of human host.24 We have now further analyzed these data to derive predictions of the thresholds as functions of recent levels of parasitemia, rather than of the age of the host. We have linked these predictions to a stochastic model that predicts parasite densities in endemic areas as a function of the pattern of transmission26,27 and fitted the model for the incidence of clinical episodes to field data from different epidemiologic settings in Ndiop, Dielmo, and Idete. The resulting model enables us to predict, for a wide range of malaria transmission settings, the occurrence of clinical episodes and to assess the likely effects of interventions on the incidence of clinical attacks.
| MATERIALS AND METHODS |
|
|
|---|
Model for clinical malaria episodes. The parasitologic simulation includes stochastic variation between individual humans in average parasite densities and also stochastic variation around that average.26 We model clinical immunity as a function of these stochastically varying parasite densities, and of a set of five parameters that are independent of the individual and of the transmission setting.
To predict the clinical outcome, for each five-day time step we draw five independent samples from the simulated parasite density distribution for each concurrent infection (to simulate potential daily changes in morbidity status) and consider only the maximum, Ymax(i,t), of the simulated densities to determine whether a clinical episode occurred. When the host is infected by several concurrent infections it is likely that one of these contributes the bulk of the parasite load, so it is logical to define Ymax(i,t) as the maximum over all infections.
A simple model is to assume that for each host there is a specific parasite density, or pyrogenic threshold, at which symptoms (e.g., fever) are triggered. Rogier and others24 considered a cohort of the inhabitants of the holoendemic village of Dielmo, Senegal and fitted a step function to the probability of fever as a function of parasite density. The parasite density at which the step occurs corresponds to the pyrogenic threshold, which was shown to vary with age.
In general, it is not realistic to assume that all individuals of the same age will have exactly the same pyrogenic threshold,2831 so it is more reasonable to expect a sigmoidal relationship between the risk of fever and the parasite density than a step function. We therefore propose a model in which the probability that an episode occurs in individual i, at time t, is related to the parasite density via a function of the following form:
![]() |
where Y*(i,t), the pyrogenic threshold for individual i at time t, is defined as the parasite density at which the probability of a clinical episode reaches 0.5, and Ymax(i,t) is the maximum density during the time interval t (note that we present only the formulae for our final choice of models).
The age pattern in the pyrogenic threshold in Dielmo, together with data derived from other study sites,29,32,33 supports the idea that the density of parasites required to stimulate acute pathology is higher in individuals who have been recently exposed to high parasite densities. This may be a result of stimulation of immune responses to toxins released at schizogony, and very likely involves physiologic tolerance of cytokines.34 The mechanism must be consistent with both rapid acquisition and rapid loss of tolerance and cannot be a simple function of antibody against toxin, which have a completely different age-pattern from that of the pyrogenic threshold.35 We model the dynamics of the pyrogenic threshold with a function of the form
![]() |
where f1(Y(i,t)) is a function describing the relationship between accrual of tolerance and the parasite density Y(i,t); f2(Y*(i,t)) describes saturation of this accrual process at high values of Y*, and the term
Y*(i,t) leads to decay of the threshold with first-order kinetics. The decay ensures that the model conforms to the epidemiologic evidence suggesting that parasite tolerance is short lived.
We define the function f1(Y(i,t)) in such a way as to ensure that the stimulus is not directly proportional to Y but rather that it asymptotically reaches a maximum at high values of Y, using
![]() |
To ensure saturation of the accrual process, we require that at high values of Y*, a higher parasite load is required to achieve the same increase by defining
![]() |
Overall therefore we propose the following dynamics for Y*:
![]() |
where
Y1* , and Y2* are constants to be estimated. To complete the specification of the model, we set the initial conditions to be Y* (i,0) = Y0* at the birth of the host, thus defining a further parameter Y0*
Data sources. We fitted the model for acute episodes to two distinct datasets. The first was published data on the age pattern of clinical episodes in the villages of Ndiop and Dielmo in Senegal.2 The village populations were visited daily to detect and treat any clinical malaria attacks (with quinine). Thus, effectively all acute episodes were thought to be treated in these villages. In the simulations of Dielmo and Ndiop we assumed that there had been no treatment of clinical malaria prior to the start of the follow-up period. To ensure that the analysis remains tractable, we approximate the patterns of transmission with recurring annual cycles (although there was variation between years in the predominant vectors and seasonality of transmission).
We also compared the predicted patterns from the simulation model for P. falciparum parasitology with those from parasitologic surveys in these two villages to evaluate its appropriateness as a basis for the predicting clinical episodes in this setting. In Dielmo, two thick blood smears were prepared per week for each individual from May 29, 1990 to September 30, 1990. In Ndiop, one thick blood smear was prepared per week for each individual from July 15, 1993 to January 15, 1994 and one per month from January 2, 1994 to July 15, 1995. Slides were only declared negative after 200 high-power fields had been scanned for parasites. Parasite densities were originally expressed as the parasite:leukocyte ratio. To adjust these densities to the same scale as that used in fitting the simulation model to other datasets, the parasite:leukocyte ratios were then multiplied by a factor of 1,416 to give a notional density in parasites/microliter of blood.26
The model was fitted to a second dataset of age-incidence rates for clinical malaria in infants less than one year of age recorded at the health center in the village of Idete, Tanzania from June 1993 to October 1994.6,25 These data were included to estimate the initial conditions (the value of Y0*) and to ensure that the model predicts the age pattern of acute episodes that is actually observed in infants. For the Idete data we used the case definitions and age groups in the paper by Vounatsou and others25 and the annual pattern of inoculations reported by Charlwood and others36 as input. We assume a common value of Y0* across all sites and therefore require data for infants from only a single transmission setting.
Implementation and fitting of the simulation model.
To obtain estimates of the five parameters
,
, Y0*, Y1*, and Y2* we fitted the model to the age-pattern of clinical malaria in all three villages (i.e., Ndiop and Dielmo in Senegal, and Idete in Tanzania) and simultaneously to the pyrogenic thresholds for Dielmo estimated by Rogier and others24 (Table 1
).
|
Simulated clinical episodes of malaria occurred with probability Pm(i,t), which was dependent on both the simulated maximum density and the current value of Y*(i,t) for each individual and each five-day time point in the 10-year follow-up period. In the simulations of Ndiop and Dielmo we simulated effective treatment of all clinical episodes within the five-day period in which they occurred. In the simulation of Idete we assumed that some proportion, Pt, of the episodes were effectively treated (i.e., the parasites were cleared within the course of one time interval), and that this proportion corresponded to the proportion of episodes reported to the village dispensary. In Idete village, simulated episodes occurring within 30 days of a preceding episode were not counted (these have been registered in the surveillance system as recrudescence, rather than new episodes). In Ndiop and Dielmo this restriction did not hold.38
For each simulated individual in each village the model thus predicted the incidence of clinical malaria as a stochastic function of the inoculation rate. These incidences were summarized over age groups and compared with the published values.2,25 Similarly, the model predicted the pyrogenic threshold, Y*(i,t), at each time point for each individual. The geometric mean of these values was calculated for each age group in the simulation of Dielmo village, and the logarithms of these values compared with the logarithms of the age-specific pyrogenic thresholds estimated by Rogier and others.24 Simulated annealing39,40 was used to identify the parameter values that minimized the residual sum of squares summed over all three villages and both outcomes for Dielmo. The Fisher information estimated from a least squares quadratic fit to the residual sum of squares was used to give approximate confidence intervals.
| RESULTS |
|
|
|---|
|
, estimated as 2.5/year, implied that in the absence of stimulation, the pyrogenic threshold decays with a half life of 0.33 years. The predicted total numbers of episodes up to age 60 were 56 for Ndiop (EIR = 20) and 53 for Dielmo (EIR = 200) compared with the published overall incidence of clinical malaria cumulative numbers of episodes up to the age of sixty of 62 and 43, respectively.2 In the simulations of both villages the age of peak incidence was a little younger than in the data predicted incidence and incidence was rather lower in the youngest individuals than the observed values, and higher in adults. The extent of the peak shift was similar in the model to the data.
Although the model was not fitted to the patterns of age prevalence and of age density in Dielmo or Ndiop, it does make predictions of these quantities, which we could therefore compare with the observed curves. The predicted age-prevalence curve for Dielmo was very similar to that observed (Figure 2a
), as were the predicted geometric mean densities in children in that village (Figure 3a
). In adults the model predicted rather higher densities than those observed in Dielmo, while for adults in Ndiop the model predicted higher prevalence in adults (Figure 2b
) but lower densities (Figure 3b
) than those observed. This would be expected if the burden of malaria is concentrated in a smaller proportion of individuals in Ndiop than in the dataset to which the parasitologic model was fitted. A reasonably good fit was obtained for the average pyrogenic threshold, but the model did not give a very good fit to the age-trend in Y*(i,t), predicting that the peak was at a greater age than the estimates of Rogier and others24 (Figure 4
).
|
|
|
| DISCUSSION |
|
|
|---|
This good fit was obtained despite the use of a parasitologic model that only crudely reproduces within-host parasite dynamics, since we fitted it to cross-sectional data.26 Day-to-day variation in parasite densities may be critical in determining levels of tolerance, and our model, based on five-day time steps, did not aim to simulate this accurately. This may explain why the density of patent parasitemia did not appear to be very important, and may also be the explanation of why we could not obtain a better fit for the age-pattern of the pyrogenic threshold. It is possible that the important variations in density and levels of tolerance are much more rapid than our model could capture, especially if they involve physiologic tolerance of cytokines.34
The relatively poor prediction of parasite prevalence and density in adults in Ndiop is possibly because the model assumes the degree of within-village heterogeneity in transmission to be the same in each village. Focality of transmission leads to lower prevalence, but higher densities in those who are infected, because of increased levels of superinfection. Based on these criteria, transmission in Ndiop appears to be more focal than that in the villages to which we fitted the parasitologic model.26 Within our model effects of focal transmission on incidence of clinical episodes should be only of secondary importance because there is little interaction between concurrent co-infections. Thus, at any given level of immunity the incidence of clinical episodes depends primarily on the overall force of infection and not on how the infections are distributed between individuals.
It was not possible to obtain a better fit for infants in Ndiop because the number of infection events predicted for this age group by our model of infection27 is less than the number of clinical episodes. We have assumed all episodes are immediately treated so that no more than one episode can occur for any one infection event, but this was not necessarily always the case. We assumed mosquito biting to be proportional to body surface area, using Tanzanian anthropometric data to estimate age-specific surface areas.41 Different patterns of human growth or mosquito behavior may account for some of the discrepancies. Selection effects that might arise because of differential mortality of susceptible individuals are an additional factor that we did not take into account.
Our model assumes particular functional forms for the relationships between the pyrogenic threshold and the risk of clinical episodes and the pyrogenic threshold and the parasite density itself. Exploratory analyses indicated that the fitted age-incidence relationships are not very sensitive to the exact functional forms used for these relationships. Empirical relationships between parasite density and risk of illness depend on how the cases are detected. In the studies in Senegal parasitemia and fever were monitored daily, so episodes were generally detected early and this may account for the abrupt pyrogenic thresholds reported by Rogier and others.24 More usually, fever cases are detected when they report to a health facility, as in the study in Idete.6 The arrival of the cases at the health facility is at varying intervals after the beginning of the episode and this tends to blur the relationship between fever risk and parasite density. If fever episodes are detected at household visits, which are carried out at intervals of more than a few days at times that are unrelated to the onset of disease, then the relationship between parasite densities and fever risk is weaker (e.g.,22,42).
In areas endemic for P. falciparum malaria, the incidence of clinical attacks is highly age dependent, with the peak incidence occurring at younger ages the higher the transmission. Such peak shifts are not only characteristic for malaria but also for many other infectious diseases.43 A superficially similar shift is also seen in patterns of age prevalence for P. falciparum,44 but the peak in prevalence is generally reached at an older age than that of acute morbidity. Unlike the pattern for clinical episodes, reduction in transmission is associated with reduction in infection prevalence over almost all of the age range, (although in some age groups there may be a small increase). The peak shift in clinical attacks is more pronounced than that in prevalence, and the incidence of acute malaria attacks in older children and adults can be substantially greater at low transmission levels than at high ones (Figure 1
).
The observation that reduction in transmission may lead to an increased incidence of disease in P. falciparum has been attributed to the phenomenon of endemic stability observed with many veterinary pathogens.45 For endemic stability to occur there must be at least two processes accounting for the age-incidence curves, one leading to an increase in incidence with age in the youngest age groups, and the other to a decreasing incidence in older individuals. Coleman and others45 suggest that the first of these conditions must be satisfied by a worsening of the outcome of infection with age over at least part of the age range. Our model demonstrates that this assumption is not necessary, for we explain the initial increase in morbidity with age as a consequence of increase in exposure to mosquitoes as the host grows in body surface area.27 Idete, Dielmo, and Ndiop are all villages with stable endemic P. falciparum. Although Ndiop is an example of mesoendemicity, contrasting with the holoendemic transmission in Idete and Dielmo, it still experiences a much higher EIR than areas of unstable transmission. The theory of endemic stability therefore needs adapting for the analysis of the case of endemic malaria.
We propose to use the sub-model of equations 15 as part of a comprehensive model for examining the likely consequences of a wide range of interventions, including vaccination. The incidence of acute illness is only one of these consequences, which can include severe life-threatening disease, chronic anemia, and indirect mortality. Even an intervention that leads to an increase in the incidence of uncomplicated illness in some age groups might lead to a reduction in mortality or severe disease. We know that parasite tolerance and anti-parasitic immunity have different dynamics, and conjecture that they make differential contributions to uncomplicated and severe disease, respectively.
Received September 18, 2005. Accepted for publication February 7, 2006.
Acknowledgments: We thank Professor Klaus Dietz for his helpful comments on an earlier draft of this report and Dr. Timothy Haley for assistance with statistical analyses.
Financial support: The mathematical modeling study was supported by the Program for Appropriate Technology in Health (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.
* Address correspondence to Thomas Smith, Swiss Tropical Institute, Socinstrasse 57, PO Box, CH-4002, Basel, Switzerland. E-mail: Thomas-A.Smith{at}unibas.ch ![]()
Authors addresses: Thomas Smith, Amanda Ross and Nicolas Maire, 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{at}unibas.ch, amanda.ross{at}unibas.ch, and nicolas.maire{at}unibas.ch. Christophe Rogier, Institut de Médecine Tropicale du Service de Santé des Armées, BP46, Parc du Pharo, 13998 Marseille-Armées, France; Telephone: 33-4-91-15-01-50/52, Fax: 33-4-91-01-64, E-mail: christophe.rogier{at}wanadoo.fr. Jean-François Trape, Institut de Recherche pour le Développement, 77 Paludologie Afrotropicale, BP 1386, CP 18524 Dakar, Senegal, Telephone: 221-849-3582, Fax: 221-832-4307, E-mail: Jean-Francois.Trape{at}ird.sn. Louis Molineaux, Peney-Dessus, CH-1242 Satigny, Geneva, Switzerland.
Reprints requests: Thomas Smith, Swiss Tropical Institute, Socinstrasse 57, PO Box, CH-4002, Basel, Switzerland.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
T. SMITH, G. F. KILLEEN, N. MAIRE, A. ROSS, L. MOLINEAUX, F. TEDIOSI, G. HUTTON, J. UTZINGER, K. DIETZ, and M. TANNER 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, August 1, 2006; 75(2_suppl): 1 - 10. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. ROSS, N. MAIRE, L. MOLINEAUX, and T. SMITH AN EPIDEMIOLOGIC MODEL OF SEVERE MORBIDITY AND MORTALITY CAUSED BY PLASMODIUM FALCIPARUM. Am J Trop Med Hyg, August 1, 2006; 75(2_suppl): 63 - 73. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. TEDIOSI, N. MAIRE, T. SMITH, G. HUTTON, J. UTZINGER, A. ROSS, and M. TANNER AN APPROACH TO MODEL THE COSTS AND EFFECTS OF CASE MANAGEMENT OF PLASMODIUM FALCIPARUM MALARIA IN SUB-SAHARAN AFRICA. Am J Trop Med Hyg, August 1, 2006; 75(2_suppl): 90 - 103. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. MAIRE, J. J. APONTE, A. ROSS, R. THOMPSON, P. ALONSO, J. UTZINGER, M. TANNER, and T. SMITH MODELING A FIELD TRIAL OF THE RTS,S/AS02A MALARIA VACCINE. Am J Trop Med Hyg, August 1, 2006; 75(2_suppl): 104 - 110. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. MAIRE, F. TEDIOSI, A. ROSS, and T. SMITH PREDICTIONS OF THE EPIDEMIOLOGIC IMPACT OF INTRODUCING A PRE-ERYTHROCYTIC VACCINE INTO THE EXPANDED PROGRAM ON IMMUNIZATION IN SUB-SAHARAN AFRICA. Am J Trop Med Hyg, August 1, 2006; 75(2_suppl): 111 - 118. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. TEDIOSI, G. HUTTON, N. MAIRE, T. A. SMITH, A. ROSS, and M. TANNER PREDICTING THE COST-EFFECTIVENESS OF INTRODUCING A PRE-ERYTHROCYTIC MALARIA VACCINE INTO THE EXPANDED PROGRAM ON IMMUNIZATION IN TANZANIA. Am J Trop Med Hyg, August 1, 2006; 75(2_suppl): 131 - 143. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |