- 1Department of Otolaryngology, Linyi Central Hospital, Linyi, China
- 2Department of Otolaryngology, Head and Neck Surgery, The Second Hospital of Weihai City, Weihai, China
- 3Department of Otolaryngology, Head and Neck Surgery, The Second Hospital of Dalian Medical University, Dalian, China
Background: Laryngeal cancer (LC) is a prevalent head and neck tumor, yet few studies quantify its disease burden and trends in risk factors of the age subgroup of middle-aged and older adult people (MAOP).
Methods: Based on Global Burden of Disease 2021 (GBD 2021) data, this study comprehensively assessed trends of disease burden in the age subgroup of LC in MAOP across 204 countries and regions from 1990 to 2021, covering incidence, prevalence, mortality, and disability-adjusted life years (DALYs). The estimated annual percentage change (EAPC) and joinpoint regression model were used to analyze the trends of disease burden and the changes in disease burden attributed to tobacco, alcohol, and occupational risk factors.
Results: Worldwide, the number of LC in MAOP prevalence, incidence, mortality, and DALYs totaled 887,500, 163,300, 97,500, and 2,284,800 in 2021. Over the previous 32 years, the rates of these indicators declined, with EAPC of −0.62 (−0.69 to −0.55), −0.95 (−1.04 to −1.87), −1.52 (−1.6 to −1.43), and −1.69 (−1.79 to −1.59), respectively. Regions with high Sociodemographic Index (SDI) showed high prevalence and morbidity. Still, they had low mortality and DALY rates, whereas areas with low SDI, such as Sub-Saharan Africa, experienced high mortality and DALY rates. Central Europe, a high-middle SDI region, was a “hotspot” for laryngeal cancer. Before the age of 75, males showed a more significant burden of disease indicators, with a peak incidence rate of 28.05 (23.93, 30.47) per 100,000 individuals between the ages of 85 and 90. However, the 85–89 and 90–94 age groups saw rising prevalence rates. The main risk factor for LC-related deaths and DALY in MAOP was tobacco exposure, particularly in low-middle SDI regions.
Conclusion: In the last 32 years, the burden of LC in MAOP has declined globally, yet age-related disparities and SDI stratification persist. Men aged before 75 and after 85 need to be targeted as key groups for prevention and control. Tobacco exposure was the leading risk factor, especially in low-middle SDI regions. These findings highlight the critical necessity for tailored public health strategies and policies.
1 Introduction
Due to its unique anatomical features, laryngeal cancer, a cancer affecting the head and neck region, represents a significant global health challenge by threatening patients’ respiratory, swallowing, and speaking abilities (1). Epidemiological statistics from the GBD 2019 study revealed that the total number of new cases of LC worldwide was 209,000, with males comprising 181,000 cases and females accounting for 28,500 instances (2). The epidemiologic profile of laryngeal cancer demonstrates significant demographic and geographic heterogeneity, with countries exhibiting relatively low SDI and MAOP male demographics experiencing the highest incidence of laryngeal cancer. Age-standardized mortality rates in developing countries (such as Pakistan, Seychelles, and Cuba) exceed 5 per 100,000, while those in developed countries (such as Japan, Singapore, and New Zealand) are below 0.5 per 100,000 (3). The disease burden is exceptionally high in MAOP. At present, the aging process is accelerating significantly around the world. According to the report by the World Health Organization, the population over 60 will grow from 1 billion to 1.4 billion from 2020 to 2030 (4).
However, regional differences still exist in the current prevention and control of LC. Statistically, one-third of the world’s poorest people have limited surgical coverage (3.5%), and there is a lack of consensus on the definition of surgery and surgical care (5). Numerous LMICs struggle with inadequate infrastructure, physical resources, and workforce necessary for delivering essential surgical care, challenging the idea that surgical issues are swiftly resolved. Furthermore, LMICs encounter substantial disparities in access to cancer treatment drugs. In many wealthy countries, patients can access cancer medications, including new treatments, without paying a substantial amount out of pocket. On the other hand, 40% of conventional chemotherapy medicines are only accessible to patients in LMICs (6).
As the immune system undergoes a gradual decline, compounded by long-term risk factors like tobacco use and alcohol intake (7), along with the significant impact of human papillomavirus (HPV) infection (8), the MAOP has become a high-risk group for LC. Despite noteworthy advancements in laryngeal cancer research within the biomedical field and the emergence of sophisticated analytical techniques, this malignancy continues to impose a substantial burden on the health of humans, particularly the older adult population. Middle-aged and older adult patients, often afflicted with underlying chronic diseases, frequently exhibit unsatisfactory prognoses, with a 5-year relative survival rate of approximately 61% (9). However, the burden of LC and its trends in MAOP, especially across distinct age subgroups, remain subjects of current understudies.
The GBD 2021 database is a global public health research project resource with a wide range of impact and scale. It systematically integrates data on 371 diseases, injuries, and risk factors from more than 200 countries and territories (10). The GBD database covers incidence, prevalence, mortality, DALYs for diseases and injuries, and attribution analysis of associated risk factors. The project provides a critical scientific basis for global public health policy, disease prevention and control development, and healthcare resource distribution, which helps to promote global health equity.
Therefore, based on GBD 2021 data, this study provides an in-depth analysis of epidemiologic characteristics, risk factors, projections, and trends in the burden of disease of LC in MAOP, especially in the age subgroups. The study is essential for formulating strategies for preventing and controlling laryngeal cancer and promoting the equitable distribution of medical resources. The study aims to achieve the following primary objectives: (1) to quantitatively analyze the spatial and temporal distributions and age subgroup distribution characteristics of the prevalence, incidence, mortality, and DALY rates of LC in MAOP from 1990 to 2021; (2) to evaluate the impact of significant risk factors across various age groups, genders, and regions; and (3) to propose stratified prevention and control strategies for LC in MAOP.
2 Materials and methods
2.1 Data sources
This research utilized the GBD 2021 database, which covers the epidemiological information of 371 diseases and injuries and 88 risk factors in 204 countries and regions worldwide between 1990 and 2021 (10). The database details, such as the construction method, index definition, and data collection channels, have been systematically discussed and explained in relevant academic literature. This study extracted epidemiologic data on LC in MAOP from this database for the world, five SDI regions, 21 GBD-defined regions, and 204 countries for people aged 55 years and older. The study covered four leading disease indicators: prevalence, morbidity, mortality, and DALY values and rate values, and included secondary risk factors: alcohol use, occupational risks, and tobacco use. The study’s dataset was sourced from the Global Health Data Exchange query tool.1
2.2 Research type and design
Based on the type of secondary data analysis, this study was divided into subgroups according to age. The disease burden indicators and attribution risk factors of LC in MAOP at the global, regional, and national levels in the GBD 2021 database were intensely discussed. The research subjects are individuals aged 55 and older, further divided into nine age subgroups (55–59, 60–64, 65–69, 70–74, 75–79, 80–84, 85–89, 90–94, 95+). Additionally, this study analyzes the association between SDI and disease burden.
2.3 Statistical analysis
2.3.1 Descriptive analysis
This study used descriptive analytic metrics covering numbers and rates of incidence, prevalence, mortality, and DALY for LC in MAOP. Mortality and DALY rates due to secondary attributable risk factors for LC in MAOP were also explored. The trends of disease burden indicators, including incidence, mortality, and DALY, were evaluated by estimating the EAPC, 95% confidence intervals (CIs), and their corresponding p-values. The formula for calculating EAPC is: EAPC = 100 × [exp(β) − 1] (11). Statistically significant increasing burden levels were identified when the EAPC point estimate and the lower confidence boundary exceeded zero. Conversely, an EAPC value with its upper CI limit below zero demonstrated a statistically significant decreasing pattern. A 95% CI range of 0 is considered to be relatively stable.
2.3.2 Correlation analysis
The SDI is a comprehensive indicator for measuring a country or region’s social and economic development. Based on SDI values, the world is divided into five levels: low SDI, low-middle SDI, middle SDI, high-middle SDI, and high SDI (12). This study explored the correlation between the prevalence, incidence, mortality, DALY rate, and SDI of LC in MAOP. Correlation coefficient (r) values were used to quantify the degree of linear correlation between SDI and the burden indicators of LC in MAOP. A positive association is evidenced when the r value exceeds the zero threshold, meaning the disease burden increases with rising SDI. A negative correlation is evidenced when the r value is below zero, meaning the disease burden increases with rising SDI.
2.3.3 Joinpoint regression analysis
The joinpoint regression model was used to analyze the temporal evolution trends of risk factors attributed to mortality and DALY rates of LC in MAOP. Differences among regions, genders, and age groups were analyzed. This study implemented significance testing for joinpoint configurations through the Monte Carlo permutation method. The evolutionary trend of the time series was analyzed using the average annual percentage change (AAPC) and annual percentage change (APC) indicators (13). APC metric quantifies temporal dynamics within a specific period, whereas the AAPC can evaluate the overall change trend during the entire observation period. When the values of APC and AAPC are above zero, they suggest a rising pattern; conversely, when these indicators are under zero, they denote a declining tendency.
When the p-value falls below 0.05, it implies that the difference has statistical significance. The analyses were executed through R (version 4.4.2), the ggplot2 package, and Joinpoint software.
3 Results
3.1 Global level
An overall decreasing trend in the rate of prevalence, incident, mortality, and DALYs of laryngeal cancer related to middle-aged and older adults was reported globally. Compared to the prevalence rate and incidence rate, the overall trend of mortality and DALYs rates was more positive (Figure 1). The prevalent cases increased from 466,713 cases (442,851 to 491,527) in 1990 to 887,580 cases (832,463 to 952,557) in 2021. The incident cases increased from 94,791 (90,077 to 99,538) in 1990 to 163,316 (151,965 to 175,890) in 2021. The rates of deaths and DALYs in 2021 were 6.56 (6.12, 7.7) per 100,000 population and 153.76 (143.31 to 165.63) per 100,000 population, with EAPC of −1.52(−1.6 to −1.43), and −1.69(−1.79 to −1.59), respectively (Figure 1 and Table 1).

Figure 1. Disease trends by 5 SDI regions and gender for the prevalence rate, incidence rate, death rate, and DALYs rate from 1990 to 2021. DALYs, disability-adjusted life-year; SDI, sociodemographic index.

Table 1. The Incidence and DALYs of LC in MAOP in 1990 and 2021, and the trends from 1990 through 2021.
3.2 SDI regional level
The period from 1990 to 2021 showed a general decrease in prevalence, incidence, mortality, and DALY rates in the five SDI regions (Figure 1 and Table 1). Notably, the low-middle SDI region showed the most minor reductions in all four measures. In contrast, the high-middle SDI region exhibited significantly more significant reductions in prevalence, mortality, and DALY rates than the other regions. In 2021, the high SDI region demonstrated the highest prevalence case and morbidity rates, with 268,129 (253,063 to 279,715) and 77.72 (73.35 to 81.07) per 100,000, respectively; in High-middle SDI, the incidence rate peaked at 12.15 (10.96 to 13.4) per 100,000. The region exhibiting the highest mortality and DALY rates was identified as low-middle SDI, at 9.9 (8.91 to 11.17) and 239.34 (215.3 to 271) cases for every 100,000 people, correspondingly, as shown in Figure 1 and Table 1. Furthermore, the middle SDI region showed the most significant cases of incidence, deaths, and DALYs, with respective figures of 44,598 (39,812 to 49,987), 29,787 (26,773 to 33,198), and 694,235 (623,988 to 773,960), respectively.
3.3 GBD regional levels
During the last 32 years, among the 21 regions, mortality and DALY rates dropped, while prevalence and morbidity demonstrated an upward trend in certain regions. Specifically, seven regions exhibited increasing trends in prevalence: East Asia, the Caribbean, Southeast Asia, Central Europe, North Africa and the Middle East, Tropical Latin America, and Western Sub-Saharan Africa (Figure 2 and Table 1). Notably, increases in prevalence were exclusively observed in the Caribbean, East Asia, and Southeast Asia, with EAPCs of 0.43 (0.32 to 0.53), 0.19 (0.04 to 0.34), and 0.07 (0.03 to 0.11), respectively. These regions also demonstrated increasing prevalence trends (Figure 2 and Table 1).

