Next Article in Journal
Bayesian Network-Based Risk Assessment of Single-Phase Grounding Accidents of Power Transmission Lines
Next Article in Special Issue
Statistical Modelling of the Effects of Weather Factors on Malaria Occurrence in Abuja, Nigeria
Previous Article in Journal
Smoking Prevalence and Secondhand Smoke Exposure during Pregnancy and Postpartum—Establishing Risks to Health and Human Rights before Developing a Tailored Programme for Smoking Cessation
Previous Article in Special Issue
Effect of Free Healthcare Policy for Children under Five Years Old on the Incidence of Reported Malaria Cases in Burkina Faso by Bayesian Modelling: “Not only the Ears but also the Head of the Hippopotamus”
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Malaria Case Fatality Rate among Children under Five in Burkina Faso: An Assessment of the Spatiotemporal Trends Following the Implementation of Control Programs

by
Mady Ouédraogo
1,2,3,
David Tiga Kangoye
4,
Sékou Samadoulougou
5,6,
Toussaint Rouamba
1,7,
Philippe Donnen
1,8 and
Fati Kirakoya-Samadoulougou
1,*
1
Centre de Recherche en Epidémiologie, Biostatistiques et Recherche Clinique, Ecole de Santé Publique, Université Libre de Bruxelles, 1070 Brussels, Belgium
2
Institut de Recherche Santé et Sociétés, Faculté de Santé Publique, Université catholique de Louvain, 1200 Brussels, Belgium
3
Institut National de la Statistique et de la Démographie [INSD], 01 BP 374 Ouagadougou 01, Burkina Faso
4
Centre National de Recherche et de Formation sur le Paludisme [CNRFP], 01 BP 2208 Ouagadougou 101, Burkina Faso
5
Evaluation Platform on Obesity Prevention, Quebec Heart and Lung Institute, Quebec, QC G1V 4G5, Canada
6
Centre for Research on Planning and Development (CRAD), Université Laval, Quebec, QC G1V 0A6, Canada
7
Unité de Recherche Clinique de Nanoro, Institut de Recherche en Sciences de la Santé, Centre National de la Recherche Scientifique et Technologique, 42 Avenue Kumda-Yonre, Ouagadougou, Kadiogo 11 BP 218 Ouagadougou CMS 11, Burkina Faso
8
Centre de Recherche en Politiques et systèmes de santé-Santé internationale, École de Santé Publique Université Libre de Bruxelles, 1070 Brussels, Belgium
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(6), 1840; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17061840
Submission received: 30 January 2020 / Revised: 4 March 2020 / Accepted: 9 March 2020 / Published: 12 March 2020
(This article belongs to the Special Issue Geo-Epidemiology of Malaria)

Abstract

:
Reducing the 2015 level of malaria mortality by 90% by 2030 is a goal set by the World Health Organization (WHO). In Burkina Faso, several malaria control programs proven to be effective were implemented over the last decade. In parallel, the progressive strengthening of the health surveillance system is generating valuable data, which represents a great opportunity for analyzing the trends in malaria burden and assessing the effect of these control programs. Complementary programs were rolled out at different time points and paces, and the present work aims at investigating both the spatial and temporal pattern of malaria case fatality rate (mCFR) by considering the effect of combining specific and unspecific malaria control programs. To this end, data on severe malaria cases and malaria deaths, aggregated at health district level between January 2013 and December 2018, were extracted from the national health data repository (ENDOS-BF). A Bayesian spatiotemporal zero-inflated Poisson model was fitted to quantify the strength of the association of malaria control programs with monthly mCFR trends at health district level. The model was adjusted for contextual variables. We found that monthly mCFR decreased from 2.0 (95% IC 1.9–2.1%) to 0.9 (95% IC 0.8–1.0%) deaths for 100 severe malaria cases in 2013 and 2018, respectively. Health districts with high mCFR were identified in the northern, northwestern and southwestern parts of the country. The availability of malaria rapid diagnosis tests (IRR: 0.54; CrI: 0.47, 0.62) and treatment (IRR: 0.50; CrI: 0.41, 0.61) were significantly associated with a reduction in the mCFR. The risk of dying from malaria was lower in the period after the free healthcare policy compared with the period before (IRR: 0.47; CrI: 0.38, 0.58). Our findings highlighted locations that are most in need of targeted interventions and the necessity to sustain and strengthen the launched health programs to further reduce the malaria deaths in Burkina Faso.

1. Introduction

Malaria remains a major public health issue in sub-Saharan Africa (SSA) and considerably contributes to child morbidity and mortality [1]. In 2018, SSA accounted for 94% of world malaria deaths, of which 67% were children under five [1]. Malaria contributed 5% and 18% of under five deaths worldwide and in SSA, respectively [2]. There is, however, some noteworthy progress that has been made. From 2010 to 2018, the malaria mortality in children under five has approximately halved in the world, and this drop in mortality is attributable to intensified control measures including vector control and drug-based interventions. This drop in mortality was even most marked in SSA [1]. In most SSA countries, this progress is, nonetheless, hampered by the challenge of achieving an early diagnostic and appropriate treatment of malaria cases to prevent the occurrence of complications and death [3].
Burkina Faso is one of the eleven highest-burden countries in the world with a 6% contribution to the 94% malaria deaths recorded in SSA in 2018 [1]. According to the statistics of the Ministry of Health, the same year, malaria was responsible for 53% of the morbidity and 66% of the deaths in children under five seen at health facilities [4]. Like other SSA countries, Burkina Faso implemented several control strategies (preventive and curative) to curb the burden of malaria as the WHO recommendations were being issued. The intermittent preventive treatment in pregnant women (with sulfadoxine–pyrimethamine in replacement of chloroquine) was implemented as of 2005. The same year, the use of artemisinin-based combination therapies (ACTs) for uncomplicated malaria was adopted in lieu of monotherapies. Four years later, the systematic rapid diagnostic testing of suspected malaria cases was implemented to promote a rational prescription of ACT and improve the quality of surveillance data. As of 2010, mass distributions campaigns of free-of-charge long-lasting insecticide-treated bed nets were conducted every three year to achieve universal coverage, the latest having been delivered in 2019. The use of injectable artemisinin derivatives, as alternatives to injectable quinine, for severe malaria treatment was introduced in 2012. The seasonal malaria chemoprophylaxis in children under five was gradually deployed from 2014 to 2019 to achieve full coverage. The indoor residual spraying with long-lasting insecticides, interrupted in 2013, was resumed in 2017 in three health districts out of 70. The IRS was actually being previously conducted in only one health district for three years [5]. Vector control and drug-based preventive measures were provided free of charge contrary to clinical malaria control measures, and this was a potential limitation in the mCFR reduction endeavor. In addition to these measures the activities of community-based health workers were intensified in 2016. The negative effects of the “pay and use health services” policy, in use in sub-Saharan African countries, on populations access to healthcare have been consistently documented [6,7]. As a consequence, the debate on the financing of healthcare has shifted from the “pay and use” to the “free of charge” service approach in favor of children under five and pregnant women who are the most vulnerable populations to malaria [8,9,10,11,12,13,14,15,16]. As of 2016, a free healthcare policy (FHC) for children under five and pregnant women was then launched [17]. Although not specific to malaria, the FHC policy is expected to serve as a catalyst especially for drug-based malaria control interventions by increasing the use of health services. Indeed, quasi-experimental assessments conducted in some sub-Saharan countries found that the free-of-charge healthcare (FHC) policy was effective in increasing access to healthcare especially for the most vulnerable groups. The launch of this policy was the starting point for a free-of-charge care for malaria cases in Burkina Faso.
The analysis of the effect of these overlapping control programs on the trends in malaria indicators is based on routine surveillance data, and for malaria case fatality rate (mCFR), their current exploitation among children under five is limited to the estimation of an unadjusted mCFR [18,19]. The robustness of this kind of analysis is limited since the data used have features that makes classic statistical methods inappropriate. Indeed, classic statistical methods do not allow reliable adjustment for intra-national heterogeneity of malaria case fatality rate. Recently developed and validated hierarchical Bayesian spatiotemporal models were implemented on malaria surveillance data and produced robust estimates [20,21,22,23,24,25,26,27,28]. These models proved useful in addressing complex issues in the data structure.
The present work aims at investigating both the spatial and temporal pattern of malaria case fatality (mCFR) by considering the effect of combining specific and unspecific malaria control strategies.

2. Materials and Methods

2.1. Study Settings

Burkina Faso is a 272,967 km2 landlocked country located in the middle of the West African region and was home to an estimated population of 20 870 060 inhabitants in 2019 [29]. The health system has three levels [4]. On the administrative viewpoint, from top to bottom, we have the central level (global data processing, decision-making), the intermediate level (regional coordination and supervision) with 13 regional health directorates and the district level (field operations) with 70 health districts [4]. From the technical viewpoint, the bottom has two layers consisting of 1896 public primary health facilities that refer to 45 district hospitals with surgical capacity. Those district hospitals, in turn, refer to the eight regional hospitals, which themselves, refer to the six teaching hospitals, which are the third level [4].
The climate is tropical with a single four month-long rainy season from June to September. The health districts located in the southwestern part of the country have a higher annual rainfall (above 950mm) compared with those located in the northern part (599–715 mm) [30]. The transmission of malaria is generally stable and occurs throughout the year with a peak in the rains. Plasmodium falciparum is responsible of more than 90% of the infections [31]. The percentage in 2014 of the population living under the poverty line was estimated at 40.1%, with regional variations from 9.6% in the Centre Region to 70.4% in the North Region [32]. This rampant poverty exists alongside a low-performance health system with high mortality rates (129 for 1000 children in 2010) [33]. In addition, the proportion of caregivers promptly seeking (within 24 h of onset) healthcare for feverish children is still low in the Sahel (28%), Boucle du Mouhoun (31%) and Southwest regions (42.8%) [34]. The performance rate of the malaria rapid diagnostic test (mRDT) is also low in the Boucle du Mouhoun (30.3%), Sahel (40.1%) and Hauts Bassins (37.9%) regions [34]. In 2017, the percentage of uncomplicated malaria cases treated with ACT among feverish children for whom care was sought at public primary health facilities was 79% overall [34].

2.2. Data Sources

2.2.1. Malaria Case Fatality Data

The monthly count data on severe malaria cases (confirmed and presumed cases) and malaria deaths in children under five from January 2013 to December 2018, aggregated at district level, were extracted from the national health statistics repository (ENDOS-BF), a DHIS2 (developed by Health Information Systems Programme, University of Oslo, Norvège) web-based platform operational as of 2013 [35]. This data repository has substantially improved the quality of the data by improving its timeliness and completeness. The data discrepancy index between the source and the reported data improved from 61% in 2013 to 80% in 2018 [36].

2.2.2. Malaria Control Programs Data

To assess the impact of the FHC policy [17] on mCFR in children under five, the period from 2013 to 2018 was considered and split into two time intervals: before FHC policy (from 2013 to 2015) and after FHC policy (from 2016 to 2018). This exposure variable was defined as a binary one: 0 for “before FHP” and 1 for “after FHP”.
The data on the availability of mRDT and ACT at health facilities were extracted from the results of the Service Availability and Readiness Assessment (SARA) surveys conducted in 2012, 2014 and 2016 [37]. These SARA surveys, sponsored by the Ministry of Health, are cross-sectional surveys that collected data on a random sample of at least 687 health facilities each year [38]. The health facilities were selected using a stratified two-stage cluster sampling method based on the administrative division of the country in regions.
The data on malaria case management with ACT were extracted from the results of the Malaria Indicators Surveys (MIS) conducted in 2014 and 2017. The MIS, sponsored by the National Statistics and Demography Institute (INSD) in collaboration with a global consulting company (ICF International), are cross-sectional surveys that collected data on a random sample of 6552 and 6332 households in 2014 and 2017, respectively [34,37]. The households were selected using a stratified two-stage cluster sampling method based on the admirative division of the country in regions.

2.2.3. Control Variables

Acute malnutrition is a debilitating condition that impairs the immune response to pathogens [39], including malaria parasites and, even more, the response to some antimalarial drugs [40]. It has been linked to an increased risk of childhood death [41]. At population level, the burden of acute undernutrition in children under five is best measured by the prevalence of global acute malnutrition (GAM), which informs on their short-term nutritional history as well the general population nutritional status. The GAM index is classified into five categories from very low to very high, and an emergency threshold was set at 10% to indicate public health immediate actions exigency [42]. We have then adjusted our model with a dichotomous variable for GAM, based on the 10% UNHCR standard threshold (0 if GAM < 10%; 1 if GAM ≥ 10%).
Two additional variables, the distance to the nearest health facility (binary: 0 for “< 5 km”, 1 for “≥5 km”) and the number of nurses for 5000 inhabitants (binary: 0 for “< 5000”; 1 for “≥5000”), were derived from health indicators (mean theoretical health facility operational range and number of inhabitants for 1 nurse, respectively) in the health statistics annual yearbooks of 2013 to 2018 [4,18,19,43,44,45]. These two variables were recoded according to the WHO standards [46].
The geographic location of the health districts and the year were used to merge the different datasets.

2.3. Statistical Analysis of Malaria CFR

2.3.1. Estimating District-Level Interventions Coverage

The data on the confounders of interest (mRDT and ACT availability, ACT use), obtained from the MIS and SARA surveys, were collected from a random sample of health districts to derive estimates at regional and national levels. These data were, therefore, not appropriate to compute reliable estimates for each of the 70 health districts. To address this issue, we fitted a Bayesian binomial model to the aggregated data on these confounders to estimate them at the district level. The details on the modelling of the estimates at the district level are provided in the Appendix A1.

2.3.2. Descriptive Analysis of Temporal Trends

The malaria case (severe) fatality rate was estimated by dividing the monthly total number of malaria deaths in children under five by the monthly total number of severe malaria cases in children under five. This estimation was done for each health district [47]. A time plot was drawn to describe the temporal evolution of the crude mCFR in Burkina Faso on a monthly basis between January 2013 and December 2018.

2.3.3. Spatiotemporal Modelling

The occurrence of malaria death in children under five on a monthly basis at the health district level was a rare event. Moreover, we fit spatio-temporal modeling with zero-inflated Poisson (ZIP) and negative binomial (ZNIB) distribution and used the deviance criterion (DIC) to determine the best spatio-temporal model. Results (Appendix A2) indicate a Bayesian ZIP regression model fit the data to describe the spatiotemporal evolution of mCFR in children under five living in Burkina Faso [48]. If y it and n it represent the number of malaria deaths and the number of severe malaria cases in children under five, respectively, in the health district i = (1, 2, …, 70) at time t = (1, 2, …, 12), the malaria deaths can be modelled as follows:
y i t Z I P ( n i t µ i t )
The risk of death µ i t can be modelled as follows:
log ( μ it ) = + s i + v t + X T β + δ it + ε it
where:
β is the regression coefficient,
s i is the spatial variance,
v t is the temporal variance,
δ i t space-time interaction.
A non-informative Gaussian prior distribution is attributed to the vector of the regression coefficient β (β N(0, 1000)) [48].
The term δ i t allows for each health district to have its own trend. This last term indicates the additional variability in the data, not explained by the other components of the model. For the spatial random effect term s i , we assign the Besag–York–Mollié (BYM) model [49]. This spatial variance was decomposed into the sum of a structured spatial random effect ( u i ) and an unstructured spatial random effect ( ζ i ) . The combination u i + ζ i accounts for the spatial dependency in the modelling process. ζ i follows a Gaussian distribution, and u i is modelled via a conditional autoregressive regression model (CAR) with a neighbor matrix W of size N × N, where the diagonal component w i i = 0 and the non-diagonal component w i j = 1 if districts i and j share a common border, and w i j = 0 otherwise. The prior posterior CAR distribution, defined for the spatial random effect, assumes that adjacent districts tend to have similar risks of mCFR in children under five. The temporal variance v t was decomposed into a sum of a structured temporal random effect ( γ t ) modelled by a AR   ( 1 ) ,   that   is γ 1 ~ N ( 0 , σ 2 1 ρ 2   )   a n d   γ l ~ N ( ρ γ l 1 , σ 2   ) ,   l = 2 , , 67 , with the autocorrelation parameter ρ quantifying the degree of dependency between three consecutive months, and an unstructured temporal random effect ( ϕ t ) modelled by a Gaussian ϕ t N ( 0 , 1 τ ϕ ) with an overdispersion parameter ϵ i t N (   σ ϵ 2   ) . As suggested by Gelman et al., we assume that all the standard deviations of random effects 1 τ ϕ   a n d   σ ε follow a strictly positive prior semi-Gaussian distribution N + ( 0 , 10 ) [50].
The Integrated Nested Laplace Approximation (INLA) package in R was used to implement the models, which yielded an a posteriori distribution of the marginal effects of the model parameters and the interventions effects. These parameters were summarized using the case fatality rate risk (CFR) and the 95% credible interval. The Richardson classification [51] was used to estimate the excess risk and classify the health districts into three categories according to the national malaria control threshold in year 2018 (1.01%) as follows: (a) high risk of mCFR if the posterior probability distribution for a child to die from severe malaria is above 0.8, (b) low risk if the probability is below 0.2 and (c) moderate risk if the probability is between 0.2 and 0.8. The results (risk of mCFR in children under five) were exported to and mapped with QGIS version 2.18.14 (developed by Open Source Geospatial Foundation, Chicago, IL, USA) [52]. The deviance information criterion (DIC) was used to compare the null model with the full model, and the model with the smaller DIC (best fitting model) was kept.

2.3.4. Ethical Considerations

The present study used data collected from public health facilities and recorded in the national health data repository (ENDOS-BF), which was accessed through the directorate in charge of the national health statistics management (DGESS) with the authorisation of the Ministry of Health. An ethical approval was not requested because the data used were aggregated at district level with no personal information accessible.

3. Results

3.1. Summary of Crude Estimates of Severe Malaria Case Fatality Rate in Children under Five

There was a decreasing trend in mCFR between 2013 and 2018. Indeed, the mCFR went from 2.0 (95% IC 1.9–2.1%) to 0.9 (95% IC 0.8–1.0%) deaths for 100 severe malaria cases in 2013 and 2018, respectively, that being a more than 55% drop (Figure 1). However, it was noticeable that the decrease appeared to stall from September 2016 onwards.

3.2. Spatio-Temporal Dynamics of Crude and Predicted Rate of Malaria Case Fatality Rate by Quarter and by District

Figure A1 (in additional file) shows the spatiotemporal evolution by quarter of the crude mCFR over the period 2013–2018, whereas Figure 2 shows the evolution of adjusted posterior mCFR. These figures show a drop in the mCFR over the period. Thus, although variable between quarters, the number of health districts with a fatality rate of malaria cases above the threshold of the national malaria control program (1.01% in 2018) fell from 40 in the fourth quarter of 2013 to 23 in fourth quarter of 2018. The adjusted fatality rate indicated the same trend with 58 health districts with a mCFR above the threshold of the national malaria control program (1.01% in 2018) in the fourth quarter of 2013 and 27 in the fourth quarter of 2018. Before July 2016, the majority of health districts had an mCFR of between 1 and 2 deaths per 100 severe malaria cases, whereas between July 2016 to December 2018, the majority of health districts exhibited an mCFR less than 1 deaths per 100 severe malaria cases. The health districts in the western (Sindou, Boromo) and northern (Gorom-Gorom) part of the country exhibited higher values with an mCFR of more than 4% between 2013 and 2016.

3.3. Spatial Pattern of Risk of Death Due to Severe Malaria among Children under Five

Figure 3 presents the probabilities of excess risk of death due to severe malaria based on the NMCP threshold for year 2018, which was estimated at 1.01%. The map indicated a heterogeneity in the risk between health districts and displayed the relative location of 28 identified high risk health districts. The high-risk (probabilities of excess risk) areas were predominantly in the health districts located in the northern, northwestern and southwestern parts of the country.

3.4. Association between Health Programs and Malaria Case Fatality Rate

Table 1 shows that the BYM full model with a temporally and spatially structured component (AR1) was more suitable than the null model to describe the spatiotemporal distribution of mCFR in children under five living in Burkina Faso. Indeed, the DIC of the full model (DIC = 17,888.29) was 108.80 points lower than that of the null model (DIC = 17,997.09). From the full model, the overall temporal variance (other the entire period) of mCFR was estimated at 1.209 (95% BCI 1.063–1.955), and the spatial variance was 1.028 (95% BCI 1.027–1.029). The space–time interaction effect was significant at 1.170 (95% BCI 1.169–1.172), indicating a heterogeneity in mCFR across the health districts over the study period.
The results of the multivariable analysis showed that a reduced risk of severe malaria death at health district level was significantly associated with the FHC policy for children under five, with the availability of mRDT and ACT and with the use of ACT for malaria treatment at health facilities. Indeed, from our findings the probability for a child under five to die from severe malaria at health facilities decreased by 50% (95% BCI 40–58%) between 2013 and 2018. In addition, a 100% increase of the proportion of health facilities with mRDT in stock was associated with a 46% (95% BCI 38–53%) decrease of mCFR, and a 100% increase in the proportion of malaria cases treated with ACT was associated with a 46% (95% BCI 38–53%) decrease of malaria case (severe) fatality rate in children under five. We also found that a 100% increase of the proportion of health facilities with ACT in stock was associated with a 50% (95% BCI 39–59%) decrease of mCFR. The distance to the nearest health facility was also significantly associated with a decrease of mCFR in children under five. However, the global acute malnutrition index and the ratio nurse/population were not significantly associated with mCFR in children under five.

4. Discussion

The present study is, to the best of our knowledge, the first attempt to use routinely collected malaria surveillance data to assess the contribution of control programs to mCFR changes in children under five at country level, with intra-national space–time details. Data on severe malaria cases and malaria deaths were extracted from the national health data repository (ENDOS-BF). Associations between the combined health programs and mCFR have rarely been studied previously using routinely collected data. Most studies in this area have rather focused on the effect of LLIN, the “test and treat” policy and the use of ACT for malaria case management on the morbidity/mortality of malaria [23,53,54]. Our results show a significant drop in mCFR of more than 47% in children under five between 2013 and 2018. Such a drop was previously reported in a recent study in which the aim was to map Plasmodium falciparum mortality in Africa between 1990 and 2015 [55]. In our study, we found that this drop was significantly associated with the launched health programs for children under five (Table 1).
It is important to note that before the implementation of FHC, several non-specific (results-based financing project, the health promotion project) or malaria-specific (seasonal chemoprevention, indoor residual spraying) control programs have been implemented for the improvement of health indicators in some health districts. These strategies have certainly contributed to the reduction of mortality among the population of children under five years of age. With the FHC policy, the test and treat components of WHO’s T3 strategy [56] are more optimal. Indeed, there was an improvement in the rate of attendance at health facilities [34] and the rate of performance of rapid diagnostic tests, enabling the detection of true malaria cases [26]. Thus, prompt diagnosis of true cases followed by effective treatment may reduce the likelihood of the progression of uncomplicated malaria cases to the severe forms responsible of case fatalities. In addition, compliance with the national malaria control program guidelines regarding the use of ACT for malaria treatment could contribute to explain this result [54].
The FHC policy may decrease mCFR by enhancing the other health interventions assessed through an increase in the use of health services, the accessibility to an early diagnostic and treatment of malaria [56,57,58,59], given that the caregivers (generally the child’s parent) are able to identify the signs of malaria in children and bring them to the nearest health facility within the 24 h of illness onset [6]. Our findings are similar to those of Posnar et al., who showed that an increased use of health services decreased malaria mortality in children under five [16]. However, studies conducted in Ghana and Cameroun did not find any significant reduction in malaria mortality associated to the FHC policy [60,61]. The authors explained their findings by the under-reporting of malaria deaths combined with a poor source document maintenance before the implementation of the FHC policy [61]. Indeed, registers were not available in many health facilities in their study. In Burkina Faso, though malaria deaths may be underestimated due to treatment-seeking behavior factors [34], underreporting is likely limited. In addition to the health facility monthly activity reports conveyed to the health statistics unit of the ministry of health, a mandatory weekly “telegram” reporting cases of core diseases of public health importance under surveillance, including malaria, is sent to the surveillance unit of the same ministry. In the end, the figures of these parallel reports must match each other. In the same study conducted in Ghana and Cameroun, the authors indicated, in addition, that an increased use of health services could potentially decrease the quality of healthcare (e.g., increased waiting time). The FHC policy should, therefore, be combined with measures to sustain healthcare quality at the risk of generating a reverse effect.
We also found that a reduction of the distance to the nearest health facility would decrease mCFR in health districts as evidenced by other studies suggesting that the shorter the distance to the health facility, the earlier malaria can be diagnosed and adequately treated [6,21].
We have also found that the effect of malaria control programs was heterogeneous across health districts. This heterogeneity could be explained by a number of logistical, medicotechnical and individual factors, varying by health district and with time. National surveys have revealed that around 21% of children who had a fever in the two weeks preceding the survey were not treated promptly [34]. This lack or delay of treatment can lead to complications resulting in death [62]. Indeed, it has been reported that in severe malaria, neurological manifestations appear within the first 12 h of admission [63] and is fatal without treatment [64]. Severe malaria treatment implies, in addition, supportive care to control complications of especially poor prognosis factors such as severe anemia, respiratory distress, bacterial co-infections and acute renal failure [65,66,67,68,69,70]. Blood, oxygen, dialysis and mechanical ventilation [71] may be necessary for an appropriate treatment of severe malaria. The health districts are uneven vis-à-vis the availability of these means for supportive treatment. Finally, other poor prognostic factors such as acute malnutrition [62,72] cannot be suppressed in an emergency. In addition, the disparities between health districts observed in the present study could be partly explained by other health programs. Indeed, the community-based health workers program was implemented in 2015 countrywide. Community-based interventions were reported to increase the coverage of most malaria control interventions [73].
Most studies assessing the impact of combined effect of health interventions (FHC, “test and treat” and use of ACT for malaria case management) on mCFR in children under five do not examine the spatiotemporal variability of the risk at the intra-national level [16,23,53,54,60,61]. It is often postulated that the efficacy of national malaria control programs relies essentially on the availability of reliable local health statistics to guide the decision-making process [74]. Our results confirm this assertion and have implications for the control of mCFR in children under five. Indeed, the measurement of the association with the combined effect of health interventions (FHC, “test and treat” and use of ACT for malaria case management), the identification of health districts with intermediate to excess risk and CFR could guide the planning of control interventions by the NMCP. In Burkina Faso, health districts are the operational level for the implementation of health policies, and one of their missions is to implement malaria control interventions [75]. Such interventions include ensuring availability of mRDT at health facilities and at the community level, the availability of ACT and injectable artesunate for the treatment of uncomplicated and severe malaria, the provision of health facilities with LLIN for free distribution to target populations as well as population information, education and communication (IEC) for an optimal use of these LLIN [53,56,57,58]. In the short term, the NMCP could reinforce the IEC for malaria control at the community level in health districts with moderate to high risk. Such a targeted and timely response from the NMCP relies on its capacity to collect, analyze data and generate situation maps in real time.
In the present study, we did not find an association between global acute undernutrition and mCFR. This result is in keeping with the findings of some prior studies [76,77,78,79]. The association between nutritional status and malaria, in the current state of evidence, remains quite fuzzy. Clinical malnutrition is described as a heterogeneous group of nutritional disorders, and its interplay with immunity and infection is complex [80]. Impaired immunity is likely a consequence as well as a cause of malnutrition [81]. While an increased risk of infection can be observed in some nutrient deficiencies, a protective effect is observed in others and, at last, an increased risk or protective effect can be observed in some other deficiencies depending on the context [82]. Earlier studies reported or suggested a protective effect of undernutrition against malaria, especially refeeding in famine context [83,84,85,86], or reported no association between malaria and undernutrition. However, there is now a strong body of evidence and general agreement that undernutrition, though it may not have a strong effect on malaria infection or incidence, can increase the illness severity and worsen severe malaria outcome [87]. The risk of death was associated with severe acute undernutrition in a number of studies including hospital-based [88,89] ones and randomized trials [90]. Some of these studies did not distinguish between wasting and stunting [91]. In a large multisite, multi-country hospital-based study, the z-score was linked to early, intermediate and late malaria mortality globally, and severe acute malnutrition was associated with mortality at six out of eight sites [92]. The wide array and disparities in the methods used limited the feasibility of conclusive meta-analyses.
The present study has limits. First, the number of severe malaria cases and malaria deaths could have been underestimated because of the low use of health services. In addition, the CAR models we used are subject to estimation biases due to ecological error [93]. However, the models we used took into account the uncertainty in the estimates and the biases related to the data sources allowing us to generate robust estimates.

5. Conclusions

The present study assessed the association between malaria control programs (free-of-charge healthcare, test and treat, use of ACT) and mCFR in children under five living in Burkina Faso. Our findings showed a decline in mCFR and an association with control interventions. Moreover, our results showed geographical differences in the mCFR decline over time in Burkina Faso. Malaria case fatality rate tended to increase in some health districts. If this tendency is maintained in the same localities in subsequent evaluations, they would then be worthy of thorough investigation to identify and address the factors impeding the policy. Moreover, in a resource-limited context, our findings could guide the NMCP in hierarchizing the health districts for the implementation of appropriate control interventions. In addition, our findings highlighted the necessity to sustain and strengthen the health programs in effect to further reduce malaria deaths in Burkina Faso. Finally, the model could be useful to other countries for the assessment of the effect of malaria control interventions and to Burkina Faso for periodic re-evaluation of these programs.

6. Research Data and Supplementary Materials

The dataset of severe malaria cases and malaria death cases aggregated at the health district level is available at the Directorate General of Studies and Sectoral Statistics of the Ministry of Health in Burkina Faso.

Author Contributions

Conceptualization: F.K.-S., S.S., M.O.; Methodology: M.O., F.K.-S., S.S.; Validation: F.K.-S., S.S., D.T.K., T.R., P.D.; Formal analysis: M.O., S.S.; Data curation: M.O., S.S.; Writing—original draft: M.O.; Writing—review & editing: D.T.K., T.R., S.S., P.D., F.K.-S.; Funding acquisition: F.K.-S., Project administration: F.K.-S.; Supervision: F.K.-S. All authors read and approved the final manuscript.

Funding

MO received a scholarship from the Conseil de l’Action Internationale (Ref. ADRI/CD/CAI-2016-NR 33); TR and DTK received a scholarship from the Académie de recherche et d’enseignement supérieur.

Acknowledgments

We thank the Directorate General of Studies and Sectoral Statistics of the Ministry of Health of Burkina Faso for access to the ENDOS-BF database.

Conflicts of Interest

The authors declare that they have no competing interests. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

BCIBayesian credible intervals
DICDeviance information criterion
DHIS2District Health Information System
INLAIntegrated Nested Laplace Approximation
mCFRmalaria case fatality risk
RDTRapid Diagnostic Test

Appendix A1. Estimating District-Level Interventions Coverage

The data on the availability of kits for the rapid diagnostic test (mRDT) for malaria, the availability of treatment for the management of malaria cases and the use of ACTs for the treatment of malaria used in this study come from periodic Survey on Availability and Readiness Assessment (2012, 2014, 2016) and Malaria Indicators surveys (MIS) for year 2014 and 2017. These data come from national survey designed to have a regionally representative estimate. A conditional autoregressive model (CAR) formulated in a Bayesian framework with a binomial distribution was performed to estimate indicators at the district level.
Let Z j be the number of health centers with mRDT kits in district j = 1, …, 70 and M j , the total number of health centers sampled and interviewed in district j. We assume that Z j follows a binomial distribution, namely Z j | M j , π ( j ) ~ Bin ( M j , π ( j ) ) ∀ j = 1, …, 70, where π ( j ) is the proportion of health center having a kit for the test in district j.
logit ( π ( j ) ) = β 0 + ω j
where β 0 is a constant and ω j , ∀ j = 1, …, 70 are modeled via a CAR prior. Each ω j conditioned by the neighbor ω k follows a normal distribution with an average equal to the average of the neighboring districts ω k and a variance inversely proportional to the number of neighboring districts m j , i.e., ω j | ω k ~ N ( γ l δ j   ω k ,   σ ω 2 m j ) , where γ quantifies the magnitude of the spatial correlation present in the data, σ ω 2 measures the spatial variance. ω j and ω k are adjacent districts in the set of all adjacent districts δ j of district j, and m j is the number of adjacent districts.
We assumed noninformative Gaussian distributions prior for β 0 ~ N (0, 102). A gamma distribution prior for σ ω 2 ~ Ga (0.1, 0.001).
We have use similar formulations for the availability of treatment in health centers and the use of artemisinin-based combination therapies (ACTs) for treatment in households.
Figure A1. Spatiotemporal trends in the crude mCFR in children under five (January 2013–December 2018. T1, T2, T3 and T4 represents the first, second, third and fourth trimester of each year respectively.
Figure A1. Spatiotemporal trends in the crude mCFR in children under five (January 2013–December 2018. T1, T2, T3 and T4 represents the first, second, third and fourth trimester of each year respectively.
Ijerph 17 01840 g0a1aIjerph 17 01840 g0a1bIjerph 17 01840 g0a1cIjerph 17 01840 g0a1d

Appendix A2. Choice of ZIP Model or ZNIB for Modelling

We performed statistic investigations to raise the over dispersion in the data. We concluded that there was an over dispersion in the data. Since the overdispersion was due to an excess of zero, a ZIP or ZNIB was suitable to handle this issue. To choose between ZIP or ZNIB we proceeded as follows.
It appeared from the two analyzes that the DIC of the ZIP was lower (DIC = 17,888.29) than the ZNIB model (DIC = 18,343.52). We then selected the ZIP model for our study.

References

  1. World Health Organization. World Malaria Report 2019. Available online: https://www.who.int/publications-detail/world-malaria-report-2019 (accessed on 21 November 2018).
  2. Rowe, A.K.; Rowe, S.Y.; Snow, R.W.; Korenromp, E.L.; Schellenberg, J.R.; Stein, C.; Nahlen, B.L.; Bryce, J.; Black, R.E.; Steketee, R.W. The burden of malaria mortality among African children in the year 2000. Int. J. Epidemiol. 2006, 35, 691–704. [Google Scholar] [CrossRef] [Green Version]
  3. Walker, K. A review of control methods for African malaria vectors. Environ. Health Proj. 2002, 2, 618–627. [Google Scholar]
  4. Ministère de la Santé. Annuaire Statistique 2018. Available online: http://www.cns.bf/IMG/pdf/annuaire_ms_2018.pdf (accessed on 21 September 2019).
  5. Ministère de la Santé. Plan Stratégique National de Lutte Contre le Paludisme du Burkina Faso 2016–2020. Available online: http://onsp-sante.bf/sites/default/files/publications/166/PSN%20%20%20%20%202016-2020_Paludisme_20_02_2017.pdf (accessed on 21 October 2019).
  6. Chuma, J.; Okungu, V.; Molyneux, C. Barriers to prompt and effective malaria treatment among the poorest population in Kenya. Malar. J. 2010, 9, 144. [Google Scholar] [CrossRef] [Green Version]
  7. Halwindi, H.; Siziya, S.; Magnussen, P.; Olsen, A. Factors perceived by caretakers as barriers to health care for under-five children in Mazabuka District, Zambia. Int. Sch. Res. Not. 2013, 2013, 905836. [Google Scholar] [CrossRef]
  8. Sieleunou, I.; Betsi, E.; Kouontchou, J.; Nguela, A.T.; Yumo, H.A. Removing the user fees for the under-fives’ simple malaria treatment in Cameroon effect on the health services utilization. J. Community Med. Health Educ. 2015, 5, 383. [Google Scholar] [CrossRef] [Green Version]
  9. Chuma, J.; Musimbi, J.; Okungu, V.; Goodman, C.; Molyneux, C. Reducing user fees for primary health care in Kenya: Policy on paper or policy in practice? Int. J. Equity Health 2009, 8, 15. [Google Scholar] [CrossRef] [Green Version]
  10. Lagarde, M.; Palmer, N. The impact of user fees on health service utilization in low- and middle-income countries: How strong is the evidence? Bull. World Health Organ. 2008, 86, 839–848. [Google Scholar] [CrossRef] [PubMed]
  11. McKinnon, B.; Harper, S.; Kaufman, J.S.; Bergevin, Y. Removing user fees for facility-based delivery services: A difference-in-differences evaluation from ten sub-Saharan African countries. Health Policy Plan. 2015, 30, 432–441. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Atchessi, N.; Ridde, V.; Zunzunegui, M.V. User fees exemptions alone are not enough to increase indigent use of healthcare services. Health Policy Plan. 2016, 31, 674–681. [Google Scholar] [CrossRef] [PubMed]
  13. Heinmüller, R.; Dembélé, Y.A.; Jouquet, G.; Haddad, S.; Ridde, V. Free healthcare provision with an NGO or by the Malian government. Impact on health center attendance by children under five. Field Actions Sci. Rep. 2013. Available online: http://journals.openedition.org/factsreports/1731 (accessed on 21 September 2019).
  14. Nyonator, F.; Kutzin, J. Health for some? The effects of user fees in the Volta Region of Ghana. Health Policy Plan. 1999, 14, 329–341. [Google Scholar] [CrossRef] [Green Version]
  15. Penfold, S.; Harrison, E.; Bell, J.; Fitzmaurice, A. Evaluation of the delivery fee exemption policy in ghana: Population estimates of changes in delivery service utilization in two regions. Ghana Med. J. 2007, 41, 100–109. [Google Scholar]
  16. Ponsar, F.; Van Herp, M.; Zachariah, R.; Gerard, S.; Philips, M.; Jouquet, G. Abolishing user fees for children and pregnant women trebled uptake of malaria-related interventions in Kangaba, Mali. Health Policy Plan. 2011, 26 (Suppl 2), ii72–ii83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Organisation International. Décret 2016-311_PRES/PM/MS/MATDSI/MINEFID Portant Gratuité de Soins au Profit des Femmes et des Enfants de Moins de Cinq ans Vivant au Burkina Faso. JO N°22.du 02 juin 2016. Available online: http://www.ilo.org/dyn/natlex/docs/ELECTRONIC/104122/126889/F331618513/BFA-104122.pdf (accessed on 21 September 2019).
  18. Ministère de la Santé. Annuaire Statistique 2016. Available online: http://www.cns.bf/IMG/pdf/annuaire_2016_ms_signe.pdf (accessed on 21 September 2019).
  19. Ministère de la Santé. Annuaire Statistique 2017. Available online: http://www.cns.bf/IMG/pdf/annuaire_ms_2017.pdf (accessed on 21 September 2019).
  20. Alegana, V.A.; Atkinson, P.M.; Lourenco, C.; Ruktanonchai, N.W.; Bosco, C.; Erbach-Schoenberg, E.Z.; Didier, B.; Pindolia, D.; Menach, A.L.; Katokele, S.; et al. Advances in mapping malaria for elimination: Fine resolution modelling of Plasmodium falciparum incidence. Sci. Rep. 2016, 6, 29628. [Google Scholar] [CrossRef] [PubMed]
  21. Ouedraogo, M.; Samadoulougou, S.; Rouamba, T.; Hien, H.; Sawadogo, J.E.M.; Tinto, H.; Alegana, V.A.; Speybroeck, N.; Kirakoya-Samadoulougou, F. Spatial distribution and determinants of asymptomatic malaria risk among children under 5 years in 24 districts in Burkina Faso. Malar. J. 2018, 17, 460. [Google Scholar] [CrossRef] [PubMed]
  22. Samadoulougou, S.; Maheu-Giroux, M.; Kirakoya-Samadoulougou, F.; De Keukeleire, M.; Castro, M.C.; Robert, A. Multilevel and geo-statistical modeling of malaria risk in children of Burkina Faso. Parasites Vectors 2014, 7, 350. [Google Scholar] [CrossRef] [Green Version]
  23. Ssempiira, J.; Kissa, J.; Nambuusi, B.; Kyozira, C.; Rutazaana, D.; Mukooyo, E.; Opigo, J.; Makumbi, F.; Kasasa, S.; Vounatsou, P. The effect of case management and vector-control interventions on space-time patterns of malaria incidence in Uganda. Malar. J. 2018, 17, 162. [Google Scholar] [CrossRef]
  24. Besag, J.; York, J.; Mollié, A. Bayesian image restoration, with two applications in spatial statistics. Ann. Inst. Stat. Math. 1991, 43, 1–20. [Google Scholar] [CrossRef]
  25. Stern, H. Inference for extremes in disease mapping. In Disease Mapping and Risk Assessment for Public Health; Lawson, A.B., Biggeri, A., Bohning, D., Lesaffre, E., Viel, J.F., Bertollini, R., Eds.; John Wiley & Son: Hoboken, NJ, USA, 1999; pp. 63–84. [Google Scholar]
  26. Ouedraogo, M.; Rouamba, T.; Samadoulougou, S.; Kirakoya-Samadoulougou, F. Effect of Free Healthcare Policy for Children under Five Years Old on the Incidence of Reported Malaria Cases in Burkina Faso by Bayesian Modelling: “Not only the Ears but also the Head of the Hippopotamus”. Int. J. Environ. Res. Public Health 2020, 17, 417. [Google Scholar] [CrossRef] [Green Version]
  27. Rouamba, T.; Samadoulougou, S.; Tinto, H.; Alegana, V.A.; Kirakoya-Samadoulougou, F. Bayesian Spatiotemporal Modeling of Routinely Collected Data to Assess the Effect of Health Programs in Malaria Incidence During Pregnancy in Burkina Faso. Sci. Rep. 2020, 10, 2618. [Google Scholar] [CrossRef] [Green Version]
  28. Rouamba, T.; Samadoulougou, S.; Tinto, H.; Alegana, V.A.; Kirakoya-Samadoulougou, F. Severe-malaria infection and its outcomes among pregnant women in Burkina Faso health-districts: Hierarchical Bayesian space-time models applied to routinely-collected data from 2013 to 2018. Spat. Spatio-Temporal Epidemiol. 2020, 33, 100333. [Google Scholar] [CrossRef]
  29. Institut National de la Statistique et de la Démographie. Projection Démographique de 2007–2020. 2009. Available online: http://www.cns.bf/IMG/pdf/projections_demographiques_sous_nationales.pdf (accessed on 21 September 2019).
  30. Ministère de la Santé. Carte_Sanitaire_2010. Available online: http://www.cns.bf/IMG/pdf/carte_sanitaire_2010.pdf (accessed on 21 September 2019).
  31. Ministère de la Santé. Plan National de Développement Sanitaire. 2011, pp. 3–22. Available online: https://www.uhc2030.org/fileadmin/uploads/ihp/Documents/Country_Pages/Burkina_Faso/Burkina_Faso_National_Health_Strategy_2011-2020_French.pdf (accessed on 21 September 2019).
  32. Institut national de la Statistique et de la Démographie. Enquête Multisectorielle Continue sur les Conditions de vie des Ménages. 2014. Available online: http://www.insd.bf/n/contenu/enquetes_recensements/Enq_EMC/Profil_de_pauvrete_et_d_inegalite_en_2014.pdf (accessed on 21 October 2019).
  33. Institut National de la Statistique et de la Démographie. Enqête Démographique et de Santé et à Indicateurs Multiples (EDSBF-MICS IV) 2010. 2012. Available online: https://dhsprogram.com/pubs/pdf/FR256/FR256.pdf (accessed on 21 October 2019).
  34. Institut National de la Statistique et de la Démographie. Rapport de L’enquête Nationale sur les Indicateurs du Paludisme 2017–2018. Available online: http://www.insd.bf/n/contenu/enquetes_recensements/enquete_palu/EIPBF_2018.pdf (accessed on 21 October 2019).
  35. Ministère de la Santé. Entrepôt de Données Sanitaires du Burkina Faso -Endos-BF, Burkina Faso 2013. Available online: https://burkina.dhis2.org/dhis/dhis-web-commons/security/login.action2013 (accessed on 21 October 2019).
  36. Ministère de la Santé. Rapport D’évaluation du Système de Gestion et de la Qualité des Données de Routine du SNIS; Direction Générale des Etudes et des Statistique Sectorielle: Ouagadougou, Burkina Faso, 2018. [Google Scholar]
  37. Institut national de la Statistique et de la Démographie. Rapport de L’enquête Nationale sur les Indicateurs du Paludisme. 2014. Available online: http://www.cns.bf/IMG/pdf/mis19.pdf (accessed on 21 November 2019).
  38. Ministère de la Santé/Direction Générale des Etudes et des Statistiques Sectorielles. Enquête Nationale sur les Prestations des Services de Santé et la Qualité des Données Sanitaires. 2012. Available online: http://apps.who.int/healthinfo/systems/datacatalog/index.php/catalog/50/reports (accessed on 21 November 2019).
  39. Rytter, M.J.; Kolte, L.; Briend, A.; Friis, H.; Christensen, V.B. The immune system in children with malnutrition--a systematic review. PLoS ONE 2014, 9, e105017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Chotsiri, P.; Denoeud-Ndam, L.; Baudin, E.; Guindo, O.; Diawara, H.; Attaher, O.; Smit, M.; Guerin, P.J.; Doumbo, O.K.; Wiesner, L.; et al. Severe Acute Malnutrition Results in Lower Lumefantrine Exposure in Children Treated With Artemether-Lumefantrine for Uncomplicated Malaria. Clin. Pharmacol. Ther. 2019, 106, 1299–1309. [Google Scholar] [CrossRef] [PubMed]
  41. Rice, A.L.; Sacco, L.; Hyder, A.; Black, R.E. Malnutrition as an underlying cause of childhood deaths associated with infectious diseases in developing countries. Bull. World Health Organ. 2000, 78, 1207–1221. [Google Scholar] [PubMed]
  42. UNHCR. Acute Malnutrition Threshold. UNHCR Emergency Handbook 2020. Available online: https://emergency.unhcr.org/entry/32604/acute-malnutrition-threshold (accessed on 22 January 2020).
  43. Ministère de la Santé. Annuaire Statistique 2013. Available online: http://www.cns.bf/IMG/pdf/annuaire_sante_2013.pdf (accessed on 21 November 2019).
  44. Ministère de la Santé. Annuaire Statistique 2014. Available online: http://www.cns.bf/IMG/pdf/annuaire_2014_du_ms.pdf (accessed on 21 November 2019).
  45. Ministère de la Santé. Annuaire Statistique 2015. Available online: http://www.cns.bf/IMG/pdf/annuaire_ms_2015_signe.pdf (accessed on 21 November 2019).
  46. Nkollo, B.; Claude, W. (Eds.) La nécessité d’une Gestion Optimale des Ressources Humaines. In Proceedings of the Conférence Régionale sur la Crise des Ressources Humaines en Santé en Afrique Subsaharienne, Douala, Cameroon, June 2007; Available online: https://www.who.int/workforcealliance/Cameroun.pdf (accessed on 21 November 2019).
  47. Ministère de la Santé. Métadonnées des Indicateurs du Système National D’information Sanitaire (SNIS). 2015. Available online: http://www.cns.bf/IMG/Metadonnees/Meta_donnees_SNIS.pdf (accessed on 21 November 2019).
  48. McCullagh, P.; Nelder, J. Generalized Linear Models; Monographs on Statistics and Applied Probability 37; Chapman and Hall: London, UK, 1989. [Google Scholar]
  49. Lawson, A.B. Bayesian Disease Mapping: Hierarchical Modeling in Spatial Epidemiology; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  50. Gelman, A. Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper). Math. Stat. Online 2006, 1, 515–534. [Google Scholar] [CrossRef]
  51. Richardson, S.; Thomson, A.; Best, N.; Elliott, P. Interpreting posterior relative risk estimates in disease-mapping studies. Environ. Health Perspect. 2004, 112, 1016–1025. [Google Scholar] [CrossRef]
  52. Quantum Gis, Development Team. QGIS Geographic Information System, Version 2.18. 2017. Available online: https://qgis.org/en/site/forusers/download.html (accessed on 21 October 2019).
  53. Dolan, C.B.; BenYishay, A.; Grépin, K.A.; Tanner, J.C.; Kimmel, A.D.; Wheeler, D.C.; McCord, G.C. The impact of an insecticide treated bednet campaign on all-cause child mortality: A geospatial impact evaluation from the Democratic Republic of Congo. PLoS ONE 2019, 14, e0212890. [Google Scholar] [CrossRef]
  54. Valle, D.; Millar, J.; Amratia, P. Spatial heterogeneity can undermine the effectiveness of country-wide test and treat policy for malaria: A case study from Burkina Faso. Malar. J. 2016, 15, 513. [Google Scholar] [CrossRef] [Green Version]
  55. Gething, P.W.; Casey, D.C.; Weiss, D.J.; Bisanzio, D.; Bhatt, S.; Cameron, E.; Battle, K.E.; Dalrymple, U.; Rozier, J.; Rao, P.C.; et al. Mapping Plasmodium falciparum Mortality in Africa between 1990 and 2015. N. Engl. J. Med. 2016, 375, 2435–2445. [Google Scholar] [CrossRef] [Green Version]
  56. Zurovac, D.; Githinji, S.; Memusi, D.; Kigen, S.; Machini, B.; Muturi, A.; Otieno, G.; Snow, R.W.; Nyandigisi, A. Major improvements in the quality of malaria case-management under the “test and treat” policy in Kenya. PLoS ONE 2014, 9, e92782. [Google Scholar] [CrossRef]
  57. O’Meara, W.P.; Mangeni, J.N.; Steketee, R.; Greenwood, B. Changes in the burden of malaria in sub-Saharan Africa. Lancet Infect. Dis. 2010, 10, 545–555. [Google Scholar] [CrossRef]
  58. Fatungase, K.O.; Amoran, O.E.; Alausa, K.O. The effect of health education intervention on the home management of malaria among the caregivers of children aged under 5 years in Ogun State, Nigeria. Eur. J. Med. Res. 2012, 17, 11. [Google Scholar] [CrossRef] [Green Version]
  59. Kenangalem, E.; Poespoprodjo, J.R.; Douglas, N.M.; Burdam, F.H.; Gdeumana, K.; Chalfein, F.; Prayoga; Thio, F.; Devine, A.; Martfurt, J.; et al. Malaria morbidity and mortality following introduction of a universal policy of artemisinin-based treatment for malaria in Papua, Indonesia: A longitudinal surveillance study. PLoS Med. 2019, 16, e1002815. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Nkwenti, H.E.; Ngowe, M.N.; Fokam, P.; Fonyuy, J.N.; Atanga, S.N.; Nkfusai, N.C.; Yankam, B.M.; Tsoka-Gwegweni, J.M.; Cumber, S.N. The effect of subsidized malaria treatment among under-five children in the Buea Health District, Cameroon. Pan Afr. Med. J. 2019, 33, 152. [Google Scholar] [CrossRef] [PubMed]
  61. Ridde, V.; Robert, E.; Meessen, B. A literature review of the disruptive effects of user fee exemption policies on health systems. BMC Public Health 2012, 12, 289. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Faye, O.; Correa, J.; Camara, B.; Dieng, T.; Dieng, Y.; Gaye, O.; Bah, I.B.; N’Dir, O.; Fall, M.; Diallo, S. Malaria lethality in Dakar pediatric environment: Study of risk factors. Med. Trop. 1998, 58, 361–364. [Google Scholar]
  63. Newton, C.R.; Taylor, T.E.; Whitten, R.O. Pathophysiology of fatal falciparum malaria in African children. Am. J. Trop. Med. Hyg. 1998, 58, 673–683. [Google Scholar] [CrossRef]
  64. Idro, R.; Marsh, K.; John, C.C.; Newton, C.R.J. Cerebral Malaria: Mechanisms of Brain Injury and Strategies for Improved Neurocognitive Outcome. Pediatric Res. 2010, 68, 267–274. [Google Scholar] [CrossRef]
  65. Sypniewska, P.; Duda, J.F.; Locatelli, I.; Althaus, C.R.; Althaus, F.; Genton, B. Clinical and laboratory predictors of death in African children with features of severe malaria: A systematic review and meta-analysis. BMC Med. 2017, 15, 147. [Google Scholar] [CrossRef] [Green Version]
  66. Conroy, A.L.; Hawkes, M.; Elphinstone, R.E.; Morgan, C.; Hermann, L.; Barker, K.R.; Namasopo, S.; Opoka, R.O.; John, C.C.; Liles, W.C.; et al. Acute Kidney Injury Is Common in Pediatric Severe Malaria and Is Associated with Increased Mortality. Open Forum Infect. Dis. 2016, 3, ofw046. [Google Scholar] [CrossRef] [Green Version]
  67. Berkley, J.; Mwarumba, S.; Bramham, K.; Lowe, B.; Marsh, K. Bacteraemia complicating severe malaria in children. Trans. R. Soc. Trop. Med. Hyg. 1999, 93, 283–286. [Google Scholar] [CrossRef] [Green Version]
  68. Church, J.; Maitland, K. Invasive bacterial co-infection in African children with Plasmodium falciparum malaria: A systematic review. BMC Med. 2014, 12, 31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Sagbo, G.G.; Alihonou, F.; Alao, M.J.; Tohodjèdé, Y.; Zohoun, L.; Bognon, G.; Agossou, J.; Noudamadjo, A.; Orou-Guidou, A. The Determinants of Fatal Outcomes during Severe Malaria in Children at the HKM University Teaching Hospital of Cotonou-Benin. J. Pediatrics 2017, 7, 245. [Google Scholar] [CrossRef] [Green Version]
  70. Camara, B.; Diagne, N.; Faye, P.; Fall, M.; Ndiaye, J.; Ba, M.; Sow, H.D. Critères de gravité et facteurs pronostiques du paludisme chez l’enfant à Dakar. Med. Mal. Infect. 2011, 41, 63–67. [Google Scholar] [CrossRef] [PubMed]
  71. Pasvol, G. The treatment of complicated and severe malaria. Br. Med. Bull. 2005, 75, 29–47. [Google Scholar] [CrossRef] [Green Version]
  72. Mutombo, A.M.; Mukuku, O.; Tshibanda, K.N.; Swana, E.K.; Mukomena, E.; Ngwej, D.T.; Luboya, O.N.; Kakoma, J.B.; Wembonyama, S.O.; Geertruyden, J.P.V.; et al. Severe malaria and death risk factors among children under 5 years at Jason Sendwe Hospital in Democratic Republic of Congo. Pan Afr. Med. J. 2018, 29, 1–8. [Google Scholar] [CrossRef]
  73. Win Han, O.; Gold, L.; Moore, K.; Agius, P.A.; Fowkes, F.J.I. The impact of community-delivered models of malaria control and elimination: A systematic review. Malar. J. 2019, 18, 269. [Google Scholar] [CrossRef]
  74. Health Organization World. Global Technical Strategy for Malaria 2016–2030; Global Malaria Programme: Geneva, Switzerland, 2015; Available online: https://www.who.int/malaria/publications/atoz/9789241564991/en/ (accessed on 21 September 2019).
  75. Ministère de la Santé. Organigramme 2016. Available online: http://ilo.org/dyn/natlex/docs/ELECTRONIC/104132/126909/F1133913732/BFA-104132.pdf (accessed on 21 October 2019).
  76. Van Den Broeck, J.; Eeckels, R.; Vuylsteke, J. Influence of nutritional status on child mortality in rural Zaire. Lancet 1993, 341, 1491–1495. [Google Scholar] [CrossRef]
  77. Monjour, L.; Palminteri, R.; Froment, A.; Renault, T.; Alfred, C.; Gentilini, M.; Gouba, E. Is cell-mediated immune response related to nutritional state, but unaffected by concomitant malarial infection? Ann. Trop. Med. Parasitol. 1982, 76, 575–577. [Google Scholar] [CrossRef]
  78. Genton, B.; al-Yaman, F.; Alpers, M.P.; Mokela, D. Indicators of fatal outcome in paediatric cerebral malaria: A study of 134 comatose Papua New Guinean children. Int. J. Epidemiol. 1997, 26, 670–676. [Google Scholar] [CrossRef]
  79. Renaudin, P. Evaluation de l’état nutritionnel de l’enfant de moins de 5 ans à Moundou, Tchad: Relations avec la morbidité et la mortalité hospitalières. Med. Trop. 1997, 57, 49–54. [Google Scholar] [CrossRef]
  80. Hughes, S.; Kelly, P. Interactions of malnutrition and immune impairment, with specific reference to immunity against parasites. Parasite Immunol. 2006, 28, 577–588. [Google Scholar] [CrossRef] [PubMed]
  81. Bourke, C.D.; Berkley, J.A.; Prendergast, A.J. Immune Dysfunction as a Cause and Consequence of Malnutrition. Trends Immunol. 2016, 37, 386–398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Shankar, A.H. Nutritional modulation of malaria morbidity and mortality. J. Infect. Dis. 2000, 182 (Suppl. 1), S37–S53. [Google Scholar] [CrossRef]
  83. Murray, M.J.; Murray, A.B.; Murray, N.J.; Murray, M.B. Diet and cerebral malaria: The effect of famine and refeeding. Am. J. Clin. Nutr. 1978, 31, 57–61. [Google Scholar] [CrossRef]
  84. Murray, M.J.; Murray, A.B.; Murray, M.B.; Murray, C.J. Somali food shelters in the Ogaden famine and their impact on health. Lancet 1976, 1, 1283–1285. [Google Scholar] [CrossRef]
  85. Murray, M.J.; Murray, N.J.; Murray, A.B.; Murray, M.B. Refeeding-malaria and hyperferraemia. Lancet 1975, 1, 653–654. [Google Scholar] [CrossRef]
  86. Nacher, M.; Singhasivanon, P.; Treeprasertsuk, S.; Vannaphan, S.; Traore, B.; Looareesuwan, S.; Gay, F. Intestinal helminths and malnutrition are independently associated with protection from cerebral malaria in Thailand. Ann. Trop. Med. Parasitol. 2002, 96, 5–13. [Google Scholar] [CrossRef]
  87. Ferreira, E.; Alexandre, M.A.; Salinas, J.L.; de Siqueira, A.M.; Benzecry, S.G.; de Lacerda, M.V.; Monteiro, W.M. Association between anthropometry-based nutritional status and malaria: A systematic review of observational studies. Malar. J. 2015, 14, 346. [Google Scholar] [CrossRef] [Green Version]
  88. Mockenhaupt, F.P.; Ehrhardt, S.; Burkhardt, J.; Bosomtwe, S.Y.; Laryea, S.; Anemana, S.D.; Otchwemah, R.N.; Cramer, J.P.; Dietz, E.; Gellert, S.; et al. Manifestation and outcome of severe malaria in children in northern Ghana. Am. J. Trop. Med. Hyg. 2004, 71, 167–172. [Google Scholar] [CrossRef] [Green Version]
  89. Orimadegun, A.E.; Fawole, O.; Okereke, J.O.; Akinbami, F.O.; Sodeinde, O. Increasing burden of childhood severe malaria in a Nigerian tertiary hospital: Implication for control. J. Trop. Pediatrics 2007, 53, 185–189. [Google Scholar] [CrossRef] [Green Version]
  90. Von Seidlein, L.; Olaosebikan, R.; Hendriksen, I.C.; Lee, S.J.; Adedoyin, O.T.; Agbenyega, T.; Nguah, S.B.; Bojang, K.; Deen, J.L.; Evans, J.; et al. Predicting the clinical outcome of severe falciparum malaria in african children: Findings from a large randomized trial. Clin. Infect. Dis. 2012, 54, 1080–1090. [Google Scholar] [CrossRef]
  91. Olumese, P.E.; Sodeinde, O.; Ademowo, O.G.; Walker, O. Protein energy malnutrition and cerebral malaria in Nigerian children. J. Trop. Pediatrics 1997, 43, 217–219. [Google Scholar] [CrossRef]
  92. Kendjo, E.; Agbenyega, T.; Bojang, K.; Newton, C.R.; Bouyou-Akotet, M.; Pedross, F.; Kombila, M.; Helbok, R.; Kremsner, P.G. Mortality patterns and site heterogeneity of severe malaria in African children. PLoS ONE 2013, 8, e58686. [Google Scholar] [CrossRef]
  93. Jenkins, R.; Omollo, R.; Ongecha, M.; Sifuna, P.; Othieno, C.; Ongeri, L.; Kingora, J.; Ogutu, B. Prevalence of malaria parasites in adults and its determinants in malaria endemic area of Kisumu County, Kenya. Malar. J. 2015, 14, 263. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Evolution of malaria case fatality rate in children under five from 2013 to 2018.
Figure 1. Evolution of malaria case fatality rate in children under five from 2013 to 2018.
Ijerph 17 01840 g001
Figure 2. Spatiotemporal trends in the predicted malaria case fatality rate per 100 severe malaria cases in children under five, estimated from the Bayesian spatiotemporal model (January 2013–December 2018). T1, T2, T3 and T4 represent the first, second, third and fourth trimester (quarter) of each year respectively.
Figure 2. Spatiotemporal trends in the predicted malaria case fatality rate per 100 severe malaria cases in children under five, estimated from the Bayesian spatiotemporal model (January 2013–December 2018). T1, T2, T3 and T4 represent the first, second, third and fourth trimester (quarter) of each year respectively.
Ijerph 17 01840 g002aIjerph 17 01840 g002bIjerph 17 01840 g002c
Figure 3. Burkina Faso health districts severe malaria deaths among children under five: exceedance probability ( P r ^ ( e u i + ζ i > 1.01 / y ) ) of posterior relative risk based on Richardson’s classification. The number in bracket indicates the number of health districts within each category.
Figure 3. Burkina Faso health districts severe malaria deaths among children under five: exceedance probability ( P r ^ ( e u i + ζ i > 1.01 / y ) ) of posterior relative risk based on Richardson’s classification. The number in bracket indicates the number of health districts within each category.
Ijerph 17 01840 g003
Table 1. Association of health programs with severe malaria case fatality rate estimated from the spatio-temporal model adjusted by contextual factors.
Table 1. Association of health programs with severe malaria case fatality rate estimated from the spatio-temporal model adjusted by contextual factors.
PredictorAdjusted IRR (95% BCI)
Free of charge health care
No1
Yes0.47 (0.38, 0.58)
Availability of mRDT a0.54 (0.47, 0.62)
Availability of treatment a0.50 (0.41, 0.61)
Use ACT for treatment a0.54 (0.47, 0.62)
Global acute malnutrition
<10%1
≥10%1.06 (0.99, 1.13)
Distance (km) b
<51
≥51.12(1.03, 1.22)
Population nurse ratio
<50001
≥50001.07 (0.99, 1.15)
Spatial varianceMean [95% BCI]
Non-structural Variance1.036 (1.035, 1.037)
Structural Variance1.028 (1.027, 1.029)
Temporal varianceMean [95% BCI]
Structural variance1.209 (1.063, 1.955)
Interaction variance1.170 (1.169, 1.172)
a Coverage (availability of malaria rapid diagnostic test (mRDT) and treatment) and use of artemisinin-based combination therapy (ACT) for treatment were modelled on a scale of 0 to 1; therefore, a one unit increase in coverage corresponds to a 100% increase, which implies a shift of the current value by 100%. b Distance to the nearest health facility.

Share and Cite

MDPI and ACS Style

Ouédraogo, M.; Kangoye, D.T.; Samadoulougou, S.; Rouamba, T.; Donnen, P.; Kirakoya-Samadoulougou, F. Malaria Case Fatality Rate among Children under Five in Burkina Faso: An Assessment of the Spatiotemporal Trends Following the Implementation of Control Programs. Int. J. Environ. Res. Public Health 2020, 17, 1840. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17061840

AMA Style

Ouédraogo M, Kangoye DT, Samadoulougou S, Rouamba T, Donnen P, Kirakoya-Samadoulougou F. Malaria Case Fatality Rate among Children under Five in Burkina Faso: An Assessment of the Spatiotemporal Trends Following the Implementation of Control Programs. International Journal of Environmental Research and Public Health. 2020; 17(6):1840. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17061840

Chicago/Turabian Style

Ouédraogo, Mady, David Tiga Kangoye, Sékou Samadoulougou, Toussaint Rouamba, Philippe Donnen, and Fati Kirakoya-Samadoulougou. 2020. "Malaria Case Fatality Rate among Children under Five in Burkina Faso: An Assessment of the Spatiotemporal Trends Following the Implementation of Control Programs" International Journal of Environmental Research and Public Health 17, no. 6: 1840. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17061840

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop