Skip to main content

Distribution of malaria exposure in endemic countries in Africa considering country levels of effective treatment

Abstract

Background

Malaria prevalence, clinical incidence, treatment, and transmission rates are dynamically interrelated. Prevalence is often considered a measure of malaria transmission, but treatment of clinical malaria reduces prevalence, and consequently also infectiousness to the mosquito vector and onward transmission. The impact of the frequency of treatment on prevalence in a population is generally not considered. This can lead to potential underestimation of malaria exposure in settings with good health systems. Furthermore, these dynamical relationships between prevalence, treatment, and transmission have not generally been taken into account in estimates of burden.

Methods

Using prevalence as an input, estimates of disease incidence and transmission [as the distribution of the entomological inoculation rate (EIR)] for Plasmodium falciparum have now been made for 43 countries in Africa using both empirical relationships (that do not allow for treatment) and OpenMalaria dynamic micro-simulation models (that explicitly include the effects of treatment). For each estimate, prevalence inputs were taken from geo-statistical models fitted for the year 2010 by the Malaria Atlas Project to all available observed prevalence data. National level estimates of the effectiveness of case management in treating clinical attacks were used as inputs to the estimation of both EIR and disease incidence by the dynamic models.

Results and conclusions

When coverage of effective treatment is taken into account, higher country level estimates of average EIR and thus higher disease burden, are obtained for a given prevalence level, especially where access to treatment is high, and prevalence relatively low. These methods provide a unified framework for comparison of both the immediate and longer-term impacts of case management and of preventive interventions.

Background

The prevalence of Plasmodium falciparum infections is routinely measured in malaria indicator surveys (MIS), and as part of various health assessments and research projects. Prevalence data are therefore relatively widely available and are often used as a measure of endemicity in geographical comparisons and in evaluating the success of intervention programmes [1]. However, although prevalence is a consequence of malaria transmission and levels of exposure, these variables do not have a one-to-one relationship but rather a non-linear relationship modified by many factors such as naturally acquired immunity, malaria interventions and of heterogeneity in transmission rates [2]. These complicate the interpretation of age-patterns of infection and disease. The relationship between exposure and prevalence of infection also depends on the amount of treatment in the population because treatment truncates infections and (depending on the drug regimen) provides a few weeks of chemoprophylaxis (Fig. 1). If access to effective treatment is good, then prevalence may remain relatively low, even at high transmission levels. The amount of effective treatment also affects the relationships of exposure (or prevalence) with morbidity, and mortality rates (Fig. 1).

Fig. 1
figure 1

Illustration of impact of treatment effectiveness on interactions between infection, clinical disease and exposure. Arrows indicate causal links and double lines show where treatment has a modifying effect

Human exposure to malaria, one part of malaria transmission, is best quantified by the entomological inoculation rate (EIR: the number of infectious bites per human host, per unit time), which is more directly related to morbidity and mortality than is prevalence. However, measuring this quantity directly requires intensive entomological studies over the whole annual period of malaria transmission. Previously established empirical relationships between prevalence and EIR have illustrated the complications and diversity by site [3]. EIR data are consequently relatively sparse, and indirect methods, that ideally account for treatment effects, are needed for estimating EIR from available prevalence data [4, 5].

The comprehensive repository of geo-located malaria parasite prevalence data maintained by the Malaria Atlas Project (MAP) is the obvious starting point for estimating how many people are exposed to malaria at different intensities, in different endemic countries. Several different algorithms have been used to infer the distribution of exposure from prevalence maps. In particular, a linear relationship between prevalence and the logarithm of the EIR approximates the empirical relationship between these variables [6], and the MAP repository includes EIR surfaces and estimates of the uncertainty based on this relationship [7]. Other researchers use process models to estimate transmission rates surfaces from prevalence data [6, 8ā€“10]. These analyses do not allow for effects of treatment on prevalence. At low transmission levels, where infection events are sporadic, and superinfection relatively infrequent, this omission can be remedied using rather simple models for translating prevalence into transmission estimates, conditional on the incidence of effective treatment [4]. At higher levels of transmission, both concurrent and sequential superinfection are frequent; so mechanistic models allowing for this, as well as for treatment rates, are needed.

Estimates of the number of clinical malaria episodes at national level and continent-wide have been made from the MAP database by assuming a standard empirically determined relationship between prevalence and the incidence of clinical malaria in children [11, 12]. Using a similar methodology based on geographical stratification of risk, estimates of clinical incidence at national level are made yearly by the World Health Organization (WHO) for the World Malaria Report (WMR) for high-burden sub-Saharan countries [13]. This report also provides up-to-date assessments of malaria-related interventions and policies, attempting to quantify the impact on disease burden. Estimates of clinical incidence for each year have been made by adjusting for changing intervention coverage levels within each country, assuming effects match those seen in controlled trials [14].

These estimates of clinical incidence do not allow for levels of access to effective treatment. This affects both the true extent of pathology, and the observed clinical incidence, whether ascertained passively or actively. Depending on underlying exposure, high treatment levels create a virtuous cycle by averting further pathology and secondary cases. Estimates of worldwide and national levels of burden should, therefore, take into account effects of treatment, as well as the shifts in age patterns of prevalence [15] and of incidence that occur as a result of transmission reducing interventions.

The OpenMalaria platform supports an ensemble of models that can be used for calibrating different malariological indices against each other [16]. OpenMalaria is a stochastic, individual-based, simulation model of malaria in humans [17] linked to a deterministic model of malaria in mosquitoes [18]. The simulation model includes sub-models of infection of humans [19], blood-stage parasite densities [20], infectiousness to mosquitoes [21], incidence of morbidity including severe and hospitalisation [22, 23] and mortality [22]. An ensemble of 14 model variants is available [24] with each model including different assumptions for decay of natural immunity, greater within-host variability between infection and entomological exposure, heterogeneity in transmission and heterogeneity in susceptibility to co-morbidities.

Six of the OpenMalaria ensemble models were used in this work to compute estimates of the distribution of exposure (EIR) for each of 43 malaria endemic countries in sub-Saharan Africa as well as estimates of clinical incidence (and also incidence of severe disease and malaria mortality) for 2010 levels of malaria control. These estimates are based on the pixel-level posterior distributions of parasite prevalence in 2010 published by MAP [7]. For each country, these estimates are conditional on national level estimates of the levels of access to effective treatment for malaria fevers [25]. The resulting estimates of the distribution of transmission and of the incidence of clinical malaria provide a basis for evaluating the impacts of both preventive and curative intervention programmes allowing for the effects of existing case management on prevalence and burden of disease.

Methods

An overview of the methods in estimating malaria exposure distributions (as EIR) and resulting burden is presented in Table 1, including inputs and outputs of each method.

Table 1 Description of the EIR and burden estimation methods A and B including their inputs and outputs

Malaria prevalence data

National levels of prevalence were taken from the prevalence surfaces estimated by the Malaria Atlas Project (MAP) for Plasmodium falciparium 2010 [7]. Estimates of prevalence in children aged 2 up to children before their 10th birthday (PfPR2ā€“10) across a 5 km Ć— 5 km grid were extracted as posterior distributions from a Markov Chain Monte Carlo (MCMC) calculated via a Bayesian geostatistical model using survey data. The primary estimates are of PfPR2ā€“10 are available from the MAP website [26] as posterior densities.

National levels of effective treatment coverage

National levels of access to effective malaria treatment were collated previously [25] and are detailed in Table 2. Effective malaria treatment is treatment that results in parasitological cure. In this work effective treatment are estimates of the probability, E 14, that effective malaria treatment will be obtained during any 14-day period in which a fever occurs. Estimates were assembled at country level taking into account multiple factors for effectiveness of malaria case management, including probably of treatment-seeking, of type of care provider, of systems compliance with the recommended anti-malarial treatment, of adherence with the drug regimen, and the quality of the anti-malarial medications.

Table 2 Coverages of effective treatment and estimated transmission profiles (EIR mean, median and quartiles) for 43 sub-Saharan Africa estimated by method B assuming country level effective treatment

Relationship between parasite prevalence and EIR

Two different methods were used to estimate distributions of malaria exposure (as measured by the entomological inoculation rate, EIR) from PfPR2ā€“10 data:

Method A: statistical relationship between EIR and prevalence

A previously published statistical model transforming PfPR2ā€“10 to EIR [7] (Additional file 7 of that paper), based on an earlier empirical analysis of the relationship of measured EIR values with PfPR2ā€“10 [6]:

$$x \sim \log {\text{Normal}}\left( {\mu ,\sigma^{2} } \right)$$
(1)

where x is EIR, Ī¼ = 1.768 + 7.247p, Ļƒ = 1.281, and p, is PfPR2ā€“10. This relationship is independent of the level of access to effective treatment, E 14, and thus does not allow for the effects of case management on the prevalence-EIR relationship. This model allows for statistical uncertainty in both variables x and p (data and fitted curve shown in Fig. 2a). Scale factors can be used to obtain the EIR estimate that would be obtained with different measurement approaches (e.g. pyrethroid-spraying catches, human landing catches, or both). This method is similar to the method used in previous analyses of the global burden of clinical malaria [11].

Fig. 2
figure 2

Relationships between malaria exposure (EIR), effective coverage, and prevalence for Method A (a) and Method B (b). a Method A: plotted empirical relationship of prevalence as a function of EIR relationship [6] (Eq. 1) with data used to fit this relationship. The relationship between standardized prevalence and EIR is approximately linear-log for all the data (grey curve fitted relationship over all data). The relationship varies by study (purple colour dots a correspond to data from a single field study, purple curve the fitted relationship to those data) and by method (red dots correspond to measurements taken via pyrethroid spray catches, and purple and blue dots measurements taken by other methods. Red and blue curves correspond to the respective fits). b Method B: OpenMalaria simulations of the relationship between prevalence and EIR (model variant R0133 only, other models shown in Additional file 2: Figure S3) for discrete levels of coverage of effective treatment (points) and best fitted model to these data as Hill functions (curve for different levels of effective treatment) (Eq. 3). Colour indicates the level of effective treatment (E 14), with red 0.001 %, yellow 5 %, light blue 20 %, dark blue 40 %

Method B: dynamic model relating EIR, prevalence, and coverage of treatment

Method B uses relationships between EIR and prevalence derived from multiple transmission models of malaria epidemiology and control, incorporating the effects of treatment on the infectious reservoir. The process of translating prevalence to EIR is illustrated in Fig. 3a, b, in essence extracting prevalence distributions at each 5 by 5 km grid from MAP (detailed above, Fig. 3a) and converting to EIR by the fitted relationship from OpenMalaria for a given coverage of effective treatment (Fig. 3b).

Fig. 3
figure 3

