Intended for healthcare professionals

CCBY Open access
Research

Breast cancer mortality in 500 000 women with early invasive breast cancer diagnosed in England, 1993-2015: population based observational cohort study

BMJ 2023; 381 doi: https://doi.org/10.1136/bmj-2022-074684 (Published 13 June 2023) Cite this as: BMJ 2023;381:e074684

Linked Opinion

Risk of breast cancer death after a diagnosis of early invasive breast cancer

  1. Carolyn Taylor, professor of oncology and honorary clinical oncologist1 2,
  2. Paul McGale, statistician1,
  3. Jake Probert, statistician1,
  4. John Broggio, cancer analytical lead3,
  5. Jackie Charman, senior cancer analyst3,
  6. Sarah C Darby, professor of medical statistics1,
  7. Amanda J Kerr, systematic reviewer1,
  8. Timothy Whelan, radiation oncologist4,
  9. David J Cutter, senior clinical research fellow and clinical oncologist1 2,
  10. Gurdeep Mannu, lecturer in general surgery1,
  11. David Dodwell, senior clinical research fellow and clinical oncologist1 2
  1. 1Nuffield Department of Population Health, University of Oxford, Oxford, UK
  2. 2Oxford University Hospitals, Oxford, UK
  3. 3National Disease Registration Service (NDRS), NHS England, Birmingham, UK
  4. 4Department of Oncology, McMaster University and Juravinski Cancer Centre, Hamilton, ON Canada
  1. Correspondence to: C Taylor carolyn.taylor{at}ndph.ox.ac.uk
  • Accepted 26 April 2023

Abstract

Objectives To describe long term breast cancer mortality among women with a diagnosis of breast cancer in the past and estimate absolute breast cancer mortality risks for groups of patients with a recent diagnosis.

Design Population based observational cohort study.

Setting Routinely collected data from the National Cancer Registration and Analysis Service.

Participants All 512 447 women registered with early invasive breast cancer (involving only breast and possibly axillary nodes) in England during January 1993 to December 2015, with follow-up to December 2020.

Main outcome measures Annual breast cancer mortality rates and cumulative risks by time since diagnosis, calendar period of diagnosis, and nine characteristics of patients and tumours.

Results For women with a diagnosis made within each of the calendar periods 1993-99, 2000-04, 2005-09, and 2010-15, the crude annual breast cancer mortality rate was highest during the five years after diagnosis and then declined. For any given time since diagnosis, crude annual breast cancer mortality rates and risks decreased with increasing calendar period. Crude five year breast cancer mortality risk was 14.4% (95% confidence interval 14.2% to 14.6%) for women with a diagnosis made during 1993-99 and 4.9% (4.8% to 5.0%) for women with a diagnosis made during 2010-15. Adjusted annual breast cancer mortality rates also decreased with increasing calendar period in nearly every patient group, by a factor of about three in oestrogen receptor positive disease and about two in oestrogen receptor negative disease. Considering just the women with a diagnosis made during 2010-15, cumulative five year breast cancer mortality risk varied substantially between women with different characteristics: it was <3% for 62.8% (96 085/153 006) of women but ≥20% for 4.6% (6962/153 006) of women.

Conclusions These five year breast cancer mortality risks for patients with a recent diagnosis may be used to estimate breast cancer mortality risks for patients today. The prognosis for women with early invasive breast cancer has improved substantially since the 1990s. Most can expect to become long term cancer survivors, although for a few the risk remains appreciable.

Introduction

Worldwide, more than 2 million patients receive a diagnosis of invasive breast cancer each year.1 For most, it is their first cancer. Most have early stage disease and receive surgery as their first treatment. Outcomes following treatment for early invasive breast cancer differ substantially across different countries and across patients with different characteristics within a country. Patients with breast cancer and the clinicians who treat them need estimates of their likely prognosis to inform treatment decisions, follow-up, and prediction of event rates for groups of patients in clinical trials. These estimates require large scale, population based studies that consider the effects of multiple patient related and tumour related factors on breast cancer mortality. Survival has improved, and detailed breast cancer mortality estimates based on patients treated recently are not available. These are needed to enable clinicians to estimate prognosis for patients treated today by using characteristics such as age, tumour size, nodal status, tumour grade, receptor status, and screening status.

Most population based studies in women with early invasive breast cancer either consider just a few specific determinants of breast cancer mortality or consider selected groups of women with specific characteristics.234567891011 An additional limitation is that few of these studies have considered the extent to which breast cancer mortality has changed over the past few decades for women with early invasive breast cancer or the way in which it varies with time since diagnosis.

We conducted the first observational cohort study showing long term outcomes in all women registered in England with early invasive breast cancer as their first cancer and who were treated initially with surgery. We describe annual breast cancer mortality rates and cumulative breast cancer mortality risks in all women with a diagnosis made during 1993-2015, and the extent to which they vary with calendar period of diagnosis, time since diagnosis, whether the cancer was screen detected, and characteristics of patients and tumours. We then focus on women with a diagnosis made during 2010-15 to inform patients and clinicians about the likely absolute mortality risks for individual patients treated for breast cancer today, taking into account their tumour characteristics, including whether the cancer was screen detected. We did not aim to assess the effects of treatment, which would need data from randomised trials.12 Neither did we aim to quantify the effects of screening on changes in mortality.

Methods

Study population

The National Cancer Registration and Analysis Service (NCRAS), together with its preceding organisations, has registered all patients with a diagnosis of cancer in England for many decades.13 These registrations are routinely linked with other information at the level of the patient. For this study, all women registered in England during January 1993 to December 2015 with breast cancer as their first invasive cancer were identified and their data were de-personalised before release for analysis.

Early invasive breast cancer was defined as disease detected only in the breast or, in the case of women with node positive disease, the breast and axillary lymph nodes with no evidence of metastatic disease.14 Surgery was either breast conserving surgery or mastectomy.

Data included pathological staging (tumour size and number of positive axillary nodes in the surgical specimen), grade (low, medium, or high), and oestrogen receptor status (with oestrogen receptor positivity defined as >10% staining or Allred score 3-8). For women with a diagnosis made during 2010-15, human epidermal growth factor receptor 2 (HER2) status was also included. These variables are routinely used in the clinic to estimate prognosis,1516 to define groups of patients in clinical guidelines, to inform treatment decisions, and to assess eligibility for clinical trials. Other variables were calendar period of diagnosis, age at diagnosis, screening status, laterality, index of multiple deprivation, region of residence, and, where applicable, dates of emigration and death, together with causes of death, up to 31 December 2020.

Data were collated and checked before analysis; see supplementary text S1 and figure S1 for details. We received data on a total of 783 980 women (figure S2). We excluded women if they were aged <18 or ≥90 years at diagnosis (n=24 157), their registrations were based on a death certificate only or they had less than three months of follow-up (n=33 086), or their histology was not invasive breast cancer (n=10 798). We also excluded the 22 577 women with diagnoses of two simultaneous cancers, as patterns of breast cancer mortality among these women may differ from those of women with just one cancer. These women comprised those with a second primary non-breast cancer within three months of registration and women with bilateral breast cancer or a contralateral breast cancer within three months of registration, who would each have two sets of breast tumour related factors. Analysing data for them would require a different method from that appropriate for women with unilateral breast cancer and just one set of tumour related factors. In addition, we excluded 138 911 women because they had likely metastatic disease. A further 42 004 were recorded as receiving neoadjuvant therapy (7.6% (42 004/554 451) of women overall; 8.1% (19 996/248 176) of women with a diagnosis made in 1993-2004 and 7.2% (22 008/306 275) with a diagnosis made in 2005-15). We excluded them as the characteristics of their tumours were recorded before the neoadjuvant therapy was given and are therefore based on clinical staging. This is in contrast to other women, for whom pathological staging was recorded.

Statistical analysis

Women were included in the study from three months after diagnosis and removed on the earliest of date of death, emigration, 95th birthday, or end of follow-up. We considered calendar period in five year categories, except that we included women with diagnoses made during 1993-94 with the 1995-99 period and included women with diagnoses made during 2015 with the 2010-14 period. Within each calendar period category, we wanted the maximum possible length of follow-up to be the same for all women. Therefore, unless a woman had previously died, emigrated, or reached her 95th birthday, we removed her from the study at 21 years from diagnosis if she had a diagnosis made during 1993-99, 16 years if she had a diagnosis made during 2000-04, 11 years if she had a diagnosis made during 2005-09, and five years if she had a diagnosis made during 2010-15 (text S2).

Information was missing for some women on tumour size, number of positive nodes, tumour grade, oestrogen receptor status, laterality, and (applicable to the period 2010-15) HER2 status. Therefore, to avoid the loss of precision and possible bias that would be caused by omitting these women from the analysis or by creating separate categories for missing values, we used multiple imputation of missing values in analyses involving these characteristics (see text S2, table S1, and table S2 for details of missing data).

We did the main analyses in four steps. In the first step, we estimated crude annual breast cancer mortality rates by time since diagnosis assuming a Poisson distribution. We then smoothed these annual rates (see figure S3 for a comparison of raw and smoothed rates) and estimated cumulative risks by time since diagnosis from the smoothed rates. We made these calculations separately for categories of calendar period of diagnosis, screening status, and oestrogen receptor status and also for combinations of these three characteristics.

In the second step, which we did separately for women with oestrogen receptor positive and oestrogen receptor negative disease, we calculated adjusted annual mortality rates according to each available characteristic by Poisson regression, with time since diagnosis and all other available characteristics as categorical variables using categories defined as shown in table 1. The third step examined the extent to which both absolute and proportional changes in the mortality rate with increasing calendar period of diagnosis varied across the different characteristics. We found substantial variability on both absolute and proportional scales. This meant that the data could not be summarised in any simple quantitative model. Therefore, in the fourth step, we did analyses of just the women with diagnoses made during 2010-15 to describe breast cancer mortality rates for groups of patients with a recent diagnosis. We calculated adjusted annual mortality rates for each available characteristic by Poisson regression, as in the second step above. This was followed by tests for pairwise interactions between the characteristics. Lastly, we calculated cumulative five year risks separately for women with different combinations of the characteristics HER2 status, age at diagnosis, screening status, number of positive nodes, and tumour size and grade.

Table 1

Characteristics of 512 447 women given a diagnosis of early breast cancer in England during 1993-2015.

View this table:

We did analyses for all cause mortality along similar lines to those for breast cancer mortality. We used Stata version 17 for all calculations. Further details are in text S2 and tables S1 and S2.

Patient and public involvement

Two patients representing the organisation Independent Cancer Patients’ Voice were involved as research partners. They helped to develop the research questions. They identified the need for up-to-date information on outcomes after breast cancer diagnosis to inform prognosis for current and future patients. They advised on which analyses were important to patients with breast cancer and provided input to the interpretation. They reviewed and commented on the main findings in the manuscript. Their input was provided in face-to-face meetings, in teleconferences, and via email. They have agreed to help with dissemination of the findings.

Results

Characteristics of study population

A total of 512 447 women with early invasive breast cancer were included in the study. Breast surgery was either breast conserving surgery (60%; 307 714 women) or mastectomy (40%; 204 733 women). For axillary surgery, 154 583 (30%) had axillary dissection only, 218 313 (43%) had sentinel lymph node biopsy only, 79 559 (16%) had both, and 59 992 (12%) had no record of either. The number of women who received a diagnosis of early breast cancer each year increased over the study period from an average of 16 193 during 1993-99 to 26 056 during 2010-15 (table 1). Before 2005 only women aged 50-64 years were eligible for screening and the overall percentage of women whose cancers were screen detected was 13% (14 981/113 354) in 1993-99 and 18% (20 190/114 826) in 2000-04. These percentages increased substantially when screening was also offered to women aged 65-70 years, to 32% (41 272/127 929) during 2005-09 and 33% (51 797/156 338) during 2010-15. Overall, nearly half the cancers among women in age groups eligible for screening were screen detected.

The percentage of women with node negative disease increased from 47% (53 773/113 354) in 1993-99 to 65% (101 898/156 338) in 2010-15, and the percentage with tumours categorised as oestrogen receptor positive increased from 77% (86 965/113 354) in 1993-99 to 86% (134 398/156 338) in 2010-15. By contrast, the percentage of women with a diagnosis of low grade disease decreased from 22% (25 257/113 354) in 1993-99 to 17% (26 738/156 338) in 2010-15.

The geographical spread reflected the population of England, with more than 30 000 women from each region. Among the 512 447 women, duration of follow-up from time of diagnosis was five years or less for 211 946 (41%) women, five to 10 years for 46 840 (9%) women, 10-15 years for 123 040 (24%) women, 15-20 years for 80 051 (16%) women, and more than 20 years for 50 570 (10%) women (table 1). By end of follow-up, 77 975 (15%) of the women in the study had died from breast cancer and 155 895 (30%) from any cause.

All women: crude breast cancer mortality

To assess the pattern of breast cancer mortality by time since diagnosis in each calendar period, we analysed data on women with a diagnosis made in each calendar period separately. Regardless of calendar period of diagnosis, annual breast cancer mortality rates increased during the two years following diagnosis, peaked during the third year, and then declined (fig 1, top left). For any given time since diagnosis, cumulative risk of breast cancer mortality decreased progressively in successive calendar periods (fig 1, top right). Five year cumulative mortality risk was 14.4% (95% confidence interval 14.2% to 14.6%) for women with a diagnosis made during 1993-99 and reduced progressively with increasing calendar period to 4.9% (4.8% to 5.0%) for women with a diagnosis made during 2010-15. We observed corresponding decreases in mortality risk when we repeated the analyses including deaths from other causes as well as breast cancer (figures S23 and S24) or including all women with a diagnosis of breast cancer, even if they had metastatic disease at diagnosis or were recorded as receiving neoadjuvant therapy (figure S25).

Fig 1
Fig 1

Crude annual breast cancer mortality rates and cumulative breast cancer mortality risks in 512 447 women with early breast cancer by time since diagnosis according to calendar period of diagnosis: for all women (top panel), for women aged 50-64 years (who would all have been eligible for screening) according to whether their cancer was screen detected (middle panels), and according to oestrogen receptor (ER) status (bottom panel). Vertical lines are 95% confidence intervals. Numbers of women at risk, numbers of deaths from breast cancer, annual rates, and cumulative mortality risks are given in tables S4 to S6. Further analyses of women by ER status are in figures S4 and S5. Analyses of non-breast cancer mortality and all cause mortality are in figures S23 and S24. Analyses of women with metastatic as well as early breast cancer, and including women who received neoadjuvant therapy, are in figure S25

To assess whether the decreases in breast cancer mortality among women given a diagnosis of early invasive breast cancer resulted solely from increases in screening, we did separate analyses for women with a diagnosis made in different calendar periods and with different screening status for women aged 50-64 years at diagnosis, all of whom would have been eligible for screening throughout the period included in the study (fig 1, middle panels). For patients with a diagnosis of either screen detected or non-screen detected cancers, annual breast cancer mortality rates and cumulative breast cancer mortality risks showed similar patterns by calendar period of diagnosis to those for all women, but with lower values in women with screen detected compared with non-screen detected cancers.

To assess differences in breast cancer mortality according to oestrogen receptor status, we did separate analyses for women with oestrogen receptor positive and oestrogen receptor negative disease. Considering all calendar periods together, annual breast cancer mortality rates peaked during the third year after diagnosis in both groups but at a much higher rate in oestrogen receptor negative disease (fig 1, bottom panel). Beyond three years, annual breast cancer mortality rates decreased in both groups but more rapidly in oestrogen receptor negative than oestrogen receptor positive disease so that, beyond 10 years, the annual breast cancer mortality rate was higher in oestrogen receptor positive than oestrogen receptor negative disease. In both groups, deaths from breast cancer continued to occur until the end of follow-up, irrespective of screening status (figure S4).

All women: adjusted breast cancer mortality

We analysed data on all 512 447 women with a diagnosis made during 1993-2015 to assess the independent effects of calendar period of diagnosis and eight other characteristics of patients and tumours on breast cancer mortality in oestrogen receptor positive and oestrogen receptor negative disease. In these analyses, we calculated annual breast cancer mortality rates separately for each characteristic by oestrogen receptor status after adjusting for the other eight characteristics and for time since diagnosis.

The adjusted annual breast cancer mortality rate declined steadily with calendar period of diagnosis in both oestrogen receptor positive and oestrogen receptor negative disease, but rates were higher in oestrogen receptor negative disease (fig 2, top left). For every other characteristic, the adjusted annual breast cancer mortality rates also followed a similar pattern in oestrogen receptor positive and oestrogen receptor negative disease but with higher rates in oestrogen receptor negative disease (fig 2; figures S6 and S7). They were lowest for women in their 40s at diagnosis and highest for women in their 80s at diagnosis, and they were lower in women with screen detected cancer than in women with non-screen detected cancer or in women who were not eligible for screening. Tumour size, number of nodes, and tumour grade all had substantial independent effects, whereas the effects of cancer laterality, deprivation, and region of residence were smaller (figure S6). See figures S26-S28 for analyses of all cause mortality.

Fig 2
Fig 2

Adjusted annual breast cancer mortality rates in 512 447 women with early breast cancer with oestrogen receptor (ER) positive or ER negative disease, by various characteristics. For each characteristic, rates are adjusted for all other characteristics shown and also for time since diagnosis, breast cancer laterality, index of multiple deprivation, and region of residence. Vertical lines are 95% confidence intervals. Results for breast cancer laterality, index of multiple deprivation, and region of residence are in figure S6. Further details are in figure S7, and separate rates for years 0-5, 5-15, and ≥15 since diagnosis are in figures S8-S10. Analyses of all cause mortality are in figures S26-S28

The variation in the adjusted breast cancer mortality rates for each characteristic considered was greatest in the first five years after diagnosis, but the effects of age at diagnosis, screening status, tumour size, number of positive nodes, and tumour grade were all still strong five to 15 years after diagnosis, and the effects of tumour size, number of nodes, and, in oestrogen receptor positive disease, age at diagnosis were still clearly present at ≥15 years (figures S8-S10).

To examine whether the decreases in breast cancer mortality with more recent calendar period (fig 2, top left) took place to the same extent among all groups of women, we calculated adjusted annual breast cancer mortality rates by oestrogen receptor status and calendar period of diagnosis for women with different characteristics. The breast cancer mortality rate decreased substantially between 1993-99 and 2010-15 for all ages at diagnosis in women with oestrogen receptor positive disease and for most ages in women with oestrogen receptor negative disease, apart from those aged 80-89 at diagnosis, for whom we saw hardly any decrease (fig 3; figures S11 and S13). In younger women, breast cancer was the main cause of death so reductions in breast cancer mortality were similar to reductions in all cause mortality according to calendar period. In older women, all cause mortality was much higher than breast cancer mortality. Mortality reductions according to calendar period were greater for all cause mortality than for breast cancer mortality (figure S30).

Fig 3
Fig 3

Adjusted annual breast cancer mortality rates in 512 447 women with early breast cancer with oestrogen receptor (ER) positive or ER negative disease by calendar period of diagnosis, according to age at diagnosis. Rates are adjusted for whether cancer was screen detected, tumour size, number of positive nodes, tumour grade, time since diagnosis, breast cancer laterality, index of multiple deprivation, and region of residence. Vertical lines are 95% confidence intervals. Further details are in figure S13. Results including only five years of follow-up (so that all age groups have the same length of follow-up) are in figure S11. Results for breast cancer laterality, index of multiple deprivation and region of residence are in figure S12. Results for all cause mortality are in figures S29 and S30

For screening status, tumour size, number of positive nodes, tumour grade, breast cancer laterality, deprivation, and region of residence, the breast cancer mortality rate decreased substantially between 1993-99 and 2010-15 for every category in both oestrogen receptor positive and oestrogen receptor negative disease (fig 4; fig 5; fig 6; fig 7; figures S12 and S14-20). See figures S29 and S30 for analyses of all cause mortality.

Fig 4
Fig 4

Adjusted annual breast cancer mortality rates in 512 447 women with early breast cancer with oestrogen receptor (ER) positive or ER negative disease by calendar period of diagnosis, according to whether their cancer was screen detected. Rates are adjusted for age at diagnosis, tumour size, number of positive nodes, tumour grade, time since diagnosis, breast cancer laterality, index of multiple deprivation, and region of residence. Vertical lines are 95% confidence intervals. Further details are in figure S14

Fig 5
Fig 5

Adjusted annual breast cancer mortality rates in 512 447 women with early breast cancer with oestrogen receptor (ER) positive or ER negative disease by calendar period of diagnosis, according to tumour size. Rates are adjusted for age at diagnosis, whether cancer was screen detected, number of positive nodes, tumour grade, time since diagnosis, breast cancer laterality, index of multiple deprivation, and region of residence. Vertical lines are 95% confidence intervals. Further details are in figure S15

Fig 6
Fig 6

Adjusted annual breast cancer mortality rates in 512 447 women with early breast cancer with oestrogen receptor (ER) positive or ER negative disease by calendar period of diagnosis, according to number of positive nodes. Rates are adjusted for age at diagnosis, whether cancer was screen detected, tumour size, tumour grade, time since diagnosis, breast cancer laterality, index of multiple deprivation, and region of residence. Vertical lines are 95% confidence intervals. Further details are in figure S16

Fig 7
Fig 7

Adjusted annual breast cancer mortality rates in 512 447 women with early breast cancer with oestrogen receptor (ER) positive or ER negative disease by calendar period of diagnosis, according to tumour grade. Rates are adjusted for age at diagnosis, whether cancer was screen detected, tumour size, number of positive nodes, time since diagnosis, breast cancer laterality, index of multiple deprivation, and region of residence. Vertical lines are 95% confidence intervals. Further details are in figure S17

Although a decrease in breast cancer mortality rate in more recent calendar periods was apparent for all groups of women, apart from those aged 80-89 with oestrogen receptor negative disease, the absolute magnitude of the decrease varied substantially. For example, in oestrogen receptor positive disease, the adjusted annual rate decreased by 2.99 percentage points (from 5.80% (5.24% to 6.41%) in 1993-99 to 2.81% (2.57% to 3.08%) in 2010-15) among women with ≥10 positive nodes, but by only 0.75 percentage points (from 1.11% (1.05% to 1.19%) in 1993-99 to 0.36% (0.34% to 0.38%) in 2010-15) for women with no positive nodes (figure S16). Similarly, in oestrogen receptor negative disease, the adjusted rate decreased by 3.30 percentage points (from 8.20% (7.21% to 9.32%) in 1993-99 to 4.90% (4.34% to 5.54%) in 2010-15) among women with ≥10 positive nodes, but by only 0.68 percentage points (from 1.56% (1.35% to 1.81%) in 1993-99 to 0.88% (0.81% to 0.94%) in 2010-15) among women with no positive nodes. We also saw substantial variations in the absolute magnitude of the decrease across the different levels of many other characteristics (figures S13-15, S17-S20). The magnitude of the decreases with calendar period of diagnosis also varied between different levels of many characteristics when considered on a proportional scale (table S7). Therefore, as women who had a diagnosis made in the most recent calendar period are most relevant to women with diagnoses of breast cancer today, further analysis considered these women in greater detail.

Women with diagnosis made 2010-15: breast cancer mortality

Among the 156 338 women with a diagnosis made during 2010-15, the adjusted annual breast cancer mortality rates varied substantially across categories of age, screening status, tumour size, number of positive nodes, and tumour grade in both oestrogen receptor positive and oestrogen receptor negative disease, whereas for HER2 status, which was available for women with diagnoses made in this time period, the adjusted annual mortality rate was lower in HER2 positive than HER2 negative disease in both oestrogen receptor positive and oestrogen receptor negative disease (figure S21). In addition, we found evidence that some of these factors did not combine in a multiplicative fashion in oestrogen receptor positive disease (table S8). Therefore, we calculated separate estimates of cumulative five year breast cancer mortality risk for all the 576 patient groups with different combinations of these factors (fig 8; table S9).

Fig 8
Fig 8

Cumulative five year breast cancer mortality risks in 156 338 women with early breast cancer diagnosed during 2010-15 by categories of tumour grade, tumour size and number of positive nodes in women with oestrogen receptor (ER) positive or ER negative disease. Figures are split by human epidermal growth factor receptor 2 (HER2) status, age and screening status. Vertical lines are 95% confidence intervals. Points are plotted for groups of women with data on ≥40 women and include 153 006/156 338 (97.9%) of women. Further details, including values of plotted points, are given in tables S9 and S10. Results on a square root scale are in figure S22. Results for all cause mortality are in figure S31

The largest such group was the 15 533 women aged 50-70 years at diagnosis with screen detected cancer that was HER2 negative, oestrogen receptor positive, medium grade, size 1-20 mm, and node negative (fig 8, second row, left side; table S9b). Women in this group had an average cumulative five year breast cancer mortality risk of 0.5% (0.4% to 0.7%). We found, however, substantial variability in the five year cumulative risk between the different combinations of characteristics. It took values <1% for 34.2% (52 257/153 006) of women, <3% for 62.8% (96 085/153 006), <5% for 73.7% (112 814/153 006), and <10% for 88.0% (134 655/153 006), whereas it was ≥20% for 4.6% (6962/153 006), ≥40% for 0.7% (1085/153 006), and ≥60% for 0.2% (238/153 006) of women (table S10).

Discussion

We have identified all 512 447 women registered with early invasive breast cancer in England during 1993 to 2015 who were treated initially with surgery and followed them to 31 December 2020. We started by documenting associations between the breast cancer mortality rate and nine individual characteristics, providing insight into the prognostic value of each. We showed that the breast cancer mortality rate has decreased with calendar period of diagnosis over the period in our study. Although decreases occurred across nearly all groups of patients, the magnitude of the decrease in the mortality rate varied substantially between women with different characteristics, both on an absolute scale and on a proportional scale. This indicates that the associations between the variables we have studied and breast cancer mortality have been changing over time. Therefore, although the long term trends in women whose cancer was diagnosed during 1993-2009 are of some relevance for long term predictions, it is the mortality of women with breast cancer diagnosed in the most recent period that is the most relevant to women receiving a diagnosis today. Breast cancer mortality is highest during the five years following diagnosis and, although deaths from breast cancer will continue to occur beyond this, the risk during each subsequent five year period is likely to be lower than during the first five years. For these reasons, we presented the five year breast cancer mortality risks for women with a diagnosis made during 2010-15 with different combinations of age, screening status, oestrogen receptor status, number of positive nodes, and tumour size and grade. Focusing on 2010-15 also enabled us to consider HER2 status, which was not available for earlier calendar periods. Clinicians can use this information to estimate breast cancer mortality risks for patients today.

Strengths and limitations of study

This is the first large, population based, long term study of women with early invasive breast cancer describing associations between breast cancer mortality and multiple characteristics of patients and tumours. Women with bilateral, previous, or metastatic cancer were excluded, as were women who did not receive initial surgery. Our results therefore derive from the large majority of women with a diagnosis of breast cancer made during the study period. Most existing studies consider only a few characteristics of patients and tumours, whereas we included all those available from NCRAS, including screening status.

Before analysis, the data underwent rigorous checking by statisticians with input from a surgeon, a pathologist, and three oncologists. We queried inconsistent or missing values with NCRAS staff and corrected them wherever possible. For any remaining missing values, we used multiple imputation. This method enables all women to be included in the analysis, taking account of any correlations between the missing variable and the variables that are known. It also allows for the uncertainty arising from the missing values in confidence intervals and significance tests.

Our study provides an accurate picture of breast cancer mortality in a complete population of women for up to 20 years. During this period, the biology of the disease may have changed owing to hormonal changes arising from obesity, the use of hormone replacement therapy, and reproductive factors.17 New systemic therapies, including aromatase inhibitors and taxanes, have become established in clinical use,18 and surgery and radiotherapy have become better targeted. These changes in treatment undoubtedly account for some of the decrease in breast cancer mortality seen over the study period. The use of trastuzumab in HER2 positive disease has also increased since the mid-2000s, and our results indicate that women with HER2 positive disease had lower breast cancer mortality than those with HER2 negative disease after adjustment for other factors, probably owing to the use of trastuzumab. Other factors, such as increased breast awareness, screening, and more sensitive breast imaging techniques, are also likely to have contributed to reducing breast cancer mortality but may in some cases have led just to earlier diagnosis and apparently longer survival without, in fact, changing the clinical course of the disease. Notably, however, the improvements in breast cancer mortality seen in women with screen detected cancers were paralleled by improvements in those whose cancers were not screen detected. Therefore, increases in screening cannot solely explain the decreases in breast cancer mortality that we observed. To account for the possible effect of stage migration, we also repeated key analyses in all women with unilateral breast cancer, including those with metastatic disease or who received neoadjuvant treatments. In these analyses, breast cancer mortality still decreased according to calendar period of diagnosis. Therefore, stage migration is not solely responsible for the observed decreases in breast cancer mortality. Neither can the decreases be attributed only to changes in tumour size, number of positive nodes, or tumour grade, as breast cancer mortality still decreased with calendar period of diagnosis, even after adjustment for these factors. Beyond these considerations, this observational study cannot determine the specific causes of these reductions in mortality. However, the main aim of our study was not to quantify the causal role played by different factors in the observed decreases in the breast cancer mortality rate. Instead, we provide information for clinicians to use when estimating breast cancer mortality risks for patients today, taking into account the characteristics of their tumour, including whether the cancer was screen detected. Our study provides, for the first time, estimates of prognosis based on all patients recently receiving a diagnosis of breast cancer in the whole of England.

A limitation of our study was that data on cancer recurrence were not available. Also, we were unable to include women who received neoadjuvant therapy, as pathological tumour size and number of positive axillary nodes at diagnosis were unavailable for them, and we could not include women simultaneously with a diagnosis of more than one cancer in these analyses.

Comparison with other studies

In recent years, the age standardised annual incidence rate of female breast cancer in the whole UK population rose from 137 per 100 000 in 1993-95 to 169 in 2016-18, whereas the corresponding breast cancer mortality rate in the whole UK population decreased steadily from 53 per 100 000 in 1993-95 to 34 in 2016-18.19 Therefore, the finding of a decrease in breast cancer specific mortality in our study of women with a diagnosis of breast cancer is not unexpected. However, national summary statistics provide little insight into the nature of the decrease—for example, whether decreases have occurred for women with early invasive breast cancer as well as those with metastatic disease or whether they have occurred in oestrogen receptor negative as well as oestrogen receptor positive disease. A literature search identified 10 studies published since 2000, each including at least 5000 women, that examined breast cancer mortality in patients with early invasive breast cancer in population based cancer registries or randomised trials (table S11).234567891011 Most studies focused on just a few specific factors, and the extent to which adjustment was made for other factors was variable. Up-to-date population based estimates of breast cancer mortality for groups of patients with individual tumour characteristics and with a recent diagnosis were unavailable. Several studies reported lower breast cancer mortality for women with more recent diagnoses,34689 and individual studies reported mortality patterns similar to those we have reported for oestrogen receptor positive versus oestrogen receptor negative disease,510 age at diagnosis,579 tumour size,211 number of positive nodes,711 and tumour grade.5711 One study showed a smaller increase in survival according to calendar period for women aged over 70 years at diagnosis than for younger women.9 In our data, the lack of any decrease in mortality for women aged 80-89 years with oestrogen receptor negative disease may be because these women do not usually receive adjuvant systematic therapy or radiotherapy,20 so any improvements in these treatments would not have affected mortality in this patient group. Patients aged less than 40 years at diagnosis had higher breast cancer mortality risk than patients in their 40s at diagnosis, after adjustment for patient and tumour related factors.79 This may be because breast cancers in younger women are intrinsically more aggressive than those in older women. Several studies have built models that focus on estimating the effects of treatment and so are not comparable with our study.2122

Implications of findings

Our findings illustrate the substantial improvement in prognosis for women with a diagnosis of early invasive breast cancer that has been made since the 1990s. More than six in 10 of the women given the diagnosis in England during 2010-15 had a risk of dying from breast cancer during the following five years of 3% or less, and nearly nine in 10 had a risk below 10%. Our findings have several implications.

Firstly, these breast cancer mortality risks inform patients today about their likely prognosis. For example, the estimated five year breast cancer mortality risk for a woman aged 60 at diagnosis, with a screen detected tumour, size <20 mm, low grade, oestrogen receptor positive, HER2 negative, and node negative would be 0.2% (table S9b). This is likely to provide her with reassurance about her prognosis. Such estimates are also relevant to decisions about treatment. Most women with breast cancer are cured by surgery, but undetected microscopic disease remains after surgery in some women, and adjuvant therapies such as chemotherapy, radiotherapy, and endocrine therapy can eradicate this, reducing long term breast cancer mortality risk. Many patients receive several adjuvant treatments. Randomised trials have shown that most adjuvant breast cancer therapies reduce breast cancer mortality by similar proportional amounts in both women with good prognosis and those with a poor prognosis.232425 By contrast, the absolute benefits of a given treatment will vary according to prognosis; for a patient with poor prognosis, the absolute mortality benefit from a particular treatment will be greater than for a patient with good prognosis. The decreases in breast cancer mortality during the past few decades suggest that the absolute benefit conferred by each individual adjuvant treatment being considered is likely to be smaller now than in the past.

Secondly, low event rates may make achieving adequate power in randomised trials difficult. Current and planned breast cancer trials may therefore need to recruit large numbers of women to meet their objectives. The likely event rates can be estimated using our data. For example, for a trial including women with tumour size 1-20 mm, grade 1 or 2, node negative, and oestrogen receptor positive disease, five year breast cancer mortality and its confidence interval may be estimated by combining the relevant rows in table S9. In this example, combining rows provides an estimated five year breast cancer mortality of 0.66% (95% confidence interval 0.59% to 0.73%). Details of the method are given in section 6.3 of text S2.

Thirdly, these decreases in breast cancer mortality have implications for decision aids, which are used in the clinic to estimate the absolute benefits of adjuvant treatments for individual patients. These models combine proportional treatment effects in past randomised trials with breast cancer mortality risks in observational data to estimate absolute treatment effects for individual patients. These models need to include breast cancer mortality risks in women with a recent diagnosis, while taking into account the long term mortality patterns seen in women with diagnoses made in the more distant past.

Finally, our findings can be used to reassure most women treated for early breast cancer that they are likely to become long term survivors. They can also be used to identify the groups of women for whom the risk of breast cancer mortality remains substantial.

What is already known on this topic

  • The risk of mortality from breast cancer after a diagnosis of early invasive breast cancer has decreased during the past few decades

  • The extent of the decrease in mortality is unknown, as is whether it is limited to patients with certain characteristics or whether it applies to all patients

  • Detailed, population based estimates of breast cancer mortality risks according to routinely available factors are unavailable for patients treated recently

What this study adds

  • Since the 1990s, the five year risk of death from breast cancer has decreased from 14.4% to 4.9% overall, with reductions seen in nearly all patient groups

  • Five year risks of death from breast cancer in patients with a recent diagnosis varied widely and were under 3% for 62.8% of women but 20% or higher for 4.6% of women

Ethics statements

Ethical approval

The study was approved by the National Cancer Registration and Analysis Service’s Office for Data Release (reference ODR1718_390). The study makes use of routinely collected data that was de-personalised before release. Informed consent from individual participants was therefore not required.

Data availability statement

De-personalised study data may be made available on request to accredited researchers who submit a proposal that is approved by NHS England’s Data Access Request Service (DARS).

Acknowledgments

We thank Jem Rashbass, former director of the National Disease Registration Service at Public Health England, for his help in initiating this study, Kezia Gaitskell for her work in reviewing the pathological information for the women in the study, and patient contributors Hilary Stobart and Mairead MacKenzie on behalf of Independent Cancer Patients’ Voice for helpful discussions. This work uses data provided by patients and collected by the NHS as part of their care and support. The data were collated, maintained, and quality assured by NCRAS, which is part of the National Disease Registration Service.

Footnotes

  • Contributors: CT, PM, GM, DJC, and SCD contributed to the study conceptualisation and design. PM, JP, JB, JC, DJC, GM, and DD contributed to the data collection and collation. PM, JP, and SCD contributed to the data analysis. CT, SD, and AJK contributed to the literature review. CT, PM, JP, SCD, TW, DJC, GM, and DD participated in data interpretation and manuscript preparation. All authors read and approved the final version. CT, PM, and JP contributed equally (joint first authors). DJC, GM, and DD contributed equally (joint last authors). CT and PM are the guarantors. The corresponding author attests that all listed authors meet authorship criteria, and that no others meeting the criteria have been omitted.

  • Funding: Funding was provided by Cancer Research UK (grants C8225/A21133, PRCRPG-Nov21\100001), the National Institute for Health Research Oxford Biomedical Research Centre, the UK Medical Research Council, and the University of Oxford. These funding bodies had no role in study design; in the collection, analysis, and interpretation of data; in the writing of the report; or in the decision to submit the article for publication.

  • Competing interests: All authors have completed the ICMJE uniform disclosure form at https://www.icmje.org/disclosure-of-interest/ and declare: support from Cancer Research UK, the National Institute for Health Research Oxford Biomedical Research Centre, the UK Medical Research Council, and the University of Oxford for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.

  • The lead authors affirm that this manuscript is an honest, accurate and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned (and, if relevant, registered) have been explained.

  • Dissemination to participants and related patient and public communities: The results have already been presented to representatives of appropriate patient groups. After publication, they will be freely available both to stakeholders and to the broader public.

  • Provenance and peer review: Not commissioned; externally peer reviewed.

http://creativecommons.org/licenses/by/4.0/

This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY 4.0) license, which permits others to distribute, remix, adapt and build upon this work, for commercial use, provided the original work is properly cited. See: http://creativecommons.org/licenses/by/4.0/.

References