Global burden of 288 causes of death and life expectancy decomposition in 204 countries and territories and 811 subnational locations, 1990-2021: a systematic analysis for the Global Burden of Disease Study 2021

BACKGROUND

Regular, detailed reporting on population health by underlying cause of death is fundamental for public health decision making. Cause-specific estimates of mortality and the subsequent effects on life expectancy worldwide are valuable metrics to gauge progress in reducing mortality rates. These estimates are particularly important following large-scale mortality spikes, such as the COVID-19 pandemic. When systematically analysed, mortality rates and life expectancy allow comparisons of the consequences of causes of death globally and over time, providing a nuanced understanding of the effect of these causes on global populations.

METHODS

The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2021 cause-of-death analysis estimated mortality and years of life lost (YLLs) from 288 causes of death by age-sex-location-year in 204 countries and territories and 811 subnational locations for each year from 1990 until 2021. The analysis used 56 604 data sources, including data from vital registration and verbal autopsy as well as surveys, censuses, surveillance systems, and cancer registries, among others. As with previous GBD rounds, cause-specific death rates for most causes were estimated using the Cause of Death Ensemble model-a modelling tool developed for GBD to assess the out-of-sample predictive validity of different statistical models and covariate permutations and combine those results to produce cause-specific mortality estimates-with alternative strategies adapted to model causes with insufficient data, substantial changes in reporting over the study period, or unusual epidemiology. YLLs were computed as the product of the number of deaths for each cause-age-sex-location-year and the standard life expectancy at each age. As part of the modelling process, uncertainty intervals (UIs) were generated using the 2·5th and 97·5th percentiles from a 1000-draw distribution for each metric. We decomposed life expectancy by cause of death, location, and year to show cause-specific effects on life expectancy from 1990 to 2021. We also used the coefficient of variation and the fraction of population affected by 90% of deaths to highlight concentrations of mortality. Findings are reported in counts and age-standardised rates. Methodological improvements for cause-of-death estimates in GBD 2021 include the expansion of under-5-years age group to include four new age groups, enhanced methods to account for stochastic variation of sparse data, and the inclusion of COVID-19 and other pandemic-related mortality-which includes excess mortality associated with the pandemic, excluding COVID-19, lower respiratory infections, measles, malaria, and pertussis. For this analysis, 199 new country-years of vital registration cause-of-death data, 5 country-years of surveillance data, 21 country-years of verbal autopsy data, and 94 country-years of other data types were added to those used in previous GBD rounds.

FINDINGS

The leading causes of age-standardised deaths globally were the same in 2019 as they were in 1990; in descending order, these were, ischaemic heart disease, stroke, chronic obstructive pulmonary disease, and lower respiratory infections. In 2021, however, COVID-19 replaced stroke as the second-leading age-standardised cause of death, with 94·0 deaths (95% UI 89·2-100·0) per 100 000 population. The COVID-19 pandemic shifted the rankings of the leading five causes, lowering stroke to the third-leading and chronic obstructive pulmonary disease to the fourth-leading position. In 2021, the highest age-standardised death rates from COVID-19 occurred in sub-Saharan Africa (271·0 deaths [250·1-290·7] per 100 000 population) and Latin America and the Caribbean (195·4 deaths [182·1-211·4] per 100 000 population). The lowest age-standardised death rates from COVID-19 were in the high-income super-region (48·1 deaths [47·4-48·8] per 100 000 population) and southeast Asia, east Asia, and Oceania (23·2 deaths [16·3-37·2] per 100 000 population). Globally, life expectancy steadily improved between 1990 and 2019 for 18 of the 22 investigated causes. Decomposition of global and regional life expectancy showed the positive effect that reductions in deaths from enteric infections, lower respiratory infections, stroke, and neonatal deaths, among others have contributed to improved survival over the study period. However, a net reduction of 1·6 years occurred in global life expectancy between 2019 and 2021, primarily due to increased death rates from COVID-19 and other pandemic-related mortality. Life expectancy was highly variable between super-regions over the study period, with southeast Asia, east Asia, and Oceania gaining 8·3 years (6·7-9·9) overall, while having the smallest reduction in life expectancy due to COVID-19 (0·4 years). The largest reduction in life expectancy due to COVID-19 occurred in Latin America and the Caribbean (3·6 years). Additionally, 53 of the 288 causes of death were highly concentrated in locations with less than 50% of the global population as of 2021, and these causes of death became progressively more concentrated since 1990, when only 44 causes showed this pattern. The concentration phenomenon is discussed heuristically with respect to enteric and lower respiratory infections, malaria, HIV/AIDS, neonatal disorders, tuberculosis, and measles.

INTERPRETATION

Long-standing gains in life expectancy and reductions in many of the leading causes of death have been disrupted by the COVID-19 pandemic, the adverse effects of which were spread unevenly among populations. Despite the pandemic, there has been continued progress in combatting several notable causes of death, leading to improved global life expectancy over the study period. Each of the seven GBD super-regions showed an overall improvement from 1990 and 2021, obscuring the negative effect in the years of the pandemic. Additionally, our findings regarding regional variation in causes of death driving increases in life expectancy hold clear policy utility. Analyses of shifting mortality trends reveal that several causes, once widespread globally, are now increasingly concentrated geographically. These changes in mortality concentration, alongside further investigation of changing risks, interventions, and relevant policy, present an important opportunity to deepen our understanding of mortality-reduction strategies. Examining patterns in mortality concentration might reveal areas where successful public health interventions have been implemented. Translating these successes to locations where certain causes of death remain entrenched can inform policies that work to improve life expectancy for people everywhere.

FUNDING

Bill & Melinda Gates Foundation.

Burden of disease scenarios for 204 countries and territories, 2022-2050: a forecasting analysis for the Global Burden of Disease Study 2021

BACKGROUND

Future trends in disease burden and drivers of health are of great interest to policy makers and the public at large. This information can be used for policy and long-term health investment, planning, and prioritisation. We have expanded and improved upon previous forecasts produced as part of the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) and provide a reference forecast (the most likely future), and alternative scenarios assessing disease burden trajectories if selected sets of risk factors were eliminated from current levels by 2050.

METHODS

Using forecasts of major drivers of health such as the Socio-demographic Index (SDI; a composite measure of lag-distributed income per capita, mean years of education, and total fertility under 25 years of age) and the full set of risk factor exposures captured by GBD, we provide cause-specific forecasts of mortality, years of life lost (YLLs), years lived with disability (YLDs), and disability-adjusted life-years (DALYs) by age and sex from 2022 to 2050 for 204 countries and territories, 21 GBD regions, seven super-regions, and the world. All analyses were done at the cause-specific level so that only risk factors deemed causal by the GBD comparative risk assessment influenced future trajectories of mortality for each disease. Cause-specific mortality was modelled using mixed-effects models with SDI and time as the main covariates, and the combined impact of causal risk factors as an offset in the model. At the all-cause mortality level, we captured unexplained variation by modelling residuals with an autoregressive integrated moving average model with drift attenuation. These all-cause forecasts constrained the cause-specific forecasts at successively deeper levels of the GBD cause hierarchy using cascading mortality models, thus ensuring a robust estimate of cause-specific mortality. For non-fatal measures (eg, low back pain), incidence and prevalence were forecasted from mixed-effects models with SDI as the main covariate, and YLDs were computed from the resulting prevalence forecasts and average disability weights from GBD. Alternative future scenarios were constructed by replacing appropriate reference trajectories for risk factors with hypothetical trajectories of gradual elimination of risk factor exposure from current levels to 2050. The scenarios were constructed from various sets of risk factors: environmental risks (Safer Environment scenario), risks associated with communicable, maternal, neonatal, and nutritional diseases (CMNNs; Improved Childhood Nutrition and Vaccination scenario), risks associated with major non-communicable diseases (NCDs; Improved Behavioural and Metabolic Risks scenario), and the combined effects of these three scenarios. Using the Shared Socioeconomic Pathways climate scenarios SSP2-4.5 as reference and SSP1-1.9 as an optimistic alternative in the Safer Environment scenario, we accounted for climate change impact on health by using the most recent Intergovernmental Panel on Climate Change temperature forecasts and published trajectories of ambient air pollution for the same two scenarios. Life expectancy and healthy life expectancy were computed using standard methods. The forecasting framework includes computing the age-sex-specific future population for each location and separately for each scenario. 95% uncertainty intervals (UIs) for each individual future estimate were derived from the 2·5th and 97·5th percentiles of distributions generated from propagating 500 draws through the multistage computational pipeline.

FINDINGS

In the reference scenario forecast, global and super-regional life expectancy increased from 2022 to 2050, but improvement was at a slower pace than in the three decades preceding the COVID-19 pandemic (beginning in 2020). Gains in future life expectancy were forecasted to be greatest in super-regions with comparatively low life expectancies (such as sub-Saharan Africa) compared with super-regions with higher life expectancies (such as the high-income super-region), leading to a trend towards convergence in life expectancy across locations between now and 2050. At the super-region level, forecasted healthy life expectancy patterns were similar to those of life expectancies. Forecasts for the reference scenario found that health will improve in the coming decades, with all-cause age-standardised DALY rates decreasing in every GBD super-region. The total DALY burden measured in counts, however, will increase in every super-region, largely a function of population ageing and growth. We also forecasted that both DALY counts and age-standardised DALY rates will continue to shift from CMNNs to NCDs, with the most pronounced shifts occurring in sub-Saharan Africa (60·1% [95% UI 56·8-63·1] of DALYs were from CMNNs in 2022 compared with 35·8% [31·0-45·0] in 2050) and south Asia (31·7% [29·2-34·1] to 15·5% [13·7-17·5]). This shift is reflected in the leading global causes of DALYs, with the top four causes in 2050 being ischaemic heart disease, stroke, diabetes, and chronic obstructive pulmonary disease, compared with 2022, with ischaemic heart disease, neonatal disorders, stroke, and lower respiratory infections at the top. The global proportion of DALYs due to YLDs likewise increased from 33·8% (27·4-40·3) to 41·1% (33·9-48·1) from 2022 to 2050, demonstrating an important shift in overall disease burden towards morbidity and away from premature death. The largest shift of this kind was forecasted for sub-Saharan Africa, from 20·1% (15·6-25·3) of DALYs due to YLDs in 2022 to 35·6% (26·5-43·0) in 2050. In the assessment of alternative future scenarios, the combined effects of the scenarios (Safer Environment, Improved Childhood Nutrition and Vaccination, and Improved Behavioural and Metabolic Risks scenarios) demonstrated an important decrease in the global burden of DALYs in 2050 of 15·4% (13·5-17·5) compared with the reference scenario, with decreases across super-regions ranging from 10·4% (9·7-11·3) in the high-income super-region to 23·9% (20·7-27·3) in north Africa and the Middle East. The Safer Environment scenario had its largest decrease in sub-Saharan Africa (5·2% [3·5-6·8]), the Improved Behavioural and Metabolic Risks scenario in north Africa and the Middle East (23·2% [20·2-26·5]), and the Improved Nutrition and Vaccination scenario in sub-Saharan Africa (2·0% [-0·6 to 3·6]).

INTERPRETATION

Globally, life expectancy and age-standardised disease burden were forecasted to improve between 2022 and 2050, with the majority of the burden continuing to shift from CMNNs to NCDs. That said, continued progress on reducing the CMNN disease burden will be dependent on maintaining investment in and policy emphasis on CMNN disease prevention and treatment. Mostly due to growth and ageing of populations, the number of deaths and DALYs due to all causes combined will generally increase. By constructing alternative future scenarios wherein certain risk exposures are eliminated by 2050, we have shown that opportunities exist to substantially improve health outcomes in the future through concerted efforts to prevent exposure to well established risk factors and to expand access to key health interventions.

FUNDING

Bill & Melinda Gates Foundation.

Global fertility in 204 countries and territories, 1950-2021, with forecasts to 2100: a comprehensive demographic analysis for the Global Burden of Disease Study 2021

BACKGROUND

Accurate assessments of current and future fertility-including overall trends and changing population age structures across countries and regions-are essential to help plan for the profound social, economic, environmental, and geopolitical challenges that these changes will bring. Estimates and projections of fertility are necessary to inform policies involving resource and health-care needs, labour supply, education, gender equality, and family planning and support. The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2021 produced up-to-date and comprehensive demographic assessments of key fertility indicators at global, regional, and national levels from 1950 to 2021 and forecast fertility metrics to 2100 based on a reference scenario and key policy-dependent alternative scenarios.

METHODS

To estimate fertility indicators from 1950 to 2021, mixed-effects regression models and spatiotemporal Gaussian process regression were used to synthesise data from 8709 country-years of vital and sample registrations, 1455 surveys and censuses, and 150 other sources, and to generate age-specific fertility rates (ASFRs) for 5-year age groups from age 10 years to 54 years. ASFRs were summed across age groups to produce estimates of total fertility rate (TFR). Livebirths were calculated by multiplying ASFR and age-specific female population, then summing across ages 10-54 years. To forecast future fertility up to 2100, our Institute for Health Metrics and Evaluation (IHME) forecasting model was based on projections of completed cohort fertility at age 50 years (CCF50; the average number of children born over time to females from a specified birth cohort), which yields more stable and accurate measures of fertility than directly modelling TFR. CCF50 was modelled using an ensemble approach in which three sub-models (with two, three, and four covariates variously consisting of female educational attainment, contraceptive met need, population density in habitable areas, and under-5 mortality) were given equal weights, and analyses were conducted utilising the MR-BRT (meta-regression-Bayesian, regularised, trimmed) tool. To capture time-series trends in CCF50 not explained by these covariates, we used a first-order autoregressive model on the residual term. CCF50 as a proportion of each 5-year ASFR was predicted using a linear mixed-effects model with fixed-effects covariates (female educational attainment and contraceptive met need) and random intercepts for geographical regions. Projected TFRs were then computed for each calendar year as the sum of single-year ASFRs across age groups. The reference forecast is our estimate of the most likely fertility future given the model, past fertility, forecasts of covariates, and historical relationships between covariates and fertility. We additionally produced forecasts for multiple alternative scenarios in each location: the UN Sustainable Development Goal (SDG) for education is achieved by 2030; the contraceptive met need SDG is achieved by 2030; pro-natal policies are enacted to create supportive environments for those who give birth; and the previous three scenarios combined. Uncertainty from past data inputs and model estimation was propagated throughout analyses by taking 1000 draws for past and present fertility estimates and 500 draws for future forecasts from the estimated distribution for each metric, with 95% uncertainty intervals (UIs) given as the 2·5 and 97·5 percentiles of the draws. To evaluate the forecasting performance of our model and others, we computed skill values-a metric assessing gain in forecasting accuracy-by comparing predicted versus observed ASFRs from the past 15 years (2007-21). A positive skill metric indicates that the model being evaluated performs better than the baseline model (here, a simplified model holding 2007 values constant in the future), and a negative metric indicates that the evaluated model performs worse than baseline.

FINDINGS

During the period from 1950 to 2021, global TFR more than halved, from 4·84 (95% UI 4·63-5·06) to 2·23 (2·09-2·38). Global annual livebirths peaked in 2016 at 142 million (95% UI 137-147), declining to 129 million (121-138) in 2021. Fertility rates declined in all countries and territories since 1950, with TFR remaining above 2·1-canonically considered replacement-level fertility-in 94 (46·1%) countries and territories in 2021. This included 44 of 46 countries in sub-Saharan Africa, which was the super-region with the largest share of livebirths in 2021 (29·2% [28·7-29·6]). 47 countries and territories in which lowest estimated fertility between 1950 and 2021 was below replacement experienced one or more subsequent years with higher fertility; only three of these locations rebounded above replacement levels. Future fertility rates were projected to continue to decline worldwide, reaching a global TFR of 1·83 (1·59-2·08) in 2050 and 1·59 (1·25-1·96) in 2100 under the reference scenario. The number of countries and territories with fertility rates remaining above replacement was forecast to be 49 (24·0%) in 2050 and only six (2·9%) in 2100, with three of these six countries included in the 2021 World Bank-defined low-income group, all located in the GBD super-region of sub-Saharan Africa. The proportion of livebirths occurring in sub-Saharan Africa was forecast to increase to more than half of the world's livebirths in 2100, to 41·3% (39·6-43·1) in 2050 and 54·3% (47·1-59·5) in 2100. The share of livebirths was projected to decline between 2021 and 2100 in most of the six other super-regions-decreasing, for example, in south Asia from 24·8% (23·7-25·8) in 2021 to 16·7% (14·3-19·1) in 2050 and 7·1% (4·4-10·1) in 2100-but was forecast to increase modestly in the north Africa and Middle East and high-income super-regions. Forecast estimates for the alternative combined scenario suggest that meeting SDG targets for education and contraceptive met need, as well as implementing pro-natal policies, would result in global TFRs of 1·65 (1·40-1·92) in 2050 and 1·62 (1·35-1·95) in 2100. The forecasting skill metric values for the IHME model were positive across all age groups, indicating that the model is better than the constant prediction.

INTERPRETATION

Fertility is declining globally, with rates in more than half of all countries and territories in 2021 below replacement level. Trends since 2000 show considerable heterogeneity in the steepness of declines, and only a small number of countries experienced even a slight fertility rebound after their lowest observed rate, with none reaching replacement level. Additionally, the distribution of livebirths across the globe is shifting, with a greater proportion occurring in the lowest-income countries. Future fertility rates will continue to decline worldwide and will remain low even under successful implementation of pro-natal policies. These changes will have far-reaching economic and societal consequences due to ageing populations and declining workforces in higher-income countries, combined with an increasing share of livebirths among the already poorest regions of the world.

FUNDING

Bill & Melinda Gates Foundation.

Global incidence, prevalence, years lived with disability (YLDs), disability-adjusted life-years (DALYs), and healthy life expectancy (HALE) for 371 diseases and injuries in 204 countries and territories and 811 subnational locations, 1990-2021: a systematic analysis for the Global Burden of Disease Study 2021

BACKGROUND

Detailed, comprehensive, and timely reporting on population health by underlying causes of disability and premature death is crucial to understanding and responding to complex patterns of disease and injury burden over time and across age groups, sexes, and locations. The availability of disease burden estimates can promote evidence-based interventions that enable public health researchers, policy makers, and other professionals to implement strategies that can mitigate diseases. It can also facilitate more rigorous monitoring of progress towards national and international health targets, such as the Sustainable Development Goals. For three decades, the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) has filled that need. A global network of collaborators contributed to the production of GBD 2021 by providing, reviewing, and analysing all available data. GBD estimates are updated routinely with additional data and refined analytical methods. GBD 2021 presents, for the first time, estimates of health loss due to the COVID-19 pandemic.

METHODS

The GBD 2021 disease and injury burden analysis estimated years lived with disability (YLDs), years of life lost (YLLs), disability-adjusted life-years (DALYs), and healthy life expectancy (HALE) for 371 diseases and injuries using 100 983 data sources. Data were extracted from vital registration systems, verbal autopsies, censuses, household surveys, disease-specific registries, health service contact data, and other sources. YLDs were calculated by multiplying cause-age-sex-location-year-specific prevalence of sequelae by their respective disability weights, for each disease and injury. YLLs were calculated by multiplying cause-age-sex-location-year-specific deaths by the standard life expectancy at the age that death occurred. DALYs were calculated by summing YLDs and YLLs. HALE estimates were produced using YLDs per capita and age-specific mortality rates by location, age, sex, year, and cause. 95% uncertainty intervals (UIs) were generated for all final estimates as the 2·5th and 97·5th percentiles values of 500 draws. Uncertainty was propagated at each step of the estimation process. Counts and age-standardised rates were calculated globally, for seven super-regions, 21 regions, 204 countries and territories (including 21 countries with subnational locations), and 811 subnational locations, from 1990 to 2021. Here we report data for 2010 to 2021 to highlight trends in disease burden over the past decade and through the first 2 years of the COVID-19 pandemic.

FINDINGS

Global DALYs increased from 2·63 billion (95% UI 2·44-2·85) in 2010 to 2·88 billion (2·64-3·15) in 2021 for all causes combined. Much of this increase in the number of DALYs was due to population growth and ageing, as indicated by a decrease in global age-standardised all-cause DALY rates of 14·2% (95% UI 10·7-17·3) between 2010 and 2019. Notably, however, this decrease in rates reversed during the first 2 years of the COVID-19 pandemic, with increases in global age-standardised all-cause DALY rates since 2019 of 4·1% (1·8-6·3) in 2020 and 7·2% (4·7-10·0) in 2021. In 2021, COVID-19 was the leading cause of DALYs globally (212·0 million [198·0-234·5] DALYs), followed by ischaemic heart disease (188·3 million [176·7-198·3]), neonatal disorders (186·3 million [162·3-214·9]), and stroke (160·4 million [148·0-171·7]). However, notable health gains were seen among other leading communicable, maternal, neonatal, and nutritional (CMNN) diseases. Globally between 2010 and 2021, the age-standardised DALY rates for HIV/AIDS decreased by 47·8% (43·3-51·7) and for diarrhoeal diseases decreased by 47·0% (39·9-52·9). Non-communicable diseases contributed 1·73 billion (95% UI 1·54-1·94) DALYs in 2021, with a decrease in age-standardised DALY rates since 2010 of 6·4% (95% UI 3·5-9·5). Between 2010 and 2021, among the 25 leading Level 3 causes, age-standardised DALY rates increased most substantially for anxiety disorders (16·7% [14·0-19·8]), depressive disorders (16·4% [11·9-21·3]), and diabetes (14·0% [10·0-17·4]). Age-standardised DALY rates due to injuries decreased globally by 24·0% (20·7-27·2) between 2010 and 2021, although improvements were not uniform across locations, ages, and sexes. Globally, HALE at birth improved slightly, from 61·3 years (58·6-63·6) in 2010 to 62·2 years (59·4-64·7) in 2021. However, despite this overall increase, HALE decreased by 2·2% (1·6-2·9) between 2019 and 2021.

INTERPRETATION

Putting the COVID-19 pandemic in the context of a mutually exclusive and collectively exhaustive list of causes of health loss is crucial to understanding its impact and ensuring that health funding and policy address needs at both local and global levels through cost-effective and evidence-based interventions. A global epidemiological transition remains underway. Our findings suggest that prioritising non-communicable disease prevention and treatment policies, as well as strengthening health systems, continues to be crucially important. The progress on reducing the burden of CMNN diseases must not stall; although global trends are improving, the burden of CMNN diseases remains unacceptably high. Evidence-based interventions will help save the lives of young children and mothers and improve the overall health and economic conditions of societies across the world. Governments and multilateral organisations should prioritise pandemic preparedness planning alongside efforts to reduce the burden of diseases and injuries that will strain resources in the coming decades.

FUNDING

Bill & Melinda Gates Foundation.

Global age-sex-specific mortality, life expectancy, and population estimates in 204 countries and territories and 811 subnational locations, 1950-2021, and the impact of the COVID-19 pandemic: a comprehensive demographic analysis for the Global Burden of Disease Study 2021

BACKGROUND

Estimates of demographic metrics are crucial to assess levels and trends of population health outcomes. The profound impact of the COVID-19 pandemic on populations worldwide has underscored the need for timely estimates to understand this unprecedented event within the context of long-term population health trends. The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2021 provides new demographic estimates for 204 countries and territories and 811 additional subnational locations from 1950 to 2021, with a particular emphasis on changes in mortality and life expectancy that occurred during the 2020-21 COVID-19 pandemic period.

METHODS

22 223 data sources from vital registration, sample registration, surveys, censuses, and other sources were used to estimate mortality, with a subset of these sources used exclusively to estimate excess mortality due to the COVID-19 pandemic. 2026 data sources were used for population estimation. Additional sources were used to estimate migration; the effects of the HIV epidemic; and demographic discontinuities due to conflicts, famines, natural disasters, and pandemics, which are used as inputs for estimating mortality and population. Spatiotemporal Gaussian process regression (ST-GPR) was used to generate under-5 mortality rates, which synthesised 30 763 location-years of vital registration and sample registration data, 1365 surveys and censuses, and 80 other sources. ST-GPR was also used to estimate adult mortality (between ages 15 and 59 years) based on information from 31 642 location-years of vital registration and sample registration data, 355 surveys and censuses, and 24 other sources. Estimates of child and adult mortality rates were then used to generate life tables with a relational model life table system. For countries with large HIV epidemics, life tables were adjusted using independent estimates of HIV-specific mortality generated via an epidemiological analysis of HIV prevalence surveys, antenatal clinic serosurveillance, and other data sources. Excess mortality due to the COVID-19 pandemic in 2020 and 2021 was determined by subtracting observed all-cause mortality (adjusted for late registration and mortality anomalies) from the mortality expected in the absence of the pandemic. Expected mortality was calculated based on historical trends using an ensemble of models. In location-years where all-cause mortality data were unavailable, we estimated excess mortality rates using a regression model with covariates pertaining to the pandemic. Population size was computed using a Bayesian hierarchical cohort component model. Life expectancy was calculated using age-specific mortality rates and standard demographic methods. Uncertainty intervals (UIs) were calculated for every metric using the 25th and 975th ordered values from a 1000-draw posterior distribution.

FINDINGS

Global all-cause mortality followed two distinct patterns over the study period: age-standardised mortality rates declined between 1950 and 2019 (a 62·8% [95% UI 60·5-65·1] decline), and increased during the COVID-19 pandemic period (2020-21; 5·1% [0·9-9·6] increase). In contrast with the overall reverse in mortality trends during the pandemic period, child mortality continued to decline, with 4·66 million (3·98-5·50) global deaths in children younger than 5 years in 2021 compared with 5·21 million (4·50-6·01) in 2019. An estimated 131 million (126-137) people died globally from all causes in 2020 and 2021 combined, of which 15·9 million (14·7-17·2) were due to the COVID-19 pandemic (measured by excess mortality, which includes deaths directly due to SARS-CoV-2 infection and those indirectly due to other social, economic, or behavioural changes associated with the pandemic). Excess mortality rates exceeded 150 deaths per 100 000 population during at least one year of the pandemic in 80 countries and territories, whereas 20 nations had a negative excess mortality rate in 2020 or 2021, indicating that all-cause mortality in these countries was lower during the pandemic than expected based on historical trends. Between 1950 and 2021, global life expectancy at birth increased by 22·7 years (20·8-24·8), from 49·0 years (46·7-51·3) to 71·7 years (70·9-72·5). Global life expectancy at birth declined by 1·6 years (1·0-2·2) between 2019 and 2021, reversing historical trends. An increase in life expectancy was only observed in 32 (15·7%) of 204 countries and territories between 2019 and 2021. The global population reached 7·89 billion (7·67-8·13) people in 2021, by which time 56 of 204 countries and territories had peaked and subsequently populations have declined. The largest proportion of population growth between 2020 and 2021 was in sub-Saharan Africa (39·5% [28·4-52·7]) and south Asia (26·3% [9·0-44·7]). From 2000 to 2021, the ratio of the population aged 65 years and older to the population aged younger than 15 years increased in 188 (92·2%) of 204 nations.

INTERPRETATION

Global adult mortality rates markedly increased during the COVID-19 pandemic in 2020 and 2021, reversing past decreasing trends, while child mortality rates continued to decline, albeit more slowly than in earlier years. Although COVID-19 had a substantial impact on many demographic indicators during the first 2 years of the pandemic, overall global health progress over the 72 years evaluated has been profound, with considerable improvements in mortality and life expectancy. Additionally, we observed a deceleration of global population growth since 2017, despite steady or increasing growth in lower-income countries, combined with a continued global shift of population age structures towards older ages. These demographic changes will likely present future challenges to health systems, economies, and societies. The comprehensive demographic estimates reported here will enable researchers, policy makers, health practitioners, and other key stakeholders to better understand and address the profound changes that have occurred in the global health landscape following the first 2 years of the COVID-19 pandemic, and longer-term trends beyond the pandemic.

FUNDING

Bill & Melinda Gates Foundation.

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.

Impacts on Human Movement in Australian Cities Related to the COVID-19 Pandemic

No studies have yet examined high-resolution shifts in the spatial patterns of human movement in Australia throughout 2020 and 2021, a period coincident with the repeated enactment and removal of varied governmental restrictions aimed at reducing community transmission of SARS-CoV-2. We compared overlapping timeseries of COVID-19 pandemic-related restrictions, epidemiological data on cases and vaccination rates, and high-resolution human movement data to characterize population-level responses to the pandemic in Australian cities. We found that restrictions on human movement and/or mandatory business closures reduced the average population-level weekly movement volumes in cities, as measured by aggregated travel time, by almost half. Of the movements that continued to occur, long movements reduced more dramatically than short movements, likely indicating that people stayed closer to home. We also found that the repeated lockdowns did not reduce their impact on human movement, but the effect of the restrictions on human movement waned as the duration of restrictions increased. Lastly, we found that after restrictions ceased, the subsequent surge in SARS-CoV-2 transmission coincided with a substantial, non-mandated drop in human movement volume. These findings have implications for public health policy makers when faced with anticipating responses to restrictions during future emergency situations.

Mapping tuberculosis prevalence in Ethiopia using geospatial meta-analysis

BACKGROUND

Reliable and detailed data on the prevalence of tuberculosis (TB) with sub-national estimates are scarce in Ethiopia. We address this knowledge gap by spatially predicting the national, sub-national and local prevalence of TB, and identifying drivers of TB prevalence across the country.

METHODS

TB prevalence data were obtained from the Ethiopia national TB prevalence survey and from a comprehensive review of published reports. Geospatial covariates were obtained from publicly available sources. A random effects meta-analysis was used to estimate a pooled prevalence of TB at the national level, and model-based geostatistics were used to estimate the spatial variation of TB prevalence at sub-national and local levels. Within the MBG Plugin Framework, a logistic regression model was fitted to TB prevalence data using both fixed covariate effects and spatial random effects to identify drivers of TB and to predict the prevalence of TB.

RESULTS

The overall pooled prevalence of TB in Ethiopia was 0.19% [95% confidence intervals (CI): 0.12%-0.28%]. There was a high degree of heterogeneity in the prevalence of TB (I2 96.4%, P <0.001), which varied by geographical locations, data collection periods and diagnostic methods. The highest prevalence of TB was observed in Dire Dawa (0.96%), Gambela (0.88%), Somali (0.42%), Addis Ababa (0.28%) and Afar (0.24%) regions. Nationally, there was a decline in TB prevalence from 0.18% in 2001 to 0.04% in 2009. However, prevalence increased back to 0.29% in 2014. Substantial spatial variation of TB prevalence was observed at a regional level, with a higher prevalence observed in the border regions, and at a local level within regions. The spatial distribution of TB prevalence was positively associated with population density.

CONCLUSION

The results of this study showed that TB prevalence varied substantially at sub-national and local levels in Ethiopia. Spatial patterns were associated with population density. These results suggest that targeted interventions in high-risk areas may reduce the burden of TB in Ethiopia and additional data collection would be required to make further inferences on TB prevalence in areas that lack data.

An archetypes approach to malaria intervention impact mapping: a new framework and example application

BACKGROUND

As both mechanistic and geospatial malaria modeling methods become more integrated into malaria policy decisions, there is increasing demand for strategies that combine these two methods. This paper introduces a novel archetypes-based methodology for generating high-resolution intervention impact maps based on mechanistic model simulations. An example configuration of the framework is described and explored.

METHODS

First, dimensionality reduction and clustering techniques were applied to rasterized geospatial environmental and mosquito covariates to find archetypal malaria transmission patterns. Next, mechanistic models were run on a representative site from each archetype to assess intervention impact. Finally, these mechanistic results were reprojected onto each pixel to generate full maps of intervention impact. The example configuration used ERA5 and Malaria Atlas Project covariates, singular value decomposition, k-means clustering, and the Institute for Disease Modeling's EMOD model to explore a range of three-year malaria interventions primarily focused on vector control and case management.

RESULTS

Rainfall, temperature, and mosquito abundance layers were clustered into ten transmission archetypes with distinct properties. Example intervention impact curves and maps highlighted archetype-specific variation in efficacy of vector control interventions. A sensitivity analysis showed that the procedure for selecting representative sites to simulate worked well in all but one archetype.

CONCLUSION

This paper introduces a novel methodology which combines the richness of spatiotemporal mapping with the rigor of mechanistic modeling to create a multi-purpose infrastructure for answering a broad range of important questions in the malaria policy space. It is flexible and adaptable to a range of input covariates, mechanistic models, and mapping strategies and can be adapted to the modelers' setting of choice.

Evaluating COVID-19-Related Disruptions to Effective Malaria Case Management in 2020-2021 and Its Potential Effects on Malaria Burden in Sub-Saharan Africa

The COVID-19 pandemic has led to far-reaching disruptions to health systems, including preventative and curative services for malaria. The aim of this study was to estimate the magnitude of disruptions in malaria case management in sub-Saharan Africa and their impact on malaria burden during the COVID-19 pandemic. We used survey data collected by the World Health Organization, in which individual country stakeholders reported on the extent of disruptions to malaria diagnosis and treatment. The relative disruption values were then applied to estimates of antimalarial treatment rates and used as inputs to an established spatiotemporal Bayesian geostatistical framework to generate annual malaria burden estimates with case management disruptions. This enabled an estimation of the additional malaria burden attributable to pandemic-related impacts on treatment rates in 2020 and 2021. Our analysis found that disruptions in access to antimalarial treatment in sub-Saharan Africa likely resulted in approximately 5.9 (4.4-7.2 95% CI) million more malaria cases and 76 (20-132) thousand additional deaths in the 2020-2021 period within the study region, equivalent to approximately 1.2% (0.3-2.1 95% CI) greater clinical incidence of malaria and 8.1% (2.1-14.1 95% CI) greater malaria mortality than expected in the absence of the disruptions to malaria case management. The available evidence suggests that access to antimalarials was disrupted to a significant degree and should be considered an area of focus to avoid further escalations in malaria morbidity and mortality. The results from this analysis were used to estimate cases and deaths in the World Malaria Report 2022 during the pandemic years.

A novel statistical framework for exploring the population dynamics and seasonality of mosquito populations

Understanding the temporal dynamics of mosquito populations underlying vector-borne disease transmission is key to optimizing control strategies. Many questions remain surrounding the drivers of these dynamics and how they vary between species-questions rarely answerable from individual entomological studies (that typically focus on a single location or species). We develop a novel statistical framework enabling identification and classification of time series with similar temporal properties, and use this framework to systematically explore variation in population dynamics and seasonality in anopheline mosquito time series catch data spanning seven species, 40 years and 117 locations across mainland India. Our analyses reveal pronounced variation in dynamics across locations and between species in the extent of seasonality and timing of seasonal peaks. However, we show that these diverse dynamics can be clustered into four 'dynamical archetypes', each characterized by distinct temporal properties and associated with a largely unique set of environmental factors. Our results highlight that a range of environmental factors including rainfall, temperature, proximity to static water bodies and patterns of land use (particularly urbanicity) shape the dynamics and seasonality of mosquito populations, and provide a generically applicable framework to better identify and understand patterns of seasonal variation in vectors relevant to public health.

Spatial codistribution of HIV, tuberculosis and malaria in Ethiopia

Background

HIV, tuberculosis (TB) and malaria are the three most important infectious diseases in Ethiopia, and sub-Saharan Africa. Understanding the spatial codistribution of these diseases is critical for designing geographically targeted and integrated disease control programmes. This study investigated the spatial overlap and drivers of HIV, TB and malaria prevalence in Ethiopia.

Methods

HIV, TB and malaria data were obtained from different nationwide prevalence surveys, and geospatial covariates were obtained from publicly available sources. A Bayesian model-based geostatistical framework was applied to each survey leveraging the strength of high-resolution spatial covariates to predict continuous disease-specific prevalence surfaces and their codistribution.

Results

The national prevalence was 1.54% (95% CI 1.40 to 1.70) for HIV, 0.39% (95% CI 0.34 to 0.45) for TB and 1.1% (95%CI 0.95 to 1.32) for malaria. Substantial subnational variation was predicted with the highest HIV prevalence estimated in Gambela (4.52%), Addis Ababa (3.52%) and Dire Dawa (2.67%) regions. TB prevalence was highest in Dire Dawa (0.96%) and Gambela (0.88%), while malaria was highest in Gambela (6.1%) and Benishangul-Gumuz (3.8%). Spatial overlap of their prevalence was observed in some parts of the country, mainly Gambela region. Spatial distribution of the diseases was significantly associated with healthcare access, demographic, and climatic factors.

Conclusions

The national distribution of HIV, TB and malaria was highly focal in Ethiopia, with substantial variation at subnational and local levels. Spatial distribution of the diseases was significantly associated with healthcare access, demographic and climatic factors. Spatial overlap of HIV, TB and malaria prevalence was observed in some parts of the country. Integrated control programmes for these diseases should be targeted to these areas with high levels of co-endemicity.

A simulation study of disaggregation regression for spatial disease mapping

Disaggregation regression has become an important tool in spatial disease mapping for making fine-scale predictions of disease risk from aggregated response data. By including high resolution covariate information and modeling the data generating process on a fine scale, it is hoped that these models can accurately learn the relationships between covariates and response at a fine spatial scale. However, validating these high resolution predictions can be a challenge, as often there is no data observed at this spatial scale. In this study, disaggregation regression was performed on simulated data in various settings and the resulting fine-scale predictions are compared to the simulated ground truth. Performance was investigated with varying numbers of data points, sizes of aggregated areas and levels of model misspecification. The effectiveness of cross validation on the aggregate level as a measure of fine-scale predictive performance was also investigated. Predictive performance improved as the number of observations increased and as the size of the aggregated areas decreased. When the model was well-specified, fine-scale predictions were accurate even with small numbers of observations and large aggregated areas. Under model misspecification predictive performance was significantly worse for large aggregated areas but remained high when response data was aggregated over smaller regions. Cross-validation correlation on the aggregate level was a moderately good predictor of fine-scale predictive performance. While these simulations are unlikely to capture the nuances of real-life response data, this study gives insight into the effectiveness of disaggregation regression in different contexts.

Mapping the endemicity and seasonality of clinical malaria for intervention targeting in Haiti using routine case data

Towards the goal of malaria elimination on Hispaniola, the National Malaria Control Program of Haiti and its international partner organisations are conducting a campaign of interventions targeted to high-risk communities prioritised through evidence-based planning. Here we present a key piece of this planning: an up-to-date, fine-scale endemicity map and seasonality profile for Haiti informed by monthly case counts from 771 health facilities reporting from across the country throughout the 6-year period from January 2014 to December 2019. To this end, a novel hierarchical Bayesian modelling framework was developed in which a latent, pixel-level incidence surface with spatio-temporal innovations is linked to the observed case data via a flexible catchment sub-model designed to account for the absence of data on case household locations. These maps have focussed the delivery of indoor residual spraying and focal mass drug administration in the Grand’Anse Department in South-Western Haiti.

Global economic costs due to vivax malaria and the potential impact of its radical cure: A modelling study

BACKGROUND

In 2017, an estimated 14 million cases of Plasmodium vivax malaria were reported from Asia, Central and South America, and the Horn of Africa. The clinical burden of vivax malaria is largely driven by its ability to form dormant liver stages (hypnozoites) that can reactivate to cause recurrent episodes of malaria. Elimination of both the blood and liver stages of the parasites ("radical cure") is required to achieve a sustained clinical response and prevent ongoing transmission of the parasite. Novel treatment options and point-of-care diagnostics are now available to ensure that radical cure can be administered safely and effectively. We quantified the global economic cost of vivax malaria and estimated the potential cost benefit of a policy of radical cure after testing patients for glucose-6-phosphate dehydrogenase (G6PD) deficiency.

METHODS AND FINDINGS

Estimates of the healthcare provider and household costs due to vivax malaria were collated and combined with national case estimates for 44 endemic countries in 2017. These provider and household costs were compared with those that would be incurred under 2 scenarios for radical cure following G6PD screening: (1) complete adherence following daily supervised primaquine therapy and (2) unsupervised treatment with an assumed 40% effectiveness. A probabilistic sensitivity analysis generated credible intervals (CrIs) for the estimates. Globally, the annual cost of vivax malaria was US$359 million (95% CrI: US$222 to 563 million), attributable to 14.2 million cases of vivax malaria in 2017. From a societal perspective, adopting a policy of G6PD deficiency screening and supervision of primaquine to all eligible patients would prevent 6.1 million cases and reduce the global cost of vivax malaria to US$266 million (95% CrI: US$161 to 415 million), although healthcare provider costs would increase by US$39 million. If perfect adherence could be achieved with a single visit, then the global cost would fall further to US$225 million, equivalent to $135 million in cost savings from the baseline global costs. A policy of unsupervised primaquine reduced the cost to US$342 million (95% CrI: US$209 to 532 million) while preventing 2.1 million cases. Limitations of the study include partial availability of country-level cost data and parameter uncertainty for the proportion of patients prescribed primaquine, patient adherence to a full course of primaquine, and effectiveness of primaquine when unsupervised.

CONCLUSIONS

Our modelling study highlights a substantial global economic burden of vivax malaria that could be reduced through investment in safe and effective radical cure achieved by routine screening for G6PD deficiency and supervision of treatment. Novel, low-cost interventions for improving adherence to primaquine to ensure effective radical cure and widespread access to screening for G6PD deficiency will be critical to achieving the timely global elimination of P. vivax.

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.

Spatiotemporal mapping of malaria prevalence in Madagascar using routine surveillance and health survey data

Malaria transmission in Madagascar is highly heterogeneous, exhibiting spatial, seasonal and long-term trends. Previous efforts to map malaria risk in Madagascar used prevalence data from Malaria Indicator Surveys. These cross-sectional surveys, conducted during the high transmission season most recently in 2013 and 2016, provide nationally representative prevalence data but cover relatively short time frames. Conversely, monthly case data are collected at health facilities but suffer from biases, including incomplete reporting and low rates of treatment seeking. We combined survey and case data to make monthly maps of prevalence between 2013 and 2016. Health facility catchment populations were estimated to produce incidence rates from the case data. Smoothed incidence surfaces, environmental and socioeconomic covariates, and survey data informed a Bayesian prevalence model, in which a flexible incidence-to-prevalence relationship was learned. Modelled spatial trends were consistent over time, with highest prevalence in the coastal regions and low prevalence in the highlands and desert south. Prevalence was lowest in 2014 and peaked in 2015 and seasonality was widely observed, including in some lower transmission regions. These trends highlight the utility of monthly prevalence estimates over the four year period. By combining survey and case data using this two-step modelling approach, we were able to take advantage of the relative strengths of each metric while accounting for potential bias in the case data. Similar modelling approaches combining large datasets of different malaria metrics may be applicable across sub-Saharan Africa.

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.

Global burden of 87 risk factors in 204 countries and territories, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019

BACKGROUND

Rigorous analysis of levels and trends in exposure to leading risk factors and quantification of their effect on human health are important to identify where public health is making progress and in which cases current efforts are inadequate. The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2019 provides a standardised and comprehensive assessment of the magnitude of risk factor exposure, relative risk, and attributable burden of disease.

METHODS

GBD 2019 estimated attributable mortality, years of life lost (YLLs), years of life lived with disability (YLDs), and disability-adjusted life-years (DALYs) for 87 risk factors and combinations of risk factors, at the global level, regionally, and for 204 countries and territories. GBD uses a hierarchical list of risk factors so that specific risk factors (eg, sodium intake), and related aggregates (eg, diet quality), are both evaluated. This method has six analytical steps. (1) We included 560 risk-outcome pairs that met criteria for convincing or probable evidence on the basis of research studies. 12 risk-outcome pairs included in GBD 2017 no longer met inclusion criteria and 47 risk-outcome pairs for risks already included in GBD 2017 were added based on new evidence. (2) Relative risks were estimated as a function of exposure based on published systematic reviews, 81 systematic reviews done for GBD 2019, and meta-regression. (3) Levels of exposure in each age-sex-location-year included in the study were estimated based on all available data sources using spatiotemporal Gaussian process regression, DisMod-MR 2.1, a Bayesian meta-regression method, or alternative methods. (4) We determined, from published trials or cohort studies, the level of exposure associated with minimum risk, called the theoretical minimum risk exposure level. (5) Attributable deaths, YLLs, YLDs, and DALYs were computed by multiplying population attributable fractions (PAFs) by the relevant outcome quantity for each age-sex-location-year. (6) PAFs and attributable burden for combinations of risk factors were estimated taking into account mediation of different risk factors through other risk factors. Across all six analytical steps, 30 652 distinct data sources were used in the analysis. Uncertainty in each step of the analysis was propagated into the final estimates of attributable burden. Exposure levels for dichotomous, polytomous, and continuous risk factors were summarised with use of the summary exposure value to facilitate comparisons over time, across location, and across risks. Because the entire time series from 1990 to 2019 has been re-estimated with use of consistent data and methods, these results supersede previously published GBD estimates of attributable burden.

FINDINGS

The largest declines in risk exposure from 2010 to 2019 were among a set of risks that are strongly linked to social and economic development, including household air pollution; unsafe water, sanitation, and handwashing; and child growth failure. Global declines also occurred for tobacco smoking and lead exposure. The largest increases in risk exposure were for ambient particulate matter pollution, drug use, high fasting plasma glucose, and high body-mass index. In 2019, the leading Level 2 risk factor globally for attributable deaths was high systolic blood pressure, which accounted for 10·8 million (95% uncertainty interval [UI] 9·51-12·1) deaths (19·2% [16·9-21·3] of all deaths in 2019), followed by tobacco (smoked, second-hand, and chewing), which accounted for 8·71 million (8·12-9·31) deaths (15·4% [14·6-16·2] of all deaths in 2019). The leading Level 2 risk factor for attributable DALYs globally in 2019 was child and maternal malnutrition, which largely affects health in the youngest age groups and accounted for 295 million (253-350) DALYs (11·6% [10·3-13·1] of all global DALYs that year). The risk factor burden varied considerably in 2019 between age groups and locations. Among children aged 0-9 years, the three leading detailed risk factors for attributable DALYs were all related to malnutrition. Iron deficiency was the leading risk factor for those aged 10-24 years, alcohol use for those aged 25-49 years, and high systolic blood pressure for those aged 50-74 years and 75 years and older.

INTERPRETATION

Overall, the record for reducing exposure to harmful risks over the past three decades is poor. Success with reducing smoking and lead exposure through regulatory policy might point the way for a stronger role for public policy on other risks in addition to continued efforts to provide information on risk factor harm to the general public.

FUNDING

Bill & Melinda Gates Foundation.

Mapping trends in insecticide resistance phenotypes in African malaria vectors

Mitigating the threat of insecticide resistance in African malaria vector populations requires comprehensive information about where resistance occurs, to what degree, and how this has changed over time. Estimating these trends is complicated by the sparse, heterogeneous distribution of observations of resistance phenotypes in field populations. We use 6,423 observations of the prevalence of resistance to the most important vector control insecticides to inform a Bayesian geostatistical ensemble modelling approach, generating fine-scale predictive maps of resistance phenotypes in mosquitoes from the Anopheles gambiae complex across Africa. Our models are informed by a suite of 111 predictor variables describing potential drivers of selection for resistance. Our maps show alarming increases in the prevalence of resistance to pyrethroids and DDT across sub-Saharan Africa from 2005 to 2017, with mean mortality following insecticide exposure declining from almost 100% to less than 30% in some areas, as well as substantial spatial variation in resistance trends.

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 123 million neonatal, infant and child deaths between 2000 and 2017

Since 2000, many countries have achieved considerable success in improving child survival, but localized progress remains unclear. To inform efforts towards United Nations Sustainable Development Goal 3.2-to end preventable child deaths by 2030-we need consistently estimated data at the subnational level regarding child mortality rates and trends. Here we quantified, for the period 2000-2017, the subnational variation in mortality rates and number of deaths of neonates, infants and children under 5 years of age within 99 low- and middle-income countries using a geostatistical survival model. We estimated that 32% of children under 5 in these countries lived in districts that had attained rates of 25 or fewer child deaths per 1,000 live births by 2017, and that 58% of child deaths between 2000 and 2017 in these countries could have been averted in the absence of geographical inequality. This study enables the identification of high-mortality clusters, patterns of progress and geographical inequalities to inform appropriate investments and implementations that will help to improve the health of all populations.

Identifying residual hotspots and mapping lower respiratory infection morbidity and mortality in African children from 2000 to 2017

Lower respiratory infections (LRIs) are the leading cause of death in children under the age of 5, despite the existence of vaccines against many of their aetiologies. Furthermore, more than half of these deaths occur in Africa. Geospatial models can provide highly detailed estimates of trends subnationally, at the level where implementation of health policies has the greatest impact. We used Bayesian geostatistical modelling to estimate LRI incidence, prevalence and mortality in children under 5 subnationally in Africa for 2000–2017, using surveys covering 1.46 million children and 9,215,000 cases of LRI. Our model reveals large within-country variation in both health burden and its change over time. While reductions in childhood morbidity and mortality due to LRI were estimated for almost every country, we expose a cluster of residual high risk across seven countries, which averages 5.5 LRI deaths per 1,000 children per year. The preventable nature of the vast majority of LRI deaths mandates focused health system efforts in specific locations with the highest burden.

Geospatial modelling shows an overall decline in morbidity and mortality due to lower respiratory infections in Africa from 2000 to 2017, but also identifies subnational areas with residual high risk.

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.

Pareto rules for malaria super-spreaders and super-spreading

Heterogeneity in transmission is a challenge for infectious disease dynamics and control. An 80-20 "Pareto" rule has been proposed to describe this heterogeneity whereby 80% of transmission is accounted for by 20% of individuals, herein called super-spreaders. It is unclear, however, whether super-spreading can be attributed to certain individuals or whether it is an unpredictable and unavoidable feature of epidemics. Here, we investigate heterogeneous malaria transmission at three sites in Uganda and find that super-spreading is negatively correlated with overall malaria transmission intensity. Mosquito biting among humans is 90-10 at the lowest transmission intensities declining to less than 70-30 at the highest intensities. For super-spreaders, biting ranges from 70-30 down to 60-40. The difference, approximately half the total variance, is due to environmental stochasticity. Super-spreading is thus partly due to super-spreaders, but modest gains are expected from targeting super-spreaders.

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 origins and relatedness structure of mixed infections vary with local prevalence of P. falciparum malaria

Individual malaria infections can carry multiple strains of Plasmodium falciparum with varying levels of relatedness. Yet, how local epidemiology affects the properties of such mixed infections remains unclear. Here, we develop an enhanced method for strain deconvolution from genome sequencing data, which estimates the number of strains, their proportions, identity-by-descent (IBD) profiles and individual haplotypes. Applying it to the Pf3k data set, we find that the rate of mixed infection varies from 29% to 63% across countries and that 51% of mixed infections involve more than two strains. Furthermore, we estimate that 47% of symptomatic dual infections contain sibling strains likely to have been co-transmitted from a single mosquito, and find evidence of mixed infections propagated over successive infection cycles. Finally, leveraging data from the Malaria Atlas Project, we find that prevalence correlates within Africa, but not Asia, with both the rate of mixed infection and the level of IBD.

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.

Mapping diphtheria-pertussis-tetanus vaccine coverage in Africa, 2000-2016: a spatial and temporal modelling study

BACKGROUND

Routine childhood vaccination is among the most cost-effective, successful public health interventions available. Amid substantial investments to expand vaccine delivery throughout Africa and strengthen administrative reporting systems, most countries still require robust measures of local routine vaccine coverage and changes in geographical inequalities over time.

METHODS

This analysis drew from 183 surveys done between 2000 and 2016, including data from 881 268 children in 49 African countries. We used a Bayesian geostatistical model calibrated to results from the Global Burden of Diseases, Injuries, and Risk Factors Study 2017, to produce annual estimates with high-spatial resolution (5 ×    5 km) of diphtheria-pertussis-tetanus (DPT) vaccine coverage and dropout for children aged 12-23 months in 52 African countries from 2000 to 2016.

FINDINGS

Estimated third-dose (DPT3) coverage increased in 72·3% (95% uncertainty interval [UI] 64·6-80·3) of second-level administrative units in Africa from 2000 to 2016, but substantial geographical inequalities in DPT coverage remained across and within African countries. In 2016, DPT3 coverage at the second administrative (ie, district) level varied by more than 25% in 29 of 52 countries, with only two (Morocco and Rwanda) of 52 countries meeting the Global Vaccine Action Plan target of 80% DPT3 coverage or higher in all second-level administrative units with high confidence (posterior probability ≥95%). Large areas of low DPT3 coverage (≤50%) were identified in the Sahel, Somalia, eastern Ethiopia, and in Angola. Low first-dose (DPT1) coverage (≤50%) and high relative dropout (≥30%) together drove low DPT3 coverage across the Sahel, Somalia, eastern Ethiopia, Guinea, and Angola.

INTERPRETATION

Despite substantial progress in Africa, marked national and subnational inequalities in DPT coverage persist throughout the continent. These results can help identify areas of low coverage and vaccine delivery system vulnerabilities and can ultimately support more precise targeting of resources to improve vaccine coverage and health outcomes for African children.

FUNDING

Bill & Melinda Gates Foundation.

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.

A stakeholder workshop about modelled maps of key malaria indicator survey indicators in Madagascar

The Demographic and Health Surveys (DHS) Program has supported three household Malaria Indicator Surveys (MIS) in Madagascar. The results of 13 key malaria indicators from these surveys have been mapped as continuous surfaces using model-based geostatistical methods. The opportunities and limitations of these mapped outputs were discussed during a workshop in Antananarivo, Madagascar in July 2018, attended by 15 representatives from various implementation, policy and research stakeholder institutions in Madagascar. Participants evaluated the findings from the maps, using these to develop figures and narratives to support their work in the control of malaria in Madagascar.

Growing evidence of Plasmodium vivax across malaria-endemic Africa

Effective malaria control strategies require an accurate understanding of the epidemiology of locally transmitted Plasmodium species. Compared to Plasmodium falciparum infection, Plasmodium vivax has a lower asexual parasitaemia, forms dormant liver-stages (hypnozoites), and is more transmissible. Hence, treatment and diagnostic policies aimed exclusively at P. falciparum are far less efficient against endemic P. vivax. Within sub-Saharan Africa, malaria control programmes justly focus on reducing the morbidity and mortality associated with P. falciparum. However, the recent emphasis on malaria elimination and increased accessibility of more sensitive diagnostic tools have revealed greater intricacies in malaria epidemiology across the continent. Since 2010, the number of studies identifying P. vivax endemic to Africa has expanded considerably, with 88 new scientific reports published since a review of evidence in 2015, approximately doubling the available data. There is evidence of P. vivax in all regions of Africa, apparent from infected vectors, clinical cases, serological indicators, parasite prevalence, exported infections, and P. vivax-infected Duffy-negative individuals. Where the prevalence of microscopic parasitaemia is low, a greater proportion of P. vivax infections were observed relative to P. falciparum. This evidence highlights an underlying widespread presence of P. vivax across all malaria-endemic regions of Africa, further complicating the current practical understanding of malaria epidemiology in this region. Thus, ultimate elimination of malaria in Africa will require national malaria control programmes to adopt policy and practice aimed at all human species of malaria.

Global, regional, and national age-sex-specific mortality and life expectancy, 1950-2017: a systematic analysis for the Global Burden of Disease Study 2017

BACKGROUND

Assessments of age-specific mortality and life expectancy have been done by the UN Population Division, Department of Economics and Social Affairs (UNPOP), the United States Census Bureau, WHO, and as part of previous iterations of the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD). Previous iterations of the GBD used population estimates from UNPOP, which were not derived in a way that was internally consistent with the estimates of the numbers of deaths in the GBD. The present iteration of the GBD, GBD 2017, improves on previous assessments and provides timely estimates of the mortality experience of populations globally.

METHODS

The GBD uses all available data to produce estimates of mortality rates between 1950 and 2017 for 23 age groups, both sexes, and 918 locations, including 195 countries and territories and subnational locations for 16 countries. Data used include vital registration systems, sample registration systems, household surveys (complete birth histories, summary birth histories, sibling histories), censuses (summary birth histories, household deaths), and Demographic Surveillance Sites. In total, this analysis used 8259 data sources. Estimates of the probability of death between birth and the age of 5 years and between ages 15 and 60 years are generated and then input into a model life table system to produce complete life tables for all locations and years. Fatal discontinuities and mortality due to HIV/AIDS are analysed separately and then incorporated into the estimation. We analyse the relationship between age-specific mortality and development status using the Socio-demographic Index, a composite measure based on fertility under the age of 25 years, education, and income. There are four main methodological improvements in GBD 2017 compared with GBD 2016: 622 additional data sources have been incorporated; new estimates of population, generated by the GBD study, are used; statistical methods used in different components of the analysis have been further standardised and improved; and the analysis has been extended backwards in time by two decades to start in 1950.

FINDINGS

Globally, 18·7% (95% uncertainty interval 18·4-19·0) of deaths were registered in 1950 and that proportion has been steadily increasing since, with 58·8% (58·2-59·3) of all deaths being registered in 2015. At the global level, between 1950 and 2017, life expectancy increased from 48·1 years (46·5-49·6) to 70·5 years (70·1-70·8) for men and from 52·9 years (51·7-54·0) to 75·6 years (75·3-75·9) for women. Despite this overall progress, there remains substantial variation in life expectancy at birth in 2017, which ranges from 49·1 years (46·5-51·7) for men in the Central African Republic to 87·6 years (86·9-88·1) among women in Singapore. The greatest progress across age groups was for children younger than 5 years; under-5 mortality dropped from 216·0 deaths (196·3-238·1) per 1000 livebirths in 1950 to 38·9 deaths (35·6-42·83) per 1000 livebirths in 2017, with huge reductions across countries. Nevertheless, there were still 5·4 million (5·2-5·6) deaths among children younger than 5 years in the world in 2017. Progress has been less pronounced and more variable for adults, especially for adult males, who had stagnant or increasing mortality rates in several countries. The gap between male and female life expectancy between 1950 and 2017, while relatively stable at the global level, shows distinctive patterns across super-regions and has consistently been the largest in central Europe, eastern Europe, and central Asia, and smallest in south Asia. Performance was also variable across countries and time in observed mortality rates compared with those expected on the basis of development.

INTERPRETATION

This analysis of age-sex-specific mortality shows that there are remarkably complex patterns in population mortality across countries. The findings of this study highlight global successes, such as the large decline in under-5 mortality, which reflects significant local, national, and global commitment and investment over several decades. However, they also bring attention to mortality patterns that are a cause for concern, particularly among adult men and, to a lesser extent, women, whose mortality rates have stagnated in many countries over the time period of this study, and in some cases are increasing.

FUNDING

Bill & Melinda Gates Foundation.

Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980-2017: a systematic analysis for the Global Burden of Disease Study 2017

BACKGROUND

Global development goals increasingly rely on country-specific estimates for benchmarking a nation's progress. To meet this need, the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2016 estimated global, regional, national, and, for selected locations, subnational cause-specific mortality beginning in the year 1980. Here we report an update to that study, making use of newly available data and improved methods. GBD 2017 provides a comprehensive assessment of cause-specific mortality for 282 causes in 195 countries and territories from 1980 to 2017.

METHODS

The causes of death database is composed of vital registration (VR), verbal autopsy (VA), registry, survey, police, and surveillance data. GBD 2017 added ten VA studies, 127 country-years of VR data, 502 cancer-registry country-years, and an additional surveillance country-year. Expansions of the GBD cause of death hierarchy resulted in 18 additional causes estimated for GBD 2017. Newly available data led to subnational estimates for five additional countries-Ethiopia, Iran, New Zealand, Norway, and Russia. Deaths assigned International Classification of Diseases (ICD) codes for non-specific, implausible, or intermediate causes of death were reassigned to underlying causes by redistribution algorithms that were incorporated into uncertainty estimation. We used statistical modelling tools developed for GBD, including the Cause of Death Ensemble model (CODEm), to generate cause fractions and cause-specific death rates for each location, year, age, and sex. Instead of using UN estimates as in previous versions, GBD 2017 independently estimated population size and fertility rate for all locations. Years of life lost (YLLs) were then calculated as the sum of each death multiplied by the standard life expectancy at each age. All rates reported here are age-standardised.

FINDINGS

At the broadest grouping of causes of death (Level 1), non-communicable diseases (NCDs) comprised the greatest fraction of deaths, contributing to 73·4% (95% uncertainty interval [UI] 72·5-74·1) of total deaths in 2017, while communicable, maternal, neonatal, and nutritional (CMNN) causes accounted for 18·6% (17·9-19·6), and injuries 8·0% (7·7-8·2). Total numbers of deaths from NCD causes increased from 2007 to 2017 by 22·7% (21·5-23·9), representing an additional 7·61 million (7·20-8·01) deaths estimated in 2017 versus 2007. The death rate from NCDs decreased globally by 7·9% (7·0-8·8). The number of deaths for CMNN causes decreased by 22·2% (20·0-24·0) and the death rate by 31·8% (30·1-33·3). Total deaths from injuries increased by 2·3% (0·5-4·0) between 2007 and 2017, and the death rate from injuries decreased by 13·7% (12·2-15·1) to 57·9 deaths (55·9-59·2) per 100 000 in 2017. Deaths from substance use disorders also increased, rising from 284 000 deaths (268 000-289 000) globally in 2007 to 352 000 (334 000-363 000) in 2017. Between 2007 and 2017, total deaths from conflict and terrorism increased by 118·0% (88·8-148·6). A greater reduction in total deaths and death rates was observed for some CMNN causes among children younger than 5 years than for older adults, such as a 36·4% (32·2-40·6) reduction in deaths from lower respiratory infections for children younger than 5 years compared with a 33·6% (31·2-36·1) increase in adults older than 70 years. Globally, the number of deaths was greater for men than for women at most ages in 2017, except at ages older than 85 years. Trends in global YLLs reflect an epidemiological transition, with decreases in total YLLs from enteric infections, respiratory infections and tuberculosis, and maternal and neonatal disorders between 1990 and 2017; these were generally greater in magnitude at the lowest levels of the Socio-demographic Index (SDI). At the same time, there were large increases in YLLs from neoplasms and cardiovascular diseases. YLL rates decreased across the five leading Level 2 causes in all SDI quintiles. The leading causes of YLLs in 1990-neonatal disorders, lower respiratory infections, and diarrhoeal diseases-were ranked second, fourth, and fifth, in 2017. Meanwhile, estimated YLLs increased for ischaemic heart disease (ranked first in 2017) and stroke (ranked third), even though YLL rates decreased. Population growth contributed to increased total deaths across the 20 leading Level 2 causes of mortality between 2007 and 2017. Decreases in the cause-specific mortality rate reduced the effect of population growth for all but three causes: substance use disorders, neurological disorders, and skin and subcutaneous diseases.

INTERPRETATION

Improvements in global health have been unevenly distributed among populations. Deaths due to injuries, substance use disorders, armed conflict and terrorism, neoplasms, and cardiovascular disease are expanding threats to global health. For causes of death such as lower respiratory and enteric infections, more rapid progress occurred for children than for the oldest adults, and there is continuing disparity in mortality rates by sex across age groups. Reductions in the death rate of some common diseases are themselves slowing or have ceased, primarily for NCDs, and the death rate for selected causes has increased in the past decade.

FUNDING

Bill & Melinda Gates Foundation.

Measuring progress from 1990 to 2017 and projecting attainment to 2030 of the health-related Sustainable Development Goals for 195 countries and territories: a systematic analysis for the Global Burden of Disease Study 2017

BACKGROUND

Efforts to establish the 2015 baseline and monitor early implementation of the UN Sustainable Development Goals (SDGs) highlight both great potential for and threats to improving health by 2030. To fully deliver on the SDG aim of "leaving no one behind", it is increasingly important to examine the health-related SDGs beyond national-level estimates. As part of the Global Burden of Diseases, Injuries, and Risk Factors Study 2017 (GBD 2017), we measured progress on 41 of 52 health-related SDG indicators and estimated the health-related SDG index for 195 countries and territories for the period 1990-2017, projected indicators to 2030, and analysed global attainment.

METHODS

We measured progress on 41 health-related SDG indicators from 1990 to 2017, an increase of four indicators since GBD 2016 (new indicators were health worker density, sexual violence by non-intimate partners, population census status, and prevalence of physical and sexual violence [reported separately]). We also improved the measurement of several previously reported indicators. We constructed national-level estimates and, for a subset of health-related SDGs, examined indicator-level differences by sex and Socio-demographic Index (SDI) quintile. We also did subnational assessments of performance for selected countries. To construct the health-related SDG index, we transformed the value for each indicator on a scale of 0-100, with 0 as the 2·5th percentile and 100 as the 97·5th percentile of 1000 draws calculated from 1990 to 2030, and took the geometric mean of the scaled indicators by target. To generate projections through 2030, we used a forecasting framework that drew estimates from the broader GBD study and used weighted averages of indicator-specific and country-specific annualised rates of change from 1990 to 2017 to inform future estimates. We assessed attainment of indicators with defined targets in two ways: first, using mean values projected for 2030, and then using the probability of attainment in 2030 calculated from 1000 draws. We also did a global attainment analysis of the feasibility of attaining SDG targets on the basis of past trends. Using 2015 global averages of indicators with defined SDG targets, we calculated the global annualised rates of change required from 2015 to 2030 to meet these targets, and then identified in what percentiles the required global annualised rates of change fell in the distribution of country-level rates of change from 1990 to 2015. We took the mean of these global percentile values across indicators and applied the past rate of change at this mean global percentile to all health-related SDG indicators, irrespective of target definition, to estimate the equivalent 2030 global average value and percentage change from 2015 to 2030 for each indicator.

FINDINGS

The global median health-related SDG index in 2017 was 59·4 (IQR 35·4-67·3), ranging from a low of 11·6 (95% uncertainty interval 9·6-14·0) to a high of 84·9 (83·1-86·7). SDG index values in countries assessed at the subnational level varied substantially, particularly in China and India, although scores in Japan and the UK were more homogeneous. Indicators also varied by SDI quintile and sex, with males having worse outcomes than females for non-communicable disease (NCD) mortality, alcohol use, and smoking, among others. Most countries were projected to have a higher health-related SDG index in 2030 than in 2017, while country-level probabilities of attainment by 2030 varied widely by indicator. Under-5 mortality, neonatal mortality, maternal mortality ratio, and malaria indicators had the most countries with at least 95% probability of target attainment. Other indicators, including NCD mortality and suicide mortality, had no countries projected to meet corresponding SDG targets on the basis of projected mean values for 2030 but showed some probability of attainment by 2030. For some indicators, including child malnutrition, several infectious diseases, and most violence measures, the annualised rates of change required to meet SDG targets far exceeded the pace of progress achieved by any country in the recent past. We found that applying the mean global annualised rate of change to indicators without defined targets would equate to about 19% and 22% reductions in global smoking and alcohol consumption, respectively; a 47% decline in adolescent birth rates; and a more than 85% increase in health worker density per 1000 population by 2030.

INTERPRETATION

The GBD study offers a unique, robust platform for monitoring the health-related SDGs across demographic and geographic dimensions. Our findings underscore the importance of increased collection and analysis of disaggregated data and highlight where more deliberate design or targeting of interventions could accelerate progress in attaining the SDGs. Current projections show that many health-related SDG indicators, NCDs, NCD-related risks, and violence-related indicators will require a concerted shift away from what might have driven past gains-curative interventions in the case of NCDs-towards multisectoral, prevention-oriented policy action and investments to achieve SDG aims. Notably, several targets, if they are to be met by 2030, demand a pace of progress that no country has achieved in the recent past. The future is fundamentally uncertain, and no model can fully predict what breakthroughs or events might alter the course of the SDGs. What is clear is that our actions-or inaction-today will ultimately dictate how close the world, collectively, can get to leaving no one behind by 2030.

FUNDING

Bill & Melinda Gates Foundation.

Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks for 195 countries and territories, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017

BACKGROUND

The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2017 comparative risk assessment (CRA) is a comprehensive approach to risk factor quantification that offers a useful tool for synthesising evidence on risks and risk-outcome associations. With each annual GBD study, we update the GBD CRA to incorporate improved methods, new risks and risk-outcome pairs, and new data on risk exposure levels and risk-outcome associations.

METHODS

We used the CRA framework developed for previous iterations of GBD to estimate levels and trends in exposure, attributable deaths, and attributable disability-adjusted life-years (DALYs), by age group, sex, year, and location for 84 behavioural, environmental and occupational, and metabolic risks or groups of risks from 1990 to 2017. This study included 476 risk-outcome pairs that met the GBD study criteria for convincing or probable evidence of causation. We extracted relative risk and exposure estimates from 46 749 randomised controlled trials, cohort studies, household surveys, census data, satellite data, and other sources. We used statistical models to pool data, adjust for bias, and incorporate covariates. Using the counterfactual scenario of theoretical minimum risk exposure level (TMREL), we estimated the portion of deaths and DALYs that could be attributed to a given risk. We explored the relationship between development and risk exposure by modelling the relationship between the Socio-demographic Index (SDI) and risk-weighted exposure prevalence and estimated expected levels of exposure and risk-attributable burden by SDI. Finally, we explored temporal changes in risk-attributable DALYs by decomposing those changes into six main component drivers of change as follows: (1) population growth; (2) changes in population age structures; (3) changes in exposure to environmental and occupational risks; (4) changes in exposure to behavioural risks; (5) changes in exposure to metabolic risks; and (6) changes due to all other factors, approximated as the risk-deleted death and DALY rates, where the risk-deleted rate is the rate that would be observed had we reduced the exposure levels to the TMREL for all risk factors included in GBD 2017.

FINDINGS

In 2017, 34·1 million (95% uncertainty interval [UI] 33·3-35·0) deaths and 1·21 billion (1·14-1·28) DALYs were attributable to GBD risk factors. Globally, 61·0% (59·6-62·4) of deaths and 48·3% (46·3-50·2) of DALYs were attributed to the GBD 2017 risk factors. When ranked by risk-attributable DALYs, high systolic blood pressure (SBP) was the leading risk factor, accounting for 10·4 million (9·39-11·5) deaths and 218 million (198-237) DALYs, followed by smoking (7·10 million [6·83-7·37] deaths and 182 million [173-193] DALYs), high fasting plasma glucose (6·53 million [5·23-8·23] deaths and 171 million [144-201] DALYs), high body-mass index (BMI; 4·72 million [2·99-6·70] deaths and 148 million [98·6-202] DALYs), and short gestation for birthweight (1·43 million [1·36-1·51] deaths and 139 million [131-147] DALYs). In total, risk-attributable DALYs declined by 4·9% (3·3-6·5) between 2007 and 2017. In the absence of demographic changes (ie, population growth and ageing), changes in risk exposure and risk-deleted DALYs would have led to a 23·5% decline in DALYs during that period. Conversely, in the absence of changes in risk exposure and risk-deleted DALYs, demographic changes would have led to an 18·6% increase in DALYs during that period. The ratios of observed risk exposure levels to exposure levels expected based on SDI (O/E ratios) increased globally for unsafe drinking water and household air pollution between 1990 and 2017. This result suggests that development is occurring more rapidly than are changes in the underlying risk structure in a population. Conversely, nearly universal declines in O/E ratios for smoking and alcohol use indicate that, for a given SDI, exposure to these risks is declining. In 2017, the leading Level 4 risk factor for age-standardised DALY rates was high SBP in four super-regions: central Europe, eastern Europe, and central Asia; north Africa and Middle East; south Asia; and southeast Asia, east Asia, and Oceania. The leading risk factor in the high-income super-region was smoking, in Latin America and Caribbean was high BMI, and in sub-Saharan Africa was unsafe sex. O/E ratios for unsafe sex in sub-Saharan Africa were notably high, and those for alcohol use in north Africa and the Middle East were notably low.

INTERPRETATION

By quantifying levels and trends in exposures to risk factors and the resulting disease burden, this assessment offers insight into where past policy and programme efforts might have been successful and highlights current priorities for public health action. Decreases in behavioural, environmental, and occupational risks have largely offset the effects of population growth and ageing, in relation to trends in absolute burden. Conversely, the combination of increasing metabolic risks and population ageing will probably continue to drive the increasing trends in non-communicable diseases at the global level, which presents both a public health challenge and opportunity. We see considerable spatiotemporal heterogeneity in levels of risk exposure and risk-attributable burden. Although levels of development underlie some of this heterogeneity, O/E ratios show risks for which countries are overperforming or underperforming relative to their level of development. As such, these ratios provide a benchmarking tool to help to focus local decision making. Our findings reinforce the importance of both risk exposure monitoring and epidemiological research to assess causal connections between risks and health outcomes, and they highlight the usefulness of the GBD study in synthesising data to draw comprehensive and robust conclusions that help to inform good policy and strategic health planning.

FUNDING

Bill & Melinda Gates Foundation.

Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017

BACKGROUND

The Global Burden of Diseases, Injuries, and Risk Factors Study 2017 (GBD 2017) includes a comprehensive assessment of incidence, prevalence, and years lived with disability (YLDs) for 354 causes in 195 countries and territories from 1990 to 2017. Previous GBD studies have shown how the decline of mortality rates from 1990 to 2016 has led to an increase in life expectancy, an ageing global population, and an expansion of the non-fatal burden of disease and injury. These studies have also shown how a substantial portion of the world's population experiences non-fatal health loss with considerable heterogeneity among different causes, locations, ages, and sexes. Ongoing objectives of the GBD study include increasing the level of estimation detail, improving analytical strategies, and increasing the amount of high-quality data.

METHODS

We estimated incidence and prevalence for 354 diseases and injuries and 3484 sequelae. We used an updated and extensive body of literature studies, survey data, surveillance data, inpatient admission records, outpatient visit records, and health insurance claims, and additionally used results from cause of death models to inform estimates using a total of 68 781 data sources. Newly available clinical data from India, Iran, Japan, Jordan, Nepal, China, Brazil, Norway, and Italy were incorporated, as well as updated claims data from the USA and new claims data from Taiwan (province of China) and Singapore. We used DisMod-MR 2.1, a Bayesian meta-regression tool, as the main method of estimation, ensuring consistency between rates of incidence, prevalence, remission, and cause of death for each condition. YLDs were estimated as the product of a prevalence estimate and a disability weight for health states of each mutually exclusive sequela, adjusted for comorbidity. We updated the Socio-demographic Index (SDI), a summary development indicator of income per capita, years of schooling, and total fertility rate. Additionally, we calculated differences between male and female YLDs to identify divergent trends across sexes. GBD 2017 complies with the Guidelines for Accurate and Transparent Health Estimates Reporting.

FINDINGS

Globally, for females, the causes with the greatest age-standardised prevalence were oral disorders, headache disorders, and haemoglobinopathies and haemolytic anaemias in both 1990 and 2017. For males, the causes with the greatest age-standardised prevalence were oral disorders, headache disorders, and tuberculosis including latent tuberculosis infection in both 1990 and 2017. In terms of YLDs, low back pain, headache disorders, and dietary iron deficiency were the leading Level 3 causes of YLD counts in 1990, whereas low back pain, headache disorders, and depressive disorders were the leading causes in 2017 for both sexes combined. All-cause age-standardised YLD rates decreased by 3·9% (95% uncertainty interval [UI] 3·1-4·6) from 1990 to 2017; however, the all-age YLD rate increased by 7·2% (6·0-8·4) while the total sum of global YLDs increased from 562 million (421-723) to 853 million (642-1100). The increases for males and females were similar, with increases in all-age YLD rates of 7·9% (6·6-9·2) for males and 6·5% (5·4-7·7) for females. We found significant differences between males and females in terms of age-standardised prevalence estimates for multiple causes. The causes with the greatest relative differences between sexes in 2017 included substance use disorders (3018 cases [95% UI 2782-3252] per 100 000 in males vs s1400 [1279-1524] per 100 000 in females), transport injuries (3322 [3082-3583] vs 2336 [2154-2535]), and self-harm and interpersonal violence (3265 [2943-3630] vs 5643 [5057-6302]).

INTERPRETATION

Global all-cause age-standardised YLD rates have improved only slightly over a period spanning nearly three decades. However, the magnitude of the non-fatal disease burden has expanded globally, with increasing numbers of people who have a wide spectrum of conditions. A subset of conditions has remained globally pervasive since 1990, whereas other conditions have displayed more dynamic trends, with different ages, sexes, and geographies across the globe experiencing varying burdens and trends of health loss. This study emphasises how global improvements in premature mortality for select conditions have led to older populations with complex and potentially expensive diseases, yet also highlights global achievements in certain domains of disease and injury.

FUNDING

Bill & Melinda Gates Foundation.

Estimates of the global, regional, and national morbidity, mortality, and aetiologies of lower respiratory infections in 195 countries, 1990-2016: a systematic analysis for the Global Burden of Disease Study 2016

BACKGROUND

Lower respiratory infections are a leading cause of morbidity and mortality around the world. The Global Burden of Diseases, Injuries, and Risk Factors (GBD) Study 2016, provides an up-to-date analysis of the burden of lower respiratory infections in 195 countries. This study assesses cases, deaths, and aetiologies spanning the past 26 years and shows how the burden of lower respiratory infection has changed in people of all ages.

METHODS

We used three separate modelling strategies for lower respiratory infections in GBD 2016: a Bayesian hierarchical ensemble modelling platform (Cause of Death Ensemble model), which uses vital registration, verbal autopsy data, and surveillance system data to predict mortality due to lower respiratory infections; a compartmental meta-regression tool (DisMod-MR), which uses scientific literature, population representative surveys, and health-care data to predict incidence, prevalence, and mortality; and modelling of counterfactual estimates of the population attributable fraction of lower respiratory infection episodes due to Streptococcus pneumoniae, Haemophilus influenzae type b, influenza, and respiratory syncytial virus. We calculated each modelled estimate for each age, sex, year, and location. We modelled the exposure level in a population for a given risk factor using DisMod-MR and a spatio-temporal Gaussian process regression, and assessed the effectiveness of targeted interventions for each risk factor in children younger than 5 years. We also did a decomposition analysis of the change in LRI deaths from 2000-16 using the risk factors associated with LRI in GBD 2016.

FINDINGS

In 2016, lower respiratory infections caused 652 572 deaths (95% uncertainty interval [UI] 586 475-720 612) in children younger than 5 years (under-5s), 1 080 958 deaths (943 749-1 170 638) in adults older than 70 years, and 2 377 697 deaths (2 145 584-2 512 809) in people of all ages, worldwide. Streptococcus pneumoniae was the leading cause of lower respiratory infection morbidity and mortality globally, contributing to more deaths than all other aetiologies combined in 2016 (1 189 937 deaths, 95% UI 690 445-1 770 660). Childhood wasting remains the leading risk factor for lower respiratory infection mortality among children younger than 5 years, responsible for 61·4% of lower respiratory infection deaths in 2016 (95% UI 45·7-69·6). Interventions to improve wasting, household air pollution, ambient particulate matter pollution, and expanded antibiotic use could avert one under-5 death due to lower respiratory infection for every 4000 children treated in the countries with the highest lower respiratory infection burden.

INTERPRETATION

Our findings show substantial progress in the reduction of lower respiratory infection burden, but this progress has not been equal across locations, has been driven by decreases in several primary risk factors, and might require more effort among elderly adults. By highlighting regions and populations with the highest burden, and the risk factors that could have the greatest effect, funders, policy makers, and programme implementers can more effectively reduce lower respiratory infections among the world's most susceptible populations.

FUNDING

Bill & Melinda Gates Foundation.

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.

Variation in Childhood Diarrheal Morbidity and Mortality in Africa, 2000-2015

BACKGROUND

Diarrheal diseases are the third leading cause of disease and death in children younger than 5 years of age in Africa and were responsible for an estimated 30 million cases of severe diarrhea (95% credible interval, 27 million to 33 million) and 330,000 deaths (95% credible interval, 270,000 to 380,000) in 2015. The development of targeted approaches to address this burden has been hampered by a paucity of comprehensive, fine-scale estimates of diarrhea-related disease and death among and within countries.

METHODS

We produced annual estimates of the prevalence and incidence of diarrhea and diarrhea-related mortality with high geographic detail (5 km2) across Africa from 2000 through 2015. Estimates were created with the use of Bayesian geostatistical techniques and were calibrated to the results from the Global Burden of Diseases, Injuries, and Risk Factors Study 2016.

RESULTS

The results revealed geographic inequality with regard to diarrhea risk in Africa. Of the estimated 330,000 childhood deaths that were attributable to diarrhea in 2015, more than 50% occurred in 55 of the 782 first-level administrative subdivisions (e.g., states). In 2015, mortality rates among first-level administrative subdivisions in Nigeria differed by up to a factor of 6. The case fatality rates were highly varied at the national level across Africa, with the highest values observed in Benin, Lesotho, Mali, Nigeria, and Sierra Leone.

CONCLUSIONS

Our findings showed concentrated areas of diarrheal disease and diarrhea-related death in countries that had a consistently high burden as well as in countries that had considerable national-level reductions in diarrhea burden. (Funded by the Bill and Melinda Gates Foundation.).

Mapping local variation in educational attainment across Africa

Educational attainment for women of reproductive age is linked to reduced child and maternal mortality, lower fertility and improved reproductive health. Comparable analyses of attainment exist only at the national level, potentially obscuring patterns in subnational inequality. Evidence suggests that wide disparities between urban and rural populations exist, raising questions about where the majority of progress towards the education targets of the Sustainable Development Goals is occurring in African countries. Here we explore within-country inequalities by predicting years of schooling across five by five kilometre grids, generating estimates of average educational attainment by age and sex at subnational levels. Despite marked progress in attainment from 2000 to 2015 across Africa, substantial differences persist between locations and sexes. These differences have widened in many countries, particularly across the Sahel. These high-resolution, comparable estimates improve the ability of decision-makers to plan the precisely targeted interventions that will be necessary to deliver progress during the era of the Sustainable Development Goals.

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.

How long do rapid diagnostic tests remain positive after anti-malarial treatment?

BACKGROUND

Rapid diagnostic tests (RDTs) are increasingly becoming a paradigm for both clinical diagnosis of malaria infections and for estimating community parasite prevalence in household malaria indicator surveys in malaria-endemic countries. The antigens detected by RDTs are known to persist in the blood after treatment with anti-malarials, but reports on the duration of persistence (and the effect this has on RDT positivity) of these antigens post-treatment have been variable.

METHODS

In this review, published studies on the persistence of positivity of RDTs post-treatment are collated, and a bespoke Bayesian survival model is fit to estimate the number of days RDTs remain positive after treatment.

RESULTS

Half of RDTs that detect the antigen histidine-rich protein II (HRP2) are still positive 15 (5-32) days post-treatment, 13 days longer than RDTs that detect the antigen Plasmodium lactate dehydrogenase, and that 5% of HRP2 RDTs are still positive 36 (21-61) days after treatment. The duration of persistent positivity for combination RDTs that detect both antigens falls between that for HRP2- or pLDH-only RDTs, with half of RDTs remaining positive at 7 (2-20) days post-treatment. This study shows that children display persistent RDT positivity for longer after treatment than adults, and that persistent positivity is more common when an individual is treated with artemisinin combination therapy than when treated with other anti-malarials.

CONCLUSIONS

RDTs remain positive for a highly variable amount of time after treatment with anti-malarials, and the duration of positivity is highly dependent on the type of RDT used for diagnosis. Additionally, age and treatment both impact the duration of persistence of RDT positivity. The results presented here suggest that caution should be taken when using RDT-derived diagnostic outcomes from cross-sectional data where individuals have had a recent history of anti-malarial treatment.

Associated patterns of insecticide resistance in field populations of malaria vectors across Africa

The development of insecticide resistance in African malaria vectors threatens the continued efficacy of important vector control methods that rely on a limited set of insecticides. To understand the operational significance of resistance we require quantitative information about levels of resistance in field populations to the suite of vector control insecticides. Estimation of resistance is complicated by the sparsity of observations in field populations, variation in resistance over time and space at local and regional scales, and cross-resistance between different insecticide types. Using observations of the prevalence of resistance in mosquito species from the Anopheles gambiae complex sampled from 1,183 locations throughout Africa, we applied Bayesian geostatistical models to quantify patterns of covariation in resistance phenotypes across different insecticides. For resistance to the three pyrethroids tested, deltamethrin, permethrin, and λ-cyhalothrin, we found consistent forms of covariation across sub-Saharan Africa and covariation between resistance to these pyrethroids and resistance to DDT. We found no evidence of resistance interactions between carbamate and organophosphate insecticides or between these insecticides and those from other classes. For pyrethroids and DDT we found significant associations between predicted mean resistance and the observed frequency of kdr mutations in the Vgsc gene in field mosquito samples, with DDT showing the strongest association. These results improve our capacity to understand and predict resistance patterns throughout Africa and can guide the development of monitoring strategies.

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.