|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
The entomologic inoculation rate (EIR), defined as the number of infective mosquitoes biting each human host in a unit time interval, is a key measurement of the exposure of human populations to malaria. A number of studies have considered how clinical outcomes relate to the EIR, since this determines the public health impact of any intervention that works by reducing malaria transmission. In this report, we focus on two such studies, an intensive study of clinical malaria in Senegal1 and a literature survey on the relationship between mortality in infants and children less than five years old and EIR.2 We do not investigate here the relationship between transmission intensity and severe malaria disease because of potential confounding by factors such as access to treatment.
In inter-age and inter-population comparisons, the EIR can be a poor indicator of parasitologic challenge for two reasons. First, the standard method for determining the EIR is to use human landing catches of mosquitoes made by adult male mosquito collectors. Since the biting rate of malaria vectors per host depends on his or her biomass35 (Figure 1
), this overestimates the challenge faced by small children. Second, since the probability that an inoculation is effective (the success probability) decreases as the EIR increases, the infection rate is not proportional to the EIR6 (Figure 2
). There is a similar relationship between P. falciparum prevalence and EIR when different study sites are compared7 and relationships similar to those in Figure 2
have also been found in a study of infection incidence in infants in Tanzania8 who have little prior exposure to malaria. This attenuation phenomenon therefore probably represents an effect of innate, not acquired, immunity.
|
|
| METHODS |
|
|
|---|
Let wi be the proportion of the total weight of people sleeping under the net contributed by person i. The expected proportion of mosquitoes feeding on that person, f(wi), (equivalent to the feeding likelihood9) is then a function of wi satisfying the constraints wi = 0
f(wi) = 0 and wi < wj
f(wi) < f(wj). These constraints are implied by the fact that there is no threshold weight, below which the host will not be bitten, however both the weight and biting rate must be positive, and the estimated biting rate should increase with the biomass of the host (Figure 1
). Moreover, two hosts of equivalent weight sleeping together are equally likely to be bitten so that wi = 0.5
f(wi) = 0.5.
To obtain estimates of f(wi) as a function of wi we fitted a Bayesian monotone regression model to the original published data. This allowed the fitted relationship to be as flexible as possible within the constraints listed above (Appendix 1).
The ratio between the number of bites received by individual i and its companion(s) is f(wi)/[1 f (wi)]. This quantity corresponds to the expected ratio between the inoculation rates of any individuals whose relative weights are wi : 1 wi. We used this relationship to convert the values of f(wi) to relative EIR values for different ages by computing the mean growth in weight for a typical African rural population (from Tanzania).10 Defining Wa as the mean weight in age group a, we estimated Ea, the adjusted EIR, for age-group a from
![]() |
where the age group 20<25 years, a typical age for mosquito collectors, is used as a standard for both weight and EIR.
Dependence of force of infection on EIR.
Beier and others6 studied the reinfection of 21 cohorts of between 43 and 50 children six months to less than six years old, each followed up during four successive two-week periods, following 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 14-day period.
We used the program WinBUGS version 1.3 (Medical Research Council Biostatistics Unit, Cambridge, United Kingdom) to fit a parametric relationship between the probability of infection, in children aged 0.56 years in a given day, h0.56, and EIR of the form
![]() |
where g and ß are both parameters describing the host susceptibility. We allow for imprecision in the EIR assessment by assuming Poisson variation in the counts of sporozoite positive mosquitoes.
To predict ha, the incidence of infection in age group a, from a given measured EIR (E2025), we use Equation 1 to estimate Ea and E0.56 and then substituted the relative number of inoculations into Equation 2, i.e.,
![]() |
Studies of outcomes in relation to incidence of infection.
To examine effects of age of the host and malaria transmission intensity on the incidence of clinical malaria, we reanalyzed data from Ndiop (EIR = 20) and Dielmo (EIR = 200) in Senegal1 (Figure 3
). We made estimates of the age-specific incidence of infection, ha, for each site and month of the year separately, using the published seasonal patterns of transmission1,11 as the source for monthly EIR figures. We then compared the age-incidence patterns reported in the original publication with estimates of the ratio of incidence of clinical malaria to ha, i.e., of the number of clinical episodes for each new infection event.
|
|
| RESULTS |
|
|
|---|
The relationship that we estimated between the infection incidence in children, h0.56 and the EIR values for the same cohorts in Kenya (E2025) shows a good fit to the data (Figure 2
) with estimates of g = 0.073 (SE = 0.011) and ß = 0.097 (SE = 0.025) (with time units of one day). This curve is similar to that reported in a study of infants in Tanzania,8 but with rather higher estimates of incidence of infection.
Analysis of clinical malaria incidence in Senegal.
The raw age-specific incidence data from Dielmo (Figure 3a
) show the typical age-incidence pattern for malaria in very high transmission areas with the peak of malaria incidence in very young children, and with very few clinical episodes attributable to malaria in older children and adults. Ndiop shows the typical pattern for a lower-transmission area, with the peak of incidence in older age groups.
The values of ha estimated for Dielmo and Ndiop from Equation 4 illustrate the age-dependence that we estimate (Figure 3b
). The steep increase with age in ha because of increasing body mass results in very different age profiles of the estimated number of episodes per new infection (Figure 3c
), from those in Figure 3a
. In both villages, the peak incidence of disease relative to that of infection is in younger children, one year of age in the case of Dielmo, and less than one year of age in Ndiop. However the number of episodes of disease in those less than one year of age in Ndiop is even higher than the infection rate.
The lower success rate of inoculations at high EIR than at low EIR (attenuation) also has substantial effects. In adults, the estimated infection rate in Ndiop is about one-fourth of that for Dielmo, despite the 10-fold difference in EIR.
Analysis of mortality data.
The infant mortality rate (IMR) shows a strong and statistically significant positive relationship with EIR2 (Figure 4a
). The reported EIR values span more than three orders of magnitude (thus, we plotted them on a logarithmic scale) and the curve fitted in the original publication2 parameterizes the IMR as a function of the logarithm of the EIR, i.e.:
![]() |
where IMR1 is equal to the IMR at an EIR of one infectious bite per year (approximately equal to the lowest level of endemic transmission), which is estimated to average 62 per 1,000 for the given study sites. However, there is considerable scatter between the points, making it very uncertain what shape of curve, if any, should be fitted.
The estimated duration of the transmission season for the 20 localities studied ranged from 4 to 12 months (mean = 6.8 months, SD = 2.0). At low EIR values, Equation 3 is equivalent to
![]() |
i.e., h is roughly proportional to the EIR. The seasonal distribution of the inoculations then makes little difference to h0 because the same constant of proportionality, gß E2025/E0.56, applies to the data for each month. However, when there are a large number of inoculations in a single month, h0 for that month tends towards the upper limit, g, implied in Equation 4. This has the consequence that at high EIR levels, the extent of attenuation depends crucially on the duration of the transmission season, with the success probability of an inoculation higher when the bites are distributed more evenly over the year. In the particular dataset available, several of the high EIR study sites have a transmission season that is short but intense, so that the estimated value of h0 for each month of transmission is close to g. This has the effect of bringing the annualized estimates of the incidence of new infections (ha) closer together at high EIR than they would be if they fell on a smooth curve drawn through the average value of ha at each EIR (Figure 4c
), and the ranking of the values of ha for these points differs from the ranking of the EIR.
This means that, if anything, it is even more difficult to discern any pattern in the IMR/h0 relationship at high values of h0 (Figure 4b
) than it is to discern pattern in the IMR/EIR relationship (Figure 4a
). The mortality rates for children 1259 months of age in the same studies were all much lower than mortality rates in infancy and there were no obvious relationships between these mortality rates and either EIR2 or ha.
| DISCUSSION |
|
|
|---|
Our attempt to make age-specific estimates of h from published EIR data is of necessity crude. It is not obvious that we have chosen the best measure of host body size, nor that the same calibration of biting rates by body size should be used for all transmission settings and sleeping patterns. Similarly, we make a strong assumption by applying Equation 3 for all age groups with the same parameters g and ß.
The very high ratio of clinical attacks to estimated new infections in the youngest age group in Ndiop also suggests that our model does not account for all the factors leading to age-dependence in infection incidence. Since all these clinical attacks were treated, repeated episodes of disease caused by the same infections cannot explain the high clinical attack rate, so we infer that the success rate of inoculations in infants is actually higher than our algorithm suggests. One possible explanation for this would be that pre-erythrocytic immunity (which is not allowed for in our calculations) also plays a role in determining the age-dependence of h. This would predict that the h increases less with age than we have estimated.
Our estimates of ha cannot therefore be considered definitive, but they nevertheless contribute to our understanding of epidemiologic patterns by clarifying which relationships can, and which cannot, be accounted for by effects of host biomass and of saturation of ha at high EIR.
There is strong evidence both from molecular typing13 and from patterns of seasonality in morbidity14 that clinical malaria is normally caused by newly invading parasites, and the most severe symptoms generally accompanying the first peak of parasite density after infection. It follows that in the short-term, any reduction in EIR will decrease the incidence of clinical episodes in proportion to the effect on the force of infection.
Reductions in transmission intensity, however, also reduce immunologic stimulation, and this may have longer-term effects, in particular resulting in shifts of the peak in the age-incidence profiles to older ages, as seen in the comparison of Dielmo and Ndiop. Consequently, it appears that a reduction in EIR will have little effect on the total number of episodes in a lifetime. However, without a reliable model for clinical immunity, the long-term effects of any intervention cannot be predicted with confidence. One model is that P. falciparum consists of a number of distinct strains, and that clinical episodes result only on the first occasion that a strain infects a host.15 Within this model, each infection confers long-term clinical immunity so that subsequent infections with the same strain do not result in clinical attacks. However, analysis of exposure-incidence relationships in Tanzanian infants suggested that the main determinant of clinical immunity is the duration of exposure to blood-stage parasites, and not the number of inoculations per se.16 There is also evidence for a role of short-term effects that depend on cross-protection between co-infecting parasites.17,18
None of these models alone can explain why the peak incidence is not in the youngest children. An explanation for this is provided by Figure 3c
, which indicates that small size is likely to be the main explanation for why incidence of clinical episodes is lower in infants. The surprisingly low incidence in infants in Dielmo could be due to passive immunity during the first 36 months of life, which means that infections then only rarely result in clinical episodes.19 In contrast, our analysis suggests that in the youngest children in Ndiop effectively all new infections result in clinical malaria attacks. Indeed, when considered relative to the infection rate, the surprise is the very high incidence in infants in Ndiop and low incidence in Dielmo. Part of the explanation could be due to the lack of allowance for acquired immunity against pre-erythrocytic stages of the parasite (see above), meaning that incidence of infection in infants in Ndiop may actually be higher than we have estimated.
While it is of some concern that reducing malaria transmission in endemic areas does not necessarily substantially reduce the incidence of clinical episodes, the public health impact of malaria depends mainly on its effects on mortality. The impact on mortality appears to show very different relationships with transmission intensity from those on acute morbidity.
The relationships of malaria-specific child mortality rates with EIR are unclear.20,21 Part of the explanation for this is that verbal autopsies for malaria are known to perform poorly,2225 leading to very imprecise estimates of malaria-specific mortality. Because malaria can be effectively treated the case-fatality rate also varies enormously depending on the effectiveness of the health services. It can be especially high in epidemics, when a high incidence rate during a short period results in social disruption.
The analysis of malaria-specific mortality is further complicated by the abundant evidence that malaria causes many deaths indirectly that are classified as due to other causes, but that would have been avoided in the absence of malaria.26,27 Much of this indirect morbidity and mortality is accounted for by acute respiratory infections, and co-infection with other pathogens worsens the outcome in malaria,28 so many deaths are clearly multi-factorial. Clear evidence for this came from the decreases in all-cause mortality in malaria control trials in the 1950s to the 1970s, which were greater than could be explained from the malaria rates.27,29 A particularly telling example is the convergence of death rates in malarious districts of Sri Lanka with those in districts without malaria, when effective control was introduced. The decreases in mortality were much greater than the mortality rates that had previously been attributed to malaria.27 A second good example is the disappearance of epidemics of acute respiratory infections from sugar plantations in coastal Guyana when the malaria was eradicated.30 Prior to eradication, these epidemics were synchronized with those of malaria.
Ultimately, it matters little whether the primary cause of a death is malaria, so for most purposes it is more relevant to analyze all-cause mortality than to concentrate only on diagnosed malaria. Randomized controlled trials of insecticide-treated bed nets and curtains have demonstrated efficacy of 1442% in reducing all-cause mortality in children one month to five years old3136 across a range of transmission intensities. This is overwhelming evidence for a benefit of transmission reduction. The strong increase in IMR with EIR (Figure 4a
) reinforces this conclusion and confirms the substantial public health benefit in reducing malaria transmission. However, the number of data points relating IMR to EIR is limited, and there is uncertainty about what shape of curve, if any, should be fitted through them to predict the effects of any specific reduction in transmission.37 This motivated us to compute h0, since we hoped to test whether the saturation effect suggested by Figure 4a
can be accounted for by the saturation in the relationship between h0 and the EIR (Figure 2
). Since h0 itself tends to a limit as the EIR increases, it would be surprising if mortality, a secondary consequence of infection, failed to plateau at high EIR, but it is less obvious what relationship might be anticipated between IMR and h0 (Figure 4c
). In fact, by demonstrating that the available points covered only a limited range of values of h0, our analysis served less to clarify this relationship than to emphasize the inadequacy of the available data for estimating it. There is clearly a need for further investigations of these relationships, and measurements of malaria transmission intensity currently in progress at a number of demographic surveillance sites across Africa will have an important role in extending the dataset.
For children 1259 months old, the limited data currently available do not suggest any particular relationship between the mortality rates and either EIR2 or ha. In the trials of insecticide-treated bed nets, the benefits of transmission reduction were also clearest in the youngest children, and at high levels of transmission (in Ghana) the benefit was limited to children less than two years of age at the start of the trial.38 Unlike the effects of transmission reduction on clinical episodes, long-term follow-up of two of these trials after the control group was given the intervention indicated that the increased survival in very young children has no price in higher mortality later on38 (Diallo DA, unpublished data). This is in keeping with findings from trials of indoor residual spraying.29,39
The finding that the main benefit of reducing transmission on mortality is in the youngest children is also at variance with effects on clinical episodes, since the latter are infrequent in the first few months of life. It is likely that the effects on infants largely relate to indirect mortality. Life-threatening acute respiratory infections and diarrhea and thus co-infection are most frequent in infants, and the mortality risk from severe anemia is highest in the second half of the first year of life. Sick infants have the greatest difficulty communicating their distress to care givers. Perhaps most importantly, low birth weight caused by malaria in pregnancy probably leads to a very substantial component of infant mortality.40 In endemic areas, malaria prophylaxis during pregnancy is estimated to correspond to lead to an average increase of 0.101 kg in the birth weight of first-born children, corresponding to a reduction of 11 per 1,000 in their neonatal mortality rate.41 Unfortunately no studies have examined the relationship between either birth weight or placental malaria and the EIR, so it is not possible to dissect out how these specific risk factors are affected by incomplete malaria control.
These indirect effects may well be most important in determining the burden of malaria in Africa and it is here that very substantial benefits are likely to be seen when levels of transmission are reduced.
| APPENDIX 1 |
|
|
|---|
To ensure that f(wi) was an increasing function of wi, we calculated it from f(wi) = f(wi1) +
i where f(w0) = 0 and the
i are estimated using a Markov Chain Monte Carlo algorithm using the program WinBUGS 1.3.42 The
i values were constrained to be positive by assigning Uniform (0, 1) prior distributions to a set of parameters
'i and computing
i from
. This transformation also constrained f(wi) to be <0.5 for the people included in the analysis.
Received August 21, 2003. Accepted for publication January 9, 2004.
Acknowledgments: We thank Don de Savigny, Louis Molineaux, Klaus Dietz, Penelope Vounatsou, and Amanda Ross for their helpful comments on an earlier draft of this manuscript.
Authors address: T. Smith, G. Killeen, C. Lengeler, and M. Tanner, Swiss Tropical Institute, Socinstrasse 57, PO Box CH-4002, Basel, Switzerland, Telephone: 41-61-284-8273, Fax: 41-61-271-7951, E-mail: Thomas-A.Smith{at}unibas.ch.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
W. P. O'Meara, T. W. Mwangi, T. N. Williams, F. E. McKenzie, R. W. Snow, and K. Marsh Relationship Between Exposure, Clinical Malaria, and Age in an Area of Changing Transmission Intensity Am J Trop Med Hyg, August 1, 2008; 79(2): 185 - 191. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Dicko, I. Sagara, D. Diemert, M. Sogoba, M. B. Niambele, A. Dao, G. Dolo, D. Yalcouye, D. A. Diallo, A. Saul, et al. Year-to-Year Variation in the Age-Specific Incidence of Clinical Malaria in Two Potential Vaccine Testing Sites in Mali With Different Levels of Malaria Transmission Intensity Am J Trop Med Hyg, December 1, 2007; 77(6): 1028 - 1033. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. J. Dunachie, M. Walther, J. E. Epstein, S. Keating, T. Berthoud, L. Andrews, R. F. Andersen, P. Bejon, N. Goonetilleke, I. Poulton, et al. A DNA Prime-Modified Vaccinia Virus Ankara Boost Vaccine Encoding Thrombospondin-Related Adhesion Protein but Not Circumsporozoite Protein Partially Protects Healthy Malaria-Naive Adults against Plasmodium falciparum Sporozoite Challenge. Infect. Immun., October 1, 2006; 74(10): 5933 - 5942. [Abstract] [Full Text] [PDF] |
||||
![]() |
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] |
||||
![]() |
T. SMITH, N. MAIRE, K. DIETZ, G. F. KILLEEN, P. VOUNATSOU, L. MOLINEAUX, and M. TANNER RELATIONSHIP BETWEEN THE ENTOMOLOGIC INOCULATION RATE AND THE FORCE OF INFECTION FOR PLASMODIUM FALCIPARUM MALARIA. Am J Trop Med Hyg, August 1, 2006; 75(2_suppl): 11 - 18. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. SMITH, A. ROSS, N. MAIRE, C. ROGIER, J.-F. TRAPE, and L. MOLINEAUX AN EPIDEMIOLOGIC MODEL OF THE INCIDENCE OF ACUTE ILLNESS IN PLASMODIUM FALCIPARUM MALARIA. Am J Trop Med Hyg, August 1, 2006; 75(2_suppl): 56 - 62. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Lundquist, L. K. Nielsen, A. Jafarshad, D. SoeSoe, L. H. Christensen, P. Druilhe, and M. H. Dziegiel Human Recombinant Antibodies against Plasmodium falciparum Merozoite Surface Protein 3 Cloned from Peripheral Blood Leukocytes of Individuals with Immunity to Malaria Demonstrate Antiparasitic Properties. Infect. Immun., June 1, 2006; 74(6): 3222 - 3231. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. G. BREMAN, M. S. ALILIO, and A. MILLS CONQUERING THE INTOLERABLE BURDEN OF MALARIA: WHAT'S NEW, WHAT'S NEEDED: A SUMMARY Am J Trop Med Hyg, August 1, 2004; 71(2_suppl): 1 - 15. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |