Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Oncol., 06 January 2026

Sec. Cancer Epidemiology and Prevention

Volume 15 - 2025 | https://doi.org/10.3389/fonc.2025.1687320

Shifting burden of nasopharyngeal carcinoma: global patterns and forecasts to 2050 from the GBD 2021

Enhui Zhou,,&#x;Enhui Zhou1,2,3†Feifei Xu,,&#x;Feifei Xu1,2,3†Tianjiao Zhou,,&#x;Tianjiao Zhou1,2,3†Jingyu Zhang,,Jingyu Zhang1,2,3Fan Song,,Fan Song1,2,4Jianxiang Li*Jianxiang Li5*Hongliang Yi,,Hongliang Yi1,2,3Qingliang Wang*Qingliang Wang6*Weijun Huang,,*Weijun Huang1,2,3*
  • 1Department of Otorhinolaryngology Head and Neck Surgery, Shanghai Sixth People’s Hospital Affiliated to Shanghai Jiao Tong University School of Medicine, Shanghai, China
  • 2Shanghai Key Laboratory of Sleep Disordered Breathing, Shanghai, China
  • 3Otolaryngology Institute of Shanghai Jiao Tong University, Shanghai, China
  • 4Department of Otorhinolaryngology Head and Neck Surgery, Shanghai General Hospital Affiliated to Shanghai Jiao Tong University School of Medicine, Shanghai, China
  • 5Department of Otolaryngology, Shanghai Sixth People’s Hospital Fujian, Jinjiang, China
  • 6Department of Hepatobiliary-Pancreatic-Splenic Surgery, The Third Affiliated Hospital, Sun Yat-sen University, Guangzhou, China

Background: Nasopharyngeal carcinoma (NPC) exhibits pronounced geographical variation, with a high burden in specific regions. We assessed global, regional, and national trends in NPC burden from 1990 to 2021 and projected estimates to 2050.

Methods: We analyzed data from the Global Burden of Disease (GBD) 2021 study for 204 countries and territories across 21 regions. Age-standardized incidence rate (ASIR), age-standardized death rate (ASDR), and age-standardized disability-adjusted life-years (DALYs) rate were estimated, and temporal trends were assessed using estimated annual percentage change (EAPC) by socio-demographic index (SDI), sex, and age group. Future burden from 2022 to 2050 was projected using a Bayesian age–period–cohort model.

Results: In 2021, there were 118,878 incident cases, 75,359 deaths, and 2,490,191 DALYs due to NPC globally. The global ASIR was 1.38 per 100,000 population (EAPC −1.06), the ASDR was 0.87 (EAPC −2.17), and the age-standardized DALYs rate was 28.91 (EAPC −2.29). Incidence rate peaked at ages 65–69 years (4.23 per 100,000 population), with consistently higher rates in males than in females. East Asia had the highest regional ASIR, and Japan bore the highest age-specific disease burden. Projections indicate continued declines in apparent ASDR and DALYs rates globally to 2050, but rising ASIR in males, particularly in East Asia and high-middle SDI regions.

Conclusions: Over the past three decades, global mortality and disability burden from NPC have decreased, whereas incidence has increased in selected populations. The projected rise in ASIR among males in high-risk regions highlights the need for targeted interventions, equitable resource allocation, and sustained surveillance to mitigate future burden.

Background

Nasopharyngeal carcinoma (NPC) is an epithelial cancer of the head and neck that originates in the mucosal lining of the nasopharynx. In 2020, NPC caused approximately, 133,000 cases and 80,000 deaths globally, with over 75% occurring in Southeast Asia and southern China (1), and men affected three times more often than women (2). In addition to its clear association with Epstein–Barr virus (EBV) infection, environmental factors and genetic susceptibility are thought to contribute to a higher risk of developing NPC (3), including well-documented dietary exposures such as salted fish in endemic regions (4). According to the World Health Organization (WHO) pathological classification, NPC is divided into keratinizing squamous cell carcinoma and non-keratinizing carcinoma. While early-stage NPC can be effectively treated, the response to treatment for advanced NPC remains unsatisfactory, with an overall five-year survival rate of 38%−63%, largely due to the high rate of early misdiagnosis (over 70% are initially diagnosed at stage III or IV) (5). Therefore, a thorough understanding of cross-regional socio-demographic patterns of NPC can inform medical policies and enhance governmental awareness to promote early diagnosis and management. Previous studies on the global burden of NPC have been limited by sociological factors, such as population structure, sex disparities, and regional economic variations (6). Although several studies have examined national NPC burden and highlighted policy needs for high-risk regions and younger populations, the data in these studies were only updated to 2019 (7, 8). These constraints hinder a thorough understanding of the current NPC global impact.

The Global Burden of Disease Study (GBD) 2021 database provides a valuable resource for assessing the global burden of diseases and injuries (9), using standardized methodologies to ensure consistency and comparability across populations and temporal trends (10). More comprehensive global statistics may be necessary to further verify the true burden of NPC. Disparities in NPC incidence and mortality persist worldwide, highlighting the need for governments and policymakers to raise awareness about prevention and ensure judicious allocation of medical resources. This study used the latest GBD data from 1990 to 2021 to assess the global burden and temporal trends of NPC, offering insights for tailored approaches to alleviate the global impact of the disease.

Methods

Data source

Data relevant to the keywords were retrieved and downloaded from the GBD Results Tool on the Global Burden of Disease (GBD) database website (https://ghdx.healthdata.org/gbd-2021). We obtained data on the sex, age, incidence, deaths, and disability-adjusted life-years (DALYs) for nasopharyngeal carcinoma (NPC) from 1990 to 2021 globally, across 21 GBD regions and 204 countries and territories. As this study was a secondary analysis of a publicly available database, patient and public involvement was not applicable.

Measures and definitions

Age-standardized rates (ASRs) were calculated for per 100,000 population. Temporal trends were analyzed these trends using a log-linear regression model, with the natural logarithm of the ASR as the dependent variable and calendar year as the independent variable (11):

y=α+βx+ϵ

where y=ln(ASR) and x=calendar year. The estimated annual percentage change (EAPC) was calculated from the β coefficient using the formula (12):

EAPC=100×(eβ1)

An EAPC value greater than zero indicates an increasing trend in age-standardized rates, whereas an EAPC value less than zero indicates a decreasing trend. In addition to estimating EAPC using log-linear regression, we further applied joinpoint regression to statistically identify significant changes in temporal trends. In addition to estimating EAPC using log-linear regression, we further applied joinpoint regression to statistically identify significant changes in temporal trends. For each segment, the annual percent change (APC) was computed, and the average annual percent change (AAPC) summarized the overall trend across the full study period. A total of 204 countries or territories and 21 GBD regions were selected in this study. The socio-demographic index (SDI) was estimated as a geometric mean of 0–1 indices of lag distributed income per capita, mean education of those ages 15 and older, and total fertility rate in females under 25 years (6). The SDI was scaled from 0 to 1, with higher values indicating higher socioeconomic levels, according to which regions were categorized into five SDI levels (low, low-middle, middle, high-middle, and high). For age group, we categorized the age groups into 5-year intervals starting from <5 to ≥95 years. We defined the aging-related ratio as the ratio of disease burden in individuals aged over 65 years to that under 65 years. Future trends (2022-2050) in age-standardized incidence, mortality, and DALY rates were projected using a Bayesian age-period-cohort (BAPC) modeling framework (13). This approach integrates temporal patterns in age, calendar period, and birth cohort while accounting for uncertainty through Bayesian inference. The model was implemented using the BAPC and INLA packages in R, following forecasting strategies established in recent GBD analyses. 95% uncertainty intervals were visualized as shaded bands around the projected trajectories in the figures.

Statistical analysis

The statistical analyses were conducted using software R (version 4.4.1). Joinpoint analyses were conducted using the Joinpoint Regression Program (Version 5.4.0; National Cancer Institute), with differences between groups evaluated using the “Advanced analysis” – “Pairwise comparison procedure”. The data visualization was performed by using packages including ggplot2, heatmap, and so on. Data on incidence, deaths, DALYs, age-standardized incidence rate (ASIR), age-standardized death rate (ASDR), and age-standardized DALYs rate were presented with corresponding 95% uncertainty intervals (UIs). EAPCs were calculated with 95% confidence intervals (CIs). Values were displayed using bar charts and categorized by age, sex, and region subgroups. Annual trends were visualized using line graphs, and regional differences were depicted using heatmaps. Fitted curves were used to illustrate the associations and annual trends between SDI and disease burden indicators. Trends were described as “increasing” or “decreasing” when the slope was statistically significant; otherwise, they were described as “stable.” Statistical significance was defined as a two-sided p < 0.05.

Results

Global burden and temporal trends of NPC

Globally, the incident cases, deaths, and DALYs of NPC were 118,878 (95% UI: 104,836–135,884), 75,359 (95% UI: 67,515–83,706), and 2,490,191 (95% UI: 2,213,998–2,778,564), respectively, with an age-standardized incidence rate (ASIR), age-standardized death rate (ASDR), and age-standardized DALYs rate of 1.38 (1.22–1.58), 0.87 (0.78–0.97), and 28.91 (25.69–32.24), respectively (Table 1; Supplementary Figures S1-S3). However, the EAPCs of ASIR, ASDR, and age-standardized DALYs rate were −1.06 (-1.28 to -0.84), -2.17 (-2.37 to -1.98), and -2.29 (-2.49 to -2.09), respectively.

Table 1
www.frontiersin.org

Table 1. Global burden of nasopharyngeal carcinoma and temporal trends from 1990 to 2021 by gender and regions.

As shown in Table 1, the ASIR was highest in East Asia at 3.41 (2.77–4.19), whereas the highest ASDR was observed in Southeast Asia at 1.60 (1.41–1.81). The highest EAPCs of ASIR, ASDR, and age-standardized DALYs rate were observed in the Caribbean at 0.56 (0.44–0.68), 0.39 (0.29–0.50), and 0.43 (0.32–0.54), respectively. Figure 1 presents the global distribution of ASIR, ASDR, age-standardized DALYs rate, and their corresponding EAPCs. At the national level, Malaysia exhibited the highest ASIR at 6.09 (5.21–7.19), ASDR at 4.76 (4.10–5.62), and age-standardized DALYs rate at 158.23 (134.38–187.14). The highest EAPCs for ASIR, ASDR, and age-standardized DALYs rate were observed in Cape Verde at 5.07 (3.55–6.62), 4.75 (3.22–6.32), and 4.58 (3.02–6.15), respectively (Supplementary Table S1).

Figure 1
World maps show changes in health metrics across countries. Panels A and B depict age-standardized incidence rates (ASIR) and their percentage changes, highlighting increases in Northern Europe and parts of Africa. Panels C and D show age-standardized death rates (ASDR) and their changes, with decreases noted in some countries. Panels E and F illustrate age-standardized DALY rates and their changes, with notable variation across regions. Insets detail specific countries like Malaysia, Ghana, and Cape Verde, with numerical values indicating changes. Color gradients indicate intensity of changes.

Figure 1. Global distribution and temporal trends of nasopharyngeal carcinoma burden in 2021. (A) Age-standardized incidence rate (ASIR); (B) estimated annual percentage change (EAPC) of ASIR; (C) age-standardized death rate (ASDR); (D) EAPC of ASDR; (E) age-standardized DALY rate; (F) EAPC of age-standardized DALYs rate.

Burden and temporal trends of NPC by SDI regions

In 2021, the highest global ASIR (2.38, 1.94–2.93) and age-standardized DALYs rate (34.84, 29.05–41.64) for NPC were observed in high-middle SDI regions in 2021(Figure 2). Middle SDI regions had the highest ASDR (1.04, 0.92–1.19) but the lowest EAPCs for ASDR and age-standardized DALYs rate at −2.97 (−3.21 to −2.72) and −3.09 (−3.34 to −2.83), respectively. By contrast, high-SDI regions exhibited the lowest ASIR, ASDR, and age-standardized DALYs rate in 2021 at 0.71 (0.66-0.76), 0.35 (0.33–0.38), and 11.11 (10.34–11.92), respectively.

Figure 2
Nine line charts display trends from 1990 to 2021 for ASIR, ASDR, and age-standardized DALY rates globally, separated by male and female. Lines representdifferent SDI levels: high, high-middle, middle, low-middle, low, and global. Each chartshows overall declines across all SDI categories, with variations in steepness, especiallyprominent in high and high-middle SDI groups.

Figure 2. Trends in age-standardized rates of nasopharyngeal carcinoma by socio-demographic index (SDI) level and sex, 1990–2021.

Across the SDI spectrum, East Asia consistently presented higher ASIRs. Age-standardized DALYs rates showed a significant negative correlation with SDI (R = −0.16, p = 0.024) (Supplementary Figure S4).

Burden and temporal trends of NPC by gender and age

In 2021, the ASIR, ASDR, and age-standardized DALYs rate in males were 2.06 (1.76-2.42), 1.30 (1.14-1.47), and 42.26 (36.70-48.46), respectively, with EAPCs of -0.85 (-1.07 to -0.63), -2.00 (-2.18 to -1.81), and -2.10 (-2.29 to -1.91). In females, the ASIR, ASDR, and age-standardized DALYs rate were 0.74 (0.64–0.87), 0.48 (0.43–0.54), and 16.12 (14.34–18.38), respectively, with EAPCs of −1.60 (−1.91 to −1.39), −2.60 (−2.82 to −2.38), and −2.71 (−2.93 to −2.49) (Table 1).

Significant sex differences in AAPCs were observed across all measures within each SDI region (p < 0.001) (Supplementary Table S2; Supplementary Figures S5-S7). The lowest ASIR, ASDR, and age-standardized DALYs rates occurred in high-SDI regions for both sexes. However, high-middle SDI regions showed a clear upward incidence trend, with an APC of 1.64 (95% CI: 1.30–1.99, p < 0.01) in females during 2013–2021 and APC 2.74 (95% CI: 2.57–2.92, p < 0.01) in males during 2012–2021 (Figure 2; Supplementary Figures S5-S7; Supplementary Table S3). ASDR and age-standardized DALYs rates in females showed significant negative correlations with SDI (R = −0.23 and −0.25, respectively; both p < 0.001) (Supplementary Figure S8).

Globally, the highest incidence rate was observed in individuals aged 65–69 years at 4.23 (3.66–4.90), the highest death rate was in those aged 85–89 years at 4.17 (3.55–4.66), and the highest DALYs rate in individuals aged 55–59 years at 90.17 (80.03–100.91) (Supplementary Table S4). From 1990 to 2021, the global rates of NPC incidence, death, and DALYs rates generally declined across age groups (Supplementary Figure S9).

When stratified by sex, males showed the highest incidence rate at ages 65–69 years (6.61, 5.53–7.97), the highest death rate at ages 85–95 years (6.75, 5.88–7.64), and the highest DALYs rate at ages 55–59 years (134.92, 117.18–155.93). In females, the highest incidence rate was observed at ages 70–74 years (2.07, 1.78–2.40), the highest death rate at ages ≥95 years (3.05, 2.17–3.65), and the highest DALYs rate at ages 55–59 years (46.82, 41.11–53.65) (Figure 3).

Figure 3
Three side-by-side population pyramids show incidence, deaths, and DALYs rates by age and gender. Each graph displays male data in blue and female data in red, with corresponding age groups from under five to over 95 years. Rates are shown as percentage and number per 100,000 individuals, featuring shaded areas and error bars for precision.

Figure 3. Number and rate for incidence, deaths, and disability-adjusted life years (DALYs) of nasopharyngeal carcinoma across different gender and age groups in 2021.

From 1990 to 2021, the ratio of male-to-female ratio of NPC incidence increased significantly among individuals aged 50–54 years (p < 0.05), while the ratios for deaths and DALYs increased significantly among those aged 25–54 years, except for ages 40–44 years (p < 0.05) (Supplementary Figure S10). Furthermore, Japan exhibited the highest aging-related ratios for incidence (3.94), deaths (5.11 times), and DALYs (4.44 times) between 1990 and 2021 (Supplementary Figure S11).

Predictions of NPC from 2022 to 2050

We projected ASIR, ASDR, and age-standardized DALYs rate for NPC over the next 30 years by using a Bayesian age–period–cohort (BAPC) model. The results indicate that the global ASIR is expected to reach 1.53 per 100,000 population by 2050, representing an annual increase of 0.43% (Figure 4). This increasing trend is primarily driven by males, whose ASIR is projected to increase by male ASIR, which increasing 0.6% annually, whereas the annual change in ASIR among females is projected to be −0.1% in female (Supplementary Figure S12).

Figure 4
Three line charts depict global trends from 1990 to 2050 for age-standardized rates. The first chart shows incidence rates slightly increasing with an annual change of 0.43 percent after 2022. The second chart displays death rates with a slight decrease, having an annual change of -0.14 percent. The third chart indicates a significant decline in DALYs (Disability-Adjusted Life Years) with an annual change of -19.43 percent post-2022. Each chart includes projections beyond 2022, illustrated with a shaded confidence interval.

Figure 4. Observed and projected trends in nasopharyngeal carcinoma burden, 1990–2050. Green dots indicate global estimates from 1990 to 2021; the dotted line marks the start of projections (2022); shaded areas represent the 95% uncertainty intervals (UI).

In contrast, the ASDR and age-standardized DALYs rates are estimated to decline to 0.75 and 23.09 per 100,000 population, respectively, by 2050, with annual decreases of −0.14% and −1.94%, respectively.

Discussion

Using the latest GBD 1990–2021 data, this study offers a comprehensive global assessment of NPC burden across SDI, age, and sex, identifies key epidemiological shifts, and projects trends through 2050. Although the ASIR of NPC in 2021 appears to have decreased significantly compared with 1990, previous studies consistently reported similar declines across children, young adults, middle-aged, and older populations (7, 14). Our results indicate that the ASIR trend is predominantly driven by males, which may contribute to a continued global increase in the future. Consistent with the previous studies, the ASDR and age-standardized DALYs rate of NPC have declined over the past decades (2, 15). Furthermore, the present research projected the continuation of these decreasing trends in the future. Overall, the significant decline is likely attributed to advancements in diagnostic and therapeutic technologies. In part, although the EBV infection remains widespread, advances in detection technology have substantially improved the sensitivity and predictive value of tests for anti-EBV antibodies and plasma or nasopharyngeal EBV DNA, including through fragment size and methylation-based analysis, even in asymptomatic period (16, 17). Moreover, substantial investment in EBV-based screening and heightened awareness within health systems have enabled earlier detection of NPC (17). Additionally, the oral–nasopharyngeal microenvironment plays an increasingly recognized role in NPC development (18). Alterations in the oral–nasopharyngeal microbiome may influence local immune responses and facilitate EBV reactivation, thereby elevating NPC risk (19). Conversely, improvements in oral health and microbial balance may help attenuate these pathways, providing a plausible mechanism linking environmental changes to reductions in NPC incidence (20). High-SDI regions tend to have a lower burden of oral disorders, potentially contributing to reduced NPC incidence (21). At the 2022 World Health Assembly, member states adopted a global strategy targeting universal oral health coverage for all populations by 2030 (22). Together, these efforts create a critical opportunity to strengthen care in a historically overlooked field and support broader NPC prevention (23). Furthermore, advances in diagnostic imaging technology, accompanied by intensity-modulated radiotherapy (IMRT) and optimization of chemotherapy regimens, have significantly improved early diagnosis and survival rates for NPC (24, 25).

Our findings highlight enduring sex disparities in NPC, with males experiencing a disproportionately higher morbidity globally and in future estimates. Furthermore, the decreasing burden trend in males was far slower than in females (26). After adjustment for major risk factors, previous studies have consistently documented sex disparities in cancer incidence and survival (27). These differences may reflect sex-specific risk exposures, including older age at diagnosis, more advanced stage at presentation, more aggressive tumor biology, and suboptimal treatment (28, 29). Risky behaviors such as tobacco and alcohol use (30), which are more prevalent among males, may further contribute to these disparities. Moreover, sex hormones and sex chromosomes are potential mechanisms driving female survival advantages (31).

A comprehensive analysis of 2021 GBD-observed NPC still showed regional variations (8). Significant differences in the cancer spectrum among regions can be attributed to living habits and unhealthy diets. Findings from the International Tobacco Control (ITC) Project indicate that Malaysia has a high level of cigarette tax avoidance or evasion (39.7%) (32), which is higher than in other Asian countries and may in turn contribute to higher smoking rates (37.5% in 15–64 years males) (33). Malaysia has introduced measures to reduce tobacco use and secondhand smoke exposure in line with the Framework Convention on Tobacco Control, ratified in 2005 (34). Greenland faces comparable challenges, with harmful alcohol use and rising obesity driven by dietary shifts from marine to domestic animal products and decreased physical activity (35, 36). The government has addressed these public health issues through its preventive programs, Inuuneritta I and Inuuneritta II (37).

Moreover, despite overall declines in morbidity and mortality, NPC incidence and mortality remain disproportionately high in East Asia. Regions with high ASIR and high-ASDR regions are therefore more likely to prioritize disease research and attract greater governmental and health-system attention, facilitating public awareness and advancement of screening technologies. Similar measures may be needed in the Caribbean, which is projected to experience an increasing NPC burden by 2021. In this region, much of the cancer burden is driven by tobacco (including SHS) and alcohol use, and a steady increase in metabolic diseases (30, 38). To address these challenges, the first edition of the Latin America and Caribbean Code Against Cancer was developed in collaboration with the Pan-American Health Organization. This code includes evidence-based, individual-level cancer prevention recommendations targeted at the general population (39). Cape Verde has exhibited a continuous increasing trend in NPC prevalence, however, this phenomenon has not attracted widespread attention. As in many African countries, the absence of a national screening program, insufficient attention, and lack of dedicated cancer treatment facilities (unavailability and high cost of high sensitivity screening tests), make cancer screening far more difficult than it should be, contributing to the region’s high mortality rates (40, 41). With changes in EBV exposure patterns, including improved oral and upper respiratory hygiene, lower childhood infection risk (42), and overall improvements in living and socioeconomic conditions (43), several historically high-incidence regions, such as Singapore, have experienced declines in NPC rates. Beyond smoking and alcohol, long-term dietary and environmental improvements—such as reduced consumption of Cantonese-style salted fish following widespread refrigeration and better air quality—may also have contributed to declining NPC incidence among younger cohorts in endemic regions (44, 45).

Although early diagnosis, behavioral change, and technological advances have reduced the burden of NPC, recurrence or metastasis (R/M) and population ageing remain major barriers to cure, with 20–30% of patients with advanced disease experiencing treatment failure due to R/M (46). Immunotherapy, which enhances antitumor immunity by modulating immune effector cells, has shown safety and efficacy in recent clinical studies (47). Furthermore, a higher age-related ratio indicates that older adults may have a more significant impact on the incidence and mortality rate of NPC. This age-driven burden has become more prominent in Japan and Singapore. Japan has undergone pronounced population ageing over recent decades, reflected in a substantial rise in life expectancy—from 79.4 years in 1990 to 85.2 years in 2021 (48). Singapore faces a similar ageing challenge (49). Population aging may increase cancer-related mortality and complicate treatment because of coexisting conditions, posing substantial challenges to both medical and community care systems. Preparation for an aging—and increasingly super-aged—society is therefore crucial, not only from healthcare and resource-planning perspectives but also from broader socioeconomic perspectives.

The present study uses the latest datasets to update and project evidence on the global NPC burden of NPC from 1990 to 2050. These findings provide valuable insights for estimating future disease trends and informing strategies aimed at improving health outcomes.

This study also has limitations. First, GBD estimates rely on modeling assumptions and data harmonization, which may introduce uncertainty (50). Second, cancer registry coverage remains incomplete in several low-SDI regions, potentially affecting estimate accuracy despite correction procedures (6). Third, long-term projections based on the BAPC inherent uncertainty, as future demographic and epidemiological changes may diverge from historical patterns. These factors should be considered when interpreting the results.

Conclusion

Although global incidence, mortality, and disability burden from nasopharyngeal carcinoma have declined since 1990, the past decade has been marked by a rising ASIR, particularly among males. East Asia, high-middle socio-demographic index regions, and aging populations continue to experience disproportionately high burden. Strengthening standardized management strategies and implementing targeted prevention efforts in these settings are essential to mitigate future increases in NPC burden.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Author contributions

EZ: Formal analysis, Writing – original draft, Visualization. FX: Writing – original draft, Methodology. TZ: Writing – original draft, Data curation, Methodology. JZ: Data curation, Writing – original draft. FS: Data curation, Writing – original draft. JL: Supervision, Writing – review & editing. HY: Resources, Writing – original draft, Supervision. QW: Writing – review & editing. WH: Supervision, Conceptualization, Writing – review & editing, Writing – original draft.

Funding

The author(s) declared financial support was received for this work and/or its publication. This study was supported by grants from the National STI2030-Major Projects of China (Grant No. 2021ZD0201900), Fundamental Research Funds for the Central Universities (Grant No. YG2025QNA44 and YG2025QNB51), National Natural Science Foundation of China (Grant Nos. 81970869, 82171125, and 82371131), Hospital Level Scientific Research Fund Program of Shanghai Sixth People’s Hospital (Grant No. ynts202404).

Acknowledgments

We sincerely appreciate the exceptional cooperation of our study team members and gratefully acknowledge the significant contributions of the Global Burden of Disease Study.

Conflict of interest

The authors declared that this work 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 Generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.1687320/full#supplementary-material

Supplementary Table 1 | Global burden of nasopharyngeal carcinoma and temporal trends from 1990 to 2021 across 204 countries or territories.

Supplementary Table 2 | Average annual percent change (AAPC) estimates derived from joinpoint regression using group comparison across the five socio-demographic index (SDI) regions, stratified by sex.

Supplementary Table 3 | Age-standardized incidence, mortality, and disability-adjusted life years (DALYs) rates across SDI regions, 1990–2023, with annual percentage changes (APC) by sex.

Supplementary Table 4 | Peak age groups and corresponding incidence, deaths, and disability-adjusted life years (DALYs) rates in 2021 across major regions.

Supplementary Figure 1 | Annual percent change (APC) and average annual percent change (AAPC) in age-standardized incidence rate (ASIR) at the global level and across five SDI regions, both sexes combined, 1990–2021. ASIR, age-standardized incidence rate; APC, annual percent change; AAPC, average annual percent change. *, P value <0.01.

Supplementary Figure 2 | Annual percent change (APC) and average annual percent change (AAPC) in age-standardized death rate (ASDR) for the global level and five SDI regions, both sexes combined, 1990–2021. ASDR, age-standardized death rate; APC, annual percent change; AAPC, average annual percent change. *, P value <0.01.

Supplementary Figure 3 | Annual percent change (APC) and average annual percent change (AAPC) in age-standardized DALYs rate for the global level and five SDI regions, both sexes combined, 1990–2021. DALYs, Disability-adjusted life year; APC, annual percent change; AAPC, average annual percent change. *, P value <0.01.

Supplementary Figure 4 | Disease burden of nasopharyngeal carcinoma by socio-demographic index (SDI) across 204 countries or territories in 2021. ASIR, age-standardized incidence rate; ASDR, age-standardized death rate; DALYs, disability-adjusted life years; SDI, Socio-demographic Index.

Supplementary Figure 5 | Annual percent change (APC) and average annual percent change (AAPC) in age-standardized incidence rate across five SDI regions, by sex, 1990–2021. Joinpoint regression was conducted using a sex-stratified group comparison approach, with male and female trends modelled concurrently. * P value of APC <0.01; ** P value of AAPC between female and males <0.01. APC, annual percent change; AAPC, average annual percent change

Supplementary Figure 6 | Annual percent change (APC) and average annual percent change (AAPC) in age-standardized death rate across five SDI regions, by sex, 1990–2021. Joinpoint regression was conducted using a sex-stratified group comparison approach, with male and female trends modelled concurrently. * P value of APC <0.01; ** P value of AAPC between female and males <0.01. APC, annual percent change; AAPC, average annual percent change

Supplementary Figure 7 | Annual percent change (APC) and average annual percent change (AAPC) in age-standardized DALYs rate across five SDI regions, by sex, 1990–2021. Joinpoint regression was conducted using a sex-stratified group comparison approach, with male and female trends modelled concurrently. * P value of APC <0.01; ** P value of AAPC between female and males <0.01. DALYs, Disability-adjusted life year; APC, annual percent change; AAPC, average annual percent change.

Supplementary Figure 8 | Disease burden of nasopharyngeal carcinoma by different gender across 204 countries or territories and socio-demographic index level in 2021. ASIR, age-standardized incidence rate; ASDR, age-standardized death rate; DALYs, disability-adjusted life years; SDI, Socio-demographic Index.

Supplementary Figure 9 | Number and rate for incidence, death, and disability-adjusted life years (DALYs) of nasopharyngeal carcinoma across different age groups in 2021. DALYs, disability-adjusted life years

Supplementary Figure 10 | Rate ratio for male to female in incidence, death, and disability-adjusted life years (DALYs) of nasopharyngeal carcinoma across different age groups in 1990, 2000, 2010, and 2021. DALYs, disability-adjusted life years.

Supplementary Figure 11 | Aging-related ratio for incidence, death, and disability-adjusted life years (DALYs) of nasopharyngeal carcinoma by gender in 204 countries or territories. DALYs, disability-adjusted life years.

Supplementary Figure 12 | The temporal and predicted trends for disease burden of nasopharyngeal carcinoma by different gender from 1990 to 2050. Blue spots represent the male levels from 1990 to 2021. Red spots represent the female levels from 1990 to 2021. The dotted line represents the beginning of the prediction (from 2022). Shading indicates the upper and lower limits of the 95% UI. DALYs, disability-adjusted life years

Abbreviations

NPC, nasopharyngeal cancer; EBV, Epstein-Barr virus; WHO, World Health Organization; GBD, Global Burden of Diseases; DALYs, disability-adjusted life-years; EAPC, estimated annual percentage change; SDI, socio-demographic index; BAPC, Bayesian age-period-cohort; ASR, age-standardized rate; ASIR, age-standardized incidence rate; ASDR, age-standardized death rate; UI, uncertainty interval; CI, confidence interval; IMRT, intensity-modulated radiotherapy; SHS, secondhand smoke exposure; R/M, recurrence and/or metastasis; APC, annual percent change; AAPC; average annual percent change.

References

1. Wong KC, Hui EP, Lo K-W, Lam WK, Johnson D, Li L, et al. Nasopharyngeal carcinoma: an evolving paradigm. Nat Rev Clin Oncol. (2021) 18:679–95. doi: 10.1038/s41571-021-00524-x

PubMed Abstract | Crossref Full Text | Google Scholar

2. Zhang Y, Rumgay H, Li M, Cao S, and Chen W. Nasopharyngeal cancer incidence and mortality in 185 countries in 2020 and the projected burden in 2040: population-based global epidemiological profiling. JMIR Public Health Surveill. (2023) 9:e49968. doi: 10.2196/49968

PubMed Abstract | Crossref Full Text | Google Scholar

3. Renaud S, Lefebvre A, Mordon S, Moralès O, and Delhem N. Novel therapies boosting T cell immunity in epstein barr virus-associated nasopharyngeal carcinoma. Int J Mol Sci. (2020) 21:4292. doi: 10.3390/ijms21124292

PubMed Abstract | Crossref Full Text | Google Scholar

4. Yu MC and Henderson BE. Intake of Cantonese-style salted fish as a cause of nasopharyngeal carcinoma. IARC Sci Publ. (1987) 84:547–9.

PubMed Abstract | Google Scholar

5. Zhang Q, Wu G, Yang Q, Dai G, Li T, Chen P, et al. Survival rate prediction of nasopharyngeal carcinoma patients based on MRI and gene expression using a deep neural network. Cancer Sci. (2023) 114:1596–605. doi: 10.1111/cas.15704

PubMed Abstract | Crossref Full Text | Google Scholar

6. GBD 2019 Cancer Risk Factors Collaborators. The global burden of cancer attributable to risk factors, 2010-19: a systematic analysis for the Global Burden of Disease Study 2019. Lancet. (2022) 400:563–91. doi: 10.1016/S0140-6736(22)01438-6

PubMed Abstract | Crossref Full Text | Google Scholar

7. Lei S, Chen L, Ji P, Li K, Li Q, Huang C, et al. Global burdens of nasopharyngeal carcinoma in children and young adults and predictions to 2040. Oral Oncol. (2024) 155:106891. doi: 10.1016/j.oraloncology.2024.106891

PubMed Abstract | Crossref Full Text | Google Scholar

8. Wei X, Chen B, Wang Z, Zhao P, and Duan X. Nasopharyngeal cancer risk assessment by country or region worldwide from 1990 to 2019. BMC Public Health. (2024) 24:1931. doi: 10.1186/s12889-024-19228-9

PubMed Abstract | Crossref Full Text | Google Scholar

9. Shen Y, Zhou T, Zou W, Zhang J, Yan S, Ye H, et al. Global, regional, and national burden of hearing loss from 1990 to 2021: findings from the 2021 global burden of disease study. Ann Med. (2025) 57:2527367. doi: 10.1080/07853890.2025.2527367

PubMed Abstract | Crossref Full Text | Google Scholar

10. Zhou T, Huang W, Wang X, Zhang J, Zhou E, Tu Y, et al. Global burden of head and neck cancers from 1990 to 2019. iScience. (2024) 27:109282. doi: 10.1016/j.isci.2024.109282

PubMed Abstract | Crossref Full Text | Google Scholar

11. Xu C, Zhou E, Shen Y, Xiao L, Huang W, Zhou T, et al. Assessment and future projection of brain and central nervous system cancer burden using a modified quality care index: evidence from the Global Burden of Disease 2021. Eur J Oncol Nurs. (2025) 79:102997. doi: 10.1016/j.ejon.2025.102997

PubMed Abstract | Crossref Full Text | Google Scholar

12. Zhou T, Wang X, Zhang J, Zhou E, Xu C, Shen Y, et al. Global burden of thyroid cancer from 1990 to 2021: a systematic analysis from the Global Burden of Disease Study 2021. J Hematol Oncol. (2024) 17:74. doi: 10.1186/s13045-024-01593-y

PubMed Abstract | Crossref Full Text | Google Scholar

13. Wang X, Shang A, Chen H, Li J, Jiang Y, Wang L, et al. Global, regional, and national disease burden of multidrug-resistant tuberculosis without extensive drug resistance, 1990-2021: Findings from the Global Burden of Disease Study 2021. Drug Resist Update. (2025) 82:101265. doi: 10.1016/j.drup.2025.101265

PubMed Abstract | Crossref Full Text | Google Scholar

14. 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

PubMed Abstract | Crossref Full Text | Google Scholar

15. Wu T, Miao W, Qukuerhan A, Alimu N, Feng J, Wang C, et al. Global, regional, and national burden of nasopharyngeal carcinoma from 1990 to 2021. Laryngoscope. (2025) 135:1409–18. doi: 10.1002/lary.31939

PubMed Abstract | Crossref Full Text | Google Scholar

16. Kwong DL, Pow EH, Sham JS, McMillan AS, Leung LH, Leung WK, et al. Intensity-modulated radiotherapy for early-stage nasopharyngeal carcinoma: a prospective study on disease control and preservation of salivary function. Cancer. (2004) 101:1584–93. doi: 10.1002/cncr.20552

PubMed Abstract | Crossref Full Text | Google Scholar

17. Miller JA, Le Q-T, Pinsky BA, and Wang H. Cost-effectiveness of nasopharyngeal carcinoma screening with epstein-barr virus polymerase chain reaction or serology in high-incidence populations worldwide. J Natl Cancer Inst. (2021) 113:852–62. doi: 10.1093/jnci/djaa198

PubMed Abstract | Crossref Full Text | Google Scholar

18. Hao Y, Zeng Z, Peng X, Ai P, Han Q, Ren B, et al. The human oral - nasopharynx microbiome as a risk screening tool for nasopharyngeal carcinoma. Front Cell Infect Microbiol. (2022) 12:1013920. doi: 10.3389/fcimb.2022.1013920

PubMed Abstract | Crossref Full Text | Google Scholar

19. Liao Y, Zhang J-B, Lu L-X, Jia Y-J, Zheng M-Q, Debelius JW, et al. Oral microbiota alteration and roles in epstein-barr virus reactivation in nasopharyngeal carcinoma. Microbiol Spectr. (2023) 11:e0344822. doi: 10.1128/spectrum.03448-22

PubMed Abstract | Crossref Full Text | Google Scholar

20. Chen Y, Li W, Chang ET, Debelius JW, Manoharan L, Zheng Y, et al. Oral fungal profiling and risk of nasopharyngeal carcinoma: a population-based case-control study. EBioMedicine. (2023) 96:104813. doi: 10.1016/j.ebiom.2023.104813

PubMed Abstract | Crossref Full Text | Google Scholar

21. Bernabe E, Marcenes W, Hernandez CR, Bailey J, Abreu LG, Alipour V, et al. Global, regional, and national levels and trends in burden of oral conditions from 1990 to 2017: A systematic analysis for the global burden of disease 2017 study. J Dent Res. (2020) 99:362–73. doi: 10.1177/0022034520908533

PubMed Abstract | Crossref Full Text | Google Scholar

22. Da Cunha AR, Compton K, Xu R, Mishra R, Drangsholt MT, Antunes JL, et al. The global, regional, and national burden of adult lip, oral, and pharyngeal cancer in 204 countries and territories: A systematic analysis for the global burden of disease study 2019. JAMA Oncol. (2023) 9:1401–16. doi: 10.1001/jamaoncol.2023.2960

PubMed Abstract | Crossref Full Text | Google Scholar

23. Benzian H, Guarnizo-Herreño CC, Kearns C, Muriithi MW, and Watt RG. The WHO global strategy for oral health: an opportunity for bold action. Lancet. (2021) 398:192–4. doi: 10.1016/S0140-6736(21)01404-5

PubMed Abstract | Crossref Full Text | Google Scholar

24. Huang H, Miao J, Xiao X, Hu J, Zhang G, Peng Y, et al. Impact on xerostomia for nasopharyngeal carcinoma patients treated with superficial parotid lobe-sparing intensity-modulated radiation therapy (SPLS-IMRT): A prospective phase II randomized controlled study. Radiother Oncol. (2022) 175:1–9. doi: 10.1016/j.radonc.2022.07.006

PubMed Abstract | Crossref Full Text | Google Scholar

25. Chen Y-P, Chan AT, Le Q-T, Blanchard P, Sun Y, and Ma J. Nasopharyngeal carcinoma. Lancet. (2019) 394:64–80. doi: 10.1016/S0140-6736(19)30956-0

PubMed Abstract | Crossref Full Text | Google Scholar

26. Bray F, Laversanne M, Sung H, Ferlay J, Siegel RL, Soerjomataram I, et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2024) 74:229–63. doi: 10.3322/caac.21834

PubMed Abstract | Crossref Full Text | Google Scholar

27. Radkiewicz C, Johansson AL, Dickman PW, Lambe M, and Edgren G. Sex differences in cancer risk and survival: A Swedish cohort study. Eur J Cancer. (2017) 84:130–40. doi: 10.1016/j.ejca.2017.07.013

PubMed Abstract | Crossref Full Text | Google Scholar

28. Li W-Z, Lv S-H, Liu G-Y, Liang H, Xia W-X, and Xiang Y-Q. Age-dependent changes of gender disparities in nasopharyngeal carcinoma survival. Biol Sex Differ. (2021) 12:18. doi: 10.1186/s13293-021-00361-8

PubMed Abstract | Crossref Full Text | Google Scholar

29. Micheli A, Ciampichini R, Oberaigner W, Ciccolallo L, de VE, Izarzugaza I, et al. The advantage of women in cancer survival: an analysis of EUROCARE-4 data. Eur J Cancer. (2009) 45:1017–27. doi: 10.1016/j.ejca.2008.11.008

PubMed Abstract | Crossref Full Text | Google Scholar

30. Reynales-Shigematsu LM, Barnoya J, Cavalcante T, Aburto TC, Romieu I, Stern MC, et al. Latin America and the Caribbean Code Against Cancer 1st edition: Tobacco and nicotine-related products, secondhand smoke, and alcohol and cancer. Cancer Epidemiol. (2023) 86 Suppl 1:102413. doi: 10.1016/j.canep.2023.102413

