Skip to main content

Breast cancer-specific survival by molecular subtype in different age groups of women in Scotland

Abstract

Background

Age and molecular subtypes are important prognostic factors in breast cancer (BC). Here, we explore how age and molecular subtypes influence BC survival in Scotland.

Methods

We analysed data from 71,784 women diagnosed with invasive BC in Scotland between 1997 and 2016, with follow-up until 31st December 2018 (median follow-up time = 5.5 years). Cox models estimated Hazard Ratios (HR) for BC-specific death by age group (with women of screening age, 50–69 years old, as the reference) within each molecular subtype, adjusting for prognostic factors. The cumulative incidence function was plotted to account for competing risks.

Results

During the study period, 37% of women died, with 53% of deaths attributed to BC. Women aged 70 + years had increased BC-specific death compared to women aged 50 to 69 years with the same subtype. HRs (95% CI) were 1.49 (1.23–1.80) for luminal A, 1.39 (1.14 to 1.69) for luminal B tumours and 1.49 (1.15 to 1.94) for triple negative breast cancer (TNBC). Women aged < 50 years had lower risk of BC death in luminal A subtype only, with HR of 0.66 (0.51–0.86) compared to women aged 50 to 69 years. Competing risks analysis showed higher cumulative incidence of death from non-BC causes, particularly for women aged 70 + years with hormone positive subtypes. Stage, treatment, and molecular subtype were the strongest prognostic factors for BC-specific mortality across all ages.

Conclusions

Age influences BC-specific mortality particularly within luminal subtypes. In contrast, other tumour characteristics and treatment are key prognostic factors for non-luminal subtypes. Future studies should investigate other markers of BC mortality particularly among over 70-year-olds, who account for 60% of BC deaths in the UK.

Background

Breast cancer (BC) is the leading cause of cancer mortality among women in over 100 countries and is currently the most common cancer globally, having overtaken lung cancer [1]. Age is an important risk factor for BC, with nearly half of cases in the United Kingdom (UK) diagnosed in women aged 50 to 70 years and a third of cases in women over 70 years of age [2]. Further, 90% of BC deaths in the UK occur in women over 50 years at the time of death, with 60% of these deaths in women over 70 years of age [2]. Given the increasing life expectancy and projections that by the middle of the 21st century almost a third of the UK population will be over 70 years of age [3], this age group is an increasingly important demographic for public health.

Despite the role of age in BC incidence and prognosis, chronological age alone should not be used as a reason to withhold specific treatments [4] since the benefit from chemotherapy is independent of age [5]. Further, age is not the only important prognostic factor in BC, multiple previous studies have also highlighted the role of the BC molecular subtypes [6], with luminal subtypes having more favourable prognosis than triple negative BC (TNBC) or human epidermal growth factor receptor (HER2) enriched subtypes [7,8,9,10,11,12]. However, fewer studies have investigated how age influences BC prognosis for the different molecular subtypes, and those who had lacked adjustment for important confounding factors, such as socio-economic status, comorbidities and competing risks.

The Scottish cancer registry is an excellent resource to address this gap, as it includes long-term data on oestrogen receptor (ER) status (since 1997) and progesterone receptor (PR) and HER2 status (since 2009)- almost a decade earlier than other UK national registries [6]. Further, the cancer registry data can be linked to other health electronic records to enable adjustment for important confounders. We have previously shown that age is an important predictor of incidence of BC tumour subtype in Scotland [6]. Here, we aimed to investigate the role of age on BC survival across the different molecular subtypes. Specifically, we assessed (1) how survival outcomes vary by molecular subtypes across different age groups, (2) whether age modifies the prognostic impact of the molecular subtypes and (3) the association of other tumour characteristics, screening, treatment, socio-economic status, comorbidities and competing risks on survival.

Methods

Data and cohort definition

In this retrospective cohort study, women aged 20 years or older diagnosed with a first primary invasive BC (defined using the International Classification of Diseases, 10th revision code [ICD10] of C50) between 1997 and 2016 in Scotland were identified from the Scottish cancer registry data [13]. The Scottish Cancer registry was established in 1958, and electronic data was available from 1981. The registry covers all Scottish residents registered with a general practitioner (GP) reason why the coverage of the registry is really high, estimated at 98% for BC cases independently of age [14]. Women with missing postcodes, living outside Scotland, or who were older than 99 years (n = 55), with unknown vital status (n = 154), with BC diagnosed by death certificate only (n = 125) or who had the same date of incidence as date of death (n = 99) were excluded (0.6% of the total sample). The final study population consisted of 71,784 women.

For each woman, information on year of BC diagnosis (in 5- year calendar groups), age at BC diagnosis, Scottish region defined by Scottish geographical regions (North, South East and West), tumour characteristics (stage, tumour grade, tumour size, nodal status, method of detection, ER, PR and HER2 status) and binary categories of treatment (surgery, radiotherapy, chemotherapy and hormone therapy) was ascertained by the cancer registry from medical and pathology records. Information on targeted treatment for HER2 tumours was not available.

An area-based measure of socio-economic status at diagnosis, the Scottish Index of Multiple Deprivation (SIMD), was derived for all women within the cancer registry with a Scottish postcode. The SIMD deprivation measure is based on 7 domains: income, employment, health, education, crime, access to services and housing [15]. SIMD is often expressed in quintiles and we used SIMD quintile to compare women in the 20% most deprived areas (quintile 1) with women in the 20% least deprived areas (quintile 5). The Charlson index of comorbidity, a measure that consists on a weighted index that takes into account the number of comorbidities and the severity of each comorbidity [16], was linked to the registry data from hospital inpatient records.

ER and PR status are measured by immunohistochemistry (IHC). HER2 status was assessed using a combination of IHC with fluorescent in-situ hybridisation for equivocal (2+) cases and the three markers were defined as positive, negative or unknown. BC subtypes were derived from the combination of ER, PR and HER2 based on the definition adopted in the St. Gallen 2011 consensus [12] and served as surrogates for the four intrinsic molecular subtypes of BC. As the registry lacks information for Ki67, tumour grade was used as a marker for cell proliferation. ER + and/or PR + and HER2- (grade 1 and 2) tumours were defined as luminal A, ER + and/or PR + and HER2- (grade 3) and ER + and/or PR + and HER2 + as luminal B, ER- and PR- and HER2 + as HER2-enriched and ER- and PR- and HER2- as TNBC.

Age at diagnosis was stratified in three categories based on the age at which women are routinely invited for mammography screening in Scotland: women aged less than 50 years, women aged 50 to 69 years (current age range for routine screening) and women aged 70 years or older (who only receive screening on request). From 1997 to 2004 the age range for routine breast screening was 50–64 years. Tumour grade was defined as low grade or well differentiated (grade I), medium grade or moderately well differentiated (grade II) and high grade or poorly differentiated (grade III). Information on method of detection was used to categorise each diagnosed tumour as screen-detected or not screen-detected (which included all remaining options: clinical examination, incidental finding and unknown).

Stage was derived from individual tumour (T), nodes (N) and metastasis (M) clinical and pathological stages. Clinical TNM stage was available during the entire study period but pathological TNM stage was only available from 2005. As a general principle, pathological TNM stage was prioritised over clinical stage as it tends to be more accurate. There were some exceptions: (1) for M stage except when clinical M stage was unknown or pathological M stage indicated metastasis (stage IV) and clinical stage did not, (2) if the woman had neoadjuvant therapy [radiotherapy, chemotherapy or hormone therapy (HT)] at least 4 weeks before surgery and (3) if clinical T stage was T4 indicating primary tumour involvement of chest wall or skin which is often obvious at clinical examination. Following the rules above and using pathological tumour size and the number of positive nodes to complete missing pathological T and N variables a final TNM stage variable was derived and categorised as I, II, III and IV following the American Joint Committee on Cancer (AJCC) 8th Edition Cancer Staging manual [17].

Statistical methods

Breast cancer specific survival (BCSS) was the primary outcome of the analysis. Primary and secondary causes of death are derived from death records linked to the cancer registry [18]. ICD9 174 and ICD10 C50 codes from primary cause of death were used to derive BC specific death using the approach described by Skyrud et al. [18]. Other primary causes of death were regarded as censored observations for the calculation of BCSS. Duration of follow-up was defined as time from date of diagnosis of BC to the earliest date of death, 31st December 2018 for women still alive at the end of the study period or embarkation date if women moved from Scotland.

Non- parametric Kaplan-Meier (KM) estimates [19] were used to estimate BCSS at 5 years by ER status and for each IHC defined molecular subtype for all women and by age group, addressing the question of how age and tumour subtype influences BCSS. Cox proportional hazards models [20] were fitted to investigate the association between the molecular subtypes and age category (as main predictors) and BC-specific mortality over time after controlling for other covariates and to assess which factors were more important for prognosis. Preliminary models included the two main exposures (ER status or the IHC defined molecular subtypes and age category at diagnosis), year of diagnosis and Scottish NHS region (model 1). Model 2 included model 1 plus tumour characteristics. Model 3 included model 2 plus treatment regimens and model 4 included model 3 plus deprivation and comorbidity measures. We also tested for an interaction between ER status and age category and between IHC defined subtypes and age category by comparing models with and without the interaction. The models were restricted to those women with information for all covariates (n = 51,140, 71% of the total). Sensitivity analyses was also conducted to estimate interactions between the model covariates with time and are presented (for the ER models only). Time varying covariates were used when the PH assumption from visual Schoenfeld residuals did not hold (this assumption is often violated in the presence of markers for BC) [21, 22].

Given that a significant interaction between molecular subtypes (either defined using ER status or the St. Gallen’s criteria) and age was found and that the PH assumption did not hold for the molecular subtypes, the main models presented were stratified by molecular subtypes to assess the effect of age (unadjusted and adjusted for other covariates) on BC-specific death within each molecular subtype. Hazard Ratio (HR) with 95% CI are presented for all models.

Competing risks analysis was performed to investigate whether the risk of early death from other causes differed between the molecular subtypes and age groups by plotting the cumulative incidence function [23] stratified by age and the molecular subtypes, with deaths from causes other than BC as the competing risk. The cumulative incidence function partitioned the contributions to mortality from BC-specific death and other causes of death and was estimated using Fine and Gray subdistribution hazard function [24]. All tests were 2-sided and deemed significant at the 5% level. Analyses were performed in Stata version 15/IC [25].

