Trends in treatment-seeking for fever in children under five years old in 151 countries from 1990 to 2020

Access to medical treatment for fever is essential to prevent morbidity and mortality in individuals and to prevent transmission of communicable febrile illness in communities. Quantification of the rates at which treatment is accessed is critical for health system planning and a prerequisite for disease burden estimates. In this study, national data on the proportion of children under five years old with fever who were taken for medical treatment were collected from all available countries in Africa, Latin America, and Asia (n = 91). We used generalised additive mixed models to estimate 30-year trends in the treatment-seeking rates across the majority of countries in these regions (n = 151). Our results show that the proportions of febrile children brought for medical treatment increased steadily over the last 30 years, with the greatest increases occurring in areas where rates had originally been lowest, which includes Latin America and Caribbean, North Africa and the Middle East (51 and 50% increase, respectively), and Sub-Saharan Africa (23% increase). Overall, the aggregated and population-weighted estimate of children with fever taken for treatment at any type of facility rose from 61% (59-64 95% CI) in 1990 to 71% (69-72 95% CI) in 2020. The overall population-weighted average for fraction of treatment in the public sector was largely unchanged during the study period: 49% (42-58 95% CI) sought care at public facilities in 1990 and 47% (44-52 95% CI) in 2020. Overall, the findings indicate that improvements in access to care have been made where they were most needed, but that despite rapid initial gains, progress can plateau without substantial investment. In 2020 there remained significant gaps in care utilisation that must be factored in when developing control strategies and deriving disease burden estimates.

Quantifying malaria acquired during travel and its role in malaria elimination on Bioko Island

BACKGROUND

Malaria elimination is the goal for Bioko Island, Equatorial Guinea. Intensive interventions implemented since 2004 have reduced prevalence, but progress has stalled in recent years. A challenge for elimination has been malaria infections in residents acquired during travel to mainland Equatorial Guinea. The present article quantifies how off-island contributes to remaining malaria prevalence on Bioko Island, and investigates the potential role of a pre-erythrocytic vaccine in making further progress towards elimination.

METHODS

Malaria transmission on Bioko Island was simulated using a model calibrated based on data from the Malaria Indicator Surveys (MIS) from 2015 to 2018, including detailed travel histories and malaria positivity by rapid-diagnostic tests (RDTs), as well as geospatial estimates of malaria prevalence. Mosquito population density was adjusted to fit local transmission, conditional on importation rates under current levels of control and within-island mobility. The simulations were then used to evaluate the impact of two pre-erythrocytic vaccine distribution strategies: mass treat and vaccinate, and prophylactic vaccination for off-island travellers. Lastly, a sensitivity analysis was performed through an ensemble of simulations fit to the Bayesian joint posterior probability distribution of the geospatial prevalence estimates.

RESULTS

The simulations suggest that in Malabo, an urban city containing 80% of the population, there are some pockets of residual transmission, but a large proportion of infections are acquired off-island by travellers to the mainland. Outside of Malabo, prevalence was mainly attributable to local transmission. The uncertainty in the local transmission vs. importation is lowest within Malabo and highest outside. Using a pre-erythrocytic vaccine to protect travellers would have larger benefits than using the vaccine to protect residents of Bioko Island from local transmission. In simulations, mass treatment and vaccination had short-lived benefits, as malaria prevalence returned to current levels as the vaccine's efficacy waned. Prophylactic vaccination of travellers resulted in longer-lasting reductions in prevalence. These projections were robust to underlying uncertainty in prevalence estimates.

CONCLUSIONS

The modelled outcomes suggest that the volume of malaria cases imported from the mainland is a partial driver of continued endemic malaria on Bioko Island, and that continued elimination efforts on must account for human travel activity.

Indirect effects of the COVID-19 pandemic on malaria intervention coverage, morbidity, and mortality in Africa: a geospatial modelling analysis

BACKGROUND

Substantial progress has been made in reducing the burden of malaria in Africa since 2000, but those gains could be jeopardised if the COVID-19 pandemic affects the availability of key malaria control interventions. The aim of this study was to evaluate plausible effects on malaria incidence and mortality under different levels of disruption to malaria control.

METHODS

Using an established set of spatiotemporal Bayesian geostatistical models, we generated geospatial estimates across malaria-endemic African countries of the clinical case incidence and mortality of malaria, incorporating an updated database of parasite rate surveys, insecticide-treated net (ITN) coverage, and effective treatment rates. We established a baseline estimate for the anticipated malaria burden in Africa in the absence of COVID-19-related disruptions, and repeated the analysis for nine hypothetical scenarios in which effective treatment with an antimalarial drug and distribution of ITNs (both through routine channels and mass campaigns) were reduced to varying extents.

FINDINGS

We estimated 215·2 (95% uncertainty interval 143·7-311·6) million cases and 386·4 (307·8-497·8) thousand deaths across malaria-endemic African countries in 2020 in our baseline scenario of undisrupted intervention coverage. With greater reductions in access to effective antimalarial drug treatment, our model predicted increasing numbers of cases and deaths: 224·1 (148·7-326·8) million cases and 487·9 (385·3-634·6) thousand deaths with a 25% reduction in antimalarial drug coverage; 233·1 (153·7-342·5) million cases and 597·4 (468·0-784·4) thousand deaths with a 50% reduction; and 242·3 (158·7-358·8) million cases and 715·2 (556·4-947·9) thousand deaths with a 75% reduction. Halting planned 2020 ITN mass distribution campaigns and reducing routine ITN distributions by 25%-75% also increased malaria burden to a total of 230·5 (151·6-343·3) million cases and 411·7 (322·8-545·5) thousand deaths with a 25% reduction; 232·8 (152·3-345·9) million cases and 415·5 (324·3-549·4) thousand deaths with a 50% reduction; and 234·0 (152·9-348·4) million cases and 417·6 (325·5-553·1) thousand deaths with a 75% reduction. When ITN coverage and antimalarial drug coverage were synchronously reduced, malaria burden increased to 240·5 (156·5-358·2) million cases and 520·9 (404·1-691·9) thousand deaths with a 25% reduction; 251·0 (162·2-377·0) million cases and 640·2 (492·0-856·7) thousand deaths with a 50% reduction; and 261·6 (167·7-396·8) million cases and 768·6 (586·1-1038·7) thousand deaths with a 75% reduction.

INTERPRETATION

Under pessimistic scenarios, COVID-19-related disruption to malaria control in Africa could almost double malaria mortality in 2020, and potentially lead to even greater increases in subsequent years. To avoid a reversal of two decades of progress against malaria, averting this public health disaster must remain an integrated priority alongside the response to COVID-19.

FUNDING

Bill and Melinda Gates Foundation; Channel 7 Telethon Trust, Western Australia.

Global estimation of anti-malarial drug effectiveness for the treatment of uncomplicated Plasmodium falciparum malaria 1991-2019

BACKGROUND

Anti-malarial drugs play a critical role in reducing malaria morbidity and mortality, but their role is mediated by their effectiveness. Effectiveness is defined as the probability that an anti-malarial drug will successfully treat an individual infected with malaria parasites under routine health care delivery system. Anti-malarial drug effectiveness (AmE) is influenced by drug resistance, drug quality, health system quality, and patient adherence to drug use; its influence on malaria burden varies through space and time.

METHODS

This study uses data from 232 efficacy trials comprised of 86,776 infected individuals to estimate the artemisinin-based and non-artemisinin-based AmE for treating falciparum malaria between 1991 and 2019. Bayesian spatiotemporal models were fitted and used to predict effectiveness at the pixel-level (5 km × 5 km). The median and interquartile ranges (IQR) of AmE are presented for all malaria-endemic countries.

RESULTS

The global effectiveness of artemisinin-based drugs was 67.4% (IQR: 33.3-75.8), 70.1% (43.6-76.0) and 71.8% (46.9-76.4) for the 1991-2000, 2006-2010, and 2016-2019 periods, respectively. Countries in central Africa, a few in South America, and in the Asian region faced the challenge of lower effectiveness of artemisinin-based anti-malarials. However, improvements were seen after 2016, leaving only a few hotspots in Southeast Asia where resistance to artemisinin and partner drugs is currently problematic and in the central Africa where socio-demographic challenges limit effectiveness. The use of artemisinin-based combination therapy (ACT) with a competent partner drug and having multiple ACT as first-line treatment choice sustained high levels of effectiveness. High levels of access to healthcare, human resource capacity, education, and proximity to cities were associated with increased effectiveness. Effectiveness of non-artemisinin-based drugs was much lower than that of artemisinin-based with no improvement over time: 52.3% (17.9-74.9) for 1991-2000 and 55.5% (27.1-73.4) for 2011-2015. Overall, AmE for artemisinin-based and non-artemisinin-based drugs were, respectively, 29.6 and 36% below clinical efficacy as measured in anti-malarial drug trials.

CONCLUSIONS

This study provides evidence that health system performance, drug quality and patient adherence influence the effectiveness of anti-malarials used in treating uncomplicated falciparum malaria. These results provide guidance to countries' treatment practises and are critical inputs for malaria prevalence and incidence models used to estimate national level malaria burden.

Modelling the impact of climate change on the distribution and abundance of tsetse in Northern Zimbabwe

Background

Climate change is predicted to impact the transmission dynamics of vector-borne diseases. Tsetse flies (Glossina) transmit species of Trypanosoma that cause human and animal African trypanosomiasis. A previous modelling study showed that temperature increases between 1990 and 2017 can explain the observed decline in abundance of tsetse at a single site in the Mana Pools National Park of Zimbabwe. Here, we apply a mechanistic model of tsetse population dynamics to predict how increases in temperature may have changed the distribution and relative abundance of Glossina pallidipes across northern Zimbabwe.

Methods

Local weather station temperature measurements were previously used to fit the mechanistic model to longitudinal G. pallidipes catch data. To extend the use of the model, we converted MODIS land surface temperature to air temperature, compared the converted temperatures with available weather station data to confirm they aligned, and then re-fitted the mechanistic model using G. pallidipes catch data and air temperature estimates. We projected this fitted model across northern Zimbabwe, using simulations at a 1 km × 1 km spatial resolution, between 2000 to 2016.

Results

We produced estimates of relative changes in G. pallidipes mortality, larviposition, emergence rates and abundance, for northern Zimbabwe. Our model predicts decreasing tsetse populations within low elevation areas in response to increasing temperature trends during 2000–2016. Conversely, we show that high elevation areas (> 1000 m above sea level), previously considered too cold to sustain tsetse, may now be climatically suitable.

Conclusions

To our knowledge, the results of this research represent the first regional-scale assessment of temperature related tsetse population dynamics, and the first high spatial-resolution estimates of this metric for northern Zimbabwe. Our results suggest that tsetse abundance may have declined across much of the Zambezi Valley in response to changing climatic conditions during the study period. Future research including empirical studies is planned to improve model accuracy and validate predictions for other field sites in Zimbabwe.

Environmental temperature and growth faltering in African children: a cross-sectional study

BACKGROUND

Child growth faltering persists in sub-Saharan Africa despite the scale-up of nutrition, water, and sanitation interventions over the past 2 decades. High temperatures have been hypothesised to contribute to child growth faltering via an adaptive response to heat, reduced appetite, and the energetic cost of thermoregulation. We did a cross-sectional study to assess whether child growth faltering is related to environmental temperature in sub-Saharan Africa.

METHODS

Data were extracted from 52 Demographic and Heath Surveys, dating from 2003 to 2016, that recorded anthropometric data in children aged 0-5 years, and were linked with remotely sensed monthly mean daytime land surface temperature for 2000-16. The odds of stunting (low height-for-age), wasting (low weight-for-height), and underweight (low weight-for-age) relative to monthly mean daytime land surface temperature were determined using multivariable logistic regression.

FINDINGS

The study population comprised 656 107 children resident in 373 012 households. Monthly mean daytime land surface temperature above 35°C was associated with increases in the odds of wasting (odds ratio 1·27, 95% CI 1·16-1·38; p<0·0001), underweight (1·09, 1·02-1·16; p=0·0073), and concurrent stunting with wasting (1·23, 1·07-1·41; p=0·0037), but a reduction in stunting (0·90, 0·85-0·96; p=0·00047) compared with a monthly mean daytime land surface temperature of less than 30°C.

INTERPRETATION

Children living in hotter parts of sub-Saharan Africa are more likely to be wasted, underweight, and concurrently stunted and wasted, but less likely to be stunted, than in cooler areas. Studies are needed to further investigate the relationship between temperature and child growth, including whether there is a direct effect not mediated by food security, regional wealth, and other environmental variables. Rising temperature, linked to anthropogenic climate change, might increase child growth faltering in sub-Saharan Africa.

FUNDING

UK Medical Research Council and UK Global Challenges Research Fund.

Housing and child health in sub-Saharan Africa: A cross-sectional analysis

BACKGROUND

Housing is essential to human well-being but neglected in global health. Today, housing in Africa is rapidly improving alongside economic development, creating an urgent need to understand how these changes can benefit health. We hypothesised that improved housing is associated with better health in children living in sub-Saharan Africa (SSA). We conducted a cross-sectional analysis of housing conditions relative to a range of child health outcomes in SSA.

METHODS AND FINDINGS

Cross-sectional data were analysed for 824,694 children surveyed in 54 Demographic and Health Surveys, 21 Malaria Indicator Surveys, and two AIDS Indicator Surveys conducted in 33 countries between 2001 and 2017 that measured malaria infection by microscopy or rapid diagnostic test (RDT), diarrhoea, acute respiratory infections (ARIs), stunting, wasting, underweight, or anaemia in children aged 0-5 years. The mean age of children was 2.5 years, and 49.7% were female. Housing was categorised into a binary variable based on a United Nations definition comparing improved housing (with improved drinking water, improved sanitation, sufficient living area, and finished building materials) versus unimproved housing (all other houses). Associations between house type and child health outcomes were determined using conditional logistic regression within surveys, adjusting for prespecified covariables including age, sex, household wealth, insecticide-treated bed net use, and vaccination status. Individual survey odds ratios (ORs) were pooled using random-effects meta-analysis. Across surveys, improved housing was associated with 8%-18% lower odds of all outcomes except ARI (malaria infection by microscopy: adjusted OR [aOR] 0.88, 95% confidence intervals [CIs] 0.80-0.97, p = 0.01; malaria infection by RDT: aOR 0.82, 95% CI 0.77-0.88, p < 0.001; diarrhoea: aOR 0.92, 95% CI 0.88-0.97, p = 0.001; ARI: aOR 0.96, 95% CI 0.87-1.07, p = 0.49; stunting: aOR 0.83, 95% CI 0.77-0.88, p < 0.001; wasting: aOR 0.90, 95% CI 0.83-0.99, p = 0.03; underweight: aOR 0.85, 95% CI 0.80-0.90, p < 0.001; any anaemia: aOR 0.87, 95% CI 0.82-0.92, p < 0.001; severe anaemia: aOR 0.89, 95% CI 0.84-0.95, p < 0.001). In comparison, insecticide-treated net use was associated with 16%-17% lower odds of malaria infection (microscopy: aOR 0.83, 95% CI 0.78-0.88, p < 0.001; RDT: aOR 0.84, 95% CI 0.79-0.88, p < 0.001). Drinking water source and sanitation facility alone were not associated with diarrhoea. The main study limitations are the use of self-reported diarrhoea and ARI, as well as potential residual confounding by socioeconomic position, despite adjustments for household wealth and education.

CONCLUSIONS

In this study, we observed that poor housing, which includes inadequate drinking water and sanitation facility, is associated with health outcomes known to increase child mortality in SSA. Improvements to housing may be protective against a number of important childhood infectious diseases as well as poor growth outcomes, with major potential to improve children's health and survival across SSA.

Mapping malaria seasonality in Madagascar using health facility data

BACKGROUND

Many malaria-endemic areas experience seasonal fluctuations in case incidence as Anopheles mosquito and Plasmodium parasite life cycles respond to changing environmental conditions. Identifying location-specific seasonality characteristics is useful for planning interventions. While most existing maps of malaria seasonality use fixed thresholds of rainfall, temperature, and/or vegetation indices to identify suitable transmission months, we construct a statistical modelling framework for characterising the seasonal patterns derived directly from monthly health facility data.

METHODS

With data from 2669 of the 3247 health facilities in Madagascar, a spatiotemporal regression model was used to estimate seasonal patterns across the island. In the absence of catchment population estimates or the ability to aggregate to the district level, this focused on the monthly proportions of total annual cases by health facility level. The model was informed by dynamic environmental covariates known to directly influence seasonal malaria trends. To identify operationally relevant characteristics such as the transmission start months and associated uncertainty measures, an algorithm was developed and applied to model realisations. A seasonality index was used to incorporate burden information from household prevalence surveys and summarise 'how seasonal' locations are relative to their surroundings.

RESULTS

Positive associations were detected between monthly case proportions and temporally lagged covariates of rainfall and temperature suitability. Consistent with the existing literature, model estimates indicate that while most parts of Madagascar experience peaks in malaria transmission near March-April, the eastern coast experiences an earlier peak around February. Transmission was estimated to start in southeast districts before southwest districts, suggesting that indoor residual spraying should be completed in the same order. In regions where the data suggested conflicting seasonal signals or two transmission seasons, estimates of seasonal features had larger deviations and therefore less certainty.

CONCLUSIONS

Monthly health facility data can be used to establish seasonal patterns in malaria burden and augment the information provided by household prevalence surveys. The proposed modelling framework allows for evidence-based and cohesive inferences on location-specific seasonal characteristics. As health surveillance systems continue to improve, it is hoped that more of such data will be available to improve our understanding and planning of intervention strategies.

Mapping Geospatial Processes Affecting the Environmental Fate of Agricultural Pesticides in Africa

The application of agricultural pesticides in Africa can have negative effects on human health and the environment. The aim of this study was to identify African environments that are vulnerable to the accumulation of pesticides by mapping geospatial processes affecting pesticide fate. The study modelled processes associated with the environmental fate of agricultural pesticides using publicly available geospatial datasets. Key geospatial processes affecting the environmental fate of agricultural pesticides were selected after a review of pesticide fate models and maps for leaching, surface runoff, sedimentation, soil storage and filtering capacity, and volatilization were created. The potential and limitations of these maps are discussed. We then compiled a database of studies that measured pesticide residues in Africa. The database contains 10,076 observations, but only a limited number of observations remained when a standard dataset for one compound was extracted for validation. Despite the need for more in-situ data on pesticide residues and application, this study provides a first spatial overview of key processes affecting pesticide fate that can be used to identify areas potentially vulnerable to pesticide accumulation.

Mapping the global endemicity and clinical burden of Plasmodium vivax, 2000-17: a spatial and temporal modelling study

BACKGROUND

Plasmodium vivax exacts a significant toll on health worldwide, yet few efforts to date have quantified the extent and temporal trends of its global distribution. Given the challenges associated with the proper diagnosis and treatment of P vivax, national malaria programmes-particularly those pursuing malaria elimination strategies-require up to date assessments of P vivax endemicity and disease impact. This study presents the first global maps of P vivax clinical burden from 2000 to 2017.

METHODS

In this spatial and temporal modelling study, we adjusted routine malariometric surveillance data for known biases and used socioeconomic indicators to generate time series of the clinical burden of P vivax. These data informed Bayesian geospatial models, which produced fine-scale predictions of P vivax clinical incidence and infection prevalence over time. Within sub-Saharan Africa, where routine surveillance for P vivax is not standard practice, we combined predicted surfaces of Plasmodium falciparum with country-specific ratios of P vivax to P falciparum. These results were combined with surveillance-based outputs outside of Africa to generate global maps.

FINDINGS

We present the first high-resolution maps of P vivax burden. These results are combined with those for P falciparum (published separately) to form the malaria estimates for the Global Burden of Disease 2017 study. The burden of P vivax malaria decreased by 41·6%, from 24·5 million cases (95% uncertainty interval 22·5-27·0) in 2000 to 14·3 million cases (13·7-15·0) in 2017. The Americas had a reduction of 56·8% (47·6-67·0) in total cases since 2000, while South-East Asia recorded declines of 50·5% (50·3-50·6) and the Western Pacific regions recorded declines of 51·3% (48·0-55·4). Europe achieved zero P vivax cases during the study period. Nonetheless, rates of decline have stalled in the past five years for many countries, with particular increases noted in regions affected by political and economic instability.

INTERPRETATION

Our study highlights important spatial and temporal patterns in the clinical burden and prevalence of P vivax. Amid substantial progress worldwide, plateauing gains and areas of increased burden signal the potential for challenges that are greater than expected on the road to malaria elimination. These results support global monitoring systems and can inform the optimisation of diagnosis and treatment where P vivax has most impact.

FUNDING

Bill & Melinda Gates Foundation and the Wellcome Trust.

Mapping the global prevalence, incidence, and mortality of Plasmodium falciparum, 2000-17: a spatial and temporal modelling study

BACKGROUND

Since 2000, the scale-up of malaria control interventions has substantially reduced morbidity and mortality caused by the disease globally, fuelling bold aims for disease elimination. In tandem with increased availability of geospatially resolved data, malaria control programmes increasingly use high-resolution maps to characterise spatially heterogeneous patterns of disease risk and thus efficiently target areas of high burden.

METHODS

We updated and refined the Plasmodium falciparum parasite rate and clinical incidence models for sub-Saharan Africa, which rely on cross-sectional survey data for parasite rate and intervention coverage. For malaria endemic countries outside of sub-Saharan Africa, we produced estimates of parasite rate and incidence by applying an ecological downscaling approach to malaria incidence data acquired via routine surveillance. Mortality estimates were derived by linking incidence to systematically derived vital registration and verbal autopsy data. Informed by high-resolution covariate surfaces, we estimated P falciparum parasite rate, clinical incidence, and mortality at national, subnational, and 5 × 5 km pixel scales with corresponding uncertainty metrics.

FINDINGS

We present the first global, high-resolution map of P falciparum malaria mortality and the first global prevalence and incidence maps since 2010. These results are combined with those for Plasmodium vivax (published separately) to form the malaria estimates for the Global Burden of Disease 2017 study. The P falciparum estimates span the period 2000-17, and illustrate the rapid decline in burden between 2005 and 2017, with incidence declining by 27·9% and mortality declining by 42·5%. Despite a growing population in endemic regions, P falciparum cases declined between 2005 and 2017, from 232·3 million (95% uncertainty interval 198·8-277·7) to 193·9 million (156·6-240·2) and deaths declined from 925 800 (596 900-1 341 100) to 618 700 (368 600-952 200). Despite the declines in burden, 90·1% of people within sub-Saharan Africa continue to reside in endemic areas, and this region accounted for 79·4% of cases and 87·6% of deaths in 2017.

INTERPRETATION

High-resolution maps of P falciparum provide a contemporary resource for informing global policy and malaria control planning, programme implementation, and monitoring initiatives. Amid progress in reducing global malaria burden, areas where incidence trends have plateaued or increased in the past 5 years underscore the fragility of hard-won gains against malaria. Efforts towards elimination should be strengthened in such areas, and those where burden remained high throughout the study period.

FUNDING

Bill & Melinda Gates Foundation.

The contribution of non-malarial febrile illness co-infections to Plasmodium falciparum case counts in health facilities in sub-Saharan Africa

Background

The disease burden of Plasmodium falciparum malaria illness is generally estimated using one of two distinct approaches: either by transforming P. falciparum infection prevalence estimates into incidence estimates using conversion formulae; or through adjustment of counts of recorded P. falciparum-positive fever cases from clinics. Whilst both ostensibly seek to evaluate P. falciparum disease burden, there is an implicit and problematic difference in the metric being estimated. The first enumerates only symptomatic malaria cases, while the second enumerates all febrile episodes coincident with a P. falciparum infection, regardless of the fever’s underlying cause.

Methods

Here, a novel approach was used to triangulate community-based data sources capturing P. falciparum infection, fever, and care-seeking to estimate the fraction of P. falciparum-positive fevers amongst children under 5 years of age presenting at health facilities that are attributable to P. falciparum infection versus other non-malarial causes. A Bayesian hierarchical model was used to assign probabilities of malaria-attributable fever (MAF) and non-malarial febrile illness (NMFI) to children under five from a dataset of 41 surveys from 21 countries in sub-Saharan Africa conducted between 2006 and 2016. Using subsequent treatment-seeking outcomes, the proportion of MAF and NMFI amongst P. falciparum-positive febrile children presenting at public clinics was estimated.

Results

Across all surveyed malaria-positive febrile children who sought care at public clinics across 41 country-years in sub-Saharan Africa, P. falciparum infection was estimated to be the underlying cause of only 37.7% (31.1–45.4, 95% CrI) of P. falciparum-positive fevers, with significant geographical and temporal heterogeneity between surveys.

Conclusions

These findings highlight the complex nature of the P. falciparum burden amongst children under 5 years of age and indicate that for many children presenting at health clinics, a positive P. falciparum diagnosis and a fever does not necessarily mean P. falciparum is the underlying cause of the child’s symptoms, and thus other causes of illness should always be investigated, in addition to prescribing an effective anti-malarial medication. In addition to providing new large-scale estimates of malaria-attributable fever prevalence, the results presented here improve comparability between different methods for calculating P. falciparum disease burden, with significant implications for national and global estimation of malaria burden.

Electronic supplementary material

The online version of this article (10.1186/s12936-019-2830-y) contains supplementary material, which is available to authorized users.

Human mobility patterns and malaria importation on Bioko Island

Malaria burden on Bioko Island has decreased significantly over the past 15 years. The impact of interventions on malaria prevalence, however, has recently stalled. Here, we use data from island-wide, annual malaria indicator surveys to investigate human movement patterns and their relationship to Plasmodium falciparum prevalence. Using geostatistical and mathematical modelling, we find that off-island travel is more prevalent in and around the capital, Malabo. The odds of malaria infection among off-island travelers are significantly higher than the rest of the population. We estimate that malaria importation rates are high enough to explain malaria prevalence in much of Malabo and its surroundings, and that local transmission is highest along the West Coast of the island. Despite uncertainty, these estimates of residual transmission and importation serve as a basis for evaluating progress towards elimination and for efficiently allocating resources as Bioko makes the transition from control to elimination.

Mapping changes in housing in sub-Saharan Africa from 2000 to 2015

Access to adequate housing is a fundamental human right, essential to human security, nutrition and health, and a core objective of the United Nations Sustainable Development Goals1,2. Globally, the housing need is most acute in Africa, where the population will more than double by 2050. However, existing data on housing quality across Africa are limited primarily to urban areas and are mostly recorded at the national level. Here we quantify changes in housing in sub-Saharan Africa from 2000 to 2015 by combining national survey data within a geostatistical framework. We show a marked transformation of housing in urban and rural sub-Saharan Africa between 2000 and 2015, with the prevalence of improved housing (with improved water and sanitation, sufficient living area and durable construction) doubling from 11% (95% confidence interval, 10-12%) to 23% (21-25%). However, 53 (50-57) million urban Africans (47% (44-50%) of the urban population analysed) were living in unimproved housing in 2015. We provide high-resolution, standardized estimates of housing conditions across sub-Saharan Africa. Our maps provide a baseline for measuring change and a mechanism to guide interventions during the era of the Sustainable Development Goals.

Heterogeneous exposure and hotspots for malaria vectors at three study sites in Uganda

Background: Heterogeneity in malaria transmission has household, temporal, and spatial components. These factors are relevant for improving the efficiency of malaria control by targeting heterogeneity. To quantify variation, we analyzed mosquito counts from entomological surveillance conducted at three study sites in Uganda that varied in malaria transmission intensity. Mosquito biting or exposure is a risk factor for malaria transmission.

Methods: Using a Bayesian zero-inflated negative binomial model, validated via a comprehensive simulation study, we quantified household differences in malaria vector density and examined its spatial distribution. We introduced a novel approach for identifying changes in vector abundance hotspots over time by computing the Getis-Ord statistic on ratios of household biting propensities for different scenarios. We also explored the association of household biting propensities with housing and environmental covariates.

Results: In each site, there was evidence for hot and cold spots of vector abundance, and spatial patterns associated with urbanicity, elevation, or other environmental covariates. We found some differences in the hotspots in rainy vs. dry seasons or before vs. after the application of control interventions. Housing quality explained a portion of the variation among households in mosquito counts.

Conclusion: This work provided an improved understanding of heterogeneity in malaria vector density at the three study sites in Uganda and offered a valuable opportunity for assessing whether interventions could be spatially targeted to be aimed at abundance hotspots which may increase malaria risk. Indoor residual spraying was shown to be a successful measure of vector control interventions in Tororo, Uganda.  Cement walls, brick floors, closed eaves, screened airbricks, and tiled roofs were features of a house that had shown reduction of household biting propensity. Improvements in house quality should be recommended as a supplementary measure for malaria control reducing risk of infection.

malariaAtlas: an R interface to global malariometric data hosted by the Malaria Atlas Project

Background

The Malaria Atlas Project (MAP) has worked to assemble and maintain a global open-access database of spatial malariometric data for over a decade. This data spans various formats and topics, including: geo-located surveys of malaria parasite rate; global administrative boundary shapefiles; and global and regional rasters representing the distribution of malaria and associated illnesses, blood disorders, and intervention coverage. MAP has recently released malariaAtlas, an R package providing a direct interface to MAP’s routinely-updated malariometric databases and research outputs.

Methods and results

The current paper reviews the functionality available in malariaAtlas and highlights its utility for spatial epidemiological analysis of malaria. malariaAtlas enables users to freely download, visualise and analyse global malariometric data within R. Currently available data types include: malaria parasite rate and vector occurrence point data; subnational administrative boundary shapefiles; and a large suite of rasters covering a diverse range of metrics related to malaria research. malariaAtlas is here used in two mock analyses to illustrate how this data may be incorporated into a standard R workflow for spatial analysis.

Conclusions

malariaAtlas is the first open-access R-interface to malariometric data, providing a new and reproducible means of accessing such data within a freely available and commonly used statistical software environment. In this way, the malariaAtlas package aims to contribute to the environment of data-sharing within the malaria research community.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2500-5) contains supplementary material, which is available to authorized users.

Mapping child growth failure in Africa between 2000 and 2015

Insufficient growth during childhood is associated with poor health outcomes and an increased risk of death. Between 2000 and 2015, nearly all African countries demonstrated improvements for children under 5 years old for stunting, wasting, and underweight, the core components of child growth failure. Here we show that striking subnational heterogeneity in levels and trends of child growth remains. If current rates of progress are sustained, many areas of Africa will meet the World Health Organization Global Targets 2025 to improve maternal, infant and young child nutrition, but high levels of growth failure will persist across the Sahel. At these rates, much, if not all of the continent will fail to meet the Sustainable Development Goal target—to end malnutrition by 2030. Geospatial estimates of child growth failure provide a baseline for measuring progress as well as a precision public health platform to target interventions to those populations with the greatest need, in order to reduce health disparities and accelerate progress.

Geospatial estimates of child growth failure in Africa provide a baseline for measuring progress and a precision public health platform to target interventions to those populations with the greatest need.

Heterogeneous exposure and hotspots for malaria vectors at three study sites in Uganda

Background: Heterogeneity in malaria transmission has household, temporal, and spatial components. These factors are relevant for improving the efficiency of malaria control by targeting heterogeneity. To quantify variation, we analyzed mosquito counts from entomological surveillance conducted at three study sites in Uganda that varied in malaria transmission intensity. Methods: Using a Bayesian zero-inflated negative binomial model, validated via a comprehensive simulation study, we quantified household differences in malaria vector density and examined its spatial distribution. We introduced a novel approach for identifying changes in malaria hotspots over time by computing the Getis-Ord statistic on ratios of household biting propensities for different scenarios. We also explored the association of household biting propensities with housing and environmental covariates. Results: In each site, there was evidence for hot and cold spots, spatial patterns associated with urbanicity, elevation, or other environmental covariates. We found some differences in the hotspots in rainy vs. dry seasons or before vs. after the application of control interventions. Housing quality explained a portion of the variation among households in mosquito counts. Conclusion: This work provided an improved understanding of heterogeneity in malaria vector density at the three study sites in Uganda and offered a valuable opportunity for assessing whether interventions could be spatially targeted to be aimed at hotspots of malaria risk. Indoor residual spraying was shown to be a successful measure of vector control interventions in Tororo, Uganda.  Cement walls, brick floors, closed eaves, screened airbricks, and tiled roofs were features of a house that had shown protective effects towards malaria risk. Improvements in house quality should be recommended as a supplementary measure for malaria control.

Spatio-temporal mapping of Madagascar's Malaria Indicator Survey results to assess Plasmodium falciparum endemicity trends between 2011 and 2016

BACKGROUND

Reliable measures of disease burden over time are necessary to evaluate the impact of interventions and assess sub-national trends in the distribution of infection. Three Malaria Indicator Surveys (MISs) have been conducted in Madagascar since 2011. They provide a valuable resource to assess changes in burden that is complementary to the country's routine case reporting system.

METHODS

A Bayesian geostatistical spatio-temporal model was developed in an integrated nested Laplace approximation framework to map the prevalence of Plasmodium falciparum malaria infection among children from 6 to 59 months in age across Madagascar for 2011, 2013 and 2016 based on the MIS datasets. The model was informed by a suite of environmental and socio-demographic covariates known to influence infection prevalence. Spatio-temporal trends were quantified across the country.

RESULTS

Despite a relatively small decrease between 2013 and 2016, the prevalence of malaria infection has increased substantially in all areas of Madagascar since 2011. In 2011, almost half (42.3%) of the country's population lived in areas of very low malaria risk (<1% parasite prevalence), but by 2016, this had dropped to only 26.7% of the population. Meanwhile, the population in high transmission areas (prevalence >20%) increased from only 2.2% in 2011 to 9.2% in 2016. A comparison of the model-based estimates with the raw MIS results indicates there was an underestimation of the situation in 2016, since the raw figures likely associated with survey timings were delayed until after the peak transmission season.

CONCLUSIONS

Malaria remains an important health problem in Madagascar. The monthly and annual prevalence maps developed here provide a way to evaluate the magnitude of change over time, taking into account variability in survey input data. These methods can contribute to monitoring sub-national trends of malaria prevalence in Madagascar as the country aims for geographically progressive elimination.

Mapping under-5 and neonatal mortality in Africa, 2000-15: a baseline analysis for the Sustainable Development Goals

BACKGROUND

During the Millennium Development Goal (MDG) era, many countries in Africa achieved marked reductions in under-5 and neonatal mortality. Yet the pace of progress toward these goals substantially varied at the national level, demonstrating an essential need for tracking even more local trends in child mortality. With the adoption of the Sustainable Development Goals (SDGs) in 2015, which established ambitious targets for improving child survival by 2030, optimal intervention planning and targeting will require understanding of trends and rates of progress at a higher spatial resolution. In this study, we aimed to generate high-resolution estimates of under-5 and neonatal all-cause mortality across 46 countries in Africa.

METHODS

We assembled 235 geographically resolved household survey and census data sources on child deaths to produce estimates of under-5 and neonatal mortality at a resolution of 5 × 5 km grid cells across 46 African countries for 2000, 2005, 2010, and 2015. We used a Bayesian geostatistical analytical framework to generate these estimates, and implemented predictive validity tests. In addition to reporting 5 × 5 km estimates, we also aggregated results obtained from these estimates into three different levels-national, and subnational administrative levels 1 and 2-to provide the full range of geospatial resolution that local, national, and global decision makers might require.

FINDINGS

Amid improving child survival in Africa, there was substantial heterogeneity in absolute levels of under-5 and neonatal mortality in 2015, as well as the annualised rates of decline achieved from 2000 to 2015. Subnational areas in countries such as Botswana, Rwanda, and Ethiopia recorded some of the largest decreases in child mortality rates since 2000, positioning them well to achieve SDG targets by 2030 or earlier. Yet these places were the exception for Africa, since many areas, particularly in central and western Africa, must reduce under-5 mortality rates by at least 8·8% per year, between 2015 and 2030, to achieve the SDG 3.2 target for under-5 mortality by 2030.

INTERPRETATION

In the absence of unprecedented political commitment, financial support, and medical advances, the viability of SDG 3.2 achievement in Africa is precarious at best. By producing under-5 and neonatal mortality rates at multiple levels of geospatial resolution over time, this study provides key information for decision makers to target interventions at populations in the greatest need. In an era when precision public health increasingly has the potential to transform the design, implementation, and impact of health programmes, our 5 × 5 km estimates of child mortality in Africa provide a baseline against which local, national, and global stakeholders can map the pathways for ending preventable child deaths by 2030.

FUNDING

Bill & Melinda Gates Foundation.

Estimating Geographical Variation in the Risk of Zoonotic Plasmodium knowlesi Infection in Countries Eliminating Malaria

BACKGROUND

Infection by the simian malaria parasite, Plasmodium knowlesi, can lead to severe and fatal disease in humans, and is the most common cause of malaria in parts of Malaysia. Despite being a serious public health concern, the geographical distribution of P. knowlesi malaria risk is poorly understood because the parasite is often misidentified as one of the human malarias. Human cases have been confirmed in at least nine Southeast Asian countries, many of which are making progress towards eliminating the human malarias. Understanding the geographical distribution of P. knowlesi is important for identifying areas where malaria transmission will continue after the human malarias have been eliminated.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 439 records of P. knowlesi infections in humans, macaque reservoir and vector species were collated. To predict spatial variation in disease risk, a model was fitted using records from countries where the infection data coverage is high. Predictions were then made throughout Southeast Asia, including regions where infection data are sparse. The resulting map predicts areas of high risk for P. knowlesi infection in a number of countries that are forecast to be malaria-free by 2025 (Malaysia, Cambodia, Thailand and Vietnam) as well as countries projected to be eliminating malaria (Myanmar, Laos, Indonesia and the Philippines).

CONCLUSIONS/SIGNIFICANCE

We have produced the first map of P. knowlesi malaria risk, at a fine-scale resolution, to identify priority areas for surveillance based on regions with sparse data and high estimated risk. Our map provides an initial evidence base to better understand the spatial distribution of this disease and its potential wider contribution to malaria incidence. Considering malaria elimination goals, areas for prioritised surveillance are identified.

Predicting the geographical distributions of the macaque hosts and mosquito vectors of Plasmodium knowlesi malaria in forested and non-forested areas

BACKGROUND

Plasmodium knowlesi is a zoonotic pathogen, transmitted among macaques and to humans by anopheline mosquitoes. Information on P. knowlesi malaria is lacking in most regions so the first step to understand the geographical distribution of disease risk is to define the distributions of the reservoir and vector species.

METHODS

We used macaque and mosquito species presence data, background data that captured sampling bias in the presence data, a boosted regression tree model and environmental datasets, including annual data for land classes, to predict the distributions of each vector and host species. We then compared the predicted distribution of each species with cover of each land class.

RESULTS

Fine-scale distribution maps were generated for three macaque host species (Macaca fascicularis, M. nemestrina and M. leonina) and two mosquito vector complexes (the Dirus Complex and the Leucosphyrus Complex). The Leucosphyrus Complex was predicted to occur in areas with disturbed, but not intact, forest cover (> 60% tree cover) whereas the Dirus Complex was predicted to occur in areas with 10-100% tree cover as well as vegetation mosaics and cropland. Of the macaque species, M. nemestrina was mainly predicted to occur in forested areas whereas M. fascicularis was predicted to occur in vegetation mosaics, cropland, wetland and urban areas in addition to forested areas.

CONCLUSIONS

The predicted M. fascicularis distribution encompassed a wide range of habitats where humans are found. This is of most significance in the northern part of its range where members of the Dirus Complex are the main P. knowlesi vectors because these mosquitoes were also predicted to occur in a wider range of habitats. Our results support the hypothesis that conversion of intact forest into disturbed forest (for example plantations or timber concessions), or the creation of vegetation mosaics, will increase the probability that members of the Leucosphyrus Complex occur at these locations, as well as bringing humans into these areas. An explicit analysis of disease risk itself using infection data is required to explore this further. The species distributions generated here can now be included in future analyses of P. knowlesi infection risk.

Treatment-seeking rates in malaria endemic countries

Background

The proportion of individuals who seek treatment for fever is an important quantity in understanding access to and use of health systems, as well as for interpreting data on disease incidence from routine surveillance systems. For many malaria endemic countries (MECs), treatment-seeking information is available from national household surveys. The aim of this paper was to assemble sub-national estimates of treatment-seeking behaviours and to predict national treatment-seeking measures for all MECs lacking household survey data.

Methods

Data on treatment seeking for fever were obtained from Demographic and Health Surveys, Malaria Indicator Surveys and Multiple Cluster Indicator Surveys for every MEC and year that data were available. National-level social, economic and health-related variables were gathered from the World Bank as putative covariates of treatment-seeking rates. A generalized additive mixed model (GAMM) was used to estimate treatment-seeking behaviours for countries where survey data were unavailable. Two separate models were developed to predict the proportion of fever cases that would seek treatment at (1) a public health facility or (2) from any kind of treatment provider.

Results

Treatment-seeking data were available for 74 MECs and modelled for the remaining 24. GAMMs found that the percentage of pregnant women receiving prenatal care, vaccination rates, education level, government health expenditure, and GDP growth were important predictors for both categories of treatment-seeking outcomes. Treatment-seeking rates, which varied both within and among regions, revealed that public facilities were not always the primary facility type used.

Conclusions

Estimates of treatment-seeking rates show how health services are utilized and help correct reported malaria case numbers to obtain more accurate measures of disease burden. The assembled and modelled data demonstrated that while treatment-seeking rates have overall increased over time, access remains low in some malaria endemic regions and utilization of government services is in some areas limited.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-1048-x) contains supplementary material, which is available to authorized users.