Skip to main content

National and subnational burden of female and male breast cancer and risk factors in Iran from 1990 to 2019: results from the Global Burden of Disease study 2019

A Correction to this article was published on 15 June 2023

This article has been updated

Abstract

Background

Breast cancer (BC) is one of the most burdensome cancers worldwide. Despite advancements in diagnostic and treatment modalities, developing countries are still dealing with increasing burdens and existing disparities. This study provides estimates of BC burden and associated risk factors in Iran at the national and subnational levels over 30 years (1990–2019).

Methods

Data on BC burden for Iran were retrieved from the Global Burden of Disease (GBD) study from 1990 to 2019. GBD estimation methods were applied to explore BC incidence, prevalence, deaths, disability-adjusted life years (DALYs), and attributable burden to risk factors based on the GBD risk factors hierarchy. Moreover, decomposition analysis was performed to find the contribution of population growth, aging, and cause-specific incidence in the total incidence change. Age-standardized rates (per 100,000 population) and 95% uncertainty intervals (UI) were reported based on sex, age, and socio-demographic index (SDI).

Results

Age-standardized incidence rate (ASIR) increased from 18.8 (95% UI 15.3–24.1)/100,000 in 2019 to 34.0 (30.7–37.9)/100,000 in 2019 among females and from 0.2/100,000 (0.2–0.3) to 0.3/100,000 (0.3–0.4) among males. Age-standardized deaths rate (ASDR) increased slightly among females from 10.3 (8.2–13.6)/100,000 in 1990 to 11.9 (10.8–13.1)/100,000 in 2019 and remained almost the same among males—0.2/100,000 (0.1–0.2). Age-standardized DALYs rate also increased from 320.2 (265.4–405.4) to 368.7 (336.7–404.3) among females but decreased slightly in males from 4.5 (3.5–5.8) to 4.0 (3.5–4.5). Of the 417.6% increase in total incident cases from 1990–2019, 240.7% was related to cause-specific incidence. In both genders, the BC burden increased by age, including age groups under 50 before routine screening programs, and by SDI levels; the high and high-middle SDI regions had the highest BC burden in Iran. Based on the GBD risk factors hierarchy, high fasting plasma glucose (FPG) and alcohol were estimated to have the most and the least attributed DALYs for BC among females, respectively.

Conclusions

BC burden increased from 1990 to 2019 in both genders, and considerable discrepancies were found among different provinces and SDI quintiles in Iran. These increasing trends appeared to be associated with social and economic developments and changes in demographic factors. Improvements in registry systems and diagnostic capacities were also probably responsible for these growing trends. Raising general awareness and improving screening programs, early detection measures, and equitable access to healthcare systems might be the initial steps to tackle the increasing trends.

Background

Cancer is responsible for a majority of deaths due to non-communicable diseases (NCDs) and is among the main barriers to increasing life expectancy and declining mortality rates in numerous nations worldwide [1, 2]. Breast cancer (BC) is one of the most burdensome cancers globally, responsible for most cancer deaths globally. BC is estimated to be the most frequently diagnosed cancer, considering both sexes, and it is the first leading cause of cancer death among females [3]. Although BC in males accounts for approximately 1% of all BCs globally (diagnosed in 1:1000 of males), its incidence seems to be increasing [4, 5].

Among all malignancies diagnosed, BC is ranked first in Iranian females, with an estimated age-standardized incidence rate (ASIR) of 44.0 (95% uncertainty interval [UI]: 36.4–52.0) 100,000 population in 2016, based on the National and Subnational Burden of Diseases (NASBOD) project report on BC [6]. Estimates from prior studies have revealed increasing trends in BC incidence and mortality rates in Iran and the globe. Recent studies predicted an increase of about 63% in the new BC cases in Iran in the year 2025 compared with 2016, as the leading cancer among Iranian females [7]. Based on previous reports, BC seems to involve Iranian males with a higher frequency than the world’s average (approximately 2.8% of primary BC cases in the country) [8]. Besides, the 5-year survival rate of BC was estimated to be relatively lower than in developed countries [9]. As a developing country, Iran has been in a transitional period and has faced numerous challenges due to modernization over recent decades [10]. The World Bank categorized Iran as a lower-middle-income and upper-middle-income country in 1990 and 2019, respectively [11]. Economic developments and demographic changes over recent years are believed to be the major causes of the growing trends in the burden of cancers [2, 12,13,14]. Geographic and climate diversities, cultural variations, genetic factors, lifestyles, socioeconomic status, and population aging and growth might also be potential factors leading to increasing trends and disparities in burden among provinces in Iran [2, 15].

Accordingly, it is worthwhile to estimate and analyze BC trends in Iran at national and subnational levels to have a clearer vision of the burden of BC on the Iranian population and public health. The most updated estimates around BC in Iran were limited to a recent study by Ataenia et al. up to 2016. The authors reported incidence, mortality, and years of life lost (YLLs) for females at the national and subnational levels [16]. In the current study, we aimed to report the most up-to-date estimates of BC burden in Iran according to the Global Burden of Disease (GBD) study 2019 by estimating all measures related to the burden of disease, including incidence, prevalence, mortality, YLLs, years lived with disability (YLDs), and disability-adjusted life year (DALYs), at the national and subnational level for both females and males, based on age groups and socio-demographic index (SDI) and to report the burden attributable to known risk factors over thirty years from 1990 to 2019. This analysis will also provide policymakers and national authorities with a broader vision and a whole picture, highlighting possible disparities in different regions of Iran as an upper-middle income country in the middle-east region, which helps them tailor up-to-date strategies and allocate resources more efficiently.

Methods

Overview of the GBD study

The GBD study provides comprehensive and comparable health estimates around the burden of diseases and injuries worldwide. The GBD study 2019 (the most up-to-date iteration of the GBD study) estimates the burden of 369 causes of death and disability and 87 risk factors and groups of risk factors at the global and regional levels for 204 countries and territories among females and males [17]. We obtained publicly available data from the Global Health Data Exchange (GHDx) query online tool (http://ghdx.healthdata.org/gbd-results-tool) regarding the burden of BC, including annual all-ages numbers and ASIR, age-standardized deaths rate (ASDR) and prevalence rate, YLLs, YLDs, and DALYs over 30 years, from 1990 to 2019, in Iran at the national and subnational levels. (All the provinces of Iran are shown on the national map in the Additional file 1: Fig. 1.) The source of data from scientific papers and disease registry systems that have been used by the GBD study 2019 to estimate the burden of BC in Iran and subnational districts is provided in Additional file 2: Table 1. International Classification of Diseases 10 (ICD-10) codes (C50-C50.629, C50.7, C50.8-C50.929, D05-D05.92, D24.0-D24.9, D48.6-D48.62, D49.3, Z12.3-Z12.39, Z80.3, Z85.3, and Z86.000) were utilized to identify BC, which is defined as the code B.1.14 in the GBD 2019 list of causes of death and disability. This study follows the Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER) [18].

Fig. 1
figure 1

Time trends of breast cancer age-standardized incidence, prevalence, deaths and disability-adjusted life years (DALYs), rates (per 100,000 population) and numbers and 95% UIs in Iran, from 1990 to 2019 for both sexes

Table 1 Breast cancer incidence, prevalence, deaths, disability-adjusted life years (DALYs), years of life lost (YLLs), and years lived with disability (YLDs) by rates (per 100,000 population) and numbers with percent changes in Iran between 1990 and 2019, for both sexes, females, and males

Cancer burden estimation

The GBD study employs various estimation and modeling tools to provide comparable results around the global burden of diseases and injuries. The detailed process of burden estimations regarding primary values, including incidence, prevalence, mortality, YLLs, YLDs, and DALYs in the GBD 2019 study, has been described in detail elsewhere [17]. Here, we briefly describe the estimation methods of these primary values by the GBD study.

In some regions, cancer mortality data are scarce, while cancer incidence data are more accessible through population-based cancer registries. So, to maximize data informing mortality models, the GBD study transforms incidence data into cancer mortality estimates using modeled mortality-to-incidence ratios (MIRs). First, incidence and mortality data obtained from cancer registries (inputs) were processed to generate crude MIRs. Final cancer cause-specific MIRs were estimated using a spatiotemporal Gaussian process regression (ST-GPR) by utilizing Healthcare Access and Quality (HAQ) index, age, and sex as covariates. The MIR estimates using the ST-GPR model were then multiplied with incidence data to generate crude mortality estimates. Subsequently, MIR-transformed mortality data were pooled with mortality data obtained from vital registries and verbal autopsies and were used as inputs in cancer-specific Cause of Death Ensemble models (CODEm) [19,20,21]. CODEm produced models of deaths for each cause, sex and age group, and time point within the GBD study by selecting models based on out-of-sample predictive validity. Eventually, the predicted mortality estimates were adjusted to align with independently modeled estimates for each sex-age group-location-year [17].

Incidence estimates were generated from modeled cancer mortality estimates using BC MIRs for each sex, age group, location, and year. Prevalence of BC was modeled using incidence, mortality background, and estimated relative survival curves and their correlation with modeled MIRs. For the estimation of YLDs, the estimated prevalence was multiplied by the corresponding BC-specific disability weight. Disability weight for each cause describes the severity of health loss with that specific cause, ranging from 0 (full health) to 1 (death) [17]. Accordingly, YLLs were calculated by multiplying the difference between the number of deaths due to BC in a specific age group by the remaining standard life expectancy at the age of death. DALYs also represent the sum of YLDs and YLLs [17, 22].

Socio-demographic index (SDI)

We categorized Iranian subnational districts by SDI into five quintiles (low, low-middle, middle, high-middle, and high SDI). We calculated this index according to the previously described methods based on the lag distributed income per capita, average educational attainment (schooling years) for individuals over 15 years, and the fertility rate in females under 25 years, which could be found elsewhere [23].

Decomposition analysis

We performed decomposition analysis in this study to find the attribution of population growth, population aging, and BC incidence rate in the total change of BC incident cases between 1990 and 2019. We performed this analysis based on the previously described method, which could be found elsewhere [24]. Briefly, two steps were taken to analyze this attribution: in the first step, the age structure, sex structure, and BC rates of 1990 were applied to the total population of 2019. The difference between this step and the total number of new cases in 1990 was attributed to population growth. In the second step, the BC incidence rate from 1990 was applied to the age structure, sex structure, and population size of 2019. The difference between these two steps was attributed to population aging. Eventually, differences between the total number of incident cases in 2019 and the second step were attributed to true changes in the BC-specific incidence rate.

Mortality-to-incidence ratio (MIR)

To report the age-standardized MIRs in this study, we divided ASDR by ASIR, obtained from the GBD results for BC, for each sex-age group-location-year combination as follows [25]:

$$\mathrm{MIR} \, (\mathrm{Mortality\!-\!to\!-\!Incidence\,Ratio})=\frac{\mathrm{Age\!-\!standardized\,deaths\,rate} \, (\mathrm{ASDR})}{\mathrm{Age\!-\!standardized\,incidence\,rate} \, (\mathrm{ASIR})}$$

To calculate the MIRs, we used the final estimated values for incidence and mortality (outputs) after the modeling systems utilized by the original GBD estimation processes described earlier in the methods section. The MIR calculation was validated previously for BC in the country using local data as part of the NASBOD project and at the global scale using GBD 2019 data [6, 26].

Estimation of the breast cancer burden attributable to risk factors

The GBD study has incorporated a comparative risk assessment framework comprising a risk factor hierarchy categorized into three main groups (environmental and occupational, behavioral, and metabolic risks) and four levels to assess the attributable burden to different risk factors. Counting all specific risk factors and aggregates, the GBD study 2019 includes 87 risks or clusters of risks [27]. Risk factors for BC were included in this study based on the GBD risk factor hierarchy at the 4th level, as follows: alcohol use, smoking, secondhand smoke, high fasting plasma glucose (FPG) [28], diet high in red meat, high body-mass index (BMI), and low physical activity. These risk factors have been included in the last iteration of the GBD study (2019) as the risk factors known to be associated with BC; for more details on the inclusion process of the risk factors, refer to the original GBD 2019 risk factors summary paper [27].

Detailed methods of estimating risk factors attributed burden could be found elsewhere [27]. Briefly, the GBD comparative risk assessment framework followed a six-step approach to compute the fraction of cancer-specific burden attributable to risk factors in the GBD study 2019. In the first step, risk factors with convincing or probable evidence for a causal association were identified using the World Cancer Research Fund criteria [29]. In the next step, relative risks (RRs) were estimated for each risk-outcome pair as a function of exposure. The meta-analytic approach was updated for a selected set of continuous risk factors in the GBD study 2019, using the meta-regression-Bayesian, regularized, trimmed (MR-BRT) mixed-effects model. In the third step, exposure levels and distributions of risk factors were modeled for each sex-age-location-year combination using existing data from published studies, household surveys, censuses, administrative data, ground monitor data, or remote sensing data. The GBD study modeled risk factor exposure levels using either a Bayesian meta-regression modeling tool (DisMod-MR 2.1), a flexible approach incorporating age- and sex-specific data, or ST-GPR, which was the preferred approach when exposure was stable across different age groups. In the fourth step, the theoretical minimum risk exposure level (TMREL) was identified for each risk factor. In the next step, the population attributable fraction (PAF) was computed for each risk-outcome pair for each sex-age-location-year based on the RRs, exposure levels, and the TMREL. Then, correction for PAF overestimation, which might occur if the independence of certain risk factors was assumed, was performed using a mediation matrix. In the last step, to estimate the BC burden attributable to each risk factor, YLLs, YLDs, and deaths were multiplied by the corresponding PAF of each risk factor for each sex-age-location-year combination. The sum of YLLs and YLDs was used to estimate the DALYs attributable to risk factors [27, 30, 31].

Statistical analysis

Age groups for BC were defined from 15 to + 80 years on a 5-year basis. Crude all age numbers for the primary values have been reported in this study, and the rates have been reported per 100,000 population. Age-standardized rates were estimated using the direct method of standardization and the GBD world population standard and were reported per 100,000 population [22]. We calculated and reported percent changes (PCs) for all the measures between 1990 and 2019. Annual percent changes (APC) were also generated to measure the trend of the age-standardized rates changes over the study period. To this end, we applied a regression line to the natural logarithm of the age-standardized rates, and 95% confidence intervals were also reported. Besides, we performed Pearson’s correlation tests to evaluate the correlation between different burden measures and SDI and between females and males in different provinces regarding age-standardized rates for each measure. We performed all analyses and data visualizations using Stata version 13 (StataCorp LLC, TX, USA) and R version 3.5.0; statistical codes are publicly available elsewhere [17]. For all estimates, the 95% UI is reported by obtaining 1000 samples of the posterior distribution and marking the 25th and 975th ranked values across all 1000 draws, which covers the proportion of data point mean values that fall between the 2.5th and 97.5th percentile of the draws of the fit values [17].

Results

Breast cancer incidence and prevalence in Iran

Overall, BC incidence case number among females increased by 5.2-fold from 2,835 (95% UI: 2352–3575) in 1990 to 14,743 (13,248–16,469) in 2019 in Iran. Similarly, ASIR for females increased from 18.8 (15.3–24.1)/100,000 to 34.0 (30.7–37.9)/100,000, with 81.2% (34.6–130.5) PC and 1.87 (1.69–2.05) APC in the same period (Table 1 and Additional file 3: Table 2). Among males, the number of incident cases also increased from 36 (27–48) to 120 (98–143) over this period, with 28.6% (− 13.2–80.4) PC and 0.39 (0.05–0.74) APC in their ASIR. Males comprised approximately 1.2% of the total incidence in 1990 (36 [27–48]/2871 [2389–3608]), while approximately 0.8% of all incident cases in 2019 (120 [98–143]/14,863 [13,372–16,590]) were male BC. Annual trends illustrated in Fig. 1 show ASIR trends and the all-ages number of incidence from 1990 to 2019. The overall age-standardized prevalence rate of BC in Iran was estimated as 312.e1 (281.5–345.4)/100,000 in 2019 for females, which was remarkably higher than its rate in 1990, estimated as 178.6 (149.7–216.4)/100,000. For males, this rate increased from 2.0 (1.6–2.6)/100,000 to 2.8 (2.3–3.3)/100,000 over this period.

Table 2 Decomposition analysis on changes in breast cancer incident cases in Iran and its 31 provinces, for both sexes, females, and males, between 1990 and 2019 (provinces have been sorted alphabetically)

At the subnational level, all provinces in Iran experienced a remarkable increase in ASIR for both females and males (Fig. 2A, B, and C), and the average ASIR, considering both sexes, almost doubled over the study period—Iran’s provinces are illustrated on the national map in Additional file 1: Fig. 1. In 1990, the highest and lowest ASIR for females among all provinces were 30.3 (21.0–43.0)/100,000 (Tehran) and 10.4 (6.6–16.2)/100,000 (Kohgiluyeh and Boyer-Ahmad), respectively, while in 2019, the highest and lowest rates were estimated as 41.8 (33.3–52.6)/100,000 (Alborz) and 20.0 (15.5–25.7)/100,000 (Sistan and Baluchistan), respectively. Taking both sexes into consideration, the maximum and minimum ASIR were 14.9 (10.4–21.2)/100,000 (Tehran) and 5.1 (3.3–7.9)/100,000 (Kohgiluyeh and Boyer-Ahmad) in 1990 and 20.8 (16.5–25.8)/100,000 (Mazandaran) and 10.1 (7.8–12.9)/100,000 (Sistan and Baluchistan) in 2019, respectively, see Additional file 4: Fig. 2A, Additional file 5: Fig. 2B, and Additional file 6: Fig. 2C. Almost all provinces experienced this remarkable increase in ASIR for females and males (Additional file 7: Table 3).

Fig. 2
figure 2figure 2

Breast cancer age-standardized incidence, prevalence, deaths, and disability-adjusted life years (DALYs) rates (per 100,000) and their changes between 1990 and 2019 in Iran at the subnational level for A both sexes, B females, and C males

In order to evaluate changes in ASIR from 1990 to 2019, PCs were assessed in two major periods, based on the two main periods over the original GBD study 2019 [17]: 1) between 1990 and 2010; 2) between 2010 and 2019 (Fig. 3A and Additional file 8: Table 4). In Tehran, the capital of Iran, the total (both sexes) PC of ASIR in the second period was almost sixfold higher than the first period, with the overall ASIR PC of 35.5% (− 12.9–111.0) from 1990 to 2019—not statistically significant increase. Regarding sexes, no province at the subnational level had a negative ASIR PC in females in both periods, while some provinces had negative ASIR PCs in males during the second period, illustrated in Fig. 3A. The highest negative ASIR PC for males was − 7.1% (− 47.9–75.2; Sistan and Baluchistan). On the other hand, the highest positive ASIR PC for females was 131.8% (30.5–310.2; Kohgiluyeh and Boyer-Ahmad), and Ilam had the highest APC in incidence rate of 3.31 (3.11–3.51) among all provinces between 1990 and 2019 (Additional file 3: Table 2).

Fig. 3
figure 3

A Breast cancer age-standardized incidence, prevalence, deaths, and disability-adjusted life years (DALYs) rates percent changes (%) over two major periods (X-axis: the first period, from 1990 to 2010; Y-axis: the second period, from 2010 to 2019) in Iran and its 31 provinces for both sexes, females, and males. The oblique line in each graph shows the equal percent change in the second period as the first period. B Breast cancer age-standardized incidence, prevalence, deaths, and disability-adjusted life years (DALYs) rates (per 100,000 population) based on the subnational socio-demographic index (SDI) quintiles in 1990 and 2019 for both sexes. The dotted line shows the mean rate within each SDI quintile

Moreover, the highest age-standardized prevalence rates at the subnational level in 2019 were estimated at 383.7 (311.1–470.1)/100,000 (Alborz) for females and 6.9 (5.0–10.0)/100,000 (Yazd) for males, while the lowest age-standardized prevalence rates for females and males were 179.6 (146.6–223.7)/100,000 and 1.8 (1.3–2.4)/100,000, respectively (both in Sistan and Baluchistan); see Additional file 7: Table 3.

Table 3 Breast cancer incidence, prevalence, deaths, disability-adjusted life years (DALYs), years of life lost (YLLs), and years lived with disability (YLDs) by rates (per 100,000 population) and numbers with percent changes in Iran between 1990 and 2019, for both sexes, females, and males, based on the age groups < 50 (15–49), 50–69, and + 70

Regarding SDI quintiles at the subnational level, ASIR increased with SDI in both 1990 and 2019. The mean ASIR for high and high-middle SDI was approximately doubled from 1990 to 2019. Low, low-middle, and middle SDI quintile regions also experienced this remarkable increase, while mean ASIRs in these regions were overall lower than that of high-middle and high-SDI regions in 1990 and 2019 (Fig. 3B).

Decomposition analysis revealed that 133.0% and 44.0% of the total 417.6% change in incident cases from 1990 to 2019 were attributable to age structure change and population growth, respectively. Thus, it was estimated that 240.7% of the total increase in new cases was related to the true increase in the BC-specific incidence over this period. All the decomposition analyses for both females and males at the national and subnational levels are presented in Table 2.

The breast cancer burden of deaths and DALYs in Iran

BC led to 4704 (4306–5192) deaths among females in 2019, almost 3.3-fold more than 1990 with 1413 (1160–1796) deaths at the national level. ASDR for females also increased slightly over these 30 years from 10.3 (8.2–13.6)/100,000 to 11.9 (10.8–13.1)/100,000, which had a PC of 14.9% (− 15.1–47.7) and 0.50 (0.34–0.66) APC. Moreover, the age-standardized DALYs rate of females was 368.7 (336.7–404.3)/100,000 in 2019, 15.2% (− 11.7–42.5) higher than of 1990 (320.2 [265.4–405.4]/100,000), 0.51 (0.36–0.65) APC; however, the overall PCs in age-standardized DALYs and ASDR in males decreased over this period (-10.7% [− 36.1–20.7] and − 11.4% [− 37.0–22.6], respectively). It is worth mentioning that 93.8% of the overall age-standardized DALYs rate in both sexes was related to YLLs in 2019, while this proportion was 96.0% in 1990 (Table 1).

At the subnational level, ASDR and age-standardized DALYs rate increased remarkably in most provinces of Iran; these increasing rates were mostly estimated for central and northern parts of Iran, as illustrated in Figs. 2A, B, and C. Females had positive ASDR PCs in most provinces, while males mostly had negative ASDR PCs in Iran (Fig. 3A). Changes in age-standardized DALYs rate also revealed almost a similar pattern to ASDR in most provinces for both sexes (Fig. 3A). Based on Fig. 3A, two provinces (Ilam and Zanjan) had the most remarkable deviation from the others and the overall PCs in age-standardized DALYs and ASDRs in the country and had the highest positive PCs in the first period; however, they had PCs lower than average in the second period. Ilam had the highest APCs in age-standardized DALYs and ASDR between 1990 and 2019, 1.75 (1.63–1.87) and 1.80 (1.63–1.96), respectively (Additional file 3: Table 2).

All the provinces’ rankings in terms of ASDR for both sexes and females and males separately are presented in Additional file 4: Fig. 2A, Additional file 5: Fig. 2B, and Additional file 6: Fig. 2C. Among females, the highest estimated ASDR in 1990 was 14.3 (9.9–20.8)/100,000 (Tehran). The first-ranked province in 2019 had, similarly, an ASDR of 14.3 (11.6–17.2)/100,000 (Qom); however, the lowest ASDR for females was estimated at 7.5 (5.7–9.7)/100,000 (Chahar Mahaal and Bakhtiari) in 2019 compared to 6.4 (4–10.2)/100,000 (Kohgiluyeh and Boyer-Ahmad) in 1990 (Additional file 5: Fig. 2B). Males had ASDRs of approximately 0.5/100,000 or lower in 1990 and 2019 (Additional file 6: Fig. 2C). Considering YLLs and YLDs, both rates increased remarkably in Iran and almost all provinces; however, only limited provinces experienced decreasing YLLs while having increasing YLDs from 1990 to 2019, such as Tehran and Chahar Mahaal and Bakhtiari (Fig. 4A).

Fig. 4
figure 4

A Breast cancer age-standardized years of life lost (YLLs) and years lived with disability (YLDs) rates (per 100,000 population), in 1990 and 2019 for both sexes. (left: for Iran and its 31 provinces; right: based on age groups). B Breast cancer incidence, prevalence, deaths, and disability-adjusted life years (DALYs) rates (per 100,000 population) in 1990 and 2019 based on age groups by sex (red: females; blue: males) in Iran

Based on SDI quintiles, the ASDR and the age-standardized DALYs rate in all quintiles were higher in 2019 compared to 1990. Overall, high and high-middle SDI quintiles had slightly higher mean ASDRs and age-standardized DALYs rate than the others (Fig. 3B).

Age-related breast cancer burden

Based on age groups, both sexes’ incidence and deaths rates were higher in almost all age groups in 2019 than in 1990 at the national level. Overall, the incidence rate increased by age in both 1990 and 2019, and the highest incidence rate was estimated in the 75–79 age group in both 1990 and 2019. Females were diagnosed with BC at younger age groups than males who were diagnosed with BC, mostly over 40–45 years. Females over 80 also experienced the highest increase in the incidence rate over this period (Fig. 4B). Although ASDR among males did not change remarkably in different age groups, it increased considerably in females of all age groups over this period (Fig. 4B). DALYs rate also increased among Iranian females of all ages from 1990 to 2019, while it had a reverse U-pattern association with age, reaching a peak at around 1,144/100,000 in the 55–59 age group in 2019. The highest age-standardized YLDs rate was estimated for the 75–79 age group, followed by 65–69 years (Fig. 4A). Of note, men had a negative change regarding DALYs rate in most age groups from 1990 to 2019. We found that the + 70-age group had highest ASIR and ASDR compared to < 50 (15–49) and 50–69 age groups, while the 50–69 age group had the highest age-standardized DALYs rate in both 1990 and 2019. It is worth noting that the most remarkable increases in the age-standardized rates were found < 50 years (15–49 age group) compared to age groups 50–69 and + 70 years since the 15–49 age group had highest PCs between 1990 and 2019 despite the relatively lower rates in this age group (Table 3).

At the subnational level, almost all 31 provinces had a similar increasing pattern to the national trend of Iran for ASIR and ASDR by age groups, with the highest rates estimated for age groups over 70 years. Besides, the highest ASDR was estimated for the + 80-age group in all provinces. In almost all provinces, DALYs rate followed a similar pattern to Iran, with the highest rates estimated for middle-aged groups (Additional file 9: Fig. 3, Additional file 10: Fig. 4, Additional file 11: Fig. 5 and Additional file 12: Fig. 6).

Correlation between measures of burden and SDI

We found a high positive correlation between age-standardized rates of all burden measures and SDI in females (Additional file 13: Table 5), and the highest Pearson’s correlation coefficients (r) were calculated for incidence, prevalence, and YLDs (r > 0.80). On the contrary, males had low correlations between the measures and SDI, and the lowest negative correlation in males was found for death rates (r = − 0.05).

Correlation between females and males in different provinces regarding BC burden

We found positive correlations between females and males with reference to ASIR, ASDR, and age-standardized DALYs rate in different provinces in the year 2019 with r of 0.54, 0.43, and 0.40 respectively. However, there were no meaningful positive correlations between females and males for the age-standardized rates in 1990 (Fig. 5).

Fig. 5
figure 5

Pearson’s correlation coefficient (r) between males and females of all provinces of Iran regarding measures of BC burden

Mortality-to-incidence ratio (MIR)

Overall, age-standardized MIR for BC revealed a downward trend in Iran and all provinces for both females and males (Fig. 6). Generally, males had a higher MIR than females in Iran, but they followed almost a similar pattern in all provinces. The overall age-standardized MIR for BC in Iran was 0.35 in 2019 compared to 0.55 in 1990. At the subnational level, the highest and the lowest MIR in 2019 were estimated as 0.48 and 0.30, respectively (Sistan and Baluchestan and Mazandaran). In 1990, the highest and the lowest MIR were also recorded in the same provinces (0.73 and 0.48, respectively).

Fig. 6
figure 6

Age-standardized mortality-to-incidence ratio (MIR) trend in Iran and all provinces from 1990 to 2019

The breast cancer burden attributable to risk factors

Due to the scarce data around risk factors attributed burden in males, we have only provided the attributable burden to risk factors in females. Overall, deaths and DALYs rate attributable to each risk factor had upward trends in Iran and all SDI quintiles from 1990 to 2019; high and high-middle SDI quintiles generally had the highest deaths and DALYs rate attributable to risk factors, except for high BMI (Fig. 7A and B). High FPG exerted the highest burden among risk factors, followed by secondhand smoke, low physical activity, diet high in red meat, smoking, and alcohol use (Table 4). High FPG-related ASDR for females was 1.1 (0.2–2.4)/100,000 in 2019 with 117.1% (52.0–210.1) increment compared to 1990, and its age-standardized DALYs rate was 28.1 (5.5–61.0)/100,000 in 2019 (PC of 117.1% [60.6–194.2] from 1990 to 2019). Alcohol use, on the other hand, had the least attributable burden, with the age-standardized DALYs rate of 1.3 (0.9–1.8)/100,000 among Iranian females, followed by high-BMI with the attributed age-standardized DALYs rate of 3.2 (-11.3–16.5)/100,000 (Table 4). BC burden attributable to each risk factor at the subnational level can be found in Additional file 14: Table 6.

Fig. 7
figure 7

Time trend of all-ages breast cancer A deaths and B disability-adjusted life years (DALYs) rates (per 100,000) attributable to risk factors are shown for Iran and its subnational socio-demographic index (SDI) quintiles from 1990 to 2019 in females

Table 4 Breast cancer deaths and disability-adjusted life years (DALYs) by rates (per 100,000 population) and numbers with percent changes attributable to risk factors in Iran between 1990 and 2019 for females and both sexes

A total DALYs (all-ages numbers) of 5563 were attributable to all the risk factors mentioned above in 1990, with a fraction of about 11.0% of the total DALYs. The fraction of attributable burden to all the included risk factors was 13.8% in 2019 (22,593/163,091); see Tables 1 and 4.

Discussion

This study revealed that BC incidence, deaths, and DALYs generally had upward trends in Iran from 1990 to 2019. Although ASIR for males was increasing almost constantly over the study period, their ASDR and age-standardized DALYs had negative PCs from 1990 to 2019 in most age groups. In contrast, these trends increased for females of all ages during this period. Decomposition analyses revealed that 133.0% and 44.0% of the total 417.6% change in incident cases from 1990 to 2019 were attributable to age structure change and population growth, respectively, and 240.7% of the total increase in new cases was related to the true increase in the BC-specific incidence over this period. Besides, females faced BC at a younger age than males; however, the incidence of BC increased remarkably with aging in both sexes. Age-standardized MIRs had a downward trend from 1990 to 2019 in Iran and all provinces, while men had higher MIRs than females over the study period. Regarding associated risk factors, high FPG was found to exert the most attributed burden in Iran, contrary to alcohol as the least.

Concerning age and sexes, our study found that the incidence rate in females increased by age, and the highest incidence rate was estimated for the 75–79 age group. Despite the less significant estimates on the burden of BC in younger age groups, previous studies estimated that BC still affects 4–6% of females under 40 years [32]. We also found that the chance of BC accelerates with aging. Also, we found that the younger population under 50 years experienced the sharpest increase in the incidence rate of BC compared to older age groups, specifically women under 50 years had the highest magnitude in ASIR increase over the study period. This finding was in line with previous studies on female and male BC at the global level. Z Chen et al. found in a recent study that the magnitude of the increase in ASIR decrease by age, and females under 50 years had the highest increase in their rates despite the relatively lower ASIR among the population [33]. So, special attention should be paid to younger age groups focusing on improving general awareness, developing screening programs and sticking to them, and modulating preventable risk factors, which might help control the increasing burden of BC among younger population as well as advanced ages. Concerning the remarkable increase in BC incidence among younger age groups, special attention should also be paid to other possible risk factors in the nation, such as stress, socioeconomic situation, and lifestyle and environmental changes. It is worth investigating by future studies to explore the possible effect of such factors on the increasing burden of BC among younger age groups more in depth. Besides, male BC has been mostly underrepresented in studies and clinical trials due to the meager incidence rate of BC among men. Limited studies on male BC revealed that men were diagnosed with BC mostly aged > 40 years, and the median age at diagnosis in males was approximately five years higher than females, 67 and 62 years, respectively, and presented at more advanced stages with a lower survival rate [8, 34, 35], concordant with our results showing that BC was mostly diagnosed in men over 40–44 years. Our study showed that although the estimated percentage of new male BC cases out of all diagnosed cases decreased from 1.2% in 1990 to 0.8% in 2019, ASIR still increased in males over this period. This study also found that MIR was higher in males than in females. This might be due to the evidence that males are mostly diagnosed at more advanced ages and stages than females. Previous studies also reported a lower survival rate in males than females, which could explain the higher MIRs in males based on our findings [35]. So, males should not be neglected, and effective strategies should be tailored to diagnose male BC at earlier ages and stages and manage them appropriately. Further studies are needed focusing on the early detection measures and treatment modalities to tackle high MIRs in males. Identifying males at higher risk for BC and developing screening programs for such individuals might be beneficial based on previous studies’ evidence [36]. MIR is associated with the clinical outcome of cancers and could be utilized as an indicator of the quality of care provided by the healthcare system [37, 38]. Based on previous studies, the global MIR for female BC decreased from 0.4 in 1990 to 0.3 in 2016 [25]. In Iran, we also found a similar downward pattern for the MIR of BC for both females and males. We inferred that these decreasing trends might be due to improvements in early diagnosis and treatment capacities over these 30 years, leading to a better survival rate.

Regarding SDI quintiles, high and high-middle SDI regions in Iran were found to have greater incidence, deaths, and DALYs rate than the others. In contrast, previous studies revealed that ASDR and DALYs rate of BC decreased globally and in high and high-middle SDI regions of the world over time, and the low-middle SDI regions surpassed the high SDI over the 2010–2017 period [39]. Improvements in therapeutic approaches, including up-to-date surgical technologies, radiation therapies, and the use of novel systemic agents based on BC subtypes [40] in developed countries, might have led to the decreased ASDR and DALYs rate in such regions. Iran is still a middle-income developing country facing limitations in diagnostic and screening capacities and access to novel therapeutic modalities. Still, better diagnostic capacities, access to healthcare facilities, and higher levels of awareness in high-SDI regions of Iran seem to be responsible for higher estimates of BC burden in these regions than in lower SDIs. Besides, lifestyle-associated and environmental and occupational risk factors might be more significant in high-SDI regions due to modernization, such as diet high in red meat and low in vegetables and organic products, and toxic environmental pollutants. Studies proposed that industrial development has led to the soil, surface water, and food, as a consequence, pollution by heavy metal salts, which have been found to stimulate BC progression and reduce its sensitivity to treatment [41]. Results of the current study also found that the impact of PCs in ASIR and ASDR is declining over time (from 2010–2019 compared to the 1990–2010 period) in wealthier provinces of Iran, while PCs in more deprived provinces in recent years were still as high as previous decades. These findings were in line with the results of a recent investigation by Rahimzadeh et al [42], on the impact of geographical and socioeconomic inequalities in female BC incidence and mortality in Iran, suggesting that a possible reversal might occur in the upcoming years, with the most deprived provinces possessing the highest ASDRs in the country. So, more attention should be paid to low SDI regions in Iran to resolve possible inequalities and control the impact of increasing mortality rates and the burden of BC in these regions. Iran is a vast country, longitudinally and latitudinally, consisting of a wide range of ethnicities. In addition to the possible factors mentioned above, disparities in lifestyles, individuals’ behaviors, cultures, genetic backgrounds, the age structure of the population in each region, and geographical differences might correspond to the divergence in rates in different country regions [43,44,45,46]. It necessitates future studies aiming at close observations and evaluations to find out possible underlying causes of disparities and any possible inequalities in access to healthcare services and the quality of care.

According to the GBD risk factor hierarchy related to BC, we found that high-FPG exerted the most attributed deaths and DALYs among Iranian females, and its attributed burden had an upward trend over the study period. Globally, alcohol consumption seems to be the most important contributing factor to the BC burden [12, 47]. However, we found that alcohol use was a minor attributable factor among the Iranian population. It is worth noting that alcohol consumption is illegal in Iran, and the data regarding this risk factor might not be accurate. Previous investigations also reported that in countries with religious-based policies to keep the rate of alcohol consumption low, the lowest PAFs were recorded regarding this factor [48].

Regarding tobacco exposure, some studies found that secondhand (passive) smoking puts individuals at a higher risk of BC than active smoking [49, 50]. Accordingly, in the present study, we found that secondhand smoking was associated with deaths and DALYs rate more than active smoking. Besides, the association between BC and BMI is controversial and still under investigation [51]. Some studies showed a strong association between increased BMI and higher BC incidence, especially in postmenopausal women who are not taking exogenous hormones, and obesity has been found as a potential risk factor for BC. By contrast, in premenopausal women, high BMI is negatively associated with BC risk [52]. In the present study, we found that high-BMI-related deaths rate had an increasing trend over the study period, while a downward trend was observed regarding DALYs. This trend could be explained by the diverse role of obesity in BC among different age groups.

Nevertheless, lifestyle-associated risk factors should be regulated in the population through preventing obesity and not having a high-fat diet. Studies also revealed that daily and routine exercise and vegetable consumption could lower the risk of BC development [50]. In a recent meta-analysis concerning the association between nutrition and the risk of BC, the authors found that fruit, vegetables, and soy consumption had an inverse association with BC. However, red meat and processed meat both had a positive association with BC risk, while processed meat consumption showed a higher impact than red meat [53]. Therefore, adopting a healthy lifestyle, including a healthy diet, physical activity, and non-smoking behavior, might help prevent BC in the community setting [54, 55].

Estimates from the GBD study provide a good vision of the burden of diseases worldwide for clinicians and policymakers. Accordingly, the present study results could raise the health information, which can be helpful in adapting appropriate strategies to control the increasing burden of BC in Iran and similar developing countries worldwide. Previous studies suggested that developing a National Cancer Control Program (NCCP) might immensely improve the management of cancers in each country by mapping the existing problems and considering possible implementation solutions, especially in developing and low-income countries, which have been reported to have no effective NCCP [56]. Besides, recent studies in Iran revealed that around 61% of Iranian women know about BC and its screening programs, while only 17% stick to these programs [57]. So, mass education of the population, especially for individuals at a higher risk of BC, based on the geographical, genetic factors, and findings from burden studies could raise general awareness, leading to better prevention and early detection of cancers.

Limitations

Despite providing high-quality estimates of the burden associated with a wide range of diseases and injuries by the GBD study, some limitations might still exist [17]. The GBD study primarily relies on data acquired from vital and cancer registries and other epidemiologic studies worldwide, and the results mainly depend on the out-of-sample predictive validity of modeling efforts [17]. In many regions of the world, especially in low-income regions, cancer registries might not be fully developed and thus reliable enough, which could influence the final results [58, 59]. Early detection of new cases requires updated systems and advanced modalities that might not be effectively accessible in low-income regions and provinces. Based on previous reports, cancer registries began in Iran in 1999, while it did not cover all cancer data from pathology laboratories and departments all over the country during the first years, and it took years to be developed to a higher standard level [60]. So, data from that period in the GBD study mostly rely on epidemiological studies and systematic reviews. The GBD study adopted advanced estimation models and methods to tackle the scarcity of data in such regions as much as possible, as described briefly in the methods of the present study. Recent investigations by the NASBOD project also attempted to estimate and resolve possible incompleteness of cancer registry systems in Iran by adapting data from the Social Security Organization Cancer Registry (SSOCR), since before receiving cancer medications, all cancer patients should be registered in the SSOCR database [16]. Results from the NASBOD project estimating the burden of cancers could be accessed using online data visualization tools at www.vizit.report website.

There were also limitations regarding males, especially concerning the burden attributable to risk factors. There were very few BC cases among males in Iran over the study period. So, the interpretations of results for males in this study should be made and applied cautiously. We were also not able to investigate the burden of BC attributable to risk factors in males due to the limited number of cases in each category. So, we could report attributable burden to risk factors only for females. It is highly recommended to investigate the male BC burden and the potential risk factors in future studies more precisely.

Concerning existing disparities among different subnational regions of Iran, it is also vital to investigate access to routine screening programs and explore some key points, such as which regions have or do not have access to screening programs, to what extent the national screening programs for BC are performed in each region, and what percentage of the population among younger and older age groups are screened. Future studies could provide valuable information and help better policy making in the nation regarding BC burden control by answering such questions require, which requires precise investigations at the national and subnational level.

Some limitations also arise from the Iranian population, culture, and genetic background. For example, it is worth noting that alcohol consumption is illegal in Iran based on culture and religion. So, the reported data might be underestimated regarding alcohol consumption. This might influence the real burden attributable to alcohol use among the Iranian population. Possible genetic diversities among different ethnicities might influence the existing disparities regarding BC burden in different provinces, which could not be covered and investigated based on the GBD study results. Further studies are needed to focus on this aspect and investigate the role of genetic diversity in the future. We could investigate the BC burden attributable to risk factors based on the GBD study 2019 risk factors hierarchy. In this study, we found that only 11.0% of the total DALYs in 1990 and 13.8% in 2019 were attributable to the included risk factors, and about 86.2% of the total DALYs in 2019 were associated with other factors. There are still some potential risk factors known to be associated with BC, such as genetic background and gene mutations, hormonal patterns, and hazardous environmental pollutants, to name but a few. The GBD study constantly investigates the burden attributable to more risk factors in each updated iteration. So, we suggest that future iterations of the GBD study incorporate more known risk factors of BC to investigate their attributed burden.

Besides, the quality of care provided by healthcare systems is another vital public health concern. It is suggested that future studies investigate the quality of care provided in different regions of Iran to BC patients. The novel quality of care index (QCI) has been recently introduced in the literature, which seems to be a reliable tool for the quality-of-care assessment for a specific cause in each sex-age group-location-year [61,62,63].

Although economic sanctions against Iran in recent decades did not directly target healthcare systems and access to medication, it seems that they led to severe issues for the healthcare system and patients with cancer, exerting a severe economic burden on the whole system and limiting access to quality care. Future studies might investigate the possible burden associated with sanctions as a public health concern in Iran in recent years more in-depth.

Conclusions

In this study, we reported for the first time the most up-to-date estimates regarding BC burden in Iran for females and males at the national and subnational levels, based on age and SDI over thirty years from 1990 to 2019. We also performed decomposition analysis to find the attribution of population growth, population aging, and BC incidence rate in the total change of BC incident cases between 1990 and 2019 and reported MIRs—as an indicator of the quality of care—to provide a more comprehensive viewpoint on the current status of BC burden in Iran for the policymakers and health authorities. Besides, we assessed and reported the burden of BC attributable to known risk factors based on the GBD 2019 risk factor hierarchy. This study also will provide insights for other similar regions of the world to Iran to adopt the methodology of this study, which will help them reach comparable results and a broad vision concerning the burden of this lethal and prevalent cancer.

Still, more attention must be paid to Iran, at national and subnational levels, especially in low SDI regions. The policy-making should focus on systematic screening programs, especially for high-risk individuals, mass education about BC risk factors, self-examination, screening measures, healthy diet, lifestyle modifications, improving diagnostic and early detection tools, increasing access to the healthcare system and diagnostic facilities, developing comprehensive cancer registries, and cost-effective therapies to control the increasing trends. Developing an efficient NCCP could be a crucial step toward controlling the BC burden in Iran and other developing regions of the world. No doubt, further studies are needed to investigate the reasons behind the existing disparities in BC burden in different parts of Iran, explore more risk factors that might affect this burden, investigate the role of sex, hormones, genetic, and ethnic factors, assess the burden of BC in males more accurately, and investigate possible effective strategies to tackle the increasing burden of BC in the future more in-depth.

Availability of data and materials

Data regarding the burden of BC for Iran, at the national and subnational levels, were obtained from the Global Health Data Exchange (GHDx) query online tool, over 30 years from 1990 to 2019, which is publicly available at http://ghdx.healthdata.org/gbd-results-tool.

Change history

References

  1. Organization WH. Global Health Observatory. Geneva: World Health Organization; 2018. p. 2018.

    Google Scholar 

  2. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  PubMed  Google Scholar 

  3. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  4. Giordano SH. Breast cancer in men. N Engl J Med. 2018;378(24):2311–20.

    Article  PubMed  Google Scholar 

  5. Cui X. The prevalence and death risk of male breast cancer: a study based on the surveillance, epidemiology, and end results database. Am J Mens Health. 2022;16(1):15579883221074818.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Ataeinia B, Saeedi Moghaddam S, Shabani M, Gohari K, Sheidaei A, Rezaei N, et al. National and subnational incidence, mortality, and years of life lost due to breast cancer in iran: trends and age-period-cohort analysis since 1990. Front Oncol. 2021;11: 561376.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Roshandel G, Ferlay J, Ghanbari-Motlagh A, Partovipour E, Salavati F, Aryan K, et al. Cancer in Iran 2008 to 2025: recent incidence trends and short-term predictions of the future burden. Int J Cancer. 2021;149(3):594–605.

    Article  CAS  PubMed  Google Scholar 

  8. Jazayeri SB, Saadat S, Ramezani R, Kaviani A. Incidence of primary breast cancer in Iran: ten-year national cancer registry data report. Cancer Epidemiol. 2015;39(4):519–27.

    Article  PubMed  Google Scholar 

  9. Rezaianzadeh A, Jalali M, Maghsoudi A, Mokhtari AM, Azgomi SH, Dehghani SL. The overall 5-year survival rate of breast cancer among Iranian women: a systematic review and meta-analysis of published studies. Breast Dis. 2017;37(2):63–8.

    Article  PubMed  Google Scholar 

  10. Farzadfar F, Naghavi M, Sepanlou SG, Moghaddam SS, Dangel WJ, Weaver ND, et al. Health system performance in Iran: a systematic analysis for the Global Burden of Disease Study 2019. Lancet. 2022;399(10335):1625–45.

    Article  Google Scholar 

  11. The World by Income and Region. The World Bank. https://datatopics.worldbank.org/world-development-indicators/the-world-by-income-and-region.html. 2022.

  12. Li N, Deng Y, Zhou L, Tian T, Yang S, Wu Y, et al. Global burden of breast cancer and attributable risk factors in 195 countries and territories, from 1990 to 2017: results from the Global Burden of Disease Study 2017. J Hematol Oncol. 2019;12(1):140.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Enayatrad M, Amoori N, Salehiniya H. Epidemiology and trends in breast cancer mortality in Iran. Iran J Public Health. 2015;44(3):430.

    PubMed  PubMed Central  Google Scholar 

  14. Rafiemanesh H, Rajaei-Behbahani N, Khani Y, Hosseini S. Incidence trend and epidemiology of common cancers in the center of Iran. Global J Health Sci. 2016;8(3):146.

    Article  Google Scholar 

  15. Danaei G, Farzadfar F, Kelishadi R, Rashidian A, Rouhani OM, Ahmadnia S, et al. Iran in transition. Lancet. 2019;393(10184):1984–2005.

    Article  PubMed  Google Scholar 

  16. Ataeinia B, Saeedi Moghaddam S, Shabani M, Gohari K, Sheidaei A, Rezaei N, et al. National and subnational incidence, mortality, and years of life lost due to breast cancer in Iran: trends and age-period-cohort analysis since 1990. Front Oncol. 2021. https://doi.org/10.3389/fonc.2021.561376/full.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Vos T, Lim SS, Abbafati C, Abbas KM, Abbasi M, Abbasifard M, et al. Global burden of 369 diseases and injuries in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet. 2020;396(10258):1204–22.

    Article  Google Scholar 

  18. Stevens GA, Alkema L, Black RE, Boerma JT, Collins GS, Ezzati M, et al. Guidelines for accurate and transparent health estimates reporting: the GATHER statement. PLoS Med. 2016;13(6): e1002056.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Foreman KJ, Lozano R, Lopez AD, Murray CJ. Modeling causes of death: an integrated approach using CODEm. Popul Health Metrics. 2012;10(1):1–23.

    Article  Google Scholar 

  20. Sharma R, Abbasi-Kangevari M, Abd-Rabu R, Abidi H, Abu-Gharbieh E, Acuna JM, et al. Global, regional, and national burden of colorectal cancer and its risk factors, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet Gastroenterol Hepatol. 2022. https://doi.org/10.1016/S2468-1253(22)00044-9.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Roth GA, Abate D, Abate KH, Abay SM, Abbafati C, Abbasi N, et al. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018;392(10159):1736–88.

    Article  Google Scholar 

  22. Wang H, Abbas KM, Abbasifard M, Abbasi-Kangevari M, Abbastabar H, Abd-Allah F, et al. Global age-sex-specific fertility, mortality, healthy life expectancy (HALE), and population estimates in 204 countries and territories, 1950–2019: a comprehensive demographic analysis for the Global Burden of Disease Study 2019. Lancet. 2020;396(10258):1160–203.

    Article  Google Scholar 

  23. Mokdad AH, Mensah GA, Krish V, Glenn SD, Miller-Petrie MK, Lopez AD, et al. Global, regional, national, and subnational big data to inform health equity research: perspectives from the Global Burden of Disease Study 2017. Ethn Dis. 2019;29(Suppl 1):159.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Fitzmaurice C, Allen C, Barber RM, Barregard L, Bhutta ZA, Brenner H, et al. Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 32 cancer groups, 1990 to 2015: a systematic analysis for the global burden of disease study. JAMA Oncol. 2017;3(4):524–48.

    Article  PubMed  Google Scholar 

  25. Sharma R. Breast cancer incidence, mortality and mortality-to-incidence ratio (MIR) are associated with human development, 1990–2016: evidence from Global Burden of Disease Study 2016. Breast Cancer. 2019;26(4):428–45.

    Article  PubMed  Google Scholar 

  26. Azadnajafabad S, Saeedi Moghaddam S, Keykhaei M, Shobeiri P, Rezaei N, Ghasemi E, et al. Expansion of the quality of care index on breast cancer and its risk factors using the Global Burden of Disease Study 2019. Cancer Med. 2022. https://doi.org/10.1002/cam4.4951.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Murray CJ, Aravkin AY, Zheng P, Abbafati C, Abbas KM, Abbasi-Kangevari M, et al. Global burden of 87 risk factors in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. The Lancet. 2020;396(10258):1223–49.

    Article  Google Scholar 

  28. Boyle P, Boniol M, Koechlin A, Robertson C, Valentini F, Coppens K, et al. Diabetes and breast cancer risk: a meta-analysis. Br J Cancer. 2012;107(9):1608–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Fund WCR, Research AIfC. Food, nutrition, physical activity, and the prevention of cancer: a global perspective. Am Inst Cancer Res. 2007.

  30. James SL, Abate D, Abate KH, Abay SM, Abbafati C, Abbasi N, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. The Lancet. 2018;392(10159):1789–858.

    Article  Google Scholar 

  31. Griswold MG, Fullman N, Hawley C, Arian N, Zimsen SR, Tymeson HD, et al. Alcohol use and burden for 195 countries and territories, 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016. The Lancet. 2018;392(10152):1015–35.

    Article  Google Scholar 

  32. Radecka B, Litwiniuk M. Breast cancer in young women. Ginekol Pol. 2016;87(9):659–63.

    Article  PubMed  Google Scholar 

  33. Chen Z, Xu L, Shi W, Zeng F, Zhuo R, Hao X, et al. Trends of female and male breast cancer incidence at the global, regional, and national levels, 1990–2017. Breast Cancer Res Treat. 2020;180(2):481–90.

    Article  PubMed  Google Scholar 

  34. Giordano SH, Cohen DS, Buzdar AU, Perkins G, Hortobagyi GN. Breast carcinoma in men: a population-based study. Cancer: Interdisciplinary. Int J Am Cancer Soc. 2004;101(1):51–7.

    Google Scholar 

  35. Nafissi N, Khayamzadeh M, Zeinali Z, Pazooki D, Hosseini M, Akbari ME. Epidemiology and histopathology of breast cancer in Iran versus other Middle Eastern countries. Middle East J Cancer. 2018;9(3):243–51.

    Google Scholar 

  36. Gao Y, Goldberg JE, Young TK, Babb JS, Moy L, 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(2):282–91.

    Article  PubMed  Google Scholar 

  37. Asadzadeh Vostakolaei F, Karim-Kos HE, Janssen-Heijnen ML, Visser O, Verbeek AL, Kiemeney LA. The validity of the mortality to incidence ratio as a proxy for site-specific cancer survival. European J Public Health. 2011;21(5):573–7.

    Article  Google Scholar 

  38. Chen W-J, Huang C-Y, Huang Y-H, Wang S-C, Hsieh T-Y, Chen S-L, et al. Correlations between mortality-to-incidence ratios and health care disparities in testicular cancer. Int J Environ Res Public Health. 2020;17(1):130.

    Article  Google Scholar 

  39. Li N, Deng Y, Zhou L, Tian T, Yang S, Wu Y, et al. Global burden of breast cancer and attributable risk factors in 195 countries and territories, from 1990 to 2017: results from the Global Burden of Disease Study 2017. J Hematol Oncol. 2019;12(1):1–12.

    Article  Google Scholar 

  40. Waks AG, Winer EP. Breast cancer treatment: a review. JAMA. 2019;321(3):288–300.

    Article  CAS  PubMed  Google Scholar 

  41. Romaniuk A, Lyndin M, Sikora V, Lyndina Y, Romaniuk S, Sikora K. Heavy metals effect on breast cancer progression. J Occupat Med Toxicol. 2017;12(1):1–9.

    Google Scholar 

  42. Rahimzadeh S, Burczynska B, Ahmadvand A, Sheidaei A, Khademioureh S, Pazhuheian F, et al. Geographical and socioeconomic inequalities in female breast cancer incidence and mortality in Iran: A Bayesian spatial analysis of registry data. PLoS ONE. 2021;16(3): e0248723.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Hiller TW, O’Sullivan DE, Brenner DR, Peters CE, King WD. Solar Ultraviolet radiation and breast cancer risk: a systematic review and meta-analysis. Environ Health Perspect. 2020;128(1): 016002.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Esmaeimzadeh N, Salahi-Moghaddam A, Khoshdel A. Geographic distribution of important cancers in Iran. Hormozgan Med J. 2015;19(2):66–76.

    Google Scholar 

  45. Keyghobadi N, Rafiemanesh H, Mohammadian-Hafshejani A, Enayatrad M, Salehiniya H. Epidemiology and trend of cancers in the province of Kerman: southeast of Iran. Asian Pac J Cancer Prev. 2015;16(4):1409–13.

    Article  PubMed  Google Scholar 

  46. Semnani S, Sadjadi A, Fahimi S, Nouraie M, Naeimi M, Kabir J, et al. Declining incidence of esophageal cancer in the Turkmen Plain, eastern part of the Caspian Littoral of Iran: a retrospective cancer surveillance. Cancer Detect Prev. 2006;30(1):14–9.

    Article  PubMed  Google Scholar 

  47. Shield KD, Soerjomataram I, Rehm J. Alcohol use and breast cancer: a critical review. Alcoholism Clin Exper Res. 2016;40(6):1166–81.

    Article  CAS  Google Scholar 

  48. Rumgay H, Shield K, Charvat H, Ferrari P, Sornpaisarn B, Obot I, et al. Global burden of cancer in 2020 attributable to alcohol consumption: a population-based study. Lancet Oncol. 2021;22(8):1071–80.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Chen Z, Shao J, Gao X, Li X. Effect of passive smoking on female breast cancer in China: a meta-analysis. Asia Pacific J Public Health. 2015;27(2):NP58–64.

    Article  Google Scholar 

  50. Shamshirian A, Heydari K, Shams Z, Aref AR, Shamshirian D, Tamtaji OR, et al. Breast cancer risk factors in Iran: a systematic review & meta-analysis. Hormone Mole Biol Clin Investig. 2020. https://doi.org/10.1515/hmbci-2020-0021/html.

    Article  Google Scholar 

  51. Argolo DF, Hudis CA, Iyengar NM. The impact of obesity on breast cancer. Curr Oncol Rep. 2018;20(6):47.

    Article  PubMed  Google Scholar 

  52. Engin A. Obesity-associated breast cancer: analysis of risk factors. Obesity and lipotoxicity. Springer; 2017. p. 571–606.

  53. Kazemi A, Barati-Boldaji R, Soltani S, Mohammadipoor N, Esmaeilinezhad Z, Clark CC, et al. Intake of various food groups and risk of breast cancer: A systematic review and dose-response meta-analysis of prospective studies. Adv Nutr. 2021;12(3):809–49.

    Article  PubMed  Google Scholar 

  54. Michels KB, Mohllajee AP, Roset-Bahmanyar E, Beehler GP, Moysich KB. Diet and breast cancer: a review of the prospective observational studies. Cancer Interdiscip Int J Am Cancer Soc. 2007;109:2712–49.

    CAS  Google Scholar 

  55. Tamimi RM, Spiegelman D, Smith-Warner SA, Wang M, Pazaris M, Willett WC, et al. Population attributable risk of modifiable and nonmodifiable breast cancer risk factors in postmenopausal breast cancer. Am J Epidemiol. 2016;184(12):884–93.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Fitzmaurice C, Abate D, Abbasi N, Abbastabar H, Abd-Allah F, Abdel-Rahman O, et al. Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2017: a systematic analysis for the global burden of disease study. JAMA Oncol. 2019;5(12):1749–68.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Mohammadi E, Aminorroaya A, Fattahi N, Azadnajafabad S, Rezaei N, Farzi Y, et al. Epidemiologic pattern of cancers in Iran; current knowledge and future perspective. J Diabetes Metabol Disorders. 2020. https://doi.org/10.1007/s40200-020-00654-6.

    Article  Google Scholar 

  58. Flaxman AD, Vahdatpour A, James SL, Birnbaum JK, Murray CJ. Direct estimation of cause-specific mortality fractions from verbal autopsies: multisite validation study using clinical diagnostic gold standards. Popul Health Metrics. 2011;9(1):1–10.

    Google Scholar 

  59. Mikkelsen L, Phillips DE, AbouZahr C, Setel PW, De Savigny D, Lozano R, et al. A global assessment of civil registration and vital statistics systems: monitoring data quality and progress. Lancet. 2015;386(10001):1395–406.

    Article  PubMed  Google Scholar 

  60. Mohammadi G, Akbari ME, Mehrabi Y, Ghanbari Motlagh A. Quality assessment of the national cancer registry in Iran: completeness and validity. Iran J Cancer Prevent. 2016;9(6).

  61. Aryannejad A, Tabary M, Ebrahimi N, Mohammadi E, Fattahi N, Roshani S, et al. Global, regional, and national survey on the burden and quality of care of pancreatic cancer: a systematic analysis for the global burden of disease study 1990–2017. Pancreatology. 2021;21(8):1443–50.

    Article  PubMed  Google Scholar 

  62. Azadnajafabad S, Saeedi Moghaddam S, Mohammadi E, Rezaei N, Ghasemi E, Fattahi N, et al. Global, regional, and national burden and quality of care index (QCI) of thyroid cancer: a systematic analysis of the Global Burden of Disease Study 1990–2017. Cancer Med. 2021;10(7):2496–508.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Mohammadi E, Ghasemi E, Azadnajafabad S, Rezaei N, Saeedi Moghaddam S, Ebrahimi Meimand S, et al. A global, regional, and national survey on burden and Quality of Care Index (QCI) of brain and other central nervous system cancers; global burden of disease systematic analysis 1990–2017. PLoS ONE. 2021;16(2): e0247120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the Institute for Health Metrics and Evaluation (IHME) for providing these invaluable and high-quality epidemiological estimations of the burden associated with diseases and injuries globally, almost in all regions and countries of the world, helping to reduce and eliminate poverty of knowledge and information worldwide. Also, the authors sincerely thank all the GBD collaborators who contributed to this study by providing data or valuable and critical feedback through precise reviewing of the manuscript.

Funding

The authors declare that they received no specific funding for this work. The Global Burden of Diseases (GBDs) study received funding from the Bill and Melinda Gates Foundation.

Author information

Authors and Affiliations

Authors

Consortia

Contributions

Please see Additional file 15 for more detailed information about individual author contributions to the research, divided into the following categories: providing data or critical feedback on data sources; developing methods or computational machinery; providing critical feedback on methods or results; drafting the work or revising is critically for important intellectual content; and managing the estimation or publications process.

Corresponding author

Correspondence to Farshad Farzadfar.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the ethical committee of EMRI under the reference code of IR.TUMS.EMRI.REC.1400.031.

Consent for publication

Not Applicable.

Competing interests

The authors declare no conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

This article has been updated to add the contributions of the GBD Iran Breast Cancer Collaborators.

Supplementary Information

Additional file 1

. Fig. 1 The geographical location of Iran in the globe (shown in red color) and its subnational administrative divisions.

Additional file 2

. Table 1 Scientific papers and data registry systems that were used by the GBD study 2019 to model and estimate the burden of breast cancer in Iran at national and subnational levels.

Additional file 3

. Table 2 Annual percent changes of breast cancer age-standardized incidence, prevalence, deaths, disability-adjusted life years (DALYs), years of life lost (YLLs) and years lived with disability (YLDs) rates from 1990 to 2019 at the national and subnational levels.

Additional file 4

. Fig. 2A Breast cancer age-standardized incidence and deaths rates (per 100,000 population) rankings in Iran and its 31 provinces in 1990 and 2019, for A) both sexes, B) females, and C) males.

Additional file 5

. Fig. 2B Breast cancer age-standardized incidence and deaths rates (per 100,000 population) rankings in Iran and its 31 provinces in 1990 and 2019, for A) both sexes, B) females, and C) males.

Additional file 6

. Fig. 2C Breast cancer age-standardized incidence and deaths rates (per 100,000 population) rankings in Iran and its 31 provinces in 1990 and 2019, for A) both sexes, B) females, and C) males.

Additional file 7

. Table 3 Breast cancer age-standardized incidence, prevalence, deaths, disability-adjusted life years (DALYs), years of life lost (YLLs) and years lived with disability (YLDs) rates (per 100,000 population), numbers and percent changes in all 31 provinces of Iran between 1990 and 2019, for both sexes, females, and males (provinces have been sorted alphabetically).

Additional file 8

. Table 4 Percent changes (%) in breast cancer age-standardized rates and all-ages numbers of incidence, prevalence, deaths, disability-adjusted life years (DALYs), years of life lost (YLLs), and years lived with disability (YLDs) over two major time periods, between 1990 and 2010 and between 2010 and 2019, in Iran and all 31 provinces for both sexes, females, and males.

Additional file 9

. Fig. 3 Breast cancer incidence rates (per 100,000 population) in 1990 and 2019 based on age groups by sex (red: female; blue: males) in Iran and its 31 provinces.

Additional file 10

. Fig. 4 Breast cancer prevalence rates (per 100,000 population) in 1990 and 2019 based on age groups by sex (red: female; blue: males) in Iran and its 31 provinces.

Additional file 11

. Fig. 5 Breast cancer deaths rates (per 100,000 population) in 1990 and 2019 based on age groups by sex (red: female; blue: males) in Iran and its 31 provinces.

Additional file 12

. Fig. 6 Breast cancer disability-adjusted life years (DALYs) rates (per 100,000 population) in 1990 and 2019 based on age groups by sex (red: female; blue: males) in Iran and its 31 provinces.

Additional file 13

. Table 5 Pearson’s correlation coefficient (r) between measures of burden and socio-demographic index (SDI)

Additional file 14

. Table 6 Breast cancer deaths and disability-adjusted life years (DALYs) rate (per 100,000 population), numbers and percent changes attributable to risk factors in all 31 provinces of Iran between 1990 and 2019 for both sexes, females, and males (provinces have been sorted alphabetically).

Additional file 15

. Contributions by the GBD Iran Breast Cancer Collaborators.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Aryannejad, A., Saeedi Moghaddam, S., Mashinchi, B. et al. National and subnational burden of female and male breast cancer and risk factors in Iran from 1990 to 2019: results from the Global Burden of Disease study 2019. Breast Cancer Res 25, 47 (2023). https://doi.org/10.1186/s13058-023-01633-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13058-023-01633-4

Keywords