key: cord-0964767-rt2flfoo
authors: Liu, Jinli; Bai, Ruhai; Chai, Zhonglin; Cooper, Mark E.; Zimmet, Paul Z.; Zhang, Lei
title: Low- and middle-income countries demonstrate rapid growth of type 2 diabetes: an analysis based on Global Burden of Disease 1990–2019 data
date: 2022-05-19
journal: Diabetologia
DOI: 10.1007/s00125-022-05713-6
sha: 29e0f66ded023f50d0797e70758c47947f464105
doc_id: 964767
cord_uid: rt2flfoo

AIMS/HYPOTHESIS: The study aims to quantify the global trend of the disease burden of type 2 diabetes caused by various risks factors by country income tiers. METHODS: Data on type 2 diabetes, including mortality and disability-adjusted life years (DALYs) during 1990–2019, were obtained from the Global Burden of Disease Study 2019. We analysed mortality and DALY rates and the population attributable fraction (PAF) in various risk factors of type 2 diabetes by country income tiers. RESULTS: Globally, the age-standardised death rate (ASDR) attributable to type 2 diabetes increased from 16.7 (15.7, 17.5)/100,000 person-years in 1990 to 18.5 (17.2, 19.7)/100,000 person-years in 2019. Similarly, age-standardised DALY rates increased from 628.3 (537.2, 730.9)/100,000 person-years to 801.5 (670.6, 954.4)/100,000 person-years during 1990–2019. Lower-middle-income countries reported the largest increase in the average annual growth of ASDR (1.3%) and an age-standardised DALY rate (1.6%) of type 2 diabetes. The key PAF attributing to type 2 diabetes deaths/DALYs was high BMI in countries of all income tiers. With the exception of BMI, while in low- and lower-middle-income countries, risk factors attributable to type 2 diabetes-related deaths and DALYs are mostly environment-related, the risk factors in high-income countries are mostly lifestyle-related. CONCLUSIONS/INTERPRETATION: Type 2 diabetes disease burden increased globally, but low- and middle-income countries showed the highest growth rate. A high BMI level remained the key contributing factor in all income tiers, but environmental and lifestyle-related factors contributed differently across income tiers. DATA AVAILABILITY: To download the data used in these analyses, please visit the Global Health Data Exchange at http://ghdx.healthdata.org/gbd-2019. GRAPHICAL ABSTRACT: [Image: see text] SUPPLEMENTARY INFORMATION: The online version of this article (10.1007/s00125-022-05713-6) contains peer-reviewed but unedited supplementary material.

Diabetes mellitus is a fast-growing ongoing global health emergency of the twenty-first century [1, 2] . Studies show that diabetes alone accounted for 3.96 million deaths worldwide in 2010 [3] , and this number increased to 4.20 million among adults in 2019 [1] . Diabetes has emerged as the fourth leading cause of disability globally, and the number of disabilityadjusted life years (DALYs) caused by diabetes was 66.3 million globally in 2019 [4] . Type 2 diabetes is the most common type of diabetes and accounts for around 90% of all diabetic cases worldwide [1, 2, 5] . The economic development and improvement in healthcare in individual countries have a considerable impact on the disease burden of type 2 diabetes. In many developed countries, including Scotland [6] , Canada [7] , UK [8] , Denmark [9] , Sweden [10] , Australia [11] and USA [12] , type 2 diabetes-related mortality has declined steadily. People living with type 2 diabetes are at higher risk of chronic diseases, viral hepatitis and even COVID-19 severity [13] [14] [15] [16] . However, large developing economies such as China reported an ongoing and increasing trend of diabetes mellitus, with diabetes mortality increasing from 5.3 deaths/100,000 people in 1990 to 10.9 deaths/ 100,000 people in 2017 [17] . In developing countries with limited resources, inaccessibility to type 2 diabetes medications and other treatments [18] is an important cause of type 2 diabetes-related mortality and disability. Many resourcelimited countries employed intensive blood glucose control, which unexpectedly increased the risk of hypoglycaemia and potential death in its population [19] . An epidemiological study suggested that during 2010-2030, the number of adults with diabetes mellitus in developing countries would increase by 69%, which is more than threefold the 20% predicted increase in developed countries [20] .