PubMed Abstract | Crossref Full Text | Google Scholar

31. Clocchiatti A, Cora E, Zhang Y, and Dotto GP. Sexual dimorphism in cancer. Nat Rev Cancer. (2016) 16:330–9. doi: 10.1038/nrc.2016.30

PubMed Abstract | Crossref Full Text | Google Scholar

32. Guindon GE, Driezen P, Chaloupka FJ, and Fong GT. Cigarette tax avoidance and evasion: findings from the International Tobacco Control Policy Evaluation (ITC) Project. Tob Control. (2014) 23:i13–22. doi: 10.1136/tobaccocontrol-2013-051074

PubMed Abstract | Crossref Full Text | Google Scholar

33. Tan QY, Zomer E, Owen AJ, Chin KL, and Liew D. Impact of tobacco use on health and work productivity in Malaysia. Tob Control. (2020) 29:111–7. doi: 10.1136/tobaccocontrol-2018-054677

PubMed Abstract | Crossref Full Text | Google Scholar

34. Abdul Mutalib RN, Abd Rani NL, Zulkifli A, Abd Latif NH, Dobson R, Engku Ibrahim TA, et al. Knowledge, beliefs, and behaviors related to secondhand smoke and smoking in the home: A qualitative study with men in Malaysia. Nicotine Tob Res. (2023) 25:821–7. doi: 10.1093/ntr/ntac239

