- Open Access
Interpreting malaria age-prevalence and incidence curves: a simulation study of the effects of different types of heterogeneity
Malaria Journal volume 9, Article number: 132 (2010)
Individuals in a malaria endemic community differ from one another. Many of these differences, such as heterogeneities in transmission or treatment-seeking behaviour, affect malaria epidemiology. The different kinds of heterogeneity are likely to be correlated. Little is known about their impact on the shape of age-prevalence and incidence curves. In this study, the effects of heterogeneity in transmission, treatment-seeking and risk of co-morbidity were simulated.
Simple patterns of heterogeneity were incorporated into a comprehensive individual-based model of Plasmodium falciparum malaria epidemiology. The different types of heterogeneity were systematically simulated individually, and in independent and co-varying pairs. The effects on age-curves for parasite prevalence, uncomplicated and severe episodes, direct and indirect mortality and first-line treatments and hospital admissions were examined.
Different heterogeneities affected different outcomes with large effects reserved for outcomes which are directly affected by the action of the heterogeneity rather than via feedback on acquired immunity or fever thresholds. Transmission heterogeneity affected the age-curves for all outcomes. The peak parasite prevalence was reduced and all age-incidence curves crossed those of the reference scenario with a lower incidence in younger children and higher in older age-groups. Heterogeneity in the probability of seeking treatment reduced the peak incidence of first-line treatment and hospital admissions. Heterogeneity in co-morbidity risk showed little overall effect, but high and low values cancelled out for outcomes directly affected by its action. Independently varying pairs of heterogeneities produced additive effects. More variable results were produced for co-varying heterogeneities, with striking differences compared to independent pairs for some outcomes which were affected by both heterogeneities individually.
Different kinds of heterogeneity both have different effects and affect different outcomes. Patterns of co-variation are also important. Alongside the absolute levels of different factors affecting age-curves, patterns of heterogeneity should be considered when parameterizing or validating models, interpreting data and inferring from one outcome to another.
A myriad of differences between individuals in a community affect the epidemiology of malaria. These differences arise in roughly four ways: heterogeneity of transmission; biological heterogeneity of the host in susceptibility and response to malaria infection[1–4]; heterogeneity in host behaviour, including quality of housing, use and knowledge of protective measures such as insecticide-treated nets (ITN) and treatment-seeking strategies; and heterogeneity in the risk of co-morbidity and malnutrition.
The different sources of heterogeneity are unlikely to be independent. A factor, such as socio-economic status (SES), may be associated with several different heterogeneities. SES can influence transmission intensity through quality of housing, knowledge of protection measures and use of ITNs [5–8]. Risks of co-morbidity and malnutrition may be associated with SES [9, 10]. SES may influence treatment-seeking behaviour . Amongst those caring for children, treatment-seeking [7, 11–13] and knowledge of danger signs [10, 11] were worse in poorer families. Such inequalities tend to persist over time . Thus, covariance between different types of heterogeneity is plausible although the degree to which it occurs is likely to vary between sites.
Differences among individuals can lead to important patterns at the community level [15, 16]. Different types of heterogeneity may have implications for both the interpretation of field data and for formulating and calibrating models of the effects of interventions, or alternatively may have little impact and could be safely ignored. Co-variation may describe sub-groups within a population, who perhaps respond differently to an intervention (for example ) or who are reached neither by the health system nor by any of several interventions.
Assessing the effects of different kinds of heterogeneity using field datasets is not practical due to difficulties in collecting the relevant data and isolating the effects of each of the heterogeneities. Therefore, heterogeneity has been investigated using mathematical models. There are few examples of heterogeneity modelling studies specific to malaria (for example [15, 18–25]) and very few incorporating multiple types of heterogeneity [23, 24]. One reason may have been the models themselves. To study the effects of multiple types of heterogeneity, a model must be sufficiently comprehensive, dynamic to allow secondary effects, and be able to easily incorporate different kinds of heterogeneity.
A recently published model of Plasmodium falciparum malaria epidemiology  satisfies these criteria. It includes processes for infection, parasitaemia, acquired immunity, infectivity, morbidity and mortality and case-management. It is dynamic, allowing feedback effects such as the effects of high treatment coverage on transmission intensity. It is also individual-based, which provides a flexible framework to conveniently incorporate heterogeneity [27–29].
This model is adopted as a base to investigate the effects of three different types of heterogeneity: transmission (including host behaviours such as bed net use and housing), co-morbidity (as a trigger for severe malaria and indirect mortality and reflecting the individual's general health and nutrition) and treatment-seeking (the chances of obtaining effective treatment). These heterogeneities are simulated singly, and in independent and co-varying pairs. This paper focuses on the effects of heterogeneity on fundamental measures of the malaria burden, age-prevalence and incidence curves. These curves are known to vary by the severity of the outcome, degree of seasonality and the overall levels of transmission intensity (reviewed in ).
The stochastic model of malaria epidemiology has been described in detail elsewhere . Briefly, a simulated population of humans are updated at each five-day time-step via components representing new infections, parasite densities, acquired immunity, uncomplicated and severe episodes, direct and indirect malaria mortality, and infectivity to mosquitoes. The course of parasite densities over an infection are described by averaged empirical data . Immunity to asexual parasites is derived from a combination of cumulative exposure to both inoculations and parasite densities, and maternal immunity . The inclusion of acquired immunity allows us to model potential effects of heterogeneity on immunity through exposure being averted or increased. The probability of a clinical attack of malaria depends on the current parasite density and a pyrogenic threshold . The pyrogenic threshold responds dynamically to recent parasite load, increasing or saturating through exposure to parasites and decaying with time, and thus is individual- and time-specific. Severe malaria can arise in two ways, either as a result of overwhelming parasite density or through uncomplicated malaria with concurrent non-malaria co-morbidity . Treatment of malaria is based on the case management model of Tediosi and colleagues  which involves implementing case management trees for severe cases, and uncomplicated cases with and without recent history of treatment. In the current simulations, we assumed first-line treatment of uncomplicated episodes to be an effective short-acting drug such as artemisinin-combination therapy (ACT) with an average probability of treatment at any one five-day time step of 20%. This value is higher than previously used [34–37]; it was raised arbitrarily to allow the effects of treatment heterogeneity to show. The treatment probability for severe cases warranting hospital admission was 48% . Mortality in this model can be either direct (following severe malaria) or indirect (uncomplicated malaria in conjunction with co-morbidity, or during the neonatal period as a result of maternal infection) . This co-morbidity is simulated as occurring with an age-dependent probability described by a hyperbolic function .
The parameter values for the base model were estimated by fitting to data from a total of 61 malaria field studies of various different aspects of malaria epidemiology , and are reported elsewhere .
Types of heterogeneity
Some individual differences are inherent to the base model and arise due to host age and an individual's propensity for high parasite densities . In addition, stochasticity occurs throughout the model.
We extended the base model by introducing an additional three types of heterogeneity: (i) transmission, (ii) co-morbidity and (iii) treatment-seeking probability. This study does not consider heterogeneity in the biological response of the host (except for the host dependent heterogeneity in parasite densities incorporated into the base model). The model processes are shown in Figure 1. Each simulated individual is assigned a status for each of the three kinds of heterogeneity at birth, which they carry throughout their life. The structure of the heterogeneities is that 50% of the population are assigned to each of the high and low status categories, with baseline values multiplied by either 1.8 or 0.2. Hence, when assuming heterogeneity in transmission, the age-adjusted entomological inoculation rate (EIR) for individual i at time t takes the value: E L (i,t) = 0.2E a (i,t) for the low exposure half of the population, where E a (i, t) is the age-adjusted EIR in the base model, and E H (i,t) = 1.8E a (i,t) in the high exposure half of the population. Similarly, in models of treatment heterogeneity, the probability that either an uncomplicated or severe attack is treated during any one simulation step takes a value of 1.8 (or 0.2) times that in the base model; in models of heterogeneous co-morbidity, two age-dependent probabilities are multiplied by 1.8 (or 0.2), the probability that a clinical attack becomes severe due to co-morbidity in any one five-day time step and the probability that, following a malaria episode which weakens the host, a non-malaria illness occurs leading to an indirect death.
In simulations of co-varying heterogeneities, the same 50% of the population are at greater risk on each of the three dimensions. Thus low access to treatment is paired with high numbers of infective bites and high risk of co-morbidity.
The present analyses make use of estimates of the remaining parameters in the model derived assuming no heterogeneity in the three variables that are the subject of this study . Results of re-estimation of these parameters conditional on different assumptions about heterogeneity will be reported elsewhere.
A reference scenario was defined using the input values in Table 1. This study does not report simulations of interventions other than case-management. The three types of heterogeneity were added to the base model and their effects systematically simulated singly, in independent pairs, and in co-varying pairs. The outcomes were parasite prevalence, uncomplicated and severe clinical episodes, direct and indirect mortality, first-line treatments and hospital admissions. The age-prevalence and incidence curves were described using the age at, and height of, the peak. To investigate the effect of transmission intensity on the effects of heterogeneity, this set of simulations was repeated for settings with 1, 6 and 200 infectious bites per person per year. The seasonality in all scenarios followed that of Namawala, Tanzania  (Table 1). Average annual transmission levels were held constant so that effects were not simply due to changes in transmission.
The effect of single heterogeneities
Heterogeneity in transmission affected all of the outcomes (Additional File 1: Table S1). The predicted prevalence of parasitaemia was reduced (Figure 2a). For the other outcomes, the incidence rates were reduced in younger age groups compared to the reference scenario and increased for older ages (Figure 2b).
There was no apparent overall effect of heterogeneity in co-morbidity risk on the age-curves. However, taken separately, low and high co-morbidity levels did show large effects on some outcomes suggesting that the combined effects cancel out in the scenario with heterogeneity (Figure 3). The outcomes affected were those directly affected by the action of co-morbidity risk: severe episodes (and hence direct mortality), hospital admissions, and indirect mortality. Co-morbidity does not influence acquired immunity in the model and so there were no feedback effects.
Probability of seeking treatment
The strongest effect of heterogeneity in treatment-seeking behaviour was a reduction in the incidence of both first-line treatments and hospital admissions (Figure 4). Taken separately, low and high levels of treatment-seeking affected all outcomes: low levels lead to increased prevalence and increased incidence of episodes and deaths. When combined in the scenario with heterogeneity in treatment-seeking levels however, the overall effect was limited, with slight increases in the peak for parasite prevalence and uncomplicated episodes (Additional File 1: Table S1).
In general, the strongest effects of the single heterogeneities were reserved for outcomes directly affected by the action of the heterogeneity. Feedback via acquired immunity produced more modest effects. For some outcomes, the effects of high and low values of the heterogeneity appeared to cancel out at the community level.
Independent and co-varying pairs of heterogeneities
The results suggest that the two independent types of heterogeneities together produce additive effects on the age-curves (Additional File 2: Table S2).
The results for sources of heterogeneity which co-vary were more variable. Changes in comparison to the independently varying case occurred only if the outcome measure was affected by both kinds of heterogeneity individually, either apparent or masked by high and low values cancelling out. Changes were observed in outcomes which were directly affected by the action of the heterogeneity. Striking differences compared to independent heterogeneities were produced for hospital admissions and direct mortality (co-morbidity and treatment-seeking: Figure 5) and first-line and hospital admissions (transmission and treatment-seeking: Figure 6). More modest changes were produced for other outcomes (Additional File 2: Table S2).
There were no changes for outcomes which were not affected by both heterogeneities individually. Prevalence, uncomplicated episodes and first-line treatments were unaffected by heterogeneity in co-morbidity and no changes were observed for the pairs co-morbidity and treatment-seeking, and co-morbidity and transmission. There were also no apparent changes to the age-curves for prevalence, uncomplicated and severe episodes, and direct and indirect mortality (transmission and treatment-seeking). Although both types of heterogeneity either affected these outcomes or masked effects at high and low levels, the effect of treatment is via dynamic feedback on acquired immunity. Outcomes which are not directly influenced by the action of the heterogeneities are less affected in our simulations.
The effect of heterogeneity in all three variables was simulated for completeness (Additional File 1: Table S1), however interpreting the results is difficult. In general, age-curves for outcomes affected by only two of the heterogeneities were similar to those produced by the pair of heterogeneities. Where all three had individual effects, the independent and co-varying triple heterogeneities could have effects which were very different to each other.
The effect of heterogeneity for different transmission intensities
These simulations were repeated for different transmission intensities (EIR 1, 6 and 200) keeping the same seasonal pattern. The age-curves shift due to the effects of transmission intensity, but the direction and strength of the impact of the different types and combinations of heterogeneity were similar to those of the reference transmission of 21 infectious bites per person per year.
This study describes the simulated effects of heterogeneity in transmission intensity, co-morbidity risk and treatment-seeking behaviour on the age-curves of a range of measures of the burden of malaria. Although this study is specific to malaria, our findings add more generally to studies of infectious diseases indicating that it is not only the absolute level of a variable, but also the presence of heterogeneity [15, 22, 25, 40] and co-variation [23, 28], which can produce important patterns at the community level.
The results show that the different types of heterogeneity have effects on different outcomes with large effects reserved for outcomes directly affected by the action of the heterogeneity, rather than dynamic feedback via acquired immunity. Transmission heterogeneity affected the age-curves for all outcomes. The peak parasite prevalence was reduced and all age-incidence curves crossed those of the reference scenario with a lower incidence in younger children and higher in older age-groups. Heterogeneity in the probability of seeking treatment reduced the peak incidence of first-line treatment and hospital admissions. Heterogeneity in co-morbidity risk showed little overall effect, but high and low values cancelled out for outcomes influenced by its action. Independently varying pairs of heterogeneities produced additive effects. More variable results were produced for co-varying pairs, with striking differences compared to independent pairs for some outcomes which were affected by both heterogeneities individually.
These findings have implications for both the interpretation of age-curves and their use in analysis and modelling. In pointing to where the different types of heterogeneity change the shape of the age-curves, this study also indicates where there is no need for concern. Since the greatest effects of single heterogeneities were reserved for outcomes directly affected or subsequent to their action, only modest effects would be expected for outcomes less directly linked unless transmission intensity is affected. Outcomes which were unaffected by single heterogeneities do not change substantially when there are concurrent heterogenous variables, either independent or co-varying.
The results illustrate the effect of heterogeneity on the ability to infer from one outcome to another. The utility of passive case detection to estimate the burden of clinical disease depends on assumptions about the underlying distribution of treatment-seeking behaviour. Similarly, the use of clinical episodes to infer effects on mortality depends on heterogeneity in the risk of co-morbidity and treatment-seeking. Patterns of heterogeneity are also important when estimating transmission from prevalence or clinical data.
A logical consequence of the findings is that estimating parameter values by fitting a model to data from a single field site falsely assuming either homogeneity or independently varying variables is liable to produce incorrect values, as found elsewhere [15, 23, 25].
This study has several limitations. The adopted representation of the patterns of heterogeneity was very crude. Whilst the strength of the impact will differ, the conclusions are not dependent on the values chosen for the high and low risk groups, on the distribution of individuals between groups or on whether there is a gradual increase in risk across individuals as long as the overall values are equal in the heterogeneity and comparison scenarios. Thus, although the scenarios do not exhibit realistic patterns, the findings nevertheless provide insights on the substantive effects of different types of heterogeneity. Likewise, there was no intention to realistically portray SES. There are complex pathways between poverty, economic activity, health-care seeking and malaria [7, 13, 41]. The scope of this study is limited to the incorporation of simple heterogeneity in three variables which are plausibly related to SES into a model of malaria epidemiology. It was also assumed that the individual's relative levels for transmission, co-morbidity risk and treatment-seeking are constant throughout their life. This simplification ignores the impact of malaria on poverty [13, 42] and mobility . More sophisticated simulations may address the need to disaggregate SES  and identify interventions of particular benefit to the poor (such as ).
Assumptions were made about the likely direction of co-variation, matching high transmission to a high risk of co-morbidity and low probability of seeking treatment. In different settings and ecotypes, these directions may be reversed. Individuals are assumed to be distributed spatially at random. Spatial patterns could be better incorporated by grouping individuals into households.
Objective criteria are lacking for determining whether two age-curves differ. However, the patterns observed were clear-cut. In this study, heterogeneity in transmission had the strongest effects on the age-curves. It is plausible that this would also be the case in many real settings. However, this study does not use use a realistic patterns and degrees of heterogeneity, and so conclusions cannot be drawn about the relative importance of the three heterogeneities on the shape of the age-curves.
The average annual level of transmission was fixed to remain constant so that the effects observed were solely due to heterogeneity when compared to the reference scenario. However, heterogeneity may also alter the absolute level of transmission intensity via dynamic feedback to the infectious reservoir.
The base model used is comprehensive and individual-based providing a flexible framework for unravelling the effects of different heterogeneities. Limitations of the model components are discussed elsewhere [26, 31–33, 46–49]. Some assumptions are especially relevant to this study. Non-malaria co-morbidity is assumed to prompt an acute episode to lead to either a severe episode or an indirect death, and the risk of co-morbidity is assumed to be age-dependent. This is reflected in the age-curves for both the reference scenario and the inclusion of heterogeneity in co-morbidity risk. The base model assumes constant probabilities for seeking treatment within a five-day period for all individuals with either acute or severe episodes. However, treatment-seeking may be age-dependent, due to differences in the recognition of fevers or perceived need for treatment in adults, children and infants and also to the tendency for different symptoms to manifest at different ages such as severe malarial anaemia in young children . The five-day time step constrains the model components for both treatment-seeking and case-management to be very simple. The time steps are currently being shortened to one day and a more sophisticated case-management model is in development.
Effects of heterogeneity on outcomes other than age-curves were beyond the scope of this paper. Heterogeneity is likely to have important effects on receptivity and elimination [15, 19–21], on individual differences and equality, and on the impact of different interventions [22–24]. Some interventions may lessen the differences between individuals in a population, but the effects of reducing heterogeneity are not yet known. In addition to further simulation studies, future work should also extend the analysis of available field data to describe the pattern of heterogeneity in different settings. Innovative methods are needed to estimate sources of heterogeneity that are difficult to measure, such as for micro-transmission [25, 51, 52], and multivariate analysis estimating the degree and patterns of co-variation would enable more realistic scenarios to be considered.
Heterogeneity is not a single entity: different types and combinations of heterogeneity lead to different effects. Interpretation and use of age-prevalence or age-incidence curves should involve consideration of the nature and extent of different heterogeneities in the populations being analysed.
Hill AV, Allsopp CE, Kwiatkowski D, Anstey NM, Twumasi P, Rowe P, Bennett S, Brewster D, McMichael AJ, Greenwood BM: Common west African HLA antigens are associated with protection from severe malaria. Nature. 1991, 352: 595-600. 10.1038/352595a0.
Mackinnon M, Mwangi T, Snow R, Marsh K, Williams T: Heritability of malaria in Africa. PLoS Med. 2005, 2: 340-10.1371/journal.pmed.0020340.
Modiano D, Petrarca V, Sirima BS, Nebie I, Diallo D, Esposito F, Coluzzi M: Different response to Plasmodium falciparum malaria in West African sympatric ethnic groups. Proc Natl Acad Sci USA. 1996, 93: 13206-13211. 10.1073/pnas.93.23.13206.
Williams TN, Mwangi TW, Wambua S, Alexander N, Kortok M, Snow RW, Marsh K: Sickle cell trait and the risk of Plasmodium falciparum malaria and other childhood diseases. J Inf Dis. 2005, 192: 178-186. 10.1086/430744.
Nathan R, Masanja H, Mshinda H, Schellenberg JA, de Savigny D, Lengeler C, Tanner M, Victora CG: Mosquito nets and the poor: can social marketing redress inequities in access?. Trop Med Int Health. 2004, 9: 1121-1126. 10.1111/j.1365-3156.2004.01309.x.
Keating J, Macintyre K, Mbogo C, Githure J, Beier JC: Self-reported malaria and mosquito avoidance in relation to household risk factors in a Kenyan coastal city. J Biosoc Sci. 2005, 37: 767-771. 10.1017/S0021932005007182.
Worrall E, Basu S, Hanson K: Is malaria a disease of poverty? A review of the literature. Trop Med Int Health. 2005, 10: 1047-1059. 10.1111/j.1365-3156.2005.01476.x.
Onwujekwe O, Malik el-F, Mustafa S, Mnzavaa A: Do malaria preventive interventions reach the poor? Socioeconomic inequities in expenditure on and use of mosquito control tools in Sudan. Health Policy Plan. 2005, 21: 10-16. 10.1093/heapol/czj004.
Masanja H, Schellenberg J, de Savigny D, Mshinda H, Victora C: Impact of Integrated Management of Childhood Illness on inequalities in child health in rural Tanzania. Health Policy Plan. 2005, i77-i84. 10.1093/heapol/czi054. Suppl 1
World Health Organisation Western Pacific Region: Integrating poverty and health into health programmes: a sourcebook for health professionals: module on malaria. 2005
Armstrong Schellenberg J, Victora C, Mushi A, de Savigny D, Schellenberg D, Mshinda H, Bryce J, for the Tanzania IMCI MCE baseline household survey study group: Inequalities among the very poor: health care for children in rural southern Tanzania. Lancet. 2003, 361: 561-566. 10.1016/S0140-6736(03)12515-9.
Njau J, Goodman C, Kachur S, Palmer N, Khatib R, Abdulla S, Mills A, Bloland P: Fever treatment and household wealth: the challenge posed for rolling out combinations therapy for malaria. Trop Med Int Health. 2006, 11: 299-313. 10.1111/j.1365-3156.2006.01569.x.
Chuma J, Gilson L, Molyneux C: Treatment-seeking behaviour, cost burdens and coping strategies among rural and urban households in Coastal Kenya: an equity analysis. Trop Med Int Health. 2007, 12: 673-686.
Power C, Matthews S: Origins of health inequalities in a national population sample. Lancet. 1997, 350: 1584-1589. 10.1016/S0140-6736(97)07474-6.
Anderson RM, May RM: Infectious diseases of humans: dynamics and control. 1991, Oxford: Oxford University Press
Grimm V: Ten years of individual-based modelling in ecology: what have we learned and what could we learn in future?. Ecolog Modelling. 1999, 115: 129-148. 10.1016/S0304-3800(98)00188-4.
Mbaye A, Richardson K, Balajo B, Dunyo S, Shulman C, Milligan P, Greenwood B, Walraven G: A randomised, placebo-controlled trial of intermittent preventive treatment with sulphadoxine-pyrimethamine in Gambian multigravidae. Trop Med Int Health. 2006, 11: 992-1002. 10.1111/j.1365-3156.2006.01649.x.
Smith DL, McKenzie FE, Snow RW, Hay SI: Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biol. 2007, 5: e42-10.1371/journal.pbio.0050042.
McKenzie FE: Why model malaria?. Parasitol Today. 2000, 16: 511-516. 10.1016/S0169-4758(00)01789-0.
Woolhouse ME, Dye C, Etard JF, Smith T, Charlwood JD, Garnett GP, Hagan P, Hii J, Ndhlovu PD, Quinnell RJ, Watts CH, Chandiwana SK, Anderson RM: Heterogeneities in the transmission of infectious agents: implications for the design of control programs. Proc Natl Acad Sci USA. 1997, 94: 338-342. 10.1073/pnas.94.1.338.
Schmid-Hempel P, Koella J: Variability and its implications for host-parasite interactions. Parasitol Today. 1994, 10: 98-102. 10.1016/0169-4758(94)90007-8.
Dye C, Hasibeder G: Population dynamics of mosquito-borne disease: effects of flies which bite some people more frequently than others. Trans R Soc Trop Med Hyg. 1986, 80: 69-77. 10.1016/0035-9203(86)90199-9.
Koella JC: On the use of mathematical models of malaria transmission. Acta Trop. 1991, 49: 1-25. 10.1016/0001-706X(91)90026-G.
White M, Griffin J, Drakeley C, Ghani A: Heterogeneity in malaria exposure and vaccine reponse: implications for the interpretation of vaccine efficacy trials. Malar J. 2010, 9: 82-10.1186/1475-2875-9-82.
Smith D, Dushoff J, Snow R, Hay S: The entomological inoculation rate and Plasmodium falciparum infection in African children. Nature. 2005, 438: 492-495. 10.1038/nature04024.
Smith T, Killeen GF, Maire N, Ross A, Molineaux L, Tediosi F, Hutton G, Utzinger J, Dietz K, Tanner M: 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 (Suppl 2): 1-10.
Grimm V, Railsback SF: Individual-based modeling and ecology. 2005, Princeton: Princeton University Press
Sabatelli L, Ghani A, Rodrigues L, Hotez P, Brooker S: Modelling heterogeneity and the impact of chemotherapy and vaccination against human hookworm. J R Soc Interface. 2008, 5: 1329-41. 10.1098/rsif.2007.1255.
Epstein J: Modelling to contain pandemics. Nature. 2009, 460: 687-10.1038/460687a.
Carneiro I, Roca-Feltrer A, Griffin J, Smith L, Tanner M, Armstrong Schellenberg J, Greenwood B, Schellenberg D: Age-patterns of malaria vary with severity, transmission intensity and seasonality in sub-Saharan Africa: a systematic review and pooled analysis. PLoS ONE. 2010, 5: e8988-10.1371/journal.pone.0008988.
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 (Suppl 2): 19-31.
Smith T, Ross A, Maire N, Rogier C, Trape JF, Molineaux L: An epidemiological model of the incidence of acute illness in Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006, 75 (Suppl 2): 56-62.
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 (Suppl 2): 63-73.
Tediosi F, Maire N, Smith T, Hutton G, Utzinger J, Ross A, Tanner M: 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 (Suppl 2): 90-103.
Penny M, Maire N, Studer A, Schapira A, Smith T: What should vaccine developers ask? Simulation of the effectiveness of malaria vaccines. PLoS ONE. 2008, 3: e3193-10.1371/journal.pone.0003193.
Maire N, Tediosi F, Ross A, Smith T: 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. 2006, 75 (Suppl 2): 111-118.
Ross A, Penny M, Maire N, Studer A, Carneiro I, Schellenberg D, Greenwood B, Tanner M, Smith T: Modelling the epidemiological impact of intermittent preventive treatment against malaria in infants. PLoS ONE. 2008, 3: e2661-10.1371/journal.pone.0002661.
Smith T, Maire N, Ross A, Penny M, Chitnis N, Schapira A, Studer A, Genton B, Lengeler C, Tediosi F, de Savigny D, Tanner M: Towards a comprehensive simulation model of malaria epidemiology and control. Parasitology. 2008, 135: 1507-1516. 10.1017/S0031182008000371.
Smith T, Charlwood J, Kihonda J, Mwankusye S, Billingsley P, Meuwissen J, Lyimo E, Takken W, Teuscher T, Tanner M: Absence of seasonal variation in malaria parasitaemia in an area of intense seasonal transmission. Acta Trop. 1993, 54: 55-72. 10.1016/0001-706X(93)90068-M.
Smith D, Dushoff J, McKenzie F: The risk of a mosquito-borne infection in a heterogeneous environment. PLoS Biology. 2004, 2: e368-10.1371/journal.pbio.0020368.
Barat L, Palmer N, Basu S, Worrall E, Hanson K, Mills A: Do malaria control interventions reach the poor? A view through the equity lens. Am J Trop Med Hyg. 2004, 71 (Suppl 2): 174-178.
Onwujekwe O, Hanson K, Uzochukwu B, Ichoku H, Ike E, Onwughalu B: Are malaria treatment expenditures catastrophic to different socio-economic and geographic groups and how do they cope with payment? A study in southeast Nigeria. Trop Med Int Health. 2010, 15: 18-25.
Stoddard S, Morrison A, Vazquez-Prokopec G, Paz Soldan V, Kochel T, Kitron U, Elder P, Scott T: The role of human movement in the transmission of vector-borne pathogens. PLoS Negl Trop Dis. 2009, 3: e481-10.1371/journal.pntd.0000481.
Bates I, Fenton C, Gruber J, Lalloo D, Medina Lara A, Squire S, Theobald S, Thomson R, Tolhurst R: Vulnerability to malaria, tuberculosis, and HIV/AIDS infection and disease. Part 1: determinants operating at individual and household level. Lancet. 2004, 4: 267-277. 10.1016/S1473-3099(04)01002-3.
Kidane G, Morrow R: Teaching mothers to provide home treatment of malaria in Tigray, Ethiopia: a randomised trial. Lancet. 2004, 356: 550-555. 10.1016/S0140-6736(00)02580-0.
Smith T, Maire N, Dietz K, Killeen G, Vounatsou P, Molineaux L, Tanner M: Relationship between the entomologic inoculation rate and the force of infection for Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006, 75 (Suppl 2): 11-18.
Ross A, Killeen G, Smith T: Relationships between host infectivity to mosquitoes and asexual parasite density in Plasmodium falciparum. Am J Trop Med Hyg. 2006, 75 (Suppl 2): 32-37.
Killeen G, Ross A, Smith T: Infectiousness of malaria-endemic human populations to vectors. Am J Trop Med Hyg. 2006, 75 (Suppl 2): 38-45.
Ross A, Smith T: The effect of malaria transmission intensity on neonatal mortality in endemic areas. Am J Trop Med Hyg. 2006, 75 (Suppl 2): 74-81.
Marsh K, Snow R: Malaria transmission and morbidity. Parassitologia. 1999, 41: 241-246.
Smith T: Estimation of heterogeneity in malaria transmission by stochastic modelling of apparent deviation from mass action kinetics. Malar J. 2007, 11: 7-12.
Bejon P, Warimwe G, Mackintosh C, Mackinnon M, Kinyanjui S, Musyoki J, Bull P, Marsh K: Analysis of immunity to febrile malaria in children that distinguishes immunity from lack of exposure. Infect Imm. 2009, 77: 1917-1923. 10.1128/IAI.01358-08.
INDEPTH Network: Population, health and survival at INDEPTH sites. 2002, Ottowa, Ontario, Canada: International Development Research Centre
We are grateful for comments from Allan Schapira, Nicolas Maire and members of the Technical Advisory Group: Michelle Gatton, Azra Ghani, Simon Hay, Steve Lindsay, Kevin Marsh and David Schellenberg. We thank Diggory Hardy and Nicolas Maire for developing the simulation software architecture and the many volunteers who donated computing power via malariacontrol.net. This study was funded by The Bill and Melinda Gates Foundation grant number 39777.01
The authors declare that they have no competing interests.
AR conceived the study, designed the simulations and drafted the manuscript. TS edited the manuscript. Both authors have read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: Summary of predicted age-prevalence and incidence curves using age and height at peak.(DOC 116 KB)
Additional file 2: Table S2: Effects of pairs of heterogeneities on predicted age-prevalence and incidence curves.(DOC 39 KB)
Authors’ original submitted files for images
About this article
Cite this article
Ross, A., Smith, T. Interpreting malaria age-prevalence and incidence curves: a simulation study of the effects of different types of heterogeneity. Malar J 9, 132 (2010). https://doi.org/10.1186/1475-2875-9-132
- Severe Malaria
- Parasite Density
- Transmission Intensity
- Parasite Prevalence