Results

Characteristics of cohort

Of the 71,784 women diagnosed with BC between 1997 and 2016 in Scotland with a median follow-up of 5.5 years, 26,280 (37%) died during the study period (Table 1). Of those who died, 53% had BC as their underlying cause of death. Differences in survival were observed by age group. Over 75% of women aged less than 70 years were alive at the end of the follow-up whereas only 35% of women 70 years and older survived. Among women younger than 50 years of age at BC diagnosis who died, 86% died of BC. Proportions of other causes of death were 43% and 57% respectively among women aged 50 to 69 years and older than 70 years at BC diagnosis.

Table 1 Individual characteristics, tumour characteristics, treatment regimens and mortality figures for the population of women diagnosed with BC in Scotland (1997 to 2016) by age group

Over 75% of tumours in the study population were ER+ (N = 55,003) with half of these tumours diagnosed in women of screening age 50 to 69 years (N = 28,385). Younger women (< 50 years) and older women (> 70 years) were slightly less likely to be diagnosed with an ER + tumour but the vast majority of tumours were hormone positive in all age groups (73% and 74% respectively compared to 80% for 50–69-year-old women) (Table 1). The distribution of molecular subtypes differed between age groups with women younger than 50 years having higher proportions of luminal B and TNBC tumours (36% and 13% respectively) and a smaller proportion of luminal A tumours (39%) than women aged 50 years or older. The proportion of HER2-enriched tumours was similar across age groups, accounting for approximately 4% of all tumours.

Other tumour characteristics also differed by age group. Women of screening age (50 to 69 years) had almost half of their tumours diagnosed through screening and had characteristics associated with less aggressive disease, including lower tumour stage, grade and size and lower positive nodal status when compared to women younger than 50 years. Women aged 70 years or older were more likely to have missing tumour characteristics than women younger than 70 years.

Treatment differences were also observed between age groups (Table 1). Women aged 70 years or older were less likely to have surgery, radiation and chemotherapy and more likely to have hormone therapy and neoadjuvant hormone therapy than women younger than 70 years.

Breast cancer-specific survival by molecular subtypes across different age groups

BCSS was highest amongst women aged 50 to 69 years, particularly for those diagnosed with ER + tumours (92.1% survival at 5 years), and lowest amongst women aged 70 years or older diagnosed with ER- tumours (55.4% survival at 5 years) (Table 2). Women younger than 50 years had only slightly lower BCSS compared to those aged 50–69, while women aged 70 years or older had lower survival across all tumour subtypes, particularly for HER2-enriched and TNBC tumours. Among all IHC-defined subtypes, TNBC had the lowest survival (72.6% survival at 5 years), regardless of age.

Table 2 Breast cancer specific survival estimates (in %) at 5 years after diagnosis (with 95% CI) by ER status and IHC defined subtypes for women of all ages (total) diagnosed in Scotland and by age group

Data from the preliminary fully adjusted models (Supplementary Tables 1 and 2) showed that age and molecular subtypes were independent significant prognostic factors for BC-specific death and supported an interaction between age and molecular subtype. Table 3 presents the effect on BCSS of each age group (compared to the reference 50–69 years old group) stratified by molecular subtype. After adjusting for covariates, women younger than 50 years were less likely to die from BC-specific death compared to women aged 50 to 69 years, particularly for ER + tumours [HR = 0.89 (95%CI: 0.82 to 0.96)], luminal A tumours [HR = 0.66 (95%CI: 0.51 to 0.86)] and HER2-enriched tumours, although the latter association with age was not statistically significant.

Table 3 Unadjusted and adjusted Cox models for breast cancer specific mortality stratified by molecular subtype showing the effect of each age group (compared to the reference group of women aged 50 to 69 years) in each molecular subtype

For women aged 70 years or older, the risk of BC-specific death was higher for all tumour subtypes compared to the 50 to 69 age group, for both ER + and ER- tumours (34% and 30% increased risk, respectively). This elevated risk remained after adjusting for all covariates, except for HER-2 enriched tumours (Table 3). The highest relative risk was observed among women aged 70 + with luminal A and TNBC tumours (HR = 1.49 for both).

Tumour characteristics, screening, treatment, socio-economic status and comorbidities in BCSS

Women with ER- tumours had worse survival outcomes compared to women with ER + tumours, while luminal B, HER2-enriched and TNBC subtypes were associated with significantly higher risks of BC death (adjusted HRs of 2.04, 1.95 and 3.93 respectively, compared to luminal A tumours) (Supplementary Tables 1 and 2).

Apart from molecular subtype, other tumour characteristics -specifically, TNM stage and grade- as well as method of detection were the strongest prognostic factors. Adjusting for these factors (model 2) attenuated the HR for age and molecular subtype. Treatments (model 3) were also significant prognostic factors for BC-specific death, and their adjustment further reduced the HR for age and molecular subtype. However, adjustment for socio-economic status and comorbidities (model 4) had limited additional effect on the associations between BC-specific death with age and molecular subtype. A directed acyclic graph (DAG) of the effect of the main covariates in BCSS is presented in Supplementary Fig. 1.

Time-varying effects of prognostic factors

To assess whether key prognostic factors changed over time, we present the results for the Cox model analysis incorporating time-varying covariates (Supplementary Table 3) with ER status as the main predictor. We observed a statistically significant interaction between time and ER status [HR for time*ER = 0.83 (95% CI: 0.81 to 0.85)], indicating that the HR for the comparison of ER- and ER + tumours decreased over time from diagnosis. HR estimates at different time points after diagnosis (1, 3, 5 and 10 years) are presented in Supplementary Table 4, which may help inform clinical decisions for improving survival.

Competing risks of BC death amongst the different age groups and by subtype

Apart from BC, cardiovascular diseases (CVD), other cancers, chronic obstructive pulmonary disease (COPD) and Alzheimer’s disease or dementia were the most common primary causes of death in Scottish women diagnosed with BC (Supplementary Table 5). There was a clear relationship between age and distribution of cause of death, with a higher proportion of women aged 70 years or older dying from causes other than BC than among women aged less than 70 years. For example, CVD accounted for 22% of all deaths in women aged 70 years or older compared to 10% of total deaths in women aged 50 to 69 years and 2% in women aged less than 50 years at BC diagnosis (Supplementary Table 5).

The cumulative incidence function showed that the probability of dying from BC was highest in women aged 70 years or older with ER- tumours (Fig. 1b), particularly within the first 5 years after BC diagnosis. The probability of dying from BC for women with an ER- tumour aged less than 50 years was also slightly higher than for those aged 50 to 69 years with the same subtype (Fig. 1a). Competing risks of death (other than BC) were particularly high for the 70 years age group diagnosed with an ER + tumour for which the probability was similar to that of BC death in the first 5 years of follow-up and considerably higher after 5 years of follow-up. Younger women under 50 years had the lowest probability of dying from causes other than BC irrespective of ER status. Cumulative incidence functions by molecular subtypes show the highest probability of BC death amongst women aged 70 years or older with a TNBC or HER2-enriched subtype (Fig. 1e and f). Competing causes of death were more likely than BC death only for women aged 70 years or older diagnosed with a luminal A tumour (Fig. 1c), who contributed about 30% of all luminal A tumours. Women younger than 50 years had a slightly higher risk of BC death if diagnosed with a luminal tumour than women aged 50 to 69 years but were less likely to die from other causes.

Fig. 1
figure 1

Cumulative incidence graph of BC death (breast: red line) and other cause of death (other: blue line) by age group for (A) ER + tumours (n = 55,003), (B) ER- tumours (n = 11,698), (C) luminal A tumours (n = 15,424), (D) luminal B tumours (n = 7,928), (E) HER2-enriched tumours (n = 1,285) and (F) TNBC tumours (n = 2,895)

Discussion

In this large population-based study, including over 70,000 women diagnosed with BC in Scotland between 1997 and 2016- more than 51,000 of whom had complete data- we observed that BCSS survival and mortality differed by molecular subtype and age. Women with ER- tumours had nearly 20% lower BCSS than women with ER + tumours and a 44% higher risk of BC mortality, even after adjusting for other individual and tumour characteristics, treatments and comorbidities.

Our analysis of IHC-defined molecular subtypes showed that women with luminal A tumours had the highest five-year survival (91%), while women with TNBC subtypes had the worst BCSS (73% at five years). Additionally, we observed a 2-fold increased risk of BC-specific death for luminal B and HER2-enriched tumours and 4-fold increase for TNBC tumours when compared to luminal A tumours. These results are similar to those observed in the US [7, 26], Canada [27], Italy [28], Germany [8], Spain [29], Switzerland [30] and Norway [9]. However, direct international comparisons are challenging due to the differences in health care systems, breast screening programmes, access to targeted therapies and differences in the incidence of non-luminal subtypes.

BCSS by molecular subtypes across different age groups

We observed differences in BCSS by age, with women aged 70 years or older having consistently higher risk of BC death compared to women of screening age (50 to 69 years), whereas women younger than 50 years had similar survival outcomes to those in the screening age group. Further, our data found evidence of survival outcomes differing by molecular subtype and age. A recent study from Norway of 21,384 women diagnosed from 2005 to 2015, showed younger age (< 40 years) was associated with higher breast cancer specific mortality among women with luminal A-like BC subtype, while old age was associated with increased death across all subtypes [9]. Notably, the 40 -to-49-year age group showed different associations with luminal A status, with a less aggressive disease characteristics than women < 40 years old (who had higher proportion of high grade HER2 + tumours). Our findings were consistent in showing women older than 70 years having worse BCSS for all tumour subtypes (except for HER2-enriched) compared to those 50–69 years. However, the magnitude of the age-related differences was smaller in Scotland than in Norway [9]. In contrast, our results showed that women younger than 50 years had lower risk of BC death compared to women of screening age among those with ER + tumours and luminal A tumours after adjusting for other covariates. Due to the limitations of our data and the small sample size for age groups < 40 and 40–49 years, we were unable to further stratify by these specific age groups. Future research should provide more granular association results by subtypes in women younger than 40 years old.

Associations of tumour characteristics, screening, treatment, socio-economic status and comorbidities in BCSS

Tumour subtype, tumour characteristics (particularly stage and grade), screening and treatment attenuated the association of age with BCSS, showing that these make an important contribution to the association. In contrast, socio-economic status or comorbidities had little influence in BCSS. Women aged 70 years or older and women younger than 50 years had tumour characteristics of a more aggressive disease (higher proportions of late stage and higher-grade tumours) compared to women aged 50–69 years. Further, women aged 70 years or older were less likely to undergo surgery compared to women of 50–69 years of age. In general, the association of age category with BCSS was less strong than for subtype and tumour characteristics. Our data provide further evidence that tumour characteristics, screening and treatment are key factors in predicting BCSS.

Our analysis shows lower proportions of women aged 70 years and older treated with surgery, chemotherapy and radiotherapy than for women younger than 70 years. Although presence of comorbidities or frailty may influence treatment decisions, recent studies in the UK suggest that older women could benefit from more extensive use of chemotherapy and radiotherapy after primary surgical treatment in terms of survival and recurrence after diagnosis of BC [31]. However, almost half of women aged 70 + years diagnosed with BC who died did so from causes other than BC with CVD being a major contributor. Previous studies have indicated that some BC treatments increase CVD risk, such as radiotherapy to the left side of the chest and specific chemotherapy agents [32, 33]. Further research is needed to investigate risk of CVD among women who have received newer treatments to help identify the balance of risk and benefits of treatment for different populations of women with BC.

Older women are underrepresented in clinical trials [34]. Age is a key risk factor for BC incidence and there is an increasing number of studies investigating age-related biologic measures [35]. However, it is not clear which measures are the most important and whether their influence differs by tumour subtypes. Some studies have suggested frailty is an important factor in decisions about the appropriateness of chemotherapy [36] as a proxy for biological ageing [37] that could help identify those most likely to benefit from treatments. Interestingly, treatment with radiation or chemotherapy may accelerate ageing particularly among older patients [38, 39]. In our data, we observed missing stage and other tumor characteristics are likely to be higher in woman older than 70 which is likely due to these women being less likely to have surgery. Wyld et al. reported that ~ 20% of women older than 70 will have primary endocrine therapy in the UK [40]. This is likely due to clinicians having to manage frail or comorbid conditions in older women, which are more prevalent in this age demographic, where the morbidity risks of surgery are higher and life expectancy is likely to be relatively short [41]. Across different countries there is considerable variation in treatment strategies. However, at least across Europe little differences in survival were noted [42]. Although some trials have been aimed at improving treatment paradigms for older women [31] the evidence base is still limited for informed clinical decision making.

Data in Scottish women suggest that screening is an important indicator for BC prognosis regardless of molecular subtype [6, 43]. These results are in line with a recent study in Sweden reporting that women who participated in screening compared to those not participating in screening had a 41% reduction in the risk of BC death in the 10 years following diagnosis [44]. The potential for selection bias remains as women who accept invitations to screening are likely to differ from women that do not attend screening in ways that may influence survival.

Strengths and limitations

The strengths of this study are the high quality of the Scottish cancer registry data [14] and linkage to mortality records that provides the opportunity to assess survival and mortality by BC subtype. As data on ER have been collected from 1997 and on PR and HER2 from 2009 this gives one of the longest follow-up periods of European cohorts. Further, possible confounders, such as method of detection, comorbidities and deprivation are not always recorded in other European and North American cancer registries and were included in the survival models and used to assess the effect of age in BCSS by subtype. Particularly, our study is further adjusted for socioeconomic factors and comorbidities using the SIMD and the Charlson score, although its inclusion only attenuated the association between age and BCSS for HER2 + tumours, suggesting that other factors contribute to the associations between age and BCSS or all-cause mortality. Other major strength of this study is the investigation of the proportional hazards assumption using time-varying models and the use of competing risks analyses.

A limitation of our study is the possibility of unmeasured confounders such as lifestyle factors (alcohol and tobacco consumption, physical activity), reproductive factors or anthropometric factors (such as body mass index). Further, there was no detailed treatment data available, particularly the registry lacked information on the type of surgery women had (mastectomy or lumpectomy), targeted treatment for HER2 + BCs which has likely affected survival of women with this tumour subtype in the last decades, nor data on BC recurrence. Additionally, we note that missingness of important confounders, which is particularly true for tumour characteristics in women older than 70 years -since these women are less likely to have surgery- might have an effect on the estimates. Another limitation is the use of subgroup analyses, which reduce statistical power and increase the risk of type I and type II errors.

Although IHC-defined subtypes are reasonably good proxies for the gold-standard of RNA-expression defined subtypes, the lack of mRNA expression assays for the classification of the molecular subtypes constitutes a clear limitation. Future research should aim to explain the differences observed in survival by mRNA expression molecular subtypes and age groups.

Finally, survival can be affected by lead time and length biases usually caused by the introduction of a national screening programme during the period of study [14]. These biases can explain survival improvements that are related to screening and not to changes in the number of deaths that are prevented or delayed. Most cancer registries do not record data about method of detection of a cancer, however the availability in the Scottish cancer registry allowed us to investigate the effect of screening in BCSS.

Conclusion

Our study, using high quality population-based data from Scotland, highlights that age and molecular subtype are independent prognostic factors. Further, other tumour characteristics and treatments are also key prognostic factors for BCSS, regardless of age. Older women face a higher risk of breast cancer death even after adjusting for tumour characteristics, method of detection (screening), treatment, socio-economic factors and comorbidities. The increased mortality risk for women aged 70 years or older suggests that both biological factors and potential undertreatment may contribute to poorer outcomes. Additionally, competing causes of death, particularly cardiovascular diseases, were more prominent among older women (especially if they had a luminal A subtype) emphasizing the need for a comprehensive approach to patient management. These findings underscore the importance of age-specific strategies in breast cancer screening, treatment, and survivorship care to optimize outcomes for all patients. Scottish data have identified that many women with BC diagnosed over 70 years of age could potentially benefit from chemotherapy which require more systematic study [45]. In particular, it is not clear what measures of ageing, such as frailty, may inform who may or may not benefit from treatment. Future efforts need to further our understanding of factors that affect risk and benefits of treatments for different molecular subtypes to improve recommendations and outcomes (including patient reported outcomes) for the growing number of women with BC diagnosed over 70 years of age.

Data availability

The data that support the findings of this study are available through application to electronic Data Research and Innovation Service (eDRIS), a part of the Information Services Division of NHS Scotland. Restrictions apply to the availability of these data and any outputs presented here have been assessed and cleared for statistical disclosure by ISD Scotland. Outputs are however available from the authors upon reasonable request and with permission of ISD Scotland.

References

  1. Bray F, Laversanne M, Sung H, Ferlay J, Siegel RL, Soerjomataram I, Jemal A. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2024;74(3):229–63. https://doiorg.publicaciones.saludcastillayleon.es/10.3322/caac.21834.

  2. Cancer Research UK. Breast cancer statistics. CRUK, London. [cited 2025 March 15]; Available from: https://www.cancerresearchuk.org/health-professional/cancer-statistics/statistics-by-cancer-type/breast-cancer

  3. Coombs N, Storey A, Giddings R. Living longer: is age 70 the new age 65. Office for National Statistics. 2019. pp. 1–10. https://backup.ons.gov.uk/wp-content/uploads/sites/3/2019/11/Living-longer_-is-age-70-the-new-age-65_-.pdf

  4. Biganzoli L, Battisti NM, Wildiers H, McCartney A, Colloca G, Kunkler IH, Cardoso MJ, Cheung KL, De Glas NA, Trimboli RM, Korc-Grodzicki B. Updated recommendations regarding the management of older patients with breast cancer: a joint paper from the European Society of Breast Cancer Specialists (EUSOMA) and the International Society of Geriatric Oncology (SIOG). Lancet Oncol. 2021;22(7):e327–40.

  5. Early Breast Cancer Trialists’ Collaborative Group. Comparisons between different polychemotherapy regimens for early breast cancer: meta-analyses of long-term outcome among 100 000 women in 123 randomised trials. Lancet. 2012;379(9814):432–44.

  6. Mesa-Eguiagaray I, Wild SH, Rosenberg PS, Bird SM, Brewster DH, Hall PS, Cameron DA, Morrison D, Figueroa JD. 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(5):852–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41416-020-0938-z.

  7. Dawood S, Hu R, Homes MD, Collins LC, Schnitt SJ, Connolly J, Colditz GA, Tamimi RM. Defining breast cancer prognosis based on molecular phenotypes: results from a large cohort study. Breast Cancer Res Treat. 2011;126:185–92. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10549-010-1113-7.

  8. Hennigs A, Riedel F, Gondos A, Sinn P, Schirmacher P, Marmé F, Jäger D, Kauczor HU, Stieber A, Lindel K, Debus J. Prognosis of breast cancer molecular subtypes in routine clinical care: a large prospective cohort study. BMC Cancer. 2016;16:1–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12885-016-2766-3.

  9. Johansson AL, Trewin CB, Hjerkind KV, Ellingjord‐Dale M, Johannesen TB, Ursin G. Breast cancer‐specific survival by clinical subtype after 7 years follow‐up of young and elderly women in a nationwide cohort. Int J Cancer. 2019;144(6):1251–61. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/ijc.31950.

  10. Nielsen TO, Parker JS, Leung S, Voduc D, Ebbert M, Vickery T, Davies SR, Snider J, Stijleman IJ, Reed J, Cheang MC. A comparison of PAM50 intrinsic subtyping with immunohistochemistry and clinical prognostic factors in tamoxifen-treated estrogen receptor–positive breast cancer. Clin Cancer Res. 2010;16(21):5222–32. https://doiorg.publicaciones.saludcastillayleon.es/10.1158/1078-0432.CCR-10-1282.

  11. O’Brien KM, Cole SR, Tse CK, Perou CM, Carey LA, Foulkes WD, Dressler LG, Geradts J, Millikan RC. Intrinsic breast tumor subtypes, race, and long-term survival in the Carolina Breast Cancer Study. Clin Cancer Res. 2010;16(24):6100–10. https://doiorg.publicaciones.saludcastillayleon.es/10.1158/1078-0432.CCR-10-1533.

  12. Goldhirsch A, Wood WC, Coates AS, Gelber RD, Thürlimann B, Senn HJ. Strategies for subtypes—dealing with the diversity of breast cancer: highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2011. Ann Oncol. 2011;22(8):1736–47. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/annonc/mdr304.

  13. Public Health Scotland. Scottish Cancer Registry- How data are collected. 2020 [cited March]. Available from: https://www.isdscotland.org/Health-Topics/Cancer/Scottish-Cancer-Registry/How-data-are-collected/

  14. Brewster DH, Stockton DL. Ascertainment of breast cancer by the Scottish Cancer Registry: an assessment based on comparison with five independent breast cancer trials databases. Breast. 2008;17(1):104–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.breast.2007.07.044.

    Article  PubMed  Google Scholar 

  15. GDP Team-NHS National Services Scotland, Deprivation Guidance for analysts, Public Health & Intelligence, Editor. 2019. Available from: https://www.isdscotland.org/Products-and-Services/GPD-Support/Deprivation/_docs/PHI-Deprivation-Guidance-version-3-2.pdf

  16. Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis. 1987;40(5):373–83.

  17. Byrd DR, Brookland RK, Washington MK, Gershelwald JE, Compton CC, Hess KR, Sullivan DC, Jessup JM. AJCC cancer staging manual. In: Amin MB, Edge SB, Greene FL, editors. New York: Springer, 2017.

  18. Skyrud KD, Bray F, Møller B. A comparison of relative and cause‐specific survival by cancer site, age and time since diagnosis. Int J Cancer. 2014;135(1):196–203. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/ijc.28645.

  19. Kaplan EL, Meier P. Nonparametric estimation from incomplete observations. J Am Stat Assoc. 1958;53(282):457–81. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/01621459.1958.10501452.

  20. Cox DR. Regression models and life‐tables. J Royal Stat Soc Ser B-Statistical Methodol. 1972;34(2):187–202. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.2517-6161.1972.tb00899.x.

  21. Gore SM, Pocock SJ, Kerr GR. Regression models and non‐proportional hazards in the analysis of breast cancer survival. J R Stat Soc Ser C Appl Stat. 1984;33(2):176–95. https://doiorg.publicaciones.saludcastillayleon.es/10.2307/2347444.

  22. Bellera CA, MacGrogan G, Debled M, De Lara CT, Brouste V, Mathoulin-Pélissier S. Variables with time-varying effects and the Cox model: some statistical concepts illustrated with a prognostic factor study in breast cancer. BMC Med Res Methodol. 2010;10:1–2. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/1471-2288-10-20.

  23. Kalbfleisch JD, Prentice RL. The statistical analysis of failure time data. John Wiley & Sons; 2002.

  24. Fine JP, Gray RJ. A proportional hazards model for the subdistribution of a competing risk. J Am Stat Assoc. 1999;94(446):496–509. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/01621459.1999.10474144.

  25. Stata. Stata statistical software release 15. StataCorp LLC.: College Station, TX; 2017.

    Google Scholar 

  26. Howlader N, Cronin KA, Kurian AW, Andridge R. Differences in breast cancer survival by molecular subtypes in the United States. Cancer Epidemiol Biomarkers Prev. 2018;27(6):619–26. https://doiorg.publicaciones.saludcastillayleon.es/10.1158/1055-9965.EPI-17-0627.

  27. Fallahpour S, Navaneelan T, De P, Borgo A. Breast cancer survival by molecular subtype: a population-based analysis of cancer registry data. CMAJ. 2017;5(3):E734–9. https://doiorg.publicaciones.saludcastillayleon.es/10.9778/cmajo.20170030.

  28. Caldarella A, Crocetti E, Bianchi S, Vezzosi V, Urso C, Biancalani M, Zappa M. Female breast cancer status according to ER, PR and HER2 expression: a population based analysis. Pathol Oncol Res. 2011;17:753–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s12253-011-9381-z.

  29. Puig-Vives M, Sánchez MJ, Sánchez-Cantalejo J, Torrella-Ramos A, Martos C, Ardanaz E, Chirlaque MD, Perucha J, Díaz JM, Mateos A, Machón M. Distribution and prognosis of molecular breast cancer subtypes defined by immunohistochemical biomarkers in a Spanish population-based study. Gynecol Oncol. 2013;130(3):609–14. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.ygyno.2013.05.039.

  30. Ess SM, Herrmann C, Bouchardy C, Neyroud I, Rapiti E, Konzelmann I, Bordoni A, Ortelli L, Rohrmann S, Frick H, Mousavi M. Impact of subtypes and comorbidities on breast cancer relapse and survival in population-based studies. Breast. 2018;41:151–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.breast.2018.07.011.

  31. Kunkler IH, Williams LJ, Jack WJ, Cameron DA, Dixon JM. Breast-conserving surgery with or without irradiation in women aged 65 years or older with early breast cancer (PRIME II): a randomised controlled trial. Lancet Oncol. 2015;16(3):266–73.

  32. Hooning MJ, Botma A, Aleman BM, Baaijens MH, Bartelink H, Klijn JG, Taylor CW, Van Leeuwen FE. Long-term risk of cardiovascular disease in 10-year survivors of breast cancer. J Natl Cancer Inst. 2007;99(5):365–75. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jnci/djk064.

  33. Zagar TM, Cardinale DM, Marks LB. Breast cancer therapy-associated cardiovascular disease. Nat Rev Clin Oncol. 2016;13(3):172–84. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nrclinonc.2015.171.

  34. Bertagnolli MM, Singh H. Treatment of older adults with cancer—addressing gaps in evidence. N Engl J Med. 2021;385(12):1062–5. https://doiorg.publicaciones.saludcastillayleon.es/10.1056/NEJMp2106089.

  35. Srour B, Hynes LC, Johnson T, Kühn T, Katzke VA, Kaaks R. Serum markers of biological ageing provide long-term prediction of life expectancy—a longitudinal analysis in middle-aged and older German adults. Age Ageing. 2022;51(2):afab271. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/ageing/afab271.

  36. Mandelblatt JS, Cai L, Luta G, Kimmick G, Clapp J, Isaacs C, Pitcher B, Barry W, Winer E, Sugarman S, Hudis C. Frailty and long-term mortality of older breast cancer patients: CALGB 369901 (Alliance). Breast Cancer Res Treat. 2017;164:107–17. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10549-017-4222-8.

  37. Mak JK, Kuja-Halkola R, Bai G, Hassing LB, Pedersen NL, Hägg S, Jylhävä J, Reynolds CA. Genetic and environmental influences on longitudinal frailty trajectories from adulthood into old age. J Gerontol A. 2023;78(2):333–41. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/gerona/glac197.

  38. Kresovich JK, O’Brien KM, Xu Z, Weinberg CR, Sandler DP, Taylor JA. Changes in methylation-based aging in women who do and do not develop breast cancer. J Natl Cancer Inst. 2023;115(11):1329–36. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jnci/djad117.

  39. Rentscher KE, Bethea TN, Zhai W, Small BJ, Zhou X, Ahles TA, Ahn J, Breen EC, Cohen HJ, Extermann M, Graham DM. Epigenetic aging in older breast cancer survivors and noncancer controls: preliminary findings from the Thinking and Living with Cancer Study. Cancer. 2023;129(17):2741–53. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/cncr.34818.

  40. Pepping RM, Portielje JE, Van De Water W, De Glas NA. Primary endocrine therapy in older women with breast cancer. Curr Geriatr Rep. 2017;6:239–46. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s13670-017-0223-z.

  41. Reed MW, Audisio RA, editors. Management of breast cancer in older women. Springer Science & Business Media; 2010.

  42. Derks MG, Bastiaannet E, Kiderlen M, Hilling DE, Boelens PG, Walsh PM, Van Eycken E, Siesling S, Broggio J, Wyld L, Trojanowski M. Variation in treatment and survival of older patients with non-metastatic breast cancer in five European countries: a population-based cohort study from the EURECCA Breast Cancer Group. Br J Cancer. 2018;119(1):121–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41416-018-0090-1.

  43. Mesa-Eguiagaray I, Wild SH, Bird SM, Williams LJ, Brewster DH, Hall PS, Figueroa JD. Breast cancer incidence and survival in Scotland by socio-economic deprivation and tumour subtype. Breast Cancer Res Treat. 2022;194(2):463–73. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10549-022-06632-1.

  44. Duffy SW, Tabár L, Yen AM, Dean PB, Smith RA, Jonsson H, Törnberg S, Chen SL, Chiu SY, Fann JC, Ku MM. Mammography screening reduces rates of advanced and fatal breast cancers: results in 549,091 women. Cancer. 2020;126(13):2971–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/cncr.32859.

  45. Castresana-Aguirre M, Johansson A, Matikas A, Foukakis T, Lindström LS, Tobin NP. Clinically relevant gene signatures provide independent prognostic information in older breast cancer patients. Breast Cancer Res. 2024;26(1):38. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13058-024-01797-7.

