- 1Department of Gastroenterology, The First Affiliated Hospital of Shandong First Medical University & Shandong Provincial Qianfoshan Hospital, Jinan, Shandong, China
- 2Primary Care Medicine Department, Second Affiliated Hospital of Anhui Medical University, Hefei, Anhui, China
- 3Department of Infectious Disease, Second Affiliated Hospital of Anhui Medical University, Hefei, Anhui, China
- 4Shandong University of Traditional Chinese Medicine, Jinan, Shandong, China
Background: Studies on global epidemiological patterns regarding the burden of pancreatic cancer (PC) attributable to smoking are limited. This study aimed to analyze the latest disease burden of PC attributable to smoking between 1990 and 2021, further analyze the main factors, and predict the trend in the coming period.
Methods: Data from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2021 was extracted and analyzed by different geographical levels, age, sex, and socio-demographic index (SDI). Key measures included age-standardized disability-adjusted life years (DALYs) rates (ASDR), age-standardized mortality rates (ASMR), and average annual percentage changes (AAPCs). Further analyses were conducted using the decomposition analysis and Bayesian Age-Period-Cohort (BAPC) model.
Results: Globally, the ASDR and ASMR of the disease demonstrated a decreasing trend. The burden was significantly higher among males, the middle-aged, and the elderly population. A positive correlation with SDI levels across most regions was exhibited. Join-point analysis indicated a decreasing trend in disease burden among high SDI regions while an increasing trend among low-middle SDI regions. Decomposition analysis indicated that proactive epidemiological changes had played a positive role in reducing the burden in High SDI regions. Projection analysis estimated that the disease burden in East Asia, showing a significant upward trend, warranted particular focus.
Conclusions: Despite ongoing tobacco control efforts, the global burden of the disease, which has declined only slightly, remains a significant concern, particularly in high-income areas and among men. Based on stronger tobacco control measures, increased emphasis on disease prevention, early screening, and treatment is essential.
Highlights
● This study comprehensively documents the global burden of PC attributable to smoking using the latest GBD 2021 database as well as projects the trend in the coming period.
● Despite the global decline in smoking rates, the burden of PC attributable to smoking has not decreased significantly.
● Greater consideration should be given to high-income countries and low- and middle-income countries, especially male populations, in the design and enhancement of global planning for tobacco control and health promotion.
Introduction
Pancreatic cancer (PC) represents a group of highly aggressive malignant tumors of the digestive system originating from pancreatic duct epithelial cells and acinar cells (1). The disease typically presents with an insidious onset, making early diagnosis difficult. Patients often lack opportunities for curative surgery and show poor tolerance to conventional cancer treatments, resulting in an abysmal prognosis. s (2),. According to data from the Surveillance, Epidemiology, and End Results project database of the National Cancer Institute’s data for PC in both sexes and all races, during 2014-2020, more than half (51%) of all cases were diagnosed at the distant stage with a 5-year survival rate of 3.1%, and even taking all the stages into account, the 5-year survival ratesare only 12.8% (3). According to the latest estimate based on global cancer statistics by world region for the year 2022 from the International Agency for Research on Cancer, this disease ranked as the sixth leading cause of cancer mortality, accounting for nearly 5% of all cancer deaths globally (4). The GBD 2021 reveals that the number of incident and death cases of PC in both sexes increased more than one time from 207905 and 211613 to 508532 and 505752 (5). All these reports have indicated a gradual increase in the number of incident cases and deaths from PC, which would not only pose a serious threat to human life and health but also constitute a significant burden on public health.
According to the most recent overview of the established modifiable risk factors for PC, cigarette smoking, obesity, alcohol use, new-onset diabetes, long-standing diabetes, family history of PC, and pancreatitis were the well-established risks factors for PC, with cigarette smoking remaining a leading cause worldwide (6). A uniquely large collaborative pooled analysis analyzing data from the International PC Case–Control Consortium, confirmed that compared with never smokers, the OR was 1.2 (95% confidence interval [CI] 1.0 to 1.3) for former smokers and 2.2 (95% CI 1.7 to 2.8) for current cigarette smokers, with a significant increasing trend in risk with increasing number of cigarettes among current smokers (odds ratios [OR] = 3.4 for ≥ 35 cigarettes per day, P for trend <0.0001) and duration of cigarette smoking (OR = 2.4 for ≥ 40years of smoking) (7). In addition to being a high-risk factor for the incidence, cigarette smoking has also been associated with a reduction in survival among patients with PC. Through survival analysis, Chen Yuan et al. observed that compared to never-smokers, the multivariable-adjusted hazard ratio (HR) for death in current smokers was 1.37 (95% CI, 1.11 to 1.69, P = 0.003). Additionally, as the number of pack-years of smoking increased, the survival rate showed a statistically significant negative correlation (P trend =0.008) a (8). In addition, tobacco exposure has been shown through basic experimental studies to initiate and accelerate the onset and progression of PC at the mechanistic level (9). In addition, recent basic research uncovered that cigarette smoke extract (CSE) exposure can not only significantly promote PC attributable to smoking progression by upregulating of Chromobox protein homolog 3 (CBX3) via promoting the overexpression of Y-box-binding protein 1 (YBX1) but also the invasion and metastasis of the cancer cells by inducing histone deacetylase 4 (HDAC4) (10, 11).
Through a 15-year period joint effect of the WHO Framework Convention on Tobacco Control (WHO FCTC) together with its technical package–MPOWER has reduced the global prevalence of smoking from 22.8% in 2007 to 17.0% in 2021, more than 8 million people still die each year from tobacco-related diseases (12). Smoking prevalence varies widely by region, age group, and sex, potentially affecting the distribution of PC attributable to smoking. However, studies on the global epidemiological patterns regarding the burden of PC attributable to smoking are relatively limited. In our study, for the first time, we used the latest data from the publicly available Global Burden of Disease database 2021 to analyze the distribution and changes in the global burden of PC attributable to smoking between 1990 and 2021 (13). Unlike previous studies that primarily focused on specific regions or diseases, our analysis not only covers the global burden of PC but also takes into account differences across regions, age groups, genders, and SDI levels. Additionally, this study delves into the impact of changes in age structure, population growth, and epidemiological trends on disease burden in different regions, and provides predictions for future trends. This offers new perspectives and scientific evidence for developing more targeted and cost-effective tobacco control policies.
Method
Data source
Estimating years lived with disability (YLDs), years of life lost (YLLs), disability-adjusted life-years (DALYs), and healthy life expectancy (HALE) for 371 diseases and injuries in seven super-regions, 21 regions, 204 countries and territories (including 21 countries with subnational locations), and 811 subnational locations, from 1990 to 2021, the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2021 aims to provide comprehensive and current global, regional, and national data on the burden of disease, injury, and risk factors to quantify the state of global health. These data were stratified by geographical location, socio-economic development level, gender, and different age groups (5). Based on 2021 Socio-demographic Index (SDI) values, a composite indicator of background social and economic conditions that influence health outcomes in each location, countries were divided into five quintiles, termed low, low-middle, middle, high-middle, and high. The SDI values for each country and their corresponding classifications can be accessed from the Institute for Health Metrics and Evaluation (IHME) website [http://ghdx.healthdata.org]. Pancreatic cancer has a poor prognosis and a short survival period for patients. The mortality can accurately reflect the contribution of smoking to the mortality burden of pancreatic cancer. Meanwhile, the DALYs can comprehensively consider both mortality and disease on health life expectancy, providing a more holistic quantification of the overall health loss caused by smoking. We extracted raw data from the GBD, including estimates of DALYs and deaths in absolute numbers, as well as age-standardized rates (ASRs) to conduct analysis (14). The data for this study were sourced from a publicly available database, requiring no ethical approval or informed consent. No ethical approval and informed consent were required because of the public availability of GBD and no identifiable information was included in the analyses.
Definition of PC and smoking
GBD 2021 provided age-sex-location and yearly estimates of PC prevalence, which were derived from pooled analyses across only one domain: the prevalence of PC. The smoking population in this study refers to individuals who use any tobacco product including all smoked tobacco products—eg. cigarettes, pipes, cigars, shisha, bidis, kreteks, and other locally smoked tobacco products daily or occasionally, including those who have used any smoking product in the past. While risks from smokeless tobacco, electronic cigarettes (also known as e-cigarettes), vaping products, heated tobacco products, chewing tobacco, and second-hand smoke are outside the scope of this study. For current smokers, exposure levels were estimated using two continuous measures: the number of cigarettes smoked per day and the cumulative number of cigarettes smoked per year. For former smokers, exposure was assessed based on the number of years since cessation (15).
Disease burden indicators
YLLs were the product of the number of deaths and standard life expectancy at each age of death, and YLDs were the product of the prevalence of each sequela and its corresponding disability weight. DALYs were calculated by summing YLDs and YLLs. The mortality numbers and the mortality rates are also included in the analysis. We also reported estimates generated as age-standardized results (e.g. ASDR and ASMR) to allow comparison of estimates made for populations with different age structures. Due to the potential heterogeneity in PC attributable to smoking, and burden resulting from variations in population age structures, we employed ASR and AAPC to quantify the trends in DALYs and mortality rates associated with PC. AAPC refers to the average annual percentage change across multiple periods, which is used to describe the overall trend over a longer period. It is calculated by fitting a regression line to the natural logarithm of the rates and provides a summary measure of a trend over time. We revealed the large differences and trends in smoking-attributable PC’s DALYs and deaths between sexes, ages, regions, countries/territories, and SDIs. 95% uncertainty intervals (UIs) were generated for all final estimates as the 25th and 95th percentile values of 500 draws. The AAPC, along with its 95% confidence intervals (CIs), was determined using a linear regression model to quantify changes over time.
Disease burden analysis
For descriptive Analysis, the first step in a more in-depth analysis of the study, the study described the distribution of the disease burden at different geographical levels (including global, the 21 geographic super-regions, and 204 countries and territories) based on the burden of disease indicators described above and created a global map of disease burden distribution, highlighting the regions with a high or a clear upward trend of disease burden. Additionally, stratifying age into five-year intervals, pyramid charts were drawn to show the differences in DALY number, DALY rate, mortality, and mortality rate between males and females across different age groups in each year. Furthermore, the distribution of disease burden across different SDI regions is also described. Based on the differences in disease burden across different SDI regions, Spearman’s rank-order correlation analysis was used to measure the intensity and direction of the correlations between the ASDR, ASMR, their AAPCs, and the SDI, and a two-tailed p-value was calculated. A p-value less than 0.05 was regarded as significant (16).
Join-point analysis
Join-pint regression model is first proposed by Kim in 2000 (17). Based on the temporal characteristics of the disease distribution, this model builds a piecewise regression and identifies significant changes in trends, referred to as joint points. Model results can be summarized using the metrics of annual percentage change (APC) and AAPC and their 95%CI to evaluate the trend of independent intervals of piecewise functions and assess the average trend over the entire study interval. Trends were described as increasing or decreasing when the AAPC or APC was statistically significant according to a two-sided p-value < 0.05, and as stable otherwise.
Decomposition analysis
We employed the decomposition method invented by Das Gupta (18) to quantify the impact of age structure, population growth, and epidemiological changes on the DALYs and mortality number of PC attributable to smoking. Age structure refers to the distribution of different age groups within a population. Population growth refers to changes in the total number of people in a region or country. Epidemiological changes refer to shifts or trends in the distribution and determinants of health-related states in populations, including the dynamics of disease occurrence, prevalence, and mortality, which can be influenced by various factors, including but not limited to environmental factors, Behavioral Factors, Socioeconomic Status, Demographic Shifts, Genetic Factors, Healthcare Infrastructure and Policies (19). This methodological approach quantifies the impact of each factor on total change thereby informing targeted interventions and effective public health strategies to reduce the disease burden.
Predictive analytics
To better forecast the future trends of DALYs and mortality rate, We utilized the Bayesian Age-Period-Cohort (BAPC) model, the probabilistic forecasts obtained by which are well calibrated and not too wide, demonstrating better predictive performance than other methods. The model, trained on historical data up to 2021, capturing the underlying patterns and dynamics of disease burden indicators over time, enables us to anticipate the evolving epidemiological landscape of PC.
Statistical analysis
Variables were expressed in numbers, percentages, and ratios. Pearson correlation analyses were conducted to evaluate the associations between ASR and SDI in 2021, AAPC, and the mean SDI of 204 countries and territories. All statistical analyses and visualizations for this study were carried out using R 4.1.2 software. A significance level of P < 0.05 was used to determine statistical significance.
Result
Global burden of PC attributable to smoking from 1990 to 2021
Globally, the absolute burden of PC attributable to smoking has increased compared to 1990 (Table 1). The total DALYs rose from 1030311.87(95% UI,935957.37 - 1124871.5) in 1990 to 1789502.84(95% UI,1567220.76 - 2042056.59) in 2021. Similarly, the number of deaths increased from 38430.87(95% UI,34757.93 - 42214.02)in 1990 to 72170.03(95% UI,62852.94 - 82937.45) in 2021. However, after adjusting for time and age structure effects, the overall burden of PC attributable to smoking showed a certain degree of downward trend. The ASDR per 100,000 people decreased from 24.98(95% UI,22.66 -27.32) in 1990 to20.38(95% UI,17.84-23.25)in 2021 (AAPC,-0.65 [95% CI,-0.72 –0.58], P = 0) and the ASMR per 100,000 people decreased from 0.98(95% UI,0.89-1.08)in 1990 to 0.83(95% UI,0.73-0.96)in 2021(AAPC,-0.51[95% CI,-0.62 - -0.4], P = 0). Overall, on a global scale, despite noticeable increases in the total DALYs and the number of deaths compared to 30 years ago, after adjusting for the effects of time and age structure, downward trends have been observed both in ASDR and ASMR(Table 1).

Table 1. Global DALY and mortality of PC attributable to smoking and their AAPC by gender, SDI level, and region.
The regional burden of PC attributable to smoking from 1990 to 2021
Regionally, the burden of PC attributable to smoking varied significantly across the 21 geographic super-regions(Table 1, Figure 1). East Asia and Oceania each recorded the highest and the lowest number of DALYs in both 1990 and 2021, respectively. However, the highest and the lowest records of the ASDR in both 1990 and 2021 are held by Central Europe and Western Sub-Saharan Africa, respectively, with Central a downward trend in Europe and a notable upward trend in Western Sub-Saharan Africa. Rising from 1.57 per 100,000 (95% UI,1.29 -1.86) in 1990 to 2.17 per 100,000(95% UI, 1.73-2.66)in 2021, with an AAPC of 1.06 ([95% CI, 1.01 - 1.11], p=0), Western Sub-Saharan Africa experienced the most significant increase in ASDR during this period. Conversely, decreasing from 15.48 per 100,000(95% UI, 13.97 - 17.11) in 1990 to 9.77 per 100,000(95% UI, 8.22 - 11.35)in 2021, with an AAPC of -1.52([95% CI,-1.88 - -1.16], p=0), Central Latin America experienced the largest declines in ASDR. Deaths and ASMR were generally higher across all regions, with East Asia recording the highest number of mortality and Central Europe recording the ASMR in 2021. It is worth mentioning that, similar to its ASDR, although Western Sub-Saharan Africa showed the lowest ASMR in both 1990 and 2021, it exhibited a significant upward trend in ASMR with an AAPC of 1.08 ([95% CI,1.03 - 1.14], p=0), second only to Central Asia, which had the highest increase with an AAPC of 1.14 ([95% CI, 0.63 - 1.66], p=0). Similar to ASDR, the largest decline in ASMR occurred in Central Latin America, decreasing from 0.63(95% UI, 0.56-0.7)in 1990 to 0.39(95% UI, 0.33-0.46), with an AAPC of -1.54 ([95% CI,-1.85 to - 1.23], p=0).

Figure 1. The burdens of PC attributable to smoking among 204 countries and territories in 2021. (A) Regional distribution of the ASDR of PC attributable to smoking in 2021. (B) Regional distribution of trends in ASDR from 1990 to 2021. (C) Regional distribution of the ASMR of PC attributable to smoking in 2021. (D) Regional distribution of trends in ASMR from 1990 to 2021. (PC, pancreatic cancer; ASDR, age-standardized disability-adjusted life year; ASMR, age-standardized mortality rate).
National burden of PC attributable to smoking from 1990 to 2021
Nationally, the burden of PC attributable to smoking varied significantly across different countries and territories. In 2021, the ASDR and ASMR per 100,000 people ranged from 0.47(95% UI, 0.34 - 0.62)and 0.02(95% UI, 0.01 - 0.03)in Nigeria to 83.2(95% UI, 66.05 - 102.38) and 3.15(95% UI, 2.46 - 4) in Greenland. In the rankings for ASDR and ASMR in 2021, countries Greenland, Greece, Montenegro, and Armenia are all among the top 5(Figures 1A, C). Despite the highest ASDR or ASMR globally, China and the United States of America had the highest absolute numbers of DALYs and deaths, driven by their large populations. For a similar reason, the Russian Federation, Japan, Germany, and India are all ranked in the top six for both of the aforementioned indicators. From 1990 to 2021, the burden of PC attributable to smoking has shifted significantly across different countries (Figures 1B, D). San Marino and Ireland experienced the most substantial decreases in both ASDR (AAPCSan Marino = -2.49 [95% CI, -2.89 - -2.08], P < 0.001, AAPCIreland = -2.43 [95% CI, -3.05 - -1.81], P < 0.001) and ASMR (AAPCSan Marino = -2.61[95% CI, -2.82 - -2.4), P < 0.001, AAPCIreland = -2.43 [95% CI, -2.89 - -1.98, P < 0.001), while Turkmenistan and Cabo Verde experienced the most substantial increases in both ASDR (AAPCTurkmenistan = 8.03 [95% CI, 5.52 - 10.6], P < 0.001, AAPCCabo Verde = 7.33 [95% CI, 6.26 - 8.4], P < 0.001)) and ASMR (AAPCTurkmenistan = 8.17 [95% CI, 5.66 - 10.74), P < 0.001, AAPCCabo Verde = 7.46 [95% CI, 6.12 - 8.81], P < 0.001). In contrast, Fiji, Germany, and the United States Virgin Islands exhibited relatively stable ASDR and ASMR, with the 95%CI for their AAPCs clustering around zero, indicating minimal changes in the disease burden(Supplementary Table S1). Overall, the trends in ASDR closely mirrored the trends in ASMR, reflecting consistent patterns in the burden of smoking-related PC across nations.
PC burden attributable to smoking by sex and age group
In the sex subgroup analysis, the burden of PC attributable to smoking was significantly higher in males compared to females (Figure 2). In the age subgroup analysis in 2021, DALY rates rose with age, peaking at 70-74 years for both sexes, at 199.4/100,000 (95% UI, 170.18 - 234.02) for males and 47.5/100,000 (95% UI, 39.67 - 56.89)for females, and then declined. The trend in mortality rates was similar for both sexes, with a peak at 90-94 years 14.82/100,000 (95% UI, 11.41 - 18.81) for males and 4.81/100,000 (95% UI, 3.25 - 6.42) for females (Supplementary Table S2). From 1990 to 2021, although DALYs for PC attributable to smoking rose to some extent each year, the ASDR reduced significantly for both sexes, with a more pronounced decrease in females (AAPCASDR = -1.01 [95% CI, -1.11 to -0.9], P < 0.001). Likewise, the number of deaths rose significantly for both sexes, but the ASMR experienced a decline to a certain extent for both sexes, with a more significant decrease in females (AAPCASMR = -0.94 [95% CI, -1.04 - -0.84], P < 0.001) (Figure 2, Table 1).

Figure 2. The burdens of PC attributable to smoking of different age groups from 1990 to 2021. The DALYs (A), ASDR (B), Death (C), mortality rate (D) of different age groups from 1990 to 2021. (PC, pancreatic cancer; ASDR, age-standardized disability-adjusted life year; ASMR, age-standardized mortality rate).
PC burden attributable to smoking by SDI
According to the SDI classification in 2021, both ASDR and ASMR exhibited some degree of positive correlation with SDI levels across most regions, with the high SDI category displaying the highest ASDR and ASMR (Supplementary Figure S1). Although joint-point analysis indicated that this positive correlation persisted in most areas, further analysis revealed significant differences in the relationship between ASR and SDI across the 21 geographic regions and 204 countries (Figures 3A–D). At the 21 geographic regions level, contrary to the majority of other regions, the ASR exhibited a negative correlation with the SDI level in high-income Asia Pacific, High−income North America, Western Europe, and Australasia (Figures 3A, B). At the national level (Figures 3C, D), most countries fulfilled the condition of a positive correlation between ASRs and SDI levels, but there were still a few countries with significant deviations. From 1990 to 2021, changes in ASRs for PC attributable to smoking have been inconsistent across SDI subregions (Table 1). In most regions with different SDI levels, except in the Low-middle SDI region, both ASDR and ASMR showed a decreasing trend. Particularly, the decline in ASDR intensified with increasing SDI levels, with the most significant decrease in ASDR occurring in the High SDI region, reducing from 40.38/100,000 (95% UI, 36.80 - 43.96) in 1990 to 31.07/100,000 (95% UI, 27.49 - 34.78) in 2021, AAPC = -0.86 [95% CI, -0.98 - -0.74]. Similarly, the most significant decrease in ASMR also occurred in the High SDI region, from 1.6/100,000 (95% UI, 1.44 - 1.76) in 1990 to 1.28/100,000 (95% UI, 1.11 - 1.47) in 2021, with an AAPC = -0.73 [95% CI, -0.85 - -0.6]. Notably, the low SDI region also showed a certain degree of reduction in ASMR, more pronounced compared to other regions, with an AAPC of -0.1 [95% CI, -0.21 - 0] (Table 1). The above AAPC trends were also reflected at the national level(Figures 3E, F).

Figure 3. The correlation between SDI and ASR in 21 regions and 204 countries. The correlation between SDI and ASDR (A) or ASMR (B) in 21 regions. The correlation between SDI and ASDR (C) or ASMR (D) in 204 countries. The correlation between mean SDI and the AAPC of ASDR (E) or the AAPC of ASMR (F) in 204 countries. (ASDR, age-standardized disability-adjusted life year; ASMR, age-standardized mortality rate; AAPC, average annual percentage change).
Trends in ASDR and ASMR across time, sex, and SDI
Globally, the overall burden of PC attributable to smoking showed a certain degree of downward trend from 1990 to 2021 (Table 1). Based on different SDI categories, apart from a more pronounced downward trend in ASDR and ASMR shown in the High SDI region and a more pronounced upward trend in the Low-middle SDI region, the trends in other SDI regions are almost uniform, remaining stable or slightly decreasing(Supplementary Figure S1, Supplementary Table S3). Further analyzing the changing trends of disease burden by gender, we observed that in the High SDI region, the declines in ASDR and ASMR for males were more pronounced than for females (Figure 4). Maintaining a relatively stable level, compared to females, the ASDR and ASMR in High-middle and Middle SDI regions for males are relatively higher, which are higher than in Low-middle and Low SDI regions to a greater extent than females (Figure 4). The distinct levels of burden between genders are indicated.

Figure 4. The join-point analysis of ASR change by sex and SDI, 1990 to 2021. The join-point analysis for PC attributable to smoking in ASDR by SDI of male (A), female (B). The join-point analysis for PC attributable to smoking in ASMR by SDI of male (C), female (D). (PC, pancreatic cancer; ASDR, age-standardized disability-adjusted life year; ASMR, age-standardized mortality rate; SDI, Socio-demographic Index).
Decomposition analysis
Our decomposition analysis provided insights into the relative contributions of aging, population growth, and epidemiological changes to DALYs and mortality of PC attributable to smoking according to various regions, intending to clearly highlight the primary driving factor (Figure 5). According to five SDI regions, High-middle SDI and Middle SDI regions had obvious increases in DALYs and death numbers, primarily driven by population growth followed by aging (Figures 5A, B). In contrast, though facing population growth and aging of almost the same magnitude as the High-middle SDI and Middle SDI regions, the increase in DALYs and mortality numbers in High SDI regions is less, largely due to the significant positive impact of the epidemiological changes. According to the 21 GBD regions, in addition to the major distributions of population growth and aging, unfavorable epidemiological changes also contributed to the noticeable increase in the above indicators in East Asia. Conversely, epidemiological changes in High-income Asia Pacific, High-income North America, and Western Europe played a role in mitigating the increase in disease burden to some extent (Figures 5C, D). When it comes to the national level, we were primarily concerned with the top six nations with the largest increases in DALYs and death numbers, where the significant increases were mainly due to the substantial negative impact of epidemiological changes (Figures 5E, F). In terms of gender, there were no significant differences between males and females.

Figure 5. The decomposition analysis of PC burden attributable to smoking, 1990 to 2021. Changes in PC attributable to smoking DALYs (A) and mortality (B) according to population-level determinants of ageing, population growth, and epidemiological change from 1990 to 2021 at the global level by SDI quintile. Changes in PC attributable to smoking DALYs (C) and mortality (D) according to population-level determinants of ageing, population growth, and epidemiological change from 1990 to 2021 at the 21 regions level. Changes in PC attributable to smoking DALYs (E) and mortality (F) according to population-level determinants of ageing, population growth, and epidemiological change from 1990 to 2021 of the top six countries with the most DALYs and mortality increased. The black dot represents the overall value of change contributed by all 3 components. For each component, the magnitude of a positive value indicates a corresponding increase in PC attributable to smoking attributed to the component; the magnitude of a negative value indicates a corresponding decrease in PC attributable to smoking attributed to the related component. (PC, pancreatic cancer; DALYs, Disability-adjusted life-years).
Projection analysis
To predict trends in the age-standardized DALY rate and mortality rate from 2020 to 2035, BAPC was performed. Globally, both ASDR and ASMR are projected to continue a similarly pronounced downward trend, reflecting the consistency between the two indicators (Figure 6A, Supplementary Table S4). Specifically for five regions, the trends of ASDR and ASMR are also similar. The decreasing trends of ASDR and ASMR are more pronounced in high and high-middle SDI regions, especially in high SDI regions, while there is a certain degree of increase in middle SDI areas, which needs to be brought to our attention, and remains relatively stable in other regions (Figures 6B, C, Supplementary Table S5). To better draw attention to the priority geographic super-regions and facilitate the implementation of effective measures at a more specific level, we have further predicted the trend of disease burden changes in subsequent years in the seven geographic super-regions with the highest disease burden in 2021. The results indicate that although the overall ASDR and ASMR levels in these 7 regions are significantly higher than the global average, it is encouraging to see that the majority of them are experiencing a downward trend over the coming years. However, it is concerning that the disease burden continues to rise in East Asia, with projections estimating the ASDR per 100,000 people at 60.70 and the ASMR per 100,000 people at 2.42 by 2035, making it the region with the highest disease burden (Figures 6D, E, Supplementary Table S6). It is worth noting that in some regions with sparse data or low predicted values, the lower limit of the 95% uncertainty interval is slightly below zero. This reflects the statistical uncertainty in the model estimates rather than implying negative values in reality.

Figure 6. The predictive analysis of PC burden attributable to smoking,1990 to 2035. The global predictive analysis of PC attributable to smoking in ASDR and ASMR (A). The predictive analysis of PC attributable to smoking of ASDR (B) and ASMR (C) in different SDI regions. The predictive analysis of PC attributable to smoking of ASDR (D) and ASMR (E) in the seven geographic super-regions with the highest disease burden in 2021.
Discussion
As one of the leading causes of cancer-related mortality, PC remains a significant global health issue. Smoking continues to be a potentially modifiable lifestyle risk factor with the strongest correlation with PC. It is necessary and meaningful to conduct disease studies on the burden of disease of PC attributable to smoking. In our study, for the first time, we have detailed the global burden and evolving trends of PC attributable to smoking by leveraging the latest data from the GBD 2021. We not only reveal the differences in smoking-attributable PC DALYs and mortality between sexes, ages, regions, countries/territories, and SDIs but also further analyze the main factors affecting the burden of disease and predict the trend of the burden of disease in the coming period.
Smoking remains the most detrimental form of tobacco use. Previous basic research has identified that smoking plays a role in initiating the onset of PC and accelerating its progression, invasion, and metastasis (10, 11). In addition, this finding is corroborated by previous clinical studies showing a higher incidence rate and a lower survival rate after onset of the disease in smokers (7, 8). Since the implementation of effective tobacco control programs and policies, such as the WHO FCTC, the global prevalence of smoking has significantly declined (20, 21) despite the uneven distribution among regions, and smoke-free environments are increasingly prevalent worldwide (22).
Globally, despite the increasing absolute burden of the disease, after adjusting for time and age structure effects, the burden of the disease has shown a downward trend from 1990 to 2021, suggesting that epidemiological change had played a positive role. Around 2010, 5years following the implementation of the WHO FCTC, noticeable decreasing trends in both the global ASDR and ASMR were indicated by the Join-point regression model, which was likely related to the period of fastest decrease in the prevalence of smoking tobacco use across the largest number of countries. It is worth mentioning that the APCs of ASDR and ASMR in recent years are higher than in the past, suggesting that the progress made in tobacco control in recent years has played a remarkably effective role in the prevention and control of PC attributable to smoking (23). This has been confirmed in our decomposition analysis section. Encouragingly, our predictive analysis indicates a continued decline in ASDR and ASMR for the disease globally over the next few years. The realization of this prediction is contingent on the premise that continued implementation of strong tobacco control measures will be taken to create favorable epidemiological conditions that can neutralize the adverse effects of population aging and demographic changes on the burden of disease.
Despite the variations of the impact of smoking on the PC burden by regions and countries/territories, males bear a significantly greater burden of disease than females, and the trend of burden reduction is more pronounced in females than in males, which should be considered in cancer prevention and early screening. According to a systematic analysis of the prevalence of smoking tobacco use in 2019, this is highly likely due to the lower prevalence, shorter duration, and lower intensity of smoking in females than in males as well as the more significant reduction in global age-standardised prevalence of smoking tobacco use from 1990 among females (15). Moreover, the increased activation of androgen receptor expression in the PC cell lines and their surrounding tissues might be involved in the carcinogenesis and cancer development of PC and further lead to a higher burden of disease in males (24).
According to the prevalence of tobacco use from the GBD 2019, East Asia is the super-region with the largest number of tobacco smokers in 2019 due to its large population base (15). The highest number of DALYs and mortality in East Asia in our study is likely related to this situation. After adjusting for time and age structure effects, with its high age-standardized prevalence of smoking tobacco use (15), Central Europe has the highest ASDR and ASMR, consistent with patterns observed in other smoking-related diseases. With the large decreases in age-standardized prevalence of smoking tobacco use in Colombia and Costa Rica, the countries have continuously introduced anti-smoking laws that have severe penalties and high costs associated with consumption, the largest declines in ASDR and ASMR were observed in Central Latin America, showing the potential of the effective tools to operate to greatly reduce the prevalence of smoking tobacco use and save millions of lives over the coming decades. With almost the lowest tobacco consumption per person, the ASDR and ASMR of smoking tobacco in some economically underdeveloped regions, such as sub-Saharan Africa, have always maintained the lowest levels, but with the largest relative increases in the number of smokers during the period (15). The manifestation that Western Sub-Saharan Africa is experiencing almost the most significant increase in ASDR and ASMR warnings us that it is high time that strong tobacco control policies should be implemented and enforced to slow down the rate of increase in the burden of disease in this area with the economic recovery and growth.
With almost the lowest and the highest age-standardized prevalence of current use of smoking tobacco, Nigeria and Greenland have the corresponding levels of ASDR and ASMR. With 341 million (30%) of 1·14 billion tobacco smokers globally, China is the country with the largest number of tobacco smokers in 2019 likely because the reduction in prevalence has not kept pace with population growth. Sufficient decreases in the prevalence of smoking are needed to offset the large population base and the demographic force of population growth, further reducing the absolute burden of disease. Respectively approving the WHO FCTC in 2004 and 2005, becoming two of the early parties to the “convention”, and have implemented a series of policies to reduce tobacco consumption and related health risks, such as smoke-free environments, tobacco tax policies, and restrictions on tobacco advertising and promotion, San Marino and Ireland experienced the most substantial decreases in both ASDR and ASMR. Reiterating, proactive tobacco control policies play a significant role in reducing the burden of this disease.
Focus on the trends of the past 5-10 years, the trends in PC attributable to smoking have varied significantly across SDI regions, characterized by declining ASDR and ASMR in high SDI regions, while observing a mild increase or relative stability in low and low-middle SDI regions and a gradual increase in medium SDI regions. This divergence can be attributed to several factors: Countries with high SDI typically enforce stringent tobacco regulations (25), promote physical activity, optimize dietary habits, and possess advanced healthcare systems with robust screening capabilities, and this is also one of the reasons why the burden of disease in high-income regions such as Asia-Pacific and North America is inversely proportional to the level of SDI. Conversely, In low and low-middle SDI regions, the current level of tobacco control policy implementation is far from sufficient. Most of the countries that have yet to adopt a single MPOWER measure at the highest level of achievement belong to low and low-middle SDI regions. Of those 20 countries where cigarettes have become more affordable, 17 belong to low and low-middle SDI regions (26). The affordability of cigarettes and the lack of necessary restrictive measures have led to a high smoking prevalence in the regions, and coupled with its relatively backward medical level, the burden of disease has gradually increased. The higher burden of disease in medium SDI regions may be linked to improving economic conditions, which have led to increased smoking prevalence and the gradual enhancement of physical examination. The notable regional disparities in the burden of disease across SDI levels underscore socio-spatial inequities in the prevention and management of PC. Future efforts should focus on not only taking effective tobacco control measures such as raising taxes but also bolstering healthcare capacity in low and middle-income regions, alongside interregional collaborations aimed at promoting smoking cessation.
The main drivers behind changes in the burden of disease are threefold: age structure, population growth, and epidemiological changes. Each factor affects the burden of disease differently across regions and periods. The preliminary conclusion that a higher incidence rate is among the elderly enlightens us to observe the weight of each factor and highlight the primary driving factors. Our study indicated that the increases in DALYs and death numbers among most regions and countries were primarily driven by population growth followed by aging, and the variations in the growth rates of burdens of disease across different regions were influenced by the distinct roles that epidemiological changes had played. The epidemiological changes in High SDI areas have shown a positive impact in mitigating the increase in the absolute burden of disease, while the nations experiencing the most rapid burden growth have shown negative effects. Within the framework of the decomposition analysis model, “epidemiological change” refers to shifts in disease incidence, prevalence, and mortality rates, influenced by various factors such as pathogen variations (e.g. the emergence of highly virulent strains), the effectiveness of public health interventions (e.g., vaccination and antibiotic use), lifestyle changes (e.g. diet and exercise habits), environmental factors (e.g. climate change), and socioeconomic changes (e.g. access to and quality of healthcare resources). Epidemiological changes affect the spread and severity of diseases, thereby altering the burden of disease. Smoking as a potentially modifiable lifestyle risk factor with the strongest correlation with PC is an important component of epidemiological changes. Thus, under the general trend of population growth and aging, taking active measures to change epidemiological factors, such as restricting tobacco use, is critical to slow down the growth of the absolute burden of disease.
It is projected that over the next 15 years, most regions will continue their current trends in burden of disease changes. However, it is a cause for concern that East Asia, which is currently among the top 7 countries with the highest burden of disease and facing severe population growth and aging issues, will continue to show a significant upward trend. If effective intervention measures are not taken in time, it is anticipated that by 2035, this region will become the area with the highest burden of disease. Therefore, strengthening disease prevention efforts by restricting the use of tobacco and increasing investment in healthcare resources for PC attributable to smoking are imminent.
What needs to be clarified is that because the original purpose of DALYs was to emphasize the greater impact of death or disability in younger age groups on social productivity, it is inevitable that the burden of disease in younger individuals may be overinflated compared to it in older individuals, based on the calculation formula for DALYs mentioned earlier. Therefore, when describing the burden of disease, we introduce the ASDR indicator to adjust the time and age structure effects. When assessing the relationship between the burden of disease and age, similar to some top-tier articles that also stratify age into five-year intervals, we simply presented DALYs by age group, without comparing DALYs between groups with large age differences because of a certain degree of unfairness. Presenting DALYs by age group has its corresponding value. On the one hand, under the same number of disabilities, the DALYs calculated for younger populations should be higher than for older groups. However, our study data indicated that the DALYs in younger groups were significantly lower than in older groups, which further suggested that the burden of PC attributable to smoking in young people was significantly lower. On the other hand, the DALYs between two age groups with close ages are reasonably comparable. According to our data, the inter-group DALYs differences among the 45-49, 50-54, and 55-59 age groups were more pronounced than the differences amongother adjacent age groups, reflecting a clear trend of the increasing burden of disease with age in the middle-aged population. This is not only related to the age characteristics of the peak incidence of PC itself but also the cumulative increase in tobacco use over the years, according to current research which has shown that the risk of PC in smokers significantly increases with the amount of daily smoking, the number of years smoked, and the cumulative number of packs smoked (27).
Several limitations should be noted. Firstly, GBD data, including its complex modeling and data input, significant methodological changes with each revision, and the potential for assumptions and knowledge gaps to influence estimates. Additionally, the GBD database only controls for certain covariates, meaning that unadjusted covariates (such as dietary habits, lifestyle, medication use, metabolic and physiological changes, etc.) may influence the study results (28). Secondly, based on the insidious onset inherent characteristics of PC, there may be misdiagnosed and underreporting of cases in some regions and countries/territories with underdeveloped healthcare systems, especially in low-income areas, leading to a lack of high-quality and detailed data, potentially leading to an underestimation of the burden of disease. Thirdly, the exposure to tobacco may be underestimated to some extent. On the one hand, the data were sourced from public databases, and information regarding smoking habits relied on self-reporting, potentially resulting in an underestimation of smoking prevalence, particularly among females and younger individuals. On the other hand, the databases assessed only traditional smoking tobacco products, excluding emerging tobacco products like e-cigarettes, whose health implications are increasingly recognized (29). Consequently, reported cases of PC attributed to smoking in public health records may be lower than the actual prevalence. Prospective studies utilizing real-world data are warranted to assess the risk of PC among smokers. Such studies would provide more accurate assessments of the burden of PC attributable to smoking and inform targeted interventions to mitigate its public health impact.
Conclusions
We analyzed the global burden of PC attributable to smoking and its variations across different regions, countries/territories, sexes, and age groups. Despite the decline in global smoking prevalence, PC attributable to smoking remains a significant issue, especially in high and high-middle-income regions and among males. Additionally, the progressive increase burden of disease in low-middle-income regions also deserves our attention. In the context of the general trend of population growth and aging, epidemiological factors have played an important role in influencing the trends in the burden of disease. Therefore, targeted tobacco control measures are essential to mitigate this persistent health burden, focusing on regions experiencing rising smoking prevalence and associated disease impacts.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: The dataset supporting the conclusions of this article is available in the Global Health Data Exchange repository. (https://vizhub.healthdata.org).
Author contributions
RD: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. YW: Formal analysis, Validation, Visualization, Writing – review & editing. MP: Validation, Writing – review & editing. JZ: Validation, Writing – review & editing. YZ: Validation, Writing – review & editing. CZ: Validation, Writing – review & editing. CL: Resources, Supervision, Writing – review & editing. YG: Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Cultivating Fund for NSFC of Shandong Provincial Qianfoshan Hospital (QYPY2019NSFC0603) and the Natural Science Foundation of Shandong Province (ZR2016HL03). The funders had no role in the study design, data collection, analysis, decision to publish, or manuscript preparation.
Acknowledgments
We appreciate the works of the Global Burden of Diseases, Injuries, and Risk Factors Study collaborators.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2025.1547029/full#supplementary-material
Abbreviations
PC, pancreatic cancer; GBD, the Global Burden of Diseases, Injuries, and Risk Factors Study; SDI, socio-demographic index; DALYs, disability-adjusted life years; ASDR, age-standardized disability-adjusted life years rates; ASMR, age-standardized mortality rates; AAPCs, average annual percentage changes; BAPC, the Bayesian Age-Period-Cohort; SEER, the Surveillance, Epidemiology, and End Results; IARC, the International Agency for Research on Cancer; MDSCs, myeloid-derived suppressor cells; HB-EGF, heparin-binding epidermal growth factor-like growth factor; CSE, cigarette smoke extract; CBX3, Chromobox protein homolog 3; YBX1, Y-box-binding protein 1; HDAC4, histone deacetylase 4; WHO, FCTC the WHO Framework Convention on Tobacco Control; YLDs, years lived with disability; YLLs, years of life lost; HALE, healthy life expectancy; IHME, the Institute for Health Metrics and Evaluation; ASRs, age-standardized rates; UIs, uncertainty intervals; CIs, confidence intervals; APC, annual percentage change.
References
1. Hu J-X, Zhao C-F, Chen W-B, Liu Q-C, Li Q-W, Lin Y-Y, et al. Pancreatic cancer: A review of epidemiology, trend, and risk factors. World J Gastroenterol. (2021) 27:4298–321. doi: 10.3748/wjg.v27.i27.4298
2. Kleeff J, Korc M, Apte M, La Vecchia C, Johnson CD, Biankin AV, et al. Pancreatic cancer. Nat Rev Dis Primer. (2016) 2:16022. doi: 10.1038/nrdp.2016.22
3. NATIONAL CANCER INSTITUTE Surveillance, Epidemiology, and End Results Program. Cancer Stat Facts: Pancreatic Cancer. Available online at: https://seer.cancer.gov/statfacts/html/pancreas.html (Accessed September 19, 2025).
4. American Cancer Society Journals. Wiley Online Library. Available online at: https://acsjournals.onlinelibrary.wiley.com/ (Accessed September 21, 2024).
5. Ferrari AJ. Global incidence, prevalence, years lived with disability (YLDs), disability-adjusted life-years (DALYs), and healthy life expectancy (HALE) for 371 diseases and injuries in 204 countries and territories and 811 subnational locations, 1990–2021: a systematic analysis for the Global Burden of Disease Study 2021. Lancet Lond Engl. (2024) 403:2133–61. doi: 10.1016/S0140-6736(24)00757-8
6. Klein AP. Pancreatic cancer epidemiology: understanding the role of lifestyle and inherited risk factors. Nat Rev Gastroenterol Hepatol. (2021) 18:493–502. doi: 10.1038/s41575-021-00457-x
7. Bosetti C, Lucenteforte E, Silverman DT, Petersen G, Bracci PM, Ji BT, et al. Cigarette smoking and pancreatic cancer: an analysis from the International Pancreatic Cancer Case-Control Consortium (Panc4). Ann Oncol. (2012) 23:1880–8. doi: 10.1093/annonc/mdr541
8. Yuan C, Morales-Oyarvide V, Babic A, Clish CB, Kraft P, Bao Y, et al. Cigarette smoking and pancreatic cancer survival. J Clin Oncol. (2017) 35:1822–8. doi: 10.1200/JCO.2016.71.2026
9. Kumar S, Torres MP, Kaur S, Rachagani S, Joshi S, Johansson SL, et al. Smoking accelerates pancreatic cancer progression by promoting differentiation of MDSCs and inducing HB-EGF expression in macrophages. Oncogene. (2015) 34:2052–60. doi: 10.1038/onc.2014.154
10. Zhang H, Yu H, Ren D, Sun Y, Guo F, Cai H, et al. CBX3 regulated by YBX1 promotes smoking-induced pancreatic cancer progression via inhibiting SMURF2 expression. Int J Biol Sci. (2022) 18:3484–97. doi: 10.7150/ijbs.68995
11. Yang J, Chheda C, Lim A, Hauptschein D, Zayou L, Tang J, et al. HDAC4 mediates smoking-induced pancreatic cancer metastasis. Pancreas. (2022) 51:190–5. doi: 10.1097/MPA.0000000000001998
12. WHO report on the global tobacco epidemic, 2023: protect people from tobacco smoke. Available online at: https://www.who.int/publications/i/item/9789240077164 (Accessed August 15, 2024).
13. Murray CJL. The Global Burden of Disease Study at 30 years. Nat Med. (2022) 28:2019–26. doi: 10.1038/s41591-022-01990-1
14. GBD Results. Institute for Health Metrics and Evaluation. Available online at: https://vizhub.healthdata.org/gbd-results (Accessed October 8, 2024).
15. Spatial, temporal, and demographic patterns in prevalence of smoking tobacco use and attributable disease burden in 204 countries and territories, 1990–2019: a systematic analysis from the Global Burden of Disease Study 2019. Lancet Lond Engl. (2021) 397:2337–60. doi: 10.1016/S0140-6736(21)01169-7
16. Zou J, Sun T, Song X, Liu Y-M, Lei F, Chen M-M, et al. Distributions and trends of the global burden of COPD attributable to risk factors by SDI, age, and sex from 1990 to 2019: a systematic analysis of GBD 2019 data. Respir Res. (2022) 23:90. doi: 10.1186/s12931-022-02011-y
17. Kim HJ, Fay MP, Feuer EJ, and Midthune DN. Permutation tests for joinpoint regression with applications to cancer rates. Stat Med. (2000) 19:335–51. doi: 10.1002/(sici)1097-0258(20000215)19:3<335::aid-sim336>3.0.co;2-z
18. Decomposition of the difference between two rates and its consistency when more than two populations are involved: Mathematical Population Studies. Math Popul Stud. (1991) 3:105–25. doi: 10.1080/08898489109525329
19. Santosa A, Wall S, Fottrell E, Högberg U, and Byass P. The development and experience of epidemiological transition theory over four decades: a systematic review. Glob Health Action. (2014) 7:23574. doi: 10.3402/gha.v7.23574
20. Puska P and Daube M. Impact assessment of the WHO Framework Convention on Tobacco Control: introduction, general findings and discussion. Tob Control. (2019) 28:s81–3. doi: 10.1136/tobaccocontrol-2018-054429
21. Lee K, Egbe CO, Bianco E, and Arora M. The 20th anniversary of the WHO Framework Convention on Tobacco Control: hard won progress amid evolving challenges. The Lancet. (2023) 402:592–4. doi: 10.1016/S0140-6736(23)01080-2
22. Semple S, Dobson R, O’Donnell R, Zainal Abidin E, Tigova O, Okello G, et al. Smoke-free spaces: a decade of progress, a need for more? Tob Control. (2022) 31:250–6. doi: 10.1136/tobaccocontrol-2021-056556
23. Global progress in implementation of the WHO FCTC: report by the Convention Secretariat (2021). Available online at: https://iris.who.int/handle/10665/368586 (Accessed September 21, 2024). Control WFC on T, Parties C of the, Session N, Geneva, Switzerland, November 2021 8-13.
24. Kanda T, Jiang X, and Yokosuka O. Androgen receptor signaling in hepatocellular carcinoma and pancreatic cancers. World J Gastroenterol WJG. (2014) 20:9229–36. doi: 10.3748/wjg.v20.i28.9229
25. Reitsma MB, Gupta V, Ng M, and Gakidou E. The effects of tobacco control policies on global smoking prevalence. Nat Med. (2021) 27:239–43. doi: 10.1038/s41591-020-01210-8
26. WHO reports progress in the fight against tobacco epidemic. Available online at: https://www.who.int/news/item/27-07-2021-who-reports-progress-in-the-fight-against-tobacco-epidemic (Accessed September 11, 2024).
27. Lugo A, Peveri G, Bosetti C, Bagnardi V, Crippa A, Orsini N, et al. Strong excess risk of pancreatic cancer for low frequency and duration of cigarette smoking: A comprehensive review and meta-analysis. Eur J Cancer. (2018) 104:117–26. doi: 10.1016/j.ejca.2018.09.007
28. Tichenor M and Sridhar D. Metric partnerships: global burden of disease estimates within the World Bank, the World Health Organisation and the Institute for Health Metrics and Evaluation. Wellcome Open Res. (2020) 4:35. doi: 10.12688/wellcomeopenres.15011.2
Keywords: global health, smoking, pancreatic cancer, disability-adjusted life years, average annual percentage changes
Citation: Du R, Wang Y, Pan M, Zhu J, Zhao Y, Zhang C, Liu C and Gao Y (2025) Global, regional, and national burdens of pancreatic cancer attributable to smoking from 1990 to 2021 and the projections to 2035:a systematic analysis from the global burden of disease study 2021. Front. Oncol. 15:1547029. doi: 10.3389/fonc.2025.1547029
Received: 17 December 2024; Accepted: 05 May 2025;
Published: 30 May 2025.
Edited by:
Khadijah A. Mitchell, Fox Chase Cancer Center, United StatesReviewed by:
Sudhanshu Kumar Bharti, Patna University, IndiaManfred Neuberger, Medical University of Vienna, Austria
Copyright © 2025 Du, Wang, Pan, Zhu, Zhao, Zhang, Liu and Gao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Changhong Liu, bGl1Y2hhbmdob25nNTY3OUBvdXRsb29rLmNvbSA=; Yuan Gao, Z2FveXVhbkBzZGhvc3BpdGFsLmNvbS5jbg==
†These authors have contributed equally to this work
‡These authors share first authorship