|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
Community effects of vaccines or of other interventions targeting the human host (herd immunity) result either directly or indirectly from changes in infectivity. To understand these, mathematical models are needed. Infectivity of the human host is a key determinant of the intensity of malaria transmission, but in contrast to field data, models of malaria transmission dynamics have generally assumed that the probability of transmission to a vector is equal across all infected humans and throughout the course of infections,10 or that complete transmission-blocking immunity is acquired as a result of repeated exposure.11
As one component of a project to develop simulation models for the potential impact of a malaria vaccine, we now develop a new model for infectivity. We avoid directly incorporating the poorly understood processes of gametocytogenesis and the subsequent fate of gametocytes. Instead, we implicitly model gametocytogenesis by estimating infectivity from analyses of the association between lagged asexual blood stage densities and infectivity. We use data from deliberate infection of human subjects with Plasmodium falciparum between 1940 and 1963 as a treatment of tertiary syphilis7,12 and fit a statistical model to these data to estimate the relationship between recent asexual parasitemia and infectivity.
| MATERIALS AND METHODS |
|
|
|---|
(i,t) as the density of asexual parasites in patient i that gave rise to the gametocytes present at time t. Our stochastic simulation model of malaria epidemiology13,14 uses five-day time steps. To incorporate transmission to the vector into this model, we propose that
(i,t) should be estimated as a weighted sum of the asexual parasite densities measured 10, 15, and 20 days prior to the feeding experiment. Thus if ß1, ß2 and ß3 are constants, then:
![]() |
This quantity is constrained to be zero when there is no history of patent infection in the last 20 days, or when patent infection first arose in the last 5-day period. We justify the choice of the lag periods of 10, 15, and 20 days by observing that previous analyses of gametocyte dynamics in the same malariatherapy patients suggested that gametocytes spend an average of 7.4 days sequestered before appearing in the peripheral circulation and have a subsequent mean circulation time of 6.4 days.15
We considered a number of models for infectivity. In our preferred model, the ratio of functional female gametocytes per blood meal to
(i,t) is log normally distributed, with some geometric mean,
, so that the density of functional female gametocytes in the host blood, yg(i,t), is related to
(i,t) by
![]() |
For a mosquito to become infected with P. falciparum, it must imbibe at least one female and one male gametocyte in the same blood meal. We assume that functional gametocytes are those able to give rise to sporozoites. We define yg* to be the density of female gametocytes necessary for infection of the mosquito, corresponding to 1 functional gametocyte per blood meal. The probability that at least one functional female gametocyte is taken up is then Pr(yg(i,t) > yg*) and it follows that:
![]() |
where
is the percentile of the cumulative standard normal distribution and
* = (ln(
) ln(yg* ))/
g is a constant which depends on the volume of the blood meal, the viability of the gametocytes, and the variability in the system. If the volume of the blood meal is assumed to be 3 µL16 (and there is assumed to be no concentration of erythrocytes in these species17,18), then the ratio of the functional female gametocyte density to
(i,t) is given by
/3.
If the sex ratio of functional male to female gametocytes is 1:1, then the probability that at least one functional male gametocyte is taken up is also Pr(yg(i,t) > yg*), and the probability that a mosquito is infected with both male and female gametocytes is Pr(yg(i,t) > yg*)2. The number of mosquitoes out of a batch of nfed(i,t) mosquitoes fed on patient i, at time t, that are infected follows a Binomial distribution
![]() |
Our estimate of the proportion of mosquitoes feeding on individual i at time point t that would become infected, Im(i,t), is then the expected proportion from this binomial
![]() |
Data and model fitting. We fitted the model given by equations 15 to archive data provided by Dr. W. Collins (Centers for Disease Control and Prevention, Atlanta, GA) collected by the United States Public Health Service in South Carolina and Georgia between 1940 and 1963. At that time, malaria therapy was a recommended treatment of neurosyphilis.12 We obtained archive data for 392 patients inoculated with varying species, strains, and parasite stages of Plasmodium. In each case, densities of asexual parasitemia and of gametocytemia were recorded on a daily basis. Batches of between 5 and 60 caged Anopheles quadrimaculatus or An. albimanus mosquitoes were allowed to feed on gametocytemic patients and then kept to allow oocysts to develop in the midgut before dissection to record infection status. Details of the methods have been previously published.12
We excluded patients inoculated or co-inoculated with Plasmodium species other than P. falciparum and those who were inoculated for a second time or who were treated with any of the antimalarial drugs used and those with insufficient data. We initially explored various possible predictors of the proportion of the mosquitoes that became infected. These comprised day since first parasitemia, fever, the concurrent asexual parasite density, and the history of parasitemia for the preceding 40 days. To maintain compatibility with our population model of asexual parasitemia,13 we simplified the history by considering parasite density only at five-day intervals.
The model described by equations 15 was selected as the most appropriate on the basis of these exploratory analyses. We treated it as a hierarchical Bayesian model, and fitted it using the Metropolis-Hastings algorithm in the program WinBUGS version 1.4.19 ß1 was fixed at a value of 1 and we used imprecise log normal prior distributions for ß2 and ß3, and a gamma prior for 1/
g2
Age-pattern of infectivity. Field estimates of the relative contribution to the infectious reservoir by different age groups have been made for settings in Liberia,9 The Gambia,8 Tanzania,8 western Kenya,20 Papua New Guinea,21 and Cameroon.22 The estimates from Africa are based on feeding of insectary An. gambiae, while those from Papua New Guinea used a colony of An. farauti.
To validate our model against these data, we simulated the epidemiology of malaria in field settings with entomologic data using our stochastic simulation model.13,14,23 This model makes predictions of parasite densities using a five-day time step assuming an annually recurring stable pattern of transmission as input. We used this model to implement an individual-based simulation of a population of at least 10,000 simulated individuals exposed from birth to the local transmission patterns. For each individual and time point during the year we use equations 15 to make predictions of Im(i,t).
There is evidence that the risk of being bitten by an anopheline is approximately proportional to body size.24 To allow for this, we make an estimate of the contribution to the infectious reservoir for age group j, weighted proportionately to A(a(i,t)), the estimated body surface area of the host, where a(i,t) is the age of individual i at time t.14 The proportion of infectiousness contributed by group j is then
![]() |
where Jj(i,t) = 1 if individual i is in age group j at time t, and Jj(i,t) = 0 otherwise, and where the summations are over the whole simulated population and a whole year of follow-up time.
| RESULTS |
|
|
|---|
A total of 562 (77%) blood slides were positive for P. falciparum 10 days prior to the feeding experiments, with median density of 1,792 parasites/µL (range = 3025,990). Parasite prevalence and densities on the other days analyzed were comparable. A total of 565 (77%) of the 730 samples on the day of the feeding experiment were gametocytemic, with median gametocyte densities of 190 (range = 103,330).
Overall, 5,420 (24%) of 22,431 mosquitoes became infected. However, there was considerable variation in the observed proportions of mosquitoes infected in any one experiment (Figure 1
). For the simulation model, we aimed to describe the average relationship between asexual parasite density and infectivity as accurately as possible for all values of the parasitemia history.
|
Model for infectivity.
We rejected a number of models for infectivity before selecting the model given by equations 15. Fixed effects models that ignored variation between patients or samples in the relationship predicted very high proportions of infected mosquitoes for hosts with a history of high parasite densities. There is substantial variation in infectivity in the data even at high values of the weighted average of asexual parasitemia
(i,t) (Figure 1
), and models with fixed effects of the weighted sum of asexual parasitemia gave a better fit only if there was allowance for overdispersion in infectivity. However, allowing for overdispersion in the response led to problems of convergence, with the estimates of random effects highly correlated with those of
(i,t).
Our preferred model (equations 15) introduces random variation into the explanatory variable,
(i,t), intended to simulate random variation in the process of gametocytogenesis. The fit was good (Figure 2
and Table 1
) except at very high values of the weighted sum of asexual density (
(i,t)). There is little data at these high values, and although proportions of 100% infected mosquitoes were observed in almost one-third of the experiments with the highest 5% of values of
(i,t), there was also considerable variation, with some proportions as low as 17%. This variation coupled with the relatively small amount of data makes it difficult to specify a model that fits well in this range.
|
|
![]() |
To test this assumption of log-normality, we plotted the logarithm of the ratio of the observed density of female gametocytes yf(i,t) to the weighted sum of the asexual parasite densities (Figure 3a
). The observed distribution of this ratio is somewhat truncated because the limit of detection in the malaria therapy dataset was 10 gametocytes/µL. However, the normal quantile plot (Figure 3b
) suggests that the log-normal approximation is roughly valid, with a mean of approximately 5 for ln(yf(i,t)/
(i,t)).
|
(i,t)exp(5), and the expectation of the normal distribution in equation 6 is
![]() |
Using our best estimate of
(Table 2
), this gives an approximation for the proportion of gametocytes that are functional of yg(i,t)/yf(i,t)
exp(5)
0.5.
|
The studies of Drakeley and others8 investigated infectivity using feeds only on blood from demonstrable gametocyte carriers, and this makes it difficult to be confident in comparing the results with our population-based estimates. The assumption that individuals without detectable gametocytes do not transmit to mosquitoes is questionable. However, the estimated contributions of different age-groups were similar to those for Ifakara in settings where population-based assessments have been carried out. The general pattern in our simulations was in agreement with the data: at higher transmission (as in Kano, Kenya) the contribution of the younger age groups was increased, and at lower transmission (as in Farafenni, or in Papua New Guinea) the older age groups are predicted to be relatively more important.
| DISCUSSION |
|
|
|---|
The latter simplification allowed us to avoid any strong assumptions about the problematical relationship between observed and functional gametocytemia or about the poorly understood factors affecting gametocytogenesis.25,26 Our examination of the distribution of gametocyte densities confirmed that our model for gametocyte densities is a reasonable approximation. Our estimate for the percentage of gametocytes that are functional is low, in agreement with estimates from life tables,27,28 which indicate that many parasites are killed between the macrogamete and oocyst stages.
We assumed a 50:50 ratio of functional male:female gametocytes. The median proportion of male gametocytes recorded in the malariatherapy dataset was 0.43 with a 90% central range of 01, while other studies have found lower proportions (approximately 0.25) of male gametocytes.2931 However, male gametocytes produce more gametes than female gametocytes,25 and thus perhaps yg* for male gametocytes should be lower than for female gametocytes to compensate for any imbalance in the sex ratio. We considered the sensitivity of Im(i,t), the proportion of mosquitoes feeding on individual i at time t who became infected, to variations in the ratio of functional male:female gametocytes. Im(i,t) showed little sensitivity to variations in the assumed sex-ratio of gametocytes over a wide range. This supports the claim that the sex ratio of gametocytes is unlikely to be of major importance in determining infectivity;25 and studies to date have shown little or no effect of sex ratio on infectivity.3235 We do not consider heterogeneity between individuals or parasite clones in their innate propensity to produce gametocytes because our interest focuses on overall infectivity. Moreover, we ignore any direct effects of fever on infectivity36 and of drug treatment on gametocytemia. These effects are likely to be of epidemiologic importance mainly in low transmission areas.
To apply our model to endemic areas in Africa, we must make several further simplifications. The malariatherapy experiments used An. quadrimaculatus and An. albimanus in controlled conditions, which differ from African field conditions where the main vectors are An. funestus and the An. gambiae complex. We consider the extent to which our approach leads to an overall bias in the estimates of the proportion of mosquito bites resulting in infections of the vector in a separate paper.37
We also assume that the relationship between asexual parasitemia and infectivity in adults infected for the first time can be generalized to those with varying levels of acquired immunity, which entails assuming that there is no acquisition of transmission-blocking immunity and no effect of antibodies to gametocytes. This is supported by our finding that the model predicts age patterns that concur with those from field data. Although children have often been thought to contribute most of the malaria transmission to vectors because of their higher parasite densities, those studies that have estimated the contributions of different ages of hosts to the infectious reservoir have found that adults make a substantial contribution (Table 2
).
Our model supports the conclusions of these field studies and indeed suggests that the contribution of adults has been underestimated because of the usual assumption that each potential host is bitten equally frequently. This corroborates the assumption that transmission-blocking immunity does not markedly increase with age, and that natural acquired immune responses to gametocytes are also not of epidemiologic importance. However, these conclusions rest on a limited evidence base. There is a need for improved population-based estimates of age- and exposure- dependence in the infectivity of humans to mosquitoes in endemic areas.
Received September 18, 2005. Accepted for publication November 5, 2005.
Acknowledgments: We thank Professor Klaus Dietz and Dr. Louis Molineaux for comments on an earlier draft of this paper. 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 GlaxoSmith-Kline Biologicals S.A. for their assistance. We gratefully acknowledge the provision of the malariatherapy data from the United States Public Health Service through Bill Collins and Geoffrey Jeffery.
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.
* Address correspondence to Amanda Ross, Swiss Tropical Institute, Socinstrasse 57, Postfach, 4002 Basel, Switzerland. E-mail: amanda.ross{at}unibas.ch ![]()
Authors addresses: Amanda Ross and Thomas Smith, Swiss Tropical Institute, Socinstrasse 57, PO Box, CH-4002, Basel, Switzerland, Telephone: 41-61-284- 8273, Fax: 41-61-284-8105, E-mails: amanda.ross{at}unibas.ch and Thomas-A.Smith{at}unibas.ch. Gerry Killeen, Ifakara Health Research and Development Center, PO Box 53, Ifakara, Kilombero District, Tanzania, Telephone: 255-748-477-118, Fax: 255-23-262-5312, E-mail: gkilleen{at}ihrdc.or.tz.
Reprint requests: Amanda Ross, 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] |
||||
![]() |
G. F. KILLEEN, A. ROSS, and T. SMITH INFECTIOUSNESS OF MALARIA-ENDEMIC HUMAN POPULATIONS TO VECTORS. Am J Trop Med Hyg, August 1, 2006; 75(2_suppl): 38 - 45. [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, 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] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |