Projected health impact of post-discharge malaria chemoprevention among children with severe malarial anaemia in Africa

Children recovering from severe malarial anaemia (SMA) remain at high risk of readmission and death after discharge from hospital. However, a recent trial found that post-discharge malaria chemoprevention (PDMC) with dihydroartemisinin-piperaquine reduces this risk. We developed a mathematical model describing the daily incidence of uncomplicated and severe malaria requiring readmission among 0-5-year old children after hospitalised SMA. We fitted the model to a multicentre clinical PDMC trial using Bayesian methods and modelled the potential impact of PDMC across malaria-endemic African countries. In the 20 highest-burden countries, we estimate that only 2-5 children need to be given PDMC to prevent one hospitalised malaria episode, and less than 100 to prevent one death. If all hospitalised SMA cases access PDMC in moderate-to-high transmission areas, 38,600 (range 16,900-88,400) malaria-associated readmissions could be prevented annually, depending on access to hospital care. We estimate that recurrent SMA post-discharge constitutes 19% of all SMA episodes in moderate-to-high transmission settings.

Joint spatiotemporal modelling reveals seasonally dynamic patterns of Japanese encephalitis vector abundance across India

Predicting vector abundance and seasonality, key components of mosquito-borne disease (MBD) hazard, is essential to determine hotspots of MBD risk and target interventions effectively. Japanese encephalitis (JE), an important MBD, is a leading cause of viral encephalopathy in Asia with 100,000 cases estimated annually, but data on the principal vector Culex tritaeniorhynchus is lacking. We developed a Bayesian joint-likelihood model that combined information from available vector occurrence and abundance data to predict seasonal vector abundance for C. tritaeniorhynchus (a constituent of JE hazard) across India, as well as examining the environmental drivers of these patterns. Using data collated from 57 locations from 24 studies, we find distinct seasonal and spatial patterns of JE vector abundance influenced by climatic and land use factors. Lagged precipitation, temperature and land use intensity metrics for rice crop cultivation were the main drivers of vector abundance, independent of seasonal, or spatial variation. The inclusion of environmental factors and a seasonal term improved model prediction accuracy (mean absolute error [MAE] for random cross validation = 0.48) compared to a baseline model representative of static hazard predictions (MAE = 0.95), signalling the importance of seasonal environmental conditions in predicting JE vector abundance. Vector abundance varied widely across India with high abundance predicted in northern, north-eastern, eastern, and southern regions, although this ranged from seasonal (e.g., Uttar Pradesh, West Bengal) to perennial (e.g., Assam, Tamil Nadu). One-month lagged predicted vector abundance was a significant predictor of JE outbreaks (odds ratio 2.45, 95% confidence interval: 1.52–4.08), highlighting the possible development of vector abundance as a proxy for JE hazard. We demonstrate a novel approach that leverages information from sparse vector surveillance data to predict seasonal vector abundance–a key component of JE hazard–over large spatial scales, providing decision-makers with better guidance for targeting vector surveillance and control efforts.

Japanese encephalitis (JE) is the leading cause of viral encephalopathy in Asia with an estimated 100,000 annual cases and 25,000 deaths. However, insufficient data on the predominant mosquito vector Culex tritaeniorhynchus–a key component of JE hazard–precludes hazard estimation required to target public health interventions. Previous studies have provided limited estimates of JE hazard, often predicting geographic distributions of potential vector occurrence without accounting for vector abundance, seasonality, or uncertainty in predictions. This study details a novel approach to predict spatiotemporal patterns in JE vector abundance using a joint-likelihood modelling technique that leverages information from sparse vector surveillance data. We showed that patterns in JE vector abundance were driven by seasonality and environmental factors and so demonstrated the limitations of previously available static vector distribution maps in estimating the vector population component of JE hazard. One-month lagged vector abundance predictions showed a positive relationship with JE outbreaks, signalling the potential use of vector abundance as a proxy for JE hazard. While vector surveillance data are limited, joint-likelihood models offer a useful approach to inform vector abundance predictions. This study provides decision-makers with a more complete picture of the distribution of JE vector abundance and can be used to target vector surveillance and control efforts and enhance the allocation of resources.

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.

Dynamics of SARS-CoV-2 with waning immunity in the UK population

The dynamics of immunity are crucial to understanding the long-term patterns of the SARS-CoV-2 pandemic. Several cases of reinfection with SARS-CoV-2 have been documented 48-142 days after the initial infection and immunity to seasonal circulating coronaviruses is estimated to be shorter than 1 year. Using an age-structured, deterministic model, we explore potential immunity dynamics using contact data from the UK population. In the scenario where immunity to SARS-CoV-2 lasts an average of three months for non-hospitalized individuals, a year for hospitalized individuals, and the effective reproduction number after lockdown ends is 1.2 (our worst-case scenario), we find that the secondary peak occurs in winter 2020 with a daily maximum of 387 000 infectious individuals and 125 000 daily new cases; threefold greater than in a scenario with permanent immunity. Our models suggest that longitudinal serological surveys to determine if immunity in the population is waning will be most informative when sampling takes place from the end of the lockdown in June until autumn 2020. After this period, the proportion of the population with antibodies to SARS-CoV-2 is expected to increase due to the secondary wave. Overall, our analysis presents considerations for policy makers on the longer-term dynamics of SARS-CoV-2 in the UK and suggests that strategies designed to achieve herd immunity may lead to repeated waves of infection as immunity to reinfection is not permanent. This article is part of the theme issue 'Modelling that shaped the early COVID-19 pandemic response in the UK'.

Predicting non-state terrorism worldwide

Several thousand people die every year worldwide because of terrorist attacks perpetrated by non-state actors. In this context, reliable and accurate short-term predictions of non-state terrorism at the local level are key for policy makers to target preventative measures. Using only publicly available data, we show that predictive models that include structural and procedural predictors can accurately predict the occurrence of non-state terrorism locally and a week ahead in regions affected by a relatively high prevalence of terrorism. In these regions, theoretically informed models systematically outperform models using predictors built on past terrorist events only. We further identify and interpret the local effects of major global and regional terrorism drivers. Our study demonstrates the potential of theoretically informed models to predict and explain complex forms of political violence at policy-relevant scales.

Engagement and adherence trade-offs for SARS-CoV-2 contact tracing

Contact tracing is an important tool for allowing countries to ease lockdown policies introduced to combat SARS-CoV-2. For contact tracing to be effective, those with symptoms must self-report themselves while their contacts must self-isolate when asked. However, policies such as legal enforcement of self-isolation can create trade-offs by dissuading individuals from self-reporting. We use an existing branching process model to examine which aspects of contact tracing adherence should be prioritized. We consider an inverse relationship between self-isolation adherence and self-reporting engagement, assuming that increasingly strict self-isolation policies will result in fewer individuals self-reporting to the programme. We find that policies which increase the average duration of self-isolation, or that increase the probability that people self-isolate at all, at the expense of reduced self-reporting rate, will not decrease the risk of a large outbreak and may increase the risk, depending on the strength of the trade-off. These results suggest that policies to increase self-isolation adherence should be implemented carefully. Policies that increase self-isolation adherence at the cost of self-reporting rates should be avoided.

This article is part of the theme issue ‘Modelling that shaped the early COVID-19 pandemic response in the UK’.

Responsible modelling: Unit testing for infectious disease epidemiology

Infectious disease epidemiology is increasingly reliant on large-scale computation and inference. Models have guided health policy for epidemics including COVID-19 and Ebola and endemic diseases including malaria and tuberculosis. Yet a coding bug may bias results, yielding incorrect conclusions and actions causing avoidable harm. We are ethically obliged to make our code as free of error as possible. Unit testing is a coding method to avoid such bugs, but it is rarely used in epidemiology. We demonstrate how unit testing can handle the particular quirks of infectious disease models and aim to increase the uptake of this methodology in our field.

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.

Predicted Impact of COVID-19 on Neglected Tropical Disease Programs and the Opportunity for Innovation

Due to the COVID-19 pandemic, many key neglected tropical disease (NTD) activities have been postponed. This hindrance comes at a time when the NTDs are progressing towards their ambitious goals for 2030. Mathematical modelling on several NTDs, namely gambiense sleeping sickness, lymphatic filariasis, onchocerciasis, schistosomiasis, soil-transmitted helminthiases (STH), trachoma, and visceral leishmaniasis, shows that the impact of this disruption will vary across the diseases. Programs face a risk of resurgence, which will be fastest in high-transmission areas. Furthermore, of the mass drug administration diseases, schistosomiasis, STH, and trachoma are likely to encounter faster resurgence. The case-finding diseases (gambiense sleeping sickness and visceral leishmaniasis) are likely to have fewer cases being detected but may face an increasing underlying rate of new infections. However, once programs are able to resume, there are ways to mitigate the impact and accelerate progress towards the 2030 goals.

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 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.

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 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.

Evaluating Bayesian spatial methods for modelling species distributions with clumped and restricted occurrence data

Statistical approaches for inferring the spatial distribution of taxa (Species Distribution Models, SDMs) commonly rely on available occurrence data, which is often clumped and geographically restricted. Although available SDM methods address some of these factors, they could be more directly and accurately modelled using a spatially-explicit approach. Software to fit models with spatial autocorrelation parameters in SDMs are now widely available, but whether such approaches for inferring SDMs aid predictions compared to other methodologies is unknown. Here, within a simulated environment using 1000 generated species' ranges, we compared the performance of two commonly used non-spatial SDM methods (Maximum Entropy Modelling, MAXENT and boosted regression trees, BRT), to a spatial Bayesian SDM method (fitted using R-INLA), when the underlying data exhibit varying combinations of clumping and geographic restriction. Finally, we tested how any recommended methodological settings designed to account for spatially non-random patterns in the data impact inference. Spatial Bayesian SDM method was the most consistently accurate method, being in the top 2 most accurate methods in 7 out of 8 data sampling scenarios. Within high-coverage sample datasets, all methods performed fairly similarly. When sampling points were randomly spread, BRT had a 1-3% greater accuracy over the other methods and when samples were clumped, the spatial Bayesian SDM method had a 4%-8% better AUC score. Alternatively, when sampling points were restricted to a small section of the true range all methods were on average 10-12% less accurate, with greater variation among the methods. Model inference under the recommended settings to account for autocorrelation was not impacted by clumping or restriction of data, except for the complexity of the spatial regression term in the spatial Bayesian model. Methods, such as those made available by R-INLA, can be successfully used to account for spatial autocorrelation in an SDM context and, by taking account of random effects, produce outputs that can better elucidate the role of covariates in predicting species occurrence. Given that it is often unclear what the drivers are behind data clumping in an empirical occurrence dataset, or indeed how geographically restricted these data are, spatially-explicit Bayesian SDMs may be the better choice when modelling the spatial distribution of target species.

A generalised random encounter model for estimating animal density with remote sensor data

Wildlife monitoring technology is advancing rapidly and the use of remote sensors such as camera traps and acoustic detectors is becoming common in both the terrestrial and marine environments. Current methods to estimate abundance or density require individual recognition of animals or knowing the distance of the animal from the sensor, which is often difficult. A method without these requirements, the random encounter model (REM), has been successfully applied to estimate animal densities from count data generated from camera traps. However, count data from acoustic detectors do not fit the assumptions of the REM due to the directionality of animal signals.

We developed a generalised REM (gREM), to estimate absolute animal density from count data from both camera traps and acoustic detectors. We derived the gREM for different combinations of sensor detection widths and animal signal widths (a measure of directionality). We tested the accuracy and precision of this model using simulations of different combinations of sensor detection widths and animal signal widths, number of captures and models of animal movement.

We find that the gREM produces accurate estimates of absolute animal density for all combinations of sensor detection widths and animal signal widths. However, larger sensor detection and animal signal widths were found to be more precise. While the model is accurate for all capture efforts tested, the precision of the estimate increases with the number of captures. We found no effect of different animal movement models on the accuracy and precision of the gREM.

We conclude that the gREM provides an effective method to estimate absolute animal densities from remote sensor count data over a range of sensor and animal signal widths. The gREM is applicable for count data obtained in both marine and terrestrial environments, visually or acoustically (e.g. big cats, sharks, birds, echolocating bats and cetaceans). As sensors such as camera traps and acoustic detectors become more ubiquitous, the gREM will be increasingly useful for monitoring unmarked animal populations across broad spatial, temporal and taxonomic scales.