Schematic diagrams of the processes in estimating EIR distributions and disease burden from MAP prevalence. The figure illustrates the steps involved in estimating geographic specific EIR (Method B) and incidence levels of malaria (Methods A and B), which includes the dynamic effect of treatment on transmission, and a dynamic model of clinical incidence. a, b Illustrate the process of extracting prevalence distributions from MAP [7] by pixel (5 km by 5 km) and converting to distributions of EIR using a statistical relationship relating prevalence and EIR for given levels of effective treatment derived from OpenMalaria simulations. Method A, not illustrated, is simpler in that it does not consider the effect of treatment on transmission, and uses the WMR method for estimating EIR. c Illustrates aggregation of these EIR distributions from pixel level to a larger spatial area, such as country level. d Illustrates the process of estimating country level burden for a distribution of EIR (derived from either method A or B), namely the EIR distributions are inputs to micro-simulation OpenMalaria with outputs of incidence of clinical cases (and mortality) which are calculated for a given coverage of effective treatment E14. The Gaussian distributions in aā€“d are illustrative only

The transmission models are six model variants from the OpenMalaria stochastic individual-based model of the dynamics of P. falciparum malaria in humans [24] (Table 3), comprising a subset a previously published model ensemble [24] with each model variant including the same sub-model for pathogenesis [23] and case-management [27], but differing by assumptions concerning immunity decay or heterogeneity in transmission or co-morbidity (Table 3). The same parameterizations as used previously [24] were used to capture human demography and the seasonality of transmission. Each model variant has been parameterized by fitting to observed relationships between seasonal patterns of EIR and a range of outcomes, including parasite prevalence [19] and morbidity rates [23] in specific field sites.

Table 3 Model-specific parameters for each model variant for the statistical models fits relating OpenMalaria EIR and prevalence among 2ā€“10 year olds

A statistical relationship was fit between simulated PfPR2ā€“10, p, and EIR, x, for a given level of effective treatment, E 14, for each model in the ensemble (Fig. 2b illustrates an example of this relationship). These simulated predictions cover a wider range of EIR and prevalence used to parameterize the transmission models originally [17, 24]. The OpenMalaria simulations use a 5-day time step and effective treatment at each 5-day time step, E 5, was obtained from the 14 day estimates using a mapping based on the pattern of fevers over time in malaria-therapy data [28] (sample values shown in Additional file 1: Table S1). A Hill function was fitted by least-squares to the simulation data in order to relate PfPR2ā€“10 and EIR, namely:

$$p\left( {x,E_{14} } \right) = \frac{{p_{max} x^{{n\left( {E_{14} } \right)}} }}{{K^{{n\left( {E_{14} } \right)}} + x^{{n\left( {E_{14} } \right)}} }},$$
(2)

where, p max , K, and n are functions of E 14. The inverse of relationship Eq. (2), relating EIR to PfPR2ā€“10 is given by

$$x = K\left( {E_{14} } \right) exp\left[ {\frac{1}{{n\left( {E_{14} } \right)}}ln\left( {\frac{{p\left( {x,E_{14} } \right)}}{{p_ {max}\left( {E_{14} } \right) - p\left( {x,E_{14} } \right)}}} \right)} \right].$$
(3)

The functional forms for x, p max and K were chosen among exponential, linear, and quadratic options to give the best fit of p(xE 14) to the simulated prevalence for different levels of coverage of effective treatment. The selected functions are:

$$K\left( {E_{14} } \right) = K_{1} { exp }\left( {K_{2} E_{14} } \right)$$
(4)
$$n\left( {E_{14} } \right) = n_{1} E_{14}^{2} + n_{2} E_{14} + n_{3}$$
(5)

and

$$p_{max} \left( {E_{14} } \right) = p_{1} exp\left( {p_{2} E_{14} } \right) .$$
(6)

where K 1, K 2, n 1,n 2, n 3, p 1, and p 2 are fitted parameters. Separate parameter sets were fitted for each of the six model variants in the ensemble (values provided in Table 3).

Estimation of EIR distributions at national level

The prevalence-EIR relationships from method A and B were used to estimate a distribution of EIR for each country from the prevalence surfaces estimated by the MAP for 2010 [7]. Prevalence from the MCMC chains are weighted by each pixel-level value of population, and the percentiles of the distributions obtained by summarizing the whole set of MCMC chains. Corresponding to the PfPR2ā€“10 value for pixel j, from MAP, and MCMC iteration i, an EIR estimate, x (i) j , is obtained. For method B this is

$$\mathbf{ x_{j}^{(i)}} = K\left( {E_{14} } \right) exp \left[ {\frac{1}{{n\left( {E_{14} } \right)}}\ln \left( {\frac{{p_{j}^{(i)} }}{{p_{{max} } - p_{j}^{(i)} }}} \right)} \right].$$
(7)

The corresponding estimate of the distribution of EIR over the whole country (including non-endemic areas, with EIR = 0) is obtained by binning x (i) j into a limited number, K, of ranges \(X_{ 1} ,X_{ 2} , \ldots ,X_{K}\), and population weighting. Aggregating the estimates from the whole set of T sampled values of p (i) j from the MCMC chains, each range is assigned probability:

$${ \Pr }\left( {X_{k} } \right) = \frac{{\mathop \sum \nolimits_{i} \mathop \sum \nolimits_{j} \left( {N_{j} {\text{I}}\left({ \mathbf{{{x}}_{{j}}^{{\left( {{i}} \right)}}} \in X_{k} } \right)} \right)}}{{T\left( {\mathop \sum \nolimits_{j} N_{j} } \right)}},$$
(8)

and I(x (i) j āˆŠ X k ) is an indicator taking value 1 if x (i) j is in range X k and zero otherwise. Here N j is the population assigned to the pixel as determined by the gridded population of the world [29, 30]. For computational convenience we carried out the summation over i before summing over j.

The resulting distributions describe the proportion of each countryā€™s population that one would expect to be living at a given level of prevalence. In many of the countries analysed, a proportion of the gridded population from [30] falls outside the boundary of the area defined by MAP as being within the spatial limits of endemic malaria transmission [7]. This proportion of the population was assigned an EIR of zero.

Two different estimates of the transmission distribution per geographic area were calculated by estimation method B, to examine sensitivity to the estimated level of access to effective treatment. To capture the situation before recent scale-up of artemisinin combination treatment (ACT), a common value of access to effective treatment for all countries was used assumed and at a value previously used in OpenMalaria simulations [31]. This value equates to approximately 15 % of all malaria cases receiving treatment resulting in parasitological cure. In addition, analyses were conducted using country-specific estimates of access to effective treatment [25]. Country levels of coverage of effective treatment are listed Table 2 and illustrated by map in Additional file 2: Figure S1 and in Fig. 3b.

Burden of disease

National level estimates of the incidence of clinical malaria were projected from the EIR distributions derived from Method A and Method B using OpenMalaria simulations. These incorporate dynamic models of clinical incidence and treatment parameterized with Senegalese and Tanzanian data [23, 31] and models for severe disease and mortality [22], and hence provide clinical incidence estimates as an extension of EIR estimation (process illustrated in Fig. 3d). Separate estimates were made using the EIR estimates with Method A, those from with Method B with E14 = 0.15, and those made with Method B with country specific E14 values.

Estimated burden, via clinical incidence, derived by both methods was compared with those national level estimates of clinical malaria from the WMR. For most sub-Saharan African countries, these use a standard empirical relationship between clinical incidence and endemicity. Clinical incidence values were assigned to each endemicity level based on estimates of the numbers of events recorded in longitudinal surveys of febrile malaria episodes in children, detected either actively or passively [32ā€“34], established independently of effects of treatment rates [12]. For countries with low endemicity, WMR uses national surveillance data to estimate burden, with adjustments to allow for incomplete reporting.

Results

National level prevalence distributions

National levels of PfPR2ā€“10 aggregated at country level after extracting from MAP [7] posterior distributions at each 5 by 5 km pixel illustrate high average levels of PfPR2ā€“10 in 2010 in many African countries but for many countries also the wide distribution of prevalence levels (values summarized in Table 4 and Additional file 1: Table S3 and shown as distributions in Additional file 2: Figure S2). Regional differences, local variation, and uncertainty within areas all contribute differently to the overall distributions, with the average levels of transmission highest in West and Central Africa. Much of Namibia, Botswana, and South Africa, and also several Sahelian countries are malaria free, as are highland areas of East Africa. Some of the variation is also a result of differences in the extent of recent intervention programmes. In some countries, intervention programmes have had little impact on 2010 prevalence (e.g. Benin, Burkina Faso, CĆ“te dā€™Ivoire), while elsewhere prevalence has been considerably reduced in the last decade (e.g. Senegal, Tanzania, Zambia), or much of the population lives in areas on the margins of stable transmission (Somalia, North Sudan). The location of some major urban centres such as Nairobi and Lusaka at relatively high altitudes, with low transmission, strongly influences some of these profiles.

Table 4 Prevalence distributions, summarized for each country: estimated prevalence (mean, median and quartiles) for 43 sub-Saharan Africa estimated from MAP prevalence posteriors at 5 km by 5 km grids, aggregated to country level and weighted by population. Mean EIR estimates for Method A aand Method B assuming effective coverage (E 14) of 15ā€‰% for all countries

Modelled relationships between EIR and prevalence

Where PfPR2ā€“10 is high, the OpenMalaria models predict on average a slightly higher EIR at a given prevalence than does the empirical model (method A), with relatively little influence of effective treatment (E 14) (Fig. 2a, b). The fitted prevalence to EIR relationships for Method B are shown in Fig. 2b (model variant R133) and Additional file 2: Figure S3 (all 6 model variants). The six model variants all predict broadly similar, but nevertheless distinct, prevalence-EIR relationships. The general pattern for Method A is for prevalence to increase steeply with EIR at low transmission levels, but to saturate at higher transmission (Fig. 2a). The considerable variation around the best fitting curve for Method A, after adjusting for the different EIR measurement techniques used, is treated as random variation that contributes to uncertainty in the estimate of EIR from prevalence. This analysis does not allow for variations in the coverage or effectiveness of case management in the different studies, however such variation could account for much of this unexplained dispersion (compare with Fig. 2b). At lower transmission levels the fitted curves for Method B (OpenMalaria) vary considerably with E 14, suggesting that effectiveness of case management is a particularly important driver of prevalence in such settings, with Method B estimating lower EIR at a given prevalence than the empirical model unless E 14 is high (Fig. 2b). This is partly because Method B constrains estimated EIR to be zero at zero prevalence, while the empirically-based Method A does not capture or force this constraint.

National level EIR distributions

The differences between the two relationships for prevalence and EIR are reflected in the estimated distributions of EIR by county (Fig. 4; Additional file 2: Figure S4). The EIR distributions are generally much more highly skewed than are the prevalence distributions. The distributions obtained with the empirical model (Method A: Fig. 4; Table 4 and Additional file 1: Table S4) and with the simulation models (Method B) that assume \(E_{14} = 0.15\) (Table 4 and Additional file 1: Table S5) are similar to each other for most countries, though the estimated median EIRs are generally somewhat higher for the estimates from the simulation models (Fig. 5). Where effectiveness of case management is high, the country specific assumptions for system effectiveness make substantial differences to the estimated EIR distributions (Fig. 4; Table 2; Additional file 2: Figure S2). In these countries, notably Zambia, Tanzania, SĆ£o TomĆ© and Principe, the EIR distribution shifts to the right when the country-specific value of \({\text{E}}_{14}\) is used, reflecting lower prevalence than in a situation with the same EIR pattern, but less effective case management. The estimate of median EIR for these countries is thus much higher when country-specific effectiveness is considered. Conversely, in a few countries, where median prevalence is low, and case management is also poor, the estimated EIR distribution allowing for country-specific effectiveness is shifted slightly to the left (e.g. both South and North Sudan).

Fig. 4
figure 4

Distribution of EIR for each of the 43 countries. Distribution of EIR (including non-endemic areas which are assigned values of 0) for each of the 43 countries. Calculated from MAP using both the empirical model (Method A) (black); the simulation models (Method B) with a common value for access to care (yellow) (\(E_{14} = 0.15)\); and country-specific values of \(E_{14}\)(blue). Countries are indicated by country code

Fig. 5
figure 5

Relationship of estimated average EIR to prevalence at country level. Calculated from MAP using both the empirical model (Method A) (black); the simulation models (Method B) with a common value for access to care, (yellow) (E14 = 0.15); and country-specific values of \({\text{E}}_{14}\)(blue)

The EIR distributions are highly skewed, so that the arithmetic means are much higher than the medians (Figs. 4, 5). Except in some cases where prevalence is very low, the average EIR is higher when there is allowance for treatment (Method B), with a much larger shift in the mean than in the median of the distribution. When country specific E 14 values are used (which are mostly higher than the 15 % shown in yellow), this makes little difference to the mean EIR, but substantial differences to the medians, reflecting the stronger relationship between treatment rates and prevalence when EIR is low, than when EIR is high.

National levels of burden of disease

The OpenMalaria simulations predict that steady state clinical incidence (over all ages) increases linearly with EIR in low transmission settings, tending to plateau at high EIR (with a suggestion, driven by the specific Senegalese data used to parameterize the models for older children and adults [35] that there may be a maximum in the curve at high prevalence). The initial slope is greater when \(E_{14}\) is higher, but the plateau occurs at a similar level of incidence irrespective of effective treatment level. These patterns are a consequence of the age-specific relationships between incidence and EIR shown in Fig. 6.

Fig. 6
figure 6

Models of the relationship between EIR and clinical incidence. Incidence of clinical episodes by EIR in OpenMalaria models with light blue \(E_{14}\) = 0.15 and dark blue \(E_{14}\) = 0.45. The continuous lines indicate the mean prediction of the overall incidence. The shading around the continuous lines indicates the range of predictions made from simulations with different model variants and random number seeds. The dashed lines indicate the incidence of clinical episodes that are treated [31]

When these models are used to infer country-specific incidence of clinical malaria, there is a clear increase in incidence with average prevalence at the country level (Fig. 7), and no plateau is reached because even the countries with highest average transmission have only small populations in the very high EIR categories (Table 4 and Additional file 1: Table S2). The relationships between country-level EIR and estimated clinical malaria incidence are similar, irrespective of whether the EIR is estimated by Method A or Method B. Similarly, Method B estimates similar relationships between country-level EIR and clinical malaria incidence, irrespective of whether a common value, or a country specific estimate is used for the effectiveness of case management.

Fig. 7
figure 7

Predicted incidence of clinical events by national level average EIR. Predicted incidence of clinical events estimated from empirical model (Method A) (black); the simulation models (Method B) with a common value for access to care (yellow) (E14 = 0.15); and country-specific values of \(E_{14}\)(blue). a uncomplicated clinical episodes (cases); b severe malaria episodes; c hospitalizations; d deaths directly attributable to malaria; e hospital deaths; f all malaria deaths (including those with co-morbidities). All rates are expressed as events per 100,000 person years at risk over all ages of hosts. The model and parameters for severe disease and mortality follows Ross et al. [22], with a common hospitalization rate assumed for severe disease across all countries

Country-specific estimates of clinical incidence using country EIR distributions is compared to published malaria cases from the World Malaria Report [14] (Fig. 8; Additional file 2: Figure S5). In general, projections of incidence using EIR derived from method B produces higher predictions than using EIR from method A, but in both cases the simulation models predict substantially more episodes of malaria than the cases reported in the World Malaria Report 2013 (Fig. 8a). There is also a much less steep relationship between the incidence rate and the overall burden (Fig. 8b). This can be explained by the empirical relationships between prevalence and case incidence used by WMR [12], which refer back to field research carried out prior to the widespread use of ACT [33, 36], and therefore do not allow for level of treatment. Moreover, only clinical episodes in children under 5 years of age are considered. The effect of high levels of treatment on reducing prevalence leads to much higher ratios of case-incidence to prevalence ratio than it would be without treatment (Fig. 7). OpenMalaria correctly predicts that in low transmission countries the majority of the clinical burden is in older age groups (Fig. 6). The difference between the methods is particularly evident for low-burden countries Namibia and Botswana, for which very low case numbers are reported by the WMR, with estimates based on adjustments to surveillance data.

Fig. 8
figure 8

Numbers of episodes per annum estimated using different approaches. Numbers of clinical cases; b incidence rate (episodes or cases per 100,000 person-year). Black points: based on EIR estimates calculated using Method A; yellow points: based on EIR estimates made using Method B with a common value for access to care (\(E_{14} = 0.15)\); blue points: based on EIR estimates made using Method B with country-specific values of \(E_{14}\)(blue). The diagonal line corresponds to a 1:1 relationship; the horizontal and vertical lines represent minimum and maximum ranges

Discussion

It has been incontrovertible since Laveranā€™s first studies of the malaria parasite [37] that effective treatment of clinical malaria results in clearance of blood-stage parasites. Treatment lowers the overall prevalence associated with malaria (or other parasites [38]), the infectiousness of the human population and the transmission level [31], both of which synergize with effects of other interventions on transmission. Most immediately, effective treatment reduces the length of illness and the incidence adverse outcomes, including severe disease, neurological sequelae, and death. This reduces the burden of disease that can potentially be averted by other interventions. All these effects need to be considered in estimating current burden of disease, in analyses of the impact of case management and treatment on the burden of disease, and analyses of treatment modifies the public health impact achievable with other curative and preventive interventions. Results presented here clearly indicate that incorporating the dynamic effects of treatment is essential for valid estimation of EIR, of clinical incidence itself, and of downstream outcomes including the incidence of severe disease and mortality rates, with substantively differences in estimates when included or excluded.

Overall, the model-based method proposed in this work (method B) provides estimates of transmission intensity, as measured by EIR, that are somewhat higher than those estimated by method A [7], especially in low endemicity countries and where case management is relatively effective. The downstream country level clinical incidence estimates are also higher than those previously reported in the World Malaria Report [13]. At the country level, allowing for uncertainty in the inputs makes a substantial difference to average values of both EIRs and disease rates, as a result of the skewness of their distributions. This means that incorporating uncertainty and spatial variation into the estimation has important consequences for both burden estimates and prediction of average health impacts of interventions, which in general vary non-linearly with EIR. The OpenMalaria models also predict, as one would expect, that the effectiveness of uncomplicated malaria treatment has substantial impact on the incidence of severe disease and malaria mortality.

Preventive interventions like insecticide treated nets (ITNs), which affect prevalence only via their impact on exposure, do not change the relationship between exposure and prevalence. Consequently, coverages of preventive interventions can be useful covariates for estimating EIR or prevalence surfaces where direct measurements are sparse, but the coverage of these interventions are not directly relevant when making estimations of disease burden from prevalence. In contrast, treatment of malaria reduces the prevalence at a given level of EIR, by preventing infections from persisting, thus modifying the relationship between the two metrics (Fig. 1). So the same prevalence can result from very different average exposures depending on the level of treatment, and the effective coverage of case management (like the degree of transmission heterogeneity [2]) should be taken account in modeling the relationships between EIR and prevalence.

Nevertheless, at least in high endemicity settings, prevalence remains the best measure on which to base geographically specific models of malaria transmission. This is because prevalence data are actively collected based on representative sampling of populations, are widely available, have been compiled into publicly accessible databases [7, 39], and have been analysed using geostatistical models to produce high resolution maps of the distribution of infection in space [7, 40]. In most sub-Saharan African countries prevalence is therefore likely to remain the main metric used in deciding when and where to distribute or target interventions. In low transmission settings such as those in Asia, Latin America, and selected African countries the annual parasite index (API) rather than the prevalence is the main metric used for monitoring and evaluation, and WMR has estimated burden in these countries using an API-based algorithm [12]. Prevalence-EIR-treatment relationships in such low transmission settings can be captured by relatively simple empirical mathematical models [4]. However in areas of moderate or high transmission it is important to allow for effects of superinfection and natural immunity, and thus mechanistic models that account for dynamics of immunity are needed.

The use of simulation models that take both prevalence and treatment rates as inputs provides a generalizable way of generating national level estimates of transmission and disease burden, applicable across the range of transmission intensities. This generalizability will be important for monitoring progress as malaria is further controlled to the point where measurement of API becomes the main metric used by many more country programmes. The approach will capture in a natural way the transitions between the different metrics, and the age shifts in the pattern of disease where transmission rates change [41, 42]. The approach can be made more robust by employing a larger ensemble including other simulation models with different assumptions about transmission heterogeneity, immunity, and pathogenesis [10, 43]. For the method to provide the best estimates of malaria attributable mortality, geographical variation in access to appropriate in-patient treatment of severe disease also needs to be taken into account.

Previous methodologies for estimating burden have applied both estimates of intervention protective efficacy derived from meta-analyses of controlled trials and/or household survey data, leading to circular reasoning. Local variability has also been ignored [1], in particular variations in access, compliance, or adherence, and also the medium- and long-term dynamics resulting from intervention-induced reductions in transmission, which include shifts of disease into older age groups [41, 42]. The burden estimation procedures proposed in this paper will allow empirical analysis of the relationships between intervention coverage and burden independently of field trial results and conditional on all these factors. This will provide a basis for assessing the impacts of both preventive and curative interventions on an equivalent basis, ensuring correct attribution of the effects of different interventions. The method can be extended to give time-dependent estimates of burden by using time-period specific input data. By linking these to intervention coverage, this will provide valid estimates of intervention impacts in time and space. Although results are presented only at country level in this work, this methodology can, in principle, be applied to any level of spatial aggregation. However, applying the approach to data disaggregated in smaller spatial units would raise additional methodological issues, as the simulation models are parameterized mainly using village-level data.

This paper demonstrates the dual importance of capturing the effects of treatment when estimating disease burden based on infection prevalence: to both improve the accuracy of those estimates and to correctly quantify the impact of treatment on reduced malaria transmission and illness. These insights are currently being incorporated into a revised WHO methodology that will lead to more refined burden estimates and ultimately better information for national and international malaria control decision-making processes.

References

  1. Gething PW, Battle KE, Bhatt S, Smith DL, Eisele TP, Cibulskis RE, et al. Declining malaria in Africa: improving the measurement of progress. Malar J. 2014;13:39.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  2. Ross A, Smith T. Interpreting malaria age-prevalence and incidence curves: a simulation study of the effects of different types of heterogeneity. Malar J. 2010;9:132.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  3. Beier JC, Killeen G, Githure JI. Short report: entomologic inoculation rates and Plasmodium falciparum malaria prevalence in Africa. Am J Trop Med Hyg. 1999;61:109ā€“13.

    CASĀ  PubMedĀ  Google ScholarĀ 

  4. Yukich J, Briet O, Bretscher MT, Bennett A, Lemma S, Berhane Y, et al. Estimating Plasmodium falciparum transmission rates in low-endemic settings using a combination of community prevalence and health facility data. PLoS One. 2012;7:e42861.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  5. Stuckey EM, Smith T, Chitnis N. Seasonally dependent relationships between indicators of malaria transmission and disease provided by mathematical model simulations. PLoS Comput Biol. 2014;10:e1003812.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  6. Smith DL, Dushoff J, Snow RW, Hay SI. The entomological inoculation rate and Plasmodium falciparum infection in African children. Nature. 2005;438:492ā€“5.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  7. Gething PW, Patil AP, Smith DL, Guerra CA, Elyazar IR, Johnston GL, et al. A new world malaria map: plasmodium falciparum endemicity in 2010. Malar J. 2011;10:378.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  8. Gemperli A, Vounatsou P, Sogoba N, Smith T. Malaria mapping using transmission models: application to survey data from Mali. Am J Epidemiol. 2006;163:289ā€“97.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  9. Gemperli A, Sogoba N, Fondjo E, Mabaso M, Bagayoko M, Briet OJ, et al. Mapping malaria transmission in West and Central Africa. Trop Med Int Health. 2006;11:1032ā€“46.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  10. Griffin JT, Hollingsworth TD, Okell LC, Churcher TS, White M, Hinsley W, et al. Reducing Plasmodium falciparum malaria transmission in Africa: a model-based evaluation of intervention strategies. PLoS Med. 2010;7:e1000324.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  11. Hay SI, Okiro EA, Gething PW, Patil AP, Tatem AJ, Guerra CA, et al. Estimating the global clinical burden of Plasmodium falciparum malaria in 2007. PLoS Med. 2010;7:750.

    Google ScholarĀ 

  12. Cibulskis RE, Aregawi M, Williams R, Otten M, Dye C. Worldwide incidence of malaria in 2009: estimates, time trends, and a critique of methods. PLoS Med. 2011;8:e1001142.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  13. WHO. World Malaria Report 2013. Geneva: World Health Organization; 2014.

    Google ScholarĀ 

  14. WHO. World Malaria Report 2012. Geneva: World Health Organization; 2013.

    Google ScholarĀ 

  15. Griffin JT, Ferguson NM, Ghani AC. Estimates of the changing age-burden of Plasmodium falciparum malaria disease in sub-Saharan Africa. Nat Commun. 2014;5:3136.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  16. Stuckey EM, Smith TA, Chitnis N. Estimating malaria transmission through mathematical models. Trends Parasitol. 2013;29:477ā€“82.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  17. Smith T, Killeen GF, Maire N, Ross A, Molineaux L, Tediosi F, et al. 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. 2006;75:1ā€“10.

    PubMedĀ  Google ScholarĀ 

  18. Chitnis N, Hardy D, Smith T. A periodically-forced mathematical model for the seasonal dynamics of malaria in mosquitoes. Bull Math Biol. 2012;74:1098ā€“2024.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  19. Smith T, Maire N, Dietz K, Killeen GF, Vounatsou P, Molineaux L, et al. Relationship between the entomologic inoculation rate and the force of infection for Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006;75:11ā€“8.

    PubMedĀ  Google ScholarĀ 

  20. Maire N, Smith T, Ross A, Owusu-Agyei S, Dietz K, Molineaux L. A model for natural immunity to asexual blood stages of Plasmodium falciparum malaria in endemic areas. Am J Trop Med Hyg. 2006;75:19ā€“31.

    PubMedĀ  Google ScholarĀ 

  21. Ross A, Killeen GF, Smith T. Relationships of host infectivity to mosquitoes and asexual parasite density in Plasmodium falciparum. Am J Trop Med Hyg. 2006;75(Suppl 2):32ā€“7.

    PubMedĀ  Google ScholarĀ 

  22. Ross A, Maire N, Molineaux L, Smith T. An epidemiologic model of severe morbidity and mortality caused by Plasmodium falciparum. Am J Trop Med Hyg. 2006;75:63ā€“73.

    PubMedĀ  Google ScholarĀ 

  23. Smith T, Ross A, Maire N, Rogier C, Trape JF, Molineaux L. An epidemiologic model of the incidence of acute illness in Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006;75:56ā€“62.

    PubMedĀ  Google ScholarĀ 

  24. Smith T, Ross A, Maire N, Chitnis N, Studer A, Hardy D, et al. Ensemble modeling of the likely public health impact of a pre-erythrocytic malaria vaccine. PLoS Med. 2012;9:e1001157.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  25. Galactionova K, Tediosi F, Savigny DD, Smith TA, Tanner M. Effective coverage and systems effectiveness for malaria case management in Sub-Saharan African countries. PLoS One. 2015;10:e0127818.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  26. Malaria Atlas Project. http://www.map.ox.ac.uk/. Accessed 15 Apr 2015.

  27. Tediosi F, Hutton G, Maire N, Smith TA, Ross A, Tanner M. Predicting the cost-effectiveness of introducing a pre-erythrocytic malaria vaccine into the expanded program on immunization in Tanzania. Am J Trop Med Hyg. 2006;75:131ā€“43.

    PubMedĀ  Google ScholarĀ 

  28. Crowell V, Yukich J, Briet OJ, Ross A. A novel approach for measuring the burden of uncomplicated Plasmodium falciparum malaria: application to data from Zambia. PLoS One. 2013;8:e57297.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  29. Balk D, Deichmann U, Yetman G, Pozzi F, Hay S, Nelson A. Determining global population distribution: methods, applications and data. Adv Parasitol. 2006;62:119ā€“56.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  30. Center for International Earth Science Information Network (CIESIN) CU: Global Rural-Urban Mapping Project (GRUMP): Urban extents. New York: Palisades; 2004.

  31. Tediosi F, Maire N, Smith T, Hutton G, Utzinger J, Ross A, et al. An approach to model the costs and effects of case management of Plasmodium falciparum malaria in sub-saharan Africa. Am J Trop Med Hyg. 2006;75:90ā€“103.

    PubMedĀ  Google ScholarĀ 

  32. Snow R, Craig M, Deichmann U, Marsh K. Estimating mortality, morbidity and disability due to malaria among Africaā€™s non-pregnant population. Bull World Health Organ. 1999;77:624ā€“40.

    PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  33. Roca-Feltrer A, Carneiro I, Armstrong Schellenberg JR. Estimates of the burden of malaria morbidity in Africa in children under the age of 5Ā years. Trop Med Int Health. 2008;13:771ā€“83.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  34. WHO. World malaria report 2008. Geneva: World Health Organization; 2008.

    Google ScholarĀ 

  35. Trape JF, Rogier C. Combating malaria morbidity and mortality by reducing transmission. Parasitol Today. 1996;12:236ā€“40.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  36. Snow R, Craig M, Newton C, Steketee RW. The public health burden of Plasmodium falciparum malaria in Africa: Deriving the numbers. Bethesda: Fogarty International Center, National Institutes of Health; 2003.

    Google ScholarĀ 

  37. Laveran A. Nature parasitaire des accidents de lā€™impaludisme: description dā€™un nouveau parasite trouvĆ© dans le sang des malades atteints de fiĆØvre palustre. Paris: J. B. Bailliere; 1881.

    Google ScholarĀ 

  38. Anderson RM, May RM. Helminth infections of humans: mathematical models, population dynamics, and control. Adv Parasitol. 1985;24:1ā€“101.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  39. Colloboration MARA. Towards an Atlas of Malaria Risk in Africa. Durban: MARA/ARMA; 1998.

    Google ScholarĀ 

  40. Gosoniu L, Vounatsou P, Sogoba N, Smith T. Bayesian modelling of geostatistical malaria risk data. Geospatial Health. 2006;1:127ā€“39.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  41. Woolhouse ME. Patterns in parasite epidemiology: the peak shift. Parasitol Today. 1998;14:428ā€“34.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  42. Smith T, Hii J, Genton B, Muller I, Booth M, Gibson N, et al. Associations of peak shifts in age-prevalence for human malarias with bed net coverage. Trans R Soc Trop Med Hyg. 2001;95:1ā€“6.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  43. Eckhoff P. Mathematical models of within-host and transmission dynamics to determine effects of malaria interventions in a variety of transmission settings. Am J Trop Med Hyg. 2013;88:817ā€“27.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

