Skip to main content
  • Research article
  • Open access
  • Published:

The Clinicopathological features and survival outcomes of patients with different metastatic sites in stage IV breast cancer



The features and survival of stage IV breast cancer patients with different metastatic sites are poorly understood. This study aims to examine the clinicopathological features and survival of stage IV breast cancer patients according to different metastatic sites.


Using the Surveillance, Epidemiology, and End Results database, we restricted our study population to stage IV breast cancer patients diagnosed between 2010 to 2015. The clinicopathological features were examined by chi-square tests. Breast cancer-specific survival (BCSS) and overall survival (OS) were compared among patients with different metastatic sites by the Kaplan-Meier method with log-rank test. Univariable and multivariable analyses were also performed using the Cox proportional hazard model to identify statistically significant prognostic factors.


A total of 18,322 patients were identified for survival analysis. Bone-only metastasis accounted for 39.80% of patients, followed by multiple metastasis (33.07%), lung metastasis (10.94%), liver metastasis (7.34%), other metastasis (7.34%), and brain metastasis (1.51%). The Kaplan-Meier plots showed that patients with bone metastasis had the best survival, while patients with brain metastasis had the worst survival in both BCSS and OS (p < 0.001, for both). Multivariable analyses showed that age, race, marital status, grade, tumor subtype, tumor size, surgery of primary cancer, and a history of radiotherapy or chemotherapy were independent prognostic factors.


Stage IV breast cancer patients have different clinicopathological characteristics and survival outcomes according to different metastatic sites. Patients with bone metastasis have the best prognosis, and brain metastasis is the most aggressive subgroup.

Peer Review reports


Breast cancer (BC), the most common malignancy in women, is also the key contributor to cancer-related deaths worldwide in female. 2.08 million women were diagnosed with breast cancer in 2018, occupying 24.2% of all new cancer cases in that year, according to the 2018 global cancer statistics. In addition, 626,679 women are expected to die from this disease, making up 15% of all cancer deaths among females [1]. Despite the high morbidity of the disease, patients with breast cancer have better prognosis compared to other aggressive cancers. According to the American Cancer Society’s biennial update on female breast cancer statistics, the 5-year survival rate is 91% and 10-year survival rate is 84% [2]. However, the survival decreases greatly if the patients develop distant metastases. The overall 5-year relative survival rate is 99% for localized diseases and 86% for regional diseases, which drops to 27% for distant-stage diseases [3].

It is estimated that 20–30% of early stage breast cancers will go on to develop metastatic disease, and 6–10% of all women suffering from breast cancer in the United States were found to have stage IV disease at diagnosis (called de novo metastatic breast cancer) [4]. These patients have poor prognoses with a median survival time of 2–3 years [5], despite technology development and medical progresses that have been made nowadays. Metastatic breast cancer is a heterogeneous disease with distinct prognoses [6], which are affected by many clinicopathological features of patients, such as age, race, marital status, performance status as well as tumor size, lymph nodes status, metastatic sites, number of metastatic sites, pathological or genotype characteristics and previous medical treatments [7,8,9]. Therefore, an accurate estimation of survival may benefit patients significantly in decision-making [10]. The TNM stage is the widely accepted tool to predict the prognosis of patients, but it deals with a limited range of factors and ignores patient-specific conditions, pathological or genotype characteristics, and treatments [10]. Therefore, it is still difficult to make precise predictions about individual prognoses for metastatic breast cancer using this method. Among all the predictors of outcome mentioned above, the metastatic site may be the most important factor with the greatest significant impact on further treatment regimens.

Population-level estimates for prognosis among breast cancer patients with distant metastases are lacking, and the relationship between clinical-related factors and the exact patterns of distant metastasis is not well elaborated. And thus, the purposes of this study were to use the Surveillance, Epidemiology, and End Results (SEER) database to present clinicopathological features and survival estimates, and to investigate the prognostic impact of stage IV breast cancer patients based on their metastatic sites at the time of cancer diagnosis on a population-based level.