PubMed Abstract | Crossref Full Text | Google Scholar

35. Bjerregaard P, Olesen I, and Larsen CV. Association of food insecurity with dietary patterns and expenditure on food, alcohol and tobacco amongst indigenous Inuit in Greenland: results from a population health survey. BMC Public Health. (2021) 21:1094. doi: 10.1186/s12889-021-11123-x

PubMed Abstract | Crossref Full Text | Google Scholar

36. Andersen S, Rex KF, Noahsen P, Sørensen HC, Larsen NH, Mulvad G, et al. Forty-five year trends in overweight and obesity in an indigenous arctic Inuit Society in transition and spatiotemporal trends. Am J Hum Biol. (2014) 26:511–7. doi: 10.1002/ajhb.22556

PubMed Abstract | Crossref Full Text | Google Scholar

37. Bjerregaard P and Larsen CV. Three lifestyle-related issues of major significance for public health among the Inuit in contemporary Greenland: a review of adverse childhood conditions, obesity, and smoking in a period of social transition. Public Health Rev. (2018) 39:5. doi: 10.1186/s40985-018-0085-8

PubMed Abstract | Crossref Full Text | Google Scholar

38. Gallardo-Rincón H, Cantoral A, Arrieta A, Espinal C, Magnus MH, Palacios C, et al. Review: Type 2 diabetes in Latin America and the Caribbean: Regional and country comparison on prevalence, trends, costs and expanded prevention. Prim Care Diabetes. (2021) 15:352–9. doi: 10.1016/j.pcd.2020.10.001

