The contribution of malaria control interventions on spatio-temporal changes of parasitaemia risk in Uganda during 2009–2014

Parasites & Vectors, Sep 2017

In Uganda, malaria vector control interventions and case management with Artemisinin Combination Therapies (ACTs) have been scaled up over the last few years as a result of increased funding. Data on parasitaemia prevalence among children less than 5 years old and coverage of interventions was collected during the first two Malaria Indicator Surveys (MIS) conducted in 2009 and 2014, respectively. In this study, we quantify the effects of control interventions on parasitaemia risk changes between the two MIS in a spatio-temporal analysis. Bayesian geostatistical and temporal models were fitted on the MIS data of 2009 and 2014. The models took into account geographical misalignment in the locations of the two surveys and adjusted for climatic changes and socio-economic differentials. Parasitaemia risk was predicted over a 2 × 2 km2 grid and the number of infected children less than 5 years old was estimated. Geostatistical variable selection was applied to identify the most important ITN coverage indicators. A spatially varying coefficient model was used to estimate intervention effects at sub-national level. The coverage of Insecticide Treated Nets (ITNs) and ACTs more than doubled at country and sub-national levels during the period 2009–2014. The coverage of Indoor Residual Spraying (IRS) remained static at all levels. ITNs, IRS, and ACTs were associated with a reduction in parasitaemia odds of 19% (95% BCI: 18–29%), 78% (95% BCI: 67–84%), and 34% (95% BCI: 28–66%), respectively. Intervention effects varied with region. Higher socio-economic status and living in urban areas were associated with parasitaemia odds reduction of 46% (95% BCI: 0.51–0.57) and 57% (95% BCI: 0.40–0.53), respectively. The probability of parasitaemia risk decline in the country was 85% and varied from 70% in the North-East region to 100% in Kampala region. The estimated number of children infected with malaria declined from 2,480,373 in 2009 to 825,636 in 2014. Interventions have had a strong effect on the decline of parasitaemia risk in Uganda during 2009–2014, albeit with varying magnitude in the regions. This success should be sustained by optimizing ITN coverage to achieve universal coverage.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://parasitesandvectors.biomedcentral.com/track/pdf/10.1186/s13071-017-2393-0?site=parasitesandvectors.biomedcentral.com

The contribution of malaria control interventions on spatio-temporal changes of parasitaemia risk in Uganda during 2009–2014

Ssempiira et al. Parasites & Vectors The contribution of malaria control interventions on spatio-temporal changes of parasitaemia risk in Uganda during 2009-2014 Julius Ssempiira 0 1 3 Betty Nambuusi 0 1 3 John Kissa 2 Bosco Agaba 2 Fredrick Makumbi 3 Simon Kasasa 3 Penelope Vounatsou 0 1 0 University of Basel , Basel , Switzerland 1 Swiss Tropical and Public Health Institute , Basel , Switzerland 2 Ministry of Health , Kampala , Uganda 3 Makerere University School of Public Health , Kampala , Uganda Background: In Uganda, malaria vector control interventions and case management with Artemisinin Combination Therapies (ACTs) have been scaled up over the last few years as a result of increased funding. Data on parasitaemia prevalence among children less than 5 years old and coverage of interventions was collected during the first two Malaria Indicator Surveys (MIS) conducted in 2009 and 2014, respectively. In this study, we quantify the effects of control interventions on parasitaemia risk changes between the two MIS in a spatio-temporal analysis. Methods: Bayesian geostatistical and temporal models were fitted on the MIS data of 2009 and 2014. The models took into account geographical misalignment in the locations of the two surveys and adjusted for climatic changes and socio-economic differentials. Parasitaemia risk was predicted over a 2 × 2 km2 grid and the number of infected children less than 5 years old was estimated. Geostatistical variable selection was applied to identify the most important ITN coverage indicators. A spatially varying coefficient model was used to estimate intervention effects at sub-national level. Results: The coverage of Insecticide Treated Nets (ITNs) and ACTs more than doubled at country and sub-national levels during the period 2009-2014. The coverage of Indoor Residual Spraying (IRS) remained static at all levels. ITNs, IRS, and ACTs were associated with a reduction in parasitaemia odds of 19% (95% BCI: 18-29%), 78% (95% BCI: 67-84%), and 34% (95% BCI: 28-66%), respectively. Intervention effects varied with region. Higher socio-economic status and living in urban areas were associated with parasitaemia odds reduction of 46% (95% BCI: 0.51-0.57) and 57% (95% BCI: 0.40-0.53), respectively. The probability of parasitaemia risk decline in the country was 85% and varied from 70% in the North-East region to 100% in Kampala region. The estimated number of children infected with malaria declined from 2,480,373 in 2009 to 825,636 in 2014. Conclusions: Interventions have had a strong effect on the decline of parasitaemia risk in Uganda during 2009-2014, albeit with varying magnitude in the regions. This success should be sustained by optimizing ITN coverage to achieve universal coverage. Malaria; Malaria indicator survey; Spatio-temporal; Parasitaemia; ITNs; IRS; ACTs; Spatially varying; Bayesian kriging; Malaria interventions Background Although malaria is still a leading global health problem, its burden has been on a decline in recent years [ 1 ]. This decline which started in the early 1990s prior to the global campaign of scaling up of control interventions in mid2000s continued through the post-scale-up period [ 2 ]. The downward trend of malaria burden in the preintervention period notwithstanding, sufficient evidence from randomized trials and field settings indicate that malaria decline during the post-scale-up period has been unprecedented [ 2–5 ]. For instance in sub-Saharan Africa (SSA) parasitaemia prevalence declined from 17% in 2010 to 13% in 2015 [1]. Also, during the period 2000–2015, declines in global malaria incidence and deaths of up to 37 and 60%, respectively were reported [ 3, 6 ]. These declines were mainly attributed to the impact of Insecticide Treated Nets (ITNs) and malaria case management with Artemisinin Combination Therapies (ACTs). In spite of these higher declines in malaria at global level, slower declines were reported in the 15 most high burden countries, the majority of which are situated in SSA [ 1 ]. This region bears the heaviest burden and accounts for an estimated 90% of all malaria deaths mainly among children less than 5 years old. Uganda is ranked fourth among these high malaria burden countries and has some of the highest malaria transmission rates in the world [7]. Since 2006, Roll Back Malaria (RBM) has funded malaria control and prevention activities in the country and periodically supports the conducting of Malaria Indicator Surveys (MIS) [8]. The MIS are standardized nationally representative surveys that collect high quality data for estimating the prevalence of parasitaemia risk in children less than 5 years old and track the progress of interventions coverage. To date, two MIS have been conducted in Uganda; MIS 2009 and MIS 2014–15 [9, 10]. Findings from the first MIS revealed a high parasitaemia risk in most regions. Malaria was hyperendemic (prevalence 50–75%) in three regions, mesoendemic (prevalence 10–50%) in six, and only hypoendemic (prevalence < 10%) in one region [10]. Results of the second MIS showed tremendous improvement in the coverage of ITNs and ACTs intervention at all levels and a reduction of parasitaemia risk of 50%. Additionally, parasitaemia risk in the majority of regions had declined to mesoendemic and hypoendemic proportions [9]. The true effect of each intervention on parasitaemia reduction is not known at national and sub-national level, and yet a new framework has been adopted by the Ministry of Health (MoH) to speed up malaria control efforts. In this framework known as Uganda Malaria Reduction Strategic Plan (UMRSP) 2014–2020, ambitious targets have been set to reduce malaria mortality to near zero, morbidity to 30 cases per 1000 population, and parasite prevalence to less than 7% [8]. To achieve these targets and ensure efficient use of scarce resources and effective programming and implementation, it is vital to understand the effect that each intervention has had on parasitaemia risk decline. Declines in malaria parasitaemia risk, morbidity and mortality have been achieved in other malaria endemic countries following scaling up of control interventions. Bhatt et al. [ 3 ] reported a reduction of 50% in Plasmodium falciparum prevalence and 40% in incidence of clinical disease in endemic African countries between 2000 and 2015. Similarly, the number of malaria cases and deaths decreased by more than 50% in southern African countries after introducing interventions during 2000–2008 [11]. In the Kilifi district of Kenya, parasitaemia prevalence declined from 35 to 1% after a mass distribution of ITNs and ACTs [11]. Also, Giardina et al. [12] demonstrated that ITNs and IRS were significantly associated with parasitaemia risk reduction in Rwanda, Tanzania, Senegal, Angola, Liberia and Mozambique. Our study aims to estimate spatio-temporal trends of parasitaemia risk changes among children less than 5 years old in Uganda during 2009–2014, and to determine the effect of interventions on parasitaemia risk decline at national and subnational levels. We analyzed MIS data using Bayesian spatio-temporal geostatistical models. The results from this study provide insight on the effectiveness of interventions and can be used by MoH and Malaria Control Program (MCP) to evaluate interventions and optimize resources for achievement of objectives of UMRSP 2014–2020. Methods Country profile Uganda is located in the great lakes region in East Africa neighboring Kenya, Tanzania, Rwanda, Democratic Republic of Congo, and South Sudan. It has a population of 37.1 million, all of which are at risk of malaria. Malaria is the leading cause of morbidity and mortality in the country, accounting for 3,631,939 (4,400,000–12,000,000) cases and 5921 (5300–17,000) deaths in 2015 [13]. The most dominant malaria parasite is Plasmodium falciparum, and the major transmission vectors are Anopheles gambiae and Anopheles funestus. In recent times, vector resistance to both pyrethroid and carbamates has been reported. Data sources Parasitological and interventions data were obtained from the MIS data of 2009 and 2014–2015. The two surveys were conducted at the peak of a high malaria transmission season towards the end of the long rainy season (December 2009 and December 2014–January 2015, respectively). The MIS are nationally representative surveys which employ a two-stage stratified cluster design. The clusters also known as census enumeration areas are selected at first stage with probability proportional-to-size sampling, and households are selected at second stage using systematic sampling. The surveys are designed to provide information on key malaria control indicators, such as the proportion of households having at least one ITN, the proportion of children under 5 years of age who slept under an ITN the previous night. Also, the survey is designed to produce key indicator estimates for urban and rural strata separately, as well as for the 10 regions/domains that constitute the country. The regions are: Kampala, Central 1, Central 2, East-Central, Mid-North, Mid-Western, North-East, South-Western and West Nile. At the first stage of sampling, 170 and 210 clusters were selected in 2009 and 2014, respectively. At the second stage, 28 households were selected from each cluster in both surveys resulting in a total of 4000 and 5880 households selected in the first and second survey, respectively [9, 10]. Coverage of ITNs was defined in terms of ownership and use indicators that were generated from data captured on the survey tools using standard definitions [14]. The following ITN ownership indicators were defined; proportion of households with at least one ITN, proportion of households with one ITN for every two people, and proportion of population with access to an ITN within their household. The ITN use indicators were: proportion of children less than 5 years of age who slept under an ITN, proportion of population that slept under an ITN, and proportion of ITNs used the night preceding the survey. IRS coverage was defined as the proportion of households that were sprayed during the last 12 and 6 months in the MIS 2009 and MIS 2014–2015, respectively. The wealth index derived from household possessions was used as a socioeconomic proxy. A case management indicator was defined as the proportion of fever episodes in children of less than 5 years old during the last 2 weeks preceding the survey which were treated with any Artemisinin Combination Therapies (ACTs). In addition, information on the location of the cluster (i.e. rural/urban) was obtained from survey data and from the Global Rural-Urban Mapping Project (GRUMP) database [15]. The GRUMP database provides gridded data at 1 km2 spatial resolution. Malaria transmission depends on the environment which affects the disease distribution, seasonality, and transmission intensity. Environmental/climatic factors were extracted from Remote Sensing (RS) sources. Weekly day and night Land Surface Temperature (LST), bi-weekly Normalized Difference Vegetation index (NDVI) and land cover data were obtained from Moderate Resolution Imaging Spectroradiometer (MODIS) at 1 km2 spatial resolution. Dekadal rainfall data at 8 × 8 km2 resolution were extracted from the US Early Warning and Environmental Monitoring System (EWES). Altitude was obtained from the shuttle radar topographic mission using the digital elevation model. Also, distances from cluster centroids to major water bodies were estimated using ESRI’s ArcGIS 10.2.1 for Desktop. The high spatial resolution population data was downloaded from WorldPop [16]. Data from remote sensing sources was acquired for the 12 month period preceding the survey and the average (cumulative value for rainfall) was calculated and extracted for each cluster. The one-year period was considered long enough to capture the actual climatic conditions that affected malaria transmission throughout the year of the survey. Statistical analysis Bayesian geostatistical models were developed to predict parasitaemia risk at the two survey time points using environmental/climatic factors as predictors. Bayesian kriging was applied to obtain parasitaemia risk estimates over a 2 × 2 km2 resolution grid. Predictions were used to determine the probability of parasitaemia risk reduction between the two surveys. The number of children infected with malaria in the two surveys was estimated by combining high spatial resolution population data obtained from WorldPop (www.worldpop.org) with the predicted pixel-level malaria prevalence estimates. The number of children less than 5 years old was estimated by multiplying population counts by a factor of 17.7%, the proportion of population under 5 years of age [17]. Regional estimates of the number of infected children were computed by aggregating pixel-level estimates at regional level. The number of infected children per pixel was obtained by multiplying pixel-wise spatially explicit prevalence estimates with high spatial resolution population estimates of number of children less than 5 years old. In both surveys, the population-adjusted prevalence was estimated by summing up estimates of the number of infected children per pixel divided by the total estimated number of children less than 5 years old. The effects of interventions were estimated by modeling the change of parasitaemia risk between the two surveys on the logit scale as a function of the effect of intervention coverage at the second survey adjusted for socioeconomic status, cluster location, and the difference in environmental/climatic factors. Geographical misalignment of the locations between the two surveys was carried out by predicting parasitaemia risk of the first survey at the second survey locations. The prediction uncertainty was incorporated by fitting an error term in the model. A spatially varying coefficients model was used to estimate intervention effects at regional level and to account for potential interactions of interventions with endemicity level. A spike and slab geostatistical Bayesian variable selection procedure was applied to select the most important ITN and environmental predictors that explain maximum variation in the change in parasitaemia risk between 2009 and 2014 [18]. Variables with the highest inclusion probability in the model were selected. Descriptive analyses were carried out in STATA (StataCorp. 2015. Stata Statistical Software: Release 14. College Station, TX: StataCorp LP). Geostatistical modeling was implemented in OpenBUGS version 3.2.3 (Imperial College and Medical Research Council, London, UK). Since implementing Bayesian kriging in OpenBUGS is very slow especially for large grids, we implemented it in R statistical software using posterior estimates of the model parameters obtained from OpenBUGS. Maps were produced in ESRI’s ArcGIS 10.2.1 (http://www.esri.com/en-us/home). Parameter estimates were summarized by their posterior medians and their corresponding 95% Bayesian Credible Intervals (BCI). The effect of a predictor was considered to be statistically important if its 95% BCI did not include zero. Detailed explanations of the fitted statistical models are presented in Additional files 1 and 2. Results A summary of the survey data is given in Tables 1 and 2, and maps of survey locations are presented in Fig. 1. A higher number of clusters, households, and children were tested in the second survey (Table 1). Results show that at country level parasitaemia prevalence declined from 42.4% in 2009 to 19.0% in 2014, a decline of 50%. At regional level, the highest malaria reduction was observed in the regions of Kampala (91.8%), Central 1 (74.0%) and Mid-North (68.6%), and the lowest in the North-East region (30.2%) and EastCentral region (35.2%). Generally interventions coverage increased at country and regional levels (Additional file 3). At country level, ITN ownership (the proportion of households with at least one ITN and the proportion of households with at least one ITN for every two people) increased by fourfold. Among regions, the biggest increase in ITN ownership was reported in East-Central (six-fold), while the smallest was observed in Mid-North (two-fold). More so, the proportion of children less than 5 years old that slept under an ITN increased by more than two times at country level. The improvement in this indicator coverage was highest in Central 2 region (5.8 times) and lowest in North-East region (1.3 times). Overall, the proportion of fever episodes treated with ACTs increased by three times. The highest increase was achieved in South-Western, East-Central and West Nile regions where coverage increased by more than five times. The least gain in ACTs coverage was observed in Mid-North region where it increased by almost two times. The national IRS coverage remained static at 5% except in the Mid-North region where an increase of 41% was achieved. Spatio-temporal trends of parasitaemia risk during 2009–2014 The effects of the most important environmental factors identified through geostatistical variable selection are shown in Table 3. Results indicate that more environmental factors were related to parasitaemia risk in 2009 compared to 2014. Also, spatial correlation was stronger in 2009. Figure 2 depicts the predicted parasitaemia risk in 2009 and 2014 over a 2 × 2 km2 resolution grid based on the 2.5th, median, and the 97.5th percentile posterior Abbreviations: MIS Malaria Indicator Survey, TNs Insecticide Treated Nets, ACTs Artemisinin Combination Therapies, IRS Indoor Residual Spraying 4 1 .5 .3 .1 .8 .7 .8 .0 3 .0 le 0 7 6 2 8 7 6 7 .2 .6 7 i 2 2 9 7 8 7 7 7 1 9 6 N e raveog ITN th r/ond iisenm iiisnnm irtcvenon ItTenoN ItTenoN ItaTnoN rtaeunnd lrsohodw iirsveuong ItaTenoN irttyaahn ,ttrseAsTCA ilifrrtaaaaveeeoongCm lircavaeeenp lilfttsssaaeeuohhohodw lilfttsssaaeeuohhohodw leeooppw liifttssccaaeunohooppw lilfttttsaaeeuhnoohoppp illfrtsscyaaee5enhnhod IraTnN iiIfttsssxTeeeuhnogdpN lfrsssyaeeuohohoddp lilfttsssaaeeuohhohodw iIltttssaSe21Ronhhnm iffrrttssvaeeeeeeooddpw rtyaehnp :IiittrsccaTeeeeenddNIssTnN leaT2b iIrtcanod irtsaaaePm irrtPoonop irrtPooonp frrtyveeo irrtPoonop irrtPoonop irrtPooonp ltseeunpd irrtPoonop irrtPoonop irrtPonoop rsyyaedbp irrtPoonop iitcanoobm iirtveabobA predictive distributions. Estimates suggest a high parasitaemia risk in 2009 where in some areas the predicted prevalence was over 80%. In 2014, parasitaemia risk was much lower in most parts of the country except in some areas of the East-Central, North-East and West Nile regions where the burden still remained high. The probability of parasitaemia decline in the country was 85%. The highest decline in malaria occurred in the regions of Central 2 and Kampala while the least was estimated in the North-East region (Fig. 3). Overall, the number of infected children reduced from over 2,480,000 to less than 830,000 between 2009 and 2014 (Table 4). This translates into a reduction of over 66%. Reduction in the estimated number of infected children was achieved in all regions. The biggest reduction occurred in Kampala (86%), Central 1 (75%), Central 2 (74%), Mid-Eastern (71%) and Mid-North region (70%), whereas the least happened in North-East (44%). In both surveys, the highest and lowest numbers of infected children were estimated in the East-Central and Kampala regions, respectively. Overall, a reduction in population adjusted prevalence of over 26% was achieved. The highest reduction (39.4%) was observed in the East-Central region while the least one (5.0%) was registered in Kampala. Figure 4 further shows that the number of infected children in 2014 shrank considerably compared to 2009 in all regions except in the East-Central region. The map also depicts a strong statistically important reduction in concentration of infected children in Mid North region in 2014. Results from geostatistical variable selection (Table 5) indicate that the proportion of population with access to Range (km) 43.3 (12.2–57.8) 43.8 (36.3–48.2) Abbreviations: MIS Malaria Indicator Survey, LST Land Surface Temperature, NDVI Normalized Difference Vegetation Index aStatistically important effect bCut-offs before and after the slash (/) are for 2009 and 2014 respectively Day LSTb Night LST Area type NDVI Rainfallb < 27.84 / < 31.4 27.84–30.18 / 31.4–33.8 > = 30.19 / > = 33.8 Rural vs urban < 17.11 / < 17.14 17.11–18.49 / 17.14–18.79 > = 18.50 / > = 18.79 Altitudeb < 1098 an ITN in their household was the only indicator able to capture the effect of ITN interventions as it has the highest inclusion probability. This indicator was used to quantify the effect of ITNs on the parasitaemia odds change. Effects of interventions on parasitaemia odds decline The effects of interventions on the change of parasitaemia odds adjusted for socioeconomic status and changes in environmental conditions between the two surveys are showed in Table 6. Results demonstrate an important protective effect of interventions on the decrease of parasitaemia odds from 2009 to 2014. ITNs, IRS and ACTs were associated with a parasitaemia odds reduction of 19% (95% BCI: 18–29%), 78% (95% BCI: 67–84%), and 34% (95% BCI: 28–66%), respectively. Similarly, higher socio-economic status had a strong effect on parasitaemia odds reduction. More so, living in urban areas was associated with a decrease in malaria odds of 57% (95% BCI: 47–60%) compared to living in rural areas. On average, rainfall, day and night LST increased from 2009 to 2014, and these increases were significantly associated with increased parasitaemia odds. However, changes in the NDVI had no effect on changes in parasitaemia odds. Intervention effects on parasitaemia odds decline varied by region (Fig. 5). The effect of ITNs at regional level was significantly higher than the national effect in MidNorth and West Nile. ITNs’ effects were significantly lower in East-Central, Mid-Eastern, Mid-Western, and South-Western regions. Likewise, the effect of ACTs was significantly higher than the national average in most regions except in Central 1, Mid-North, Mid-Western, and West Nile regions. Discussion In this study we have determined the spatio-temporal trends of parasitaemia odds and the effect of control interventions on the change of parasitaemia risk in Uganda during 2009–2014. Furthermore, we estimated the probability Fig. 3 Probability of parasitaemia risk decline from 2009 to 2014 of parasitaemia risk decline and the number of infected children at the two survey time points. Our study results showed a strong ITNs effect on parasitaemia risk reduction during 2009–2014 following a two-fold increase in coverage in the 5 years. These results support findings in similar malaria endemic settings [ 3 ]. This protective effect can be attributed to the physical barrier provided by ITNs to block mosquitoes from infecting humans with Plasmodium sporozoites, thus preventing parasites from completing their development cycle [19]. Also, the insecticide in ITNs reduces the lifespan of vectors when they come into contact, thus decreasing the chances of transmission [6]. Furthermore, the high coverage and utilization registered in the country may have achieved a ‘mass effect’ that reduces the mosquito population and thus protects people in communities who are not using ITNs but live in close proximity to households with ITNs [20, 21]. The high increase in ITNs coverage can be credited to increased donor support that funded ITNs purchase and distribution through effective distribution outreach channels [8]. These channels include mass distribution campaigns, antenatal care clinics, Expanded Program for Immunization (EPI), and commercial sale of subsidized ITNs through the private sector. These distribution channels have had an immediate success of raising the proportion of households possessing at least one ITN from less than 50% to more than 90%. In spite of the high ITN coverage across the country, ITN effects on parasitaemia odds reduction varied with region. Effects were highest in regions which were initially the most high burdened in 2009. The varying effects of interventions could be explained by regional heterogeneities in malaria transmission rates [22], ecology, and access to health services [23]. Furthermore, case management with ACTs was strongly associated with parasitaemia risk reduction following a three-fold increase in coverage during the study period. Prompt treatment of malaria with ACTs suppresses and kills malaria parasites in the body which prevents progression to Abbreviations: ITNs Insecticide Treated Nets, ACTs Artemisinin Combination Therapies, IRS Indoor Residual Spraying, LST Land Surface Temperature, NDVI Normalized Difference Vegetation Index aStatistically important effect Covariate Difference in LST (day) Difference in LST (night) Difference in NDVI Difference in rainfall Area type (urban vs rural) Wealth index ITN IRS ACTs Spatial variance Range (km) OR (95% BCI) 1.10 (1.02–1.13)a 1.09 (1.03–1.18)a 1.00 (0.94–1.08) 1.14 (1.08–1.23)a 0.43 (0.40–0.53)a 0.54 (0.51–0.57)a 0.81 (0.71–0.82)a 0.22 (0.16–0.33)a 0.66 (0.34–0.72)a 0.63 (0.56–0.76) 35.4 (24.3–37.0) Probability of inclusion (%) severe disease, thus reducing transmission and subsequently parasitaemia load in the population [24]. In line with our study findings, Bhatt et al. [ 3 ] also found that ACTs together with ITNs were the most impactful interventions on malaria risk reduction in African endemic countries during 2000–2015. Also, effects of ACTs also varied with region. However, despite the two-fold increase in ACTs coverage in the 5 years, its coverage was still lower than targeted. This could possibly be attributed to supply chain constraints [25], the semi-regulated private health facilities and drug stores and the inadequate laboratory diagnostic capacity in most of the lower level facilities [8]. Indoor residual house spraying also had a very strong effect on parasitaemia odds reduction despite its coverage remaining static between 2009 and 2014. The endophilic behavior of the predominant Anopheles spp. mosquitoes makes this intervention highly effective in Uganda as vectors are killed by the insecticide as they rest on house walls after taking a blood meal [26]. The static coverage is perhaps explained by the high costs involved in IRS implementation. This prompted NMCP to roll out IRS gradually initially starting in 2009 with the 10 most high malaria burden districts located in the Mid-North region [8]. Following a significant reduction in malaria transmission in the 10 districts [27], IRS was later extended to another 14 high burden districts in the North-East, MidEastern, and East-Central regions. Effectiveness of IRS on malaria risk reduction has been reported in other studies in Uganda [27], Kenya [ 28 ], Bioko, Equatorial Guinea and Mozambique [29]. Our results further showed that urban areas were associated with a decreased parasitaemia risk compared to rural areas. This could be explained by uneven access to healthcare services between urban and rural areas in developing countries [30]. In Uganda, lower level health facilities, which are the major source of health services in rural areas, are poorly equipped and understaffed [31]. On the other hand, urban areas are served by a much bigger network of better equipped higher level facilities both public and private. Indeed urbanization is one of the reasons that has been suggested as a strong possible causal factor of the downward trend of malaria risk in the pre-intervention period [32]. The effect of urbanization on socio-economic and landscape changes mitigates the risk of malaria transmission. The inverse relationship between urbanization and malaria risk has also been reported in other malaria endemic settings [32–35]. Higher socio-economic status was strongly associated with parasitaemia odds reduction. Related to this finding, our results also showed that the highest probability of parasitaemia decline was attained in Kampala region and the lowest in the North-East region. The former is the capital city and the most developed region, while the latter is the least developed and most hard-to-reach region in Uganda. Socio-economic status affects the ability to afford healthcare services, better housing conditions, and knowledge of malaria prevention [36], which are important determinants of severity and outcome of the disease. These results are in agreement with other studies that reported a higher burden of malaria among poor countries [37] and in hard-to-reach areas [6, 37]. This finding augments evidence that malaria is a disease associated with poverty [38, 39] and low socio-economic development [39–42]. Furthermore, increased land surface temperature and rainfall between 2009 and 2014 were associated with a higher parasitaemia risk. This result is expected since malaria is a vector-borne disease sensitive to changes in climatic conditions [ 2 ]. Temperature influences the speed of development of mosquitoes and Plasmodium parasites [43]. Rainfall is the most important driver of mosquito population dynamics and malaria transmission because it provides the optimal humidity and medium for mosquito fertilization and breeding [44, 45]. Although a reduction in parasitaemia risk was achieved in all regions, nevertheless, parasitaemia risk was still high in the regions of North-East, West Nile, and East-Central compared to other regions. This disproportionately high risk in these regions in spite of the high intervention coverage might be attributed to low socio-economic development [46], and limited access to health services [23]. In the case of East-Central region, rice growing practiced in this region has been documented as a potential driver of malaria risk transmission due to the large swamps that provide a favorable habitat for mosquito breeding [47]. Similarly, other studies have reported a higher malaria risk in settings with low socio-economic status [42], poor access to health services [39], and rice paddies [48]. The strong reduction in the estimated number of malaria-infected children may also underline the effect of increases in interventions coverage [9], urbanization [49], and generally improving socio-economic conditions [50]. Conclusions Our study demonstrates that malaria control interventions have had a strong effect on the decline of parasitaemia risk in Uganda during 2009–2014, albeit with varying magnitude in the regions. This success should be sustained by optimizing ITN coverage to achieve universal coverage and by timely replacing worn-out ITNs. NMCP should sustain the malaria prevention awareness campaigns through the use of Information, Education and Communication (IEC) materials to further promote the use of ITNs. In the high burden districts where IRS implementation is on-going, efforts should be made to ensure that all households are sprayed periodically every 6 months. NMCP should address the problems limiting ACTs coverage scale-up by providing free RDTs to all healthcare providers in line with the WHO ‘Test and Treat’ campaign, and increasing supervision for private health facilities. The varying intervention effects in different regions maybe an indication that interventions work differently in different regions of the country. This therefore calls for a better understanding of the environmental and entomological conditions in each region to tailor a combination of interventions suitable to local settings that will have maximum reduction on transmission. Also, in the regions where the risk remains disproportionately high, NMCP needs to conduct specific studies to understand human and/or vector behavior responsible for this problem. In these regions, other tools should be introduced such as chemoprevention especially in the high risk group of children less than 5 years and mass drug administration to reduce the parasite load in the population. In order to maximize intervention effects and avert reversal in malaria risk reduction, government and donor funded poverty reduction programs should prioritize regions/districts where socio-economic conditions are low. In summary, the ambitious targets of UMRSP 2014–2020 can be achieved if the country commits to implementing an integrated package to cover all aspects of disease prevention, management, and health. However, this will only be possible if the current funding portfolio is increased from the contemporary less than $1 average per head per year to the recommended $4 per head per year [51] equivalent to $140 million per year. Additional files Additional file 1: Details of statistical models to estimate parasitaemia risk, effects of interventions on the change of parasitaemia risk, and spatially varying interventions effects. (DOCX 40 kb) Additional file 2: Joint posterior distributions of the fitted statistical models. (DOCX 32 kb) Additional file 3: Malaria intervention coverage in 2009 and 2014. Percentage of households with one ITN (a), percentage of households with at least 1 ITN for every two people (b), percentage of population with access to an ITN (c), percentage of population that slept under an ITN the previous night (d), percentage of children less than 5 years who slept under an ITN the previous night (e), proportion of fever episodes treated with any ACT (f). (PDF 211 kb) Abbreviations ACTs: Artemisinin combination therapies; BCI: Bayesian credible intervals; DHS: Demographic health survey; EWES: Environmental monitoring system; GRUMP: Global rural-urban mapping project; IRS: Indoor residual spraying; ITNs: Insecticide treated nets; LST: Land surface temperature; MIS: Malaria indicator surveys; MODIS: Moderate resolution imaging spectroradiometer; MoH: Ministry of Health; NDVI: Normalized difference vegetation index; NMCP: National malaria control program; PMI: President’s malaria initiative; RS: Remote sensing; SSA: Sub-Saharan Africa; UMRSP: Uganda malaria reduction strategic plan (UMRSP); WHO: World Health Organization Acknowledgments We are grateful to Uganda Ministry of Health, MCP, Uganda Bureau of Statistics (UBOS), Makerere University School of Public Health, DHS MEASURE, PMI and the Global Fund. Funding This research was supported and funded by the Swiss Programme for Research on Global Issues for Development (r4d) project no. IZ01Z0–147,286 and the European Research Council (ERC) advanced grant project no. 323180. Availability of data and materials The DHS MEASURE program prohibits researchers from redistributing data as per their “Dataset Terms of Use”. However, the data are available in the DHS MEASURE program website (www.dhsprogram.com) upon request following data access instructions (http://dhsprogram.com/data/Access-Instructions.cfm). Also, data can be requested through the following contact; Tel: (301) 572–0851, E-mail: . Authors’ contributions JS developed methodology, analyzed and synthesized data, fitted models, carried out data validation, and wrote the manuscript; BN participated in data analysis and synthesis; JK carried out data curation and participated in manuscript writing; BA carried out data curation and participated in manuscript writing; FM formulated research goals and objectives, participated in the process of acquisition of project financial support, and manuscript writing; SK formulated research goals and objectives, planned, coordinated, and executed research, and manuscript writing; PV formulated research goals and objectives, planned, coordinated, and executed research, spearheaded study methodology development, and manuscript writing. All authors read and approved the final manuscript. Ethics approval and consent to participate In this study we analyzed secondary data made available by the Demographic Health Survey (DHS) MEASURE. According to survey protocols and related documents of the two surveys, ethical approval was obtained from the Institutional Review Board of International Consulting Firm (ICF) of Calverton, Maryland, USA, and also from Makerere University School of Biomedical Sciences Higher Degrees Research and Ethics committee (SBS-HDREC), and the Uganda National Council for Science and Technology (UNCST). Details of ethical clearance are published in the Uganda MIS 2009 and MIS 2014–15 reports for the first and second survey, respectively [9, 10]. Consent for publication Not applicable. Competing interests The authors declare that they have no competing interests. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. 1. World Health Organization. World malaria report: 2016 . 2016 . Available: http://apps.who.int/iris/bitstream/10665/252038/1/9789241511711-eng.pdf. 2. Snow RW , Kibuchi E , Karuri SW , Sang G , Gitonga CW , Mwandawiro C , et al. Changing malaria prevalence on the Kenyan Coast since 1974: climate, drugs and vector control . PLoS One . 2015 ; 10 : e0128792 . 3. Bhatt S , Weiss DJ , Cameron E , Bisanzio D , Mappin B , Dalrymple U , et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015 . Nature . 2015 ; 526 : 207 - 11 . Impact of artemisinin-based combination therapy and insecticide-treated nets on malaria burden in Zanzibar . PLoS Med . 2007 ; 4 : e309 . Cochrane Database Syst Rev . 2004 ; (2):CD000363. WHO , UNICEF. WHO | Achieving the malaria MDG target: reversing the incidence of malaria 2000-2015 . In: WHO 2015. Available: http://www.who. int/malaria/publications/atoz/9789241509442/en/. Accessed 10 June 2017 . Measures of malaria burden after long-lasting insecticidal net distribution and indoor residual spraying at three sites in Uganda: A prospective observational study . PLoS Med . 2016 ; 13 : e1002167 . National Malaria Control Program . Uganda Malaria Operational Plan . 2016 . Available: https://www.pmi.gov/docs/default-source/default-documentlibrary/malaria-operational-plans/fy16/fy-2016 - uganda -malaria-operationalplan . pdf?sfvrsn=5 Uganda Bureau of Statistics and ICF International. Uganda Malaria Indicator Survey 2014 - 15 . Kampala, Uganda, and Rockville , Maryland: UBOS and ICF International 2015 . Available: https://dhsprogram.com/pubs/pdf/MIS21/ MIS21.pdf. Accessed 30 May 2016 . Uganda Bureau of Statistics and ICF International. Uganda Malaria Indicator Survey 2009 . Kampala, Uganda, and Rockville , Maryland: UBOS and ICF International . 2010 Available: https://dhsprogram.com/pubs/pdf/MIS6/MIS6. pdf. Accessed 30 May 2016 . O'Meara WP , Mangeni JN , Steketee R , Greenwood B . Changes in the burden of malaria in sub-Saharan Africa . Lancet Infect Dis . 2010 ; 10 : 545 - 55 . Giardina F , Kasasa S , Sié A , Utzinger J , Tanner M , Vounatsou P. Effects of vectorcontrol interventions on changes in risk of malaria parasitaemia in sub-Saharan Africa: a spatial and temporal analysis . Lancet Glob Health . 2014 ; 2 : e601 - 15 . WHO. World Malaria Report 2015 . Geneva: WHO; 2015 . World Health Organisation. Household Survey Indicators for Malaria Control 2013 . Available: http://www.malariasurveys. org/documents/ Household%20Survey%20Indicators%20for%20Malaria%20Control.pdf. Global Rural-Urban Mapping Project (GRUMP), v1 | SEDAC . Available: http:// sedac.ciesin.columbia.edu/data/collection/grump-v1. Accessed 3 Feb 2017 . org.uk/. Accessed 30 May 2016 . Uganda Bureau of Statistics. The National Population and Housing Census 2014 - Main Report , Kampala, Uganda 2016 . Available: http://www.ubos.org/ onlinefiles/uploads/ubos/NPHC/ NPHC%202014 %20FINAL%20RESULTS%20REPORT.pdf. Accessed 14 Sep 2016 . Modelling the geographical distribution of soil-transmitted helminth infections in Bolivia . Parasit Vectors . 2013 ; 6 : 152 . Bueno-Mari R , Jimenez-Peydro R . Global change and human vulnerability to vector-borne diseases . Lausane: Frontiers Media SA ; 2010 . Louis VR , Bals J , Tiendrebéogo J , Bountogo M , Ramroth H , De Allegri M , et al. Long-term effects of malaria prevention with insecticide-treated mosquito nets on morbidity and mortality in African children: randomised controlled trial . Tropical Med Int Health . 2012 ; 17 : 733 - 41 . Maxwell CA , Msuya E , Sudi M , Njunwa KJ , Carneiro IA , Curtis CF . Effect of community-wide use of insecticide-treated nets for 3-4 years on malarial morbidity in Tanzania . Trop Med Int Health TM IH . 2002 ; 7 : 1003 - 8 . Okello PE , Van Bortel W , Byaruhanga AM , Correwyn A , Roelants P , Talisuna A , et al. Variation in malaria transmission intensity in seven sites throughout Uganda . Am J Trop Med Hyg . 2006 ; 75 : 219 - 25 . Yeka. Malaria in Uganda: Challenges to control on the long road to elimination: I. Epidemiology and current control efforts 2012 . Available: http://www.sciencedirect.com/science/article/pii/S0001706X11000611. Accessed 28 Nov 2016 . Baird JK . Real-world therapies and the problem of vivax malaria . N Engl J Med . 2008 ; 359 : 2601 - 3 . Access to and utilisation of health services for the poor in Uganda: a systematic review of available evidence . Trans R Soc Trop Med Hyg . 2008 ; 102 : 1067 - 74 . Becker N , Petric D , Zgomba M , Boase C , Minoo M , Dahl C , et al. Mosquitoes and their control . New York: Springer; 2010 . Bukirwa H , Yau V , Kigozi R , Filler S , Quick L , Lugemwa M , et al. Assessing the impact of indoor residual spraying on malaria morbidity using a sentinel site surveillance system in Western Uganda . Am J Trop Med Hyg . 2009 ; 81 : 611 - 4 . 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 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 . Combining indoor residual spraying and insecticide-treated net interventions . Am J Trop Med Hyg . 2009 ; 81 : 519 - 24 . Dolea C. Increasing access to health workers in remote and rural areas through improved retention: global policy recommendations . Geneva: World Health Organization; 2010 . Changes in utilization of health services among poor and rural residents in Uganda: are reforms benefitting the poor? Int J Equity Health . 2009 ; 8 : 39 . Tatem AJ , Gething PW , Smith DL , Hay SI . Urbanization and the global malaria recession . Malar J . 2013 ; 12 : 133 . Omumbo JA , Guerra CA , Hay SI , Snow RW . The influence of urbanisation on measures of Plasmodium falciparum infection prevalence in East Africa . Acta Trop . 2005 ; 93 : 11 - 21 . Ramroth H , Ndugwa RP , Müller O , Yé Y , Sié A , Kouyaté B , et al. Decreasing childhood mortality and increasing proportion of malaria deaths in rural Burkina Faso . Glob Health Action . 2009 ;2: 1909 . Rapid urban malaria appraisal (RUMA) I: Epidemiology of urban malaria in Ouagadougou . Malar J . 2005 ; 4 : 43 . Yadav K , Dhiman S , Rabha B , Saikia P , Veer V . Socio-economic determinants for malaria transmission risk in an endemic primary health centre in Assam, India . Infect Dis Poverty . 2014 ; 3 : 19 . Snow RW , Marsh K. Malaria in Africa: progress and prospects in the decade since the Abuja Declaration . Lancet . 2010 ; 376 : 137 - 9 . Tanner M , de Savigny D. Malaria eradication back on the table . Bull World Health Organ . 2008 ; 86 : 82 . Malaria: progress, perils, and prospects for eradication . J Clin Invest . 2008 ; 118 : 1266 - 76 . Protopopoff N , Bortel WV , Speybroeck N , Geertruyden J-PV , Baza D , D'Alessandro U , et al. Ranking malaria risk factors to guide malaria control efforts in African highlands . PLoS One . 2009 ; 4 : e8022 . Gullan PJ , Cranston PS . The Insects: An Outline of Entomology . 5th ed. New York: Wiley-Blackwell; 2014 . Githeko AK , Ndegwa W. Predicting malaria epidemics in the Kenyan highlands using climate data: a tool for decision makers . Glob Change Hum Health . 2001 ; 2 : 54 - 63 . Kynast-Wolf G , Hammer GP , Müller O , Kouyaté B , Becher H . Season of death and birth predict patterns of mortality in Burkina Faso . Int J Epidemiol . 2006 ; 35 : 427 - 35 . Ministry of Finance. Uganda poverty status report 2014 . Available: zotero:// report/items/0_H927KNFD/html/report.html. Accessed 27 Nov 2016 . Pullan RL , Bukirwa H , Staedke SG , Snow RW , Brooker S. Plasmodium infection and its risk factors in eastern Uganda . Malar J . 2010 ; 9 : 2 . Diboulo E , Sié A , Vounatsou P. Assessing the effects of malaria interventions on the geographical distribution of parasitaemia risk in Burkina Faso . Malar J . 2016 ; 15 : 228 . Malar J. 2015 ; 14 : 374 . Tusting LS , Rek J , Arinaitwe E , Staedke SG , Kamya MR , Cano J , et al. Why is malaria associated with poverty? Findings from a cohort study in rural Uganda . Infect Dis Poverty . 2016 ; 5 : 78 . Teklehaimanot A , McCord GC , Sachs JD . Scaling up malaria control in Africa: an economic and epidemiological assessment . Am J Trop Med Hyg . 2007 ; 77 : 138 - 44 .


This is a preview of a remote PDF: https://parasitesandvectors.biomedcentral.com/track/pdf/10.1186/s13071-017-2393-0?site=parasitesandvectors.biomedcentral.com

Julius Ssempiira, Betty Nambuusi, John Kissa, Bosco Agaba, Fredrick Makumbi, Simon Kasasa, Penelope Vounatsou. The contribution of malaria control interventions on spatio-temporal changes of parasitaemia risk in Uganda during 2009–2014, Parasites & Vectors, 2017, 450, DOI: 10.1186/s13071-017-2393-0