The recent version of the Surveillance, Epidemiology and End Results (SEER) 18 registries Custom Data (with additional treatment fields) was used as the data source for the present population-based investigation. SEER*Stat Software version 8.3.5 ( (Information Management Service, Inc. Calverton, MD, USA) was used to generate the case listing. Maintained by the National Cancer Institute, the SEER program is the largest publicly available cancer dataset in the world, which consists of 18 population-based cancer registries and covers approximately 26% of the US population across several geographic regions [11]. Patients diagnosed with female breast cancer as the primary cancer from 2010 to 2015, with distant metastasis, were enrolled into the study. The tumors were classified based on their primary site of presentation and histology utilizing the International Classification of Disease for Oncology, Third Edition (ICD-O-3). Patients with unknown metastatic sites and survival were excluded. Of 19,913 women with a diagnosis of stage IV breast cancer between the year 2010 and 2015 included into the SEER Registry, 18,322 women were eligible for inclusion into the present study (Fig. 1, Flowchart). Eligible patients were grouped according to metastatic site. All procedures performed in studies involving human participants were in accordance with the 1964 Helsinki declaration and its later amendments, or comparable ethical standards. Because the SEER database is publicly accessible, this study does not require informed patient consent and was deemed exempt from review by the Ethics Committee of the First Affiliated Hospital of Xian Jiaotong University.

Fig. 1
figure 1

Flowchart for patient selection from the Surveillance, Epidemiology and End Results (SEER) database

Demographic and clinical variables

The relationship between metastatic sites and clinical characteristics, including age at diagnosis, year of diagnosis, race, marital status, tumor grade, tumor size, nodal status, subtype, and treatment was analyzed. Initial metastatic sites were registered as single or multiple and were categorized as bone-only, lung-only, liver-only, brain-only, other-only (metastatic sites other than the aforementioned sites such as skin and soft tissue) and multiple metastasis (more than one of the aforementioned distant metastatic sites). Overall survival (OS) was calculated from the date of diagnosis to the date of death due to any cause, the date of last follow-up, or December 31, 2015. Breast cancer specific survival (BCSS) was measured as the time from the date of diagnosis to the date of death attributed to breast cancer. Both overall survival and breast cancer specific survival were used as endpoints.

Statistical analysis

The baseline characteristics of patients and treatment were described using summary statistics, with continuous variables being shown as mean ± standard deviation. Differences between qualitative variables and continuous variables were analyzed using χ2 statistics and analysis of variance, respectively. Survival was estimated using the Kaplan–Meier method and compared between the different metastatic groups using log-rank test. Hazard ratios (HRs) and 95% confidence intervals (CIs) were calculated by univariable and multivariable Cox proportional hazard models to assess the relative contribution of metastatic sites and other factors to survival after diagnosis of metastatic breast cancer. The statistical analyses above were carried out using SPSS software version 22.0 (SPSS Inc.). Cumulative incidence function (CIF) was employed to show the probability of each competing event and the differences between the groups were estimated using Gray’s test. Competing risk analysis were performed using R version 3.6.1 software (The R Foundation for Statistical Computing, Vienna, Austria; with the R packages cmprsk and survival. All the tests above were 2-tailed, and a p value of less than 0.05 was considered statistically significant.


Patient characteristics

According to the above inclusion and exclusion criteria, a total of 18,322 patients were enrolled in this study. Table 1 shows baseline demographic and clinicopathological characteristics of the included patients according to metastatic sites. Median survival after diagnosis of metastatic disease was 26.0 months, with 8442 patients (46.1%) being alive at the end of the follow-up period.

Table 1 Patient characteristics

The largest subgroup was the bone metastasis, comprising 39.8% of all patients (7292), followed by multiple metastasis (33.07%, 6059), lung metastasis (10.94%, 2005), liver metastasis (7.34%, 1346), other metastasis (7.34%, 1344), and brain metastasis (1.51%, 276). The median patient age at initial diagnosis was 62 years (range 20–100 years). The mean patient age for liver metastasis was the lowest (59.0 years), while lung metastasis was the highest (66.2 years) among all metastasis groups (range 61.2–64.3 years). Poorly or undifferentiated tumors were most common in liver metastasis (48.66%), followed by lung metastasis (47.16%) and brain metastasis (42.39%). Patients with lung metastasis (50.69%) and multiple metastases (48.07%) tend to have larger tumors of T3-T4 at initial diagnosis. Furthermore, these two groups of patients also have a higher proportion of later N stage (N3), making up 29.48 and 33.94% respectively. The distribution of subtypes is significantly different between initial metastatic sites (Table 1). Bone was the predominant initial site of metastasis for the HR+/HER2- (64.4%) group and the least common site in the HR−/HER2+ group (3.24%). Of patients with brain, lung and liver metastasis, 26.9, 21.22, and 14.42%, respectively, had the triple negative subtype. Detailed patient characteristics are presented in Table 1.

Kaplan-Meier Survival analysis

Of all the 18,322 patients finally recruited, 9880 patients had died by the end of the last follow-up, 7239 of whom died of breast cancer specifically. The median survival of all patients since diagnosis of distant metastases was 26.0 months. The median survival ranged from 8.0 months (95% CI 5.65–10.35 months) for patients with brain metastasis to 36.0 months (95% CI 34.74–37.27 months) for patients with bone metastasis.

The Kaplan-Meier plots were displayed in Fig. 2 to show the survival of all populations. Figure 2a shows the overall survival (OS) of stage IV patients who were enrolled into the study according to metastatic sites. Patients with bone metastasis had the best survival, with 3-year OS rate of 50.5%, followed by patients with other metastasis, liver metastasis, and lung metastasis, (with OS rate of 41.9, 38.2, and 37.5% respectively). Patients with brain metastasis and multiple metastasis had worse OS than other subgroups: the 3-year OS rate was 19.9, and 27.4%, respectively (p < 0.001). Figure 2b shows the breast cancer specific survival (BCSS) according to metastatic sites of stage IV patients. Similar to OS, brain metastasis had the worst survival, with a 3-year BCSS rate of 50.6%. The 3-year BCSS rate with liver, lung, and multiple metastases were, 67.7, 66.1, and 62.6% respectively. Bone and other metastasis had better BCSS, for which the 3-year BCSS rate was 76.6 and 65.3%, respectively (p < 0.001).

Fig. 2
figure 2

Survival curves with the log-rank tests of overall survival (OS, a, p < 0.001) and breast cancer-specific survival (BCSS, b, p < 0.001) based on metastatic sites for breast cancer patients

Breast cancer subtype is an important prognostic factor since it affects survival significantly. The triple negative pattern tends to decrease survival dramatically compared to the other groups (Additional file 1: Figure S1). Survival estimates overall (Fig. 3) and breast cancer specific (Additional file 2: Figure S2) as stratified by subtype are graphically displayed in the Figure; HR+/HER2- (A), HR+/HER2 + (B), HR−/HER2+ (C), and TN (D). In the first three subtypes, worst survival was seen in the brain metastasis and best survival in the bone metastasis, while in triple negative subtype, patients with multiple metastases had the worst prognoses.

Fig. 3
figure 3

Survival curves with the log-rank tests of overall survival per metastatic sites according to subtype; HR+/HER2-(a), HR+/HER2 + (b), HR−/HER2 + (c), TN (d). Abbreviations: HR: Hormone receptor, HER2: Human epidermal growth receptor 2, TN: Triple negative

Cox regression analysis of survival

In order to further figure out the effect of multiple factors on BCSS and OS, the Cox proportional hazard model was applied to the analysis. Univariable analysis of BCSS and OS proved unmarried marital status, African descent, high tumor grade, large tumor size, later N stage, and the Her2 positive and triple negative cancer subtypes to be distinct risk factors for poor survival (hazard ratio [HR] > 1, p < 0.001). By contrast, married status, other race, apply of chemotherapy, radiotherapy, and surgery were found to be protective factors for better survival (hazard ratio [HR] < 1, p < 0.001) (Table 2). As for metastatic sites, the results were consistent with those from the Kaplan-Meier analysis. Patients with bone metastasis had the best BCSS and OS, followed by patients with other, liver, lung, and multiple metastases. Specifically, patients with brain metastasis exhibited the worst BCSS, with a hazard ratio of 1.708 (95% confidence interval [CI] = 1.442–2.023, p < 0.001), and OS (HR = 2.492, 95% CI =2.161–2.874, p < 0.001), when compared to bone metastasis.

Table 2 Univariate Analysis of Prognostic factors of BCSS and OS in Metastatic Breast Cancer

The variables age, race, marital status, grade, subtype, T stage, N stage, surgery, radiotherapy, chemotherapy and metastatic sites were subsequently analyzed with the multivariable Cox analysis. After adjusting for age, race/ethnicity, marital status, tumor grade, breast cancer subtype, tumor size, nodal status, surgery, radiotherapy, and chemotherapy in the analysis, metastatic site remained an independent prognostic factor of BCSS (p < 0.001) and OS (p < 0.001). The detailed results of Cox regression analysis of BCSS and OS can be found in Table 3. Compared to patients with bone metastasis, the BCSS (HR 0.994, 95% CI 0.881–1.122, p = 0.921) and OS (HR 0.994, 95% CI 0.897–1.100, p = 0.902) of patients with lung metastasis were not significantly different. Patients with other metastasis had similar BCSS (HR 0.955, 95% CI 0.827–1.103, p = 0.532) as patients with bone metastasis, but worse OS (HR 1.127, 95% CI 1.001–1.269, p = 0.048). Patients with liver (HR 1.384, 95% CI 1.208–1.586, p < 0.001; OS, HR 1.428, 95% CI 1.272–1.602, p < 0.001) and multiple metastases (BCSS, HR 1.475, 95% CI 1.361–1.599, p < 0.001; OS, HR 1.806, 95% CI 1.684–1.937, p < 0.001) had worse BCSS and OS than bone metastasis. Similarly, multivariable analysis also indicated that those patients with brain metastasis had significant inferior BCSS (HR 1.975, 95% CI 1.551–2.514, p < 0.001) and OS (HR 2.307, 95% CI 1.862–2.859, p < 0.001) over other metastatic sites. Age (p < 0.001), marital status (p < 0.001), tumor grade (p < 0.001), triple negative subtype (p < 0.001), tumor size (p < 0.001), surgery (p < 0.001), chemotherapy (p < 0.001) were also statistically significant in a multivariable analysis with forward selection.

Table 3 Multivariate Analysis of Prognostic factors of BCSS and OS in Metastatic Breast Cancer

We also did a competing risk analysis to estimate the marginal probability for each cause of death. The competing risk analysis on cause of death for the entire cohort and according to metastatic sites is reported in supplementary data (Additional file 3: Figure S3). Of the cohort, the respective 5-year estimates of breast cancer-specific mortality and other cause of mortality were 56.2 and 20.7%.


In the current study, we analyzed the clinicopathological characteristics and survival of de novo metastatic breast cancer according to metastatic sites using the SEER data. Metastatic breast cancer (MBC) is widely known to have a poor prognosis compared to non-metastatic breast cancer. Bone, liver, lung, and brain are common sites of distant metastasis in breast cancer [12]. Previous studies have shown that bone is the most common distant metastatic organ in breast cancer patients [7, 13]. Similarly, our results also show that bone metastasis is the most prevalent subgroup among the study cohort, accounting for 39.8% of the total patients.

Recently published studies proposed that the metastatic patterns may differ from one another among different biological types of breast cancer. And it has also been demonstrated to have different prognostic impact on different patterns of distant metastasis. In this study, we showed that the patients with bone metastasis mainly have the HR+/HER2- subtype, which is supported by several studies that patients with HR+ are more prone to develop bone metastasis [14, 15]. The predictive value and the relationship between visceral metastasis and breast cancer subtype are still controversial. Kennecke et al. concluded that HER2-enriched breast tumors were more likely to develop liver metastases when compared to luminal A tumors [7]. However, there are some other studies reported that liver metastasis was not associated with breast cancer subtype [16]. In this study, HR−/HER2+ tumors exhibited highest rates of liver metastasis, while HR+/HER2- exhibited highest rates of lung metastasis. The brain has been previously described as a preferred site of metastasis among triple negative cancers. Martin, et al. previously reported on the incidence proportion of breast cancer patients with brain metastases from 2010 to 2013. They found that the incidence proportion of brain metastasis was highest among patients with HR−/HER2+ and triple-negative subtypes (11.37 and 11.45%, respectively) to any distant sites at diagnosis of breast cancer [17]. Other studies based on incidence of brain involvement among patients with metastatic HER-2 positive tumor has been reported to be from 25 to 34% [18], which is significantly higher than the 3.4 and 1.9% rate for triple negative and HR−/HER2+ subtypes respectively in another study [8]. The main reason for this is that there is a very limited number of patients with brain metastasis enrolled in this study compared to other subgroups, accounting for only 1.51% of the whole cohort population, which makes it not as representative as other results. Further studies with more patients and prospective design are needed to better interpret the situation.

Some studies have shown that the survival of female breast cancer may be associated with different metastatic patterns [12, 19]. However, these studies are controversial, and there’s no population-based study focus on studying the survival differences of patients with different metastatic pattern. Gerratana et al. reported that breast cancer patients with lung as the first site of distant metastasis had the best survival outcome (58.5 months) compared with those with bone (44.4 months), liver (36.7 months), or brain (7.35 months) as the first metastasis identified [12]. But the 5-year survival elevated significantly in patients with bone metastasis when compared with other metastasis in another study [19], which is consistent with our conclusion. In our study, the Kaplan–Meier curves showed that patients with bone metastasis had the best survival in both OS and BCSS. Lung, liver and other metastasis had very similar survival according to Kaplan-Meier analysis. Women with multiple organs affected had worse survival than those with minimal metastatic disease, a finding confirmed in this study. The presence of brain metastasis was shown in this analysis as a poor prognostic factor compared to all the other distant metastases, worse than patients with multiple metastases. We got the same result in both Kaplan-Meier survival analysis and Cox regression analysis. Taking bone metastasis as the reference, univariable and multivariable analysis showed that patients with brain metastasis had the worst prognosis. The differences of the reported outcomes may be partially due to the clinicopathological and molecular characteristics specific to different subtypes of BC. This result is in line with another study based on SEER database [17] as well as some retrospective studies [20, 21]. Of note, our study may be more comprehensively offering better guidance for the daily clinical practice, given that our analysis was based on a much larger population of BC patients.

Besides de novo MBC, most MBC are developed from primary cancer, which is known as recurrent MBC. The majority of BC patients diagnosed with early stages of disease are eligible for surgical resections and adjuvant systemic therapies, leading to longer disease-free survival. While, those develop metastatic disease during or after cycles of standard care are barely curable. A recent study conducted by Rueda et al. presented a statistical model based on distinct disease stages and competing risks of mortality to predict individual risk of distant recurrence [22]. They applied the model to a population of 3240 patients and analyzed the rates and routes of metastases and their lethality, suggesting that ER-negative patients harbored more visceral disease than the ER-positive cases. It took longer time to develop bone metastases which were more common in ER-positive than in ER-negative patients (71% versus 43%). They also denoted that this model was of more predictive value in predicting brain metastases incorporating the influence of the number and site of relapses on the risk of death after recurrence. We had the same conclusion in this study, which indicates that both de novo MBC and recurrent MBC share the same molecular characteristics and some of the biological behaviors. As for survival and prognosis, previous studies have shown that recurrent MBC patients tend to have a worse prognosis relative to de novo MBC patients. Yamamura et al. studied differences in survival outcomes between de novo and recurrent MBC patients, and found that de novo MBC led to a better prognosis than recurrent MBC with a disease-free interval of less than 2 years or an AFI (interval from the end of adjuvant treatment to the first recurrence) less than 1 year [23]. The reason for this may be that de novo MBC might be more sensitive to systemic treatments for its treatment-naïve status. Recurrent MBC, however, might be more resistant to systemic treatments after receiving adjuvant therapies. Other possible biological differences besides therapeutic factors between de novo and recurrent MBC patients should also be considered.

Historically, women diagnosed with metastatic breast cancer were not treated with surgery and received only systemic therapy [24]. It is believed that surgical resection of the primary tumor was palliative and performed only to relieve symptoms such as bleeding, infection, or pain. However, value of local treatment of the primary in cases of a metastatic solid tumor has been shown for metastatic renal cell carcinoma, nonfunctioning pancreatic neuroendocrine tumors and hepatocellular carcinoma [25,26,27]. However, there hasn’t been an evidence-based consensus on surgery for patients with metastatic BC. A similar strategy is currently being explored in several ongoing studies for patients with MBC at initial presentation [28]. Some former retrospective data suggested no benefit from surgical resection of the primary tumor and metastases, since surgeries may increase the progression of distant metastases, in spite of a better local control of breast cancer. But an increasing number of studies suggested that surgical resection of primary tumor for highly selected stage IV patients who had favorable responses to systemic neoadjuvant chemotherapy should be considered [29, 30]. This is consistent with observations from our study describing a survival benefit for both BCSS and OS of the patients accepted surgery compared to those who did not undergo surgery. The conclusion is also true when it comes to radiotherapy and chemotherapy. The multivariable Cox proportional hazards model showed both radiotherapy and chemotherapy were significantly associated with improved survival of patients in this cohort. However, the study still carries the methodological defects of a retrospective analysis, so it is impossible to derive clear recommendations. Further randomized clinical trials are needed to reach conclusions on the benefits and risks of breast local treatment including surgery and radiotherapy together with systemic treatment for women diagnosed with MBC.

In this study, age at diagnosis, race, marital status and tumor grade are prognostic factors that influenced survival besides initial sites of metastases, which is in accordance with other studies investigating prognostic factors in MBC [8, 31]. The current analysis showed that married patients have better overall and breast cancer-specific survival compared to unmarried patients. This difference may be explained by the social psychological support given by a partner, which greatly benefits the patient. In addition to tumor characteristics as mentioned, menstrual history (age at menarche, cycle length, irregular menstruation, lifetime number of menstrual cycles, menopausal status at diagnosis and age at menopause) and reproductive factors (parity, number of full-term pregnancy, age at first full-term pregnancy, and breastfeeding) are among the most well-established risk factors for breast cancer. A meta-analysis revealed a decrease in risk by 4–9% for each increased year of postponement in menarche, 7% for each extra birth, and 4% for every addition of 12 months of breast feeding, and an increase in risk by 3–5% each increased year of age at first birth [32]. There are also studies investigating influence of reproductive factors on tumor characteristics and breast cancer survival in women. An association had been found between early age at menarche and reduced survival. Age at menarche was also shown to be significantly associated with tumor grade and lymph node involvement [33, 34]. Unfortunately, we are unable to do such analysis since there’s no female history data in the SEER database. Further studies with focus on this point of view are needed to better delineate the picture.

With the increase of cancer incidence all over the world, a growing number of patients suffer from cancer invasion and metastasis. Metastatic breast cancer patients could be composed of heterogeneous groups and present with totally different metastatic patterns and prognostic outcomes, thus distinct treatment strategies are required. We hope to provide deeper insights into the heterogeneity of metastatic breast cancer through this population-based study. It may help to determine whether women with metastatic breast cancer after diagnosis might benefit from extended therapy to improve survival. Specifically, the brain metastasis is the most lethal subgroup of metastatic breast cancer that should be paid more attention after diagnosis, paving the way for individualized management strategies for these high-risk patient populations.

Inevitably, there are several limitations in this study that should not be overlooked. First, as a retrospective study rather than a prospective cohort study, inherent selection biases cannot be avoided and could limit the external validity of this study. Second, information about disease recurrence or subsequent sites of disease involvement is not provided in the SEER database, so that the focus of this study is de novo metastatic breast cancer, and no conclusions can be made about patients who developed metastases later in their disease course. Third, the follow-up period was relatively short with the median follow-up time of only 14 months, as the data of HER-2 status were not available until 2010. Thus, the study only focused on the short-term prognosis of metastatic breast cancer patients. These limitations may have contributed to study bias and undermined the power of analysis.


Despite these limitations, our results provide insight into the epidemiology, clinicopathological characteristics and survival outcomes of distant metastases in patients with newly diagnosed breast cancer in the United States. Even though patients with distant metastases are all defined as advanced cancer, the prognosis varies greatly depending on metastatic site. Patients with bone metastasis and other metastatic sites other than bone, liver, lung, or brain have better prognoses, while brain metastasis is the most aggressive group and requires additional treatment options. The metastatic sites should be taken into consideration when making therapeutic strategies for patients with advanced cancer. This knowledge is useful for the treatment of metastatic breast cancer patients in this era of individualized therapy.

Availability of data and materials

The data were abstracted from an open database, the Surveillance, Epidemiology, and End Results (SEER) database (



Breast cancer


Breast cancer-specific survival


Confidence interval


Human epidermal growth factor receptor 2


Hazard ratio


Overall survival


Surveillance, epidemiology and end results


Triple negative


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

    Article  Google Scholar 

  2. DeSantis CE, Ma J, Gaudet MM, Newman LA, Miller KD, Goding Sauer A, Jemal A, Siegel RL. Breast cancer statistics, 2019. CA Cancer J Clin:2019.

  3. Noone AMHN, Krapcho M, Miller D, Brest A, Yu M, Ruhl J, Tatalovich Z, Mariotto A, Lewis DR, Chen HS, Feuer EJ, Cronin KA, editors. SEER Cancer statistics review, 1975–2015: National Cancer Institute; 2018.

  4. Howlader NNA, Krapcho M. SEER Fast Stats, 1975-2014. Stage Distrib. 2018:2005–14.

  5. Mariotto AB, Etzioni R, Hurlbert M, Penberthy L, Mayer M. Estimation of the number of women living with metastatic breast Cancer in the United States. Cancer Epidemiol Biomarkers Prev. 2017;26(6):809–15.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Peart O. Metastatic breast Cancer. Radiol Technol. 2017;88(5):519M–39M.

    PubMed  Google Scholar 

  7. Kennecke H, Yerushalmi R, Woods R, Cheang MCU, Voduc D, Speers CH, Nielsen TO, Gelmon K. Metastatic behavior of breast Cancer subtypes. J Clin Oncol. 2010;28(20):3271–7.

    Article  PubMed  Google Scholar 

  8. Lobbezoo DJ, van Kampen RJ, Voogd AC, Dercksen MW, van den Berkmortel F, Smilde TJ, van de Wouw AJ, Peters FP, van Riel JM, Peters NA, et al. Prognosis of metastatic breast cancer subtypes: the hormone receptor/HER2-positive subtype is associated with the most favorable outcome. Breast Cancer Res Treat. 2013;141(3):507–14.

    Article  CAS  PubMed  Google Scholar 

  9. Largillier R, Ferrero JM, Doyen J, Barriere J, Namer M, Mari V, Courdi A, Hannoun-Levi JM, Ettore F, Birtwisle-Peyrottes I, et al. Prognostic factors in 1,038 women with metastatic breast cancer. Ann Oncol. 2008;19(12):2012–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol. 2015;16(4):e173–80.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Wingo PAJP, Hiatt RA. Building the infrastructure for nationwide cancer surveillance and control–a comparison between the National Program of Cancer registries (NPCR) and the surveillance, epidemiology, and end results (SEER) program (United States). Cancer Causes Control. 2003;14:175–93.

    Article  PubMed  Google Scholar 

  12. Gerratana L, Fanotto V, Bonotto M, Bolzonello S, Minisini AM, Fasola G, Puglisi F. Pattern of metastasis and outcome in patients with breast cancer. Clin Exp Metastasis. 2015;32(2):125–33.

    Article  CAS  PubMed  Google Scholar 

  13. Wu SG, Sun JY, Yang LC, Tang LY, Wang X, Chen XT, Liu GH, Lin HX, Lin Q, He ZY. Patterns of distant metastasis in Chinese women according to breast cancer subtypes. Oncotarget. 2016;7(30):47975–84.

    PubMed  PubMed Central  Google Scholar 

  14. Liede A, Jerzak KJ, Hernandez RK, Wade SW, Sun P, Narod SA. The incidence of bone metastasis after early-stage breast cancer in Canada. Breast Cancer Res Treat. 2016;156(3):587–95.

    Article  CAS  PubMed  Google Scholar 

  15. Kai M, Kogawa T, Liu DD, Fouad TM, Kai K, Niikura N, Hsu L, Willey JS, Theriault RL, Valero V, et al. Clinical characteristics and outcome of bone-only metastasis in inflammatory and noninflammatory breast cancers. Clin Breast Cancer. 2015;15(1):37–42.

    Article  PubMed  Google Scholar 

  16. Park HS, Kim S, Kim K, Yoo H, Chae BJ, Bae JS, Song BJ, Jung SS. Pattern of distant recurrence according to the molecular subtypes in Korean women with breast cancer. World J Surg Oncol. 2012;10:4.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Martin AM, Cagney DN, Catalano PJ, Warren LE, Bellon JR, Punglia RS, Claus EB, Lee EQ, Wen PY, Haas-Kogan DA, et al. Brain metastases in newly diagnosed breast Cancer: a population-based study. JAMA Oncol. 2017;3(8):1069–77.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Gaedcke J, Traub F, Milde S, Wilkens L, Stan A, Ostertag H, Christgen M, von Wasielewski R, Kreipe HH. Predominance of the basal type and HER-2/neu type in brain metastasis from breast cancer. Mod Pathol. 2007;20(8):864–70.

    Article  CAS  PubMed  Google Scholar 

  19. Kast K, Link T, Friedrich K, Petzold A, Niedostatek A, Schoffer O, Werner C, Klug SJ, Werner A, Gatzweiler A, et al. Impact of breast cancer subtypes and patterns of metastasis on outcome. Breast Cancer Res Treat. 2015;150(3):621–9.

    Article  CAS  PubMed  Google Scholar 

  20. Ording AG, Heide-Jorgensen U, Christiansen CF, Norgaard M, Acquavella J, Sorensen HT. Site of metastasis and breast cancer mortality: a Danish nationwide registry-based cohort study. Clin Exp Metastasis. 2017;34(1):93–101.

    Article  PubMed  Google Scholar 

  21. Gerdan L, Segedin B, Nagy V, Khoa MT, Trang NT, Schild SE, Rades D. The number of involved extracranial organs: a new predictor of survival in breast cancer patients with brain metastasis. Clin Neurol Neurosurg. 2013;115(10):2108–10.

    Article  PubMed  Google Scholar 

  22. Rueda OM, Sammut SJ, Seoane JA, Chin SF, Caswell-Jin JL, Callari M, Batra R, Pereira B, Bruna A, Ali HR, et al. Dynamics of breast-cancer relapse reveal late-recurring ER-positive genomic subgroups. Nature. 2019;567(7748):399–404.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Yamamura J, Kamigaki S, Fujita J, Osato H, Komoike Y. The difference in prognostic outcomes between De novo stage IV and recurrent metastatic patients with hormone receptor-positive, HER2-negative breast Cancer. In Vivo. 2018;32(2):353–8.

    PubMed  PubMed Central  Google Scholar 

  24. Khan SA, Stewart AK, Morrow M. Does aggressive local therapy improve survival in metastatic breast cancer? Surgery. 2002;132(4):620–6 discussion 626-627.

    Article  PubMed  Google Scholar 

  25. Heng DY, Wells JC, Rini BI, Beuselinck B, Lee JL, Knox JJ, Bjarnason GA, Pal SK, Kollmannsberger CK, Yuasa T, et al. Cytoreductive nephrectomy in patients with synchronous metastases from renal cell carcinoma: results from the international metastatic renal cell carcinoma database consortium. Eur Urol. 2014;66(4):704–10.

    Article  PubMed  Google Scholar 

  26. Keutgen XM, Nilubol N, Glanville J, Sadowski SM, Liewehr DJ, Venzon DJ, Steinberg SM, Kebebew E. Resection of primary tumor site is associated with prolonged survival in metastatic nonfunctioning pancreatic neuroendocrine tumors. Surgery. 2016;159(1):311–8.

    Article  PubMed  Google Scholar 

  27. Abdel-Rahman O. Role of liver-directed local tumor therapy in the management of hepatocellular carcinoma with extrahepatic metastases: a SEER database analysis. Expert Rev Gastroenterol Hepatol. 2017;11(2):183–9.

    Article  CAS  PubMed  Google Scholar 

  28. Ruiterkamp J, Voogd AC, Tjan-Heijnen VC, Bosscha K, van der Linden YM, Rutgers EJ, Boven E, van der Sangen MJ, Ernst MF, Dutch breast Cancer Trialists G. SUBMIT: systemic therapy with or without up front surgery of the primary tumor in breast cancer patients with distant metastases at initial presentation. BMC Surg. 2012;12:5.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Badwe R, Hawaldar R, Nair N, Kaushik R, Parmar V, Siddique S, Budrukkar A, Mittra I, Gupta S. Locoregional treatment versus no treatment of the primary tumour in metastatic breast cancer: an open-label randomised controlled trial. Lancet Oncol. 2015;16(13):1380–8.

    Article  PubMed  Google Scholar 

  30. Soran A, Ozmen V, Ozbas S, Karanlik H, Muslumanoglu M, Igci A, Canturk Z, Utkan Z, Ozaslan C, Evrensel T, et al. A randomized controlled trial evaluating resection of the primary breast tumor in women presenting with de novo stage IV breast cancer: Turkish Study (Protocol MF07–01). J Clin Oncol. 2016;34(15_suppl):1005.

    Article  Google Scholar 

  31. Lobbezoo DJ, van Kampen RJ, Voogd AC, Dercksen MW, van den Berkmortel F, Smilde TJ, van de Wouw AJ, Peters FP, van Riel JM, Peters NA, et al. Prognosis of metastatic breast cancer: are there differences between patients with de novo and recurrent metastatic breast cancer? Br J Cancer. 2015;112(9):1445–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Cancer CGoHFiB. Breast cancer and breastfeeding: collaborative reanalysis of individual data from 47 epidemiological studies in 30 countries, including 50 302 women with breast cancer and 96 973 women without the disease. Lancet. 2002;360(9328):187–95.

    Article  Google Scholar 

  33. Orgeas CC, Hall P, Rosenberg LU, Czene K. The influence of menstrual risk factors on tumor characteristics and survival in postmenopausal breast cancer. Breast Cancer Res. 2008;10(6):R107.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Trivers KF, Gammon MD, Abrahamson PE, Lund MJ, Flagg EW, Kaufman JS, Moorman PG, Cai J, Olshan AF, Porter PL, et al. Association between reproductive factors and breast cancer survival in younger women. Breast Cancer Res Treat. 2007;103(1):93–102.

    Article  PubMed  Google Scholar 

Download references


We thank Yupeng Yang from the Department of Surgery, Zhongshan Hospital, Fudan University, for helping us with some of the data analysis. We also thank Chelsy Chen from Johns Hopkins University, for helping us with the revision of the manuscript.


Not applicable.

Author information

Authors and Affiliations



RW, JH and LN designed the study. XL1 and XL2 collected the data. RW and YZ analyzed the data. RW, YZ, and XL1 organized the manuscript. JH and LN reviewed the paper and revised the manuscript. All authors (RW, YZ, XL1, XL2, JH and LN) have read and approved the final manuscript. All authors contributed toward data analysis, drafting and revising the paper and agree to be accountable for all aspects of the work.

Corresponding author

Correspondence to Ligang Niu.

Ethics declarations

Ethics approval and consent to participate

This study used previously collected de-identified data, the need for informed consent had been waived due to the retrospective nature of the study, and was deemed exempt from review by the Ethics Committee of the First Affiliated Hospital of Xi’an Jiaotong University.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Additional files

Additional file 1: Figure S1.

Survival curves with the log-rank tests of overall survival (OS, A, p < 0.001) and breast cancer-specific survival (BCSS, B, p < 0.001) based on subtype. Abbreviations: HR: Hormone receptor, HER2: Human epidermal growth receptor 2, TN: Triple negative.

Additional file 2: Figure S2.

Survival curves with the log-rank tests of breast cancer-specific survival per metastatic sites according to subtype; HR+/HER2-(A), HR+/HER2 + (B), HR−/HER2 + (C), TN(D). Abbreviations: HR: Hormone receptor, HER2: Human epidermal growth receptor 2, TN: Triple negative.

Additional file 3: Figure S3.

Cumulative incidence curves of deaths to show the probability of each competing event in the entire cohort (A) and according to metastatic sites (B). The real line represents breast cancer-specific mortality and the dotted line represents competing mortality. Cumulative incidence curves of deaths to show the probability of death due to breast cancer according to metastatic sites (C).

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, R., Zhu, Y., Liu, X. et al. The Clinicopathological features and survival outcomes of patients with different metastatic sites in stage IV breast cancer. BMC Cancer 19, 1091 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: