Skip to main content

Development and validation of a point-based scoring system for predicting axillary lymph node metastasis and disease outcome in breast cancer using clinicopathological and multiparametric MRI features

Abstract

Background

Axillary lymph node (ALN) metastasis is used to select treatment strategies and define the prognosis in breast cancer (BC) patients and is typically assessed using an invasive procedure. Noninvasive, simple, and reliable tools to accurately predict ALN status are desirable. We aimed to develop and validate a point-based scoring system (PSS) for stratifying the ALN metastasis risk of BC based on clinicopathological and quantitative MRI features and to explore its prognostic significance.

Methods

A total of 219 BC patients were evaluated. The clinicopathological and quantitative MRI features of the tumors were collected. A multivariate logistic regression analysis was used to create the PSS. The performance of the models was evaluated using receiver operating characteristic curves, and the area under the curve (AUC) of the models was calculated. Kaplan–Meier curves were used to analyze the survival outcomes.

Results

Clinical features, including the American Joint Committee on Cancer (AJCC) stage, T stage, human epidermal growth factor receptor-2, estrogen receptor, and quantitative MRI features, including maximum tumor diameter, Kep, Ve, and TTP, were identified as risk factors for ALN metastasis and were assigned scores for the PSS. The PSS achieved an AUC of 0.799 in the primary cohort and 0.713 in the validation cohort. The recurrence-free survival (RFS) and overall survival (OS) of the high-risk (> 19.5 points) groups were significantly shorter than those of the low-risk (≤ 19.5 points) groups in the PSS.

Conclusion

PSS could predict the ALN metastasis risk of BC. A PSS greater than 19.5 was demonstrated to be a predictor of short RFS and OS.

Background

Breast cancer (BC) is the most commonly diagnosed cancer among women worldwide and has become the second leading cause of cancer-related death [1]. BC can spread to the regional lymph nodes, primarily to the axillary and internal mammary nodes and subsequently to the medial supraclavicular nodes [2]. Patients with BC frequently experience axillary lymph node (ALN) metastasis, which affects the clinical stage, therapy options, surgical approach, and patient prognosis [1, 2]. Thus, accurate identification of ALN involvement in BC patients is crucial for prognosis and treatment decisions, particularly in the current time of downgrading off axillary surgery. Currently, ALN dissection, sentinel lymph node biopsy (SLNB), and lymph node biopsy before surgery are important procedures for ALN staging. Generally, lymph node biopsy is very important in daily clinical practice for pretreatment staging in patients with suspicious lymph node metastasis. SLNB is the current best standard approach for axillary staging in patients with clinically node-negative BC [3, 4]. ALN dissection is the standard treatment for BC patients with clinically positive nodes, but it might be avoided in patients with negative SLNB, as well as in patients with one or two sentinel lymph node-positivity receiving breast radiation and systemic therapy [5, 6]. However, these are invasive methods with the risk of postoperative complications, especially for ALN dissection [7, 8].

As noninvasive approaches, physical examination, mammography, ultrasonography, and positron emission tomography–computed tomography (PET/CT) are widely used to predict ALN metastasis. However, their abilities to assess the risk of ALN metastasis are limited because of their high false-negative rates and low sensitivity [9,10,11,12]. To date, only a few studies have used preoperative MRI to predict ALN metastasis, and these studies have shown that MRI morphological features and enhancement parameters of tumors are correlated with ALN status [13,14,15]. However, these models have not been widely used due to their complexity in clinical application. Furthermore, the mere identification of ALN metastasis by imaging is insufficient to change the paradigm for axillary surgical strategy [3]. A recently developed and validated nomogram was used to predict the metastasis status of ALNs based on a combination of clinicopathologic and MRI features, in which the developed model performed well for identifying ALN metastasis [10]. However, using the nomogram might be time-consuming and hard to interpret [16]. Therefore, to simplify the eventual clinical application of the predictive model, a point-based scoring system (PSS) was proposed and has been used to predict lymphovascular invasion and differential diagnosis of tumors [16, 17]. In this study, we aimed to use the PSS to predict the risk of ALN metastasis based on clinicopathological and quantitative MRI features and investigated the prognostic effect of the PSS.

Materials and methods

Patient characteristics

This retrospective study was approved by the ethics committee of Meizhou People's Hospital and the requirement for written informed consent was waived. We reviewed a total of 457 consecutive female patients with a postoperative histopathological diagnosis of BC who had undergone preoperative dynamic contrast-enhanced MRI (DCE-MRI) and diffusion-weighted imaging (DWI) scans and axillary biopsy or SLNB (complete ALN dissection was performed if the SLNB or axillary biopsy was positive) from January 2016 to May 2019. The patient exclusion criteria were as follows: (1) patients who had received neoadjuvant therapy prior to breast MRI or surgery; (2) patients with breast cancer recurrence; (3) patients with occult breast cancer; (4) patients with a tumor smaller than 1.0 cm;(5) patients with incomplete clinicopathological data; and (6) patients with poor visualization of the tumors on breast MRI. The patient inclusion/exclusion criteria are shown in Fig. 1. Finally, a total of 219 patients were enrolled and randomly divided into a primary cohort and a validation cohort at a ratio of 7:3 according to previous studies [18, 19]. Histological results after surgery were used as the ground truth of ALN metastasis. This study's patients are part of a large retrospective breast MRI database, of which 165 patients have been reported in a previously published study [20]. A previous report evaluated BC receptor status and molecular subtypes; however, this novel study focused on ALN metastasis by developing a new PSS.

Fig. 1
figure 1

Patient selection flowchart

Clinicopathological characteristics

Clinicopathological data, including age, Ki-67 level, human epidermal growth factor receptor-2 (HER-2), progesterone receptor (PR), estrogen receptor (ER), tumor node metastasis (TNM) stage, American Joint Committee on Cancer (AJCC) stage [21], pathological status of ALNs, number of metastatic ALNs, LVI, and perineural invasion status, were obtained from the patients' medical records. A pathological positive status of ALN metastasis was defined as macrometastasis (foci > 2.0 mm) or micrometastasis (foci 0.2–2.0 mm) identified by hematoxylin–eosin staining [7]. LVI was defined as the presence of tumor cells within the space of the endothelial cells [22]. The cutoff value for the positivity of HER-2, PR, and ER was used from internationally recognized standards [23, 24].

MRI technique

All images were acquired by a 3 T MR system (Magnetom Skyra, Siemens Healthcare) using a 16-channel bilateral breast coil. Table E1 in the supplementary material summarizes the MRI acquisition parameters. The contrast agent gadopentetate dimeglumine (Bayer Pharmaceuticals Corporation) at a dose of 0.2 ml/kg was given by intravenous injection at a rate of 3.0 ml/s.

Image analysis

MR image analyses were performed on the DCE software package (Tissue 4D, version: syngo MR D13, Siemens Healthcare), and DCE-derived parametric maps were automatically generated after motion correction. Breast tumors were identified on the DCE-MRI images, as the prominent enhancement area corresponded to the high signal area in the DWI image (b = 800 s/mm2) with a low signal area on the ADC map. Two radiologists (W.F. and F.C., with more than 15 years of experience) who were both blinded to the clinicopathological data reviewed the DCE-MRI images and ADC maps independently to extract the quantitative parameters of the tumors, including the volume transfer constant (Ktrans, min−1), reverse reflux rate constant (Kep, min−1), extracellular extravascular volume fraction (Ve), rate of contrast enhancement for inflow (W-in, min−1), rate of contrast decay for outflow (W-out, min−1), time to peak enhancement after injection (TTP, min), and ADC value (× 10 −3 mm2/s). The region of interest (ROI) drawing principles were as follows (Figure E1 in the supplementary material): ROIs with a minimum area of 0.10 cm2 were manually drawn on the continuous three maximum sections with the greatest enhancement areas of the tumors, avoiding visible blood vessels, obvious bleeding, and necrotic and cystic areas, with the same position and size on both the DCE-derived parametric maps and the ADC maps. The mean values of the parameters were taken as the final values. In addition, tumor sizes, including the maximum, minimum, and effective diameter of the tumor on the largest section, were recorded on the DCE-MRI images.

Clinical treatment

All patients underwent surgical treatment after breast MRI. The surgical methods for BC include breast-conserving surgery, mastectomy, and breast reconstruction surgery, which were determined by the National Guidelines for the Diagnosis and Treatment of Breast Cancer 2022 in China (English version)[23]. Of 219 patients, 21 patients had received breast-conserving surgery, 189 patients had received mastectomy, and 9 patients had received breast reconstruction surgery. Furthermore, 171 patients received adjuvant chemotherapy after surgery, 71 received radiation therapy after surgery, and 121 received adjuvant endocrine therapy after surgery. Moreover, 43 patients had received all three treatment methods after surgery.

Follow-up

All patients underwent clinical and imaging follow-up with chest CT, breast mammography, or breast MRI until progression, followed by routine follow-up until death. According to the Chinese Anti-Cancer Association guidelines for the diagnosis and treatment of breast cancer (2021 version), all patients were postoperatively followed up with a physical examination, laboratory tests, CT, mammography, and ultrasound every 3 months for the first two years, once every 6 months during years 3 to 5, and once a year after 5 years. Additionally, the patients with breast-conserving surgery were followed up with breast MRI once a year. Recurrence-free survival (RFS) was measured in months from the date of surgery to the first date of local recurrence, distant metastasis, or the last follow-up date, whichever came first. Overall survival (OS) was measured in months from the date of surgery to the date of death or the last follow-up date. The last follow-up date was March 1, 2021.

Point-based scoring system development

A PSS was constructed using the method described by Sullivan et al. [25]. Continuous variables were transformed into binary variables according to the cutoff values determined by the largest Youden index to reduce the PSS score range and simplify its eventual clinical application. In the primary cohort, the cutoff values for age, maximum tumor diameter, minimum tumor diameter, effective tumor diameter, Ktrans, Kep, Ve, W-in, W-out, TTP, and ADC were 48.5 years,3.55 cm,1.95 cm,3.13 cm,0.298 min−1,0.982 min−1,0.253, 0.674 min−1,-0.013 min−1,0.778 min, and 0.736 × 10−3 mm2/s, respectively. All variables were included in the multivariate logistic regression analysis following the backward stepwise selection method. The point of each selected variable was assigned based on its β-coefficient in the multivariate logistic regression analysis and rounded to the nearest 0.5. The final PSS score was obtained by summing the points of each variable. The odds ratios (ORs) with 95% confidence intervals (CIs) of each variable were also calculated.

Statistical analysis

R version 3.6.0 was used for statistical analysis. Continuous variables with a normal distribution are expressed as the mean ± standard deviation, and categorical variables are expressed as counts and percentages. Data between the ALN metastasis and non-ALN metastasis groups were compared with Student's t test, chi-squared test, or Kruskal–Wallis H test, if appropriate. In addition, the intraclass correlation coefficient (ICC) was used to assess the interobserver agreement of the DCE-MRI parameters and the ADC values between the two radiologists [26, 27]. An ICC value between 0.61 and 0.80 indicates good agreement, while an ICC value ≥ 0.81 indicates excellent agreement [28]. A receiver operating characteristic (ROC) curve was created to estimate the discriminating power of the models, and the area under the curve (AUC), accuracy, sensitivity and specificity were calculated [27]. The ROCs of the models were compared by using the DelLong test. For the survival analysis, the patients were divided into low-risk and high-risk groups according to the PSS point, and the grouping threshold was determined by the X-tile method [29]. Kaplan–Meier curves for the high-risk and low-risk groups were plotted, and the prognosis of the high- and low-risk groups was compared by using the log rank test. P < 0.05 was considered statistically significant.

Results

Patient characteristics

Of 219 patients, 110 patients had non-ALN metastasis and 109 patients had ALN metastasis. The characteristics of the patients in the total, primary, and validation cohorts are listed in Table 1. The ALN metastasis positivity was 49.7% (76/153) and 50.0% (33/66) in the primary cohort and validation cohort, respectively. The patients with ALN metastasis had higher maximum, minimum, and effective tumor diameters, and higher T stages than the patients without ALN metastasis in the total, primary, and validation cohorts (all P < 0.05). Compared to the patients without ALN metastasis, higher positive rates of LVI and lower Ve and W-out values were more pronounced in the patients with ALN metastasis in the total and primary cohorts (all P < 0.05), but these differences were not confirmed in the validation cohort (all P > 0.05). In addition, the patients with ALN metastasis had a higher Ki-67 index than patients without ALN metastasis in the total cohort (P = 0.038), but this was not confirmed in the primary and validation cohorts (all P > 0.05).

Table 1 Characteristics of patients with and without axillary lymph node metastasis

Agreement between two readers

The ICC values between the two readers for Ktrans, Kep, Ve, W-in, W-out, TTP, and ADC were 0.995 (95% CI: 0.994–0.996), 0.775 (95% CI: 0.716–0.823), 0.988 (95% CI: 0.984–0.990), 0.990 (95% CI: 0.988–0.993), 0.793 (95% CI: 0.739–0.838), 0.965 (95% CI: 0.955–0.973), and 0.894 (95% CI: 0.863–0.917), respectively.

Point-based scoring system development based on the primary cohort

Based on the backward stepwise selection method in the multivariate logistic regression analysis (Table 2), eight variables were associated with ALN metastasis and they were assigned scores for the final prediction rule: AJCC stage [II (1.0 points), III (17.5 points), IV (16.0 points)], T stage [T2 (1.0 points), T3 (17.0 points), T4 (2.0 points)], ER (positive,16.0 points), HER-2 (positive,1.5 points), maximum tumor diameter(≥ 3.55 cm,1.0 points), Kep (≥ 0.982 min−1,1.0 points), Ve (≥ 0.253, -0.5 points), and TTP (≥ 0.778 min−1,1.0 points). The PSS point threshold for predicting the risk of ALN metastasis was 17.5, as obtained by the largest Youden index.

Table 2 Multivariate logistic regression analysis for predicting axillary lymph node metastasis with the assigned points

Performance of the point-based scoring system

The logistic regression model that used all of the predictors for ALN metastasis obtained an AUC of 0.802 in the primary cohort and 0.707 in the validation cohort, and the PSS achieved an AUC of 0.799 in the primary cohort and 0.713 in the validation cohort (Table 3). ROC curve comparisons between the logistic regression model and the PSS are shown in Fig. 2. There were no significant AUC differences between the logistic regression model and PSS in either the primary cohort (P = 0.645) or in the validation cohort (P = 0.572). Table 4 presents the risk of ALN metastasis according to the PSS, and the ALN metastasis risk ranged from 4.7% to 100.0%. Figure 3 and Fig. 4 represent examples of the PSS in use.

Table 3 Performance of the prediction model for axillary lymph node metastasis
Fig. 2
figure 2

Receiver operator characteristic (ROC) curve for the prediction models in the primary cohort (A) and validation cohort (B)

Table 4 Axillary lymph node metastasis risk according to the point-based scoring system
Fig. 3
figure 3

One example of the point-based scoring system for predicting breast cancer patients without axillary lymph node metastasis. A ~ F are MRI images of a female right breast cancer patient with invasive ductal carcinoma who was positive for human epidermal growth factor receptor-2, negative for estrogen receptor, AJCC stage III, and T2 stage. An axial contrast-enhanced T1WI image (A) shows a tumor with a maximum diameter of 2.3 cm located in the lower-outer quadrant. The mean values of Kep, Ve, and TTP from the corresponding pseudocolor images of Kep (B), Ve, (C), and TTP (D) were 0.760 min−1, 0.190, and 0.620 min, respectively. Breast MRI was suspicious for axillary lymph node metastasis based on the following features: axial T2WI with fat suppression imaging (E) and axial contrast-enhanced T1WI imaging (F) showing cortical thickening, an oval shape and a long-to-short axis ratio of 1.4 for a maximal lymph node (white arrow) in axillary level I (nodes lateral and inferior to pectoralis minor muscle). However, the risk of axillary lymph node metastasis assessed by the scoring system was 17.0, with an axillary lymph node metastasis probability of 26.9%. The final postoperative pathology report showed that this tumor had no axillary lymph node metastasis

Fig. 4
figure 4

One example of the point-based scoring system for predicting breast cancer patients with axillary lymph node metastasis. A ~ F are MRI images of a female left breast cancer patient with invasive ductal carcinoma who was negative for human epidermal growth factor receptor-2, positive for estrogen receptor, AJCC stage II and T2 stage. An axial contrast-enhanced T1WI image (A) shows a tumor with a maximum diameter of 3.7 cm located in the upper-outer quadrant. The mean values of Kep, Ve, and TTP from the corresponding pseudocolor images of Kep (B), Ve, (C), and TTP (D) were 1.130 min−1, 0.310, and 0.530 min, respectively. Breast MRI diagnosed no axillary lymph node metastasis for level I nodes due to the lack of all of the following features from axial T2WI with fat suppression imaging (E) and contrast-enhanced T1WI imaging (F): cortical thickening, missing fatty hilum, round shape, or a long-to-short axis ratio of less than 2 (long-to-short axis ratio was 2.2 for this patient). However, the risk of axillary lymph node metastasis assessed by the scoring system was 19.5,with an axillary lymph node metastasis probability of 81.8%.The final postoperative pathology report showed that this tumor had three axillary lymph node metastases at axillary level I (white arrow)

Predictors of survival

All patients had completed RFS follow-up, and the median RFS of all patients was 32 [interquartile range (IQR):23.0–44.0] months, of which 26 patients had tumor recurrence and 193 patients had no tumor recurrence. The recurrence rate was 19.3% in the ALN metastasis group and 4.5% in the non-ALN metastasis group. In addition, 9 patients with tumor recurrence were deceased and 2 patients without tumor recurrence were deceased. While 218 patients had completed OS follow-up, the median OS of all patients was 43.8 (IQR:32.0–55.0) months, of which 11 patients were deceased. Of note, the overall death rate was 8.3% in the ALN metastasis group and 1.8% in the non-ALN metastasis group.

The PSS point threshold for the low-risk group and high-risk group was 19.5, as obtained by the X-tile method [29]. The median RFS was 32.0 (IQR:15.5–46.0) months for the high-risk group (> 19.5 points) and 32.5 (IQR:22.9–44.3) for the low-risk group (≤ 19.5 points) in the primary cohort (the comparative P value for the two survival curves using the log-rank test was lower than 0.0001, hereinafter, Fig. 5A) and 32.0 (IQR:21.0–40.0) months for the high-risk group and 32.0 (IQR:23.5–40.0) for the low-risk group in the validation cohort (P = 0.032, Fig. 5B). The Kaplan–Meier survival analysis showed that the median OS was 43.3 (IQR:31.9–56.0) months for the high-risk group and 50.0 (IQR:29.0–55.5) for the low-risk group in the primary cohort (P = 0.042, Fig. 6A) and 43.5 (IQR:32.0–55.0) months for the high-risk group and 50.5 (IQR:45.0–55.0) for the low-risk group in the validation cohort (P = 0.029, Fig. 6B).

Fig. 5
figure 5

Kaplan–Meier recurrence-free survival (RFS) analyses of the low-risk and high-risk groups in the PSS. The high-risk group patients had a significantly worse RFS than the low-risk groups patients in both the primary cohort (A) and validation cohort (B)

Fig. 6
figure 6

Kaplan–Meier overall survival (OS) analyses of the low-risk and high-risk groups in the PSS. The high-risk group patients had a significantly worse OS than the low-risk group patients in both the primary cohort (A) and validation cohort (B)

Discussion

In the present study, we developed and validated a PSS for stratifying the ALN metastasis risk of BC based on the combination of clinical features and quantitative MRI features and investigated the prognostic effect of the PSS. Our results showed that the PSS has good discrimination for stratifying the ALN metastasis risk of BC in the primary cohort and in the validation cohort. In addition, the RFS and OS of the high-risk group (> 19.5 points) patients were significantly lower than those of the low-risk (≤ 19.5 points) group patients in the PSS.

Risk stratification scoring systems have been used to guide clinical decision-making [16, 30,31,32,33]. A recent study conducted by MY et al. developed a scoring system to stratify the risk of lymphovascular invasion in BC patients with an AUC of 0.824 [16]. Another study from Ouldame et al. used a scoring system to predict axillary response after neoadjuvant chemotherapy in initially node-positive women with BC [34]. However, little is known about the ALN based on this scoring system. In this study, to develop the risk stratification scoring system, a multivariable logistic regression analysis was used to select the risk factors for ALN metastasis. The PSS point threshold for predicting the risk of ALN metastasis is 17.5 and is very easy for clinicians to use. In terms of assessing ALN metastasis, the scoring system showed good performance in the primary (AUC = 0.799) and validation cohorts (AUC = 0.713), which was similar to the result from Yang et al. [35]. Contrary to previously developed models that might be time-consuming and hard to apply, this risk stratification scoring system provides an easy tool for surgeons to assess breast cancer patients' risk for ALN metastasis prior to surgery [7, 10, 36, 37].

In this study, eight candidates, including AJCC stage, T stage, ER status, HER-2 status, maximum tumor diameter, Kep, Ve, and TTP, were identified as risk factors for ALN metastasis from a stepwise logistic regression model. Our results showed that the selected variables in the scoring system were in line with those of the results reported in previous studies. For instance, Chen et al. found that HER-2 and tumor size were related to breast cancer ALN metastasis [38]. Yang et al. reported that ER and T stage were significantly associated with ALN status [35]. Ya et al. showed the ability of dynamic contrast-enhanced MRI parameters of Kep and TTP to predict ALN metastasis in breast cancer patients [39]. Among the risk factors in this study, clinicopathologic risk factors for AJCC stage, T stage, ER, and HER-2 were more important than the MRI quantitative parameters of maximum tumor diameter, Kep, Ve, and TTP based on their β-coefficient results in the logistic regression model. These results were consistent with previous results, in which clinicopathological risk factors were better predictors than imaging biomarkers [18, 40]. In addition, Kep plays a more important role in the PSS than the parameters of Ve and TTP. The possible explanations may be as follows: Kep may be less sensitive to the absolute value of the contrast agent concentration compared with Ve and TT. Another reason could be due to the complexity of the complex pathological characteristics of BC, which were consistent with previous results [28, 41].

The determination of ALN metastasis has important clinical significance for patients' overall recurrence and survival [3, 42]. Generally, patients with ALN metastasis have poor outcomes and a higher risk of local recurrence and distant metastasis than patients without ALN metastasis. Similar results were observed in our study: the patients with ALN metastasis had a higher recurrence rate and overall death rate than the patients without ALN metastasis (recurrence rate:19.3% vs. 4.5%, and death rate 5.0% vs. 1.8%, respectively). To explore prognostic significance of the PSS, we compared the prognosis of high- and low-risk groups in the PSS. The Kaplan–Meier survival analysis showed that the median RFS and OS of the high-risk group patients were significantly lower than those of the low-risk group patients in both the primary cohort and validation cohort, suggesting that the PSS can be routinely used in predicting the prognosis of BC patients.

This study has several limitations. First, the PSS was developed based on the combination of clinicopathological and MRI features of BC instead of the ALN itself because it is difficult to match the biopsied or dissected ALN on MRI imaging. Second, these retrospective study results were obtained in a single institution with some inherent limitations and the results should be validated in a multicenter study. Third, some patients were followed up for less than 5 years, which partially affected the strength of prognostic information.

Conclusion

In conclusion, a simple and reliable point-based scoring system can be used to stratify the ALN metastasis risk of BC. PSS greater than 19.5 was confirmed to be a predictor of short RFS and OS.

Availability of data and materials

The data cohorts used and/or analyzed in the present study are available from the corresponding authors upon reasonable request.

Abbreviations

ALN:

Axillary lymph node

BC:

Breast cancer

ER:

Estrogen receptor

PR:

Progesterone receptor

HER-2:

Human epidermal growth factor receptor-2

PSS:

Point-based scoring system

RFS:

Recurrence-free survival

OS:

Overall survival

References

  1. Poortmans PM, Weltens C, Fortpied C, Kirkove C, Peignaux-Casasnovas K, Budach V, van der Leij F, Vonk E, Weidner N, Rivera S, et al. Internal mammary and medial supraclavicular lymph node chain irradiation in stage I-III breast cancer (EORTC 22922/10925): 15-year results of a randomised, phase 3 trial. Lancet Oncol. 2020;21(12):1602–10.

    Article  PubMed  Google Scholar 

  2. Zheng X, Yao Z, Huang Y, Yu Y, Wang Y, Liu Y, Mao R, Li F, Xiao Y, Wang Y, et al. Deep learning radiomics can predict axillary lymph node status in early-stage breast cancer. Nat Commun. 2020;11(1):1236.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Chang JM, Leung JWT, Moy L, Ha SM, Moon WK. Axillary nodal evaluation in breast cancer: state of the art. Radiology. 2020;295(3):500–15.

    Article  PubMed  Google Scholar 

  4. Esposito E, Di Micco R, Gentilini OD. Sentinel node biopsy in early breast cancer. A review on recent and ongoing randomized trials. Breast. 2017;36:14–9.

    Article  PubMed  Google Scholar 

  5. Magnoni F, Galimberti V, Corso G, Intra M, Sacchini V, Veronesi P. Axillary surgery in breast cancer: An updated historical perspective. Semin Oncol. 2020;47(6):341–52.

    Article  PubMed  Google Scholar 

  6. Noguchi M, Inokuchi M, Noguchi M, Morioka E, Ohno Y, Kurita T. Axillary surgery for breast cancer: past, present, and future. Breast Cancer (Tokyo, Japan). 2021;28(1):9–15.

    Article  PubMed  Google Scholar 

  7. Zhao M, Wu Q, Guo L, Zhou L, Fu K. Magnetic resonance imaging features for predicting axillary lymph node metastasis in patients with breast cancer. Eur J Radiol. 2020;129:109093.

    Article  PubMed  Google Scholar 

  8. Zhang J, Li X, Huang R, Feng WL, Kong YN, Xu F, Zhao L, Song QK, Li J, Zhang BN, et al. A nomogram to predict the probability of axillary lymph node metastasis in female patients with breast cancer in China: A nationwide, multicenter, 10-year epidemiological study. Oncotarget. 2017;8(21):35311–25.

    Article  PubMed  Google Scholar 

  9. Zhang YN, Wang CJ, Xu Y, Zhu QL, Zhou YD, Zhang J, Mao F, Jiang YX, Sun Q. Sensitivity, specificity and accuracy of ultrasound in diagnosis of breast cancer metastasis to the axillary lymph nodes in chinese patients. Ultrasound Med Biol. 2015;41(7):1835–41.

    Article  PubMed  Google Scholar 

  10. Han L, Zhu Y, Liu Z, Yu T, He C, Jiang W, Kan Y, Dong D, Tian J, Luo Y. Radiomic nomogram for prediction of axillary lymph node metastasis in breast cancer. Eur Radiol. 2019;29(7):3820–9.

    Article  PubMed  Google Scholar 

  11. An YS, Lee DH, Yoon JK, Lee SJ, Kim TH, Kang DK, Kim KS, Jung YS, Yim H. Diagnostic performance of 18F-FDG PET/CT, ultrasonography and MRI. Detection of axillary lymph node metastasis in breast cancer patients. Nuklearmedizin. 2014;53(3):89–94.

    Article  CAS  PubMed  Google Scholar 

  12. Marino MA, Avendano D, Zapata P, Riedl CC, Pinker K. Lymph node imaging in patients with primary breast cancer: concurrent diagnostic tools. Oncologist. 2020;25(2):e231–42.

    Article  PubMed  Google Scholar 

  13. Loiselle C, Eby PR, Kim JN, Calhoun KE, Allison KH, Gadi VK, Peacock S, Storer BE, Mankoff DA, Partridge SC, et al. Preoperative MRI improves prediction of extensive occult axillary lymph node metastases in breast cancer patients with a positive sentinel lymph node biopsy. Acad Radiol. 2014;21(1):92–8.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Loiselle CR, Eby PR, DeMartini WB, Peacock S, Bittner N, Lehman CD, Kim JN. Dynamic contrast-enhanced MRI kinetics of invasive breast cancer: a potential prognostic marker for radiation therapy. Int J Radiat Oncol Biol Phys. 2010;76(5):1314–9.

    Article  PubMed  Google Scholar 

  15. Tuncbilek N, Karakas HM, Okten OO. Dynamic magnetic resonance imaging in determining histopathological prognostic factors of invasive breast cancers. Eur J Radiol. 2005;53(2):199–205.

    Article  PubMed  Google Scholar 

  16. Ni-Jia-Ti MY, Ai-Hai-Ti DL, Huo-Jia AS, Wu-Mai-Er PL, A-bu-li-zi AB, Shi Y, Rou-zi NE, Su WJ, Dai GZ, Da-mo-la MH. Development of a risk-stratification scoring system for predicting lymphovascular invasion in breast cancer. BMC Cancer. 2020;20(1):94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Lee MW, Kim GH, Kim KB, Kim YH, Park DY, Choi CI, Kim DH, Jeon TY. Digital image analysis-based scoring system for endoscopic ultrasonography is useful in predicting gastrointestinal stromal tumors. Gastric Cancer. 2019;22(5):980–7.

    Article  PubMed  Google Scholar 

  18. Chen X, Chen X, Yang J, Li Y, Fan W, Yang Z. Combining Dynamic Contrast-Enhanced Magnetic Resonance Imaging and Apparent Diffusion Coefficient Maps for a Radiomics Nomogram to Predict Pathological Complete Response to Neoadjuvant Chemotherapy in Breast Cancer Patients. J Comput Assist Tomogr. 2020;44(2):275–83.

    Article  PubMed  Google Scholar 

  19. Chen X, Yang Z, Yang J, Liao Y, Pang P, Fan W, Chen X. Radiomics analysis of contrast-enhanced CT predicts lymphovascular invasion and disease outcome in gastric cancer: a preliminary study. Cancer Imag. 2020;20(1):24.

    Article  Google Scholar 

  20. Yang Z, Chen X, Zhang T, Cheng F, Liao Y, Chen X, Dai Z, Fan W. Quantitative Multiparametric MRI as an imaging biomarker for the prediction of breast cancer receptor status and molecular subtypes. Front Oncol. 2021;11(3692):628824.

  21. Kalli S, Semine A, Cohen S, Naber SP, Makim SS, Bahl M. American Joint Committee on cancer’s staging system for breast cancer, eighth edition: what the radiologist needs to know. Radiographics. 2018;38(7):1921–33.

    Article  PubMed  Google Scholar 

  22. Hamy AS, Lam GT, Laas E, Darrigues L, Balezeau T, Guerin J, Livartowski A, Sadacca B, Pierga JY, Vincent-Salomon A, et al. Lymphovascular invasion after neoadjuvant chemotherapy is strongly associated with poor prognosis in breast carcinoma. Breast Cancer Res Treat. 2018;169(2):295–304.

    Article  PubMed  Google Scholar 

  23. Health Commission Of The People’s Republic Of China N. National guidelines for diagnosis and treatment of breast cancer 2022 in China (English version). Chin J Cancer Res=Chung kuo yen cheng yen chiu. 2022;34(3):151–75.

    Article  PubMed  Google Scholar 

  24. Liang X, Chen X, Yang Z, Liao Y, Wang M, Li Y, Fan W, Dai Z, Zhang Y. Early prediction of pathological complete response to neoadjuvant chemotherapy combining DCE-MRI and apparent diffusion coefficient values in breast Cancer. BMC Cancer. 2022;22(1):1250.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Sullivan LM, Massaro JM, D’Agostino RB Sr. Presentation of multivariate data for clinical use: The Framingham Study risk score functions. Stat Med. 2004;23(10):1631–60.

    Article  PubMed  Google Scholar 

  26. Yang Z, Chen X, Huang R, Li S, Lin D, Yang Z, Sun H, Liu G, Qiu J, Tang Y, et al. Atypical presentations of coronavirus disease 2019 (COVID-19) from onset to readmission. BMC Infect Dis. 2021;21(1):127.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Zhong M, Yang Z, Chen X, Huang R, Wang M, Fan W, Dai Z, Chen X. Readout-segmented echo-planar diffusion-weighted MR imaging improves the differentiation of breast cancer receptor statuses compared with conventional diffusion-weighted imaging. J Magnet Reson Imaging. 2022;56(3):691–9.

    Article  Google Scholar 

  28. Lai T, Chen X, Yang Z, Huang R, Liao Y, Chen X, Dai Z. Quantitative parameters of dynamic contrast-enhanced magnetic resonance imaging to predict lymphovascular invasion and survival outcome in breast cancer. Cancer Imaging. 2022;22(1):61.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Wu J, Zhang H, Li L, Hu M, Chen L, Xu B, Song Q. A nomogram for predicting overall survival in patients with low-grade endometrial stromal sarcoma: A population-based analysis. Cancer Commun (London, England). 2020;40(7):301–12.

    Article  Google Scholar 

  30. Turner BM, Gimenez-Sanders MA, Soukiazian A, Breaux AC, Skinner K, Shayne M, Soukiazian N, Ling M, Hicks DG. Risk stratification of ER-positive breast cancer patients: A multi-institutional validation and outcome study of the Rochester Modified Magee algorithm (RoMMa) and prediction of an Oncotype DX(®) recurrence score <26. Cancer Med. 2019;8(9):4176–88.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Lai J, Wang H, Peng J, Chen P, Pan Z. Establishment and external validation of a prognostic model for predicting disease-free survival and risk stratification in breast cancer patients treated with neoadjuvant chemotherapy. Cancer Manage Res. 2018;10:2347–56.

    Article  CAS  Google Scholar 

  32. Li JL, Lin XY, Zhuang LJ, He JY, Peng QQ, Dong YP, Wu JX. Establishment of a risk scoring system for predicting locoregional recurrence in T1 to T2 node-negative breast cancer patients treated with mastectomy: Implications for postoperative radiotherapy. Medicine (Baltimore). 2017;96(26):e7343.

    Article  PubMed  Google Scholar 

  33. Giuliani M, Rinaldi P, Rella R, D’Angelo A, Carlino G, Infante A, Romani M, Bufi E, Belli P, Manfredi R. A new risk stratification score for the management of ultrasound-detected B3 breast lesions. Breast J. 2018;24(6):965–70.

    Article  CAS  PubMed  Google Scholar 

  34. Ouldamer L, Chas M, Arbion F, Body G, Cirier J, Ballester M, Bendifallah S, Daraï E. Risk scoring system for predicting axillary response after neoadjuvant chemotherapy in initially node-positive women with breast cancer. Surg Oncol. 2018;27(2):158–65.

    Article  PubMed  Google Scholar 

  35. Yang J, Wang T, Yang L, Wang Y, Li H, Zhou X, Zhao W, Ren J, Li X, Tian J, et al. Preoperative Prediction of Axillary Lymph Node Metastasis in Breast Cancer Using Mammography-Based Radiomics Method. Sci Rep. 2019;9(1):4429.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Yu Y, Tan Y, Xie C, Hu Q, Ouyang J, Chen Y, Gu Y, Li A, Lu N, He Z, et al. Development and validation of a preoperative magnetic resonance imaging radiomics-based signature to predict axillary lymph node metastasis and disease-free survival in patients with early-stage breast cancer. JAMA Netw Open. 2020;3(12):e2028086.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Arefan D, Chai R, Sun M, Zuley ML, Wu S. Machine learning prediction of axillary lymph node metastasis in breast cancer: 2D versus 3D radiomic features. Med Phys. 2020;47(12):6334–42.

    Article  PubMed  Google Scholar 

  38. Chen W, Wang C, Fu F, Yang B, Chen C, Sun Y. A model to predict the risk of lymph node metastasis in breast cancer based on clinicopathological characteristics. Cancer Manag Res. 2020;12:10439–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Ya G, Wen F, Xing-Ru L, Zhuan-Zhuan G, Jun-Qiang L. Difference of DCE-MRI parameters at different time points and their predictive value for axillary lymph node metastasis of breast cancer. Acad Radiol. 2022;29(Suppl 1):S79–86.

    Article  PubMed  Google Scholar 

  40. Liu Z, Li Z, Qu J, Zhang R, Zhou X, Li L, et al. Radiomics of multiparametric MRI for pretreatment prediction of pathologic complete response to neoadjuvant chemotherapy in breast cancer: a multicenter study. Clin Cancer Res. 2019;25(12):3538–47.

    Article  CAS  PubMed  Google Scholar 

  41. Fedorov A, Fluckiger J, Ayers GD, Li X, Gupta SN, Tempany C, Mulkern R, Yankeelov TE, Fennessy FM. A comparison of two methods for estimating DCE-MRI parameters via individual and cohort based AIFs in prostate cancer: a step towards practical implementation. Magn Reson Imaging. 2014;32(4):321–9.

    Article  PubMed  PubMed Central  Google Scholar 

  42. McGale P, Taylor C, Correa C, Cutter D, Duane F, Ewertz M, Gray R, Mannu G, Peto R, Whelan T, et al. Effect of radiotherapy after mastectomy and axillary surgery on 10-year recurrence and 20-year breast cancer mortality: meta-analysis of individual patient data for 8135 women in 22 randomised trials. Lancet. 2014;383(9935):2127–35.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

National Natural Science Foundation of China (82101985); Medical Scientific Foundation of Guangdong Province (B2021280).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: Zhuozhi Dai, Weixiong Fan. Data curation: Xiaofeng Chen, Zhiqi Yang. Formal analysis: Ruibin Huang, Yue Li. Funding acquisition: Xiaofeng Chen, Zhuozhi Dai. Methodology: Yuting Liao, Guijin Li, Mengzhu Wang. Project administration: Xiangguang Chen, Weixiong Fan. Supervision: Xiangguang Chen, Zhuozhi Dai. Writing—original draft: Xiaofeng Chen, Zhiqi Yang. Writing—review & editing: Zhuozhi Dai, Weixiong Fan. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Zhuozhi Dai or Weixiong Fan.

Ethics declarations

Ethics approval and consent to participate

This retrospective study was approved by the institutional review board of Meizhou People's Hospital.

Consent for publication

The institutional review board waived the need to obtain informed consent for this retrospective study.

Competing interests

One author of the study (YT L) is a consultant for GE Healthcare and was responsible for assisting with data, and another two authors of the study (MZ W) and (GJ L) are consultants for Siemens Healthcare and provided some critical technical support in the manuscript. Those authors have not controlled any data or information that may present a conflict of interest. The other authors are not employees or consultants in the industrial fields and had control of the data and information submitted for publication.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Table E1.

Shows the MRI acquisition parametersAQ. Figure E1. Measurement of MRI quantitative parameters of tumors in a right breast cancer patient with invasive ductal carcinoma.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, X., Yang, Z., Huang, R. et al. Development and validation of a point-based scoring system for predicting axillary lymph node metastasis and disease outcome in breast cancer using clinicopathological and multiparametric MRI features. Cancer Imaging 23, 54 (2023). https://doi.org/10.1186/s40644-023-00564-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40644-023-00564-9

Keywords