Skip to main content

Clinically relevant gene signatures provide independent prognostic information in older breast cancer patients

Abstract

Background

The clinical utility of gene signatures in older breast cancer patients remains unclear. We aimed to determine signature prognostic capacity in this patient subgroup.

Methods

Research versions of the genomic grade index (GGI), 70-gene, recurrence score (RS), cell cycle score (CCS), PAM50 risk-of-recurrence proliferation (ROR-P), and PAM50 signatures were applied to 39 breast cancer datasets (N = 9583). After filtering on age ≥ 70 years, and the presence of estrogen receptor (ER) and survival data, 871 patients remained. Signature prognostic capacity was tested in all (n = 871), ER-positive/lymph node-positive (ER + /LN + , n = 335) and ER-positive/lymph node-negative (ER + /LN−, n = 374) patients using Kaplan–Meier and multivariable Cox-proportional hazard (PH) modelling.

Results

All signatures were statistically significant in Kaplan–Meier analysis of all patients (Log-rank P < 0.001). This significance remained in multivariable analysis (Cox-PH, P ≤ 0.05). In ER + /LN + patients all signatures except PAM50 were significant in Kaplan–Meier analysis (Log-rank P ≤ 0.05) and remained so in multivariable analysis (Cox-PH, P ≤ 0.05). In ER + /LN− patients all except RS were significant in Kaplan–Meier analysis (Log-rank P ≤ 0.05) but only the 70-gene, CCS, ROR-P, and PAM50 signatures remained so in multivariable analysis (Cox-PH, P ≤ 0.05).

Conclusions

We found that gene signatures provide prognostic information in survival analyses of all, ER + /LN + and ER + /LN- older (≥ 70 years) breast cancer patients, suggesting a potential role in aiding treatment decisions in older patients.

Background

Human life expectancy is predicted to increase globally by 4.4 years for both men and women in the coming two decades [1]. This will result in a larger population of older adults and as cancer is generally a disease of aging, it is estimated that 60% of newly diagnosed cancers in 2035 will come from adults aged 65 + [2]. Older cancer patients are however typically underrepresented in clinical trials [3,4,5,6] and may also be undertreated relative to younger patient populations [7, 8]. This implies that there is a lack of data on whether the tools used to guide treatment decisions in younger (< 70 years) cancer patients are also applicable to older patient populations.

Treatment decisions in early breast cancer are based on tumour size, lymph node involvement, stage and prognostic and predictive biomarkers including the estrogen, progesterone and human epithelial growth factor 2 receptors (ER, PR and HER2) and expression of the proliferation marker Ki67 [9]. A patient’s age is also recommended to be taken into consideration under the provision that it should not be used as a reason to withhold specific treatments [9]. This is in line with data from the Early Breast Cancer Trialists Collaborative Group (EBTCG) showing that the relative benefit from chemotherapy is independent of age [10], and international treatment guidelines for older breast cancer patients (> 70 years) which state that endocrine treatment should be offered to postmenopausal women irrespective of age [11]. It is important to note however that treatments should only be offered after an initial screening assessment for frailty [11].

In addition to routine pathological staging and clinical biomarkers, recent years have also seen an increase in the use of multigene signatures to aid in risk stratification of early breast cancer patients. The signature field is most mature in breast cancer with the 21-gene recurrence score (RS) and 70-gene signatures demonstrating prognostic capacity in large scale randomized clinical trials [12,13,14]. Moreover, these signatures along with others can be used to guide treatment decisions primarily in postmenopausal early breast cancer patients with node negative or positive (1–3 nodes) invasive tumours [15]. The use of gene signatures in older breast cancer patients remains controversial however as there is currently insufficient evidence to support their use [11]. As such, we aimed to perform the first comprehensive comparison of the additional prognostic capacity of clinically relevant breast cancer gene signatures beyond that of routine biomarkers in a single older breast cancer cohort. Specifically, we apply the genomic grade index (GGI), 70-gene, 21-gene recurrence score (RS), cell cycle score (CCS), PAM50 Risk of Recurrence score—Proliferation (ROR-P), and PAM50 signatures to 39 open access breast cancer datasets with a combined total of 871 patients aged 70 years or older.

Methods

Cohort description

The data for this study was extracted from the R package MetaGxBreast [16], a gene expression database of 39 open access breast cancer datasets with manually-curated and standardized clinical, pathological, survival, and treatment metadata for breast cancer totalling 9583 patients. After first selecting patients who were aged 70 years or older (n = 1399), samples were subsequently excluded on the basis of: lacking information on ER status (n = 62), lacking survival information (n = 323), or insufficient coverage of gene signature genes (n = 143, further described below), 871 patients remained in total (see CONSORT diagram in Fig. 1). We further stratified these patients into clinically relevant subgroups of ER + /LN + (N = 335) and ER + /LN− (N = 374) to analyse the prognostic capacity of gene signatures taking ER and lymph node statuses into account. Information regarding the number of positive lymph nodes was not available. Importantly, in this study we define older patients as those over the age of 69 years, in line with previous publications [17,18,19,20,21] and randomized clinical trials[22, 23]. In order to compare older patients to a younger postmenopausal breast cancer patient population we also selected ER + /LN− patients between 55 and 65 years of age (N = 478, labelled as “55–65” in subsequent analyses). This 5-year gap between 65 and 70 years was selected to have a clear separation by age and the 10 year interval was chosen in order to obtain a large enough cohort for comparison.

Fig. 1
figure 1

Consort diagram of older breast cancer patient selection from the 39 datasets in the MetaGxBreast database. Patients were excluded owing to lacking information on ER, being under the age of 70 years old, lacking survival information or having insufficient coverage of gene signature genes. * See methods for a full description

Gene expression signatures

Research versions of the Genomic Grade Index (GGI) [24], 70 gene [25], 21-gene recurrence score (RS) [26], Cell Cycle Score (CCS) [27]), PAM50 Risk of Recurrence score—Proliferation (ROR-P), and Prediction Analysis of Microarray 50 (PAM50) [28], were applied on an individual dataset basis to each of the 39 MetaGxBreast studies as described in the original publications. Signature classifications were subsequently pooled with clinical data for statistical analysis. For GGI, tumour grade was not available for all patients (missing in 38%), as such, we used a variation of the tool to compute the tumour grade [29] and then ran the conventional research version of GGI. The original RS signature cutoffs (Low < 18, Intermediate 18–31, High > 31) were used throughout the study, however the updated cutoffs from the TAILORx clinical trial (Low < 10, Intermediate 11–25, High > 26) were also assessed for the sake of completeness, as specified in the results section. PAM50 and ROR-P work most accurately if the dataset is ER status-balanced [30], however population based breast cancer datasets have a higher proportion of ER + tumours (as this tumour type is more common). To address this skew in the MetaGxBreast studies (ER proportions shown in Additional file 1: Supplementary Table 1) we used Monte Carlo sampling to have equal proportions of ER + and ER− patients per dataset. Specifically, we compute the frequencies of ER + and ER-patients and then subsampled patients from the predominant group to ensure that the frequencies of ER + and ER− cases were equal within datasets. This balancing process involved conducting 100 subsamples for each dataset, thereby facilitating the computation of the median expression levels of the probes for data centering. This methodology has been previously applied by us and others [31,32,33] and full code to reproduce our signature calls and statistical analyses can be found in the following code repository (https://bitbucket.org/tobingroup/elderly). The gene signatures were chosen owing to their relevance in randomized clinical trials such as MINDACT [13], ASTER70 [22], TAILORX [34] and also in real world evaluation [35]. This article was performed and is reported in accordance with REMARK guidelines [36].

Gene expression and gene mapping

Gene expression data extracted from MetaGxBreast comes pre-processed and normalized, a detailed description of this can be found in the original publication [16]. Probe to gene mapping was achieved by merging annotation sources from MetaGxBreast, supplementary files from the original signature publications and Bioconductor 3.15 in R. Probes mapping to the same gene were combined by averaging their expression values. MetaGxBreast combines datasets from different gene expression array platforms, the majority of which are Affymetrix. This means that all genes are not found in all datasets. The 70-gene signature is derived on an Agilent array platform, and approximately 75% of the signature’s genes are mappable to the Affymetrix platform. For this reason, we excluded datasets if less than 75% of the 70-gene signature probes were present and similarly if any non-reference genes for the RS were absent. Consequently, the median gene coverage for GGI was 100%, 100% for PAM50, 91% for the 70-gene signature, 100% for RS, and 96% for CCS.

Statistical analysis

Kaplan–Meier and multivariable Cox proportional hazard analyses were used to assess older patient survival in the context of gene expression signature subgroups. We used the R package survminer version 0.4.9 for Kaplan–Meier and the R package survival version 3.4–0 for multivariable Cox proportional hazard analyses. The latter was adjusted for ER status, lymph node status, tumour grade, tumour size and whether the patients had received hormonal therapy or not. We did not adjust for treatment with chemotherapy as few patients received it (N = 50). These methodologies were applied to All patients as well the subgroups of ER + /LN+ and ER + /LN – patients. Subgroups were adjusted for tumour size, grade, and hormonal therapy only. The clinical endpoint used was Recurrence Free Survival (RFS) defined as the time from date of curative surgery to the time recurrence (distant metastatic events and loco-regional recurrences). RFS was censored at 10 years and the median follow-up time was 6.2 years. MetaGxBreast does not provide RFS data for METABRIC, so this information was instead extracted from the supplementary data of Rueda et al. [37]. The likelihood ratio (LR) and the concordance index (c-index) were computed using univariable models as a measure of signature prognostic capacity. In general, the LR is a useful metric for comparing variables with differing numbers of subgroups (e.g. GGI has two subgroups and PAM50 has five) as it provides a single overall value alongside a χ2 test statistic. It has been previously used by us and others for head to head signature and immunohistochemical biomarker comparisons [27, 33, 38, 39]. For determination of the additional prognostic capacity of signatures beyond clinico-pathological markers we calculated the delta likelihood ratio (∆LR), by comparing the LR of a multivariable model that included the adjustment variables noted above with and without the gene expression signature. This allows us to explicitly quantify the prognostic capacity of the models when gene signatures are included and more directly compare the prognostic value of each signature. To assess if there was any statistical difference in the classification of patients by gene signatures between younger postmenopausal patients (between 55 and 65 years of age) and older patients, we performed Chi-squared tests. Tests used are indicated in table legends. All statistical tests were two-sided and a significance level ⍺ of 5% was used. All statistical analyses were performed using R statistical software version 4.1.2.

Results

Cohort clinico-pathological characteristics and gene signature distribution

To assess the prognostic capacity of gene expression signatures in older (≥ 70 years) breast cancer patients, we applied research versions of the GGI, 70 gene, RS, CCS, ROR-P, and PAM50 signatures to expression array data from 39 open access breast cancer datasets individually. Signature classifications and clinico-pathological variables were then merged into a single dataset. After limiting the cohort to an older population only and filtering on the basis of exclusion criteria (see methods), 871 patients remained (CONSORT diagram, Fig. 1). The median age of these patients was 75 years old (range 70–96) and as expected, a decrease in age frequency is readily apparent as patients tend towards the upper age range (Additional file 2: Supplementary Fig. 1). Clinico-pathological characteristics for these patients are shown in Table 1. The majority of patients were ER-positive (87%), PR-positive (47%) and HER2-negative (65%, Table 1, PR and HER2 status unknown in 23 and 27% of patients, respectively). Half of all patients were negative for lymph node metastases and the majority of tumours were of larger size (69% ≥ 2 cm) and intermediate or high grade (81%). 57% of patients received hormonal therapy but few received chemotherapy (< 6%), as anticipated given that many of the datasets used in this study are from cohorts assembled before the 2000s where treatments were less standardized. The number of patients categorized into subgroups on the basis of gene expression signatures is shown in Additional file 3: Supplementary Table 2. Binary signatures demonstrated an approximate even split into grouping patients into low and high-risk groups (see GGI and 70-gene, Additional file 3: Supplementary Table 2). No common pattern of sample distribution was noted for signatures with three levels (RS, CCS and ROR-P) but PAM50 classified 80% of patients into luminal A or B subtypes, in line with the high level of ER-positive patients in the cohort.

Table 1 Clinico-pathological characteristics of the older breast cancer cohort

Gene signatures provide independent prognostic information for older breast cancer patients

We next assessed the prognostic capacity of the GGI, 70-gene, RS, CCS, ROR-P, and PAM50 gene signatures in All (N = 871), ER + /LN + (N = 335), and ER + /LN− patients (N = 374) using Kaplan–Meier analysis. All signatures provided prognostic information in All patients (log-rank test, P < 0.05, Fig. 2). In the ER + /LN + subgroup, all signatures except PAM50 were statistically significant (log-rank test, P < 0.05, Additional file 2: Supplementary Fig. 2). In the ER + /LN— subgroup, the GGI, 70-gene, CCS, ROR-P and PAM50 signatures were statistically significant (log-rank test, P < 0.05, Additional file 2: Supplementary Fig. 3) with a non-significant trend observed for RS (log-rank test, P = 0.068, Additional file 2: Supplementary Fig. 3). In general, statistically significant gene signatures remained so in multivariable Cox proportional hazards analyses in All and ER + /LN + subgroups after adjusting for tumour size, tumour grade, ER status, lymph node status, and whether the patient received hormonal therapy or not (Cox proportional hazards modeling, P < 0.05 vs. signature reference group. Table 2). In ER + /LN− patients however, only the 70-gene, CCS, ROR-P and PAM50 signatures remained statistically significant in the same analysis (Table 2). Of note, we also analyzed ER + /LN−/HER2− patients and observed similar trends (Additional file 3: Supplementary Table 3) as for ER + /LN− patients but only the 70-gene signature and PAM50 showed statistically significant independent prognostic information (P = 0.02 for 70-gene high vs. low risk and P = 0.002 for PAM50 Her2-enriched vs Luminal A). We also applied Likelihood ratio statistics to assess the additional prognostic capacity of signatures beyond the routine clinico-pathological markers ER, lymph node status, tumour grade, tumour size and hormonal therapy received or not. This allows for a more direct comparison of signature prognostic capacity relative to each other. The largest ∆LRs were found for the RS, CCS and ROR-P signatures in All (∆LR = 14.11, 16.97, and 15.38) and ER + /LN + (∆LR = 13.6, 13.2 and 15.82) patients, respectively (LR−test, P < 0.01, Additional file 3: Supplementary Table 4). Lower ∆LRs were found for signatures in ER + /LN− patients where only the 70 gene and CCS remained statistically significant (∆LR = 6.18 and 6.47 for the 70 gene and CCS signatures respectively, LR-test, P < 0.05, Additional file 3: Supplementary Table 4).

Fig. 2
figure 2

Kaplan–Meier analysis of gene expression signatures in All patients of the older cohort. a Genomic Grade Index (GGI) b 70-gene c recurrence score (RS) d cell-cycle score (CCS) e PAM50 risk of recurrence score—proliferation (ROR-P) f prediction analysis of microarray 50 (PAM50)

Table 2 Multivariable analysis for gene signatures in All, ER + /LN + and ER + /LN− older patients

Differences in signature risk stratification in comparison to a younger postmenopausal patient cohort

The clinical utility of gene signatures is most frequently discussed in relation to ER + /LN− postmenopausal breast cancer patients, however in our analysis all signatures showed reduced prognostic capacity (∆LR) in this patient subgroup. To understand if there are differences in signature prognostic capacity between an older ER + /LN− patient cohort and a younger postmenopausal ER + /LN− cohort, we ran the same analyses but this time selecting ER + /LN− patients between 55 and 65 years old (N = 478). Patient characteristics for this subgroup are shown in Additional file 3: Supplementary Table 5. We found that all gene signatures except PAM50 provided independent prognostic information in ER + /LN− patients aged 55–65 years old (Additional file 3: Supplementary Table 6) and that the additional prognostic capacity of signatures beyond routine clinico-pathological markers were in general higher in this younger subgroup relative to older (≥ 70 years) ER + /LN− breast cancer patients (55–65 vs > 70, ∆LR; GGI: 6.23 vs 1.43; 70-gene: 10.11 vs 6.18; RS: 6.59 vs 1.84; CCS: 8.93 vs 6.47; ROR-P: 12.02 vs 5.25; PAM50: 4.58 vs 7.89, Additional file 3: Supplementary Tables 4 and 6). Of note, we also tested the TAILORx RS cutoffs in the same subgroups (55–65 and > 70 ER + /LN− patients) but found it to perform worse than the original RS cutoffs—it did not provide independent prognostic information in either subgroup (data not shown).

To assess whether the generally higher LRs in 55–65 year olds could be owing to a difference in signature risk stratification (assignment of tumours into signature subgroups) we also compared signature subgroup composition between older and younger ER + /LN− patients. We found that only the ROR-P signature showed a statistically significant difference in risk stratification between these two groups (Chi-square test, P = 0.02, Additional file 3: Supplementary Table 7). This suggests that the difference in signature prognostic capacity between older and younger patients is unlikely to be related to signature stratification and may point to a need for re-optimization of signatures and their prognostic cutoffs in older breast cancer populations.

Discussion

In this study we assessed the prognostic capacity of six gene expression signatures in a cohort of 871 older (≥ 70 years) breast cancer patients. We found that all gene signatures provided independent prognostic information in All patients. In ER + /LN + patients all but PAM50 provided independent prognostic information while in ER + /LN− patients the 70-gene, CCS, ROR-P and PAM50 signatures remained statistically significant after adjusting for routine clinico-pathological variables. LR statistics showed lower additional prognostic capacity of signatures beyond these routine variables in ER + /LN− patients relative to All and ER + /LN + patients. Further comparison of older ER + /LN− patients to a postmenopausal ER + /LN− cohort of younger age (55–65) showed higher signature ∆LR values in the younger cohort despite similar signature subgroup stratification in both groups.

This is the first comprehensive study comparing the prognostic performance of multiple clinically relevant gene expression signatures in a single older breast cancer patient cohort, adjusting for clinco-pathological tumour and patient characteristics. Two studies have, however, focused on the prognostic or treatment predictive capacity of individual signatures in older breast cancer patients. In the first, the ASTER 70 s randomized phase III clinical trial demonstrated that GGI is prognostic in older (≥ 70 years) breast cancer patients [40] with a similar age range to our study. The trial did not however find a statistically significant overall survival benefit with the addition of chemotherapy to endocrine therapy after surgery in ER + /HER2− patients whose tumours were classified as GGI high (GG3) [40]. In the second, Noordhoek et al. found that patients with high clinical risk based on the St. Gallen risk classification, but classified as ultra-low risk by the 70-gene signature, had excellent prognosis [41] in a study of 418 older patients above 70 years old (median age 78). Other studies whilst not focusing on prognosis or treatment predictive value, do provide some indications that gene signatures may be useful in older patients. A recent study of the 70-gene signature in 89 older patients showed no statistical difference in the proportion of patients classified into low and high-risk groups when compared to the MINDACT clinical trial implying that this signature could also be applied to older breast cancer patients [42]. Regarding RS, Iles et al. showed a decline in its usage with increasing age and a higher prevalence of low-risk classifications in patients above 70 years old [43]. ER + /LN−/HER2− patients derive low or no benefit from chemotherapy if a tumour is classified as RS low [34] and in line with this, Barni et al. found a 38% relative reduction of chemotherapy usage in 230 older patients on the basis of RS [44]. Whilst the above studies focus on de-escalation, age alone should not be a contraindication for chemotherapy usage [45,46,47]. This implies that running gene signatures on older breast cancer patients could also identify those who would benefit from chemotherapy treatment but may not in the past have been treated owing to advanced age. RS has proved useful for escalation of chemotherapy in older patients, with a treatment decision change of 18.6% mainly from no-chemotherapy to chemotherapy in a cohort of 237 patients [48]. Interestingly, in our study gene expression signatures classified between 25 and 53% of all tumours into high-risk groups however, only 6% of the patients received chemotherapy, none of whom were ER + /LN−.

Breast cancer clinical oncology guidelines support the use of the 70-gene (commercially MammaPrint) and RS (OncotypeDx) signatures to guide endocrine or chemotherapy treatment decisions in postmenopausal ER + /LN−/HER2− patients (or with 1–3 positive lymph nodes) and PAM50 (Prosigna) in postmenopausal ER + /LN−/HER2- patients [9, 15]. In our study, the 70-gene and PAM50 signatures provided independent prognostic information in ER + /LN− older patients (N = 374) and showed similar trends in ER + /LN−/HER2— patients (P ≤ 0.05, N = 222). However, we did note that the additional prognostic capacity of signatures beyond routine clinico-pathological markers was reduced in our older ER + /LN− patient cohort relative to those aged 55–65. The reason for this is not readily apparent but one potential concern is that the cutoffs for signature classification into prognostic groups were optimized on a younger patient population and may need to be re-optimised/changed for application to an older cohort. Related to this, Jezequel et al. noted a difference in the proportions of patients classified into good/poor prognosis groups by GGI, 70-gene and RS when comparing the age groups ≤ 40, 40–70, and ≥ 70[49]. Similarly, Kruijf et al. found differences in PAM50 subtype proportions and weaker signature prognostic capacity in older patients (≥ 65) when compared to a younger (< 65) patient population [50]. In addition, an increase in the proportion of tumours classified as luminal subtypes and a decrease in basal-like subtypes has also been previously reported in older patients [49,50,51]. Taken together, these studies imply that biomarkers and cutpoints used in younger postmenopausal breast cancer population might not be directly transferable to older patients without modification. In the current study only the ROR-P signatures showed a statistically significant difference in risk stratification when comparing older vs. younger postmenopausal patients (Additional file 3: Supplementary Table 7), but non-significant trends were noted for GGI and PAM50. As such these differences are unlikely to explain the reduced additional prognostic capacity of signatures in ER + /LN− older patients and this reduction could be owing to a difference in the biology of tumours from older populations. This is also supported by one study showing that luminal B tumours from patients over the age of 70 years were less aggressive than those from younger age groups and that this was related to differences in pathways for iron metabolism, mitochondrial oxidative phosphorylation and reactive stroma [49].

The strengths of our study are as follows. First, we provide a comprehensive analysis of the prognostic capacity of six gene expression signatures in an older breast cancer patient cohort with a median age of 75.4 years old (52% of patients were over the age of 75). This is notable as no patients above the age of 75 were included in the TAILORx or MINDACT clinical trials and only 12% of patients were above 70 years old in the RxPONDER trial [14], despite 30% of the breast cancer diagnosis occurring in patients above 70 [52]. Second, we assess the additional prognostic capacity of these signatures beyond routine clinico-pathological biomarkers—something that is currently lacking in published literature. There are also some limitations to this study. First, no clear definition of what constitutes a patient as “older” is routinely applied in a clinical setting. Even though we defined this as patients ≥ 70 years of age, the usage of chronological age may ignore the diverse ways that time affects individuals. Since cancer is a disease of aging [53], a better definition of older could be obtained using the biological age which takes multiple biological and physiological developmental factors including genetics, lifestyle, diet and comorbidities [54]. Related to this, a second limitation is that we did not adjust our analyses for patient frailty which is known to negatively impact prognosis [55], owing to no data on this being available for this metric. Third, clinical survival endpoints were not identical across the 39 independent datasets we assessed; therefore, we combined two different end points (RFS and DMFS) into a single survival metric. Of note, we chose not to use an overall survival (OS) endpoint owing to potential competing causes of death. Fourth, we relied on the research versions of gene-expression signatures in place of their commercial implementations, fifth, we lack patient numbers to assess the treatment predictive value of these signature and sixth the CCS risk stratification cutoffs were derived from the METABRIC dataset which is included in our study. Whilst we have never optimized the CCS cutoffs for prognostic capacity, there is still the potential for an overfit bias that could possibly impact our results for this signature only.

In conclusion, we show that gene expression signatures provide independent prognostic information in All, ER + /LN + and ER + /LN− patients who are over the age of 70 years, supporting the rationale of the ASTER70s clinical trial. These results suggest a potential role for gene expression signatures in aiding treatment decisions in older breast cancer patients and indicate that further investigation is warranted in prospective clinical study to elucidate their treatment predictive value.

Data availability

The data used in this study are publicly available on MetaGxBreast package in R (https://bioconductor.org/packages/release/data/experiment/html/MetaGxBreast.html). R-code to reproduce the results of this study are publicly available at https://bitbucket.org/tobingroup/elderly.

Abbreviations

GGI:

Genomic grade index

RS:

Recurrence score

CCS:

Cell cycle score

ROR-P:

PAM50 risk of recurrence score-proliferation

ER:

Estrogen receptor

LN:

Lymph node

RFS:

Recurrence free survival

PR:

Progesterone receptor

HER2:

Human epidermal growth factor receptor 2

LR:

Likelihood ratio

DMFS:

Distant metastasis free survival

OS:

Overall survival

References

  1. Foreman KJ, Marquez N, Dolgert A, Fukutaki K, Fullman N, McGaughey M, et al. Forecasting life expectancy, years of life lost, and all-cause and cause-specific mortality for 250 causes of death: reference and alternative scenarios for 2016–40 for 195 countries and territories. Lancet. 2018;392:2052–90.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Pilleron S, Sarfati D, Janssen-Heijnen M, Vignat J, Ferlay J, Bray F, et al. Global cancer incidence in older adults, 2012 and 2035: a population-based study. Int J Cancer. 2019;144:49–58.

    Article  CAS  PubMed  Google Scholar 

  3. Townsley CA, Selby R, Siu LL. Systematic review of barriers to the recruitment of older patients with cancer onto clinical trials. J Clin Oncol. 2005;23:3112–24.

    Article  PubMed  Google Scholar 

  4. Hutchins LF, Unger JM, Crowley JJ, Coltman CA Jr, Albain KS. Underrepresentation of patients 65 years of age or older in cancer-treatment trials. N Engl J Med. 1999;341:2061–7.

    Article  CAS  PubMed  Google Scholar 

  5. Lewis JH, Kilgore ML, Goldman DP, Trimble EL, Kaplan R, Montello MJ, et al. Participation of patients 65 years of age or older in cancer clinical trials. J Clin Oncol. 2003;21:1383–9.

    Article  PubMed  Google Scholar 

  6. Talarico L, Chen G, Pazdur R. Enrollment of elderly patients in clinical trials for cancer drug registration: a 7-year experience by the US Food and Drug Administration. J Clin Oncol. 2004;22:4626–31.

    Article  PubMed  Google Scholar 

  7. Bouchardy C, Rapiti E, Fioretta G, Laissue P, Neyroud-Caspar I, Schäfer P, et al. Undertreatment strongly decreases prognosis of breast cancer in elderly women. J Clin Oncol. 2003;21:3580–7.

    Article  PubMed  Google Scholar 

  8. Noon AP, Albertsen PC, Thomas F, Rosario DJ, Catto JWF. Competing mortality in patients diagnosed with bladder cancer: evidence of undertreatment in the elderly and female patients. Br J Cancer. 2013;108:1534–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Cardoso F, Kyriakides S, Ohno S, Penault-Llorca F, Poortmans P, Rubio IT, et al. Early breast cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up†. Ann Oncol. 2019;30:1194–220.

    Article  CAS  PubMed  Google Scholar 

  10. Peto R, Davies C, Godwin J, Gray R, Pan HC, et al. 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:432–44.

    Article  Google Scholar 

  11. Biganzoli L, Battisti NML, Wildiers H, McCartney A, Colloca G, Kunkler IH, et al. 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:e327–40.

    Article  PubMed  Google Scholar 

  12. Sparano JA, Gray RJ, Ravdin PM, Makower DF, Pritchard KI, Albain KS, et al. Clinical and genomic risk to guide the use of adjuvant therapy for breast cancer. N Engl J Med. 2019;380:2395–405.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Cardoso F, van’t Veer LJ, Bogaerts J, Slaets L, Viale G, Delaloge S, et al. 70-Gene signature as an aid to treatment decisions in early-stage breast cancer. N Engl J Med. 2016;375:717–29.

    Article  CAS  PubMed  Google Scholar 

  14. Kalinsky K, Barlow WE, Gralow JR, Meric-Bernstam F, Albain KS, Hayes DF, et al. 21-Gene assay to inform chemotherapy benefit in node-positive breast cancer. N Engl J Med. 2021;385:2336–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Andre F, Ismaila N, Allison KH, Barlow WE, Collyar DE, Damodaran S, et al. Biomarkers for adjuvant endocrine and chemotherapy in early-stage breast cancer: ASCO guideline update. J Clin Oncol. 2022;40:1816–37.

    Article  CAS  PubMed  Google Scholar 

  16. Gendoo DMA, Zon M, Sandhu V, Manem VSK, Ratanasirigulchai N, Chen GM, et al. MetaGxData: clinically annotated breast, ovarian and pancreatic cancer datasets and their use in generating a multi-cancer gene signature. Sci Rep. 2019;9:8770.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  17. Derks MGM, Bastiaannet E, Kiderlen M, Hilling DE, Boelens PG, Walsh PM, et al. 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:121–9.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Wildiers H, de Glas NA. Anticancer drugs are not well tolerated in all older patients with cancer. Lancet Healthy Longev. 2020;1:e43–7.

    Article  PubMed  Google Scholar 

  19. Hotta K, Ueoka H, Kiura K, Tabata M, Tanimoto M. An overview of 48 elderly-specific clinical trials of systemic chemotherapy for advanced non-small cell lung cancer. Lung Cancer. 2004;46:61–76.

    Article  PubMed  Google Scholar 

  20. Tesarova P. Breast cancer in the elderly-Should it be treated differently? Rep Pract Oncol Radiother. 2012;18:26–33.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Sanabria A, Carvalho AL, Vartanian JG, Magrin J, Ikeda MK, Kowalski LP. Comorbidity is a prognostic factor in elderly patients with head and neck cancer. Ann Surg Oncol. 2007;14:1449–57.

    Article  PubMed  Google Scholar 

  22. Brain E, Girre V, Rollot F, Bonnetain F, Debled M, Lacroix M, et al. ASTER 70s: benefit of adjuvant chemotherapy for estrogen receptor-positive HER2-negative breast cancer in women over 70 according to genomic grade—A French GERICO/UCBG UNICANCER multicenter phase III trial. J Clin Orthod. 2012;30:TPS667–TPS667.

    Google Scholar 

  23. Hughes KS, Schnaper LA, Berry D, Cirrincione C, McCormick B, Shank B, et al. Lumpectomy plus tamoxifen with or without irradiation in women 70 years of age or older with early breast cancer. N Engl J Med. 2004;351:971–7.

    Article  CAS  PubMed  Google Scholar 

  24. Sotiriou C, Wirapati P, Loi S, Harris A, Fox S, Smeds J, et al. Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis. J Natl Cancer Inst. 2006;98:262–72.

    Article  CAS  PubMed  Google Scholar 

  25. van’t Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AAM, Mao M, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002;415:530–6.

    Article  Google Scholar 

  26. Paik S, Shak S, Tang G, Kim C, Baker J, Cronin M, et al. A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N Engl J Med. 2004;351:2817–26.

    Article  CAS  PubMed  Google Scholar 

  27. Tobin NP, Lundberg A, Lindström LS, Harrell JC, Foukakis T, Carlsson L, et al. PAM50 provides prognostic information when applied to the lymph node metastases of advanced breast cancer patients. Clin Cancer Res. 2017;23:7225–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Parker JS, Mullins M, Cheang MCU, Leung S, Voduc D, Vickery T, et al. Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol. 2009;27:1160–7.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Wennmalm K, Bergh J. A simple method for assigning genomic grade to individual breast tumours. BMC Cancer. 2011;11:306.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Raj-Kumar P-K, Liu J, Hooke JA, Kovatich AJ, Kvecher L, Shriver CD, et al. PCA-PAM50 improves consistency between breast cancer intrinsic and clinical subtyping reclassifying a subset of luminal A tumors as luminal B. Sci Rep. 2019;9:7956.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  31. Johansson A, Yu NY, Iftimi A, Tobin NP, van’t Veer L, Nordenskjöld B, et al. Clinical and molecular characteristics of estrogen receptor-positive ultralow risk breast cancer tumors identified by the 70-gene signature. Int J Cancer. 2022;150:2072–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Tobin NP, Harrell JC, Lövrot J, Egyhazi Brage S, Frostvik Stolt M, Carlsson L, et al. Molecular subtype and tumor characteristics of breast cancer metastases as assessed by gene expression significantly influence patient post-relapse survival. Ann Oncol. 2015;26:81–8.

    Article  CAS  PubMed  Google Scholar 

  33. Lundberg A, Lindström LS, Harrell JC, Falato C, Carlson JW, Wright PK, et al. Gene expression signatures and immunohistochemical subtypes add prognostic value to each other in breast cancer cohorts. Clin Cancer Res. 2017;23:7512–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Sparano JA, Gray RJ, Makower DF, Pritchard KI, Albain KS, Hayes DF, et al. Adjuvant chemotherapy guided by a 21-gene expression assay in breast cancer. N Engl J Med. 2018;379:111–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Kjällquist U, Acs B, Margolin S, Karlsson E, Kessler LE, Garcia Hernandez S, et al. Real world evaluation of the prosigna/PAM50 test in a node-negative postmenopausal Swedish population: a multicenter study. Cancers. 2022. https://doi.org/10.3390/cancers14112615.

    Article  PubMed  PubMed Central  Google Scholar 

  36. McShane LM, Altman DG, Sauerbrei W, Taube SE, Gion M, Clark GM. REporting recommendations for tumour MARKer prognostic studies (REMARK). Br J Cancer. 2005;93:387–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Rueda OM, Sammut S-J, Seoane JA, Chin S-F, Caswell-Jin JL, Callari M, et al. Dynamics of breast-cancer relapse reveal late-recurring ER-positive genomic subgroups. Nature. 2019;567:399–404.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  38. Buus R, Sestak I, Barron S, Loughman T, Fender B, Ruiz CL, et al. Validation of the OncoMasTR risk score in estrogen receptor-positive/HER2-negative patients: a TransATAC study. Clin Cancer Res. 2020;26:623–31.

    Article  CAS  PubMed  Google Scholar 

  39. Sestak I, Buus R, Cuzick J, Dubsky P, Kronenwett R, Denkert C, et al. Comparison of the performance of 6 prognostic signatures for estrogen receptor-positive breast cancer: a secondary analysis of a randomized clinical trial. JAMA Oncol. 2018;4:545–53.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Brain E, Viansone AA, Bourbouloux E, Rigal O, Ferrero J-M, Kirscher S, et al. Final results from a phase III randomized clinical trial of adjuvant endocrine therapy ± chemotherapy in women ≥ 70 years old with ER HER2- breast cancer and a high genomic grade index: the Unicancer ASTER 70s trial. J Clin Oncol. 2022. https://doi.org/10.1200/jco.2022.40.16_suppl.500.

    Article  Google Scholar 

  41. Noordhoek I, Bastiaannet E, de Glas NA, Scheepens J, Esserman LJ, Wesseling J, et al. Validation of the 70-gene signature test (MammaPrint) to identify patients with breast cancer aged ≥ 70 years with ultralow risk of distant recurrence: a population-based cohort study. J Geriatr Oncol. 2022. https://doi.org/10.1016/j.jgo.2022.07.006.

    Article  PubMed  Google Scholar 

  42. Mansani F, Celinski V, Freitas-Junior R. De-escalation of chemotherapy in elderly women using a 70-gene platform: comparison of the MINDACT study with a real-world study in the Brazilian population (AGEMA-BRA). J Clin Orthod. 2022;40:e12568–e12568.

    Google Scholar 

  43. Iles K, Roberson ML, Spanheimer P, Gallagher K, Ollila DW, Strassle PD, et al. The impact of age and nodal status on variations in oncotype DX testing and adjuvant treatment. npj Breast Cancer. 2022;8:1–8.

    Article  Google Scholar 

  44. Barni S, Cognetti F, Petrelli F. Is the oncotype DX test useful in elderly breast cancer patients: a subgroup analysis of real-life Italian PONDx study. Breast Cancer Res Treat. 2022;191:477–80.

    Article  CAS  PubMed  Google Scholar 

  45. Brain EGC, Mertens C, Girre V, Rousseau F, Blot E, Abadie S, et al. Impact of liposomal doxorubicin-based adjuvant chemotherapy on autonomy in women over 70 with hormone-receptor-negative breast carcinoma: a French Geriatric Oncology Group (GERICO) phase II multicentre trial. Crit Rev Oncol Hematol. 2011;80:160–70.

    Article  PubMed  Google Scholar 

  46. Freyer G, Campone M, Peron J, Facchini T, Terret C, Berdah J-F, et al. Adjuvant docetaxel/cyclophosphamide in breast cancer patients over the age of 70: results of an observational study. Crit Rev Oncol Hematol. 2011;80:466–73.

    Article  PubMed  Google Scholar 

  47. Takabatake D, Taira N, Hara F, Sien T, Kiyoto S, Takashima S, et al. Feasibility study of docetaxel with cyclophosphamide as adjuvant chemotherapy for Japanese Breast Cancer patients. Jpn J Clin Oncol. 2009;39:478–83.

    Article  PubMed  Google Scholar 

  48. Yu J, Wu J, Huang O, He J, Li Z, Chen W, et al. Do 21-gene recurrence score influence chemotherapy decisions in T1bN0 breast cancer patients? Front Oncol. 2020;10:708.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Jézéquel P, Sharif Z, Lasla H, Gouraud W, Guérin-Charbonnel C, Campion L, et al. Gene-expression signature functional annotation of breast cancer tumours in function of age. BMC Med Genom. 2015;8:80.

    Article  Google Scholar 

  50. de Kruijf EM, Bastiaannet E, Rubertá F, de Craen AJM, Kuppen PJK, Smit VTHBM, et al. Comparison of frequencies and prognostic effect of molecular subtypes between young and elderly breast cancer patients. Mol Oncol. 2014;8:1014–25.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Jenkins EO, Deal AM, Anders CK, Prat A, Perou CM, Carey LA, et al. Age-specific changes in intrinsic breast cancer subtypes: a focus on older women. Oncologist. 2014;19:1076–83.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Gosain R, Pollock Y, Jain D. Age-related disparity: breast cancer in the elderly. Curr Oncol Rep. 2016;18:69.

    Article  PubMed  Google Scholar 

  53. Berben L, Floris G, Wildiers H, Hatse S. Cancer and aging: two tightly interconnected biological processes. Cancers. 2021. https://doi.org/10.3390/cancers13061400.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Diebel LWM, Rockwood K. Determination of biological age: geriatric assessment vs biological biomarkers. Curr Oncol Rep. 2021;23:104.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Hong C-C, Ambrosone CB, Goodwin PJ. Comorbidities and their management: potential impact on breast cancer outcomes. Adv Exp Med Biol. 2015;862:155–75.

    Article  PubMed  Google Scholar 

Download references

Funding

Open access funding provided by Karolinska Institute. This work was supported by the Iris, Stig och Gerry Castenbäcks Stiftelse for cancer research (N.P.T.); the King Gustaf V Jubilee Foundation (N.P.T.); the Stockholm Cancer Society (Cancerföreningen i Stockholm to L.S.L.); the Swedish Cancer Society (Cancerfonden, N.P.T. grant Number: 200802; L.S.L. Grant Number: 222081, 220552SIA, A.M. Grant Number 21 0277 JCIA 01 H); the Swedish Research Council (Vetenskapsrådet, Grant Number 2020-02466 to L.S.L); the Swedish Research Council for Health, Working life and Welfare, (FORTE, Grant Number 2019-00477 to L.S.L.); ALF medicine (Grant Number FoUI-974882 to L.S.L.) and the Gösta Milton Donation Fund (Stiftelsen Gösta Miltons donationsfond, to L.S.L).

Author information

Authors and Affiliations

Authors

Contributions

MCA and NPT contributed to the study concept and design. MCA and NPT contributed to the acquisition and analyses of data. All authors interpreted the data, drafted the manuscript, read, and approved the final version.

Corresponding author

Correspondence to Nicholas P. Tobin.

Ethics declarations

Ethics approval and consent to participate

Publicly available data—Not applicable.

Competing interests

N.P.T. has performed consultations to Pfizer. A.M. has performed consultations to Veracyte (no financial or any other compensation) and Roche (no financial or any other compensation). T.F. consultancy to Astra Zeneca, Gilead, Affibody, Pfizer, Novartis, Veracyte, Exact Sciences and Roche; honoraria from UpToDate; research funding to institution from Pfizer, Astra Zeneca, Novartis.

Additional information

Publisher's Note

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

Supplementary Information

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

Castresana-Aguirre, M., Johansson, A., Matikas, A. et al. Clinically relevant gene signatures provide independent prognostic information in older breast cancer patients. Breast Cancer Res 26, 38 (2024). https://doi.org/10.1186/s13058-024-01797-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13058-024-01797-7