Figure 2. Disease trends in global and 21 regions from 1990 to 2021. The EAPC of prevalence rate (A), incidence rate (B), death rate (C), and DALYs rate (D). DALYs, disability-adjusted life-year; GBD, Global Burden of Disease; EAPC, estimated annual percentage changes.
From 1990 through 2021, minor reductions in mortality and DALY rates have been observed in the Caribbean, Western Sub-Saharan Africa, and Southern Sub-Saharan Africa, with the most significant reductions in prevalence, incidence, and mortality demonstrated in Central Asia (Figure 2). The regions with the most considerable DALY rate reductions were Australia, high-income Asia, and Central Pacific Latin America, respectively (Figure 2D). Notably, Central Europe, a region with a high-moderate SDI, had the highest prevalence rate and incidence rate in 2021 at 111 (102.45 to 120.09) and 20.27 (18.55 to 22.09) per 100,000 people, as well as the second highest rate of deaths and DALYs (Figure 2 and Table 1). The number of DALYs reported in South Asia in 2021 was 709,622 (618,985 to 816,799) cases, about one-third of the global number of DALYs.
3.4 Country level
Between 1990 and 2021, the burden of LC in MAOP in most countries declined. However, the proportion of countries with an increasing burden of disease was 24% for incidence rate, 38% for prevalence rate, 14% for mortality rate, and 16% for DALY rate (Figures 3A–C; Supplementary Figure S1; Supplementary Tables S2–S5). The countries with the highest increases in mortality and DALY rates were Guinea and Lesotho, with EAPC of 1.84 (1.64 to 2.05) and 2.04 (1.68 to 2.39), respectively. The top five countries with the increase in the rate of incidence, mortality, and DALY are Guinea, Chad, Lesotho, Ghana, and Sri Lanka (Figure 3F). It is worth noting that Sri Lanka had the highest prevalence and incidence increase, with EAPC of 3.02 (2.59 to 3.47) and 2.49 (2.05 to 2.94), respectively (Figure 3F).

Figure 3. The global burdens of LC in MAOP for both sexes across 204 countries from 1990 to 2021 and the ranks of the leading locations. The EAPC of incidence rate (A), death rate (B), and DALY rate (C) in 2021 worldwide and the top six countries. The trend of incidence number (D), incidence rate (E) of in the top 6 ranked countries for LC in MAOP from 1990 to 2021. The EAPC of DALY, incidence, and mortality in the top 6 ranked countries (F).
In 2021, China, India, the United States of America, Brazil, the Russian Federation, and France ranked in the top six cases of incidence and prevalence (Figure 3D and Supplementary Tables S2, S3). China and India had the highest cases of incidence, prevalence, mortality, and DALY due to their large populations. In 2021, China’s incidence and DALY case numbers were 32,588 (25,555 to 41,417) and 386,575 (301,814 to 490,830). The Indian numbers for this number are 26,512 (23,106 to 30,559) and 543,682 (472,859 to 627,443). With rates of 45.71 (33.98 to 61.77), 39.14 (30.41 to 51.66), and 38.15 (31.59 to 46.67) per 100,000 population, respectively, Monaco, Montenegro, and Cuba were the leading three nations by incidence rate, which also occupied the top three positions in terms of prevalence rate (Figure 3E and Supplementary Figures S2, S3). Monaco saw a notable decline during these 32 years, while Montenegro and Cuba saw increases (Figure 3F). The top four countries of death rate and DALYs rate were Cuba, Montenegro, Pakistan, and Seychelles (Supplementary Figures S4, S5). That year, India, China, Brazil, Pakistan, the United States of America, and the Russian Federation were the top six nations in mortality and DALY cases. These nations especially need to improve their efforts to prevent and control disease.
3.5 Patterns by sex and age group
Globally, mortality and DALY rate in the MAOP age groups demonstrated a general tendency toward decline from 1990 through 2021 (Figures 4C,D), while prevalence rates in the 85–89 and 90–94 age groups showed an increase, with the EAPC for these age groups 0.35 (0.3 to 0.4), and 0.38 (0.34 to 0.43), respectively (Figure 4A and Table 2).

Figure 4. EAPC of prevalence rate (A), incidence rate (B), death rate (C), and DALYs rate (D) in 9 age groups in Global and 5 SDI regions from 1990 to 2021. DALYs, disability-adjusted life-year; SDI, socio-demographic index; EAPC, estimated annual percentage changes.

Table 2. Prevalence cases and rates of LC in MAOP for both sexes in 1990 and 2021, and its temporal trends from 1990 to 2021.
Over the past 32 years, global reductions in the rate of prevalence, morbidity, mortality, and DALY were slightly better for men than for women (Table 1 and Figure 1). However, in 2021, globally, the rate of these indicators was considerably more severe for male than for female in different age groups and across SDI regions (Figures 1, 4 and Supplementary Figure S6). Specifically, the case of incidence and mortality in men worldwide were 140,040 (130,074 to 151,216) and 83,699 (77,817 to 90,455), respectively; the DALYs rate was 282.12 (261.85, 305.98) per 100,000, while the morbidity rate in women was 23,276 (20,094 to 26,886) (Table 1). The number and rate of prevalence, incidence, and DALYs were superior in males among MAOP before age 75; incidence and prevalence rates increased again between ages 80 and 94 years, with a peak prevalence rate of 28.05 (23.93, 30.47) per 100,000 persons in men aged 85 to 90 years. The prevalence rate and DALYs rates peaked in males aged 70–74 years at 132.21 (122.68, 142.76) per 100,000 population and 47.78 (39.55, 57.81) for every 100,000 people, while a secondary peak was observed in the prevalence rate of males aged 85–89 years at 119.37 (105.01, 128.56) per 100,000 population. The peak number of DALY was in the 60–64 age group, with 463,638 (429,177 to 504,950). The mortality rate exhibited an age-dependent increase, with the maximum case of deaths occurring aged 65–69 years, amounting to 16,841 (15,597 to 18,459) (Figure 5; Table 2; Supplementary Tables S6–S8).

