Forecasts and alternative scenarios of COVID-19 mortality have been critical inputs into a range of policies and decision-makers need information about predictive performance. We identified n=386 public COVID-19 forecasting models and included n=8 that were global in scope and provided public, date-versioned forecasts. For each, we examined the median absolute percent error (MAPE) compared to subsequently observed mortality trends, stratified by weeks of extrapolation, world region, and month of model estimation. Models were also assessed for ability to predict the timing of peak daily mortality. The MAPE among models released in July rose from 1.8% at one week of extrapolation to 24.6% at twelve weeks. The MAPE at six weeks were the highest in Sub-Saharan Africa (34.8%), and the lowest in high-income countries (6.3%). At the global level, several models had about 10% MAPE at six weeks, showing surprisingly good performance despite the complexities of modelling human behavioural responses and government interventions. The framework and publicly available codebase presented here ( https://github.com/pyliu47/covidcompare ) can be routinely used to compare predictions and evaluate predictive performance in an ongoing fashion.
Forecasts and alternative scenarios of COVID-19 mortality have been critical inputs for pandemic response efforts, and decision-makers need information about predictive performance. We screen n = 386 public COVID-19 forecasting models, identifying n = 7 that are global in scope and provide public, date-versioned forecasts. We examine their predictive performance for mortality by weeks of extrapolation, world region, and estimation month. We additionally assess prediction of the timing of peak daily mortality. Globally, models released in October show a median absolute percent error (MAPE) of 7 to 13% at six weeks, reflecting surprisingly good performance despite the complexities of modelling human behavioural responses and government interventions. Median absolute error for peak timing increased from 8 days at one week of forecasting to 29 days at eight weeks and is similar for first and subsequent peaks. The framework and public codebase (https://github.com/pyliu47/covidcompare) can be used to compare predictions and evaluate predictive performance going forward.
Despite a long history of mosquito-borne virus epidemics in the Americas, the impact of the Zika virus (ZIKV) epidemic of 2015–2016 was unexpected. The need for scientifically informed decision-making is driving research to understand the emergence and spread of ZIKV. To support that research, we assembled a data set of key covariates for modeling ZIKV transmission dynamics in Colombia, where ZIKV transmission was widespread and the government made incidence data publically available. On a weekly basis between January 1, 2014 and October 1, 2016 at three administrative levels, we collated spatiotemporal Zika incidence data, nine environmental variables, and demographic data into a single downloadable database. These new datasets and those we identified, processed, and assembled at comparable spatial and temporal resolutions will save future researchers considerable time and effort in performing these data processing steps, enabling them to focus instead on extracting epidemiological insights from this important data set. Similar approaches could prove useful for filling data gaps to enable epidemiological analyses of future disease emergence events.
Despite a long history of mosquito-borne virus epidemics in the Americas, the impact of the Zika virus (ZIKV) epidemic of 2015-2016 was unexpected. The need for scientifically informed decision-making is driving research to understand the emergence and spread of ZIKV. To support that research, we assembled a data set of key covariates for modeling ZIKV transmission dynamics in Colombia, where ZIKV transmission was widespread and the government made incidence data publically available. On a weekly basis between January 1, 2014 and October 1, 2016 at three administrative levels, we collated spatiotemporal Zika incidence data, nine environmental variables, and demographic data into a single downloadable database. These new datasets and those we identified, processed, and assembled at comparable spatial and temporal resolutions will save future researchers considerable time and effort in performing these data processing steps, enabling them to focus instead on extracting epidemiological insights from this important data set. Similar approaches could prove useful for filling data gaps to enable epidemiological analyses of future disease emergence events.
Despite a long history of mosquito-borne virus epidemics in the Americas, the impact of the Zika virus (ZIKV) epidemic of 2015–2016 was unexpected. The need for scientifically informed decision-making is driving research to understand the emergence and spread of ZIKV. To support that research, we assembled a data set of key covariates for modeling ZIKV transmission dynamics in Colombia, where ZIKV transmission was widespread and the government made incidence data publically available. On a weekly basis between January 1, 2014 and October 1, 2016 at three administrative levels, we collated spatiotemporal Zika incidence data, nine environmental variables, and demographic data into a single downloadable database. These new datasets and those we identified, processed, and assembled at comparable spatial and temporal resolutions will save future researchers considerable time and effort in performing these data processing steps, enabling them to focus instead on extracting epidemiological insights from this important data set. Similar approaches could prove useful for filling data gaps to enable epidemiological analyses of future disease emergence events.
Despite a long history of mosquito-borne virus epidemics in the Americas, the impact of the Zika virus (ZIKV) epidemic of 2015-2016 was unexpected. The need for scientifically informed decision-making is driving research to understand the emergence and spread of ZIKV. To support that research, we assembled a data set of key covariates for modeling ZIKV transmission dynamics in Colombia, where ZIKV transmission was widespread and the government made incidence data publically available. On a weekly basis between January 1, 2014 and October 1, 2016 at three administrative levels, we collated spatiotemporal Zika incidence data, nine environmental variables, and demographic data into a single downloadable database. These new datasets and those we identified, processed, and assembled at comparable spatial and temporal resolutions will save future researchers considerable time and effort in performing these data processing steps, enabling them to focus instead on extracting epidemiological insights from this important data set. Similar approaches could prove useful for filling data gaps to enable epidemiological analyses of future disease emergence events.
BACKGROUND: The 2018–2019 Ebola virus disease (EVD) outbreak in North Kivu and Ituri provinces in the Democratic Republic of the Congo (DRC) is the largest ever recorded in the DRC. It has been declared a Public Health Emergency of International Concern. The outbreak emerged in a region of chronic conflict and insecurity, and directed attacks against health care workers may have interfered with disease response activities. Our study characterizes and quantifies the broader conflict dynamics over the course of the outbreak by pairing epidemiological and all available spatial conflict data. METHODS: We build a set of conflict variables by mapping the spatial locations of all conflict events and their associated deaths in each of the affected health zones in North Kivu and Ituri, eastern DRC, before and during the outbreak. Using these data, we compare patterns of conflict before and during the outbreak in affected health zones and those not affected. We then test whether conflict is correlated with increased EVD transmission at the health zone level. FINDINGS: The incidence of conflict events per capita is ~ 600 times more likely in Ituri and North Kivu than for the rest of the DRC. We identified 15 time periods of substantial uninterrupted transmission across 11 health zones and a total of 120 bi-weeks. We do not find significant short-term associations between the bi-week reproduction numbers and the number of conflicts. However, we do find that the incidence of conflict per capita was correlated with the incidence of EVD per capita at the health zone level for the entire outbreak (Pearson's r = 0.33, 95% CI 0.05–0.57). In the two provinces, the monthly number of conflict events also increased by a factor of 2.7 in Ebola-affected health zones (p value < 0.05) compared to 2.0 where no transmission was reported and 1.3 in the rest of the DRC, in the period between February 2019 and July 2019. CONCLUSION: We characterized the association between variables documenting broad conflict levels and EVD transmission. ...
BACKGROUND: Wolbachia-infected mosquitoes reduce dengue virus transmission, and city-wide releases in Yogyakarta city, Indonesia, are showing promising entomological results. Accurate estimates of the burden of dengue, its spatial distribution and the potential impact of Wolbachia are critical in guiding funder and government decisions on its future wider use. METHODS: Here, we combine multiple modelling methods for burden estimation to predict national case burden disaggregated by severity and map the distribution of burden across the country using three separate data sources. An ensemble of transmission models then predicts the estimated reduction in dengue transmission following a nationwide roll-out of wMel Wolbachia. RESULTS: We estimate that 7.8 million (95% uncertainty interval [UI] 1.8-17.7 million) symptomatic dengue cases occurred in Indonesia in 2015 and were associated with 332,865 (UI 94,175-754,203) lost disability-adjusted life years (DALYs). The majority of dengue's burden was due to non-severe cases that did not seek treatment or were challenging to diagnose in outpatient settings leading to substantial underreporting. Estimated burden was highly concentrated in a small number of large cities with 90% of dengue cases occurring in 15.3% of land area. Implementing a nationwide Wolbachia population replacement programme was estimated to avert 86.2% (UI 36.2-99.9%) of cases over a long-term average. CONCLUSIONS: These results suggest interventions targeted to the highest burden cities can have a disproportionate impact on dengue burden. Area-wide interventions, such as Wolbachia, that are deployed based on the area covered could protect people more efficiently than individual-based interventions, such as vaccines, in such dense environments.
BACKGROUND: Yellow fever cases are under-reported and the exact distribution of the disease is unknown. An effective vaccine is available but more information is needed about which populations within risk zones should be targeted to implement interventions. Substantial outbreaks of yellow fever in Angola, Democratic Republic of the Congo, and Brazil, coupled with the global expansion of the range of its main urban vector, Aedes aegypti, suggest that yellow fever has the propensity to spread further internationally. The aim of this study was to estimate the disease's contemporary distribution and potential for spread into new areas to help inform optimal control and prevention strategies. METHODS: We assembled 1155 geographical records of yellow fever virus infection in people from 1970 to 2016. We used a Poisson point process boosted regression tree model that explicitly incorporated environmental and biological explanatory covariates, vaccination coverage, and spatial variability in disease reporting rates to predict the relative risk of apparent yellow fever virus infection at a 5 × 5 km resolution across all risk zones (47 countries across the Americas and Africa). We also used the fitted model to predict the receptivity of areas outside at-risk zones to the introduction or reintroduction of yellow fever transmission. By use of previously published estimates of annual national case numbers, we used the model to map subnational variation in incidence of yellow fever across at-risk countries and to estimate the number of cases averted by vaccination worldwide. FINDINGS: Substantial international and subnational spatial variation exists in relative risk and incidence of yellow fever as well as varied success of vaccination in reducing incidence in several high-risk regions, including Brazil, Cameroon, and Togo. Areas with the highest predicted average annual case numbers include large parts of Nigeria, the Democratic Republic of the Congo, and South Sudan, where vaccination coverage in 2016 was estimated to be substantially less than the recommended threshold to prevent outbreaks. Overall, we estimated that vaccination coverage levels achieved by 2016 avert between 94 336 and 118 500 cases of yellow fever annually within risk zones, on the basis of conservative and optimistic vaccination scenarios. The areas outside at-risk regions with predicted high receptivity to yellow fever transmission (eg, parts of Malaysia, Indonesia, and Thailand) were less extensive than the distribution of the main urban vector, A aegypti, with low receptivity to yellow fever transmission in southern China, where A aegypti is known to occur. INTERPRETATION: Our results provide the evidence base for targeting vaccination campaigns within risk zones, as well as emphasising their high effectiveness. Our study highlights areas where public health authorities should be most vigilant for potential spread or importation events. FUNDING: Bill & Melinda Gates Foundation.
BACKGROUND: Substantial outbreaks of yellow fever in Angola and Brazil in the past 2 years, combined with global shortages in vaccine stockpiles, highlight a pressing need to assess present control strategies. The aims of this study were to estimate global yellow fever vaccination coverage from 1970 through to 2016 at high spatial resolution and to calculate the number of individuals still requiring vaccination to reach population coverage thresholds for outbreak prevention. METHODS: For this adjusted retrospective analysis, we compiled data from a range of sources (eg, WHO reports and health-service-provider registeries) reporting on yellow fever vaccination activities between May 1, 1939, and Oct 29, 2016. To account for uncertainty in how vaccine campaigns were targeted, we calculated three population coverage values to encompass alternative scenarios. We combined these data with demographic information and tracked vaccination coverage through time to estimate the proportion of the population who had ever received a yellow fever vaccine for each second level administrative division across countries at risk of yellow fever virus transmission from 1970 to 2016. FINDINGS: Overall, substantial increases in vaccine coverage have occurred since 1970, but notable gaps still exist in contemporary coverage within yellow fever risk zones. We estimate that between 393·7 million and 472·9 million people still require vaccination in areas at risk of yellow fever virus transmission to achieve the 80% population coverage threshold recommended by WHO; this represents between 43% and 52% of the population within yellow fever risk zones, compared with between 66% and 76% of the population who would have required vaccination in 1970. INTERPRETATION: Our results highlight important gaps in yellow fever vaccination coverage, can contribute to improved quantification of outbreak risk, and help to guide planning of future vaccination efforts and emergency stockpiling. FUNDING: The Rhodes Trust, Bill & Melinda Gates Foundation, the Wellcome Trust, the National Library of Medicine of the National Institutes of Health, the European Union's Horizon 2020 research and innovation programme.
BACKGROUND: Since late 2015, an epidemic of yellow fever has caused more than 7334 suspected cases in Angola and the Democratic Republic of the Congo, including 393 deaths. We sought to understand the spatial spread of this outbreak to optimise the use of the limited available vaccine stock. METHODS: We jointly analysed datasets describing the epidemic of yellow fever, vector suitability, human demography, and mobility in central Africa to understand and predict the spread of yellow fever virus. We used a standard logistic model to infer the district-specific yellow fever virus infection risk during the course of the epidemic in the region. FINDINGS: The early spread of yellow fever virus was characterised by fast exponential growth (doubling time of 5-7 days) and fast spatial expansion (49 districts reported cases after only 3 months) from Luanda, the capital of Angola. Early invasion was positively correlated with high population density (Pearson's r 0·52, 95% CI 0·34-0·66). The further away locations were from Luanda, the later the date of invasion (Pearson's r 0·60, 95% CI 0·52-0·66). In a Cox model, we noted that districts with higher population densities also had higher risks of sustained transmission (the hazard ratio for cases ceasing was 0·74, 95% CI 0·13-0·92 per log-unit increase in the population size of a district). A model that captured human mobility and vector suitability successfully discriminated districts with high risk of invasion from others with a lower risk (area under the curve 0·94, 95% CI 0·92-0·97). If at the start of the epidemic, sufficient vaccines had been available to target 50 out of 313 districts in the area, our model would have correctly identified 27 (84%) of the 32 districts that were eventually affected. INTERPRETATION: Our findings show the contributions of ecological and demographic factors to the ongoing spread of the yellow fever outbreak and provide estimates of the areas that could be prioritised for vaccination, although other constraints such as vaccine supply and delivery need to be accounted for before such insights can be translated into policy. FUNDING: Wellcome Trust.
Background Since late 2015, an epidemic of yellow fever has caused more than 7334 suspected cases in Angola and the Democratic Republic of the Congo, including 393 deaths. We sought to understand the spatial spread of this outbreak to optimise the use of the limited available vaccine stock. Methods We jointly analysed datasets describing the epidemic of yellow fever, vector suitability, human demography, and mobility in central Africa to understand and predict the spread of yellow fever virus. We used a standard logistic model to infer the district-specific yellow fever virus infection risk during the course of the epidemic in the region. Findings The early spread of yellow fever virus was characterised by fast exponential growth (doubling time of 5–7 days) and fast spatial expansion (49 districts reported cases after only 3 months) from Luanda, the capital of Angola. Early invasion was positively correlated with high population density (Pearson's r 0·52, 95% CI 0·34–0·66). The further away locations were from Luanda, the later the date of invasion (Pearson's r 0·60, 95% CI 0·52–0·66). In a Cox model, we noted that districts with higher population densities also had higher risks of sustained transmission (the hazard ratio for cases ceasing was 0·74, 95% CI 0·13–0·92 per log-unit increase in the population size of a district). A model that captured human mobility and vector suitability successfully discriminated districts with high risk of invasion from others with a lower risk (area under the curve 0·94, 95% CI 0·92–0·97). If at the start of the epidemic, sufficient vaccines had been available to target 50 out of 313 districts in the area, our model would have correctly identified 27 (84%) of the 32 districts that were eventually affected. Interpretation Our findings show the contributions of ecological and demographic factors to the ongoing spread of the yellow fever outbreak and provide estimates of the areas that could be prioritised for vaccination, although other constraints such as vaccine supply and delivery need to be accounted for before such insights can be translated into policy. Funding Wellcome Trust. ; ISSN:1473-3099 ; ISSN:1474-4457
Achieving universal health coverage (UHC) requires health financing systems that provide prepaid pooled resources for key health services without placing undue financial stress on households. Understanding current and future trajectories of health financing is vital for progress towards UHC. We used historical health financing data for 188 countries from 1995 to 2015 to estimate future scenarios of health spending and pooled health spending through to 2040.We extracted historical data on gross domestic product (GDP) and health spending for 188 countries from 1995 to 2015, and projected annual GDP, development assistance for health, and government, out-of-pocket, and prepaid private health spending from 2015 through to 2040 as a reference scenario. These estimates were generated using an ensemble of models that varied key demographic and socioeconomic determinants. We generated better and worse alternative future scenarios based on the global distribution of historic health spending growth rates. Last, we used stochastic frontier analysis to investigate the association between pooled health resources and UHC index, a measure of a country's UHC service coverage. Finally, we estimated future UHC performance and the number of people covered under the three future scenarios.
BACKGROUND: The number of individuals living with dementia is increasing, negatively affecting families, communities, and health-care systems around the world. A successful response to these challenges requires an accurate understanding of the dementia disease burden. We aimed to present the first detailed analysis of the global prevalence, mortality, and overall burden of dementia as captured by the Global Burden of Diseases, Injuries, and Risk Factors (GBD) Study 2016, and highlight the most important messages for clinicians and neurologists. METHODS: GBD 2016 obtained data on dementia from vital registration systems, published scientific literature and surveys, and data from health-service encounters on deaths, excess mortality, prevalence, and incidence from 195 countries and territories from 1990 to 2016, through systematic review and additional data-seeking efforts. To correct for differences in cause of death coding across time and locations, we modelled mortality due to dementia using prevalence data and estimates of excess mortality derived from countries that were most likely to code deaths to dementia relative to prevalence. Data were analysed by standardised methods to estimate deaths, prevalence, years of life lost (YLLs), years of life lived with disability (YLDs), and disability-adjusted life-years (DALYs; computed as the sum of YLLs and YLDs), and the fractions of these metrics that were attributable to four risk factors that met GBD criteria for assessment (high body-mass index [BMI], high fasting plasma glucose, smoking, and a diet high in sugar-sweetened beverages). FINDINGS: In 2016, the global number of individuals who lived with dementia was 43·8 million (95% uncertainty interval [UI] 37·8-51·0), increased from 20.2 million (17·4-23·5) in 1990. This increase of 117% (95% UI 114-121) contrasted with a minor increase in age-standardised prevalence of 1·7% (1·0-2·4), from 701 cases (95% UI 602-815) per 100 000 population in 1990 to 712 cases (614-828) per 100 000 population in 2016. More women than men had dementia in 2016 (27·0 million, 95% UI 23·3-31·4, vs 16.8 million, 14.4-19.6), and dementia was the fifth leading cause of death globally, accounting for 2·4 million (95% UI 2·1-2·8) deaths. Overall, 28·8 million (95% UI 24·5-34·0) DALYs were attributed to dementia; 6·4 million (95% UI 3·4-10·5) of these could be attributed to the modifiable GBD risk factors of high BMI, high fasting plasma glucose, smoking, and a high intake of sugar-sweetened beverages. INTERPRETATION: The global number of people living with dementia more than doubled from 1990 to 2016, mainly due to increases in population ageing and growth. Although differences in coding for causes of death and the heterogeneity in case-ascertainment methods constitute major challenges to the estimation of the burden of dementia, future analyses should improve on the methods for the correction of these biases. Until breakthroughs are made in prevention or curative treatment, dementia will constitute an increasing challenge to health-care systems worldwide. FUNDING: Bill & Melinda Gates Foundation. ; AA received financial support from the Department of Science and Technology, Government of India, (New Delhi, India) through the INSPIRE Faculty program. MSBS received Australian Government Research and Training Program funding for post-graduates to study at the Australian National University (Canberra, ACT, Australia). FC acknowledges support from the European Union (FEDER funds POCI/01/0145/FEDER/007728 and POCI/01/0145/FEDER/007265) and National Funds (FCT/MEC, Fundação para a Ciência e a Tecnologia and Ministério da Educação e Ciência) under the Partnership Agreements PT2020 UID/MULTI/04378/2013 and PT2020 UID/QUI/50006/2013. EC is supported by an Australian Research Council Future fellowship (FT3 140100085). AK was supported by the Miguel Servet contract financed by the CP13/00150 and PI15/00862 projects, integrated into the National R + D + I and funded by the ISCIII (General Branch Evaluation and Promotion of Health Research) and the European Regional Development fund (ISCIII-FEDER). MOO is supported by grant U54HG007479 from the National Institutes of Health. TCR is a member of the Alzheimer Scotland Dementia Research Centre (University of Edinburgh, Edinburgh, UK) and is supported by Alzheimer Scotland. RT-S was partly supported by grant number PROMETEOII/2015/021 from Generalitat Valenciana and the national grant PI17/00719 from ISCIII-FEDER. TW acknowledges academic support from University of Rajarata (Mihintale, Sri Lanka). ; Sí
BACKGROUND: Traumatic brain injury (TBI) and spinal cord injury (SCI) are increasingly recognised as global health priorities in view of the preventability of most injuries and the complex and expensive medical care they necessitate. We aimed to measure the incidence, prevalence, and years of life lived with disability (YLDs) for TBI and SCI from all causes of injury in every country, to describe how these measures have changed between 1990 and 2016, and to estimate the proportion of TBI and SCI cases caused by different types of injury. METHODS: We used results from the Global Burden of Diseases, Injuries, and Risk Factors (GBD) Study 2016 to measure the global, regional, and national burden of TBI and SCI by age and sex. We measured the incidence and prevalence of all causes of injury requiring medical care in inpatient and outpatient records, literature studies, and survey data. By use of clinical record data, we estimated the proportion of each cause of injury that required medical care that would result in TBI or SCI being considered as the nature of injury. We used literature studies to establish standardised mortality ratios and applied differential equations to convert incidence to prevalence of long-term disability. Finally, we applied GBD disability weights to calculate YLDs. We used a Bayesian meta-regression tool for epidemiological modelling, used cause-specific mortality rates for non-fatal estimation, and adjusted our results for disability experienced with comorbid conditions. We also analysed results on the basis of the Socio-demographic Index, a compound measure of income per capita, education, and fertility. FINDINGS: In 2016, there were 27·08 million (95% uncertainty interval [UI] 24·30-30·30 million) new cases of TBI and 0·93 million (0·78-1·16 million) new cases of SCI, with age-standardised incidence rates of 369 (331-412) per 100 000 population for TBI and 13 (11-16) per 100 000 for SCI. In 2016, the number of prevalent cases of TBI was 55·50 million (53·40-57·62 million) and of SCI was 27·04 million (24·98-30·15 million). From 1990 to 2016, the age-standardised prevalence of TBI increased by 8·4% (95% UI 7·7 to 9·2), whereas that of SCI did not change significantly (-0·2% [-2·1 to 2·7]). Age-standardised incidence rates increased by 3·6% (1·8 to 5·5) for TBI, but did not change significantly for SCI (-3·6% [-7·4 to 4·0]). TBI caused 8·1 million (95% UI 6·0-10·4 million) YLDs and SCI caused 9·5 million (6·7-12·4 million) YLDs in 2016, corresponding to age-standardised rates of 111 (82-141) per 100 000 for TBI and 130 (90-170) per 100 000 for SCI. Falls and road injuries were the leading causes of new cases of TBI and SCI in most regions. INTERPRETATION: TBI and SCI constitute a considerable portion of the global injury burden and are caused primarily by falls and road injuries. The increase in incidence of TBI over time might continue in view of increases in population density, population ageing, and increasing use of motor vehicles, motorcycles, and bicycles. The number of individuals living with SCI is expected to increase in view of population growth, which is concerning because of the specialised care that people with SCI can require. Our study was limited by data sparsity in some regions, and it will be important to invest greater resources in collection of data for TBI and SCI to improve the accuracy of future assessments. FUNDING: Bill & Melinda Gates Foundation. ; Bill & Melinda Gates Foundation ; We acknowledge the funding and support of the Bill & Melinda Gates Foundation. AK was supported by the Miguel Servet contract, which was financed by the CP13/00150 and PI15/00862 projects integrated into the National Research, Development, and Implementation,and funded by the Instituto de Salud Carlos III General Branch Evaluation and Promotion of Health Research and the European Regional Development Fund (ERDF-FEDER). AMS is supported by the Egyptian Fulbright Mission Program. AF acknowledges the Federal University of Sergipe (Sergipe, Brazil). AA received financial assistance from the Indian Department of Science and Technology (New Delhi, India) through the INSPIRE faculty programme. AS is supported by Health Data Research UK. DJS is supported by the South African Medical Research Council. AB is supported by the Public Health Agency of Canada. SMSI received a senior research fellowship from the Institute for Physical Activity and Nutrition, Deakin University (Waurn Ponds, VIC, Australia), and a career transition grant from the High Blood Pressure Research Council of Australia. FP and CF acknowledge support from the European Union (FEDER funds POCI/01/0145/FEDER/007728 and POCI/01/0145/FEDER/007265) and National Funds (FCT/MEC, Fundação para a Ciência e a Tecnologia, and Ministério da Educação e Ciência) under the Partnership Agreements PT2020 UID/MULTI/04378/2013 and PT2020 UID/QUI/50006/2013. TB acknowledges financial support from the Institute of Medical Research and Medicinal Plant Studies, Yaoundé, Cameroon. AM of Imperial College London is grateful for support from the Northwest London National Institute for Health Research (NIHR) Collaboration for Leadership in Applied Health Research andCare and the Imperial NIHR Biomedical Research Centre. KD is funded by a Wellcome Trust Intermediate Fellowship in Public Health and Tropical Medicine (grant number 201900). PSA is supported by an Australian National Health and Medical Research Council Early Career Fellowship. RT-S was supported in part by grant number PROMETEOII/2015/021 from Generalitat Valenciana and the national grant PI17/00719 from ISCIII-FEDER. The Serbian part of this contribution (by MJ) has been co-financed with grant OI175014 from the Serbian Ministry of Education, Science and Technological Development; publication of results was not contingent upon the Ministry's approval. MMMSM acknowledges support from the Serbian Ministry of Education, Science and Technological Development (contract 175087). MM's research was supported by the NIHR Biomedical Research Centre at Guy's and St Thomas' NHS Foundation Trust (London, UK) and King's College London. The views expressed are those of the authors and not necessarily those of the UK National Health Service, the NIHR, or the UK Department of Health. TWB was supported by the Alexander von Humboldt Foundation through the Alexander von Humboldt professor award, which was funded by the German Federal Ministry of Education and Research ; Sí