Malaria Journal volume 22, Article number: 16 (2023 ) Cite this article
To achieve malaria elimination it is essential to understand the impact of insecticide-treated net (ITNs) programmes. Here, the impact of ITN access and use on malaria prevalence in children in Malawi was investigated using Malaria Indicator Survey (MIS) data.
MIS data from 2012, 2014 and 2017 were used to investigate the relationship between malaria prevalence in children (6–59 months) and ITN use. Generalized linear modelling (GLM), geostatistical mixed regression modelling and non-stationary GLM were undertaken to evaluate trends, spatial patterns and local dynamics, respectively.
Malaria prevalence in Malawi was 27.1% (95% CI 23.1–31.2%) in 2012 and similar in both 2014 (32.1%, 95% CI 25.5–38.7) and 2017 (23.9%, 95% CI 20.3–27.4%). ITN coverage and use increased during the same time period, with household ITN access growing from 19.0% (95% CI 15.6–22.3%) of households with at least 1 ITN for every 2 people sleeping in the house the night before to 41.7% (95% CI 39.1–44.4%) and ITN use from 41.1% (95% CI 37.3–44.9%) of the population sleeping under an ITN the previous night to 57.4% (95% CI 55.0–59.9%). Both the geostatistical and non-stationary GLM regression models showed child malaria prevalence had a negative association with ITN population access and a positive association with ITN use although affected by large uncertainties. The non-stationary GLM highlighted the spatital heterogeneity in the relationship between childhood malaria and ITN dynamics across the country.
Malaria prevalence in children under five had a negative association with ITN population access and a positive association with ITN use, with spatial heterogeneity in these relationships across Malawi. This study presents an important modelling approach that allows malaria control programmes to spatially disentangle the impact of interventions on malaria cases.
Malaria is one of the most important causes of morbidity and mortality in Malawi. Since 2007, mass distribution of insecticide-treated nets (ITNs) to the demographically most vulnerable population groups has been a major part of Malawi’s vector control efforts. This, in combination with improved diagnosis and treatment of cases, has resulted in a 36% reduction in malaria cases, from an estimated 5.6 million cases in 2010 to 3.6 million in 2016 [1]. Despite this progress, continued investment in malaria control has not led to a further decrease in cases [1, 2] with an estimated 4.3 million cases still occurring in 2020.
The Malawi National Malaria Control Programme (NMCP) aimed to reduce malaria incidence by at least 50% from a 2016 baseline of 386 per 1000 population to 193 per 1000, and reduce malaria deaths by at least 50% from 23 per 100,000 population to 12 per 100,000 population by 2022 [3]. To reach these targets, it is necessary to understand the impact of current control activities and tailor future vector control to local epidemiological and entomological dynamics [4] in the context of a fast growing population. The Malawian population has grown almost 30% during from 2012 to 2022 to 20.4 million [5]. To control the malaria vectors, in 2007, the government started distributing long-lasting insecticide-treated nets (LLINs). It is generally accepted that ITN coverage has helped decrease malaria prevalence in Malawi [6,7,8,9]. However, retrospective studies have also found that a 13% increase in bed net access from 2012 to 2014 was not associated with a reduction in malaria burden in children [10]. This was corroborated by a similar study investigating overall trends within the 2012 and 2014 malaria indicator survey (MIS) data [11]. Comparison of data from 2004 and 2016 also showed that community malaria prevalence was not related to ITN use [12]. No personal, and limited community, protection from ITNs was found in a field study from 2012 [13]. These inconsistent conclusions suggest that ITN access and use may have a heterogeneous impact on malaria prevalence.
One of the main challenges for malaria elimination is the heterogeneity of the current malaria landscape [1, 14, 15]. It is not well understood why this heterogeneity has emerged, although the possible varied efficacy of ITNs could be partly responsible. ITNs can impact areas differently due to different vector population compositions and behaviours, climate variation and the presence of resistance [16,17,18]. Furthermore, gaps may exist between policy and implementation [19], with human behaviour one of the most complex variables involved in malaria transmission. The local population might accept, but not use and maintain nets, use nets for other purposes or migrate to areas with higher malaria risk [20,21,22]. Further reasons include social factors, such as autonomy in health care decisions [23], bed net integrity and insecticide degradation [24].
The best method for measuring the efficacy of ITNs directly is randomized controlled trials, ideally conducted across different parts of a country. No such randomized control trials have been conducted in Malawi due to the unethical nature of withholding nets from a proportion of the population. However, an important alternative source of information exists in the MIS. A search on the National Library of Medicine identifies more then 2,000 papers that have used MIS data in some capacity (search date 16–02-2022). These routine, large-scale household surveys are designed to produce snapshots of the malaria situation at national, regional and urban/rural levels. In Malawi, national MIS were conducted in 2010, 2012, 2014 and 2017 [25,26,27]. These MIS, together, capture the malaria situation and dynamics throughout the country and represent an essential source of information for policy development [1]. For example, MIS studies helped recognize the possible ineffectiveness of ITNs [10, 11], which led to the recent addition of pyrethroid-Piperonyl butoxide (PBO) nets to the vector control programme. The latest MIS survey was conducted in Malawi in April 2021. Data are currently being analysed and have not been released to the public.
Even with the recent addition of PBO and dual active ingredient nets to the vector control programme, traditional permethrin nets will likely remain an important part of control efforts in Malawi due to the uncertain durability of the next generation bed nets[28], their additional costs and unknown acceptance by the local population. In this study, the spatial relationship between ITN access and use, and malaria prevalence in children was investigated using MIS data from 2012, 2014 and 2017.
Malawi is a southern African country with an estimated population of 18.6 million in 2019 [29]. The country consists of the Northern, Central and Southern regions, further divided into 28 districts. Malawi has three seasons, a rainy season from November to April, a cool-dry season from May to August and a hot-dry season from August to November. This sub-tropical climate is favourable for transmission of malaria. The main malaria vector, Anopheles funestus, is present throughout the country and throughout the year. The population along Lake Malawi and in the southern lowland are especially at high risk [15], as their environment is also ideal for the secondary malaria vectors Anopheles gambiae sensu stricto (s.s.) and Anopheles arabiensis [15, 30]. Even though An. funestus and An. gambiae sensu lato (s.l.) were susceptible to all insecticides in 2007, by 2010 pyrethroid resistance was found in both species [16, 31]. At present, it is suspected that pyrethroid and carbamate resistance is widespread across the country [32]. In 2007, the Malawian government started the distribution of LLINs. In 2012, 41% of the population owned at least one ITN with 29% sleeping under it the preceding night [25]. In the same year, the aim of universal bed net coverage (one LLIN per 1.8 people) was expressed. To achieve this, in 2014, 7 million nets were distributed throughout the country, further complemented by an additional 8 million in 2016 [1][1]. During all distribution campaigns hyperendemic and endemic areas were prioritized. To supplement the ITN campaigns, ITNs are also routinely distributed through the antenatal care service. Sporadic indoor residual spray (IRS) campaigns have also been conducted in industrial agricultural estates and in several districts [25, 34]. During the mass ITN distribution campaigns IRS districts do not receive ITNs. Although IRS use has been related to a reduction in parasitaemia in Malawi [35], due to their high cost and the spread of insecticide resistance, it was scaled back to Nkhotakota district in 2012 and implemented only sporadically until 2018. A timeline of the different malaria vector control activities in Malawi can be found in Additional file 1.
The MIS are designed to provide nationally, regionally and urban/rural representative data on 14 core malaria indicators during peak malaria transmission, including malaria prevalence and malaria control [36]. A different combination of households is surveyed each survey year using a two-stage stratified cluster design [37]. Field teams ask the head of household an array of questions. If the target respondent is unavailable, a person aged 15 or over living in each household is asked the questions. The field team also conduct malaria testing and geo-reference the location [38]. Household anonymity is maintained by displacing the coordinates randomly between 0 and 10 km within the second administrative boundary. The locations of the household clusters and their urban/rural classification are visualized in Additional file 2.
For this study, household-level MIS data from 2012, 2014 and 2017 were analysed with ‘household’ as the unit of analysis. The timing of the MIS in combination with the implementation of vector control activities in Malawi are tabulated in Additional file 1. A household is defined as one person or a group of people living together in a housing unit who acknowledge one adult as the head of the household [37]. The 2010 MIS was excluded, as it was not conducted by the DHS program, not geo-referenced and did not include Rapid Diagnostic Test (RDTs) results.
In the MIS surveys, after verbal consent, children aged 6-to-59 months were tested for malaria using rapid diagnostic tests (RDTs) and microscopy of blood smears [38]. The RDT has a sensitivity of 99.5% and specificity of 98% for Plasmodium falciparum [39]. A small drop of blood from the finger or heel of a child was tested for P. falciparum in 2012 and 2014 using the SD Bioline rapid diagnostic test (Standard Diagnostic Inc., Korea) and in 2017 for P. falciparum and other Plasmodium species using the SD Bioline Malaria Ag P.f/Pan rapid diagnostic test (Standard Diagnostic Inc., Korea). The children who tested positive were offered malaria pharmacological treatment according to the current standard procedures for malaria treatment in Malawi [38]. A second blood sample was taken from the children for a thick and thin blood smear to confirm infection. The slides were checked in the laboratory for the presence of Plasmodium parasites by two independent microscopists. Both the RDT and blood smear results are shared in the MIS report. The NMCP in Malawi uses only the blood smear results for malaria prevalence calculations, while the RDT results are used for treatment of cases in the field. Correspondingly, the focus in this paper was on blood smear prevalence results.
The questionnaire respondents are asked to show the surveyor all the nets in the household. Surveyors are instructed to count and inspect these nets. From this ITN data, household ITN ownership (indicator 1: proportion of households with at least 1 ITN), household ITN access (indicator 2: proportion of households with at least 1 ITN for every 2 people sleeping in the house the night before), population ITN access (indicator 3: proportion of the population that could potentially sleep under an ITN the previous night, assuming each net could be used by 1.6 people (as suggested by [40])) and ITN use (indicator 4: proportion of the population that slept under an ITN the previous night) were calculated as proxies for ITN access and use [36]. A detailed description of the calculations is available in Additional file 3. These indicators were used to identify ownership and behavioural gaps. The ownership gap is the gap between household ITN ownership and population ITN access [36]. This gap highlights that although nets are reaching a large proportion of the population, the number of ITNs available in the householdis insufficient to cover everyone in the household (assuming one ITN provides protection for 1.6 people). The behavioural gap is the gap between population ITN access and ITN use [36]. This gap identifies if net access is sufficient to cover all people sleeping in the household (assuming one ITN provides protection for 1.6 people), yet ITNs are not being used.
The data were analysed in the R software [41,42,43]. National averages and 95% confidence intervals were calculated for each variable in each survey separately using a weighting system that adjusts for selection probability differences [37]. Paired t-tests were undertaken to compare differences in malaria prevalence between RDTs and blood smear results. Four different model types were used to investigate the relationship between childhood malaria and ITN use in Malawi (Fig. 1). Two assumptions were made: (i) population ITN access is a good representation of the accessibility of good quality ITNs for the population in the area; and (ii) ITN use the previous night is a good representation of the overall use of ITNs by the population throughout the malaria transmission seasons due to the large sample size of the MIS.
The four different model types used to investigate the relationship between ITN use and childhood malaria in Malawi
An initial investigation was undertaken using the multivariate linear regression model (model 1). As the residuals were not normally distributed, the multivariate linear regression model was extended to a binomial generalized linear model (BGLM, model 2). In this model, variable selection was based on the Akaike information criterion (AIC). The selected predictor covariates consisted of the year of survey, urban/rural areas, population ITN access, ITN use and the interaction between population ITN access and ITN use, with malaria prevalence from blood smear used as the outcome variable. Residuals were mapped and tested for autocorrelation to check if values were independent.
As the BGLM model (model 2) showed autocorrelation in the residuals, data were formally tested for spatial dependence in the variance. Spatial dependence was shown to be significant [see Additional file 4] and, therefore, spatial analysis was undertaken by fitting a Binomial generalized geostatistical linear model (BGGLM) [44] to the aggregated malaria prevalence data (model 3):
where \(\alpha\) is the intersect; \({\beta }_{1}\) the linear regression parameter for the survey year, \({\beta }_{2}\) the parameter for the rural/urban variable at location \(i;\) \({\beta }_{3}\) the parameter for ITN population access; \({\beta }_{4}\) the parameter for ITN use; \({\beta }_{5}\) the parameter for the interaction between ITN population access and ITN use; \(S\left({x}_{i}\right)\) a spatial Gaussian process; \({x}_{i}\) the ith spatial location and \({Z}_{i}\) the residual extra-binomial variation within a sampling location. The Monte Carlo maximum likelihood (MCML) estimation was used to fit the binomial model. The MCML estimation was repeated using the parameter values from the initial MCML to improve the MCML estimates [44].
There is a logical relation between the ITN population access and ITN use data variables, converging at 0 value. Further investigation of these variables indicated weak collinearity (variance inflation factor < 5), and adjustment in the models was not needed. However, as an interaction between both variables is expected, the interaction was left in the model.
As is common practice, a mixed model with a geostatistical spatial random effects term was adopted. This type of mixed model is justified where the residuals from the linear fixed effects term of the regression model are autocorrelated. The difference between urban and rural areas was accounted for by adding it as a covariate. To account for differences in data collection teams, malaria cycle periods and climate between the different survey years, the year was also included as a covariate.
The BGGLM (model 3) is a global model in that its parameters are spatially stationary. To investigate the variation of parameters across space, a non-stationary GLM was fitted (model 4). This non-stationary method of regression considers only a portion of the data around the prediction point at a time, and is sometimes referred to as ‘local mapping.’ A so-called ‘kernel’ is moved across the study area and at each location in a pre-defined grid, the regression coefficients are estimated with the GLM (model 2) as its base, using the data weighted by the kernel. The size of kernel can be adaptive, to ensure that sufficient information is included even in areas where data are sparse. In this way, the non-stationary method allows the estimated regression parameters to be mapped throughout the study area [45]. Specifically, it calculates local regression parameters by fitting a regression model for malaria prevalence \(p\) such that:
where (\({u}_{i}\) ,\({v}_{i})\) represents the coordinates at location \(i\) ; \({\alpha }_{i}\) the intersect at location \(i\) ;\(m\) the number of independent variables; \(\beta {(u}_{i}\) ,\({v}_{i})\) the local regression parameter for the \(k\) th independent variable at location \(i\) ; \({x}_{k,i}\) the \(k\) th covariate; and \({Z}_{i}\) the error at location \(i\) . This expression allows the \(\beta\) to vary with the location coordinates (\({u}_{i}\) ,\({v}_{i})\) , making the model spatially non-stationary; and allows \(\beta\) to be estimated via weighted least squares with the weights matrix obtained from a Gaussian kernel, attributing larger weights to values of predictors from more proximate locations. Analysis was undertaken using the package 'GWmodel' [46]. Standard errors were estimated using the bootstrap function, with coefficients re-estimated 1,000 times at each grid point. The size of the kernel is controlled by the bandwidth, identified using both the cross-validation and AIC corrected (AICc) approach within the package. Collinearity is minimized as described by the R package documentation.
The MIS data used for this research consisted of 10,538 households in 430 different locations collected during three different years (Table 1). Malaria prevalence calculated using RDT results was higher than when calculated using blood smears for 2012 (paired t-test P < 0.001, df = 139), 2014 (P = 0.012, df = 139) and 2017 (P < 0.001, df = 149). A comparison of malaria prevalence for both RDTs and blood smears on the same children showed more than 80% overlap [see Additional file 5].
Although malaria prevalence in the different survey years did not differ (Table 1), the number of positive cases decreased in central Malawi (Fig. 2). Malaria prevalence was lower in urban areas compared to rural areas (both RDT and blood smear P < 0.001) [see Additional file 6].
Malaria prevalence and ITN access summary results for the household clusters surveyed in 2012, 2014 and 2017. A malaria prevalence in children measured using blood smear test in the different household clusters. B Population ITN access in the different household clusters with the assumption that one ITN provides protection for 1.6 people
All ITN indicators suggested an increase in ITN access and use through time (Table 1). By 2017, household ownership had increased from 55 to 82%. Ownership was high throughout the country, with most households owning at least one ITN (Fig. 2) [see Additional file 7]. However, this high household ownership did not translate into similarly high household access. In 2017, only 42% of households had sufficient ITNs to cover all sleeping people and ITN access for the population was below 60% for the majority of areas surveyed. Net coverage in central Malawi was slightly lower than for the rest of the country, with also fewer people sleeping under an ITN [see Additional file 7]. There is evidence of an ownership gap, with ITN campaigns reaching most households across the country, but not supplying sufficient numbers for complete coverage (Fig. 3). There is no evidence of a behavioural gap in 2017, with the increase in population ITN access from 2012 to 2017 aligned with an increase in ITN use.
Scatterplots visualizing the ITN ownership gap and behavioural gap in 2012, 2014 and 2017, with colours indicating the different regions. The black lines are the trend lines representing no gap. A Comparison of household ITN ownership to household ITN access to visualize the ITN ownership gap. B Comparison of population ITN access to ITN use to visualize the behavioural gap
The BGGLM revealed a negative association between malaria prevalence and ITN population access in all regions in Malawi, with low child malaria prevalence occurring in areas with high population access to nets. In contrast, malaria prevalence was positively associated with sleeping under a net, with high use occurring in areas with high child malaria prevalence (Table 2). Both ITN variables showed heterogeneity across the different regions, suggesting a heterogeneous relationship between child malaria prevalence and ITN population access and ITN use. The interaction between ITN population access and ITN use was not significant.
There was a clear relationship between malaria prevalence and the urban/rural variable, with child malaria cases three times higher in rural areas than urban areas. The same model, run independently for each region, indicated that the difference in malaria prevalence between urban and rural areas was demonstrated mainly in central and southern Malawi, where malaria was four times higher in rural areas (Table 2). Child malaria prevalence was lower in 2012 than 2014, but did not differ between 2012 and 2017.
The non-stationary coefficient estimates for both ITN population access and ITN use for the three different years were mapped [see Additional file 8]. The parameters (i.e., relationship) between child malaria prevalence and the indicators exhibit some geographical variation, although the standard errors are large (Table 3). The north of the country shows parameters consistently significant compared to the south of the country. A temporal difference appears between 2014 and 2017, with the correlation between child malaria prevalence and both indicators reduced in 2017 compared to 2014.
This study presents an important, key modelling approach that allows malaria control programmes to spatially unravel the relationship between child malaria prevalence and ITN distribution. The non-stationary generalized linear model helped visualize variation in the relationship between malaria cases and bed net indicators. Increasing household access to and ownership of ITNs resulted in a decrease in malaria cases. However, this did not occur homogeneously across the country. The ITN coverage and use increased from 2012 to 2017, while child malaria prevalence decreased only in some areas. The BGGLM regression using MIS survey data showed that child malaria prevalence had a negative association with ITN population access and a positive association with ITN use. However, it is challenging to identify cause and effect without temporally linked data. The MIS data are currently the best information available to understand the impact of ITNs on malaria prevalence in Malawi. The large variation in odds ratio reflects the high uncertainty in the data and the need for further localized data that takes into account any potential confounder, to spatially disentangle the relationship between ITN indicators and malaria prevalence.
A negative association between malaria prevalence and ITN population access in all regions of Malawi was found. Areas with high population access to nets had lower malaria prevalence than areas with low access. This corresponds to earlier work in Malawi, where ITN ownership was found to be protective against malaria parasitaemia in children [7, 47]. The non-stationary spatial analysis highlighted heterogeneity in the relationship with some areas showing no negative association. The assumption was made that nets maintained physical integrity and bioefficacy of insecticides at least one year post-distribution. Yet, the lack of a discernible relationship could indicate that the roll out was slow, new nets from the campaigns were not used (still in packaging), nets were not used consistently or nets were ineffective. The nets are unlikely to be entirely ineffective, as even with high insecticide resistance the physical barrier of nets and sublethal effect of the insecticides results in a negative relationship between malaria prevalence and ITNs [47,48,49]. It is likely that other factors also affected prevalence, creating artefacts in the relationship between the outcome and predictor. For example, behavioural resistance, with mosquitoes biting when people are not protected by the nets, could play a role [50, 51]. The quality of the nets has also not been considered, while studies have shown that using a mosquito bed net that is more than one year old is a risk factor for malaria in Malawi [52]. Moreover, extreme droughts the year prior could have impacted overall malaria prevalence dynamics as the model does not account for key environmental factors [53].
There was a positive relationship between malaria prevalence and ITN use in all Malawian regions, with high net use occurring in areas with high child malaria prevalence. The question ‘did you sleep under a bed net last night?’ was used as a proxy for overall ITN use. The question focuses on ‘last night’ and does not capture the overall ITN use behaviour. The positive association found here between malaria prevalence and net use is likely due to the large numbers of malaria cases motivating the population to sleep under nets more regularly in the short term. The MIS household surveys are sensitive to recall bias and social desirability factors [54]. Additionally, nightly variables, such as temperature and mosquito density, could also greatly impact the use of nets [55]. For ITN use throughout the malaria season to be related to malaria prevalence, thus identifying its protective ability, temporal studies are necessary where fieldworkers independently confirm bed net use over both space and time throughout the season.
This study revealed spatial heterogeneity in the relationship between malaria prevalence and the ITN indicators, with child malaria prevalence higher in rural than urban areas. It is widely accepted that malaria dynamics differ between rural and urban areas due to differences in demographics, socioeconomics, housing, drainage and access to health care [22, 56,57,58,59]. Generally, in Malawi malaria is more prevalent in rural areas, although high case numbers have also been reported in urban areas [10, 11, 22, 56, 58, 60]. Malaria control likely has a different impact in the urban and rural habitats, due to environmental differences and population movements. Furthermore, the dichotomization of the factors areinconsistent across the research field (the definition of what an urban or rural area is differs), with a more nuanced definition of urban malaria risk and prevention efforts necessary in Malawi to control adequately for the contextual factors that drive malaria prevalence [61]. The results from this study further substantiate the complexity of urban/rural malaria dynamics [22, 56] and show the importance of including this factor in the analysis to investigate the impact of vector control for sustainable intervention measures.
Many studies have shown that the scale-up of ITNs has protected the Malawi population from malaria [6,7,8,9]. This study shows that this relationship varies across space. Heterogeneity in malaria cases has been identified previously in Malawi [14, 52, 62]. Malaria epidemiology is a complex dynamic between many factors at the individual, household and community levels [62]. Confounders can have a large impact on this heterogeneity, and they are challenging to identify. The introduction of vector control in this environment can have very different outcomes, with the decline rate in malaria prevalence known to be very different across the country [14]. One possible explanation for spatial heterogeneity in the relationship between malaria prevalence and ITN use is that it captures the heterogeneity in malaria prevalence, with net coverage homogeneously high throughout the country. A small group of households can account for a majority of cases, which results in spatial heterogeneity in even small geographical areas [63]. Another possibility is that insecticide and behavioural resistance of vector species is impacting net efficacy. Little information about the spatial distribution of malaria vectors and their resistance status is available, but reports do indicate heterogeneity in its spread with high resistance reported in the south and around Lake Malawi [3]. It is also important to note the potential effect of other malaria control interventions on the heterogeneity of the relationship between malaria prevalence and ITNs, including the use of IRS, larval source management and house improvements [34, 64] on malaria prevalence, which have not been considered in this study. Understanding the spatial dynamics will help increase the effectiveness of vector control campaigns, for example, by changing to bed nets that kill resistant mosquitoes more effectively or shifting to different control tools altogether in specific areas. The fine-scale spatial and temporal heterogeneity and their causes need to be investigated further by implementing field studies designed to answer these specific questions.
In concurrence with the recommendations from the Ministry of Health in Malawi, and in alignment with prior Malawian MIS studies [10, 11], blood smear data were used as the sole proxy for malaria prevalence. Both the RDT and blood smear malaria testing methods have strengths and limitations [65]. As shown here, the RDT generally produces higher positivity rates than blood smears, as it measures antigens that are detectable in the blood up to four weeks after parasite clearance [66]. This is contrary to blood smear tests, which measure the physical presence of malaria parasites. From a modelling perspective, instead of choosing one indicator, joint distribution modelling of RDT and blood smear results could help improve inference. Recently, Amoah et al.[67] developed a geostatistical framework to combine spatially referenced disease prevalence data from multiple diagnostics. Joint distribution models draw benefit from the combination of different diagnostics, although particular care needs to be taken for diagnostics with large discrepancies in sensitivity and specificity.
It is important to note that malaria prevalence is not the only way to measure malaria burden in a country. Studies have found mortality reduction without a decrease in malaria prevalence, and the other way around [68]. Furthermore, malaria prevalence here is focused on children aged 6-to-59 months, while school-aged children are both at higher risk of infection and asymptomatic infection [69, 70]. The study could, thus, be improved by using malaria prevalence in the entire population as an outcome variable, which is not possible with MIS data. Excluding this important risk group from the analysis could have skewed the data and results. Additionally, the MIS data are a snapshot and do not include the seasonality of malaria. Although surveys are planned during peak malaria season, as shown by Chirombo et al.[15], this peak differs yearly in Malawi. Whether data are collected during the peak season or two weeks later hugely impacts malaria prevalence estimates. If snapshots of malaria prevalence are compared between the different years without a clear understanding of the seasonal dynamics in these different years, this will greatly influence the analysis and incorrect conclusions can be made. It is important to place malaria prevalence from MIS data in the context of the country and investigate how this is linked to mortality and other malaria indicators.
The MIS data are currently the only data available in Malawi to investigate the relationship between malaria prevalence and ITN use. As the MIS has not been designed for this purpose, caution is advised when interpreting the results. Studies designed to investigate the impact of ITNs are indispensable to understand the efficacy of ITNs [68]. Until these studies are available, creative solutions are necessary to analyse MIS data. Two studies have previously used MIS data from 2012 and 2014 to investigate the relationship between child malaria prevalence and ITN use. Contrary to this study, both found that the number of bed nets per household was not significantly associated with malaria morbidity [10, 11]. Both studies focused on socio-demographic characteristics, while this study included the spatial coordinates of the data within the analysis. This study shows that the impact of ITNs differs geographically. For vector control programmes to make informed decisions about future control activities, it is essential to have access to both country-wide and spatially disaggregated analyses of ITNs impact. This spatial disaggregation is especially important with an increase in IRS activities since 2019 and a combination of PBO and dual active ingredient bed nets distributed during mass ITN campaigns in 2018 and 2021 [33].
The non-stationary generalized linear model visualized the geographically changing relation between malaria prevalence and the ITN indicators. Although the geostatistical mixed regression model presented the overall spatial relationship between child malaria prevalence and bednet indicators, it did not allow for visualization of variation of these relationships across Malawi. The maps produced by the non-stationary model can help vector control programme spatially disentangle the impact of interventions on malaria prevalence for vector control programs. The non-stationary model has some limitations. For example, the bandwidth is optimized based on accurate prediction of the response variable, not on accurate estimation of the coefficients [46]. Especially when the regression model is fitted within a small kernel or with limited data, collinearity can be a problem [71]. Nevertheless, for spatially clustered data such as the MIS data, this method can be appropriate to provide insights into how estimated relations vary across the country. Although relationships are not causal, they could highlight underlying covariates that have been left unmeasured. It is an important tool for eco-epidemiological studies [72,73,74], especially for a disease such as malaria, that is so closely related to the environment and the sociodemographic dynamics of the population. Yet, spatially non-stationary models have only rarely been used in malaria research.
Malaria prevalence in children under five had a negative association with ITN population access and a positive association with ITN use in unadjusted models, showing spatial heterogeneity in these relationships across Malawi. The non-stationary generalized linear model is an important modelling approach that helps vector control programmes visualize variation in the relationship between malaria cases and intervention methods. This study highlights the complexity of the relationship between malaria and ITNs and the clear need for spatially disaggregated data and models to inform localized control activities.
The datasets analysed during the current study are available in the DHS program repository, which can be accessed through (https://dhsprogram.com/).
National Malaria Control Programme K-WTRPaLSoHTM, Malawi Liverpool Wellcome Trust Cliinical Research Programme. Epidemiology and control profile of malaria in Malawi. 2018.
WHO. World Malaria report 2021. Geneva: World Health Organization; 2021.
NMCP. Integrated Vector Control Strategy for Malaria Control Malawi 2020–2024. Lilongwe, 2020.
Wilson AL, Courtenay O, Kelly-Hope LA, Scott TW, Takken W, Torr SJ, et al. The importance of vector control for the control and elimination of vector-borne diseases. PLoS Negl Trop Dis. 2020;14: e0007831.
Population, total - Malawi 2022. Available from: https://data.worldbank.org/indicator/SP.POP.TOTL?locations=MW.
Escamilla V, Alker A, Dandalo L, Juliano JJ, Miller WC, Kamthuza P, et al. Effects of community-level bed net coverage on malaria morbidity in Lilongwe. Malawi Malar J. 2017;6:142.
Hershey CL, Florey LS, Ali D, Bennett A, Luhanga M, Mathanga DP, et al. Malaria control interventions contributed to declines in malaria parasitemia, severe anemia, and all-cause mortality in children less than 5 years of age in Malawi, 2000–2010. Am J Trop Med Hyg. 2017;97:76–88.
Boudová S, Divala T, Mawindo P, Cohee L, Kalilani-Phiri L, Thesing P, et al. The prevalence of malaria at first antenatal visit in Blantyre, Malawi declined following a universal bed net campaign. Malar J. 2015;14:422.
Buchwald AG, Coalson JE, Cohee LM, Walldorf JA, Chimbiya N, Bauleni A, et al. Insecticide-treated net effectiveness at preventing Plasmodium falciparum infection varies by age and season. Malar J. 2017;16:32.
Zamawe CO, Nakamura K, Shibanuma A, Jimba M. The effectiveness of a nationwide universal coverage campaign of insecticide-treated bed nets on childhood malaria in Malawi. Malar J. 2016;15:505.
Zgambo M, Mbakaya BC, Kalembo FW. Prevalence and factors associated with malaria parasitaemia in children under the age of five years in Malawi: a comparison study of the 2012 and 2014 Malaria Indicator Surveys (MISs). PLoS ONE. 2017;12: e0175537.
Nkoka O, Chuang TW, Chen YH. Effects of insecticide-treated net access and use on infant mortality in Malawi: a pooled analysis of demographic health surveys. Prev Med. 2019;127: 105790.
Mathanga DP, Mwandama DA, Bauleni A, Chisaka J, Shah MP, Landman KZ, et al. The effectiveness of long-lasting, insecticide-treated nets in a setting of pyrethroid resistance: a case-control study among febrile children 6 to 59 months of age in Machinga District. Malawi Malar J. 2015;14:457.
Chipeta MG, Giorgi E, Mategula D, Macharia PM, Ligomba C, Munyenyembe A, et al. Geostatistical analysis of Malawi’s changing malaria transmission from 2010 to 2017. Wellcome Open Res. 2019;4:57.
Chirombo J, Ceccato P, Lowe R, Terlouw DJ, Thomson MC, Gumbo A, et al. Childhood malaria case incidence in Malawi between 2004 and 2017: spatio-temporal modelling of climate and non-climate factors. Malar J. 2020;19:5.
Wondji CS, Coleman M, Kleinschmidt I, Mzilahowa T, Irving H, Ndula M, et al. Impact of pyrethroid resistance on operational malaria control in Malawi. Proc Natl Acad Sci USA. 2012;109:19063–70.
Mzilahowa T, Chiumia M, Mbewe RB, Uzalili VT, Luka-Banda M, Kutengule A, et al. Increasing insecticide resistance in Anopheles funestus and Anopheles arabiensis in Malawi, 2011–2015. Malar J. 2016;15:563.
Riveron JM, Chiumia M, Menze BD, Barnes KG, Irving H, Ibrahim SS, et al. Rise of multiple insecticide resistance in Anopheles funestus in Malawi: a major concern for malaria vector control. Malar J. 2015;14:344.
Mwendera CA, de Jager C, Longwe H, Kumwenda S, Hongoro C, Phiri K, et al. Challenges to the implementation of malaria policies in Malawi. BMC Health Serv Res. 2019;19:194.
Mategula D, Ndeketa L, Gichuki J, Zimba B, Ching’ani W, Chipeta MG. Effect of bed net colour and shape preferences on bed net usage: a secondary data analysis of the 2017 Malawi Malaria Indicator Survey. Malar J. 2020;19:428.
Berthe S, Harvey SA, Lynch M, Koenker H, Jumbe V, Kaunda-Khangamwa B, et al. Poverty and food security: drivers of insecticide-treated mosquito net misuse in Malawi. Malar J. 2019;18:320.
Mathanga DP, Tembo AK, Mzilahowa T, Bauleni A, Mtimaukenena K, Taylor TE, et al. Patterns and determinants of malaria risk in urban and peri-urban areas of Blantyre. Malawi Malaria Journal. 2016;15:590.
Nkoka O, Chuang TW, Chuang KY, Chen YH. Factors associated with insecticide-treated net usage among women of childbearing age in Malawi: a multilevel analysis. Malar J. 2018;17:372.
Andronescu LR, Buchwald AG, Coalson JE, Cohee L, Bauleni A, Walldorf JA, et al. Net age, but not integrity, may be associated with decreased protection against Plasmodium falciparum infection in southern Malawi. Malar J. 2019;18:329.
Ministry of Health, NMCP, Measure DHS, ICF, International. Malawi Malaria Indicator Survey (MIS) 2012. Lilongwe, Malawi, and Calverton, USA; 2012.
Ministry of Health, NMCP, Measure DHS, ICF, Malawi Malaria Indicator Survey (MIS) 2014. Lilongwe, Malawi, and Rockville, USA; 2014.
Ministry of Health, NMCP, Measure DHS, ICF, Malawi Malaria Indicator Survey 2017. Lilongwe, Malawi, and Rockville, USA, 2018.
Gleave K, Lissenden N, Richardson M, Ranson H. Piperonyl butoxide (PBO) combined with pyrethroids in long‐lasting insecticidal nets (LLINs) to prevent malaria in Africa. Cochrane Database Syst Rev. 2017;2017:CD012776.
United Nations Population Division. World Population Prospects Revision. 2019.
Bennett A, Kazembe L, Mathanga DP, Kinyoki D, Ali D, Snow RW, et al. Mapping malaria transmission intensity in Malawi, 2000–2010. Am J Trop Med Hyg. 2013;89:840–9.
Mathanga DP, Walker ED, Wilson ML, Ali D, Taylor TE, Laufer MK. Malaria control in Malawi: current status and directions for the future. Acta Trop. 2012;121:212–7.
President’s Malaria Initiative. Malaria Operational Plan FY 2016. USAID, 2016.
Topazian HM, Gumbo A, Brandt K, Kayange M, Smith JS, Edwards JK, et al. Effectiveness of a national mass distribution campaign of long-lasting insecticide-treated nets and indoor residual spraying on clinical malaria in Malawi, 2018–2020. BMJ Glob Health. 2021;6: e005447.
Chanda E, Mzilahowa T, Chipwanya J, Mulenga S, Ali D, Troell P, et al. Preventing malaria transmission by indoor residual spraying in Malawi: grappling with the challenge of uncertain sustainability. Malar J. 2015;14:254.
Skarbinski J, Mwandama D, Wolkon A, Luka M, Jafali J, Smith A, et al. Impact of indoor residual spraying with lambda-cyhalothrin on malaria parasitemia and anemia prevalence among children less than five years of age in an area of intense, year-round transmission in Malawi. Am J Trop Med Hyg. 2012;86:997–1004.
MEASURE Evaluation, The Demographic and Health Surveys Program, President’s Malaria Initiative, Roll Back Malaria Partnership, United Nations Children’s Fund, World Health Organization. Household survey indicators for malaria control. 2018.
Croft TN, Aileen M. J. Marshall, Courtney K. Allen, et al. Guide to DHS Statistics. Rockville, USA: ICF, 2018.
ICF. Demographic and Health Surveys Standard Recode Manual for DHS7. Rockville, USA. DHS Program, 2018.
Tadesse E, Workalemahu B, Shimelis T. Diagnostic performance evaluation of the sd bioline malaria antigen AG PF/PAN TEST (05FK60) in a malaria endemic area of southern Ethiopia. Rev Inst Med Trop Sao Paulo. 2016;58:59–59.
Kilian A, Koenker H, Paintain L. Estimating population access to insecticide-treated nets from administrative data: correction factor is needed. Malar J. 2013;12:259.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. 2019; https://www.R-project.org/.
R Core Team. foreign: Read Data Stored by 'Minitab', 'S', 'SAS', 'SPSS', 'Stata', 'Systat', 'Weka', 'dBase', .... R package. version 08–71. 2018.
Lumley T. Analysis of complex survey samples. J Stat Softw. 2004;9:1–19.
Giorgi E, Diggle PJ. PrevMap: An R Package for Prevalence Mapping. J Stat Softw. 2017;78:1–29.
Brunsdon C, Fotheringham AS, Charlton ME. Geographically weighted regression: a method for exploring spatial nonstationarity. Geogr Anal. 1996;28:281–98.
Gollini I, Lu B, Charlton M, Brunsdon C, Harris P. GWmodel: an R Package for exploring spatial heterogeneity using geographically weighted models. J Stat Softw. 2015;63:1–50.
Lindblade KA, Mwandama D, Mzilahowa T, Steinhardt L, Gimnig J, Shah M, et al. A cohort study of the effectiveness of insecticide-treated bed nets to prevent malaria in an area of moderate pyrethroid resistance. Malawi Malar J. 2015;14:31–31.
Glunt KD, Coetzee M, Huijben S, Koffi AA, Lynch PA, N’Guessan R, et al. Empirical and theoretical investigation into the potential impacts of insecticide resistance on the effectiveness of insecticide-treated bed nets. Evol Appl. 2017;11:431–41.
Churcher TS, Lissenden N, Griffin JT, Worrall E, Ranson H. The impact of pyrethroid resistance on the efficacy and effectiveness of bednets for malaria control in Africa. Elife. 2016;5: e16090.
Kreppel KS, Viana M, Main BJ, Johnson PCD, Govella NJ, Lee Y, et al. Emergence of behavioural avoidance strategies of malaria vectors in areas of high LLIN coverage in Tanzania. Sci Rep. 2020;10:14527–14527.
Carrasco D, Lefèvre T, Moiroux N, Pennetier C, Chandre F, Cohuet A. Behavioural adaptations of mosquito vectors to insecticide control. Curr Opin Insect Sci. 2019;34:48–54.
Hajison PL, Feresu SA, Mwakikunga BW. Malaria in children under-five: a comparison of risk factors in lakeshore and highland areas, Zomba district. Malawi PLoS One. 2018;13:e0207207–e0207207.
Drought-Ravaged Malawi Faces Largest Humanitarian Emergency in its History: Center for Strategic and International Studies; 2016 [Accessed: 2022 16–02]. Available from: https://www.csis.org/analysis/drought-ravaged-malawi-faces-largest-humanitarian-emergency-its-history.
Battle KE, Gumbo A, Hamuza G, Kwizombe C, Banda AT, Chipeta S, et al. Consultative meeting that examined alignment and discrepancies between health facility and household survey data in Malawi. Malar J. 2019;18:411.
Koenker H, Taylor C, Burgert-Brucker CR, Thwing J, Fish T, Kilian A. Quantifying seasonal variation in insecticide-treated net use among those with access. Am J Trop Med Hyg. 2019;101:371–82.
Wilson ML, Krogstad DJ, Arinaitwe E, Arevalo-Herrera M, Chery L, Ferreira MU, et al. Urban malaria: understanding its epidemiology, ecology, and transmission across seven diverse ICEMR network sites. Am J Trop Med Hyg. 2015;93:110–23.
Hennessee I, Chinkhumba J, Briggs-Hagen M, Bauleni A, Shah MP, Chalira A, et al. Household costs among patients hospitalized with malaria: evidence from a national survey in Malawi, 2012. Malar J. 2017;16:395.
Mzilahowa T, Luka-Banda M, Uzalili V, Mathanga DP, Campbell CH Jr, Mukaka M, et al. Risk factors for Anopheles mosquitoes in rural and urban areas of Blantyre District, southern Malawi. Malawi Med J. 2016;28:154–8.
Kazembe LN, Mathanga DP. Estimating risk factors of urban malaria in Blantyre, Malawi: A spatial regression analysis. Asian Pacific J Trop Biomed. 2016;6:376–81.
Lungu EA, Biesma R, Chirwa M, Darker C. Is the urban child health advantage declining in Malawi?: evidence from demographic and health surveys and multiple indicator cluster surveys. J Urban Health. 2019;96:131–43.
Dear NF, Kadangwe C, Mzilahowa T, Bauleni A, Mathanga DP, Duster C, et al. Household-level and surrounding peri-domestic environmental characteristics associated with malaria vectors Anopheles arabiensis and Anopheles funestus along an urban-rural continuum in Blantyre. Malawi Malar J. 2018;17:229.
Chirombo J, Lowe R, Kazembe L. Using structured additive regression models to estimate risk factors of malaria: analysis of 2010 Malawi malaria indicator survey data. PLoS ONE. 2014;9: e101116.
Kabaghe AN, Chipeta MG, Gowelo S, Mburu M, Truwah Z, McCann RS, et al. Fine-scale spatial and temporal variation of clinical malaria incidence and associated factors in children in rural Malawi: a longitudinal study. Parasit Vectors. 2018;11:129.
McCann RS, Kabaghe AN, Moraga P, Gowelo S, Mburu MM, Tizifa T, et al. The effect of community-driven larval source management and house improvement on malaria transmission when added to the standard malaria control strategies in Malawi: a cluster-randomized controlled trial. Malar J. 2021;20:232.
Azikiwe CC, Ifezulike CC, Siminialayi IM, Amazu LU, Enye JC, Nwakwunite OE. A comparative laboratory diagnosis of malaria: microscopy versus rapid diagnostic test kits. Asian Pac J Trop Biomed. 2012;2:307–10.
Moody A. Rapid diagnostic tests for malaria parasites. Clin Microbiol Rev. 2002;15:66–78.
Amoah B, Diggle PJ, Giorgi E. A geostatistical framework for combining spatially referenced disease prevalence data from multiple diagnostics. Biometrics. 2020;76:158–70.
Nkumama IN, O’Meara WP, Osier FHA. Changes in malaria epidemiology in Africa and new challenges for elimination. Trends Parasitol. 2017;33:128–40.
Walldorf JA, Cohee LM, Coalson JE, Bauleni A, Nkanaunena K, Kapito-Tembo A, et al. School-age children are a reservoir of malaria infection in Malawi. PLoS ONE. 2015;10: e0134061.
Topazian HM, Gumbo A, Puerto-Meredith S, Njiko R, Mwanza A, Kayange M, et al. Asymptomatic Plasmodium falciparum malaria prevalence among adolescents and adults in Malawi, 2015–2016. Sci Rep. 2020;10:18740.
Wheeler D, Tiefelsdorf M. Multicollinearity and correlation among local regression coefficients in geographically weighted regression. J Geogr Syst. 2005;7:161–87.
Acharya BK, Cao C, Lakes T, Chen W, Naeem S, Pandit S. Modeling the spatially varying risk factors of dengue fever in Jhapa district, Nepal, using the semi-parametric geographically weighted regression model. Int J Biometeorol. 2018;62:1973–2016.
Mayfield HJ, Lowry JH, Watson CH, Kama M, Nilles EJ, Lau CL. Use of geographically weighted logistic regression to quantify spatial variation in the environmental and sociodemographic drivers of leptospirosis in Fiji: a modelling study. Lancet Planet Health. 2018;2:e223–32.
Lin CH, Wen TH. Using geographically weighted regression (GWR) to explore spatial varying relationships of immature mosquitoes and human densities with the incidence of dengue. Int J Environ Res Public Health. 2011;8:2798–815.
JT is funded by the MRC Skills Development Fellowship. LS is supported by the Wellcome Trust Seed Award (212501/Z/18/Z), Wellcome Trust NIHR–Wellcome Partnership for Global Health Research Collaborative Award, CEASE (220870/Z/20/Z), and by the Academy of Medical Sciences GCRF Networking Grant Scheme (GCRFNGR7/1329).
Vector Biology Department, Liverpool School of Tropical Medicine, Liverpool, UK
Julie-Anne A. Tangena & Donnie Mategula
Malawi-Liverpool Wellcome Trust Clinical Research Programme, Blantyre, Malawi
Lancaster Ecology and Epidemiology Group, Lancaster University, Lancaster, UK
Lancaster Environment Centre, Lancaster University, Bailrigg, Lancaster, LA1 4YR, UK
Geography and Environmental Science, University of Southampton, Highfield, Southampton, SO17 1BJ, UK
Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, 11A Datun Road, Beijing, 100101, China
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
JT and PA designed the study. JT conducted the analysis under supervision of PA and LS. During the analysis, valuable feedback was provided by DM. Results were interpreted by all authors and written up by JT. All authors read, edited and approved the final manuscript.
Correspondence to Julie-Anne A. Tangena.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Timeline of malaria indicator surveys (MIS) in combination with the main vector control activities implemented in Malawi from 2011 to 2017.
Household cluster location for each Malaria Indicator Survey (MIS) round.
Description of the insecticide-treated net indicator calculations.
Variogram showing spatial dependence in the variance.
Comparison of malaria prevalence calculated using RDT and blood smear results.
Blood smear malaria prevalence in urban and rural areas for 2012, 2014 and 2017 with median.
Insecticide treated net indicators mapped for 2012, 2014 and 2017.
The geographically weighted regression analysis results with β the estimated coefficient range and the standard error (se).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Tangena, JA.A., Mategula, D., Sedda, L. et al. Unravelling the impact of insecticide-treated bed nets on childhood malaria in Malawi. Malar J 22, 16 (2023). https://doi.org/10.1186/s12936-023-04448-y
DOI: https://doi.org/10.1186/s12936-023-04448-y
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By using this website, you agree to our Terms and Conditions, California Privacy Statement, Privacy statement and Cookies policy. Manage cookies/Do not sell my data we use in the preference centre.
© 2023 BioMed Central Ltd unless otherwise stated. Part of Springer Nature.