Download references

Acknowledgements

We thank NHS Scotland and Information Services Division for collection of data used for this analysis. We specifically would like to acknowledge Andrew Deas, Suhail Iqbal, Rita Nogueira and Ross Murdoch of NHS National Services Scotland, Jeremy Thomas, Catherine Dhaliwal and Mike Dixon for helpful discussions of these analyses. JDF is currently funded by the Intramural program of the National Cancer Institute Division of Cancer Epidemiology and Genetics.

Funding

JDF is currently funded by the Intramural program, Division of Cancer Epidemiology and Genetics the National Cancer Institute, National Institutes of Health USA. This project was funded by Wellcome Trust grant 207800/Z/17/Z.

Author information

Authors and Affiliations

Authors

Contributions

Conception & Design of the study: IME, JDF, SHW and LW. Interpretation of data: All authors. Drafting of the manuscript: IME, JDF and SHW. Revised work and provided important intellectual content: All authors. Final approval of the manuscript: All authors.

Corresponding authors

Correspondence to Ines Mesa-Eguiagaray or Jonine D. Figueroa.

Ethics declarations

Ethics approval and consent to participate

Approval from the Public Benefit and Privacy Panel for Health (PBPP) and Social Care is a requirement for data access. Our project was approved by PBPP reference number 1718-0057.

Consent for publication

Not applicable.

Competing interests

SMB holds shares in GlaxoSmithKline. Other authors declare no competing interests.

Additional information

Publisher’s note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mesa-Eguiagaray, I., Wild, S.H., Williams, L.J. et al. Breast cancer-specific survival by molecular subtype in different age groups of women in Scotland. Breast Cancer Res 27, 59 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13058-025-02012-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13058-025-02012-x

Keywords