As a relatively minimally invasive method with both high accuracy and sensitivity, SLNB has increasingly become the preferred method for assessing ALN status in patients diagnosed with breast cancer who have clinically negative ALNs. However not all patients with positive SLNB develop further lymph node metastases. In our study, non-SLN metastasis occurred in only 47.1% (280/594) of patients, which is roughly the same as in other studies [1, 27–30]. The AMAROS and OTOASOR trials showed that these patients had similar local recurrence rates, disease-free survival (DFS), and overall survival (OS) with axillary radiotherapy compared to ALND, with significantly fewer side effects than with ALND, and ultimately enhanced quality of life [11, 12]. Of course, it is important to note that approximately 50% of patients do not have non-SLN metastasis, and additional axillary RNI may also have adverse effects and increase the burden on the patient. Therefore the non-SLN metastasis model is useful to help identify the presence of non-SLN metastasis in individuals diagnosed with early breast cancer. Allowing for a 10% FNR, our new model can predict that 24.5% of SLN-positive patients do not develop further metastases. This subset of patients could avoid further ALND, and additional axillary RNI. This prediction allows patients to avoid superfluous treatment and mitigates the various complications associated with treatment.
Tevis et al. [31] found that the OncotypeDX recurrence score in clinically ALN-negative breast cancer was not predictive of lymph node burden and was not conducive to guiding decisions about the extent of axillary surgery. Currently, the most widely employed and efficacious approach involves modeling to predict non-SLN status, utilizing clinicopathological parameters, including those pertaining to the primary tumor and the SLN. In the present study, the univariate analysis revealed that 9 variables, namely tumor size, multifocality, lymphovascular invasion, ECE, ER status, number of negative SLNs, number of positive SLNs, size of the SLN metastasis, and metastatic SLN location in axilla, were significantly associated with non-SLN metastasis. Subsequent multivariate analysis demonstrated that all parameters, except for ER status, independently contributed to the risk of non-SLN metastasis. The pathological indicators of the primary tumor, including tumor size, multifocality, and lymphovascular invasion, exhibited robust predictive capabilities for non-SLN status, as validated by multiple studies [16, 28–30]. In contrast, tumor location, histological grade, ER status, and molecular subtype were rarely used as significant predictors of non-SLN metastasis, which is consistent with our results. The pathologic indicators of SLNs with significant independent predictive value in our study encompassed ECE, number of negative SLNs, number of positive SLNs, size of the SLN metastasis, and metastatic SLN location in axilla. The prevalence of ECE in SLN-positive patients has been documented to range from 30 to 37.6 percent [32–34]. In contrast, this rate was only 22.7% (135/594) in our study. Nevertheless, our research findings still indicate that ECE has a significant impact on non-SLN metastasis, which is consistent with the results of various prior studies [35–38], despite the contrary view of Wu et al [29].
Multiple studies have revealed that the number of negative SLNs and number of positive SLNs are the most common indicators applied to predict non-SLN metastasis [16, 18, 19, 39–42]. An escalation in the number of positive SLNs has been linked to a heightened risk of non-SLN metastasis, while a rise in the number of negative SLNs implies a diminished probability of non-SLN metastasis. The size of SLN metastasis includes macrometastasis, micrometastases, and ITC. Of those with SLN micrometastasis in our study, only 12.8% (5/39) had non-SLN metastases. The IBCSG 23 − 01 and AATRM trials have both indicated a 13% rate of metastasis in patients with SLN micrometastasis who underwent ALND to identify non-SLN [43, 44]. It is noteworthy that the definition of micrometastasis included ITC in IBCSG 23 − 01. Macrometastasis was found to be associated with a higher incidence of non-SLN metastasis in comparison to micrometastasis [45, 46].
In our new model, the metastatic SLN location is a novel and significant predictor of non-SLN metastasis. The positive SLN is categorized into lower and upper groups based on its relative location to the ICBN. Positive SLNs located both above and below the ICBN are also included in the upper group. Two studies by Li et al. found that MB-stained SLNs were predominantly situated under the ICBN; in 95% (56/59) of patients, positive SLNs were found to be located below the ICBN, and only 3 cases of positive SLN were located on both lower and upper sides [47, 48]. In contrast, our study showed that the lower group accounted for 59% (350/594) of cases. This may be because Li et al. used MB for SLN tracing, whereas we used MB in combination with ICG. Extensive research has demonstrated that the employment of MB alone for SLNB may result in a certain degree of FNR, and a fluorescent combined dye method can improve the detection rate of SLN [49–52]. The results suggested that positive SLN location above the ICBN was a high risk factor for non-SLN metastasis. Therefore, we should not ignore ALND when positive SLNs appear in the upper group. To our knowledge, this is the first study to include the metastatic SLN location in a non-SLN metastasis prediction model.
The MSKCC nomogram is a widely recognized model for predicting non-SLN metastasis. The model includes tumor size, tumor type and grade, number of positive SLNs, number of negative SLNs, lymphovascular invasion, multifocality, and ER status. Numerous studies have validated the MSKCC model, but the AUC varies from 0.6 to 0.8 because of differences in regions and patient populations [3, 4, 17, 20, 29, 53–59]. The SCH nomogram was the first model established using a population of Chinese individuals diagnosed with breast cancer, incorporating variables such as number of positive SLNs, number of negative SLNs, lymphovascular invasion, and size of the SLN metastasis [20], with an original AUC of 0.779. Upon validation using our database, the MSKCC and SCH models demonstrated AUC values of 0.715 and 0.755, suggesting a good discriminatory ability. However, Wu et al. [29] tested the MSKCC and SCH models in 236 Chinese individuals diagnosed with breast cancer and their AUCs were only 0.677 and 0.674. The reason for this discrepancy may be the large difference in the age distribution of our patients and those of Wu et al. In our study, 56.3% (271/481) of the patients were older than 50 years, compared to 37.7% (89/236) in theirs. These findings demonstrate that the efficacy of the models is contingent upon factors such as sample size, parameter definition, parameter inclusion, and ethnic disparities, thereby limiting their applicability to breast cancer in the Chinese population.
Our model was built based on 481 SLN-positive patients with an AUC of 0.830. Our validation cohort had an AUC of 0.785, suggesting favorable performance of the model. Comparative analysis with the conventional MSKCC and SCH models revealed superior diagnostic performance of our model. This can be attributed to our inclusion of a greater number of variables, which effectively mitigated the impact of variable distribution differences. Consequently, our nomogram demonstrated enhanced performance within the patient population. In order to avoid the damage caused by ALND in SLN-positive patients, our model was able to predict 3.3% (0%), 17.7% (< 5%), and 24.5% (< 10%) of patients who could be considered to forgo ALND with a guaranteed FNR. Notably, the MSKCC and SCH models were inferior to our model. However, the focus of the AUC is solely on the predictive accuracy of the model, which does not necessarily imply its clinical applicability. To address this limitation, DCA was employed to evaluate the clinical utility of the model by assessing its impact on clinical decision-making. By setting three risk thresholds based on the clinical non-SLN metastasis rate, we determined that our model yielded a higher net benefit compared to both the MSKCC and SCH models. These findings suggest that our model is more clinically relevant and serves as a valuable reference for determining the need for ALND.
The IBCSG23-01, AATRM, and Z0011 studies revealed no statistically significant disparities in DFS and OS between SLNB-alone and ALND groups [10, 43, 44]. Moreover, the SLNB-alone group exhibited a decrease in the incidence of postoperative upper extremity lymphedema and sensory and motor nerve injury. The AMAROS and OTOASOR studies similarly found no significant difference in prognosis between axillary radiotherapy and ALND groups, and the axillary recurrence rate was extremely low in both groups [11, 12]. Together, these findings imply that the option of disregarding ALND or selecting axillary RNI may be viable in cases of early breast cancer with a limited or moderate tumor burden in the SLNs. Current clinical trials have mostly focused on patients with 1 to 2 positive SLNs. Previous studies and our results also confirm that patients with 3 or more positive SLNs do not necessarily have non-SLN metastases, which would perhaps exempt them from ALND as well [1, 36]. Our new model predicted this subgroup and found an AUC of 0.843, which demonstrates the ability of the new model to accurately evaluate the non-SLN status of individuals diagnosed with breast cancer with 3 or more positive SLNs. Therefore, our model has the ability to precisely forecast the non-SLN status and, consequently, assess the tumor load of ALNs in Chinese women diagnosed with early breast cancer, which can provide valuable information for the subsequent decision-making regarding ALND or additional axillary RNI.
However, our study also has some limitations. Firstly, the study included only Chinese patients from a single hospital, thus lacking the diversity and representation of a larger multicenter cohort. Secondly, the metastatic SLN location has not appeared in previous models, and further validation with a larger sample size is required to establish its reliability. Lastly, patients with ITCs in SLN who do not undergo ALND were excluded from our cohort, thereby limiting the generalizability of our model to this subgroup of patients.