Figure 5. The case and rate of prevalence (A), incidence (B), death (C), and DALYs (D) in 9 age groups globally in 2021. DALYs, disability-adjusted life-year; SDI, socio-demographic index.
Over the past 32 years, an age-related worsening trend has been observed in the global and five SDI regional burden of disease indicators, including prevalence, incidence, mortality, and DALY rates (Figure 4; Table 2; Supplementary Tables S6–S9). Mortality and DALY rates for the MAOP showed a decreasing trend in other regions, except for the 80+ age group in the Low SDI region (Figures 4C,D). However, regarding the incidence rate, increasing trends were observed among the following age groups: individuals aged 80+ years in low SDI regions, 85–94 years in low-middle SDI regions, 80–85 years in high-middle SDI regions, and 90–94 years in high SDI regions, while the remaining age groups show decreasing trends (Figure 4B and Supplementary Table S7). Regarding prevalence rate, several regions, and age groups still presented increasing trends, with the high-middle SDI region showing the best-improving trends in the rate of incidence and mortality and the low SDI region demonstrating the worst (Figure 4). The 55–59 age group showed the best improvement in prevalence, incidence, mortality, and DALY. It is noteworthy that the prevalence and incidence of females in high SDI regions remained persistently elevated relative to other SDI regions (Table 1).
3.6 Indicators correlation with SDI
Across most regions, SDI was positively correlated with prevalence and incidence rate, with correlation coefficients (R-values) of 0.604 and 0.038, respectively (Supplementary Figures S7, S8). Conversely, the rate of DALY and mortality exhibited an overall negative correlation with SDI, with correlation coefficients (R-values) of −0.122 and −0.154, respectively (Supplementary Figures S9, 10). Initially, the prevalence and incidence rates increased but decreased as the SDI level continued to rise. For mortality rate and DALYs rate, these indicators demonstrated slight fluctuations with the SDI until it reached 0.7, after which they began to decline. For instance, high SDI regions (e.g., Western Europe) had higher prevalence and morbidity rates, while their mortality and DALY rates are relatively low.
3.7 Risk factor analysis
Based on GBD 2021 data, our study evaluated the secondary-level contributable risk factors linked to LC in MAOP. Three risk factors were identified for laryngeal cancer mortality and DALY, including alcohol use, tobacco use, and occupational risk. Notably, tobacco was the leading risk factor, contributing substantially to laryngeal cancer mortality and DALYs in middle-aged and older adults, followed by alcohol consumption, and occupational risk factors were relatively minor. In 2021, the global numbers of deaths associated with alcohol use, tobacco use, and occupational hazards were 11,840 (6,542 to 17,002), 66,667 (59,522 to 73,780), and 5,885 (3,674 to 8,534), respectively, with corresponding death rates of 0.8 (0.44 to 1.14), 4.49 (4.01 to 4.97), and 0.4 (0.25 to 0.57) per 100,000 people. These factors contributed to DALY rates of 19.6 (10.79 to 27.66), 106.39 (95.51 to 117.35), and 9.06 (5.64 to 13.37) per 100,000 people, respectively.
Mortality and DALY rates associated with risk factors demonstrated systematic gender/age variations in LC in MAOP associated with risk factors. Men faced a much greater burden of disease related to tobacco and alcohol use than women, illustrating gender disparities (Figure 6). The number of DALYs and mortality related to alcohol use, tobacco use, and occupational risks increased initially and then decreased with age (Figures 7A,C). Conversely, mortality rates exhibited an increasing trend with age, peaking at 90–94 years for males (Figure 7B). DALY rates related to the three risk factors increased with age up to 74 years, after which a decreasing trend was observed (Figure 7D). Compared with 1990, the 2021 data revealed a decline in the proportion of tobacco-related deaths and DALYs among the male population aged over 90, while the proportion attributable to alcohol consumption and occupational risks increased (Figures 6A–D). A slight uptick was observed in the percentage of deaths and DALYs among females due to LC from occupational factors (Figures 6E–H). The three regions with the highest burden of DALY and deaths due to smoking were Eastern Europe, Central Europe, and East Asia (Figure 8).

Figure 6. Proportion of death and DALYs by age and gender attributable to risk factors in 1990 and 2021. Death (A) and DALYs (B) in males and 1990, Death (C) and DALYs (D) in males and 2021, Death (E) and DALYs (F) in females and 1990, Death (G) and DALYs (H) in female and 2021. DALYs, disability-adjusted life-year.

Figure 7. The case and rate of mortality and DALY associated with alcohol use, tobacco use, and occupational risk in 9 age groups globally in 2021. Death number (A), Death rate (B), DALY number (C), DALYs rate (D). DALYs, Disability-adjusted life-year.

Figure 8. Proportion of laryngeal cancer death (A) and DALYs (B) attributable to alcohol, tobacco, and occupational risk among 21 GBD regions. DALYs, disability-adjusted life-year.
Since the early 1990s, the mortality and DALY rates of LC in MAOP attributed to tobacco, alcohol use, and occupational factors demonstrated a downward trend globally and across five SDI regions, yet regional disparities persist. Males in the high SDI region showed the most remarkable drop in tobacco-related mortality and DALY rates in MAOP (Figures 9A,C). Conversely, the low-middle SDI region saw the least drop among males and became the region with the most excellent tobacco-related mortality rate among men in 2021. Among females, the middle SDI region showed t the most marked reduction, while the low SDI region experienced the least reduction (Figures 9B,D). Notably, although the burden of smoking-attributable LC in MAOP was lowest in males in the high SDI region, the high SDI regions consistently had higher rates of smoking-related mortality and DALY rate in females.

Figure 9. Temporal trends of mortality and DALY rates attributed to tobacco of LC in MAOP globally and in five SDI regions. Death rate attributable to tobacco in male (A) and female (B). DALYs rate attributable to tobacco in male (C) and female (D). MAOP, middle-aged and older populations; LC, laryngeal cancer; DALYs, disability-adjusted life-year.
Regarding the mortality rate and DALY of LC in MAOP caused by alcohol consumption, the high middle SDI region had consistently had the highest male alcohol-related mortality rate and DALY rate over the previous 32 years, despite showing the most significant decline in both male and female populations (Supplementary Figure S11). The mortality and DALY rate of alcohol-related LC in MAOP were increasing in both male and female populations in the low-middle SDI and Low SDI regions. Throughout the whole study period, females in high SDI areas had consistently increased alcohol-related death rates and DALYs for LC in MAOP (Supplementary Figures S11B,D).
Both male and female groups revealed a rising tendency in the low-middle SDI area in the occupational factors related to laryngeal cancer mortality and DALY rate, with the male group experiencing the most significant decrease in the high-middle SDI region and the female group in the high SDI region (Supplementary Figure S12). Therefore, it is essential to strengthen the prevention and control of smoking and alcohol consumption in low-middle SDI areas. At the same time, women should strengthen these measures in low SDI and high SDI areas.
4 Discussion
Based on the GBD database, this study provided a detailed analysis of the epidemiological patterns of LC in MAOP from 1990 to 2021. It explored the burden of disease in different age groups. The study’s findings indicated that while the overall prevalence, incidence, mortality, and DALY of LC in MAOP decreased globally, the burden of disease indicators increased in specific regions and age-specific groups. Tobacco was notably the primary risk factor, heavily influencing laryngeal cancer mortality and DALYs in MAOP. Moreover, the forecasted global disease burden for laryngeal cancer is expected to continue following the same trends that were observed from 1990 to 2021.
The non-linear correlation between SDI and health burden indicators carried significant policy ramifications. The present study identified a positive relationship between the prevalence of the condition and the SDI, with high SDI regions exhibiting a distinctive pattern characterized by a “high prevalence rate-low mortality rate.” For example, despite higher prevalence and morbidity, Western Europe and North America had significantly lower mortality and DALY rates than low and medium SDI regions. This occurrence might be related to higher alcohol and tobacco use per person (14). At the same time, these developed countries have better early detection systems (e.g., widespread laryngoscopy technology) (15) and targeted therapies (e.g., use of immunotherapy) (16), as well as better disease registration systems (17) and longer survival, which together improve survival and prolong the time patients live with the disease. In contrast, the negatively correlated relationship between mortality and SDI highlights the role of medical technology in improving prognosis. The low SDIs (e.g., sub-Saharan Africa) are limited by inadequate medical resources, with fewer than 50% of individuals diagnosed with cancer achieving survival rates beyond 3 years post-diagnosis (18). Its persistent barriers to cancer treatment, including surgery, radiotherapy, and chemotherapy, may contribute to the observed low cancer survival rates in sub-Saharan Africa. Prioritizing the strengthening of surgical care systems, expanding oncology-specific surgical education programs, and implementing multidisciplinary diagnostic support should constitute critical objectives for healthcare development in low SDI regions.
Central Europe, a high-middle SDI region, was a “hotspot” for the laryngeal cancer burden, with the highest prevalence and incidence rates in 2021. Additionally, it had the second-highest mortality and DALY rates, with a substantial proportion of DALYs and deaths attributed to smoking. Factors contributing to this include differences in unhealthy behaviors, including smoking and drinking (19), along with challenges and inequalities in cancer screening, diagnosis, treatment, and care (20). The Caribbean is also a country with a high burden of laryngeal cancer. Its incidence and prevalence rate were on an upward trend, and the rates of death and DALYs decreased the least. Resources for cancer control are limited due to the geographic distribution of islands in the Caribbean, lack of health system resources, and increased vulnerability to climatic hazards (21), and suboptimal cancer surveillance systems. A key contributor to the Caribbean’s high mortality is the prevalence of late-stage disease at diagnosis, resulting from the factors outlined earlier. These make cancer care in the Caribbean extremely challenging (22).
Epidemiological characterization of gender, age, and risk factors in middle-aged and older laryngeal cancer patients reveals complex interactions between biological differences, environmental exposures, and social behavior patterns. The findings of this study suggest that tobacco use plays a substantial role in the development of laryngeal carcinogenesis. The observed correlation appears mediated through nicotine’s regulatory effects on neoplastic cell dissemination dynamics in head and neck malignancies. This effect is associated with nicotine-altering CHRNA5 gene expression, which accelerates cancer cell metastasis (23). Statistical analyses demonstrated a significant correlation between smoking and the propensity for developing head and neck tumors, inclusive of a dose–response relationship with smoking intensity and longevity.
In instances where tobacco use is non-existent, the connection between alcohol intake and head and neck cancer is relatively modest. This association is evident primarily in the circumstances involving high-dose alcohol consumption and the occurrence of LC (24). This is consistent with the burden of smoking and alcohol in our study. Notably, individuals who chew tobacco exhibit an increased risk of laryngeal cancer (25). Especially in India and Sri Lanka, due to their unique habit of chewing tobacco, there is an urgent need to carry out targeted etiological research (26). The daily quantity of smoking is a significant indicator of disease risk. Middle-aged and older men bear a substantially greater share of smoking-related head and neck cancer cases compared to women. This phenomenon is closely related to gender differences in tobacco consumption. In 2019, the number of men aged 30 or over who smoke 20 or more cigarettes per day reached 231 million (32.3%), compared to 27.1 million women (18.5%) (27). In terms of the geographical distribution of smoking, studies have shown a significant increase in smoking prevalence in North Africa, the Middle East, and sub-Saharan Africa since 1990, with China accounting for over a third of global tobacco demand. In 2019, countries with the highest per capita tobacco consumption were primarily located in Europe, such as Montenegro, where the health burden of head and neck cancer among middle-aged and older individuals is relatively high. In low-income countries, over 2.8 billion people remain unprotected by smoke-free environment laws, and strict restrictions on tobacco marketing cover only 18% of the global population (28). Moreover, studies indicate that adult smoking rates are higher in high-income groups of wealthier countries, while the average adult smoking rate in middle-low-income countries is nearly twice that of low-income countries (LICs) (23% vs. 12%, respectively) (29).
Studies show a correlation between the age of the onset of tobacco-related deaths and the economic level of the country. Death caused by tobacco among the older population was higher in LMICs and upper-middle-income countries (UMICs), with proportions of 56 and 51%, respectively, as compared to 74% in LICs (29). The cost of smoking-related mortality as a proportion of the total socioeconomic cost of tobacco use was 64 and 83% in Chad and Sri Lanka, respectively, highlighting the severity of the tobacco-related burden in these two countries (29). Therefore, there is an urgent need to strengthen tobacco control measures in these countries. The study found that the proportion of smoking-attributable deaths increased with age, peaking at 60–64 (22.0% of smoking-attributable deaths), and then declined in older age groups (27), which is generally consistent with our findings. In the MAOP, the relatively long duration of smoking exposure, coupled with immune senescence, i.e., immune dysfunction that occurs with age in this population, affects tumor growth and the effectiveness of anti-tumor immune reactivity in older adults (30). Alcohol, a recognized risk factor for laryngeal cancers, has a heightened risk associated with the potency and frequency of alcohol consumption. In contrast, the effect of the duration of alcohol consumption is more complex. Ethanol may contribute to cancer development, including alterations or damage to DNA, proteins, and lipids acetaldehyde (7). In addition, specific industrial exposures [including asbestos (31) and benzene] have been associated with LC development. People working in specific industries, such as construction, metals, textiles, ceramics, logging, and food, are more likely to develop laryngeal cancer (32). This suggests the need to strengthen occupational health and safety legislation 2021. The increase in DALYs and deaths due to occupational exposures in women reflects the inadequate protection of female workers in manufacturing industries. In addition to these etiologies, there is evidence that infection with high-risk HPV, particularly HPV16, is a risk factor for LC (8).
Significant variations in the burden of LC worldwide require implementing differentiated control measures based on epidemiological characteristics and development status. It is recommended to actively implement a coordinated strategy for the prevention and control of laryngeal cancer at three levels: primary prevention should focus on tobacco and alcohol control and strengthen relevant legislation, especially for low-income, lower-middle-income and high-risk countries, and occupational safety and health; secondary prevention should promote low-cost screening technology in high-risk areas; and tertiary prevention should optimize treatment protocols for older adult patients. Capacity development for cancer control in LICs should be strengthened (5). The availability of surgical infrastructure, perioperative care, surgical expertise, team building, and comprehensive services should be improved (33).
For high-risk countries, tiered intervention strategies should be implemented, especially in low SDI nations (e.g., Chad), Central Europe and low-middle SDI, where primary prevention should be prioritized. Increasing taxes on tobacco products is the most cost-effective and valid intervention (29). Actively promote HPV vaccination; widespread vaccination has a possibility for a reduction in the incidence of HPV-positive LC (1). Improve the protection of occupationally exposed people, especially in Southeast Asia and low-middle SDI, where wearing N95 masks and laryngeal protection devices should be mandatory. Meanwhile, early detection of laryngeal cancer by peripheral blood exosome extraction improves the diagnostic capacity at the grassroots level (34); a range of advanced imaging techniques (e.g., multimodal optics, micro endoscopy, hyperspectral and infrared thermography) has relatively low-cost and non-invasive characteristics that can be used for rapid diagnosis (35). To improve precision screening and early intervention in populous countries, e.g., China, India, and other populous countries should implement 2-year laryngoscopic screening for high-risk men (≥40 years old, ≥20 pack-years of smoking), particularly in individuals under 75 and over 85 years of age.
Strengthening partnerships between Europe and developing nations will improve the research potential of these regions that are undergoing development (36). Integrating artificial intelligence, including machine learning, deep learning, and robotic surgery, into the overall management of laryngeal cancer can improve the quality of medical care (37). For older adult patients with early-stage disease, minimally invasive surgical or radiation treatments, such as transoral laser microsurgery, offer the potential for organ preservation (9). This revision underscores the significance of the multidisciplinary team (MDT), particularly for individuals with locally advanced cancer. In such cases, organ-sparing surgical procedures, Chemoradiotherapy, or radiotherapy alone can provide the possibility of laryngeal preservation while maintaining overall survival rates (38).
This study was limited by the retrospective characterization of Global Burden of Disease (GBD) data, with incomplete registration of deaths in some countries, particularly in areas with a low SDI. In addition, risk factors like tobacco and alcohol were not further delineated for measurement and stratification, and the etiologic role of HPV in LC was not separately assessed, potentially underestimating the impact of biological factors on epidemiologic patterns. In addition, laryngeal cancers were not examined by subtyping, e.g., glottic and subglottic types. In the future, there is a need to deepen molecular epidemiological studies of laryngeal cancer by combining multi-omics data (e.g., genomic and environmental exposure groups) and to develop region-specific risk prediction models to optimize screening strategies.
5 Conclusion
Globally, the burden of LC in MAOP experienced a consistent decline since 1990; however, the disease exhibits significant SDI stratification and substantial age and geographic variations. The disease burden is higher in low and low-middle SDI regions. Prevention and control strategies should be the primary focus before 75 and after 85 years old. The main risk factor for LC related deaths and DALY in MAOP is tobacco exposure, particularly in low-middle SDI regions where the burden is significant. Prevention and control of smoking and alcohol consumption in low SDI and high SDI areas should be strengthened among women. To effectively address the reduction of the global burden of LC in MAOP, it is imperative to comprehensively consider sociodemographic indices, geographic differences, gender disparities, age and risk factors, and the implementation of tailored intervention measures.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found at: https://ghdx.healthdata.org/gbd-2021.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
HX: Software, Data curation, Conceptualization, Formal analysis, Writing – original draft, Visualization. HW: Software, Methodology, Data curation, Validation, Writing – review & editing. HZ: Writing – review & editing, Data curation. SW: Writing – review & editing, Supervision, Project administration. BY: Writing – review & editing, Project administration, Funding acquisition.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This project was supported by the Clinical Skills Enhancement Programme for Young Healthcare Professionals of the Second Hospital of Dalian Medical University (No. 2024LCJSYL20), the In-hospital cultivation project of the Second Hospital of Dalian Medical University (No. 1) and Chinese Nursing Association Research Fund (ZHKYQ202420).
Acknowledgments
We would like to express our sincere gratitude to the research team members for their excellent collaboration and to honor the remarkable contribution of the Global Burden of Disease study.
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 authors declare that no Gen 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/fpubh.2025.1601029/full#supplementary-material
Footnotes
References
1. Johnson, DE, Burtness, B, Leemans, CR, Lui, VWY, Bauman, JE, and Grandis, JR. Head and neck squamous cell carcinoma. Nat Rev Dis Primers. (2020) 6:92. doi: 10.1038/s41572-020-00224-3
2. GBD 2019 Respiratory Tract Cancers Collaborators. Global, regional, and national burden of respiratory tract cancers and associated risk factors from 1990 to 2019: a systematic analysis for the global burden of disease study 2019. Lancet Respir Med. (2021) 9:1030–49. doi: 10.1016/S2213-2600(21)00164-8
3. Huang, A, Wu, XL, Song, J, Wang, YT, Yao, Y, Liu, Z, et al. Global trend and risk factors of the disease burden for pharynx and larynx cancers between 1990 and 2019: a systematic analysis of the global burden of disease study 2019. BMC Public Health. (2022) 22:2192. doi: 10.1186/s12889-022-14654-z
4. World Health Organization. (2024) Ageing and health. Available online at: https://www.who.int/news-room/fact-sheets/detail/ageing-and-health (accessed May 2, 2025)
5. Shawar, YR, Shiffman, J, and Spiegel, DA. Generation of political priority for global surgery: a qualitative policy analysis. Lancet Glob Health. (2015) 3:e487–95. doi: 10.1016/s2214-109x(15)00098-4
6. Cherny, NI, Trapani, D, Galotti, M, Saar, M, Bricalli, G, Roitberg, F, et al. ESMO global consortium study on the availability, out-of-pocket costs, and accessibility of cancer medicines: 2023 update. Ann Oncol. (2025) 36:247–62. doi: 10.1016/j.annonc.2024.12.005
7. Di Credico, G, Polesel, J, Dal Maso, L, Pauli, F, Torelli, N, Luce, D, et al. Alcohol drinking and head and neck cancer risk: the joint effect of intensity and duration. Br J Cancer. (2020) 123:1456–63. doi: 10.1038/s41416-020-01031-z
8. Kava, CM, Chaturvedi, AK, Senkomago, V, Mix, JM, Markowitz, LE, Kreimer, AR, et al. The role of HPV16 in oral cavity and laryngeal cancers in the United States. J Natl Cancer Inst. (2024) 1–6. doi: 10.1093/jnci/djae320
9. Machiels, JP, René Leemans, C, Golusinski, W, Grau, C, Licitra, L, and Gregoire, V. Squamous cell carcinoma of the oral cavity, larynx, oropharynx and hypopharynx: EHNS-ESMO-ESTRO clinical practice guidelines for diagnosis, treatment and follow-up. Ann Oncol. (2020) 31:1462–75. doi: 10.1016/j.annonc.2020.07.011
10. GBD 2021 Diseases and Injuries Collaborators. 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. (2024) 403:2133–61. doi: 10.1016/s0140-6736(24)00757-8
11. Liu, Q, Wang, H, Chen, Z, Xiong, J, Huang, Y, Zhang, S, et al. Global, regional, and national epidemiology of nasopharyngeal carcinoma in middle-aged and elderly patients from 1990 to 2021. Ageing Res Rev. (2025) 104:102613. doi: 10.1016/j.arr.2024.102613
12. Qin, N, Fan, Y, Yang, T, Yang, Z, and Fan, D. The burden of gastric Cancer and possible risk factors from 1990 to 2021, and projections until 2035: findings from the global burden of disease study 2021. Biomark Res. (2025) 13:5. doi: 10.1186/s40364-024-00720-8
13. Jani, CT, Kareff, SA, Morgenstern-Kaplan, D, Salazar, AS, Hanbury, G, Salciccioli, JD, et al. Evolving trends in lung cancer risk factors in the ten most populous countries: an analysis of data from the 2019 global burden of disease study. EClinicalMedicine. (2025) 79:103033. doi: 10.1016/j.eclinm.2024.103033
14. World Health Organization. (2024) Global status report on alcohol and health and treatment of substance use disorders. Available online at: https://www.who.int/publications/i/item/9789240096745 (accessed March 20, 2025)
15. Esmaeili, N, Davaris, N, Boese, A, Illanes, A, Navab, N, Friebe, M, et al. Contact endoscopy – narrow band imaging (CE-NBI) data set for laryngeal lesion assessment. Sci Data. (2023) 10:733. doi: 10.1038/s41597-023-02629-7
16. Mell, LK, Torres-Saavedra, PA, Wong, SJ, Kish, JA, Chang, SS, Jordan, RC, et al. Radiotherapy with cetuximab or durvalumab for locoregionally advanced head and neck cancer in patients with a contraindication to cisplatin (NRG-HN004): an open-label, multicentre, parallel-group, randomised, phase 2/3 trial. Lancet Oncol. (2024) 25:1576–88. doi: 10.1016/s1470-2045(24)00507-2
17. Shephard, EA, Parkinson, MA, and Hamilton, WT. Recognising laryngeal cancer in primary care: a large case-control study using electronic records. Br J Gen Pract. (2019) 69:e127–33. doi: 10.3399/bjgp19X700997
18. Joko-Fru, WY, Bardot, A, Bukirwa, P, Amidou, S, N’da, G, Woldetsadik, E, et al. Cancer survival in sub-Saharan Africa (SURVCAN-3): a population-based study. Lancet Glob Health. (2024) 12:e947–e59. doi: 10.1016/s2214-109x(24)00130-x
19. Rumgay, H, Shield, K, Charvat, H, Ferrari, P, Sornpaisarn, B, Obot, I, et al. Global burden of cancer in 2020 attributable to alcohol consumption: a population-based study. Lancet Oncol. (2021) 22:1071–80. doi: 10.1016/s1470-2045(21)00279-5
20. Lawler, M, Davies, L, Oberst, S, Oliver, K, Eggermont, A, Schmutz, A, et al. European Groundshot—addressing Europe's cancer research challenges: a lancet oncology commission. Lancet Oncol. (2023) 24:e11–56. doi: 10.1016/s1470-2045(22)00540-x
21. Ortiz, AP, Hospedales, CJ, Méndez-Lázaro, PA, Hamilton, WM, Rolle, LD, Shepherd, JM, et al. Protecting Caribbean patients diagnosed with cancer from compounding disasters. Lancet Oncol. (2024) 25:e217–24. doi: 10.1016/s1470-2045(24)00071-8
22. Spence, D, Dyer, R, Andall-Brereton, G, Barton, M, Stanway, S, Argentieri, MA, et al. Cancer control in the Caribbean island countries and territories: some progress but the journey continues. Lancet Oncol. (2019) 20:e503–21. doi: 10.1016/s1470-2045(19)30512-1
23. Feng, C, Mao, W, Yuan, C, Dong, P, and Liu, Y. Nicotine-induced CHRNA5 activation modulates CES1 expression, impacting head and neck squamous cell carcinoma recurrence and metastasis via MEK/ERK pathway. Cell Death Dis. (2024) 15:785. doi: 10.1038/s41419-024-07178-4
24. Hashibe, M, Brennan, P, Benhamou, S, Castellsague, X, Chen, C, Curado, MP, et al. Alcohol drinking in never users of tobacco, cigarette smoking in never drinkers, and the risk of head and neck cancer: pooled analysis in the international head and neck Cancer epidemiology consortium. J Natl Cancer Inst. (2007) 99:777–89. doi: 10.1093/jnci/djk179
25. Gil, GF, Anderson, JA, Aravkin, A, Bhangdia, K, Carr, S, Dai, X, et al. Health effects associated with chewing tobacco: a burden of proof study. Nat Commun. (2024) 15:1082. doi: 10.1038/s41467-024-45074-9
26. GBD 2019 Chewing Tobacco Collaborators. Spatial, temporal, and demographic patterns in prevalence of chewing tobacco use in 204 countries and territories, 1990-2019: a systematic analysis from the global burden of disease study 2019. Lancet Public Health. (2021) 6:e482–99. doi: 10.1016/s2468-2667(21)00065-7
27. GBD 2019 Tobacco Collaborators. 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. (2021) 397:2337–60. doi: 10.1016/s0140-6736(21)01169-7
28. Odo, DB, Ayo-Yusuf, O, Dinku, Y, Mekonnen, AG, and Maddox, R. Trends in the prevalence and factors associated with indoor smoking in 24 countries party to the WHO FCTC: implications for equitable policy implementation. BMJ Glob Health. (2025) 10:e017110. doi: 10.1136/bmjgh-2024-017110
29. Mann, N, Spencer, G, Hutchinson, B, Ngongo, C, Tarlton, D, Webb, D, et al. Interpreting results, impacts and implications from WHO FCTC tobacco control investment cases in 21 low-income and middle-income countries. Tob Control. (2024) 33:s17–26. doi: 10.1136/tc-2023-058337
30. Lian, J, Yue, Y, Yu, W, and Zhang, Y. Immunosenescence: a key player in cancer development. J Hematol Oncol. (2020) 13:151. doi: 10.1186/s13045-020-00986-z
31. Wronkiewicz, SK, Roggli, VL, Hinrichs, BH, Kendler, A, Butler, RA, Christensen, BC, et al. Chrysotile fibers in tissue adjacent to laryngeal squamous cell carcinoma in cases with a history of occupational asbestos exposure. Mod Pathol. (2020) 33:228–34. doi: 10.1038/s41379-019-0332-7
32. Godono, A, Dito, A, Martini, G, Picciaiola, MV, Di Lorenzo, A, Ciocan, C, et al. Occupational benzene exposure and risk of head and neck cancer: a systematic review and meta-analysis. Environ Res. (2024) 263:120033. doi: 10.1016/j.envres.2024.120033
33. Spence, D, Argentieri, MA, Andall-Brereton, G, Anderson, BO, Duggan, C, Bodkyn, C, et al. Advancing cancer care and prevention in the Caribbean: a survey of strategies for the region. Lancet Oncol. (2019) 20:e522–34. doi: 10.1016/s1470-2045(19)30516-9
34. Bocchetti, M, Luce, A, Iannarone, C, Pasquale, LS, Falco, M, Tammaro, C, et al. Exosomes multiplex profiling, a promising strategy for early diagnosis of laryngeal cancer. J Transl Med. (2024) 22:582. doi: 10.1186/s12967-024-05396-0
35. Mahmood, H, Shaban, M, Rajpoot, N, and Khurram, SA. Artificial intelligence-based methods in head and neck cancer diagnosis: an overview. Br J Cancer. (2021) 124:1934–40. doi: 10.1038/s41416-021-01386-x
36. Izarn, F, Henry, J, Besle, S, Ray-Coquard, I, Blay, JY, and Allignet, B. Globalization of clinical trials in oncology: a worldwide quantitative analysis. ESMO Open. (2025) 10:104086. doi: 10.1016/j.esmoop.2024.104086
37. Zhong, NN, Wang, HQ, Huang, XY, Li, ZZ, Cao, LM, Huo, FY, et al. Enhancing head and neck tumor management with artificial intelligence: integration and perspectives. Semin Cancer Biol. (2023) 95:52–74. doi: 10.1016/j.semcancer.2023.07.002
38. Forastiere, AA, Ismaila, N, Lewin, JS, Nathan, CA, Adelstein, DJ, Eisbruch, A, et al. Use of larynx-preservation strategies in the treatment of laryngeal Cancer: American Society of Clinical Oncology clinical practice guideline update. J Clin Oncol. (2018) 36:1143–69. doi: 10.1200/jco.2017.75.7385
Keywords: laryngeal cancer, global burden of disease, middle-aged, older adult, risk factors
Citation: Xiao H, Wang H, Zhang H, Wu S and Yu B (2025) Global, regional, and national burden of laryngeal cancer in middle-aged and older adults from 1990 to 2021: an analysis of age and sex differences and attributable risk factors. Front. Public Health. 13:1601029. doi: 10.3389/fpubh.2025.1601029
Edited by:
Monika Rucinska, University of Warmia and Mazury in Olsztyn, PolandReviewed by:
Julio Cesar Quispe Mamani, Universidad Nacional del Altiplano, PeruMohamed Hady Ashry, New Giza University, Egypt
Copyright © 2025 Xiao, Wang, Zhang, Wu and Yu. 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: Shanying Wu, c2hhbnlpbmd3dUAxMjYuY29t; Bo Yu, eXVibzE3NzA5ODcxMjYxQGRtdS5lZHUuY24=
†These authors have contributed equally to this work