Baseline patient characteristics
Table I shows demographic and clinicopathological characteristics of patients in this study. The follow-up time ranges from 0 to 71 months, with a median time if 31 months. There were 4288 deaths or 2905 deaths caused by breast cancer during follow-up. The two groups have significantly different distribution of clinicopathological characteristics. Compared with the no PMRT group, tumors in the PMRT group were more aggressive. Specifically, PMRT group shows higher grade (III and IV, 55.4% vs. 41.0%, p < 0.001), lager tumor size (T2, 3, 4, 76.5% vs. 42.9%, p < 0.001), more lymph nodes (LN positive, 78.6% vs. 29.2%, p < 0.001), more absence of ER (26.2% vs. 21.6%, p < 0.001) and PR (37.4% vs. 32.1%, p < 0.001), HER2 overexpression (23.8% vs. 20.5%, p < 0.001) and younger age (15–59, 68.5% vs. 56.5%, p < 0.001). In addition, the PMRT group received more radical surgery (MRM, RM and ERM, 61.3% vs. 30.9%, p < 0.001) and more chemotherapy (87.2% vs. 44.4%, p < 0.001).
Comparison of survival between the PMRT and no PMRT groups
We performed Kaplan–Meier analysis to compare the survival between the PMRT and no PMRT groups. The PMRT group had significantly shorter survival in comparison of the no PMRT group both in BCSS (HR, 1.975, 95% CI, 1.835–2.126, p < 0.001) and OS (HR, 1.378, 95% CI, 1.293–1.468, p < 0.001) (Fig. 2). The five-year BCSS rates were 85.2% (95% CI, 84.3–86.1%) and 92.4% (95% CI, 92.0–92.8%) in the PMRT group and no PMRT group, respectively. Similarly, OS rates were 82.7% (95% CI, 81.7–83.6%) and 87.5% (95% CI, 87.0–88.0%) in the two groups, respectively.
In the univariate model, radiation, race, grade, tumor size, lymph node status, diagnosis year (2012–2013), age group (40–49,50–59,60–69 and 80–99), ER status, PR status, HER2 status, surgery, and chemotherapy were significantly associated with BCSS (Table II). In the multivariate model, radiation, race, grade, tumor size, lymph nodes, age group (70–79, 80–99), ER status, PR status, HER2 status, surgery and chemotherapy remained significantly related to BCSS after controlling the above factors. Compared with the no PMRT group, the PMRT group had longer BCSS (HR, 0.739, 95% CI, 0.679–0.805, p < 0.001). Several characteristics, including radiation, race, tumor size, lymph node status, diagnosis year (2014–2015), age group (40–49,70–79 and 80–99), ER status, PR status, HER2 status and surgery were significantly related to OS in the univariate model (Table III). And in the multivariate model, radiation, race, grade, tumor size, lymph nodes, age group (60–69, 70–79, 80–99), ER status, PR status, HER2 status, surgery and chemotherapy remained significantly related to OS. It is noteworthy, PMRT was a significant better prognostic predictor for OS (HR, 0.721, 95% CI, 0.670–0.777, p < 0.001).
Comparison of survival between the PMRT and no PMRT groups based on LN subgroups
We stratified all patients into three subgroups based on positive axillary lymph nodes. As presented in Fig. 3 and Fig. 4, the PMRT group showed improved BCSS and OS in the LN 1 to 3 subgroup(HR, 0.738, 95% CI, 0.639–0.853, p < 0.001 and HR, 0.684, 95% CI, 0.604–0.776, p < 0.001,respectively) and LN > = 4 subgroup (HR, 0.635, 95% CI, 0.563–0.717, p < 0.001 and HR, 0.623, 95% CI, 0.558–0.696, p < 0.001,respectively) except for LN 0 subgroup (HR, 1.146 95% CI, 0.925–1.419, p = 0.211 and HR, 1.035, 95% CI, 0.868–1.234 p = 0.704, respectively).
Comparison of survival between the PMRT and no PMRT groups in the LN 1–3 and T1-2 subgroup
In the univariate model, PMRT, grade, race, age group (80–99), ER status, PR status, HER2 status, surgery and chemotherapy were significantly related to BCSS (Table IV). After controlling the above factors, beam radiation, race (other race), grade, age group (80–99), ER status, PR status, HER2 status and surgery remained significantly related to BCSS in the multivariate model. The PMRT group had better BCSS (HR, 0.826, 95% CI, 0.688–0.992, p = 0.040) compared with the no PMRT group. Some characteristics, including radiation, race, grade, age group (60–69,70–79,80–99), ER status, PR status, HER2 status, surgery and chemotherapy were significantly related to OS in the univariate model, as is shown in Table IV. In the multivariate model, radiation, race (other race), grade, age group (70–79,80–99), ER status, PR status, HER2 status, surgery and chemotherapy remained significantly related to OS. Specially, PMRT was a significant better prognostic predictor for OS (HR, 0.751, 95% CI, 0.643–0.878, p < 0.001).
Subgroup analysis of patients with 1–3 lymph nodes and T1-2 based on T stage and molecular subtype
Since the effect on survival of PMRT may be different in different T stages and molecular subtypes, subgroup analyses of survival by stratifying T stage and molecular subtype were performed. When stratified by T stage, PMRT can significantly improve BCSS in T2 stage(HR, 0.773, 95% CI, 0.629–0.949, p = 0.014) and T4 stage (HR, 0.605, 95% CI, 0.415–0.882, p = 0.009), but not in T1 stage (HR, 0.841, 95% CI, 0.565–1.253, p = 0.395) and T3 stage (HR, 0.789, 95% CI, 0.583–1.069, p = 0.127) (Fig. 5). In terms of OS, PMRT can significantly improve survival in T2 (HR, 0.710, 95% CI, 0.594–0.848, p < 0.001),T3 (HR, 0.737, 95% CI, 0.561–0.969, p = 0.029) and T4 (HR, 0.550, 95% CI, 0.392–0.773, p = 0.001) stages except T1 stage(HR, 0.744, 95% CI, 0.535–1.033, p = 0.078) (Fig. 6).
As presented in Figs. 7 and 8, the PMRT group showed improved BCSS and OS in the HER2-/HR + subtype (HR, 0.722, 95% CI, 0.548–0.953, p = 0.021 and HR, 0.648, 95% CI, 0.518–0.811, p < 0.001,respectively), while there was no significant difference in the HER2+/HR+ (HR, 1.253, 95% CI, 0.673–2.331, p = 0.477 and HR, 0.971, 95% CI, 0.572–1.647, p = 0.912,respectively), HER2+/HR- (HR, 0.831, 95% CI, 0.415–1.662, p = 0.600 and HR, 0.651, 95% CI, 0.344–1.233, p = 0.188, respectively)and triple-negative (HR, 0.863, 95% CI, 0.646–1.153, p = 0.320 and HR, 0.843, 95% CI, 0.646–1.100, p = 0.209, respectively) subtypes.