Simulation of malaria epidemiology and control in the highlands of western Kenya
Malaria Journal volume 11, Article number: 357 (2012)
Models of Plasmodium falciparum malaria epidemiology that provide realistic quantitative predictions of likely epidemiological outcomes of existing vector control strategies have the potential to assist in planning for the control and elimination of malaria. This work investigates the applicability of mathematical modelling of malaria transmission dynamics in Rachuonyo South, a district with low, unstable transmission in the highlands of western Kenya.
Individual-based stochastic simulation models of malaria in humans and a deterministic model of malaria in mosquitoes as part of the OpenMalaria platform were parameterized to create a scenario for the study area based on data from ongoing field studies and available literature. The scenario was simulated for a period of two years with a population of 10,000 individuals and validated against malaria survey data from Rachuonyo South. Simulations were repeated with multiple random seeds and an ensemble of 14 model variants to address stochasticity and model uncertainty. A one-dimensional sensitivity analysis was conducted to address parameter uncertainty.
The scenario was able to reproduce the seasonal pattern of the entomological inoculation rate (EIR) and patent infections observed in an all-age cohort of individuals sampled monthly for one year. Using an EIR estimated from serology to parameterize the scenario resulted in a closer fit to parasite prevalence than an EIR estimated using entomological methods. The scenario parameterization was most sensitive to changes in the timing and effectiveness of indoor residual spraying (IRS) and the method used to detect P. falciparum in humans. It was less sensitive than expected to changes in vector biting behaviour and climatic patterns.
The OpenMalaria model of P. falciparum transmission can be used to simulate the impact of different combinations of current and potential control interventions to help plan malaria control in this low transmission setting. In this setting and for these scenarios, results were highly sensitive to transmission, vector exophagy, exophily and susceptibility to IRS, and the detection method used for surveillance. The level of accuracy of the results will thus depend upon the precision of estimates for each. New methods for analysing and evaluating uncertainty in simulation results will enhance the usefulness of simulations for malaria control decision-making. Improved measurement tools and increased primary data collection will enhance model parameterization and epidemiological monitoring. Further research is needed on the relationship between malaria indices to identify the best way to quantify transmission in low transmission settings. Measuring EIR through mosquito collection may not be the optimal way to estimate transmission intensity in areas with low, unstable transmission.
Rationale for work
In order to make informed decisions for malaria control, programme managers require information on the optimal mix of intervention strategies tailored to specific transmission patterns of malaria[1–3]. This information is often unavailable due to the difficulty in measuring rates of malaria transmission and determining the impact of control interventions on transmission. While the efficacy of individual malaria control interventions in reducing morbidity and mortality in western Kenya has been demonstrated by field trials[4, 5], there have been fewer studies investigating the effects across a range of transmission intensities or for combinations of interventions[6, 7].
Since 2008, a number of epidemiological and entomological studies have been carried out in Rachuonyo South, Kenya, as part of the Malaria Transmission Consortium (MTC). The availability of data from these and other studies presents an opportunity for site-specific parameterization of models of malaria transmission. The results of these model simulations can be translated into evidence-based decision making for malaria control programme managers. This project applies individual-based stochastic models of malaria to MTC sites with transmission data to simulate the impact of a range of malaria control strategies.
Rachuonyo South district is situated in Nyanza province, bordering Lake Victoria in western Kenya (Figure1) and encompasses an area of 930km2. The main MTC study site is located in the south west of the district and represents a highland “fringe” area (1,400-1,600 meters above sea level). Ethnicity in Rachuonyo South is predominantly the Luo ethnic group. Residents depend upon farming and cattle and goat herding for subsistence. Homesteads are distributed broadly across a rolling landscape intersected with small streams and rivers. Total annual rainfall in this area averages 1,200 mm per year (Figure2) while average daily temperatures range from 17-27°C. The area is characterized by generally low malaria endemicity with marked seasonal and inter-annual variations in transmission[8, 9].
The main malaria vectors in the highlands were previously recorded to be Anopheles gambiae sensu stricto, Anopheles arabiensis and Anopheles funestus[10, 11]. In recent years, there is evidence that An. gambiae s.s. is disappearing from lowlands Nyanza leaving An. arabiensis as the predominant species within the An. gambiae sensu lato complex and An. funestus as the primary Plasmodium falciparum vector (Stevenson, personal communication). These changes are most likely due to intensive targeting of malaria control interventions, but climatic factors may also have played a role[12–14].
In western Kenya indoor residual spraying (IRS) campaigns were carried out in the Kericho district in the 1940s and Nandi district in the 1950s (using dichlorodiphenyltrichloroethane (DDT) and dieldrin, respectively). It is thought that malaria transmission was largely eliminated from large portions of the highlands as a result[15, 16]. While epidemics re-emerged in the 1980s[17, 18], it was not until after the year 2000 that routine, large-scale vector control interventions were introduced in these areas.
The main control methods used today in the epidemic highland areas include mass-distribution of long-lasting insecticide-treated nets (LLINs), IRS with pyrethroids, and prompt and effective treatment of malaria[19–21]. Artemisinin-based combination therapy (ACT), specifically artemether-lumefantrine (AL) was adopted as the first line treatment drug in 2006 following a decline in efficacy of sulphadoxine-pyrimethamine (SP) and amodiaquine, the previous first and second line treatments, respectively. In 2006 and 2011 Rachuonyo South was included within the Kenyan national mass distribution LLIN campaign and distribution continues through antenatal clinics, child welfare clinics, and comprehensive care clinics for people living with HIV. Since 2005 Rachuonyo South has been targeted for universal coverage of IRS once per year in advance of the main transmission season. Different formulations of pyrethroid insecticide have been used over the years with lambdacyhalothrin (ICON) used in 2009, alphacypermethrin (FENDONA) used in 2010, ICON again in 2011, and 2012 started with ICON and then switched to deltamethrin.
OpenMalaria transmission model
A team at the Swiss Tropical and Public Health Institute (Swiss TPH) and Liverpool School of Tropical Medicine (LSTM) has developed stochastic simulation models of transmission of malaria based on the simulation of infection in individuals that are able to simulate the impact (cost-effectiveness, clinical, epidemiological and entomological) of numerous intervention strategies for malaria control[22–26]. These models form part of the OpenMalaria platform that makes the considerable code base written in C++ freely available online. Users are able to carry out predictive simulations either via a downloadable stand-alone programme or via a volunteer grid computing resource and semi-automated experiment design and analysis platform capable of handling entire experiments of 10,000-100,000 scenarios.
Individual infections are simulated by stochastic series of parasite densities, which determine an individual’s morbidity and mortality risks as well as their infectiousness to vectors[22, 27]. The simulated infections are nested within simulations of individuals in human populations, and linked to a model of transmission of malaria between humans and mosquitoes and to models of interventions[22, 23, 27]. The transmission model is based on a periodically-forced difference equation model for malaria mosquitoes feeding on, infecting and getting infected from a heterogeneous population of hosts. These dynamics are calibrated by a seasonal pattern of EIR for each mosquito species assuming that in the absence of interventions EIR seasonality is fixed across years. Simulations are run for one human life span to induce an “equilibrium” level of immunity in the population. Subsequent dynamics are used to predict available malaria outcomes, such as patterns of infection in humans or patterns of disease by age and season, which can then be compared to field data.
The details of the methods to build and parameterize the transmission model used in this project have been published elsewhere[22–26] and therefore are not covered in this paper. In this paper the model components described above are employed to an ensemble of 14 model variants for malaria in humans to address stochasticity and model uncertainty. Simulations were repeated with multiple random seeds to address parameter uncertainty.
The models included in the OpenMalaria platform were initially parameterized from published data from Namawala, Tanzania[22–28]; 61 data sets were used to optimize certain parameters[22–26]. To update the parameterization for the Rachuonyo South scenario, data collected as part of the MTC project in the study area was the first choice to use for the model parameters. A description of these studies and how they were used to parameterize the model can be found in Additional file1.
MTC field studies
A number of field studies were carried out in Kisii and Rachuonyo South districts between 2009 and 2011 with the goal of establishing an evidence base to help malaria control programme managers monitor malaria transmission and implement and adjust malaria control interventions. Data from these studies are currently being analysed and will be described in detail in forthcoming publications. For the purposes of the modelling work described in this paper, the datasets used are described in Table1.
Where data were not available from MTC surveys, parameter inputs were identified via a literature review of publications using the PubMed electronic database using the key words “Kenya, Nyanza, Rachuonyo, western Kenya, malaria, Plasmodium falciparum, transmission, antimalarials, artemether- lumfantrine, insecticide residual spraying, insecticide-treated nets, larviciding, intermittent preventive treatment, modelling, malaria incidence, treatment seeking, mosquito resting duration, extrinsic incubation period, Anopheles.” An internet review was also conducted on the websites for the Kenyan Ministry of Health, Division of Malaria Control, the National Bureaus of Statistics, and the National Demographic Health Surveys. The sources were prioritized in the following strata in order of precedence: study area districts MTC data collection, study area districts existing literature, study area provincial data, national level data, existing model parameterization. Where more than one data source was found within any one stratum the study with the closest site characteristics or, where applicable, date of data collection closest to that of the MTC studies was used.
To determine the annual average EIR, the transmission parameter in the model, seroconversion rates using the MSP-1 antigen were estimated from the July 2009 cross-sectional survey as described in Drakeley et al. 2005 and derived EIR equivalents were calculated as described in Corran et al. 2007. The average monthly EIR values used to calibrate the seasonal pattern of transmission in the scenario were calculated by separating the annual average EIR from existing literature for a neighboring district into the monthly proportion of rainfall in Rachuonyo South recorded by the Kogalo weather station so that the peak malaria transmission month corresponded to one month later than the peak rainfall month (Figure2). Because the annual average EIR is based on serology, the model incorporates the overall temperature and humidity effects but excludes the seasonality of these effects.
In practice, many of the entomological and health system parameters were based on data from elsewhere used in other modelling exercises[26–32] as they are thought to be fairly standard across anopheline species and anti-malarials. However, because several entomological parameters are sensitive to temperature, particularly the extrinsic incubation period (EIP) and mosquito resting duration[33, 34], these values were adjusted for each study area based on the average annual temperature collected by the Kogalo weather station. Also, the latest data from the study site challenges the assumption that vectors are normally predominantly endophilic and endophagic. For the purposes of this experiment, emphasis was placed on overall vector biting behaviour rather than simulating individual species. This was due to the design of the entomological field studies for which results were available at the time of model parameterization that focused on indoor/outdoor species composition and trap evaluation rather than the biting behaviour within individual species. The efficacy of LLINs and IRS were adjusted to affect the indoor mosquitoes but not the outdoor mosquitoes and the proportion of bites on a human compared to other mammals was reduced for the outdoor mosquitoes.
The monitoring measures serving as the outputs simulated by the model were chosen based on the indicators of malaria transmission measured by the field studies described above.
Before the main simulation, the scenario was run for one human life span to ensure each simulated individual acquired the expected natural immunity for his or her age. The fitting of the dynamic EIR in the transmission model to the pre-intervention calibration EIR was done during the last five years of the life span simulation. A subpopulation was considered as a cohort and received mass drug administration (MDA) at the beginning of the main simulation, to “mimic” the MTC cohort study conditions, where participants were given a course of the first-line malaria treatment upon enrollment into the study to clear any existing malaria parasites. Finally, the effect of interventions on epidemiological outcomes of malaria in the full population of the study area was simulated for two years.
Validation and sensitivity analysis
The project addresses uncertainty on three levels: stochasticity, model uncertainty, and parameter uncertainty. Each simulation was repeated by the OpenMalaria simulator on an ensemble of 14 model variants using ten random seeds in order to address model uncertainty and stochasticity. Results in the form of graphs from the ensemble of model variants were visually analysed and compared to observed data from the study areas using Stata (version 11; College Station, TX, USA). Further analysis of the scenario simulation and observed data for the selected impact measures was conducted using Stata. The proportion of simulation results falling within the 95% confidence intervals of the observed cohort data was measured in order to assess goodness of fit.
A sensitivity analysis to address parameter uncertainty was driven by results of the visual comparison of stochasticity. Elements of the model central to the epidemiology and control of malaria in this particular study area were identified based on whether there was uncertainty about parameter estimates and their potential impact on the composition and behaviour of vectors, effectiveness of interventions, and population-level monitoring. These included effectiveness of IRS, indoor versus outdoor biting behaviour of the vectors, the detection limit of the survey method used for malaria in humans, annual average EIR, and climate and weather patterns affecting vector biology and parasite development in the vector. Parameters were altered one at a time and results analysed by comparing the simulated number of cases per person per year for each scenario to those of the baseline parameterization.
Model design and baseline scenario parameterization
The following parameter estimates are based on currently unpublished data from the MTC studies described above in Table1. The key entomological feature of the scenario involves one primary malaria vector that bites and rests outdoors 62% of the time. At the time of enrollment into the cohort, 30.3% of the cohort population slept under a net the previous night and 69.3% of survey households received IRS. The IRS deployment schedule happened yearly over a period of two months as described above in the Background section. The annual mean temperature in Rachuonyo South from 2009 to 2010 was recorded as 20.3 degrees Celsius, setting the estimate of the extrinsic incubation period of An. gambiae at 14 days and the resting duration 3 days. Malaria transmission is highly variable following two distinct rainy seasons. The EIR is unstable with a last recorded value from an entomological survey of 0.4 infectious bites per person per year. This study was conducted in neighboring Kisii district before LLIN and IRS scale-up in 2006. More recent results from the July 2009 MTC cross sectional study estimate an EIR of 1.5 infectious bites per person per year based on serological data (Table2).
For OpenMalaria to simulate dynamics of the study population, code was included in the scenario to select a cohort representing 15% of the total population over one year old matching the cohort enrollment criteria, all of whom received a course of anti-malarials at the start of the survey period. The validation of the model uses the model outputs from only this cohort, while the remaining simulations represent the larger study area population of 10,000 individuals. The details of the values used to parameterize the model along with their sources can be found in Additional Files2,3,4,5,6,7.
Simulation and validation
OpenMalaria is able to simulate the seasonality and level of the EIR for the Rachuonyo South scenario with greater stochasticity in the peak months and in the scenario with observed interventions (Figure3). Simulations show prevalence between 5.58% and 10.81% in Rachuonyo South’s peak transmission month and between 2.99% and 6.04% in the lowest transmission month (Figure4).
The Figures5a and b compare the simulation of P. falciparum prevalence in the population with observed data from the MTC cohort study conducted from June 2009 – June 2010 in Rachuonyo South District as detected by a rapid diagnostic test (RDT) using EIR values derived from entomological studies (0.4 infectious bites per person per year, seasonality from neighboring district, Table2) versus serology (1.5 infectious bites per person per year, seasonality from study site weather station). The prevalence was especially high in June of 2010, possibly due to a combination of more rainfall than normal during the rainy season and rollout of IRS at a later month compared to the previous year.
While the model is able to predict the level of prevalence in both scenarios, using an EIR from serology and seasonality from weather station data represents a visually better fit with both level of overall and seasonal prevalence compared to using an EIR and seasonality from entomology. Using the proportion of simulation runs falling within the 95% confidence intervals of the observed cohort data as a benchmark for comparing simulation results, the final scenario was able to improve both the number of months (six months with more than 30% of simulations runs predicted compared to three months, n=12) and the proportion of total simulation runs (29.9% vs. 14.6%, n=1,680).
Indoor residual spraying
The two main malaria control measures in the study area are distribution of LLINs and IRS. While net use is assumed constant over the time frame of the simulation, IRS is a timed intervention that occurred between April and May of 2009 and June and July of 2010 (Figure2). To simulate the impact of IRS effectiveness at killing and deterring vectors and the rate at which the insecticide decays on model predictions, scenarios were created to simulate very high and very low IRS effectiveness (Table3).
Compared to the baseline, increasing the duration and effectiveness of IRS had the effect of greatly reducing the simulated number of patent infections (Figure6b). While prevalence is greatly reduced, transmission is never completely interrupted even in the scenario simulating highly effective IRS.
To study the effects of changes in vector diversity and biting behaviour, different scenarios of proportion of indoor vs. outdoor biting are considered. The baseline scenario assumes one primary vector species which bites outdoors 64% of the time and indoors 36% of the time. The experiment includes one scenario with increased exophagy with 74% of transmission occurring outdoors and 26% of transmission occurring indoors and a second scenario with transmission is split equally indoors and outdoors. This is modeled by reducing the effectiveness of vector control interventions.
The scenario in which the biting behaviour of a single vector species is altered and a greater proportion of EIR is due to indoor biting (increased from 36% to 50%) shows a reduction in prevalence (Figure7c). This is because the indoor mosquitoes would be affected by the IRS campaigns conducted in April – May of the first year and June – July of the second year. The scenario with a greater proportion of transmission from outdoor biting shows a similar level of transmission during the low season but greater amplitude in peak months (Figure7b).
Survey detection limit
To address the model’s sensitivity to the ability of a given test to detect a P. falciparum infection, an experiment was created to mimic the detection limits of polymerase chain reaction (PCR), skilled microscopy, and a low-quality diagnostic such as a poor-quality RDT or unskilled microscopy (Table4). The number of simulated infections decreases with higher detection limits, as does the stochasticity of the predictions (Figure8). This indicates a population that has a considerable proportion of infections occurring characterized by low parasitaemia.
EIR and climatic patterns
In order to account for differences in collection and calculation method as well as micro-variations in EIR within the study area, an experiment was conducted with varying levels of the annual average EIR while keeping the seasonal pattern the same over the baseline. This includes a scenario with a low EIR value that was measured in an neighbouring district a slightly higher altitude before large-scale control programmes were implemented in 2006, two scenarios with medium EIR (one equal to double the recorded value and one equal to the recorded value in the neighboring lowland districts), and a larger EIR. OpenMalaria is able to simulate the scenarios with EIRs of 7 and 20 with less stochasticity than the scenarios with smaller EIRs.
To examine model sensitivity to changes in the entomological parameters that could occur as a result of different climate patterns an experiment was created to simulate decreased rainfall, increased rainfall, decreased temperature, increased temperature, and two long rainy periods instead of the long and short rains the study area currently experiences. Compared to the baseline, simulating increased rainfall in the same seasonal pattern did not have as great an effect on number of patent infections as did the scenario which increased the short rainy season to match the longer rains. Simulating temperature changes by altering the extrinsic incubation period and resting duration did not have a visible impact on the predicted number of patent infections (Figure9).
Figure9 demonstrates the overall results of the one-way sensitivity analysis in relation to the baseline scenario for Rachuonyo South.
Validation and model analysis
The scenario parameterized for Rachuonyo South district is able to replicate the overall level of prevalence in the given population for the majority of the months out of the year. However, the timing of the peak transmission month is delayed and the depth of the trough in November of the first year of simulations is not captured. Thus, the number of runs simulating the number of patent infections falling within the 95% confidence intervals (CIs) of the observed number of patent infections is lower than optimal. This could be due to several factors, such as inter-annual variation in transmission in the study area.
A main challenge for transmission models is calibration and validation with data from the field. The OpenMalaria transmission model is calibrated primarily by the intensity of malaria transmission, or EIR, for each vector. There are several methods to measure EIR in the field, and the method used varies by location depending on the implementer of the study. Usually the types of surveys necessary to quantify transmission are not done on a regular basis, and in low transmission settings where mosquito densities are low, the longitudinal studies required to estimate EIR are intensive and inherently expensive. Entomological studies with the aim of identifying sporozoite-positive mosquitoes, while important for monitoring vector biting behavior, are not suitable for developing a seasonality pattern for a given total EIR in this area of low, unstable transmission. Perhaps monthly or even weekly studies measuring mosquito density and changes in vector biting behaviour over one or multiple years would be a way of determining the seasonality pattern of transmission.
Despite a clearer picture of overall annual transmission, serology is not able to characterize seasonality of transmission, and without a baseline is unable to give an indication of pre- vs. post- intervention exposure. Serology combined with a seasonality pattern from rainfall data offered a more accurate picture than entomological data alone. While this information can be approximated from weather station data and from remote sensing in areas lacking a weather station, a challenge is relating the amount and seasonality of transmission to the amount of rainfall as their relationship is not linear.
The method of evaluation used in this study was to analyse the number of simulation runs which fall within the 95% CIs of the observed data. There is not yet any consensus on how to evaluate uncertainty and goodness-of-fit for model ensembles. The merits of different methods have been discussed for models used in meteorology, climate change and macroeconomics, but questions remain on whether model averaging is appropriate and how to quantify an acceptable level of stochasticity for basing programmatic decisions on model predictions. A consensus should be achieved on these criteria if quantitative projections from such models are to become an integral part of the range of decision-making tools for malaria control.
Implications of the sensitivity analysis
The sensitivity analysis highlights the robustness of the OpenMalaria transmission model for simulating a range of entomological and epidemiological scenarios. The majority of the simulation results for extreme scenarios of the entomological and biological components of the model remain similar to the simulation results for the baseline scenario, suggesting that small changes in these parameters are unlikely to have a large impact on prevalence, while changes in EIR and effectiveness of IRS have a greater impact on the estimated prevalence in the study area (Figure9).
Pyrethroid resistance has already been documented in western Kenya and elsewhere and much depends on the effectiveness of these insecticides[41–43]. In the study area there are not yet reports of pyrethroid knockdown resistance (kdr) mutations due to the lack of presence of An. gambiae s.s., but there may be other resistance mechanisms present, for example metabolic resistance, given the high numbers of An. funestus in the study area. The results of the sensitivity analysis suggest that malaria incidence and prevalence are likely to increase as this resistance continues to rise. As noted in the background section, the Kenya DOMC has alternated deployment of different types of pyrethroid insecticides for different years. While this could potentially have the effect of discouraging resistance to any one formulation, until new insecticides are developed the continued use of only pyrethroids has the potential to encourage resistance.
Initial results of entomological surveys (Cooke, personal communication) show evidence of a shift in the relative importance of outdoor biting compared to what has been observed in the neighbouring highland district in the past[10, 44]. It is unclear whether this is a behavioral change in response to high LLIN and IRS coverage or whether there have been alterations in overall species composition. For Rachuonyo South there are no baseline data to compare this to. Evidence from lowland districts within Nyanza indicate that both composition and biting behaviour of the malaria vectors has changed over the past five years, coinciding with a substantial scaling up of vector control interventions. Entomological surveys conducted in 2009 – 2011(Stevenson, personal communication) show that An. arabiensis is now seen more frequently inside and outside dwellings than An. gambiae s.s., the previously-documented major vector in Kisii district[10, 44]. Preliminary data from the study sites also indicate that An. funestus or other species may be playing an ever increasingly important role on malaria transmission in the area.
The observed data more closely resemble the scenario with indoor/outdoor biting profile based on 2009 – 2011 data, which supports the hypothesis of a greater proportion of outdoor biting. If this is the case, there is a limit to the effectiveness of the current vector control interventions in Rachuonyo South (IRS and LLINs) at controlling P. falciparum because they target the shrinking proportion of the infective bites occurring indoors. While these interventions will still offer an important level of deterrency, interventions that have a killing effect on exophagic mosquitoes may be an appropriate addition to existing indoor interventions. Larviciding, area repellents, and even interventions targeting the human-stage of the parasite could also be taken into consideration to complement existing methods. Implementation of a number of these methods is currently being piloted in Rachuonyo South (Bousema, personal communication).
Survey detection limit
The outcome simulated in this scenario is number of patent infections as measured by a Paracheck® rapid test kit manufactured by Orchid Biomedical Systems. The 2010 WHO malaria case management guidelines recommend treatment after parasite-based diagnosis. Quality assurance measures for these tests are based on their ability to detect either 100 or 200 parasites per microliter, not because of the limitations of the RDT technology but rather because of limited accuracy and error of expert microscopy, the “gold” standard in malaria diagnosis in the absence of PCR[47, 48]. In addition, there is evidence for changes in the accuracy of diagnosis by RDTs in the East African highlands both over time and across age groups.
The implication of the sensitivity of the model to a change in survey detection limit is that if RDTs used in surveys perform poorly, whether the result of low quality manufacturing or improper storage conditions or use, according to simulation results up to 50% of infected individuals would be misclassified. When put in a broader public health context, there are a number of scenarios applicable to the study area when decision-making can be affected by detection limit. These range from a health worker deciding to administer an anti-malarial drug following malaria diagnosis in an individual, to country-wide planning in the public sector health system for estimating quantities of antimalarial drugs required for a given year, to deciding the appropriate time to change the vector control strategy if the measure of transmission is based on an estimate of prevalence.
When approaching a situation where transmission is interrupted, attention must be paid to the type of screening strategy (active vs. passive case detection) and screening method used to detect the last remaining parasitaemia in the population. In these cases the higher presence of asymptomatic, sub-patent infections representing the infectious reservoir of parasites in the population indicates the PCR method would be preferable over a less sensitive method. While molecular diagnostic tools such as PCR and loop-mediated isothermal amplification (LAMP) are both able to detect infections at a much lower parasite density than microscopy or RDTs and may be appropriate in study settings, studies show these methods are not currently suitable for routine diagnosis at a community level[50, 51]. However, even the most sensitive PCR diagnostic does not detect all infections in a population. If a large proportion of infections occur at a high parasite density the detection limit of the diagnostic would not be as important a consideration. This sensitivity analysis shows that observed prevalence depends on the method used for detection, a point relevant for study design and modelling alike.
EIR and climatic patterns
The sensitivity analysis results show that an increase in EIR corresponds to an increase in cases of malaria. The OpenMalaria transmission model is dependent on the length of the gonotrophic cycle of the vector, which is in turn affected by environmental changes. The mosquito resting duration and EIP both decrease as the ambient temperature decreases[33, 52, 53]. If the EIP duration decreases, a vector infected with P. falciparum becomes infectious more quickly. A shorter gonotrophic cycle means both increased biting frequency and increased daily mortality of the vector. The highlands of western Kenya have variable seasonal temperature and rainfall changes; for example, in the late 1990s the study area experienced a resurgence of malaria not seen for decades. Simulation results indicate that changes in temperature resulting in a change in EIP or resting duration and changes in the overall volume of rainfall resulting in a slight change in EIR are not likely to effect the impact of IRS deployment or result in a shift in P. falciparum prevalence in the population. In addition, there is preliminary evidence that in the study area increased relative humidity is associated with an increased number of anophelines (Cooke, personal communication). However, even taking into account the caveats for the relationship between malaria transmission and rainfall, changing the pattern of transmission to simulate the effect of an extension of the historically short rainy season to match the rainfall profile of the longer rainy season could result in greater amplitude of incidence in the peak months.
Many parameters in the model remain from the initial Tanzanian model parameterization, for example the parameters for the mosquito feeding cycle (Additional File2) and treatment-seeking behaviour (Additional File3), because there are not yet site-specific studies with this focus. While ample entomological data were collected in the study area, there was less available information on treatment-seeking behaviour and its consequences outside the public sector.
Coartem® was given to all MTC cohort study participants to clear any prevalent P. falciparum parasitaemia making it possible to measure malaria incidence at each follow-up. The study excluded pregnant women from the cohort due to the limited data on use in pregnancy and contraindication in the 1st trimester pregnancy of artemether-lumefantrine, the active ingredients of Coartem®. Infection with P. falciparum during pregnancy has been shown to be associated with increased parasitaemia of the mother due to a weakened immune system as well as an increased likelihood of manifestation of clinical disease in addition to adverse affects on the fetus and newborn. Although this is unlikely to have a major effect on transmission in the population as a whole, the patent infections and uncomplicated episodes in the age groups for women of childbearing age could be underestimated.
Since the OpenMalaria transmission model was developed to examine the effect of moderate to high transmission, it does not include a mechanism to account for inter-annual variation in EIR as driven by climatic factors. Thus, every year is treated as the same, which is not the case in the study area. In the western Kenyan highlands sharp increases in incidence occur every few years and are likely to be driven by climate variability; the higher than usual transmission in the cohort following heavy rains during the time of the survey provides an example of such a sharp increase. As a result of the validation using one year’s data in this area of substantial year-to-year variation, firm conclusions are unable to be drawn about the longer-term seasonal transmission in the population.
Nyanza province has the highest prevalence of HIV in Kenya at 15.1% of the population[56, 57]. HIV infection increases an individual’s susceptibility to malaria infection and severity of clinical outcomes and decreases immunity[58, 59]. A limitation of the transmission model is that it does not account for the interaction between malaria and HIV.
The models analysed here do not explicitly take spatial associations into account. Variation in proximity to breeding sites could be a factor driving the difference in epidemic profile of the study area. The parameterizations used in this study do not take into account the rate of imported cases from lowland areas, as there is frequent travel between the highland Kisii and Nyamira districts and the lowland areas of Rachuonyo North, Nyando and Kisumu districts. Heterogeneity in availability to vectors and imported cases should be taken into account in future simulations of the study area.
Individual-based stochastic simulations of malaria can be used as a tool to assist decision making for malaria control programmes by testing assumptions about the seasonal pattern of transmission, vector diversity and behavior, and intervention effectiveness in district-level settings. Efforts should be made to ensure models aiding in the understanding of site-specific transmission dynamics are more accessible to programme managers. The sensitivity analysis shows that in order to simulate malaria in the Rachuonyo South highlands, attention must be paid to vector biting behaviour, their susceptibility to IRS, and the detection method used for human surveys. These features will have an impact on predicting the impact of interventions in areas with low and/or variable P. falciparum transmission. The sensitivity analysis also demonstrates the accuracy of the model and can lend confidence to end users of these results in informing control options. New methods and tools for analysing and evaluating simulation results will enhance the usefulness of simulations for malaria control decision-making. Measuring EIR through mosquito collection may not be the optimal way to define transmission in areas with low, unstable transmission. Further research into the relationship between different measures of malaria is needed to better quantify transmission in low transmission settings.
Artemisinin Combination Therapy
Extrinsic Incubation Period
Entomological Inoculation Rate
Indoor Residual Spraying
Kenya Medical Research Institute/Centers for Disease Control and Prevention
Loop-Mediated Isothermal Amplification
Long-Lasting Insecticide-Treated Net
Mass Drug Administration
Malaria Transmission Consortium
Polymerase Chain Reaction
Rapid Diagnostic Test
- Swiss TPH:
Swiss Tropical and Public Health Institute.
Hay SI, Rogers DJ, Toomer JF, Snow RW: Annual Plasmodium falciparum entomological inoculation rates (EIR) across Africa: literature survey, Internet access and review. Trans R Soc Trop Med Hyg. 2000, 94: 113-127. 10.1016/S0035-9203(00)90246-3.
The malERA Consultative Group on Vector Control: A research agenda for malaria eradication: vector control. PLoS Med. 2011, 8: e1000401-
Chitnis N, Schapira A, Smith T, Steketee R: Comparing the effectiveness of malaria vector-control interventions through a mathematical model. Am J Trop Med Hyg. 2010, 83: 230-240. 10.4269/ajtmh.2010.09-0179.
Hawley WA, ter Kuile FO, Steketee RS, Nahlen BL, Terlouw DJ, Gimnig JE, Shi YP, Vulule JM, Alaii JA, Hightower AW, Kolczak MS, Phillips-Howard PA: Implications of the western Kenya permethrin-treated bed net study for policy, program implementation, and future research. Am J Trop Med Hyg. 2003, 68: 168-173.
Pluess B, Tanser FC, Lengeler C, Sharp BL: Indoor residual spraying for preventing malaria. Cochrane Database Syst Rev. 2010, 4: CD006657-
Guyatt HL, Corlett SK, Robinson TP, Ochola SA, Snow RW: Malaria prevention in highland Kenya: indoor residual house-spraying vs. insecticide-treated bednets. Trop Med Int Health. 2002, 7: 298-303. 10.1046/j.1365-3156.2002.00874.x.
Hamel MJ, Otieno P, Bayoh N, Kariuki S, Were V, Marwanga D, Laserson KF, Williamson J, Slutsker L, Gimnig J: The combination of indoor residual spraying and insecticide-treated nets provides added protection against malaria compared with insecticide-treated nets alone. Am J Trop Med Hyg. 2011, 85: 1080-1086. 10.4269/ajtmh.2011.10-0684.
Mabaso ML, Craig M, Ross A, Smith T: Environmental predictors of the seasonality of malaria transmission in Africa: the challenge. Am J Trop Med Hyg. 2007, 76: 33-38.
Kristan M, Abeku TA, Beard J, Okia M, Rapuoda B, Sang J, Cox J: Variations in entomological indices in relation to weather patterns and malaria incidence in East African highlands: implications for epidemic prevention and control. Malar J. 2008, 7: 231-10.1186/1475-2875-7-231.
Ndenga B, Githeko A, Omukunda E, Munyekenye G, Atieli H, Wamai P, Mbogo C, Minakawa N, Zhou G, Yan G: Population dynamics of malaria vectors in western Kenya highlands. J Med Entomol. 2006, 43: 200-206. 10.1603/0022-2585(2006)043[0200:PDOMVI]2.0.CO;2.
Okara RM, Sinka ME, Minakawa N, Mbogo CM, Hay SI, Snow RW: Distribution of the main malaria vectors in Kenya. Malar J. 2010, 9: 69-10.1186/1475-2875-9-69.
Bayoh MN, Mathias DK, Odiere MR, Mutuku FM, Kamau L, Gimnig JE, Vulule JM, Hawley WA, Hamel MJ, Walker ED: Anopheles gambiae: historical population decline associated with regional distribution of insecticide-treated bed nets in western Nyanza Province. Kenya. Malar J. 2010, 9: 62-10.1186/1475-2875-9-62.
Minakawa N, Sonye G, Mogi M, Githeko A, Yan G: The effects of climatic factors on the distribution and abundance of malaria vectors in Kenya. J Med Entomol. 2002, 39: 833-841. 10.1603/0022-2585-39.6.833.
Zhou G, Githeko AK, Minakawa N, Yan G: Community-wide benefits of targeted indoor residual spray for malaria control in the western Kenya highland. Malar J. 2010, 9: 67-10.1186/1475-2875-9-67.
White GB: The Anopheles gambiae complex and malaria transmission around Kisumu, Kenya. Trans R Soc Trop Med Hyg. 1972, 66: 572-581. 10.1016/0035-9203(72)90302-1.
John CC, Riedesel MA, Magak NG, Lindblade KA, Menge DM, Hodges JS, Vulule JM, Akhwale W: Possible interruption of malaria transmission, highland Kenya, 2007–2008. Emerg Infect Dis. 2009, 15: 1917-1924.
Malakooti MA, Biomndo K, Shanks GD: Reemergence of epidemic malaria in the highlands of western Kenya. Emerg Infect Dis. 1998, 4: 671-676. 10.3201/eid0404.980422.
Shanks GD, Biomndo K, Hay SI, Snow RW: Changing patterns of clinical malaria since 1965 among a tea estate population located in the Kenyan highlands. Trans R Soc Trop Med Hyg. 2000, 94: 253-255. 10.1016/S0035-9203(00)90310-9.
Division of Malaria Control [Ministry of Public Health and Sanitation], Ministry of Medical Services: National Guidelines for Diagnosis, Treatment and Prevention of Malaria in Kenya. 2010, Division of Malaria Control [Ministry of Public Health and Sanitation], Nairobi, Third
Fillinger U, Ndenga B, Githeko A, Lindsay SW: Integrated malaria vector control with microbial larvicides and insecticide-treated nets in western Kenya: a controlled trial. Bull World Health Organ. 2009, 87: 655-665. 10.2471/BLT.08.055632.
Division of Malaria Control [Ministry of Public Health and Sanitation], Kenya National Bureau of Statistics, and ICF Macro: 2010 Kenya Malaria Indicator Survey. 2011, DOMC, KNBS and ICF Macro, Nairobi, Kenya
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: 1-10.
Chitnis N, Smith T, Steketee R: A mathematical model for the dynamics of malaria in mosquitoes feeding on a heterogeneous host population. J Biol Dynamics. 2008, 2: 259-285. 10.1080/17513750701769857.
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, Ross A, Maire N, Chitnis N, Studer A, Hardy D, Brooks A, Penny M, Tanner M: Ensemble modeling of the likely public health impact of a pre-erythrocytic malaria vaccine. PLoS Med. 2012, 9: e1001157-10.1371/journal.pmed.1001157.
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-1124. 10.1007/s11538-011-9710-0.
Owusu-Agyei S, Smith T, Beck HP, Amenga-Etego L, Felger I: Molecular epidemiology of Plasmodium falciparum infections among asymptomatic inhabitants of a holoendemic malarious area in northern Ghana. Trop Med Int Health. 2002, 7: 421-428. 10.1046/j.1365-3156.2002.00881.x.
Drakeley CJ, Corran PH, Coleman PG, Tongren JE, McDonald SL, Carneiro I, Malima R, Lusingu J, Manjurano A, Nkya WM, Lemnge MM, Cox J, Reyburn H, Riley EM: Estimating medium- and long-term trends in malaria transmission by using serological markers of malaria exposure. Proc Natl Acad Sci USA. 2005, 102: 5108-5113. 10.1073/pnas.0408725102.
Corran P, Coleman P, Riley E, Drakeley C: Serology: a robust indicator of malaria transmission intensity?. Trends Parasitol. 2007, 23: 575-582. 10.1016/j.pt.2007.08.023.
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.
Briet OJ, Hardy D, Smith TA: Importance of factors determining the effective lifetime of a mass, long-lasting, insecticidal net distribution: a sensitivity analysis. Malar J. 2012, 11: 20-10.1186/1475-2875-11-20.
Minakawa N, Omukunda E, Zhou G, Githeko A, Yan G: Malaria vector productivity in relation to the highland environment in Kenya. Am J Trop Med Hyg. 2006, 75: 448-453.
Macdonald G: Malariology: A comprehensive survey of all aspects of this group of diseases from a global standpoint. 1957, W.B. Saunders Company, Philadelphia
Stevenson J, StLauernt B, Lobo N, Cooke M, Kahindi S, Oriango R, Harbach R, Cox J, Drakeley C: Novel vectors of malaria parasite in the western highlands of Kenya [letter]. Emerg Infect Dis. 2012, 18: 1547-1549. 10.3201/eid1809.120283.
Shililu JI, Maier WA, Seitz HM, Orago AS: Seasonal density, sporozoite rates and entomological inoculation rates of Anopheles gambiae and Anopheles funestus in a high-altitude sugarcane growing zone in Western Kenya. Trop Med Int Health. 1998, 3: 706-710. 10.1046/j.1365-3156.1998.00282.x.
Kelly-Hope LA, McKenzie FE: The multiplicity of malaria transmission: a review of entomological inoculation rate measurements and methods across sub-Saharan Africa. Malar J. 2009, 8: 19-10.1186/1475-2875-8-19.
Paaijmans KP, Wandago MO, Githeko AK, Takken W: Unexpected high losses of Anopheles gambiae larvae due to rainfall. PLoS One. 2007, 2: e1146-10.1371/journal.pone.0001146.
Knutti R: The end of model democracy? An editorial comment. Climatic Change. 2010, 102: 395-404. 10.1007/s10584-010-9800-2.
Stainforth DA, Allen MR, Tredger ER, Smith LA: Confidence, uncertainty and decision-support relevance in climate predictions. Philos Transact A Math Phys Eng Sci. 2007, 365: 2145-2161. 10.1098/rsta.2007.2074.
Ranson H, Abdallah H, Badolo A, Guelbeogo WM, Kerah-Hinzoumbe C, Yangalbe-Kalnone E, Sagnon N, Simard F, Coetzee M: Insecticide resistance in Anopheles gambiae: data from the first year of a multi-country study highlight the extent of the problem. Malar J. 2009, 8: 299-10.1186/1475-2875-8-299.
Ombok MO G, Bayoh N, Vulule J, Gimnig J, Walker E: Entomological monitoring of the indoor residual spraying (IRS) program in western Kenya. Kenya National Malaria Forum: 10-11 October 2011. 2011, Center for Disease Control and Prevention, Nairobi, Kenya
Kawada H, Dida GO, Ohashi K, Komagata O, Kasai S, Tomita T, Sonye G, Maekawa Y, Mwatele C, Njenga SM, Mwandawiro C, Minakawa N, Takagi M: Multimodal pyrethroid resistance in malaria vectors, Anopheles gambiae s.s., Anopheles arabiensis, and Anopheles funestus s.s. in western Kenya. PLoS One. 2011, 6: e22574-10.1371/journal.pone.0022574.
Githeko AK, Adungo NI, Karanja DM, Hawley WA, Vulule JM, Seroney IK, Ofulla AV, Atieli FK, Ondijo SO, Genga IO, Odada PK, Situbi PA, Oloo JA: Some observations on the biting behavior of Anopheles gambiae s.s., Anopheles arabiensis, and Anopheles funestus and their implications for malaria control. Exp Parasitol. 1996, 82: 306-315. 10.1006/expr.1996.0038.
Killeen GF, Chitnis N, Moore SJ, Okumu FO: Target product profile choices for intra-domiciliary malaria vector control pesticide products: repel or kill?. Malar J. 2011, 10: 207-10.1186/1475-2875-10-207.
WHO: Guidelines for the Treatment of Malaria. 2010, World Health Organization, Geneva, Second
WHO: The role of laboratory diagnosis to support malaria disease management: focus on the use of RDTs in areas of high transmission. 2004, WHO Technical Consultation, Geneva
WHO: Malaria Rapid Diagnostic Test Performance: Results of WHO product testing of malaria RDTs: Round 2 (2009). 2009, World Health Organization, Geneva
Abeku TA, Kristan M, Jones C, Beard J, Mueller DH, Okia M, Rapuoda B, Greenwood B, Cox J: Determinants of the accuracy of rapid diagnostic tests in malaria case management: evidence from low and moderate transmission settings in the East African highlands. Malar J. 2008, 7: 202-10.1186/1475-2875-7-202.
Han ET, Watanabe R, Sattabongkot J, Khuntirat B, Sirichaisinthop J, Iriko H, Jin L, Takeo S, Tsuboi T: Detection of four Plasmodium species by genus- and species-specific loop-mediated isothermal amplification for clinical diagnosis. J Clin Microbiol. 2007, 45: 2521-2528. 10.1128/JCM.02117-06.
Poschl B, Waneesorn J, Thekisoe O, Chutipongvivate S, Karanis P: Comparative diagnosis of malaria infections by microscopy, nested PCR, and LAMP in northern Thailand. Am J Trop Med Hyg. 2010, 83: 56-60. 10.4269/ajtmh.2010.09-0630.
Macdonald G: The Epidemiology and Control of Malaria. 1957, Oxford University Press, London
Rogers DJ, Randolph SE: Global Mapping of Infectious Diseases: Methods, Examples and Emerging Applications. Adv Parasitol. Volume 62. Edited by: Hay SI, Graham A, Rogers D. 2006, Elsevier, San Diego, 345-381.
Githeko A, Ndegwa W: Predicting malaria epidemics in the Kenyan highlands using climate data: a tool for decision makers. Global Change & Human Health. 2001, 2: 10-
Desai M, ter Kuile FO, Nosten F, McGready R, Asamoa K, Brabin B, Newman RD: Epidemiology and burden of malaria in pregnancy. Lancet Infect Dis. 2007, 7: 93-104. 10.1016/S1473-3099(07)70021-X.
Montana L, Neuman M, Mishra V, Hong R: Spatial Modeling of HIV Prevalence in Cameroon, Kenya, and Tanzania. 2005
Montana L, Neuman M, Mishra V: Spatial modeling of HIV prevalence in Kenya. Demographic and Health Research. 2007
Grimwade K, French N, Mbatha DD, Zungu DD, Dedicoat M, Gilks CF: HIV infection as a cofactor for severe falciparum malaria in adults living in a region of unstable malaria transmission in South Africa. AIDS. 2004, 18: 547-554. 10.1097/00002030-200402200-00023.
Whitworth J, Morgan D, Quigley M, Smith A, Mayanja B, Eotu H, Omoding N, Okongo M, Malamba S, Ojwiya A: Effect of HIV-1 and increasing immunosuppression on malaria parasitaemia and clinical episodes in adults in rural Uganda: a cohort study. Lancet. 2000, 356: 1051-1056. 10.1016/S0140-6736(00)02727-6.
The authors thank the families who participated in the MTC studies, the members and external advisors of the MTC for valuable discussions, and the staff of KEMRI/CDC in Kisumu, Kenya, in particular Nabie Bayoh and John Gimnig, for providing data for parameterization of the model. We acknowledge our colleagues in the Swiss TPH Health Systems Research and Dynamical Modeling Unit, particularly Michael Tarantino, for help with the models and simulations, the MalariaControl.net participants for volunteering computing resources to run OpenMalaria, and two anonymous reviewers for their helpful comments. This project is a part of the Malaria Transmission Consortium and has been funded in whole or in part with funds from the Bill and Melinda Gates Foundation under Grant No. 45114. This article has been approved by the Director of the Kenya Medical Research Institute.
The authors declare that they have no competing interests.
EMS designed the experiments, performed the literature review for model parameterization, analysed results and drafted the manuscript. JS participated in parameterization of the model and designed, supervised and conducted the MTC field studies. MC provided field implementation and sample analysis for the MTC entomological field studies. CO and EM supervised and coordinated the field collection of samples. GO was responsible for data management. DH programmed the simulation software. CD provided serological analysis for the MTC field studies. TAS and JC conceived of and designed the study. NC participated in designing the experiments, analysing the results, and drafting the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Title: Model parameterization source overview. Description: Tables containing a detailed description of the various studies in Rachuonyo South district conducted by MTC and how the data was used to parameterize the base simulation scenario. (PDF 174 KB)
Additional file 2: Title: Parameter values for the model of the mosquito feeding cycle. Description: Tables containing a detailed description of the parameter values and their source(s) for the model of the mosquito feeding cycle. (PDF 405 KB)
Additional file 3: Title: Health system parameter values. Description: Tables containing a detailed description of the parameter values and their source(s) for the model of the health system. (PDF 400 KB)
Additional file 4: Title: Description of model demographic parameters. Description: Tables containing a detailed description of the parameter values and their source(s) for the model of the demography. (PDF 259 KB)
Additional file 5: Title: Vector control intervention effective length of protection parameter values. Description: Tables containing a detailed description of the parameter values and their source(s) for effective length of protection for the model of vector control interventions. (PDF 393 KB)
Additional file 6: Title: Vector control intervention effectiveness parameter values. Description: Tables containing a detailed description of the parameter values and their source(s) for effectiveness for the model of vector control interventions. (PDF 433 KB)
Additional file 7: Title: Vector control intervention implementation parameter values. Description: Tables containing a detailed description of the parameter values and their source(s) for implementation schedule and coverage levels for the model of vector control interventions. (PDF 283 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Stuckey, E.M., Stevenson, J.C., Cooke, M.K. et al. Simulation of malaria epidemiology and control in the highlands of western Kenya. Malar J 11, 357 (2012). https://doi.org/10.1186/1475-2875-11-357