PubMed Abstract | Crossref Full Text | Google Scholar

39. Espina C, Feliu A, Maza M, Almonte M, Ferreccio C, Finck C, et al. Latin America and the Caribbean Code Against Cancer 1st Edition: 17 cancer prevention recommendations to the public and to policy-makers (World Code Against Cancer Framework). Cancer Epidemiol. (2023) 86 Suppl 1:102402. doi: 10.1016/j.canep.2023.102402

PubMed Abstract | Crossref Full Text | Google Scholar

40. Vilares AT, Ciabattoni R, Cunha TM, and Félix A. Cervical cancer in Cape Verde: reappraisal upon referral to a tertiary cancer centre. Ecancermedicalscience. (2022) 16:1471. doi: 10.3332/ecancer.2022.1471

PubMed Abstract | Crossref Full Text | Google Scholar

41. Okyere J, Ayebeng C, and Dickson KS. Factors associated with age at first screening for cervical cancer among adult Cape Verdean women: a cross-sectional study. BMC Public Health. (2024) 24:2444. doi: 10.1186/s12889-024-19917-5

PubMed Abstract | Crossref Full Text | Google Scholar

42. Klumb CE, Hassan R, De ODE, De RLM, MK Carriço, De ADJ, et al. Geographic variation in Epstein-Barr virus-associated Burkitt’s lymphoma in children from Brazil. Int J Cancer. (2004) 108:66–70. doi: 10.1002/ijc.11443

PubMed Abstract | Crossref Full Text | Google Scholar

43. Miller JA, Liu Z, Pinsky BA, Le Q-T, Li T, Yu KJ, et al. Optimization and local cost-effectiveness of nasopharyngeal carcinoma screening strategies in southern China: secondary analysis of the guangdong randomized trial. Cancer Epidemiol Biomarkers Prev. (2024) 33:884–95. doi: 10.1158/1055-9965.EPI-23-1486

PubMed Abstract | Crossref Full Text | Google Scholar

44. Jia W-H, Luo X-Y, Feng B-J, Ruan H-L, Bei J-X, Liu W-S, et al. Traditional Cantonese diet and nasopharyngeal carcinoma risk: a large-scale case-control study in Guangdong, China. BMC Cancer. (2010) 10:446. doi: 10.1186/1471-2407-10-446

PubMed Abstract | Crossref Full Text | Google Scholar

45. Luo J, Chia KS, Chia SE, Reilly M, Tan CS, and Ye W. Secular trends of nasopharyngeal carcinoma incidence in Singapore, Hong Kong and Los Angeles Chinese populations, 1973-1997. Eur J Epidemiol. (2007) 22:513–21. doi: 10.1007/s10654-007-9148-8

PubMed Abstract | Crossref Full Text | Google Scholar

