Nomogram prediction for breast cancer-specific survival in patients with ER−/PR−/HER2+ breast cancer: an analysis based on the SEER population
Highlight box
Key findings
• A robust nomogram incorporating age, lung metastasis, tumor size, and stage accurately predicts 1-, 2-, and 3-year breast cancer (BC)-specific survival in patients with estrogen receptor (ER)−/progesterone receptor (PR)−/human epidermal growth factor receptor 2 (HER2)+ BC, achieving high areas under the curve (>0.77) and good calibration in both training and validation cohorts.
• Lung metastasis emerged as a critical independent prognostic factor, underscoring the importance of metastatic screening and tailored management in this aggressive BC subtype.
What is known and what is new?
• ER−/PR−/HER2+ BC is an aggressive subtype, and factors like tumor size and grade are associated with prognosis.
• This study establishes a validated competing-risk nomogram—incorporating age, lung metastasis, tumor size, and stage—that accurately predicts 1‑ to 3‑year survival specifically for this subgroup, with lung metastasis identified as a critical novel prognostic indicator.
What is the implication, and what should change now?
• This nomogram enables personalized prognosis estimation for patients with ER−/PR−/HER2+ BC, supporting earlier identification of high-risk cases—particularly those with lung metastasis—and calls for integrating such risk tools into clinical decision-making to optimize surveillance, systemic therapy, and resource allocation.
Introduction
According to the estimates of cancer incidence and mortality released by the International Agency for Research on Cancer in GLOBOCAN 2022, breast cancer (BC) is the most common type of cancer among women, with approximately 2.3 million new cases, accounting for 11.6% of all cancers. This disease is the fourth leading cause of cancer deaths globally, with 666,000 deaths (accounting for 6.9% of all cancer deaths) (1). The treatment and prognosis of BC are linked to the expression of human epidermal growth factor receptor 2 (HER2) and hormone receptor (HR), including the estrogen receptor (ER) and progesterone receptor (PR) (2-5). Positive ER and PR suggest that cells can receive signals from estrogen and/or progesterone, and drugs through the endocrine system can inhibit tumor growth (6,7); HER2 is a transmembrane protein of which its intracellular structural region contains sites that bind to adenosine triphosphate, which can reduce the expression of cancer suppressor genes. Its positive expression indicates that cells are growing and proliferating vigorously, which is correlated with tumor invasion and metastasis. Targeted treatment is often recommended for patients with HER2-positive BC (8,9).
Given the critical correlation of ER, PR, and HER2 with treatment and prognosis in BC patients, BC can be classified into four clinical subtypes based on their expression in tumor cells. Triple negative (ER−/PR−/HER2−), HR-negative and HER2-positive (ER−/PR−/HER2+), HR-positive and HER2-positive, and HR-positive and HER2-negative (10). Based on statistics from the National Cancer Institute and the Centers for Disease Control and Prevention, the American Cancer Society found that from 2017 to 2021, 70% of BC subtypes in the United States were HR+/HER2−, 9% were HR+/HER2+, 4% were HR−/HER2+, and 10% were HR−/HER2−. Although the HR−/HER2+ subtype represents the smallest proportion, it has the second lowest 5-year survival rate at 86%, surpassed only by the triple-negative subtype, which has the poorest prognosis at 78% (11). This may be associated with the higher invasiveness and risk of locoregional recurrence of the ER−/PR−/HER2+ subtype, as well as its negative impact on survival following distant metastasis (12-15). Therefore, further in-depth research into this type of BC is still warranted.
In conventional survival analyses, such as the Kaplan-Meier method or the Cox proportional hazards model, cancer-specific mortality (CSM) is considered the outcome of interest, while deaths from other causes are censored. However, this approach ignores the competing risk posed by non-cancer deaths, which can introduce bias when estimating the cumulative incidence of CSM. In clinical practice, BC patients frequently die due to a variety of causes, and these causes compete with one another. The likelihood of CSM may be impacted by non-cancer-specific deaths, and there may be a chance that they will compete with one another. Because conventional survival analysis techniques do not account for the influence of competing risk factors on patient death, they may overestimate the risk of CSM (16). To analyze survival data with competing risks, Austin et al. (17) proposed a method called the competing risks model (CRM), which took into account the effect of other risk factors on CSM and was more consistent with clinical practice to study patient outcomes. CSM from conventional survival analysis was higher than that (18) from CRM.
Due to BC’s biological heterogeneity, the current staging systems based on Roman Numeral and tumor-node-metastasis (TNM), even integrating molecular subtyping information, are not sufficient to fully address clinical needs. In addition, other factors, including marital status, age, adjuvant therapy, and race, may affect prognosis. It has been demonstrated that the presentation of a nomogram as a prediction model outperforms traditional staging systems and shows significant strengths in the prediction of tumor recurrence, prognosis, and outcome (19).
Currently, nomograms have been developed for most BC subtypes, but because the ER−/PR−/HER2+ subtype accounts for the smallest proportion, there have been few studies on this subtype. The current study investigated the survival prognosis and clinical characteristics of patients with ER−/PR−/HER2+ BC via the Surveillance, Epidemiology, and End Results (SEER) database, and established a nomogram to accurately predict the BC-specific survival (BCSS) of patients with ER−/PR−/HER2+ BC based on independent prognostic factors. We present this article in accordance with the TRIPOD reporting checklist (available at https://tcr.amegroups.com/article/view/10.21037/tcr-2025-aw-2562/rc).
Methods
Patient selection and data acquisition
The National Cancer Institute created SEER, a database that gathers cancer patients’ information. The current study utilized SEER*Stat to acquire data on BC patients diagnosed in the period from 2010 to 2015. These cases were analyzed to include female patients with ER−/PR−/HER2+ BC and exclude patients with the following conditions: race, tumor stage, axillary lymph node (LN) positive grade I–II, local LN biopsy, local LN biopsy positive, bone/lung metastasis, tumor size unknown, no surgery or radiotherapy/status of surgery or radiotherapy unknown, survival time <1 month, unmarried or cohabiting, patients with non-primary tumors. Finally, 1,525 eligible subjects were selected for analysis (Figure 1). Because the SEER database’s information is accessible to the general public, the current study did not require patients’ informed consent or ethical approval. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments.
Variable collection
The SEER database’s demographic information incorporated year of diagnosis, age recode with single ages and 100+, marital status at diagnosis, and race recode (White, Black, American Indian/Alaska Native, Asian/Pacific Islander). Clinical indicators included primary site-labeled, International Classification of Diseases for Oncology, 3rd Edition (ICD-O-3) histological behavior, grade (through 2017), laterality, combined summary stage (2004+), derived American Joint Committee on Cancer (AJCC) stage group, 7th ed [2010–2015], derived AJCC M, 7th ed [2010–2015], derived AJCC N, 7th ed [2010–2015], derived AJCC T, 7th ed [2010–2015], RX Summ-Surg Prim Site (1998+), RX Summ-Surg/Rad Seq, radiation recode, chemotherapy recode, regional nodes positive (1988+), regional nodes examined (1988+), SEER Combined Mets at DX-lung (2010+), SEER Combined Mets at DX-liver (2010+), SEER Combined Mets at DX-bone (2010+), sequence number, first malignant primary indicator, CS tumor size [2004–2015], Total number of in situ/malignant tumors for patients.
Outcome
The current study selected BC-specific death (BCSD) as the outcome of interest, and deaths from other causes were considered as competing risk events. Survival time is defined as the period from diagnosis to death from BC. To facilitate subsequent study and analysis, the current study renamed the included indicators and assigned values to them (Table 1).
Table 1
| Factors | SEER | Assignment |
|---|---|---|
| Marriage | Marital status at diagnosis | 1: married (including common law) |
| 2: others: divorced, separated, single (never married); widowed | ||
| Age | Age recode with single ages and 100+ | 1: ≤40 years |
| 2: 41–65 years | ||
| 3: ≥66 years | ||
| Race | Race recode | 1: White |
| 2: Black | ||
| 3: Asian or Pacific Islander, American Indian/Alaska Native | ||
| Year | Year of diagnosis | 1: 2010–2013 |
| 2: 2014–2015 | ||
| Site | Primary site-labeled | 1: inner (C50.2-upper-inner quadrant of breast, C50.3-lower-inner quadrant of breast) |
| 2: outer (C50.4-upper-outer quadrant of breast, C50.5-lower-outer quadrant of breast) | ||
| 3: central (C50.0-nipple, C50.1-central portion of breast) | ||
| 4: others (C50.6-axillary tail of breast, C50.8-overlapping lesion of breast) | ||
| Histological behavior | ICD-O-3 histological behavior | 1: infiltrating duct carcinoma (8500/3: infiltrating duct carcinoma, NOS, 8521/3: infiltrating ductular carcinoma, 8522/3: infiltrating duct and lobular carcinoma, 8523/3: infiltrating duct mixed with other types of carcinoma) |
| 2: others (8500/3: infiltrating duct carcinoma, NOS, 8140/3: adenocarcinoma, NOS, 8230/3: solid carcinoma, NOS, 8401/3: apocrine adenocarcinoma, 8480/3: mucinous adenocarcinoma, 8501/3: comedocarcinoma, NOS, 8503/3: intraductal papillary adenocarcinoma with invasion, 8507/3: ductal carcinoma, micropapillary, 8510/3: medullary carcinoma, NOS, 8520/3: lobular carcinoma, NOS, 8524/3: infiltrating lobular mixed with other types of carcinoma, 8530/3: inflammatory carcinoma, 8541/3: Paget disease and infiltrating ductal carcinoma of breast, 8543/3: Paget disease and intraductal carcinoma, 8575/3: metaplastic carcinoma, NOS) | ||
| Grade | Grade (through 2017) | 1: well differentiated; grade I |
| 2: moderately differentiated; grade II | ||
| 3: poorly differentiated; grade III, undifferentiated; anaplastic; grade IV | ||
| Laterality | Laterality | 1: left, origin of primary |
| 2: right, origin of primary | ||
| Stage | Derived AJCC stage group, 7th ed [2010–2015] | 1: 0, IA, IB |
| 2: IIA, IIB | ||
| 3: IIIA, IIIB, IIIC, IV | ||
| T | Derived AJCC T, 7th ed [2010–2015] | 1: Tis, T1a, T1b, T1c, T1mic, T1NOS |
| 2: T2 | ||
| 3: T3, T4a, T4b, T4d | ||
| N | Derived AJCC N, 7th ed [2010–2015] | 1: N0, NO(i+), NO(i−), N0(mol−) |
| 2: N1, N1c, N1mi, N1NOS, N2, N2a, N2b, N2NOS, N3, N3a, N3b, N3NOS | ||
| M | Derived AJCC M, 7th ed [2010–2015] | 1: M0 |
| 2: M1 | ||
| Surgery | RX Summ-Surg Prim Site (1998+) | 1: 0.8–0.90 (unoperated or unknown) |
| 2: 20–30 (partial mastectomy) | ||
| 3: 40–75 (total mastectomy) | ||
| Radiation | RX Summ-Surg/Rad Seq | 1: intraoperative radiation, radiation after surgery, radiation before and after surgery, radiation prior to surgery |
| 2: no radiation and/or cancer-directed surgery | ||
| Rad_recode | Radiation recode | 1: beam radiation, combination of beam with implants or isotopes, radiation, NOS method or source not specified, radioactive implants (includes brachytherapy) (1988+) |
| 2: none/unknown, refused (1988+) | ||
| Chemotherapy | Chemotherapy recode | 1: yes |
| 2: no/unknown | ||
| Nodes | Regional nodes examined (1988+) | 1: 1–5 |
| 2: 6–46 | ||
| Nodes_Pos | Regional nodes positive (1988+) | 1: 0 |
| 2: 1–35 | ||
| DX_bone | SEER Combined Mets at DX-bone (2010+) | 1: yes |
| 2: no | ||
| DX_liver | SEER Combined Mets at DX-liver (2010+) | 1: yes |
| 2: no | ||
| DX_lung | SEER Combined Mets at DX-lung (2010+) | 1: yes |
| 2: no | ||
| Sequence | Sequence number | 1: one primary only |
| 2: 1st of 2 or more primaries, 2nd of 2 or more primaries, 3rd of 3 or more primaries, 4th of 4 or more primaries, 5th of 5 or more primaries | ||
| Size | CS tumor size (2004–2015) | 1: 1–10.991 mm |
| 2: 11–20.992 mm | ||
| 3: 21–30.993 mm | ||
| 4: >31.995 mm | ||
| Pri_indicator | First malignant primary indicator | 1: yes |
| 2: no | ||
| Tumors | Total number of in situ/malignant tumors for patient | 1: 1 |
| 2: >1 | ||
| Survival | Survival months | Continuous variable, no assignment required |
AJCC, American Joint Committee on Cancer; ICD-O-3, International Classification of Diseases for Oncology, 3rd Edition; NOS, not otherwise specified; SEER, Surveillance, Epidemiology, and End Results; TNM, tumor-node-metastasis.
Statistical analysis
R software (4.3.1) was employed for statistical analysis. Frequencies and percentages were utilized to express categorical variables, and the Chi-squared test was conducted to compare differences between groups. Mean ± standard deviation (SD) or interquartile range (IQR) was employed to express continuous variables. The 1,525 patients with ER−/PR−/HER2+ BC were randomized into the training (n=1,067) and validation (n=458) sets in 7:3. Multicollinearity among the candidate variables was assessed using the variance inflation factor (VIF). Variables with a VIF <10 were included in the model to ensure the absence of severe multicollinearity. Significant prognostic factors were initially identified through univariable analysis using the Fine-Gray test. Variables with a P value <0.05 were subsequently entered into a multivariable analysis. Then, a nomogram was constructed using R software to predict 1-, 2-, and 3-year cancer-specific survival (CSS) for patients with ER−/PR−/HER2+ BC. The receiver operating characteristic (ROC) curve and area under the curve (AUC)/concordance index (C-index) were adopted to quantitatively assess the model’s efficacy and predictive power. Calibration curves were utilized to validate the model and to assess the nomogram’s performance. In addition, the performance of CRM and conventional survival analysis was compared regarding 1-, 2-, and 3-year CSS. P<0.05 indicated statistical significance.
Results
Clinicopathologic and demographic characteristics
This study incorporated 1,525 ER−/PR−/HER2+ BC patients and randomized them in 7:3 into a training set (n=1,067) and a validation set (n=458). The demographic and clinical characteristics of all patients are illustrated in Table 2. Most patients had an age of 40 to 70 years, representing 75% of all patients. Among the histologic types, infiltrating duct carcinoma (IDC) was the dominant one, accounting for 94.6% of all patients. Among the grades (through 2017), patients with grade III/IV accounted for the largest proportion (72.9%). Tis/T1 (49.8%), stage N0 (60.7%), and M0 (97.2%) accounted for the largest proportion in the AJCC T/N/M staging systems, respectively. The tumor size was predominant in the 10–20 mm group, making up 28.3% of all patients. The total number of deaths was 203, of which 73 were non-specific, accounting for 36.0% of all deaths (Table 2).
Table 2
| Factors | Define | All (n=1,525), n (%) | Training (n=1,067), n (%) | Validation (n=458), n (%) | P† |
|---|---|---|---|---|---|
| Marriage | Married | 944 (61.9) | 669 (62.7) | 275 (60.0) | 0.35 |
| Others | 581 (38.1) | 398 (37.3) | 183 (40.0) | ||
| Age (years) | <40 | 136 (8.9) | 98 (9.2) | 38 (8.3) | 0.84 |
| 40–70 | 1,144 (75.0) | 799 (74.9) | 345 (75.3) | ||
| >70 | 245 (16.1) | 170 (15.9) | 75 (16.4) | ||
| Race | White | 938 (61.5) | 654 (61.3) | 284 (62.0) | 0.11 |
| Black | 232 (15.2) | 152 (14.2) | 80 (17.5) | ||
| Other | 355 (23.3) | 261 (24.5) | 94 (20.5) | ||
| Site | Inner | 300 (19.7) | 209 (19.6) | 91 (19.9) | 0.99 |
| Outer | 728 (47.7) | 509 (47.7) | 219 (47.8) | ||
| Central | 91 (6.0) | 65 (6.1) | 26 (5.7) | ||
| Others | 406 (26.6) | 284 (26.6) | 122 (26.6) | ||
| Histological behavior | IDC | 1,443 (94.6) | 1,014 (95.0) | 429 (93.7) | 0.33 |
| Others | 82 (5.4) | 53 (5.0) | 29 (6.3) | ||
| Grade | I | 22 (1.4) | 15 (1.4) | 7 (1.5) | 0.92 |
| II | 391 (25.6) | 271 (25.4) | 120 (26.2) | ||
| III/IV | 1,112 (72.9) | 781 (73.2) | 331 (72.3) | ||
| Laterality | Left | 791 (51.9) | 567 (53.1) | 224 (48.9) | 0.14 |
| Right | 734 (48.1) | 500 (46.9) | 234 (51.1) | ||
| Stage | 0/I | 615 (40.3) | 427 (40.0) | 188 (41.0) | 0.90 |
| II | 635 (41.6) | 448 (42.0) | 187 (40.8) | ||
| III/IV | 275 (18.0) | 192 (18.0) | 83 (18.1) | ||
| T stage | Tis/T1 | 759 (49.8) | 523 (49.0) | 236 (51.5) | 0.55 |
| T2 | 591 (38.8) | 423 (39.6) | 168 (36.7) | ||
| T3/T4 | 175 (11.5) | 121 (11.3) | 54 (11.8) | ||
| N stage | N0 | 925 (60.7) | 644 (60.4) | 281 (61.4) | 0.75 |
| N1/N2/N3 | 600 (39.3) | 423 (39.6) | 177 (38.6) | ||
| M stage | M0 | 1,482 (97.2) | 1,039 (97.4) | 443 (96.7) | 0.59 |
| M1 | 43 (2.8) | 28 (2.6) | 15 (3.3) | ||
| Surgery | No/unknown | 15 (1.0) | 11 (1.0) | 4 (0.9) | 0.92 |
| Partial | 767 (50.3) | 539 (50.5) | 228 (49.8) | ||
| Total | 743 (48.7) | 517 (48.5) | 226 (49.3) | ||
| Radiation | Yes | 827 (54.2) | 587 (55.0) | 240 (52.4) | 0.37 |
| No | 698 (45.8) | 480 (45.0) | 218 (47.6) | ||
| Radiation recode | Yes | 828 (54.3) | 588 (55.1) | 240 (52.4) | 0.35 |
| No | 697 (45.7) | 479 (44.9) | 218 (47.6) | ||
| Chemotherapy | Yes | 1,251 (82.0) | 869 (81.4) | 382 (83.4) | 0.40 |
| No | 274 (18.0) | 198 (18.6) | 76 (16.6) | ||
| Nodes | 1–5 | 956 (62.7) | 666 (62.4) | 290 (63.3) | 0.78 |
| 6–46 | 569 (37.3) | 401 (37.6) | 168 (36.7) | ||
| Nodes positive | 0 | 1,008 (66.1) | 697 (65.3) | 311 (67.9) | 0.35 |
| 1–35 | 517 (33.9) | 370 (34.7) | 147 (32.1) | ||
| Bone metastasis | Yes | 16 (1.0) | 10 (0.9) | 6 (1.3) | 0.70 |
| No | 1,509 (99.0) | 1,057 (99.1) | 452 (98.7) | ||
| Liver metastasis | Yes | 17 (1.1) | 12 (1.1) | 5 (1.1) | >0.99 |
| No | 1,508 (98.9) | 1,055 (98.9) | 453 (98.9) | ||
| Lung metastasis | Yes | 10 (0.7) | 7 (0.7) | 3 (0.7) | >0.99 |
| No | 1,515 (99.3) | 1,060 (99.3) | 455 (99.3) | ||
| Sequence | First | 1,193 (78.2) | 829 (77.7) | 364 (79.5) | 0.48 |
| Others | 332 (21.8) | 238 (22.3) | 94 (20.5) | ||
| Size (mm) | 1–10 | 330 (21.6) | 241 (22.6) | 89 (19.4) | 0.07 |
| 11–20 | 431 (28.3) | 284 (26.6) | 147 (32.1) | ||
| 21–30 | 361 (23.7) | 264 (24.7) | 97 (21.2) | ||
| ≥31 | 403 (26.4) | 278 (26.1) | 125 (27.3) | ||
| Primary indicator | Yes | 1,352 (88.7) | 946 (88.7) | 406 (88.6) | >0.99 |
| No | 173 (11.3) | 121 (11.3) | 52 (11.4) | ||
| Tumors | 1 | 1,212 (79.5) | 847 (79.4) | 365 (79.7) | 0.94 |
| >1 | 313 (20.5) | 220 (20.6) | 93 (20.3) |
†, represents the comparison between the training and validation groups. ER, estrogen receptor; HER2, human epidermal growth factor receptor 2; IDC, infiltrating duct carcinoma; PR, progesterone receptor; TNM, tumor-node-metastasis.
Nomogram based on CRM
Univariable analysis revealed that marital status at diagnosis, age, grade, stage, TNM stage, surgical method, number of LN biopsy, number of positive LN biopsy, tumor size, lung metastasis, and liver metastasis were independent risk factors for specific death in ER−/PR−/HER2+ BC patients. A multivariable analysis containing these factors indicated that age [<40 years as reference, hazard ratio =0.532; 95% confidence interval (CI): 0.284–0.997 for 40–70 years, hazard ratio =0.909; 95% CI: 0.409–2.018 for >70 years], stage (stage I as reference, hazard ratio =2.001; 95% CI: 0.758–5.279 for stage II, hazard ratio =4.938; 95% CI: 1.467–16.629 for stage III/IV), lung metastasis (with metastasis as reference, hazard ratio =0.198; 95% CI: 0.049–0.8 for no metastasis), and tumor size (<10 mm as reference, hazard ratio =1.963; 95% CI: 0.732–5.264 for 10–30 mm, hazard ratio =3.255; 95% CI: 1.015–10.441 for >31 mm) were predictors (Table 3). According to the results, a nomogram based on a CRM was created to predict 1-, 2-, and 3-year CSS for ER−/PR−/HER2+ BC patients (Figure 2).
Table 3
| Factors | Define | Univariable analysis | Multivariable analysis | |||
|---|---|---|---|---|---|---|
| HR (95% CI) | P | HR (95% CI) | P | |||
| Marriage | Married | REF | NA | REF | NA | |
| Others | 1.51 (1.01–2.26) | 0.044 | 1.194 (0.784–1.819) | 0.41 | ||
| Age (years) | <40 | REF | NA | REF | NA | |
| 40–70 | 0.49 (0.27–0.86) | 0.01 | 0.532 (0.284–0.997) | 0.049 | ||
| >70 | 0.73 (0.36–1.46) | 0.38 | 0.909 (0.409–2.018) | 0.81 | ||
| Race | White | REF | NA | REF | NA | |
| Black | 1.27 (0.74–2.16) | 0.39 | NA | NA | ||
| Other | 0.76 (0.45–1.28) | 0.30 | NA | NA | ||
| Site | Inner | REF | NA | REF | NA | |
| Central | 0.7 (0.24–2.07) | 0.52 | NA | NA | ||
| Others | 1.09 (0.6–1.98) | 0.79 | NA | NA | ||
| Outer | 1.04 (0.6–1.79) | 0.90 | NA | NA | ||
| Histological behavior | IDC | REF | NA | REF | NA | |
| Others | 1.13 (0.46–2.78) | 0.80 | NA | NA | ||
| Grade | I/II | REF | NA | REF | NA | |
| III/IV | 1.94 (1.11–3.36) | 0.01 | 1.404 (0.759–2.598) | 0.28 | ||
| Laterality | Left | REF | NA | REF | NA | |
| Right | 0.85 (0.56–1.28) | 0.43 | NA | NA | ||
| Stage | I | REF | NA | REF | NA | |
| II | 3.05 (1.6–5.81) | 0.0007 | 2.001 (0.758–5.279) | 0.16 | ||
| III/IV | 9.22 (4.89–17.38) | <0.001 | 4.938 (1.467–16.629) | 0.01 | ||
| T stage | Tis/T1 | REF | NA | REF | NA | |
| T2 | 2.48 (1.52–4.05) | 0.0003 | 0.742 (0.348–1.581) | 0.44 | ||
| T3/T4 | 4.71 (2.68–8.29) | <0.001 | 0.473 (0.176–1.271) | 0.14 | ||
| N stage | N0 | REF | NA | REF | NA | |
| N1/N2/N3 | 3.98 (2.54–6.22) | <0.001 | 1.192 (0.38–3.742) | 0.76 | ||
| M stage | M0 | REF | NA | REF | NA | |
| M1 | 5.39 (2.68–10.83) | <0.001 | 0.956 (0.292–3.136) | 0.94 | ||
| Surgery | No/unknown | REF | NA | REF | NA | |
| Partial | 0.2 (0.06–0.67) | 0.009 | 1.016 (0.288–3.583) | 0.98 | ||
| Total | 0.35 (0.1–1.16) | 0.08 | 1.18 (0.336–4.145) | 0.80 | ||
| Radiation | Yes | REF | NA | REF | NA | |
| No | 1.35 (0.9–2.03) | 0.14 | NA | NA | ||
| Radiation recode | Yes | REF | NA | REF | NA | |
| No | 1.36 (0.91–2.04) | 0.14 | NA | NA | ||
| Chemotherapy | Yes | REF | NA | REF | NA | |
| No | 1.18 (0.72–1.95) | 0.51 | NA | NA | ||
| Nodes | 1–5 | REF | NA | REF | NA | |
| 6–46 | 2.11 (1.41–3.17) | <0.001 | 0.686 (0.396–1.189) | 0.18 | ||
| Nodes positive | 0 | REF | NA | REF | NA | |
| 1–35 | 3.75 (2.45–5.76) | <0.001 | 1.633 (0.605–4.409) | 0.33 | ||
| Bone metastasis | Yes | REF | NA | REF | NA | |
| No | 0.39 (0.09–1.67) | 0.21 | NA | NA | ||
| Liver metastasis | Yes | REF | NA | REF | NA | |
| No | 0.15 (0.06–0.4) | <0.001 | 0.794 (0.197–3.199) | 0.75 | ||
| Lung metastasis | Yes | REF | NA | REF | NA | |
| No | 0.08 (0.04–0.2) | <0.001 | 0.198 (0.049–0.8) | 0.02 | ||
| Sequence | First | REF | NA | REF | NA | |
| Others | 1.21 (0.77–1.91) | 0.40 | NA | NA | ||
| Size (mm) | <10 | REF | NA | REF | NA | |
| 10–30 | 3.12 (1.33–7.32) | 0.008 | 1.963 (0.732–5.264) | 0.18 | ||
| >31 | 7.27 (3.11–16.97) | <0.001 | 3.255 (1.015–10.441) | 0.047 | ||
| Primary indicator | Yes | REF | NA | REF | NA | |
| No | 1.04 (0.56–1.94) | 0.90 | NA | NA | ||
| Tumors | 1 | REF | NA | REF | NA | |
| >1 | 1.2 (0.75–1.92) | 0.44 | NA | NA | ||
CI, confidence interval; CSS, cancer-specific survival; HR, hazard ratio; IDC, infiltrating duct carcinoma; NA, not applicable; REF, reference; TNM, tumor-node-metastasis.
Validation of nomogram
The current study employed C-index, ROC curve, and AUC to evaluate the predictive accuracy of nomograms, and their calibration performance was verified by calibration curves. In the training set, the C-index of nomograms was 0.780 (SE: 0.007), and the AUC for 1-, 2-, and 3-year CSS prediction was 0.825 (95% CI: 0.729–0.921), 0.794 (95% CI: 0.716–0.873), and 0.781 (95% CI: 0.724–0.838), respectively (Figure 3), demonstrating a good predictive performance of the model. In the validation set, the C-index was 0.742 (SE: 0.015), and the AUC for 1-, 2-, and 3-year CSS prediction was 0.864 (95% CI: 0.71–1), 0.793 (95% CI: 0.668–0.917), and 0.772 (95% CI: 0.667–0.876), respectively (Figure 3), indicating high predictive value and high reliability of the model in internal validation. The calibration plots indicated good consistency between the nomogram’s prediction and the actual observation in the training and validation sets. Calibration results are illustrated in Figure 4.
Discussion
In this study, a SEER-based CRM was applied to introduce nomograms to predict CSS in ER−/PR−/HER2+ BC patients at 1, 2, and 3 years, and its performance was evaluated via a validation set. The results indicated that the nomogram exhibited high reliability and significant predictive value.
In the current study, it was found that only 130 of the 1,525 ER−/PR−/HER2+ patients died of BC, and 1,395 survived or died of other diseases. Therefore, CRM was applied to predict prognostic factors for ER−/PR−/HER2+ BC.
Nomograms are a common tool for evaluating oncologic and medical outcomes. Its ability to produce a single numerical probability of a clinical event by integrating various prognostic and determinant variables satisfies the desire for integrated clinical and biological models and promotes the development of personalized medicine (19). The nomogram based on multivariable analysis has good accuracy as a comprehensive statistical model. Previous studies revealed that nomograms had advantages in the evaluation of cancer risks, selection of therapies and drugs, and prediction of survival outcomes for various cancers (19,20).
The current study indicated the important role of age in affecting the prognosis of ER−/PR−/HER2+ BC. Since the age of 40 years is a reasonable cutoff for defining “young” (21,22), this study set three age groups: <40 years, 40–70 years, and >70 years. It was found that BC patients aged >70 years had the highest risk of death. Hershman et al. also discovered that whether the endpoint was overall survival (OS) or BCSD, the prognosis was worse in older BC patients (23). The reason might be that older patients were unable to tolerate intense standard chemotherapy and appeared to be more susceptible to cardiotoxicity from chemotherapy than younger patients (23). Additionally, the incidence of other chronic diseases (e.g., diabetes mellitus) that affected survival was higher in older patients. Nomograms in the current study demonstrated that BC patients aged <40 years exhibited a higher risk of death than those aged 40–70 years. The finding was similar to that of Chen et al. (24), LWingo P et al. (25), and Yancik et al. (26). The reason might be related to the tissue type of BC patients aged <40 years. Azim et al. revealed that regardless of subtype, grade, and stage, the expression of RANK ligand, c-kit, mammary stem and luminal progenitor cells, and BRCA1 mutation signatures were higher in young patients (27). Increasing evidence also suggests that differences in the breast matrix of young patients, as well as changes during pregnancy and lactation, may contribute to differences in tumor biology that emerge later (28).
The prognosis of metastatic breast cancer (MBC) is known to be worse than that of non-metastatic BC. Common sites for distant metastasis of BC include bone, liver, lung, and brain (29). The current study indicated that liver metastasis was most common in patients with ER−/PR−/HER2+ BC, followed by bone metastasis and lung metastasis. Kennecke et al. (14) also found that compared to luminal A tumors, HER2-rich breast tumors had a higher chance of developing liver metastases. Similarly, Leone et al. (30) utilized the SEER database to reveal that liver metastasis was more common in the HR−/HER2+ subtype. Some studies suggest that CXCR4 is a chemokine receptor that is upregulated by HER2 activation and has been proposed to be involved in promoting tumor cell invasion into the liver (31). However, other studies have reported that liver metastasis is not linked to BC subtypes (32). The current study further analyzed the effect of bone, liver, and lung metastasis on prognosis and indicated that only lung metastasis had an obvious effect on the prognosis of HR−/HER2+ BC patients. However, Gerratana et al. (29) reported that compared to BC patients with metastases of bone (44.4 months), liver (36.7 months), or brain (7.35 months) as the first distant metastasis, those with lung metastasis exhibited the best survival outcomes (58.5 months). The reason for inconsistent findings might be that Gerratana et al. did not conduct separate studies of HER2-overexpressing BC and the sample size in our study was small.
Tumor size has long been considered a significant prognostic factor for BC, and it has been reviewed by various investigators via data from the Medline database (33), SEER database (34), and Norway Cancer Registry (35). Therefore, the components of the AJCC staging system include tumor size, regional LN status, and distant metastasis (36). The current study revealed that the tumor size of ER−/PR−/HER2+ BC was also an important factor in patient prognosis. When the tumor size was below 10 mm, the risk of death was lowest, while the risk of death was highest when it was above 31 mm. This was comparable to what Rosenberg et al. discovered (34,37). However, Yu et al. found that the correlation of tumor size with breast cancer-specific mortality (BCSM) was segmented in LN-negative diseases. Using tumors of 21 to 30 mm as a reference, the hazard ratio of BCSM rose as tumor size increased until it peaked at a tumor size of 41 to 50 mm, after which the increase in tumor size was unexpectedly correlated with a decrease in the risk ratio, with the lowest point at a tumor size of 61 to 80 mm (38). Furthermore, Wo et al. suggested that small tumors with four LN positives might predict a higher BCSM compared to larger tumors (39).
BC in the SEER database is classified primarily by the differentiation degree of tumor cells. Studies have shown that the classification of BC was important for the prognosis of patients (40,41). The current study illustrated that among all prognostic factors, tumor grade had the greatest impact on the prognosis of HER2-overexpressing BC patients, and the nomogram score of BC patients with poorly differentiated and undifferentiated tumor grades (III/IV) reached 100. Poorly differentiated or undifferentiated BC cells have the characteristics of obvious karyotype mutations and structural disorders. From a biological point of view, such cells tend to have stronger invasiveness and metastasis, and may also have stronger immune escape capabilities. Therefore, such tumors are more malignant and have strong drug resistance and reduced sensitivity to chemotherapy and radiotherapy, greatly affecting patient prognosis.
However, there were limitations to the current study. Firstly, it was a retrospective study, and selection bias was inevitable. Secondly, due to the limitations of public databases, it is impossible to fully study the impact of endocrine therapy, specific radiotherapy and chemotherapy regimens, lifestyle, and family history on the outcomes. For instance, the application of pertuzumab from 2010 to 2015 might have affected the prognosis of patients. Thirdly, our model is designed specifically for treatment-naïve, de novo ER−/PR−/HER2+ BC. It may not be applicable to patients with recurrent or metastatic disease. Finally, as no external validation set was included, the application of nomograms to risk groups in Chinese patient cohorts was not validated. To increase the precision and integrity of nomograms, prospective data can be collected subsequently, and more prospective validation can be included in future studies.
Conclusions
A prediction model was successfully established using the SEER database to accurately predict CSS in patients with ER−/PR−/HER2+ BC based on independent prognostic factors. This model exhibited good performance in both the training and validation sets. It can assist medical professionals in creating customized treatment and follow-up plans.
Acknowledgments
None.
Footnote
Reporting Checklist: The authors have completed the TRIPOD reporting checklist. Available at https://tcr.amegroups.com/article/view/10.21037/tcr-2025-aw-2562/rc
Peer Review File: Available at https://tcr.amegroups.com/article/view/10.21037/tcr-2025-aw-2562/prf
Funding: None.
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://tcr.amegroups.com/article/view/10.21037/tcr-2025-aw-2562/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Siegel RL, Miller KD, Fuchs HE, et al. Cancer Statistics, 2021. CA Cancer J Clin 2021;71:7-33. [Crossref] [PubMed]
- Bartlett JM, Brookes CL, Robson T, et al. Estrogen receptor and progesterone receptor as predictive biomarkers of response to endocrine therapy: a prospectively powered pathology study in the Tamoxifen and Exemestane Adjuvant Multinational trial. J Clin Oncol 2011;29:1531-8. [Crossref] [PubMed]
- Kurozumi S, Matsumoto H, Hayashi Y, et al. Power of PgR expression as a prognostic factor for ER-positive/HER2-negative breast cancer patients at intermediate risk classified by the Ki67 labeling index. BMC Cancer 2017;17:354. [Crossref] [PubMed]
- Slamon DJ, Clark GM, Wong SG, et al. Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science 1987;235:177-82. [Crossref] [PubMed]
- Press MF, Bernstein L, Thomas PA, et al. HER-2/neu gene amplification characterized by fluorescence in situ hybridization: poor prognosis in node-negative breast carcinomas. J Clin Oncol 1997;15:2894-904. [Crossref] [PubMed]
- Jensen EV, Jordan VC. The estrogen receptor: a model for molecular medicine. Clin Cancer Res 2003;9:1980-9.
- Ueno T. Biomarkers of neoadjuvant/adjuvant endocrine therapy for ER-positive/HER2-negative breast cancer. Chin Clin Oncol 2020;9:35. [Crossref] [PubMed]
- Gianni L, Eiermann W, Semiglazov V, et al. Neoadjuvant chemotherapy with trastuzumab followed by adjuvant trastuzumab versus neoadjuvant chemotherapy alone, in patients with HER2-positive locally advanced breast cancer (the NOAH trial): a randomised controlled superiority trial with a parallel HER2-negative cohort. Lancet 2010;375:377-84. [Crossref] [PubMed]
- Ross JS, Slodkowska EA, Symmans WF, et al. The HER-2 receptor and breast cancer: ten years of targeted anti-HER-2 therapy and personalized medicine. Oncologist 2009;14:320-68. [Crossref] [PubMed]
- Curigliano G, Burstein HJ, Winer EP, et al. De-escalating and escalating treatments for early-stage breast cancer: the St. Gallen International Expert Consensus Conference on the Primary Therapy of Early Breast Cancer 2017. Ann Oncol 2017;28:1700-12.
- Giaquinto AN, Sung H, Newman LA, et al. Breast cancer statistics 2024. CA Cancer J Clin 2024;74:477-95. [Crossref] [PubMed]
- Voduc KD, Cheang MC, Tyldesley S, et al. Breast cancer subtypes and the risk of local and regional relapse. J Clin Oncol 2010;28:1684-91. [Crossref] [PubMed]
- Lobbezoo DJ, van Kampen RJ, Voogd AC, 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:507-14. [Crossref] [PubMed]
- Kennecke H, Yerushalmi R, Woods R, et al. Metastatic behavior of breast cancer subtypes. J Clin Oncol 2010;28:3271-7. [Crossref] [PubMed]
- Kast K, Link T, Friedrich K, et al. Impact of breast cancer subtypes and patterns of metastasis on outcome. Breast Cancer Res Treat 2015;150:621-9. [Crossref] [PubMed]
- Meng X, Hao F, Ju Z, et al. Conditional survival nomogram predicting real-time prognosis of locally advanced breast cancer: Analysis of population-based cohort with external validation. Front Public Health 2022;10:953992. [Crossref] [PubMed]
- Austin PC, Lee DS, Fine JP. Introduction to the Analysis of Survival Data in the Presence of Competing Risks. Circulation 2016;133:601-9. [Crossref] [PubMed]
- Li Z, Shi Y, Wu L, et al. Establishment and verification of a nomogram to predict tumor-specific mortality risk in triple-negative breast cancer: a competing risk model based on the SEER cohort study. Gland Surg 2022;11:1961-75. [Crossref] [PubMed]
- Balachandran VP, Gonen M, Smith JJ, et al. Nomograms in oncology: more than meets the eye. Lancet Oncol 2015;16:e173-80. [Crossref] [PubMed]
- Iasonos A, Schrag D, Raj GV, et al. How to build and interpret a nomogram for cancer prognosis. J Clin Oncol 2008;26:1364-70. [Crossref] [PubMed]
- Fu J, Wu L, Fu W, et al. How Young Is Too Young in Breast Cancer?-Young Breast Cancer Is Not a Unique Biological Subtype. Clin Breast Cancer 2018;18:e25-39. [Crossref] [PubMed]
- Zheng H, Ge C, Lin H, et al. The impact of age on outcomes of breast cancer in different hormone receptor and HER2 groups. PLoS One 2023;18:e0280474. [Crossref] [PubMed]
- Hershman DL, Till C, Shen S, et al. Association of Cardiovascular Risk Factors With Cardiac Events and Survival Outcomes Among Patients With Breast Cancer Enrolled in SWOG Clinical Trials. J Clin Oncol 2018;36:2710-7. [Crossref] [PubMed]
- Chen HL, Zhou MQ, Tian W, et al. Effect of Age on Breast Cancer Patient Prognoses: A Population-Based Study Using the SEER 18 Database. PLoS One 2016;11:e0165409. [Crossref] [PubMed]
- Wingo PA, Ries LA, Parker SL, et al. Long-term cancer patient survival in the United States. Cancer Epidemiol Biomarkers Prev 1998;7:271-82.
- Yancik R, Wesley MN, Ries LA, et al. Effect of age and comorbidity in postmenopausal breast cancer patients aged 55 years and older. JAMA 2001;285:885-92. [Crossref] [PubMed]
- Azim HA Jr, Partridge AH. Biology of breast cancer in young women. Breast Cancer Res 2014;16:427. [Crossref] [PubMed]
- Azim HA Jr, Brohée S, Peccatori FA, et al. Biology of breast cancer during pregnancy using genomic profiling. Endocr Relat Cancer 2014;21:545-54. [Crossref] [PubMed]
- Gerratana L, Fanotto V, Bonotto M, et al. Pattern of metastasis and outcome in patients with breast cancer. Clin Exp Metastasis 2015;32:125-33. [Crossref] [PubMed]
- Leone JP, Leone BA, Zwenger AO, et al. The prognostic significance of metastatic pattern in stage IV male breast cancer at initial diagnosis: a population-based study. Breast Cancer Res Treat 2021;187:237-44. [Crossref] [PubMed]
- Li YM, Pan Y, Wei Y, et al. Upregulation of CXCR4 is essential for HER2-mediated tumor metastasis. Cancer Cell 2004;6:459-69. [Crossref] [PubMed]
- Park HS, Kim S, Kim K, et al. Pattern of distant recurrence according to the molecular subtypes in Korean women with breast cancer. World J Surg Oncol 2012;10:4. [Crossref] [PubMed]
- Mirza AN, Mirza NQ, Vlastos G, et al. Prognostic factors in node-negative breast cancer: a review of studies with sample size more than 200 and follow-up more than 5 years. Ann Surg 2002;235:10-26. [Crossref] [PubMed]
- Zheng YZ, Wang L, Hu X, et al. Effect of tumor size on breast cancer-specific survival stratified by joint hormone receptor status in a SEER population-based study. Oncotarget 2015;6:22985-95. [Crossref] [PubMed]
- Johansson ALV, Trewin CB, Fredriksson I, et al. In modern times, how important are breast cancer stage, grade and receptor subtype for survival: a population-based cohort study. Breast Cancer Res 2021;23:17. [Crossref] [PubMed]
- Boffa DJ, Greene FL. Reacting to changes in staging designations in the 7th edition of the AJCC staging manual. Ann Surg Oncol 2011;18:1-3.
- Rosenberg J, Chia YL, Plevritis S. The effect of age, race, tumor size, tumor grade, and disease stage on invasive ductal breast cancer survival in the U.S. SEER database. Breast Cancer Res Treat 2005;89:47-54.
- Yu KD, Jiang YZ, Chen S, et al. Effect of large tumor size on cancer-specific mortality in node-negative breast cancer. Mayo Clin Proc 2012;87:1171-80. [Crossref] [PubMed]
- Wo JY, Chen K, Neville BA, et al. Effect of very small tumor size on cancer-specific mortality in node-positive breast cancer. J Clin Oncol 2011;29:2619-27. [Crossref] [PubMed]
- Zhong X, Lin Y, Zhang W, et al. Predicting diagnosis and survival of bone metastasis in breast cancer using machine learning. Sci Rep 2023;13:18301. [Crossref] [PubMed]
- Huang X, Luo Z, Liang W, et al. Survival Nomogram for Young Breast Cancer Patients Based on the SEER Database and an External Validation Cohort. Ann Surg Oncol 2022;29:5772-81. [Crossref] [PubMed]