Download references

Authorsā€™ contributions

MAP, CB, NM and TAS designed the experiments and analyzed results. PPR, OJTB, DLS, and PWG contributed to the methodology and carried out data analyses. MAP and TAS drafted the manuscript. TAS, MAP and NM conceived of the study. All authors read and approved the final manuscript.

Acknowledgements

The authors would like to thank Michael Tarantino and Erin Stuckey for the help with running simulations. The authors would also like to thank the many volunteers who made their computers available to run simulations via malariacontrol.net.

Compliance with ethical guidelines

Competing interests The authors declare that they have no competing interests.

Financial disclosure

MAP and TAS acknowledge funding by the Bill and Melinda Gates Foundation (#OPP1032350) and PATH-Malaria Vaccine Initiative (MVI). PWG is a Career Development Fellow (#K00669X) jointly funded by the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement and also receives support from the Bill and Melinda Gates Foundation (#OPP1068048). DLS acknowledges funding from the Bill and Melinda Gates Foundation (#OPP1110495). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Melissa A. Penny.

Additional files

12936_2015_864_MOESM1_ESM.pdf

Additional file 1. This file includes additional Tables that support and expand some of the results in the main text, but whose inclusion would detract from the main argument.

12936_2015_864_MOESM2_ESM.pdf

Additional file 2. This file includes additional Figures of results that support and expand some of the results in the main text, but whose inclusion would detract from the main argument.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Penny, M.A., Maire, N., Bever, C.A. et al. Distribution of malaria exposure in endemic countries in Africa considering country levels of effective treatment. Malar J 14, 384 (2015). https://doi.org/10.1186/s12936-015-0864-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-015-0864-3

Keywords