46. Chen L, Zhang Y, Lai S-Z, Li W-F, Hu W-H, Sun R, et al. 10-year results of therapeutic ratio by intensity-modulated radiotherapy versus two-dimensional radiotherapy in patients with nasopharyngeal carcinoma. Oncologist. (2019) 24:e38–45. doi: 10.1634/theoncologist.2017-0577

PubMed Abstract | Crossref Full Text | Google Scholar

47. Palucka AK and Coussens LM. The basis of oncoimmunology. Cell. (2016) 164:1233–47. doi: 10.1016/j.cell.2016.01.049

PubMed Abstract | Crossref Full Text | Google Scholar

48. GBD 2021 Japan Collaborators. Three decades of population health changes in Japan, 1990-2021: a subnational analysis for the Global Burden of Disease Study 2021. Lancet Public Health. (2025) 10:e321–32. doi: 10.1016/S2468-2667(25)00044-1

PubMed Abstract | Crossref Full Text | Google Scholar

49. Tan CC, Lam CS, Matchar DB, Zee YK, and Wong JE. Singapore’s health-care system: key features, challenges, and shifts. Lancet. (2021) 398:1091–104. doi: 10.1016/S0140-6736(21)00252-X

PubMed Abstract | Crossref Full Text | Google Scholar

50. Zeng D-Y, Li J-M, Lin S, Dong X, You J, Xing Q-Q, et al. Global burden of acute viral hepatitis and its association with socioeconomic development status, 1990-2019. J Hepatol. (2021) 75:547–56. doi: 10.1016/j.jhep.2021.04.035

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: Bayesian age-period-cohort model, global burden of disease, incidence, mortality, nasopharyngeal carcinoma

Citation: Zhou E, Xu F, Zhou T, Zhang J, Song F, Li J, Yi H, Wang Q and Huang W (2026) Shifting burden of nasopharyngeal carcinoma: global patterns and forecasts to 2050 from the GBD 2021. Front. Oncol. 15:1687320. doi: 10.3389/fonc.2025.1687320

Received: 18 August 2025; Accepted: 08 December 2025; Revised: 01 December 2025;
Published: 06 January 2026.

Edited by:

Faris Lami, University of Baghdad, Iraq

Reviewed by:

Seyed Amir Ahmad Safavi-Naini, Shahid Beheshti University of Medical Sciences, Iran
Chen Dong, Chung Shan Medical University, Taiwan

Copyright © 2026 Zhou, Xu, Zhou, Zhang, Song, Li, Yi, Wang and Huang. 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: Jianxiang Li, bGlqaWFueGlhbmcyMDAyQDE2My5jb20=; Qingliang Wang, d2FuZ3FsNUBtYWlsLnN5c3UuZWR1LmNu; Weijun Huang, aGVsbG9odWFuZ3dqQDEyNi5jb20=

These authors have contributed equally to this work

Disclaimer: 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.