The Global Burden of Disease (GBD) study indicated that high BMI was the leading risk factor contributing to 41.2% type 2 diabetes mortality and 54.0% DALYs, followed by ambient particulate matter pollution risk factors (accounting for 18.0% and 18.3% of type 2 diabetes-related deaths and DALYs), and diet low in whole grains (accounting for 16.5% and 21.2% of type 2 diabetes-related deaths and DALYs) in 2017 [21] followed. However, these previous studies did not address the impact of diabetes on mortality and disability attributable to major risk factors in countries at different income levels. Understanding this knowledge will enable policymakers to tailor their policies to adjust for the various levels of economic development in these countries. This will provide informed strategies for coordinated actions to improve their health and social security systems to mitigate diabetes-related mortality and disability. In this study, we aimed to estimate type 2 diabetes-related deaths and DALYs attributable to various modifiable risk factors over the period of 1990-2019 in countries at various income tiers.

Overview The GBD 2019 estimation of attributable burden followed the general framework established for comparative risk assessment [22, 23] used in GBD since 2002 [24] . GBD 2019 estimated prevalence of exposure and attributable deaths, years of life lost (YLLs), years lived with disability (YLDs), and DALYs for 23 age groups; men, women, and both sexes combined; 87 behavioural, environmental and occupational, and metabolic risk factors; reported estimates for 369 diseases and injuries; and 203 countries (Taiwan: province of China) and territories that were grouped into 21 regions and seven super-regions [4, 25] . All available data on causes of death are standardised and pooled into a single database used to generate cause-specific mortality estimates by age, sex, year and geography [4] . In GBD 2019, DALYs were computed by adding YLLs and YLDs for each cause, location, age group, sex and year.

The attributable burden is the potential reduction in the current disease burden if past population exposure had shifted to an alternative or counterfactual distribution of risk exposure. In this study, we included 13 risk factors introduced in the GBD studies. These attributable risk factors include metabolic (i.e. high BMI), environmental and occupational (i.e. ambient particulate matter pollution and household air pollution from solid fuels), behavioural (i.e. tobacco [smoking and second-hand smoking], dietary [diet low in whole grains, diet low in fruits, diet low in nuts and seeds, diet high in sugarsweetened beverages, diet high in red meat, diet low in fibre and diet high in processed meat], and low physical activity). The definition of these risk factors has been provided in detail in previous publications [26] . The GBD 2019 study adopted a set of unified RRs for all attribution analysis globally [25] . In brief, the RRs were modelled using meta-regression with pooled data from prospective cohort studies or published literature reviews. For each risk-outcome pair, we used the expected summary exposure value (SEV) to calculate the expected population attributable fractions (PAFs) [4, 25] . The method of evaluating the impact of environmental pollutants has been described in a GBD study [25] .

Socioeconomic status was defined based on the gross national per capita income, as classified by the World Bank, i.e. low-income countries (LICs), lower-middle-income countries (LMICs), upper-middle-income countries (UMICs), and high-income countries (HICs) [27] .

Disease burden attributable to risk factors For continuous risk factors, such as BMI, the formula for the PAF is defined as:

Where RR(x) is the RR of a certain disease for exposure level x, P 1 (x) is the population distribution of the exposure, P 2 (x) is the minimum theoretical exposure distribution, and m is the maximum exposure level. The minimum theoretical exposure is the counterfactual condition of exposure and has been previously defined for each risk factor [26] . For categorical risk factors (e.g. smoking and physical inactivity), the formula for PAF is:

where i is exposure level, RR i is RR for exposure level i, and P i is the prevalence of exposure level i. For calculating the burden of multiple risk factors, validation studies have reported congruency between the true risk associated with multiple risk factors affecting the same outcome and a multiplicative aggregation of the PAFs of the individual risk factors [28] .

We calculated the joint PAF of multiple risk factors based on the following formula. The mediation factors were estimated using GBD 2019 [25] .

where J is the number of risk factors for calculating the joint effect, PAF ioast is the attributable fraction of i risk factor, MF jio is the mediation factor between risk factor i and a certain disease o through risk factor j, a was the age group, s was sex and t was year.

Indicators of annual change in the burden of disease The rate of annual growth was measured as the average annual change in type 2 diabetes burden from 1990 to 2019:

All statistical analyses were performed using the R v3.5.1 (https://www.r-project.org/). A probability value of p<0.05 was considered statistically significant.

Global type 2 diabetes over time Globally, the number of deaths attributable to type 2 diabetes increased from 0.61 (95% CI 0.57, 0.64) million in 1990 to 1.47 (1.37, 1.57) million in 2019, with an annual growth rate of 4.9% (4.6, 5.2). The age-standardised risk-attributable death rates increased from 16.7 (15.7, 17.5)/100,000 person-years to 18 (Fig. 2a,b) . The annual growth rates for ASDR and age-standardised DALY rate showed no significant differences across all three time periods in all income tiers (Fig. 2a,b) .

PAFs for type 2 diabetes Globally, the leading risk factor for type 2 diabetes-related ASDR was high BMI. The ASDR attributable to high BMI had increased from 5.0 (95% CI 3.0, 7.3)/100,000 person-years to 7.6 (5.3, 10.0)/100,000 person-years during 1990-2019. This was followed by ambient particulate matter pollution (increased from 1.6 [1.0, 2.2] to 2.5 [1.7, 3.2]/100,000 person-years) and low physical activity (increased from 1.5 [0.8, 2.5] to 1.6 [0.8, 2.7]/ 100,000 person-years) ( Table 1 ). The three risk factors collectively contributed to over half (54.1%) of all type 2 diabetes-related ASDR in 2019 globally (Fig. 3a) (Fig. 3b) . With the exception of smoking and second-hand smoke, the proportions of deaths and DALYs attributable to type 2 diabetes varied less by sex globally in 2019 (Fig. 3a, b) .

In LICs, the leading risk factor for type 2 diabetes-related ASDR was high BMI (ASDR attributable to high BMI increased from 6.2 [95% CI 2.7, 10.9]/100,000 person-years to 9.6 [5.7, 14.2]/100,000 person-years during 1990-2019). This was followed by household air pollution from solid fuels (reduced from 9.4 [5.7, 18 .5] to 7.1 [4.7, 11 .2]/100,000 person-years) and low fruits in the diet (stabilised at 2.6/ 100,000 person-years during 1990-2019) (Table 1, Fig. 3a) . The leading risk factors contributing to age-standardised DALY rates were similar (Table 1, Fig. 3b ).

In LMICs, high BMI remained the leading risk factor for type 2 diabetes-related deaths (ASDR attributable to high BMI increased from 3.9 [95% CI 1.9, 6.6] to 10 Fig. 3a) . The leading risk factors contributing to age-standardised DALY rates were similar (Table 1, Fig. 3b ).

In UMICs, high BMI remained the leading risk factor attributable to type 2 diabetes-related deaths (ASDR attributable to high BMI increased from 4.9 [95% CI 3.1, 7.1] to 6.8 [4.7, 9. 1]/ 100,000 person-years during 1990-2019). This was followed by ambient particulate matter pollution (increased from 1.5 [1.0, 2.1] to 2.5 [1.8, 3.2]/100,000 person-years) and low physical activity (increased from 1.4 [0.7, 2.2] to 1.5 [0.8, 2.5]/ 100,000 person-years) ( Table 1 , Fig. 3a ) However, we observed a slightly different contribution pattern in agestandardised DALY rates, with high BMI, ambient particulate matter pollution and smoking being the leading risk factors (Table 1, Fig. 3b ).

In HICs, despite high BMI remaining the leading risk factor for type 2 diabetes-related ASDR, its contribution was substantially lower than other income tiers and decreasing over time (ASDR attributable to high BMI reduced from 5.3 [95% CI 3.4, 7.2]/100,000 person-years to 4.7 [3.4, 6.2]/100,000 personyears during 1990-2019). This was followed by diet high in processed meat (reduced from 1.6 [1.1, 1.9] to 1.2 [0.9, 1.5]/ 100,000 person-years) and low physical activity (reduced from 1.4 [0.7, 2.3] to 1.0 [0.5, 1.7]/100,000 person-years) ( Table 1 , Fig. 3a) . For age-standardised DALY rates, high BMI, diet high in processed meat and smoking are leading risk factors, and all three factors demonstrated flattening the upward tendency over time (Table 1, Fig. 3b ).

Our study analysed the correlation between per capita gross domestic product (GDP) and the attributable burden of risk factors in 203 countries in 2019 (electronic supplementary material [ESM] Table 1 ). Per capita GDP was positively correlated with PAFs for death and DALYs caused by diet high in red meat, diet high in processed meat, diet high in sugar-sweetened beverages and low physical activity (all p<0.0001), but negatively correlated with PAFs for death and DALYs caused by household air pollution from solid fuels, smoking and diet low in fruits (all p<0.0001). However, per capita GDP showed a negative correlation with PAFs only for DALYs caused by second-hand smoke (p<0.0001).

Our study showed that type 2 diabetes-related mortality and DALYs varied substantially during 1990-2019, globally and by income tiers. The numbers of deaths and DALYs attributable to type 2 diabetes doubled from its 1990 levels to 1.47 million and 66.3 million in 2019. ASDR and age-standardised DALY rate also increased by 10.8% and 27.6% to 18.5/100,000 and 801.5/100,000, respectively, during the same period. Overall, type 2 diabetes-related ASDR and age-standardised DALY rate increased with the declining income tier of the countries. Both rates in LICs and LMICs were 2-3 times higher than those in UMICs and HICs in 2019. Further, LMICs reported the largest increase in the average annual growth of ASDR (1.3%) and an age-standardised DALY rate (1.6%) of type 2 diabetes. In contrast, the annual growth of ASDR and age-standardised DALY rate was the lowest (−0.8% and nearly zero, respectively) in HICs, during 1990-2019. BMI is the key common leading risk factor of type 2 diabetes disease burden across all income tiers. With the exception of BMI, while in low-and middle-(lower-middle-and upper-middle-) income countries, risk factors attributable to type 2 diabetes-related deaths and DALYs are mostly environment-related, the risk factors in HICs are mostly lifestyle-related.

Our study demonstrated a rapidly growing trend of type 2 diabetes burden in low-and middle-income country settings, although concerns related to type 2 diabetes burden have historically focused on HICs, which was consistent with previous studies [29, 30] . Both ASDR and age-standardised DALY rate in LICs and LMICs are consistently 2-3 times higher than those in UMICs and HICs throughout our study duration. In particular, the fact that HICs reported a negative or nearly zero annual growth in both age-standardised rates whereas LMICs have more significant positive annual growth suggests that the patients with diabetes may result in significant delays in diagnosis and treatment of diabetes, leading to a severe burden of diabetes disability and related complications. Third, PAFs for death and DALYs differed by income tiers, as indicated by our study. Exposure to environmental risk factors for diabetes is more common in lower income countries as a result of limited options of energy consumption in these countries. If the current trend continues, type 2 diabetes will likely cause an increasingly severe disease burden in lower income countries than their higher income counterparts. Other studies have also demonstrated similar associations between income tiers and type 2 diabetes-related mortality and DALYs. Lin et al [31] reported that, worldwide, the rates of mortality and disability caused by type 2 diabetes exhibited an upward trend in low-income and lower-middle-income countries. In comparison, the mortality rate due to type 2 diabetes in HICs, such as Australia [11] and Sweden [10] , have shown a declining trend. Consistently, a separate study reported a large reduction in complications related to type 2 diabetes in HICs during 2000-2015, indicating an alleviation of disability caused by type 2 diabetes in these countries [32] . Lin et al [31] reported that for type 2 diabetes, the association between mortality or DALY rates with Sociodemographic Index (SDI) demonstrated an inverse U-shaped curve with the higher rates occurring in low-middle, middle, and highmiddle SDI countries. The corresponding rates in low-SDI and high-SDI countries were substantially lower.

Our study indicates that BMI remains the leading risk factor for type 2 diabetes-related deaths and DALYs globally and across all income tiers. Obesity is a well-documented important public health issue and a key contributor to numerous chronic diseases [33, 34] . Historically, the increase in the prevalence of obesity began in HICs in the 1970s, followed by most middle-income countries, and more recently, lowincome countries [35] . The increase in obesity prevalence is likely a consequence of economic development and affects the occurrence of type 2 diabetes. However, economic development and wealth can enable better type 2 diabetes intervention. Over the past decade, lifestyle interventions in HICs, aiming to reduce the risk of overweight and obesity, have reportedly reduced the risk of type 2 diabetes, as demonstrated by the US Diabetes Prevention Program [36] and the Dutch Diabetes Prevention Study [37] . The American Cancer Society's Cancer Prevention Study also reported that intentional weight loss was associated with a 28% reduction in type 2 diabetes-related mortality and its complications (cardiovascular disease) [38] . In contrast, prevention programmes for overweight and obesity in low-and middle-income countries have commenced later than high-income countries. In 2016, the Chinese government declared the 'Healthy China 2030' initiative, aiming to facilitate appropriate diet and physical activities to reduce obesity and hence type 2 diabetes in the Chinese population. A steady growing trend of new type 2 diabetes has been documented during 1990-2017 [39] , which is closely related to the rise in prevalence of overweight and obesity in China in recent years [40] . The country with the largest number of adults with diabetes aged 20-79 years in 2021 is China, and it accounts for 26.2% of the number of people living with diabetes globally [41] .

Our study shows that both household and ambient air pollution are major risk factors for type 2 diabetes deaths and DALYs in LICs and LMICs, consistent with previous findings [42] . Pollution and poverty are closely related. An estimated 3 billion people in low-and middle-income countries, mostly in rural communities, still use solid fuels (firewood, biomass or charcoal) and traditional stoves for heating and cooking [43] . This results in large populations being exposed to household air pollution and subsequent type 2 diabetes onset and complications because air pollution is a leading cause of insulin resistance [44] . Ambient particulate matter pollution, often measured by the air density of particulate matter 2.5 with an aerodynamic diameter less than 2.5 μm (PM 2.5 ), has led to 4.2 million deaths and 103.1 million DALYs in 2015 [45] . Ambient particulate matter pollution disproportionately affects the poor and the vulnerable in low-income and middle-income countries, which account for 79% of adults living with diabetes [1, 42] . The population is often neglected in these settings where economic development is prioritised, and its health impacts on people living with diabetes are largely underestimated. Our study demonstrates significant contributions of both household and ambient air pollution to type 2 diabetes and promotes sustainable green development and pollutant reduction.

Our study also demonstrates that diet rich in processed meat is the second most important risk factor for type 2 diabetes-related deaths and DALYs in HICs. Meta-analysis studies have consistently shown that diet rich in processed meat is associated with an increased risk of type 2 diabetes in HICs [46, 47] . Over the past decade, despite some HICs having witnessed declining trends in the consumption of red and processed meat, consumption in these countries remains higher than in other income countries [48] . McMichael et al found that the global average consumption was 100 g per person per day, whereas the consumption now is as high as 200-250 g in HICs [49] . However, as global meat prices have fallen, red and processed meat has become increasingly available in low-and middle-income countries [48] . At the same time, health-motivated tax on red and processed meat is low in these countries [50] . This could potentially increase consumption of processed meat and hence the disease burden of type 2 diabetes in low-and middle-income countries.

Our study has several limitations. First, the projection relied heavily on estimates driven by GBD 2019; therefore, similar limitations in estimates of deaths, DALYs and attributable burden in the GBD study also apply to this study. Second, it deviates somewhat from the comparative risk assessment approach when the GBD study estimates the burden for total particulate matter pollution and divides this total burden proportionately between the ambient and household particulates. Third, for exposure measurement, patterns of data availability are non-uniform across geography and over time. Although the GBD study has modified this, it still causes a bias in comparison.

Conclusions Type 2 diabetes disease burden has increased globally in countries of all income tiers. High BMI, ambient particulate matter pollution and low physical activity collectively contributed to half of type 2 diabetes-related deaths and two-thirds of type 2 diabetes-related DALYs globally in 2019. BMI remains the most important contributing factor to the type 2 diabetes disease burden. Following BMI, the environmental factors are the next important contributor to the type 2 diabetes burden in low-and lower-middle-income countries, whereas dietary factors are the next important contributing factor in HICs.

International Diabetes Federation

Diabetes: a 21st century challenge

Mortality attributable to diabetes: estimates for the year 2010

Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: a systematic analysis for the global burden of disease study 2019

Diabetes mellitus statistics on prevalence and mortality: facts and fallacies

Trends in type 2 diabetes incidence and mortality in Scotland between

Trends in diabetes prevalence, incidence, and mortality in Ontario, Canada 1995-2005: a populationbased study

Mortality trends in patients with and without diabetes in Ontario, Canada and the UK from 1996 to 2009: a population-based study

Improved survival among patients with complicated type 2 diabetes in Denmark: a prospective study

Mortality and cardiovascular disease in type 1 and type 2 diabetes

Mortality trends among people with type 1 and type 2 diabetes in Australia: 1997-2010

Trends in death rates among U.S. adults with and without diabetes between 1997 and 2006: findings from the National Health Interview Survey

Incidence of medication-treated depression and anxiety associated with long-term cancer, cardiovascular disease, diabetes and osteoarthritis in community-dwelling women and men

Comparison of type 2 diabetes mellitus incidence in different phases of hepatitis B virus infection: a meta-analysis

Elevated glucose level leads to rapid COVID-19 progression and high fatality

High fasting blood glucose level with unknown prior history of diabetes is associated with high risk of severe adverse COVID-19 outcome

Mortality, morbidity, and risk factors in China and its provinces, 1990-2017: a systematic analysis for the global burden of disease study 2017

Differences in the availability of medicines for chronic and acute conditions in the public and private sectors of developing countries

Global burden of hypoglycaemia-related mortality in 109 countries, from 2000 to 2014: an analysis of death certificates

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

Global mortality, disability, and the contribution of risk factors: global burden of disease study

Comparative quantification of health risks conceptual framework and methodological issues

Measuring universal health coverage based on an index of effective coverage of health services in 204 countries and territories, 1990-2019: a systematic analysis for the global burden of disease study 2019

a systematic analysis for the global burden of disease study 2019

A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990-2010: a systematic analysis for the global burden of disease study 2010

World Bank Country and Lending Groups (2019) The World Bank

Validation of a new predictive risk model: measuring the impact of the major modifiable risks of death for patients and populations

The burden and costs of chronic diseases in low-income and middle-income countries

Availability and affordability of essential medicines for diabetes across highincome, middle-income, and low-income countries: a prospective epidemiological study

Global, regional, and national burden and trend of diabetes in 195 countries and territories: an analysis from 1990 to 2025

The changing face of diabetes complications

Metabolic mediators of the effects of body-mass index, overweight, and obesity on coronary heart disease and stroke: a pooled analysis of 97 prospective cohorts with 1·8 million participants

The global obesity pandemic: shaped by global drivers and local environments

Long-term outcomes of lifestyle intervention to prevent diabetes in American Indian and Alaska native communities: the special diabetes program for Indians diabetes prevention program

High diabetes distress among ethnic minorities is not explained by metabolic, cardiovascular, or lifestyle factors: findings from the Dutch diabetes pearl cohort

Intentional weight loss and mortality among overweight individuals with diabetes

Trends in the incidence of diabetes mellitus: results from the global burden of disease study 2017 and implications for diabetes mellitus prevention

Marriage contributes to higher obesity risk in China: findings from the China health and nutrition survey

International Diabetes Federation

The lancet commission on pollution and health

Fires, fuel, and the fate of 3 billion: the state of the energy impoverished

Effect of environmental air pollution on type 2 diabetes mellitus

Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the global burden of diseases study 2015

Red meat consumption and risk of type 2 diabetes: 3 cohorts of US adults and an updated meta-analysis

Red and processed meat consumption and risk for all-cause mortality and cardiometabolic outcomes: a systematic review and meta-analysis of cohort studies

Socioeconomic and demographic drivers of red and processed meat consumption: implications for health and environmental sustainability

Food, livestock production, energy, climate change, and health

Healthmotivated taxes on red and processed meat: a modelling study on optimal tax levels and associated health impacts

Global Burden of Disease Collaborative Network

United States of America: Institute for Health Metrics and Evaluation (IHME)

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

Acknowledgements This manuscript uses publicly available results from the Global Burden of Disease (GBD) Results Tool. The Institute for Health Metrics and Evaluation (IHME), the University of Washington (UW) and the GBD Collaborator Network were not involved in the preparation of this manuscript; the contents and views in this manuscript are those of the authors and should not be construed to represent the views or interpretation of results of IHME, UW or the GBD Collaborator Network.Data availability To download the data used in these analyses, please visit the Global Health Data Exchange at http://ghdx.healthdata.org/ gbd-2019 [51] .Funding Open Access funding enabled and organized by CAUL and its Member Institutions. This study was supported by the National Natural Science Foundation of China (Grant number: 8191101420), Outstanding Young Scholars Funding (Grant number: 3111500001), Xi'an Jiaotong University Basic Research and Profession Grant (Grant number: xtr022019003, xzy032020032) and Xi'an Jiaotong University Young Talent Support Grant (Grant number: YX6J004). The funders had no role in study design, data collection, data analysis, interpretation or writing of the manuscript.Contribution statement LZ and JL substantially contributed by developing the conceptual framework and design of the study. JL wrote the first draft of the manuscript and performed the statistical analysis. All authors contributed to the interpretation of the results and writing. LZ and PZ critically revised the manuscript for important intellectual content. All authors have approved the final version to be published. LZ is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/ .

Authors' relationships and activities The authors declare that there are no relationships or activities that might bias, or be perceived to bias, their work.