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

ORIGINAL RESEARCH article

Front. Endocrinol., 06 February 2026

Sec. Cancer Endocrinology

Volume 17 - 2026 | https://doi.org/10.3389/fendo.2026.1674679

This article is part of the Research TopicAI-Powered Insights: Predicting Treatment Response and Prognosis in Breast CancerView all 18 articles

Global trends and disease burden of elderly male breast cancer, 1990-2021: a population-based study

Shaochun Liu&#x;Shaochun Liu1†Jian Xu&#x;Jian Xu2†Han Zhang&#x;Han Zhang1†Yuhan TangYuhan Tang1Xiaoxi HanXiaoxi Han1Wenjing XiongWenjing Xiong1Linlin FanLinlin Fan1Chang SuChang Su1Zhiqi SuiZhiqi Sui1Wenhui Zhao*Wenhui Zhao1*
  • 1Department of Medical Oncology, Harbin Medical University Cancer Hospital, Harbin Medical University, Harbin, Heilongjiang, China
  • 2Department of Oncology, the Second Affiliated Hospital of Anhui Medical University, Hefei, Anhui, China

Objective: To assess global, regional, and national trends in the burden of elderly male breast cancer (EMBC) from 1990 to 2021 and to evaluate projected patterns to 2040.

Methods: Estimates were obtained from the Global Burden of Disease 2021 study. Age-standardised incidence, mortality, and disability-adjusted life year rates (ASIR, ASMR, ASDR) were analysed across all countries and Sociodemographic Index (SDI) strata. Long-term changes were quantified using the average annual percent change derived from log-linear models. Joinpoint regression identified temporal inflection points. Age–period–cohort (APC) models characterise independent temporal effects. Mortality changes were decomposed into components attributable to population growth, ageing, and epidemiological change. Inequality was assessed using slope and concentration indexes. Attributable mortality and DALYs were evaluated for alcohol use, dietary risks, and tobacco. Future rates to 2040 were estimated using a Bayesian age–period–cohort (BAPC) model.

Results: Globally, EMBC incidence, mortality, and DALYs increased from 1990 to 2021, with average annual percent changes(AAPC) of 1.8 (95% CI, confidence interval, 1.63 to 1.98), 0.58 (95% CI 0.38 to 0.77), and 0.68 (95% CI 0.43 to 0.92). East Asia recorded the steepest rise in incidence, increasing from 1.65 (95% UI, uncertainty interval, 1.16 to 2.42) to 6.65 per 100000 population (95% UI 2.77 to 9.78). The middle SDI quintile showed the largest increases in all three metrics, rising from 1.62 (95% UI 1.18 to 2.09) to 4.92 per 100000 population (95% UI 2.26 to 6.81). APC analysis indicated pronounced period and cohort effects in middle and low-middle SDI settings. Decomposition identified population growth as the dominant driver of rising burden. Alcohol use and dietary risks accounted for most increases in mortality and DALYs, while tobacco contributed minimally. Cross-country inequality was modest for incidence but more marked for mortality and DALYs. Projections suggest that age-standardised rates will decline gradually by 2040, although absolute case numbers may continue to rise in rapidly ageing regions.

Conclusion: The global burden of EMBC continues to increase, with substantial regional and socioeconomic disparities. Although age-standardised rates are projected to decline, population ageing is expected to sustain or expand absolute numbers of cases and deaths.

Introduction

Breast cancer arises from the malignant transformation of mammary epithelial cells and constitutes a leading contributor to global cancer morbidity and mortality. Classified primarily into invasive and non-invasive subtypes, it is closely linked to genetic susceptibility, hormonal exposure, and environmental determinants (1). As the most frequently diagnosed cancer worldwide, breast cancer also ranks as the leading cause of disability-adjusted life years (DALYs) across all malignancies (2).

Although male breast cancer (MBC) accounts for less than 1% of all newly diagnosed breast cancers (3), emerging evidence underscores its distinct biological behaviour and clinical manifestations compared with female breast cancer (FBC). Men typically present with larger tumour size, higher nodal involvement rates, and increased likelihood of distant metastasis at diagnosis (4), mainly attributable to delayed recognition of breast-related symptoms and the absence of gender-specific routine screening initiatives. Long-term Surveillance, Epidemiology, and End Results (SEER) data demonstrate that the age-adjusted incidence of MBC rose from 0.85 per 100,000 in 1975 to 1.43 in 2011, with a mean age at diagnosis of approximately 67 years (5). The global burden of MBC has steadily increased over recent decades (6).

This pattern reflects the rapid ageing of the global population. The risk of cancer increases markedly with advancing age (7). According to the WHO and the Lancet Commission on Ageing, individuals aged 60 years or older are internationally recognized as older adults (8). The World Health Organization projects that more than 20% of the global population will be aged 60 years or older by the mid-21st century, indicating that the cancer burden in older adults will continue to grow and contribute substantially to overall morbidity (9). The rising burden of cancer in this population poses a major public health challenge. It has direct implications for achieving Sustainable Development Goal 3.4, which targets reductions in premature mortality from non-communicable diseases through prevention and treatment (10).

In the present study, EMBC is explicitly defined as male breast cancer occurring in adults aged ≥60 years. Older men exhibit unique vulnerabilities, including a higher prevalence of BRCA2 mutations, greater multimorbidity burden, diminished physiological reserve, under-recognition of breast symptoms, and frequent misclassification as benign conditions such as gynecomastia (11, 12). Well-documented diagnostic delays in this population are associated with more advanced disease at presentation and poorer outcomes, highlighting the need for elderly-specific surveillance and intervention strategies beyond those developed for younger or all-age male populations.

Most epidemiological research on breast cancer has focused on women, and evidence specific to older men remains limited. Existing GBD-based assessments have described global male breast cancer patterns, yet they predominantly report all-age MBC and therefore obscure age-related risk structures, diagnostic pathways, and demographic pressures that shape the burden in older adults (13). The Global Burden of Disease (GBD) 2021 study offers a harmonised, methodologically consistent framework for evaluating these patterns across countries and over time (14). Within this context, a dedicated assessment of elderly male breast cancer (EMBC) is needed to determine whether temporal trajectories, regional disparities, and risk attributions diverge from those observed in younger men or in the male population as a whole, particularly as health systems confront the accelerating implications of population ageing.

To further elucidate the mechanisms underlying temporal change, the age-period-cohort (APC) model distinguishes biological ageing effects, period influences related to diagnostic or therapeutic evolution, and cohort-specific exposures shaped by generational context. This analytical structure provides a multidimensional interpretation of long-term incidence and mortality trends (15). Using GBD 2021 estimates, the present study applies APC modelling to characterise global, regional, and national trends in EMBC from 1990 to 2021, offering a comprehensive evaluation of its evolving burden in an ageing world.

Methods

Data source

Data for EMBC, defined as breast cancer occurring in men aged ≥60 years, were obtained from the GBD 2021 study. This study represents a secondary analysis of publicly available, de-identified GBD data. GBD 2021 provides harmonised estimates of incidence, mortality, and DALYs for 371 diseases and 88 risk factors across 204 countries and territories from 1990 to 2021. Estimates are generated using a hierarchical Bayesian meta-regression framework implemented in DisMod-MR 2.1, which incorporates epidemiological data, predictive covariates, and location-specific priors to produce internally consistent measures at global, regional, and national levels (16). The Sociodemographic Index (SDI), a composite indicator of income, education, and fertility, was used to stratify locations by development level (14). Detailed GBD query specifications, including cause codes, age groups, locations, extraction interface, and run date, are provided in the Supplementary Methods.

Analysis of trends globally, five SDI quintiles, and 204 countries and territories

Age-standardized incidence (ASIR), mortality (ASMR), and disability-adjusted life year rates (ASDR) were evaluated from 1990 to 2021. Long-term temporal trends were quantified using the Average Annual Percent Change (AAPC), estimated through log-linear regression of annual standardized rates (17). Joinpoint regression was used to identify inflection points and to calculate annual percent change (APC) across detected segments (18). Model specifications, including the number of joinpoints tested, permutation procedures, and approaches to uncertainty estimation, are described in the Supplementary Methods.

Age–period–cohort modelling

To assess the independent contributions of age, period, and cohort to observed temporal patterns, age–period–cohort models were fitted to GBD estimates. The selection of interval widths, reference categories, and the identifiability solution followed established methodological principles. The intrinsic estimator (IE) approach was applied to address the linear dependency among the three temporal dimensions (19). Details on model implementation, software versions, procedures for propagating GBD uncertainty intervals, and sensitivity analyses are provided in the Supplementary Methods.

Decomposition analysis

Changes in EMBC mortality between 1990 and 2021 were decomposed into components attributable to population growth, population ageing, and changes in age-specific mortality rates. Interaction effects were also quantified to capture combined demographic and epidemiological influences (20). Full mathematical expressions and a numerical example illustrating component contributions are presented in the Supplementary Methods.

Cross-country inequality analysis

Cross-national disparities in EMBC burden were evaluated using the slope index of inequality (SII) to quantify absolute inequality and the concentration index (CI) to quantify relative inequality. SII was estimated by regressing standardized outcomes on the SDI-based relative rank of each country, and CI was computed from the cumulative distribution of DALYs across ranked populations (21). Ranking procedures, weighting approaches, tie treatment, and statistical methods for uncertainty estimation are detailed in the Supplementary Methods.

Risk factor analysis

Temporal changes in ASMR and ASDR attributable to dietary risks, alcohol use, and tobacco were quantified using the risk–outcome pairs defined in the GBD comparative risk assessment framework. This framework integrates exposure distributions, relative risks, and theoretical minimum risk exposure levels to estimate the burden attributable to modifiable risk factors across locations and years.

Projection analysis

Projections of ASIR, ASMR, and ASDR for EMBC to 2040 were generated at the global level using a Bayesian age–period–cohort (BAPC) model (22). This approach applies Bayesian priors to stabilise age, period, and cohort effects and to propagate uncertainty, projecting future rates under the assumption that recently observed trends persist (23).

Results

Global trends

Table 1, Supplementary Tables S1S3 present the ASIR, ASMR, ASDR, and AAPC for EMBC across the global, five SDI quintiles, and 21 regions and 204 countries and territories in 1990 and 2021. The global trend analysis reveals that the ASIR, ASMR, and ASDR of EMBC all showed an increasing trend from 1990 to 2021, with AAPCs of 1.8 [95% CI, confidence interval, 1.63 to 1.98], 0.58 [95% CI 0.38 to 0.77], and 0.68 [95% CI 0.43 to 0.92], respectively. As shown in Figure 1, a global Joinpoint regression analysis identified a notable turning point in ASIR in 2013, with the APC dropping from 3.16 (2004–2013) to 0.24 (2013–2021). Similarly, both ASMR and ASDR experienced significant declines in 2013 and 2012, respectively, with APCs decreasing from 0.6 (2000–2013) to -0.4 (2013–2021) for ASMR and from 1.22 (2004–2012) to –0.29 (2012–2021) for ASDR.

Table 1
www.frontiersin.org

Table 1. Incident cases and ASIR of EMBC in 1990 and 2021, and its temporal trends from 1990 to 2021.

Figure 1
Three line graphs depict trends in male breast cancer data from 1990 to 2021. Graph A shows incidence rates rising sharply through 2013, then stabilizing. Graph B illustrates increasing death rates until 2013, followed by a slight decline. Graph C presents Disability-Adjusted Life Years (DALYs) with notable increases until 2012, then a decline. Each graph highlights joinpoints reflecting significant trend changes, with specific Average Percent Change (APC) values for each period. Color-coded lines represent different time segments, accompanied by legends explaining APC values.

Figure 1. Joinpoint regression analysis of EMBC for ASIR (A), ASMR (B), and ASDR (C), 1990–2021. Age definition: EMBC refers to men aged ≥60 years. Age standardization: Standardized to the GBD 2017 world standard population. Unit: per 100,000 population. EMBC, elderly male breast cancer; ASIR, age-standardized incidence rate; ASMR, age-standardized mortality rate; ASDR, age-standardized disability-adjusted life years; APC, annual percent change; AAPC, average annual percent change.

Trends across the five SDI quintiles

ASIR increased across all SDI quintiles. The most notable increase occurred in the middle SDI quintile, where it rose from 1.62 per 100,000 population (95% UI, uncertainty interval, 1.18 to 2.09) in 1990 to 4.92 per 100,000 population (95% UI 2.26–6.81) in 2021, with an AAPC of 3.71 (95% CI 3.36 to 4.07). Similarly, the most significant increases in ASMR and ASDR were observed in the middle SDI quintile, with AAPCs of 1.56 (95% CI 1.35 to 1.78) and 1.75 (95% CI 1.56 to 1.93), respectively. In contrast, the smallest increase in ASIR was seen in the high SDI quintile (AAPC 0.35 [95% CI –0.35 to 1.05]), and the high SDI quintile showed the most pronounced declines in ASDR and ASMR among the five SDI quintiles, with AAPCs of –0.47 (95% CI –0.99 to 0.06) and –0.37 (95% CI –0.93 to 0.2).

Trends across 21 regions and 204 countries and territories

By region, East Asia recorded the highest ASIR in 2021, whereas Eastern Sub-Saharan Africa had the highest ASMR and ASDR. Between 1990 and 2021, East Asia also experienced the most pronounced increase in ASIR, rising from 1.65 per 100,000 population (95% UI 1.16 to 2.42) in 1990 to 6.65 per 100,000 population (95% UI 2.77 to 9.78) in 2021, with an AAPC of 4.6 (95% CI 4.24 to 4.96). The ASMR and ASDR in Eastern Sub-Saharan Africa reached 14.41 per 100,000 population (95% UI 9.41 to 25.39) and 269.29 per 100,000 population (95% UI 174.93 to 476.59), respectively. In contrast, the ASIR in Western Sub-Saharan Africa showed no significant temporal change, with an AAPC of 0.05 (95% CI to0.02 to 0.12).

At the country level, Uganda had the highest ASIR and ASMR for elderly male breast cancer in both 1990 and 2021. The ASIR was 17.85 per 100,000 population in 1990 (95% UI 9.43 to 30.75) and increased to 24.98 per 100,000 population in 2021 (95% UI 12.48 to 45.27). The ASMR was 16.87 per 100,000 population in 1990 (95% UI 8.83 to 29.1) and reached 21.17 per 100,000 population in 2021 (95% UI 10.65 to 37.55). The country with the highest ASDR in 2021 was Malawi, with a rate of 308.26 per 100,000 population (95% UI 130.87 to 603.06).

Figure 2 presents the differences in ASIR and the corresponding AAPCs across countries and regions in 1990 and 2021.

Figure 2
Nine line graphs divided into three panels (A, B, C) each containing three graphs. The graphs show trends by age, periods, and birth cohorts, respectively. Data is grouped by SDI levels: Global, Low, Low-middle, Middle, High-middle, and High SDI, represented by different colored lines. Each graph illustrates changes in annual rates, rate ratios, or risk ratios across age groups (60-64 to 95-plus), five-year periods (1992-2021), and birth cohorts (1892-1981). Shaded areas indicate confidence intervals.

Figure 2. Global map of ASIR of EMBC in 1990 (A) and 2021 (B), and corresponding AAPC from 1990 to 2021 (C). Age definition: EMBC refers to men aged ≥60 years. Age standardization: Age-standardized to the GBD 2017 world standard population. Unit: per 100,000 population. EMBC, elderly male breast cancer; ASIR, age-standardized incidence rate; AAPC, average annual percentage change.

Age-period-cohort analysis on ASIR, ASMR, and ASDR

Figure 3 summarises the age, period and cohort patterns of incidence, mortality and DALYs for EMBC at the global level and across SDI categories. The incidence results show small positive local drifts in the younger elderly groups and values close to zero or slightly negative in advanced ages. The period effects show consistent increases, and the cohort effects rise sharply in later birth cohorts, particularly in middle- and low-middle-SDI regions. In contrast, the mortality and DALYs results show uniformly negative local drifts across ages, modest period increases in most low-SDI settings and cohort effects that decline towards the mid-century birth cohorts before levelling off or showing only mild recovery in later cohorts. High-SDI regions display the most stable incidence drifts, the smallest period increases and cohort patterns that remain close to unity for incidence and below unity for both mortality and DALYs.

Figure 3
Three line graphs labeled A, B, and C show age-standardized rates per one hundred thousand populations from 1990 to 2040. Each graph includes a solid line for observed data until 2020 and a dotted line for predicted data post-2020. A shaded area indicates data confidence levels. The trend generally peaks around 2020 and declines afterward.

Figure 3. Age–period–cohort effects on incidence (A), mortality (B), and DALYs (C) of EMBC globally and by SDI level. Age definition: EMBC refers to men aged ≥60 years. Unit: per 100,000 population. EMBC, elderly male breast cancer; SDI, Sociodemographic Index; DALYs, disability-adjusted life years.

Decomposition analysis

Figure 4 presents the decomposition analysis of ASIR, ASMR, and ASDR for breast cancer in older males across different SDI quintiles, highlighting the roles of aging, population growth, and epidemiological changes. Regarding ASIR, epidemiological changes exerted a greater influence than ASMR and ASDR, particularly in the middle and low-middle SDI quintiles. Population growth remained the predominant driver of increases in these regions, while the effect of aging remained relatively limited across all SDI levels. For ASMR, population growth emerged as the primary contributor, with the most pronounced changes observed in the middle and low SDI quintiles. The influence of aging was most significant in the high SDI quintile, though its effect was less notable elsewhere. The impact of epidemiological changes varied, with marked contributions seen in the middle SDI quintile.

Figure 4
Bar chart showing the breakdown of Disability-Adjusted Life Years (DALYs), deaths, and incidence by location and influencing factors. Locations include global, high, high-middle, middle, low-middle, and low SDI. Each bar is subdivided into percent change of aging, population, and epidemiological change (purple, teal, yellow). Black dots mark specific values within each category.

Figure 4. Decomposition of changes in incidence, mortality, and DALYs of EMBC from 1990 to 2021 by SDI quintile. DALYs, disability-adjusted life years; SDI, Sociodemographic Index.

Regarding ASDR, population growth accounted for the largest share of the increase across all SDI quintiles, particularly in middle- and low-SDI areas. The influence of aging was more prominent in the high and high-middle SDI quintiles. At the same time, the effects of epidemiological changes were generally modest, with some increase seen in the middle and low SDI quintile. In summary, the decomposition analysis underscores that population growth is the leading factor behind the rising breast cancer burden in older males, with aging playing a more prominent role in higher SDI quintiles and epidemiological changes contributing more variably, especially in middle and low SDI quintiles. Supplementary Tables S4S6 detail the changes in incidence, deaths, and DALYs attributed to population-level determinants and causes from 1990 to 2021.

Cross-country inequality analysis

There is a degree of absolute and relative inequality in SDI observed in breast cancer burden, although it is not particularly significant (Figure 5). The differences in the concentration index, which measures relative gradient inequality, between 1990 and 2021 show that while the concentration index for ASIR was 0.12 (95% CI: 0.05–0.19) in 1990 and decreased to 0.06 (95% CI: 0.01–0.10) in 2021, indicating a narrowing imbalance in the distribution of incidence rates across countries with different SDI, the inequality remains. In contrast, the concentration index for ASMR was -0.12 (95% CI: -0.21 to -0.04) in 1990 and -0.17 (95% CI: -0.23 to -0.11) in 2021, while the concentration index for ASDR was -0.12 (95% CI: -0.21 to -0.03) in 1990 and -0.16 (95% CI: -0.22 to -0.10) in 2021. These findings suggest that the distribution of mortality and DALYs burden across countries with different SDIs shows an increasing imbalance.

Figure 5
Scatter and line charts depicting health data across six panels, labeled A to F. Panels A, C, and E display crude incidence, death, and DALY rates per 100,000 against relative SDI rank for 1990 and 2021, with bubble sizes representing population. Panels B, D, and F show cumulative fractions of incidence, deaths, and DALY against population ranked by SDI, highlighting concentration indices for 1990 and 2021. Blue represents 1990 data, while red indicates 2021. The charts collectively illustrate trends and disparities in health data over time and by country, with notable emphasis on data from China.

Figure 5. SDI-related inequality in ASIR, ASMR, and ASDR of EMBC in 1990 and 2021: regression-based slope index of inequality (A, C, E) and concentration curves (B, D, F). Age definition: EMBC refers to men aged ≥60 years. Age standardization: Age-standardized to the GBD 2017 world standard population. Unit: per 100,000 population. EMBC, elderly male breast cancer; ASIR, age-standardized incidence rate; ASMR, age-standardized mortality rate; ASDR, age-standardized disability-adjusted life years; SDI, Sociodemographic Index.

Risk factor analysis

Figure 6 depicts the temporal trends in ASMR and ASDR for EMBC from 1990 to 2021, attributable to three principal risk factors at both global and regional levels. Globally, alcohol consumption emerged as the foremost contributor to increased ASMR and ASDR, with dietary risks following as the second most significant factor. From 1990 to 2015, the burden linked to alcohol and dietary risks steadily escalated before declining between 2015 and 2021.

Figure 6
Graphs showing trends in ASMR and ASDR from 1990 to 2020 across global and different SDI regions: Low, Low-middle, Middle, High-middle, and High. Each graph details the impact of alcohol use, dietary risks, and tobacco, represented by blue, green, and red areas, respectively.

Figure 6. Temporal trends in ASMR (A) and ASDR (B) of EMBC attributable to alcohol use, dietary risks, and tobacco across SDI quintiles and globally, 1990–2021. Age definition: EMBC refers to men aged ≥60 years. Age standardization: Age-standardized to the GBD 2017 world standard population. Unit: per 100,000 population. EMBC, elderly male breast cancer; ASMR, age-standardized mortality rate; ASDR, age-standardized disability-adjusted life years; SDI, Sociodemographic Index.

At the SDI quintile level, the low SDI quintile initially experienced a gradual decline in the impact of alcohol and dietary factors on ASMR and ASDR from 1990, only to see a modest resurgence around 2010. In contrast, the low‐middle quintile exhibited a continuous rise in alcohol‐driven ASMR and ASDR throughout the period, while the middle quintile showed an upward trend from 1990 that reversed in the early 2010s. The high‐middle quintile maintained relative stability from 1990 to 2021, and in the high quintile, ASMR and ASDR increased from 1990 before gradually declining in the 21st century. Notably, tobacco use had the least impact across all quintiles, although a slight upward trend was observed in the low‐middle and middle SDI groups.

BAPC prediction

Figure 7 shows that the projected global burden of EMBC declines steadily after 2021. By 2040, the predicted ASIR, ASMR, and ASDR for EMBC are 4.06 (95% UI 1.45–7.08), 1.63 (95% UI 0.77–2.49), and 32.85 (95% UI 16.15–49.55), respectively.

Figure 7
Three world maps labeled A, B, and C, showing ASIR and AAPC values by region. Colors represent different ranges: yellow, blue, red, gray, and others. Insets feature Caribbean, Central America, Persian Gulf, Balkan Peninsula, Southeast Asia, West Africa, Eastern Mediterranean, and Northern Europe.

Figure 7. Observed and projected global ASIR (A), ASMR (B), and ASDR (C) of EMBC, 1990–2040. Age definition: EMBC refers to men aged ≥60 years. Age standardization: Age-standardized to the GBD 2017 world standard population. Unit: per 100,000 population. EMBC, elderly male breast cancer; ASIR, age-standardized incidence rate; ASMR, age-standardized mortality rate; ASDR, age-standardized disability-adjusted life years.

Discussion

This study provides the first comprehensive global assessment of EMBC from 1990 to 2021. The results show that EMBC follows a temporal trajectory that is distinct from both all-age MBC and female breast cancer. Incidence, mortality and DALYs rose steadily for over two decades and then slowed after 2012–2013, consistent with the turning points detected in Joinpoint analyses. This synchronous rise contrasts with the global pattern in female breast cancer, where incidence has continued to rise while mortality has declined steadily over the same period (2426). They also differ from all-age MBC, which generally shows an earlier separation between rising incidence and declining mortality (13). This distinction confirms the added value of restricting analysis to older men and demonstrates that EMBC cannot be inferred from existing all-age estimates.

The unique epidemiology of EMBC likely reflects the intersection of ageing biology, male-specific risk profiles and health system limitations. Older men have a higher prevalence of BRCA2 mutations, multimorbidity and diminished physiologic reserve, and they commonly experience diagnostic delays because breast symptoms are often misattributed to benign causes such as gynecomastia (11, 12, 27). These factors contribute to more advanced stage at diagnosis and may explain why mortality trends in EMBC mirror incidence more closely than in female breast cancer, where long-standing investment in screening and treatment has produced substantial survival gains (25, 26).

Marked geographic heterogeneity was observed. East Asia shows the highest incidence in 2021 and the steepest long-term increase, which corresponds with rapid population ageing, lifestyle transitions, and improved registry completeness in the region (28, 29). At the same time, several African countries exhibit the highest mortality and DALYs, which likely reflects persistent barriers to diagnosis, limited pathology infrastructure, delayed presentation and constrained access to therapy. These findings align with broader evidence that many cancers in Africa continue to show limited improvements in survival despite global progress (30).

High SDI regions demonstrated earlier declines in incidence, mortality, and DALYs than other SDI quintiles. This may reflect stabilised exposure to major risk factors, broader access to diagnostic imaging and therapeutic advances. However, older men with breast cancer benefit less from endocrine therapy and targeted agents such as CDK4/6 inhibitors compared with younger patients because multimorbidity, frailty and drug tolerance frequently limit treatment intensification (6, 12). Competing risks from cardiovascular disease, diabetes, and chronic pulmonary conditions also attenuate the mortality reductions that might otherwise accompany therapeutic advances (31).

Age–period–cohort analysis provides additional insight into these patterns. Incidence increases most consistently across younger elderly groups and in middle and low-middle SDI settings, suggesting that later birth cohorts may have accumulated higher lifetime risk. The rise in cohort risk after 1962 mirrors patterns observed for hormone receptor positive female breast cancer in Scotland and elsewhere (32). Period effects demonstrate broad global increases, which likely reflect enhanced diagnostic access and improved cancer registration, particularly in evolving health systems across Asia and Latin America (29). Mortality and DALYs show negative local drifts across many age groups, although gains are not as marked as those observed for women.

The decomposition analysis highlights population growth as the dominant driver of the global rise in EMBC, particularly in middle- and low-SDI regions. Ageing contributes more substantially in high SDI regions. Epidemiological change plays a varying role but is most pronounced in settings undergoing rapid urbanisation and metabolic transition, which mirrors patterns observed in colorectal cancer and pancreatic cancer among older adults (31, 33).

Although cross-country inequality remained modest in absolute terms, its direction differed across indicators. Consistent with our inequality analysis, incidence inequality narrowed slightly from 1990 to 2021. In contrast, mortality and DALYs became increasingly concentrated in lower SDI countries, reflecting a widening fatal burden in settings with limited diagnostic and treatment capacity. These patterns align with evidence that cancer survival improvements have disproportionately benefited higher-SDI regions, while lower-resource countries continue to experience persistently higher mortality despite similar incidence levels (30, 34). Without targeted strategies, these disparities are likely to grow as populations age.

The interpretation of risk factors requires caution because most etiologic evidence for breast cancer originates from studies in women, whereas research in older men remains scarce. In our analysis, alcohol and dietary risks were the predominant contributors to EMBC mortality and DALYs, increasing steadily from 1990 to 2015 before declining modestly thereafter. In contrast, tobacco consistently accounted for the smallest proportion and showed only slight increases in low-middle and middle SDI regions (13). These temporal patterns align with findings from multinational pooling studies of male breast cancer, which report modest but significant associations between alcohol intake and cancer risk, with odds ratios typically ranging from 1.2 to 1.8 for moderate-to-high consumption (35, 36). In contrast, female breast cancer meta-analyses demonstrate a more gradual dose–response relationship, with risk increasing by approximately 7 to 10 percent per 10 g/day of alcohol and reaching relative risks of 1.3 to 1.5 only at higher intake levels (37). Beyond differences in effect magnitude, underlying mechanisms may diverge between sexes: in men, alcohol-related carcinogenesis appears more strongly linked to hepatic dysfunction and impaired oestrogen clearance, whereas studies in women emphasise direct hormonal pathways (35). The regional heterogeneity in our results, including sharper alcohol- and diet-related burdens in middle and low-middle SDI settings, likely reflects variation in exposure patterns, metabolic profiles and health system resources. These discrepancies underscore the uncertainty in extrapolating evidence from female-dominated studies to EMBC and highlight the need for male-specific epidemiologic research. The prominent role of alcohol and metabolic risks observed in our findings supports integrated non-communicable disease prevention strategies tailored to ageing male populations.

Based on projections from the BAPC model, the global burden of EMBC is expected to decline, a pattern that mirrors forecasts for all-age male breast cancer. In clear contrast, the projected incidence and mortality of female breast cancer continue to rise through 2050 (38). This divergence likely reflects the influence of multiple risk factors that drive female breast cancer globally, with shifting exposure patterns compounded by demographic transitions that enlarge the pool of individuals at risk (39). Concurrently, delayed childbearing, rising obesity and metabolic risk, increasing breast density, and accumulated hormone-related exposures in several countries further reinforce the sustained upward trajectory of female breast cancer (24). It is important to note that although age-standardised rates of EMBC are declining, the absolute number of cases and deaths may continue to increase in regions experiencing rapid population ageing, particularly in East Asia, Eastern Europe, and several middle-to-high SDI settings where the population of older men is expanding substantially (40).

Diagnostic improvements may also explain portions of the observed temporal changes. Increased utilisation of breast imaging modalities among high-risk men, including mammography, ultrasound and MRI, has enhanced case detection and contributed to greater awareness of male breast disease in clinical practice (41). Therapeutic advances in endocrine therapy, chemotherapy, and targeted agents have improved outcomes for some older men, although treatment limitations related to frailty and polypharmacy persist (5). Older patients face substantial competing mortality from cardiovascular disease, diabetes, chronic respiratory diseases, and frailty-related conditions that may attenuate survival gains even when cancer control improves (42).

This study has several limitations. First, GBD estimates depend on statistical modelling in data-sparse settings where registration systems are incomplete, which may influence regional comparisons (39). Second, the ecological association between SDI and EMBC outcomes carries inherent limitations and may not capture individual-level risk. Third, GBD does not include cancer stage, molecular subtype or treatment information, which limits interpretation of age-specific pathways. Fourth, independent validation using data sources such as SEER is limited for elderly male populations. Fifth, the large number of comparisons across SDI levels, regions and countries increases the likelihood of false-positive findings. Although consistent direction of patterns and sensitivity analyses mitigate this concern, caution is warranted when interpreting marginal differences. Finally, uncertainty intervals around incidence, mortality and DALYs should be considered when comparing country ranks and regional contrasts.

Conclusion

In conclusion, EMBC displays a distinct epidemiologic trajectory, marked by long-term increases in ASIR, ASMR, and ASDR and emerging evidence of recent stabilisation in selected high-SDI settings. BAPC projections indicate a gradual decline in global ASIR, ASMR, and ASDR through 2040. Persistent excess mortality and DALYs in low-SDI countries reflect constrained diagnostic capacity, delayed presentation, and limited access to systemic therapy. Alcohol and dietary risks remain the dominant contributors to the attributable burden, while cohort patterns point to heightened susceptibility among generations born after the early 1960s. These findings underscore an urgent need for region-tailored strategies that include awareness initiatives for older men, integration of context-appropriate early detection approaches, strengthened pathology infrastructure, and equitable access to treatment. Such investments will be essential to mitigate preventable mortality and disability from EMBC as global populations continue to age.

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

Author contributions

SL: Writing – original draft. JX: Writing – review & editing. HZ: Writing – review & editing. YT: Writing – review & editing. XH: Writing – review & editing. WX: Writing – review & editing. LF: Writing – review & editing. CS: Writing – review & editing. ZS: Writing – review & editing. WZ: Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the Beijing Medical Award Foundation (YXJL−2021−0424−0382), the Harbin Medical University Cancer Hospital Haiyan Scientific Research Fund (JJZD2022−14), and the Natural Science Foundation of Heilongjiang Province (PL2024H183).

Acknowledgments

Thanks to Xiao Ming for his work in the GBD database. The authors also thank Home for Researchers (www.home-for-researchers.com) for providing professional English language editing services.

Conflict of interest

The author(s) 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) declared that generative AI was used in the creation of this manuscript. During the preparation of this work, the authors used chatGPT to polish the language. After using this tool, the authors reviewed and edited the content as needed and took full responsibility for the content of the publication.

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

References

1. Barzaman K, Karami J, Zarei Z, Hosseinzadeh A, Kazemi MH, Moradi-Kalbolandi S, et al. Breast cancer: Biology, biomarkers, and treatments. Int Immunopharmacol. (2020) 84:106535. doi: 10.1016/j.intimp.2020.106535, PMID: 32361569

PubMed Abstract | Crossref Full Text | Google Scholar

2. Coles CE, Anderson BO, Cameron D, Cardoso F, Horton R, Knaul FM, et al. The Lancet Breast Cancer Commission: tackling a global health, gender, and equity challenge. Lancet (London England). (2022) 399:1101–3. doi: 10.1016/S0140-6736(22)00184-2, PMID: 35189077

PubMed Abstract | Crossref Full Text | Google Scholar

3. Siegel RL, Miller KD, Fuchs HE, and Jemal A. Cancer Statistics, 2021. CA: Cancer J Clin. (2021) 71:7–33. doi: 10.3322/caac.21654, PMID: 33433946

PubMed Abstract | Crossref Full Text | Google Scholar

4. Gucalp A, Traina TA, Eisner JR, Parker JS, Selitsky SR, Park BH, et al. Male breast cancer: a disease distinct from female breast cancer. Breast Cancer Res Treat. (2019) 173:37–48. doi: 10.1007/s10549-018-4921-9, PMID: 30267249

PubMed Abstract | Crossref Full Text | Google Scholar

5. Giordano SH. Breast cancer in men. New Engl J Med. (2018) 378:2311–20. doi: 10.1056/NEJMra1707939, PMID: 29897847

PubMed Abstract | Crossref Full Text | Google Scholar

6. Fox S, Speirs V, and Shaaban AM. Male breast cancer: an update. Virchows Archiv: an Int J Pathol. (2022) 480:85–93. doi: 10.1007/s00428-021-03190-7, PMID: 34458944

PubMed Abstract | Crossref Full Text | Google Scholar

7. Bisht S, Mao Y, and Easwaran H. Epigenetic dynamics of aging and cancer development: current concepts from studies mapping aging and cancer epigenomes. Curr Opin Oncol. (2024) 36:82–92. doi: 10.1097/CCO.0000000000001020, PMID: 38441107

PubMed Abstract | Crossref Full Text | Google Scholar

8. Prince MJ, Wu F, Guo Y, Gutierrez Robledo LM, O’Donnell M, Sullivan R, et al. The burden of disease in older people and implications for health policy and practice. Lancet (London England). (2015) 385:549–62. doi: 10.1016/S0140-6736(14)61347-7, PMID: 25468153

PubMed Abstract | Crossref Full Text | Google Scholar

9. Fane M and Weeraratna AT. How the ageing microenvironment influences tumour progression. Nat Rev Cancer. (2020) 20:89–106. doi: 10.1038/s41568-019-0222-9, PMID: 31836838

PubMed Abstract | Crossref Full Text | Google Scholar

10. Murthy SS, Trapani D, Cao B, Bray F, Murthy S, Kingham TP, et al. Premature mortality trends in 183 countries by cancer type, sex, WHO region, and World Bank income level in 2000-19: a retrospective, cross-sectional, population-based study. Lancet Oncol. (2024) 25:969–78. doi: 10.1016/S1470-2045(24)00274-2, PMID: 38964357

PubMed Abstract | Crossref Full Text | Google Scholar

11. Khan NAJ and Tirona M. An updated review of epidemiology, risk factors, and management of male breast cancer. Med Oncol (Northwood London England). (2021) 38:39. doi: 10.1007/s12032-021-01486-x, PMID: 33721121

PubMed Abstract | Crossref Full Text | Google Scholar

12. Shinde SD, Satpute DP, Behera SK, and Kumar D. Computational Biology of BRCA2 in Male Breast Cancer, through Prediction of Probable nsSNPs, and Hit Identification. ACS omega. (2022) 7:30447–61. doi: 10.1021/acsomega.2c03851, PMID: 36061650

PubMed Abstract | Crossref Full Text | Google Scholar

