Population and fertility by age and sex for 195 countries and territories, 1950–2017: a systematic analysis for the Global Burden of Disease Study 2017

Background: Population estimates underpin demographic and epidemiological research and are used to track progress on numerous international indicators of health and development. To date, internationally available estimates of population and fertility, although useful, have not been produced with transparent and replicable methods and do not use standardised estimates of mortality. We present single-calendar year and single-year of age estimates of fertility and population by sex with standardised and replicable methods. Methods: We estimated population in 195 locations by single year of age and single calendar year from 1950 to 2017 with standardised and replicable methods. We based the estimates on the demographic balancing equation, with inputs of fertility, mortality, population, and migration data. Fertility data came from 7817 location-years of vital registration data, 429 surveys reporting complete birth histories, and 977 surveys and censuses reporting summary birth histories. We estimated age-specific fertility rates (ASFRs; the annual number of livebirths to women of a specified age group per 1000 women in that age group) by use of spatiotemporal Gaussian process regression and used the ASFRs to estimate total fertility rates (TFRs; the average number of children a woman would bear if she survived through the end of the reproductive age span [age 10–54 years] and experienced at each age a particular set of ASFRs observed in the year of interest). Because of sparse data, fertility at ages 10–14 years and 50–54 years was estimated from data on fertility in women aged 15–19 years and 45–49 years, through use of linear regression. Age-specific mortality data came from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2017 estimates. Data on population came from 1257 censuses and 761 population registry location-years and were adjusted for underenumeration and age misreporting with standard demographic methods. Migration was estimated with the GBD Bayesian demographic balancing model, after incorporating information about refugee migration into the model prior. Final population estimates used the cohort-component method of population projection, with inputs of fertility, mortality, and migration data. Population uncertainty was estimated by use of out-of-sample predictive validity testing. With these data, we estimated the trends in population by age and sex and in fertility by age between 1950 and 2017 in 195 countries and territories. Findings: From 1950 to 2017, TFRs decreased by 49\ub74% (95% uncertainty interval [UI] 46\ub74–52\ub70). The TFR decreased from 4\ub77 livebirths (4\ub75–4\ub79) to 2\ub74 livebirths (2\ub72–2\ub75), and the ASFR of mothers aged 10–19 years decreased from 37 livebirths (34–40) to 22 livebirths (19–24) per 1000 women. Despite reductions in the TFR, the global population has been increasing by an average of 83\ub78 million people per year since 1985. The global population increased by 197\ub72% (193\ub73–200\ub78) since 1950, from 2\ub76 billion (2\ub75–2\ub76) to 7\ub76 billion (7\ub74–7\ub79) people in 2017; much of this increase was in the proportion of the global population in south Asia and sub-Saharan Africa. The global annual rate of population growth increased between 1950 and 1964, when it peaked at 2\ub70%; this rate then remained nearly constant until 1970 and then decreased to 1\ub71% in 2017. Population growth rates in the southeast Asia, east Asia, and Oceania GBD super-region decreased from 2\ub75% in 1963 to 0\ub77% in 2017, whereas in sub-Saharan Africa, population growth rates were almost at the highest reported levels ever in 2017, when they were at 2\ub77%. The global average age increased from 26\ub76 years in 1950 to 32\ub71 years in 2017, and the proportion of the population that is of working age (age 15–64 years) increased from 59\ub79% to 65\ub73%. At the national level, the TFR decreased in all countries and territories between 1950 and 2017; in 2017, TFRs ranged from a low of 1\ub70 livebirths (95% UI 0\ub79–1\ub72) in Cyprus to a high of 7\ub71 livebirths (6\ub78–7\ub74) in Niger. The TFR under age 25 years (TFU25; number of livebirths expected by age 25 years for a hypothetical woman who survived the age group and was exposed to current ASFRs) in 2017 ranged from 0\ub708 livebirths (0\ub707–0\ub709) in South Korea to 2\ub74 livebirths (2\ub72–2\ub76) in Niger, and the TFR over age 30 years (TFO30; number of livebirths expected for a hypothetical woman ageing from 30 to 54 years who survived the age group and was exposed to current ASFRs) ranged from a low of 0\ub73 livebirths (0\ub73–0\ub74) in Puerto Rico to a high of 3\ub71 livebirths (3\ub70–3\ub72) in Niger. TFO30 was higher than TFU25 in 145 countries and territories in 2017. 33 countries had a negative population growth rate from 2010 to 2017, most of which were located in central, eastern, and western Europe, whereas population growth rates of more than 2\ub70% were seen in 33 of 46 countries in sub-Saharan Africa. In 2017, less than 65% of the national population was of working age in 12 of 34 high-income countries, and less than 50% of the national population was of working age in Mali, Chad, and Niger. Interpretation: Population trends create demographic dividends and headwinds (ie, economic benefits and detriments) that affect national economies and determine national planning needs. Although TFRs are decreasing, the global population continues to grow as mortality declines, with diverse patterns at the national level and across age groups. To our knowledge, this is the first study to provide transparent and replicable estimates of population and fertility, which can be used to inform decision making and to monitor progress. Funding: Bill & Melinda Gates Foundation

Predicting the environmental suitability for onchocerciasis in Africa as an aid to elimination planning

Recent evidence suggests that, in some foci, elimination of onchocerciasis from Africa may be feasible with mass drug administration (MDA) of ivermectin. To achieve continental elimination of transmission, mapping surveys will need to be conducted across all implementation units (IUs) for which endemicity status is currently unknown. Using boosted regression tree models with optimised hyperparameter selection, we estimated environmental suitability for onchocerciasis at the 5 × 5-km resolution across Africa. In order to classify IUs that include locations that are environmentally suitable, we used receiver operating characteristic (ROC) analysis to identify an optimal threshold for suitability concordant with locations where onchocerciasis has been previously detected. This threshold value was then used to classify IUs (more suitable or less suitable) based on the location within the IU with the largest mean prediction. Mean estimates of environmental suitability suggest large areas across West and Central Africa, as well as focal areas of East Africa, are suitable for onchocerciasis transmission, consistent with the presence of current control and elimination of transmission efforts. The ROC analysis identified a mean environmental suitability index of 071 as a threshold to classify based on the location with the largest mean prediction within the IU. Of the IUs considered for mapping surveys, 502% exceed this threshold for suitability in at least one 5 × 5-km location. The formidable scale of data collection required to map onchocerciasis endemicity across the African continent presents an opportunity to use spatial data to identify areas likely to be suitable for onchocerciasis transmission. National onchocerciasis elimination programmes may wish to consider prioritising these IUs for mapping surveys as human resources, laboratory capacity, and programmatic schedules may constrain survey implementation, and possibly delaying MDA initiation in areas that would ultimately qualify.SUPPORTING INFORMATION : FIGURE S1. Data coverage by year. Here we visualise the volume of data used in the analysis by country and year. Larger circles indicate more data inputs. ‘NA’ indicates records for which no year was reported (eg, ‘pre-2000’). https://doi.org/10.1371/journal.pntd.0008824.s001FIGURE S2. Illustration of covariate values for year 2000. Maps were produced using ArcGIS Desktop 10.6. https://doi.org/10.1371/journal.pntd.0008824.s002FIGURE S3. Environmental suitability of onchocerciasis including locations that have received MDA for which no pre-intervention data are available. This plot shows suitability predictions from green (low = 0%) to pink (high = 100%), representing those areas where environmental conditions are most similar to prior pathogen detections. Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s003FIGURE S4. Environmental suitability prediction uncertainty including locations that have received MDA for which no pre-intervention data are available. This plot shows uncertainty associated with environmental suitability predictions colored from blue to red (least to most uncertain). Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s004FIGURE S5. Environmental suitability of onchocerciasis excluding morbidity data. This plot shows suitability predictions from green (low = 0%) to pink (high = 100%), representing those areas where environmental conditions are most similar to prior pathogen detections. Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s005FIGURE S6. Environmental suitability prediction uncertainty excluding morbidity data. This plot shows uncertainty associated with environmental suitability predictions colored from blue to red (least to most uncertain). Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. https://doi.org/10.1371/journal.pntd.0008824.s006FIGURE S7. Covariate Effect Curves for all onchocerciasis occurrences (measures of infection prevalence and disability). On the right set of axes we show the frequency density of the occurrences taking covariate values over 20 bins of the horizontal axis. The left set of axes shows the effect of each on the model, where the mean effect is plotted on the black line and its uncertainty is represented by the upper and lower confidence interval bounds plotted in dark grey. The figures show the fit per covariate relative to the data that correspond to specific values of the covariate. https://doi.org/10.1371/journal.pntd.0008824.s007FIGURE S8. Covariate Effect Curves for all onchocerciasis occurrences (measures of infection prevalence and disability). On the right set of axes we show the frequency density of the occurrences taking covariate values over 20 bins of the horizontal axis. The left set of axes shows the effect of each on the model, where the mean effect is plotted on the black line and its uncertainty is represented by the upper and lower confidence interval bounds plotted in dark grey. https://doi.org/10.1371/journal.pntd.0008824.s008FIGURE S9. ROC analysis for threshold. Results of the area under the receiver operating characteristic (ROC) curve analysis are presented below, with false positive rate (FPR) on the x-axis and true positive rate (TPR) on the y-axis. The red dot on the curve represents the location on the curve that corresponds to a threshold that most closely agreed with the input data. For each of the 100 BRT models, we estimated the optimal threshold that maximised agreement between occurrence inputs (considered true positives) and the mean model predictions as 0·71. https://doi.org/10.1371/journal.pntd.0008824.s009TABLE S1. Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER) checklist. https://doi.org/10.1371/journal.pntd.0008824.s010TABLE S2. Total number of occurrence data classified as point and polygon inputs by diagnostic. We present the total number of occurrence points extracted from the input data sources by diagnostic type. ‘Other diagnostics’ include: DEC Patch test; Knott’s Method (Mazotti Test); 2 types of LAMP; blood smears; and urine tests. https://doi.org/10.1371/journal.pntd.0008824.s011TABLE S3. Total number of occurrence data classified as point and polygon inputs by location. https://doi.org/10.1371/journal.pntd.0008824.s012TABLE S4. Covariate information. https://doi.org/10.1371/journal.pntd.0008824.s013TEXT S1. Details outlining construction of occurrence dataset. https://doi.org/10.1371/journal.pntd.0008824.s014TEXT S2. Covariate rationale. https://doi.org/10.1371/journal.pntd.0008824.s015TEXT S3. Boosted regression tree methodology additional details. https://doi.org/10.1371/journal.pntd.0008824.s016APPENDIX S1. Country-level maps and data results. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s017This work was primarily supported by a grant from the Bill & Melinda Gates Foundation OPP1132415 (SIH). Financial support from the Neglected Tropical Disease Modelling Consortium (https://www.ntdmodelling.org/), which is funded by the Bill & Melinda Gates Foundation (grants No. OPP1184344 and OPP1186851), and joint centre funding (grant No. MR/R015600/1) by the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement which is also part of the EDCTP2 programme supported by the European Union (MGB).The Neglected Tropical Disease Modelling Consortium which is funded by the Bill & Melinda Gates Foundation, the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement which is also part of the EDCTP2 programme supported by the European Union (MGB).http://www.plosNTDS.orgam2022Medical Microbiolog

Global, regional, and national burden of suicide mortality 1990 to 2016 : Systematic analysis for the Global Burden of Disease Study 2016

Objectives To use the estimates from the Global Burden of Disease Study 2016 to describe patterns of suicide mortality globally, regionally, and for 195 countries and territories by age, sex, and Socio-demographic index, and to describe temporal trends between 1990 and 2016. Design Systematic analysis. Main outcome measures Crude and age standardised rates from suicide mortality and years of life lost were compared across regions and countries, and by age, sex, and Socio-demographic index (a composite measure of fertility, income, and education). Results The total number of deaths from suicide increased by 6.7% (95% uncertainty interval 0.4% to 15.6%) globally over the 27 year study period to 817 000 (762 000 to 884 000) deaths in 2016. However, the age standardised mortality rate for suicide decreased by 32.7% (27.2% to 36.6%) worldwide between 1990 and 2016, similar to the decline in the global age standardised mortality rate of 30.6%. Suicide was the leading cause of age standardised years of life lost in the Global Burden of Disease region of high income Asia Pacific and was among the top 10 leading causes in eastern Europe, central Europe, western Europe, central Asia, Australasia, southern Latin America, and high income North America. Rates for men were higher than for women across regions, countries, and age groups, except for the 15 to 19 age group. There was variation in the female to male ratio, with higher ratios at lower levels of Socio-demographic index. Women experienced greater decreases in mortality rates (49.0%, 95% uncertainty interval 42.6% to 54.6%) than men (23.8%, 15.6% to 32.7%). Conclusions Age standardised mortality rates for suicide have greatly reduced since 1990, but suicide remains an important contributor to mortality worldwide. Suicide mortality was variable across locations, between sexes, and between age groups. Suicide prevention strategies can be targeted towards vulnerable populations if they are informed by variations in mortality rates. © Published by the BMJ Publishing Group Limited.Peer reviewe