Article Text

Download PDFPDF

Original research
Global and regional burden of chronic respiratory disease in 2016 arising from non-infectious airborne occupational exposures: a systematic analysis for the Global Burden of Disease Study 2016
  1. GBD 2016 Occupational Chronic Respiratory Risk Factors Collaborators
    1. Correspondence to Tim Driscoll, School of Public Health, University of Sydney, Sydney, NSW 2006, Australia; tim.driscoll{at}


    Objectives This paper presents detailed analysis of the global and regional burden of chronic respiratory disease arising from occupational airborne exposures, as estimated in the Global Burden of Disease 2016 study.

    Methods The burden of chronic obstructive pulmonary disease (COPD) due to occupational exposure to particulate matter, gases and fumes, and secondhand smoke, and the burden of asthma resulting from occupational exposure to asthmagens, was estimated using the population attributable fraction (PAF), calculated using exposure prevalence and relative risks from the literature. PAFs were applied to the number of deaths and disability-adjusted life years (DALYs) for COPD and asthma. Pneumoconioses were estimated directly from cause of death data. Age-standardised rates were based only on persons aged 15 years and above.

    Results The estimated PAFs (based on DALYs) were 17% (95% uncertainty interval (UI) 14%–20%) for COPD and 10% (95% UI 9%–11%) for asthma. There were estimated to be 519 000 (95% UI 441,000–609,000) deaths from chronic respiratory disease in 2016 due to occupational airborne risk factors (COPD: 460,100 [95% UI 382,000–551,000]; asthma: 37,600 [95% UI 28,400–47,900]; pneumoconioses: 21,500 [95% UI 17,900–25,400]. The equivalent overall burden estimate was 13.6 million (95% UI 11.9–15.5 million); DALYs (COPD: 10.7 [95% UI 9.0–12.5] million; asthma: 2.3 [95% UI 1.9–2.9] million; pneumoconioses: 0.58 [95% UI 0.46–0.67] million). Rates were highest in males; older persons and mainly in Oceania, Asia and sub-Saharan Africa; and decreased from 1990 to 2016.

    Conclusions Workplace exposures resulting in COPD, asthma and pneumoconiosis continue to be important contributors to the burden of disease in all regions of the world. This should be reducible through improved prevention and control of relevant exposures.

    • occupational exposure
    • COPD
    • occupational asthma
    • pneumoconiosis
    • work

    This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See:

    View Full Text

    Statistics from

    Key messages

    What is already known about this subject?

    • Occupational respiratory exposures have been shown to be an important cause of chronic work-related respiratory disease at national and global level.

    • The last analysis of this issue at the global level was for the year 2000—this paper provides a new analysis for 2016.

    What are the new findings?

    • Analysis of the Global Burden of Disease data set suggests that globally there were about 519,100 deaths and 13.6 million disability-adjusted life years in 2016 from chronic respiratory disease due to occupational airborne exposures.

    • The population attributable fraction for chronic obstructive pulmonary disease (COPD) was 17% and for asthma was 10%.

    • Workplace exposures resulting in COPD, asthma and pneumoconioses remain important contributors to the burden of disease in all regions of the world.

    How might this impact on policy or clinical practice in the foreseeable future?

    • These findings highlight opportunities to continue to reduce chronic respiratory disease burden worldwide by improving prevention and control of workplace airborne exposures.


    Airborne respiratory hazards (inorganic and organic particulate matter, vapours, gases and fumes) are a common exposure in occupational settings and many studies have identified resulting malignant and chronic respiratory disease as an important component of the occupational injury and disease burden at both country and global levels.1–10 Work-related respiratory diseases remain a problem even in high-income countries, as shown by incident cases of pneumoconioses that are still occurring.11 12

    The Global Burden of Disease (GBD) Comparative Risk Assessment (CRA) project was the first to consider the burden of occupational chronic respiratory disease comprehensively at a regional and global level, estimating the burden for the year 2000.13 That study included airborne exposures leading to asthma, chronic obstructive pulmonary disease (COPD), asbestosis, coal workers’ pneumoconiosis (CWP) and silicosis.14

    The new GBD initiative, conducted by the Institute of Health Metrics and Evaluation, first focused on 201015 and has been updated several times since.16–19 It provides a detailed analysis of the burden of disease and injury overall and arising from specific risk factors. One set of those risk factors comprises occupational exposures,18but no detailed analysis of the occupational risk factor results has been undertaken.

    This paper presents a detailed analysis of the global and regional burden of chronic respiratory disease arising from non-infectious occupational airborne exposures, as estimated in the GBD 2016 study. Malignant occupational respiratory disease20 and an overview of all occupational risk factors21 are considered in companion papers.20


    General approach

    The general methodology used in GBD 2016 is described elsewhere,18 as is the overall approach to occupational risk factors.21 These methods are briefly summarised here. A more detailed description is provided here of the analyses of occupational exposures to particulate matter, gases and fumes (PMGF), secondhand smoke (SHS), asthmagens and pneumoconiotic dusts and their associated outcomes.

    The burden of occupational respiratory disease for PMGF and SHS (causing COPD) and for asthmagens was estimated using the population attributable fraction (PAF), that is, the proportion of deaths or disability-adjusted life years (DALYs) that would not have occurred if exposure was at the theoretical minimum risk exposure level; this was then used to estimate attributable numbers of deaths or DALYs. The PAF requires information on the relative risk of the disease due to the exposure of interest and the proportion of the target population exposed. Pneumoconioses were estimated directly as part of the overall GBD estimates of prevalence and deaths for each included cause. Age-standardised rates (per 100 000 people) were based only on persons aged 15 and above. Results were calculated for all years from 1990 to 2016, inclusive; the 2016 findings are the focus of this paper. The socio-demographic index (SDI) is a composite indicator of development status based on total fertility rate, mean education for those aged 15 and older and lag distributed income per capita.18Region-specific, SDI-specific and global results are reported here. Country-specific information is available through the GBD Compare data visualisation.22 High-income countries were defined as countries in the Australasia, high-income North America, Western Europe and Asia Pacific regions, and low/middle-income (LMI) countries as all other countries. Employment data came from the International Labour Organization Labour Force,23 supplemented where necessary by sub-national data sources and modelling. PAFs for all carcinogens except asbestos were estimated for each age-sex-country group using the equation based on Levin24:

    Embedded Image

    where P(x) is the proportion of persons exposed at level x in the relevant population and RR(x) is the relative risk corresponding to exposure level x.

    PMGF and SHS

    Industry was used as a proxy for exposure to PMGF because we identified no suitable and valid data sources at a country or global level of exposure to PMGF, either singly or to PMGF as a group. Current industry was used as the basis of exposure estimates, but the estimates of proportions exposed (ie, workers who experienced more than trivial exposure) within each industry (nine categories—see online supplementary table S1) were designed to take into account past exposure (to estimate ever exposed), given that both past and current exposure appear to increase the risk of COPD. Estimates of proportion exposed at lower and higher levels in high income and LMI countries were based on sparse published data (see online supplementary material) and expert opinion by GBD collaborators (online supplementary table S1). Information on risk was obtained by conducting a systematic review of international literature and meta-analysis (unpublished) of relevant results.5 25 Relative risks in these studies were for COPD greater than or equal toGlobal Initiative for Chronic Obstructive Lung Disease(GOLD) stage II: defined as requiring non-reversibility after using bronchodilators for provocation, a forced expiratory volume in one second/forced vital capacity (FEV1/FVC) ratio of less than 0.70 and an FEV1 of less than 80% predicted.26 Relative risk estimates were used for an overall ‘lower’ level (RR=1.44; 95% CI 1.07–1.95) and an overall ‘higher’ level (RR=2.31; 95% CI 1.45 to 3.73) of exposure to the agents of concern (‘higher’ and ‘lower’ were based on the exposure descriptions in the papers).5 25 The reference group was persons not working and persons working in trade, finance or service industries. The prevalence of exposure to PMGF was determined using the following equation:

    Supplementary data

    Embedded Image

    where EAP = economically active population, c = country, s = sex, EA = economic activity, l = level of exposure, y =year and a = age.

    Exposure information on SHS was based on the CAREX (Carcinogen Exposure) database, which provides industry-specific information from 1990 to 1993 on the prevalence of exposure to various carcinogens in countries of Western Europe,27 as described elsewhere.20 The relevant relative risks were those used for SHS in the general GBD 2016 analysis.18


    Exposure and relative risks for asthmagens were based on the current occupation distribution (eight categories—see online supplementary table S2) because there were no suitable and valid data sources at a country or global level describing exposure to the wide range of occupational asthmagens. All relative risk information, except that for agricultural occupations, came from a study by Karjalainen and coworkers, a comprehensive national population study of incident asthma.7 8 Relative risks for agricultural occupations were based on a study by Kogevinas and coworkers,28 using a weighted average of the separate estimates for ‘farmers’ and ‘agricultural’ workers provided in the paper. This information was used because the results were thought to be more generalisable to agriculture in the rest of the world, especially for LMI regions. Separate risks were available and used for males and females (except for agricultural operations), although the sex-specific risks were similar and within the limits of random variation. The same relative risks were used for all age groups. The counterfactual was persons not working and administrative workers. Byssinosis was included as asthma for the purposes of the analysis. The prevalence of exposure to asthmagens was determined using the following equation:

    Embedded Image

    where EAP = economically active population, c = country, s = sex, OCC = occupation, y = year and a = age.

    Pneumoconiotic dusts

    As mentioned, pneumoconioses were estimated directly as part of the overall GBD estimates of prevalence and deaths for each included cause, rather than using the attributable fraction approach. The methods used are described elsewhere.29 The attributable fraction is essentially 100% because virtually all pneumoconioses arise as a result of occupational exposure. Separate estimates were available for silicosis, asbestosis and CWP, with the remaining cases grouped under an ‘other pneumoconiosis’ category.23

    Statistical approach

    The main modelling and analyses employed to produce the GBD 2016 data, and the calculation and use of 95% uncertainty intervals (95% UI), were as described elsewhere.18 21 Uncertainty intervals are primarily presented in detail in the tables to assist with the flow of the text.


    There were estimated to be about 519 000 (95% UI 441 000–609 000) deaths from chronic respiratory disease in 2016 due to occupational airborne risk factors. The vast majority (460 000 [95% UI 382 000–551 000]; 89%) of these were due to COPD arising from PMGF and SHS. The remaining deaths were from asthma (37 600 [95% UI 28 400–47 900]; 7%), due to exposure to a range of asthmagens, and from pneumoconiosis (21 500 [95% UI 17 900–25 400]; 4%), arising from exposure to pneumoconiotic dusts. Males accounted for 75% (390 000) of the deaths overall and between 69% (asthmagens) and 88% (pneumoconiotic dusts) for individual risk factors. The relative contribution of the different risk factors was similar when the burden was measured in terms of DALYs (13.6 [95% UI 11.9–15.5] million DALYs overall; COPD: 10.7 [95% UI 9.0–12.5] million; asthma: 2.3 [95% UI 1.9–2.9] million; pneumoconioses: 0.58 [95% UI 0.49–0.67] million), with 79% of the DALYs due to PMGF and SHS. Males accounted for 73% (9.9 million) of the DALYs (table 1).

    Table 1

    Global occupational-attributable deaths, DALYs and PAFs from chronic respiratory disease due to airborne exposures by risk factor and sex, 2016 (number, percent and proportion [95% uncertainty interval])*

    PMGF and SHS

    The PAF for COPD arising from occupational exposures was 17% (95% UI 14%–20%) for DALYs (16% for deaths), ranging from a low of 10% in Central sub-Saharan Africa to 21% in East Asia (table 2). The PAF was much higher in males (21%) than females (11%) and peaked at about 24% in 60–64-year-old males. The highest number of deaths and the highest rate of deaths from COPD due to occupational exposures occurred in the older age groups, often beyond usual retirement age. Males had three to four times the number and rate of deaths compared with females. The peak for number of DALYs occurred in the 65–74 year age group, but the rate of DALYs increased considerably with age and was highest in the 75–84 year group for both males and females (online supplementary figures S1-S4).

    Table 2

    Deaths, YLLs, YLDs, DALYs and PAFs from COPD due to occupational exposure to PMGF and SHS, by region, 2016 (number, percent, rate, proportion)*

    By far the highest number of deaths and DALYs from COPD occurred in East Asia and South Asia, the regions with the largest populations, which together accounted for about 71% of both measures. The highest rates of deaths were in Oceania, South Asia and East Asia, the rates in these three regions being considerably higher than elsewhere (the lowest rates were in high-income Asia Pacific and Eastern Europe).

    The same regions had the highest DALY rates (the lowest DALY rates were in Andean Latin America, high-income Asia Pacific and the Caribbean). Rates tended to be higher in low-middle and middle SDI regions, and there was considerable variation between regions. Seventy-five percent of the DALYs were due to years of life lost (YLLs); this predominance of YLLs was seen in nearly all regions (table 2). Information on the separate contribution of SHS to COPD is presented in the online supplementary material.


    The PAF for asthma from occupational exposures was estimated to be 9.9% (95% UI 9.0%–10.9%) based on DALYs (8.9% [95% UI 7.8–10.1%] based on deaths), ranging from 4.1% in Central sub-Saharan Africa to 12.0% in South Asia (table 3). The PAF was higher for males (13%) than females (7%) and peaked at around 18% between the ages of about 35 and 49 years.

    Table 3

    Deaths, YLLs, YLDs, DALYs and PAFs from asthma due to occupational exposure to asthmagens, by region, 2016 (number, percent, rate, proportion)*

    Deaths arising from occupational exposure to asthmagens occurred at all ages from 15 to 79 years, but the highest numbers occurred in persons aged 55–64 and the highest rates in persons aged 65–79. The burden was spread more evenly across age groups in terms of DALYs, with the highest number of DALYs in the 45–54 year age group and the highest rates in the 55–64 year age group (online supplementary figures S5-S8).

    The highest number of deaths occurred in South Asia and Southeast Asia, and rates were highest in the low and low-middle SDI regions, particularly Oceania, South Asia and Southeast Asia (the lowest rates were in Western Europe and Central Europe). A similar pattern was seen for DALYs (the lowest rates were in East Asia and Tropical Latin America). Overall, YLLs and years of life with disability each contributed about 50% to the DALYs. However, low and low-middle SDI regions had a much higher proportion of DALYs due to YLLs compared with high-income regions, reflecting that a higher proportion of deaths occurred at younger ages in these regions compared with the high and high-middle SDI regions (table 3).

    Pneumoconiotic dusts

    The PAF for all pneumoconioses was assumed to be 100%. Silicosis (48%) was the largest specific cause of death from pneumoconiosis, ahead of asbestosis (16%) and CWP (12%), but about one-quarter of the deaths were classified in the ‘other pneumoconiosis’ category. There was a similar distribution between pneumoconiosis categories in terms of DALYs (table 4). The number of deaths increased with age until age 85 years and over, and the age-standardised death rates were highest in the older age groups. There was a broader distribution of DALYs across age groups, and although the rates still increased with increasing age, the rate was highest in the 75–84 year age group (online supplementary figures S9-S12).

    Table 4

    Global occupational-attributable deaths and DALYs from pneumoconioses due to exposure to pneumoconiotic dusts, by region, 2016 (number and rate)*

    The highest number of deaths and DALYs overall and for silicosis and CWP occurred in East Asia, South Asia and Western Europe, with high-income North America replacing East Asia for asbestosis deaths. The age-standardised death rates were highest in high-income Asia Pacific, East Asia and Oceania (the lowest rates were in Southeast Asia and the Caribbean), and the DALY rates highest in East Asia, Oceania and Southern sub-Saharan Africa (the lowest rates were in the Caribbean and Southeast Asia). Sixty-two percent of the silicosis deaths and 36% of the CWP deaths occurred in East Asia, and 27% of the asbestosis deaths occurred in Western Europe, which also had the second-highest rate (behind East Asia) of silicosis deaths. Western Europe, South Asia and East Asia had the highest number of asbestosis deaths, and East Asia, Australasia and Western Europe had the highest rate of asbestosis deaths (table 4—the rate data for individual pneumoconioses are not shown here).

    Changes over time

    For COPD, there was little change (4% rise) in the number of deaths due to occupational exposure to PMGF and SHS between 1990 and 2016, but the (standardised) rate of death from COPD declined by 41% over this time. For asthmagens, the number of deaths due to occupational exposure increased by 7% and the rate of death declined by 36%. The number of deaths from pneumoconioses changed minimally (1%) over this period, but the rate of death from pneumoconioses declined by 41%. Changes in the numbers and rates of DALYs were similar to those seen for deaths, except for asthma, which had a 27% increase in DALYs between 1990 and 2016. The PAFs for asthma rose considerably over this time (21% for deaths; 28% for DALYs), but there was little change in the PAFs for COPD (table 5).

    Table 5

    Change in global occupational-attributable deaths and DALYs from chronic respiratory disease due to occupational exposure to asthmagens, PMGF, SHS and pneumoconiotic dusts between 1990 and 2016, number and per capita (number and percent [95% uncertainty interval])*


    This analysis of the GBD 2016 study has shown there is a considerable burden of chronic respiratory disease worldwide and in all regions arising from exposure to occupational risk factors. Chronic obstructive pulmonary disease is the primary resulting disease, in terms of both deaths and DALYs, but asthma and pneumoconioses are also important. Rates were much higher in males than females for all these disorders, but important in both. The lower female rates reflect the fact that women are less likely to be employed in tasks that involve the relevant exposures.23 The results are consistent with those from the overall GBD respiratory analysis.30 The decreases in per capita burden for most measures, and the increase for asthma DALYs, result primarily from changes in the relevant PAFs that, in turn, reflect changes in the occupation and industry distribution, which are the basis of the exposure assessments.

    PMGF, SHS and COPD

    The global estimate of the PAF for COPD arising from occupational exposure to PMGF and SHS (17% for DALYs; 16% for deaths) is consistent with most published findings for individual countries and overall. These have typically reported PAFs of the order of 10%–15%, although much higher values have been estimated, particularly for non-smokers, typically due to differences in the level or type of exposures of the included subjects or the use of different assumptions.1 10 31–35 In addition, as smoking rates diminish, the PAF for occupational risk factors will increase. In comparison, the GBD 2016 study estimated PAFs for COPD in regard to smoking and SHS of 43% and ambient particulate matter pollution of 27%.22 The CRA study (covering the year 2000) estimated 318 000 deaths and PAFs from occupational exposure of 13% based on DALYs and 12% based on deaths.14 The Burden of Obstructive Lung Disease (BOLD) study documented a direct relationship between COPD prevalence and number of years worked in dusty jobs.36

    Asthmagens and asthma

    As with COPD related to occupational exposures, the occupational asthma PAF estimates of 10% for DALYs and 9% for deaths from this study are consistent with most published findings for individual countries, which are of the order of 10%–15%3 7 8 10 28 and comparable to the PAF due to smoking (10% for DALYs; 14% for deaths).22 The CRA study, which was based on the year 2000, estimated a PAF of 11% based on DALYs and 17% based on deaths (and estimated 38 000 deaths),14 the differences primarily arising from changes in the employment distribution and slight differences in the general methodology.22

    Pneumoconiotic dusts and pneumoconioses

    Obtaining reliable global information on pneumoconiosis cases is challenging. This analysis identified silicosis as the predominant pneumoconiosis, with much lower numbers of cases of asbestosis and CWP. The increase in rates with age is consistent with the published literature,37 and the number of deaths is consistent with the publicly available data for many countries, but also differs considerably for some others for which the estimates here are notably different from the numbers reported in the WHO Mortality Database.38 The reason for this is not clear, but presumably is because of the use of different primary data sources and assumptions in the GBD modelling process. It is likely that most of the moderate proportion of pneumoconiosis deaths and DALYs (both 23%) coded in GBD 2016 as due to ‘Other pneumoconioses’ were actually due to silicosis, asbestosis or CWP, as these have always been identified as the three main pneumoconioses. The different coding is likely to have arisen due to incomplete coding in the source data and the way this was allocated to specific categories.

    Methodological considerations and limitations

    Most of the methodological issues specific to the three main outcomes of interest have already been considered in the relevant sections of the Discussion. The main general uncertainties have been considered in detail in the companion overview paper.21 Issues of particular relevance to the presented analysis included basing exposure prevalence estimates on industry (for PMGF and SHS) and occupation (for asthmagens); uncertainty in the prevalence and level of exposure to PMGF overall and in different industries; the potential for mismatch between the relative risk estimates used and the exposure circumstances to which they have been applied; not explicitly taking into account the potential effect of differences in smoking habits and environmental exposures between regions and over time; probable heterogeneity in terms of how chronic respiratory conditions are identified, diagnosed and managed worldwide; and not including some potentially relevant risk factors and outcomes such as respiratory infections,39 other occupational causes of fibrosis apart from pneumoconioses and lung disease arising from nanoparticle exposure.40 For both COPD and asthma, the extent and effect of any mismatch between the exposure and the relative risk estimate applied in LMI countries are not clear. It would be helpful to have usable information on this from LMI countries, which might allow different risk estimates to be applied in these countries if appropriate. However, currently the necessary data are not available.

    Implications and uses of the data

    The main finding of this study is that workplace exposures resulting in COPD, asthma and pneumoconioses remain important contributors to the burden of disease in all regions of the world. The relevant exposures are respiratory and it should be possible to minimise all (or most), and in some instances to essentially eliminate them, through appropriate commitment to, and implementation of, exposure control interventions to decrease the airborne exposure levels of the relevant hazards. However, it must be recognised that there are a range of PMGF implicated as increasing the risk of COPD and hundreds of known occupational asthmagens. Elimination or appropriate control of many of these exposures will take considerable resources and effort and requires continued vigilance. The study does not provide information on the cost or practicality of eliminating or better controlling the relevant exposures, and the results for COPD and pneumoconioses largely reflect past exposures. However, the high burden of COPD cases suggests the relevant exposures should be a priority in the area of occupational airborne exposures resulting in chronic respiratory disease. The findings also have implications for healthcare costs and social protection in older individuals. Finally, further investment in country-level data sources, especially in LMI countries, would help improve the accuracy and usefulness of the estimates generated by the GBD study.


    There are many respiratory conditions that can arise directly, or indirectly, from work. The results from this study indicate that non-malignant/non-infectious respiratory diseases arising from occupational exposures are an important cause of death and disability worldwide. Many of these cases should be preventable by adopting better health and safety approaches, particularly through improved engineering and working conditions.

    Supplementary data


    View Abstract


    • Collaborators Collaborators details are as follows: Tim Driscoll, Kyle Steenland, Neil Pearce, Lesley Rushton, Sally J. Hutchings, Kurt Straif, Degu Abate, Dilaram Acharya, Anurag Agrawal, Fares Alahdab, Kefyalew Addis Alene, Sofia Androudi, Mina Anjomshoa, Carl Abelardo T. Antonio, Olatunde Aremu, Zerihun Ataro, Alaa Badawi, Joseph Adel Mattar Banoub, Suzanne Lyn Barker-Collo, Neeraj Bedi, Derrick A. Bennett, Robert Bernstein214, Mircea Beuran, Krittika Bhattacharyya, Ali Bijani, Zahid A. Butt, Juan J. Carrero, Carlos A. Castañeda-Orjuela, Odgerel Chimed-Ochir, Lalit Dandona, Rakhi Dandona, Anh Kim Dang, Ahmad Daryani, Beruk Berhanu Desalegn, Samath Dhamminda Dharmaratne, Shirin Djalalinia, Eleonora Dubljanin, Soheil Ebrahimpour, Ziad El-Khatib, Mohammad Fareed, Andre Faro, Eduarda Fernandes, Florian Fischer, Takeshi Fukumoto, Silvano Gallus, Teklu Gebrehiwo Gebremichael, Kebede Embaye Gezae, Ayman Grada, Yuming Guo, Rahul Gupta, Arvin Haj-Mirzaian, Arya Haj-Mirzaian, Samer Hamidi, Mehedi Hasan, Milad Hasankhani, Simon I. Hay, Chi Linh Hoang, Michael K. Hole, H Dean Hosgood, Mihaela Hostiuc, Sorin Hostiuc, Seyed Sina Naghibi Irvani, Sheikh Mohammed Shariful Islam, Mihajlo Jakovljevic, Ravi Prakash Jha, Jost B. Jonas, Amaha Kahsay, Amir Kasaeian, Norito Kawakami, Yousef Saleh Khader, Morteza Abdullatif Khafaie, Ejaz Ahmad Khan, Mohammad Hossein Khosravi, Jagdish Khubchandani, Yun Jin Kim, Ruth W. Kimokoti, Adnan Kisa, Manolis Kogevinas215, Soewarta Kosen, Parvaiz A. Koul, Ai Koyanagi, Barthelemy Kuate Defo, G Anil Kumar, Dharmesh Kumar Lal, Arman Latifi, James Leigh, Miriam Levi, Shanshan Li, Shai Linn, Narayan Bahadur Mahotra, Marek Majdan, Reza Malekzadeh, Mohammad Ali Mansournia, Francisco Rogerlândio Martins-Melo, Benjamin Ballard Massenburg, Varshil Mehta, Addisu Melese, Mulugeta Melku, Ziad A. Memish, Walter Mendoza, Tuomo J. Meretoja, Tomislav Mestrovic, GK Mini, Erkin M. Mirrakhimov, Babak Moazen, Naser Mohammad Gholi Mezerji, Shafiu Mohammed, Ali H Mokdad, Lorenzo Monasta, Yoshan Moodley, Mahmood Moosazadeh, Ghobad Moradi, Lidia Morawska, Shane Douglas Morrison, Seyyed Meysam Mousavi, Ghulam Mustafa, Vinay Nangia, Ionut Negoi, Ruxandra Irina Negoi, Cuong Tat Nguyen, Trang Huyen Nguyen, Molly R. Nixon, Richard Ofori-Asenso, Felix Akpojene Ogbo, Andrew T. Olagunju, Bolajoko Olubukunola Olusanya, Mahesh P A, Songhomitra Panda-Jonas, Eun-Kee Park, Sanghamitra Pati, Mostafa Qorbani, Anwar Rafay, Alireza Rafiei, Fakher Rahim, Vafa Rahimi-Movaghar, Fatemeh Rajati, Robert C. Reiner, Satar Rezaei, Leonardo Roever, Luca Ronfani, Gholamreza Roshandel, Basema Saddik, Saeid Safiri, Mohammad Ali Sahraian, Abdallah M. Samy, David C. Schwebel, Sadaf G. Sepanlou, Berrin Serdar, Masood Ali Shaikh, Aziz Sheikh, Mika Shigematsu, Rahman Shiri, Reza Shirkoohi, Si Si, João Pedro Silva, Dhirendra Narain Sinha, Moslem Soofi, Joan B. Soriano, Chandrashekhar T. Sreeramareddy, Jeffrey D. Stanaway, Mark A. Stokes, Mu'awiyyah Babale Sufiyan, Ipsita Sutradhar, Rafael Tabarés-Seisdedos, Ken Takahashi, Yonatal Mesfin Tefera, Mohamad-Hani Temsah, Marcos Roberto Tovani-Palone, Bach Xuan Tran, Khanh Bao Tran, Lorainne Tudor Car, Irfan Ullah, Pascual R. Valdez, Job F. M. van Boven, Tommi Juhani Vasankari, Francesco S. Violante, Giang Thu Vu, Gregory R. Wagner, Yasir Waheed, Yuan-Pang Wang, Biruck Desalegn Yirsaw, Naohiro Yonemoto, Chuanhua Yu, Mohammad Zamani, and Stephen S. Lim (online supplementary file)

    • Contributors The draft manuscript was prepared by Tim Driscoll, with input from Sally Hutchings, Lesley Rushton, Kyle Steenland and Kurt Straif. All listed authors have contributed appropriately to the GBD project and to the review and modification of the manuscript. The final manuscript was prepared by TD following comments from co-authors and Journal reviewers and editors.

    • Funding The overall GBD study is partly funded by the Bill & Melinda Gates Foundation. The work reported in this paper was partly supported by funding from the World Health Organization. The funders had no role in the study design, data collection, data analysis, data interpretation or writing of the report. The authors had access to the data in the study and the final responsibility to submit the paper.

    • Disclaimer The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health.

    • Competing interests CATA reports personal fees from Johnson & Johnson (Philippines), Inc., outside the submitted work. NK reports personal fees from Junpukai Foundation and Softbank, Co.; and grants from Fujitsu, LTD, Fujitsu Software Technologies, LTD and Softbank, Co., outside the submitted work. JK reports grants from Merck Pharmaceuticals, outside the submitted work. TJM reports grants from Cancer Foundation Finland sr., during the conduct of the study.

    • Patient consent for publication Not required.

    • Provenance and peer review Not commissioned; externally peer reviewed.

    • Data availability statement Data are available in a public, open access repository. The data on which this analysis is based are available on the GBD Compare web site ( Some of the raw data (where data owners give permission or where it is already public access) is available on the data section of the IHME GBD web site (

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Linked Articles