13. Li Y, Huang Y, Huang H, Wei T, Zhang A, Xing L, et al. Global, regional, and national burden of male breast cancer in 204 countries and territories: a systematic analysis from the global burden of disease study, 1990-2021. EClinicalMedicine. (2025) 80:103027. doi: 10.1016/j.eclinm.2024.103027, PMID: 39831130

PubMed Abstract | Crossref Full Text | Google Scholar

14. Ong KL, Stafford LK, McLaughlin SA, Boyko EJ, Vollset SE, Smith AE, et al. Global, regional, and national burden of diabetes from 1990 to 2021, with projections of prevalence to 2050: a systematic analysis for the Global Burden of Disease Study 2021. Lancet (London England). (2023) 402:203–34. doi: 10.1016/S0140-6736(23)01301-6, PMID: 37356446

PubMed Abstract | Crossref Full Text | Google Scholar

15. Fosse E and Winship C. Analysing age-period-cohort data: A review and critique. Annu Rev Sociol. (2019) 45:467–92. doi: 10.1146/annurev-soc-073018-022616

Crossref Full Text | Google Scholar

16. Steinmetz JD, Culbreth GT, Haile LM, Rafferty Q, Lo J, Fukutaki KG, et al. Global, regional, and national burden of gout, 1990-2020, and projections to 2050: a systematic analysis of the Global Burden of Disease Study 2021. Lancet Rheumatol. (2024) 6:e507–17. doi: 10.1016/S2665-9913(24)00117-6, PMID: 38996590

PubMed Abstract | Crossref Full Text | Google Scholar

17. Iannuzzi JP, King JA, Leong JH, Quan J, Windsor JW, Tanyingoh D, et al. Global incidence of acute pancreatitis is increasing over time: A systematic review and meta-analysis. Gastroenterology. (2022) 162:122–34. doi: 10.1053/j.gastro.2021.09.043, PMID: 34571026

PubMed Abstract | Crossref Full Text | Google Scholar

18. Qiu H, Cao S, and Xu R. Cancer incidence, mortality, and burden in China: a time-trend analysis and comparison with the United States and United Kingdom based on the global epidemiological data released in 2020. Cancer Commun (London England). (2021) 41:1037–48. doi: 10.1002/cac2.12197, PMID: 34288593

PubMed Abstract | Crossref Full Text | Google Scholar

19. Huang D, Lai H, Shi X, Jiang J, Zhu Z, Peng J, et al. Global temporal trends and projections of acute hepatitis E incidence among women of childbearing age: Age-period-cohort analysis 2021. J infection. (2024) 89:106250. doi: 10.1016/j.jinf.2024.106250, PMID: 39181413

PubMed Abstract | Crossref Full Text | Google Scholar

20. Cheng X, Yang Y, Schwebel DC, Liu Z, Li L, Cheng P, et al. Population ageing and mortality during 1990–2017: a global decomposition analysis. PLoS Med. (2020) 17:. doi: 10.1371/journal.pmed.1003138, PMID: 32511229

PubMed Abstract | Crossref Full Text | Google Scholar

21. Luo H. Global burden and cross-country inequalities in six major immune-mediated inflammatory diseases from 1990 to 2021: A systemic analysis of the Global Burden of Disease Study 2021. Autoimmun Rev. (2024) 23:103639. doi: 10.1016/j.autrev.2024.103639, PMID: 39260756

PubMed Abstract | Crossref Full Text | Google Scholar

22. Riebler A and Held L. Projecting the future burden of cancer: Bayesian age-period-cohort analysis with integrated nested Laplace approximations. Biometrical J Biometrische Z. (2017) 59:531–49. doi: 10.1002/bimj.201500263, PMID: 28139001

PubMed Abstract | Crossref Full Text | Google Scholar

23. Bender RG, Sirota SB, Swetschinski LR, Dominguez RMV, Novotney A, Wool EE, et al. Global, regional, and national incidence and mortality burden of non-COVID-19 lower respiratory infections and aetiologies, 1990-2021: a systematic analysis from the Global Burden of Disease Study 2021. Lancet. (2024) 24:974–1002. doi: 10.1016/S1473-3099(24)00176-2, PMID: 38636536

PubMed Abstract | Crossref Full Text | Google Scholar

24. Sha R, Kong XM, Li XY, and Wang YB. Global burden of breast cancer and attributable risk factors in 204 countries and territories, from 1990 to 2021: results from the Global Burden of Disease Study 2021. Biomark Res. (2024) 12:87. doi: 10.1186/s40364-024-00631-8, PMID: 39183342

PubMed Abstract | Crossref Full Text | Google Scholar

25. Australian Institute of Health and Welfare. Cancer in Australia: actual incidence and mortality data from 1982 to 2007 and projections to 2010. Asia-Pacific J Clin Oncol. (2011) 7:325–38. doi: 10.1111/j.1743-7563.2011.01502.x, PMID: 22151981

PubMed Abstract | Crossref Full Text | Google Scholar

26. Renard F, Van Eycken L, and Arbyn M. High burden of breast cancer in Belgium: recent trends in incidence (1999-2006) and historical trends in mortality (1954-2006). Arch Public Health = Arch belges sante publique. (2011) 69:2. doi: 10.1186/0778-7367-69-2, PMID: 22958447

PubMed Abstract | Crossref Full Text | Google Scholar

27. Gnerlich JL, Deshpande AD, Jeffe DB, Seelam S, Kimbuende E, and Margenthaler JA. Poorer survival outcomes for male breast cancer compared with female breast cancer may be attributable to in-stage migration. Ann Surg Oncol. (2011) 18:1837–44. doi: 10.1245/s10434-010-1468-3, PMID: 21484520

PubMed Abstract | Crossref Full Text | Google Scholar

28. Yeung WJ and Lee Y. Aging in east asia: new findings on retirement, health, and well-being. J Gerontology. (2022) 77:589–91. doi: 10.1093/geronb/gbab055, PMID: 33788932

PubMed Abstract | Crossref Full Text | Google Scholar

29. Rong J, Cheng P, Li D, Wang X, and Zhao D. Global, regional, and national temporal trends in prevalence for depressive disorders in older adults, 1990-2019: An age-period-cohort analysis based on the global burden of disease study 2019. Ageing Res Rev. (2024) 100:102443. doi: 10.1016/j.arr.2024.102443, PMID: 39097004

PubMed Abstract | Crossref Full Text | Google Scholar

30. Hu K, Ding P, Wu Y, Tian W, Pan T, and Zhang S. Global patterns and trends in the breast cancer incidence and mortality according to sociodemographic indices: an observational study based on the global burden of diseases. BMJ Open. (2019) 9:e028461. doi: 10.1136/bmjopen-2018-028461, PMID: 31594871

PubMed Abstract | Crossref Full Text | Google Scholar

31. Zhao H, Zhang Y, Liu H, Wang Y, and Song Z. Age-period-cohort analysis of global, regional, and national pancreatic cancer incidence, mortality, and disability-adjusted life years, 1990-2019. BMC Cancer. (2024) 24:1063. doi: 10.1186/s12885-024-12835-0, PMID: 39198814

PubMed Abstract | Crossref Full Text | Google Scholar

32. Mesa-Eguiagaray I, Wild SH, Rosenberg PS, Bird SM, Brewster DH, Hall PS, et al. Distinct temporal trends in breast cancer incidence from 1997 to 2016 by molecular subtypes: a population-based study of Scottish cancer registry data. Br J Cancer. (2020) 123:852–9. doi: 10.1038/s41416-020-0938-z, PMID: 32555534

PubMed Abstract | Crossref Full Text | Google Scholar

33. Breekveldt ECH, Toes-Zoutendijk E, Spaander MCW, van de Schootbrugge-Vandermeer HJ, van Vuuren AJ, van Kemenade FJ, et al. Advanced-stage CRC incidence patterns following the phased implementation of the CRC screening programme in the Netherlands. Eur J Cancer (Oxford Engl. (2023) 1990) 178:60–7. doi: 10.1016/j.ejca.2022.10.016, PMID: 36403368

PubMed Abstract | Crossref Full Text | Google Scholar

34. Lu G, Li J, Wang S, Pu J, Sun H, Wei Z, et al. The fluctuating incidence, improved survival of patients with breast cancer, and disparities by age, race, and socioeconomic status by decade, 1981-2010. Cancer Manage Res. (2018) 10:4899–914. doi: 10.2147/CMAR.S173099, PMID: 30464592

PubMed Abstract | Crossref Full Text | Google Scholar

35. Cook MB, Guénel P, Gapstur SM, van den Brandt PA, Michels KB, Casagrande JT, et al. Tobacco and alcohol in relation to male breast cancer: an analysis of the male breast cancer pooling project consortium. Cancer epidemiology Biomarkers Prev. (2015) 24:520–31. doi: 10.1158/1055-9965.EPI-14-1009, PMID: 25515550

PubMed Abstract | Crossref Full Text | Google Scholar

36. Guénel P, Cyr D, Sabroe S, Lynge E, Merletti F, Ahrens W, et al. Alcohol drinking may increase risk of breast cancer in men: a European population-based case-control study. Cancer causes control: CCC. (2004) 15:571–80. doi: 10.1023/B:CACO.0000036154.18162.43, PMID: 15280636

PubMed Abstract | Crossref Full Text | Google Scholar

37. Hamajima N, Hirose K, Tajima K, Rohan T, Calle EE, Heath CW Jr., et al. Alcohol, tobacco and breast cancer--collaborative reanalysis of individual data from 53 epidemiological studies, including 58,515 women with breast cancer and 95,067 women without the disease. Br J Cancer. (2002) 87:1234–45. doi: 10.1038/sj.bjc.6600596, PMID: 12439712

PubMed Abstract | Crossref Full Text | Google Scholar

38. Xu Y, Gong M, Wang Y, Yang Y, Liu S, and Zeng Q. Global trends and forecasts of breast cancer incidence and deaths. Sci Data. (2023) 10:334. doi: 10.1038/s41597-023-02253-5, PMID: 37244901

PubMed Abstract | Crossref Full Text | Google Scholar

39. Roheel A, Khan A, Anwar F, Akbar Z, Akhtar MF, Imran Khan M, et al. Global epidemiology of breast cancer based on risk factors: a systematic review. Front Oncol. (2023) 13:1240098. doi: 10.3389/fonc.2023.1240098, PMID: 37886170

PubMed Abstract | Crossref Full Text | Google Scholar

40. Zou S and Li C. Global and regional trends in male breast cancer burden: a comparative analysis of incidence, prevalence, mortality, and disability-adjusted life years in China and worldwide (1990-2021). Gland Surg. (2025) 14:818–33. doi: 10.21037/gs-2025-10, PMID: 40546833

PubMed Abstract | Crossref Full Text | Google Scholar

41. Gao Y, Goldberg JE, Young TK, Babb JS, Moy L, and Heller SL. Breast cancer screening in high-risk men: A 12-year longitudinal observational study of male breast imaging utilization and outcomes. Radiology. (2019) 293:282–91. doi: 10.1148/radiol.2019190971, PMID: 31526252

PubMed Abstract | Crossref Full Text | Google Scholar

42. Zhang H, Lin W, Chen D, Wang K, Tu W, Lin H, et al. Cardiovascular and other competing causes of death in male breast cancer patients: A population-based epidemiologic study. Clin Interventions Aging. (2021) 16:1393–401. doi: 10.2147/CIA.S314689, PMID: 34321871

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: breast cancer, elderly male, global burden, incidence, mortality

Citation: Liu S, Xu J, Zhang H, Tang Y, Han X, Xiong W, Fan L, Su C, Sui Z and Zhao W (2026) Global trends and disease burden of elderly male breast cancer, 1990-2021: a population-based study. Front. Endocrinol. 17:1674679. doi: 10.3389/fendo.2026.1674679

Received: 28 July 2025; Accepted: 12 January 2026; Revised: 05 December 2025;
Published: 06 February 2026.

Edited by:

Antonino Belfiore, University of Catania, Italy

Reviewed by:

Muhammad Imran Khan, Pak-Austria Fachhochschule Institute of Applied Sciences and Technology, Pakistan
Yiyuan Liu, The second Affiliated Hospital of Shantou University, China

Copyright © 2026 Liu, Xu, Zhang, Tang, Han, Xiong, Fan, Su, Sui and Zhao. 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: Wenhui Zhao, MDE4N0BocmJtdS5lZHUuY24=

These authors have contributed equally to this work and share first authorship

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.