Skip to main content

Automated Breast Ultrasound (ABUS)-based radiomics nomogram: an individualized tool for predicting axillary lymph node tumor burden in patients with early breast cancer

Abstract

Objectives

Preoperative evaluation of axillary lymph node (ALN) status is an essential part of deciding the appropriate treatment. According to ACOSOG Z0011 trials, the new goal of the ALN status evaluation is tumor burden (low burden, < 3 positive ALNs; high burden, ≥ 3 positive ALNs), instead of metastasis or non-metastasis.

We aimed to develop a radiomics nomogram integrating clinicopathologic features, ABUS imaging features and radiomics features from ABUS for predicting ALN tumor burden in early breast cancer.

Methods

A total of 310 patients with breast cancer were enrolled. Radiomics score was generated from the ABUS images. Multivariate logistic regression analysis was used to develop the predicting model, we incorporated the radiomics score, ABUS imaging features and clinicopathologic features, and this was presented with a radiomics nomogram. Besides, we separately constructed an ABUS model to analyze the performance of ABUS imaging features in predicting ALN tumor burden. The performance of the models was assessed through discrimination, calibration curve, and decision curve.

Results

The radiomics score, which consisted of 13 selected features, showed moderate discriminative ability (AUC 0.794 and 0.789 in the training and test sets). The ABUS model, comprising diameter, hyperechoic halo, and retraction phenomenon, showed moderate predictive ability (AUC 0.772 and 0.736 in the training and test sets). The ABUS radiomics nomogram, integrating radiomics score with retraction phenomenon and US-reported ALN status, showed an accurate agreement between ALN tumor burden and pathological verification (AUC 0.876 and 0.851 in the training and test sets). The decision curves showed that ABUS radiomics nomogram was clinically useful and more excellent than US-reported ALN status by experienced radiologists.

Conclusions

The ABUS radiomics nomogram, with non-invasive, individualized and precise assessment, may assist clinicians to determine the optimal treatment strategy and avoid overtreatment.

Peer Review reports

Background

Breast cancer is the most common malignant tumor and the main cause of cancer-related death among women [1]. The axillary lymph node (ALN) status determines the need for systemic therapy, the extent of surgery, reconstruction options, and the need for radiation therapy after mastectomy. Accurate assessment of ALN status plays an important role in breast cancer treatment and prognosis.

Based on preoperative assessment, patients without suspicious ALN metastasis by ultrasound (US) should have a sentinel lymph node biopsy (SLNB). Patients with positive ALN metastasis on US are candidates for ultrasound-guided fine needle aspiration (FNA) or core needle biopsy (CNB). For a negative FNA/CNB, SLNB is indicated. If FNA/CNB is positive, axillary lymph node dissection (ALND) is indicated unless neoadjuvant therapy is given. ALND is associated with several complications, including lymphedema (found in up to 25% of women after surgery), infection, shoulder motion restriction, and major vascular and nerve damage [2]. Besides, clinical practice has demonstrated that an important number of patients undergo a secondary ALND when SLNB displays major lymph node involvement [3]. Therefore, accurate preoperative noninvasive assessment of ALN status is essential.

According to the findings of the Z0011 trials, the American Society of Clinical Oncology updated clinical practice guidelines, ALND can be omitted in patients with breast cancer with 1–2 positive SLNs without a decrease in disease-free survival or overall survival [4, 5]. So, the new goal of the ALN status evaluation is tumor burden (low burden, < 3 positive ALNs; high burden, ≥ 3 positive ALNs), instead of metastasis or non-metastasis.

Ultrasound is a widely-used tool in assessing the ALN status preoperatively as it is non-invasive, radiation-free, real-time, rapid, and convenient. Previous studies have proved that axillary ultrasound (AUS) may provide valuable information relevant to ALN status in breast cancer [6]. However, AUS mainly obtains visual image information and focuses on the qualitative analysis of lesions, the diagnostic performances of axillary ultrasound to detect ALN involvement highly depends on the experience of radiologists [3]. Ahmed et al. [7] showed that 43.2% of patients with positive ALN metastasis on US had a low lymph node tumor burden (ALN metastasis < 3). It means that almost half of patients with ALN metastasis positive US assessments are exempt from ALND, which some researchers believe may constitute overtreatment for these patients [8]. AUS may not be a reliable predictor for nodal metastasis [6, 9]. In the age of precision medicine, a more effective and individualized method is urgent to resolve this problem.

As an emerging three-dimensional imaging technique, Automated Breast Ultrasound (ABUS), addresses the limitations of conventional handheld ultrasound (HHUS) and automatically scans the breast based on special high-frequency broadband sensors [10]. Several recent studies [11, 12] have shown that some unique features of ABUS, although they are also visual assessment and qualitative descriptions, may provide additional information for breast lesions. Specifically, retraction phenomenon appears as a satellite model around the lesion, with high sensitivity (80%–89%) and specificity (96%–100%) for breast cancer [13,14,15]. Radiomics extracts high-throughput quantitative features that may not be directly observable with the naked eye from single or multiple medical images. Radiomics has been more recently applied to distinguish benign malignant breast lesions [16], predict lymph node status [17, 18], and even evaluate treatment response [19]. According to the radiomics quality score proposed by Lambin et al. [20], ABUS images with standardized, repeatable, and high-resolution characteristics would be fit for radiomics analysis. However, to our knowledge, there is no ABUS-based radiomics study to differentiate ALN tumor burden. For patients with high tumor burden can assist in identifying what kind of initial axillary surgery can overlook the SLNB and undergo ALND specifically and assist in employing neoadjuvant chemotherapy or individualized adjuvant radiotherapy. For patients with low ALN tumor burden, unnecessary treatment, and potential complications due to surgery can be avoided.

The purpose of this study was to develop a radiomics nomogram integrating clinicopathologic features, ABUS imaging features and radiomics features from ABUS for predicting ALN tumor burden in early breast cancer.

Methods

Patient selection

Ethical approval for this retrospective study was obtained from our institutional review board, and informed consent was canceled. All the patients with breast cancer confirmed by pathology in our institution from November 2018 to January 2021 were selected. The inclusion criteria were as follows: (1) axillary US and ABUS examinations were performed before biopsy or resection; (2) the ALN status of the patients was clearly verified by pathology after SLNB/ALND; (3) breast lesion with a diameter less than 5 cm (stages T1 and T2).

The exclusion criteria were as follows: (1) no complete clinicopathological data or axillary US and ABUS images; (2) the patient had undergone anticancer therapy (radiotherapy or chemotherapy); (3) the ABUS image quality was poor with artifacts.

All patients included in this study were randomly divided into the training and test sets at a ratio of 7:3.

Data acquisition

ABUS examination was performed by two well-trained technologists using the Invenia™ automatic breast ultrasound system with automatic 6 ~ 14 MHz linear broadband transformers (covering a volume of 15.4 × 17.0 × 5.0 cm) (Invenia™ ABUS, automatic breast ultrasound system, Ge Healthcare, Sunnyvale, the United States). The converter could be automatically moved in the scanning box. The thickness of each frame was 0.5 mm, and 330 images were collected axial. Patients were supine position with their arms fully raised to expose the breast. Meanwhile, a wedge-shaped cushion was placed under one side of the body that helps keep the breast stable with the nipple pointing towards the ceiling. A hypoallergenic lotion was distributed evenly over the breast with an additional amount over the area of the nipple. A disposable membrane was used to aid coupling and to uniformly compress the entire breast, enabling greater penetration, improving detail resolution at depth, and eliminating the creation of artifacts at the periphery. After the image collection was completed, all images were sent to the workstation for 3D reconstruction to obtain crown, horizontal and sagittal surface images. Radiologists can read all the image information in the workstation at any time for diagnosis. For the target tumor, diameter was measured as the largest diameter found on the axial plane of ABUS. ABUS imaging features were synthetically analyzed in three planes (axial, coronal and sagittal) using the Breast Imaging-Reporting and Data System (BI-RADS) lexicon, including margin (smooth, spiculated, angular or circumscribed), shape (regular or irregular), echo pattern (hypoechoic or complex), posterior acoustic feature (no change, enhance or decrease), calcification (no, macro or micro), orientation (horizontal or vertical), hyperechoic halo, retraction phenomenon. The above assessment was done by two radiologists (radiologist 1 and radiologist 2, with 9 and 15 years of experience in breast US) who were blinded the clinical and pathological information, and any differences were resolved through consultation.

The US-reported ALN status was obtained from the US reports, and axillary images including important features of suspicious lymph nodes were documented into the Picture Archiving and Communication Systems (PACS). It was retrospectively reviewed and verified by two radiologists (radiologist 1 and radiologist 2, with 9 and 15 years of experience in breast US). Axillary US features of lymph nodes used to assess suspicion for malignancy were as follows: (1) cortical thickness of 3 mm or greater; (2) longest/shortest axis ratio < 2; (3) absence of fatty hilum [21].

Tumor segmentation and radiomics feature extraction

The three-dimensional region of interest (3D-ROI) was manually drawn around the boundary of the mass on the axial ABUS images by a radiologist with 5 years of experience in breast US using the SEG3D2 software (https://www.sci.utah.edu/cibc-software/seg3d.html).

The first-order statistics, textural and wavelet features were extracted automatically from each ABUS image by pyradiomics (https://pyradiomics.readthedocs.io/en/latest/index.html). All the extracted features were in concordance with the standard set by the Imaging Biomarker Standardization Initiative (IBSI) [22]. Next, all radiomic features were rescheduled using Z-score normalization to facilitate subsequent statistical analysis. Figure 1 shows the flowchart of the radiomics score workflow and study flowchart.

Fig. 1
figure 1

Radiomics score workflow and study flowchart. Tumor was manually drawn ROI from the axis ABUS images by using the SEG3D2 software. Next, 837 features, including first-order statistics, textural and wavelet features, were extracted by pyradiomics. ICC > 0.75, correlation coefficient ≥ 0.9, tenfold cross-test, and the LASSO regression were applied to data dimension reduction and select the most significant ALN tumor burden related radiomics features. Univariable and multivariate logistic regression analysis was used to develop the predicting model. We incorporated all independent predictors, and this was presented as ABUS radiomics nomogram

To assess the inter-observer reproducibility, another radiologist (with 3 years of experience in breast US) drew 3D-ROI from 60 randomly chosen images. Intraclass correlation coefficient (ICC) was used to assess the inter-observer agreement, which was graded as excellent (0.90 to 1.00), good (0.75 to 0.90), moderate (0.50 to 0.75),or poor (< 0.50) [23]. The stable features with ICC > 0.75 were selected to adopt different segmentations.

Radiomics score

We used Pearson correlation coefficient (normally distributed data) or Spearman’s rank correlation coefficient (nonnormal or rank data) to evaluate the redundancy of the features and eliminated redundant features with correlation coefficient ≥ 0.9, with only the one left for next analysis. Then, the least absolute shrinkage and selection operator (LASSO) regression using tenfold cross-test was applied to select the most significant ALN tumor burden related radiomics features.

Model construction

Univariate and multivariate logistic regression analyses were performed to select the significant factors for ALN tumor burden. In univariate analysis, factors having P values < 0.10 were included in the multivariate analysis. Then, factors having P values < 0.05 were considered independent predictors after the multivariate analysis. Finally, ABUS model and radiomics nomogram was developed by incorporating these independent predictors.

Model validation

In this study, the validity of the prediction model was assessed through Receiver operating characteristic (ROC) curve, calibration curve, and decision curve.

ROC curves were plotted to assess the performance of the prediction model for ALN tumor burden in the training and test sets. The relevant metrics, including areas under the curve (AUC), sensitivity, specificity, accuracy, positive predictive value (PPV), and negative predictive value (NPV), were also calculated.

Calibration curves were plotted to explore the predictive accuracy of the radiomics nomogram in the training and test sets. Besides, the goodness-of-fit was evaluated with the Hosmer–Lemeshow test.

To demonstrate the clinical value of the radiomics nomogram decision curves were drawn.

Statistic analysis

Statistical analyses were performed using R software (version 4.1.2, https://www.r-project.org/). Categorical variables were compared using the chi-square test, and continuous variables were compared using the t-test or Man-Whitney U test to evaluate the consistency of the factors in the training and test sets. The reported statistical significance levels were all two-sided, and p values less than 0.05 were considered statistically significant. R software was used to construct and assess the radiomics score and the prediction model (details shown in Additional file 1).

Results

Basic information

A total of 310 patients were enrolled and randomly divided into a training set (n = 213) and a test set (n = 97) in this study. Table 1 showed the characteristics of breast cancer in training and test sets. There were 41(19.2%) and 17(17.5%) patients with ALN high tumor burden of early breast cancer in the training and test sets. There were no significant differences between the sets in the characteristics (p > 0.05).

Table 1 Characteristics of breast cancer in training and test sets

Radiomics score

A total of 837 features (Additional file 2) were extracted from each patient’s axial ABUS images, including first-order statistics (n = 18), texture features (n = 75) and wavelet features (n = 738). First, in the evaluation of reproducibility, 508 showing ICC more than 0.75 were used for subsequent analysis. Second, we eliminated redundant features with correlation coefficient ≥ 0.9, and only one was chosen. Third, after the LASSO regression and tenfold cross-test, 13 radiomics features with nonzero coefficients were chosen (Fig. 2). Then, we got the final formula for radiomics score.

Fig. 2
figure 2

Radiomics features selection using tenfold cross-test and LASSO regression. Tuning parameter(λ) selection in the LASSO regression used tenfold cross-test based on the minimum criterion. Dotted vertical lines indicated the optimal values using the minimum criteria and the1-SE criteria. A λ value of 0.024(log(λ) = -3.715) was chosen (minimum criterion) according to tenfold cross-test. B A coefficient profile plot was produced against the log(λ) sequence. Dotted vertical lines indicated the value obtained by the above tenfold cross-test, which resulted in 13 radiomics features with nonzero coefficients. C After tenfold cross-test and LASSO regression, the name and coefficient of selected radiomics features showed by a bar diagram

$$\text{radiomics score= -1.57486+wavelet.LLH\_glcm\_Imc1* 0.06655}{ } + \text{wavelet.LLH\_gldm\_DependenceVariance*-0.12317} \, \text{+wavelet.LHL\_glcm\_JointEnergy*-0.05780}$$
$$\text{+wavelet.LHH\_glcm\_JointEnergy*-0.24198}$$
$$\text{+wavelet.LHH\_glcm\_MaximumProbability *-0.09588}$$
$$\text{+wavelet.HLL\_firstorder\_Kurtosis *0.23917}$$
$$\text{+wavelet.HLL\_glcm\_Imc2*0.06358}$$
$$\text{+wavelet.HLL\_glszm\_HighGrayLevelZoneEmphasis*0.60089}$$
$$\text{+wavelet.HLL\_gldm\_SmallDependenceHighGrayLevelEmphasis*0.03775}$$
$$\text{+wavelet.HHL\_firstorder\_Skewness *0.21118}$$
$$\text{+wavelet.HHL\_glcm\_SumEntropy*-0.28907}$$
$$\text{+wavelet.HHH\_glcm\_ClusterShade*-0.09944}$$
$$\text{+wavelet.LLL\_glrlm\_RunEntropy *0.09735}$$

The radiomics score of the high and low tumor burden was calculated using 13 radiomic features, respectively. The violin plot (Fig. 3) showed that radiomics score in the high tumor burden was significantly higher than low tumor burden (P < 0.001).

Fig. 3
figure 3

Distribution of radiomics score in high and low tumor burden patients. The patients with high tumor burden had significantly higher score than those with low tumor burden (P < 0.001)

Model construction

Univariable logistic regression analysis results showed in Table 2. On the multivariable analysis (Table 3) that included ABUS imaging features (the ABUS model), diameter (odds ratio [OR], 2.03, 95% CI: [1.292,3.256], P = 0.002), hyperechoic halo (OR, 2.48,95% CI: [1.149,5.388], P = 0.020), and retraction phenomenon (OR, 6.53, 95% CI: [2.972,15.542], P < 0.001) were independently associated with ALN tumor burden. On the multivariable analysis (Table 3) that included both clinicopathologic features, ABUS imaging features and radiomics score (the radiomics model), retraction phenomenon (OR, 3.34, 95% CI: [1.413,8.229], P = 0.006), US-reported ALN status (OR, 7.62, 95% CI: [2.429,33.904], P = 0.002), and radiomics score (OR, 3.82, 95% CI: [2.175.7.400], P < 0.001) were independently associated with ALN tumor burden.

Table 2 Results of univariable analysis in the training set
Table 3 Comparison of the multivariable models for ALN tumor burden in the training set

The radiomics nomogram that incorporated the above independent predictors was developed and presented as the nomogram (Fig. 4).

Fig. 4
figure 4

Developed ABUS radiomics nomogram. The ABUS radiomics nomogram was developed in the training set, incorporating the radiomics score, US-reported ALN status and retraction phenomenon. US, ultrasound; ALN, axillary lymph node

Model validation

Figure 5 A and B showed the performance of radiomics score, ABUS model, and radiomics nomogram for discriminating ALN status. The results of AUCs for radiomics score, ABUS model and radiomics nomogram were 0.794(95%CI:0.709,0.879), 0.772(95%CI:0.677,0.868),0.876(95%CI:0.815,0.937) in the training set and 0.789(95%CI:0.657,0.921),0.736(95%CI:0.602,0.870),0.851(95%CI:0.738,0.964) in the test set. Moreover, the performance of the radiomics nomogram was significantly more excellent than radiomics score and ABUS model both in the training and test sets (DeLong test P < 0.05). The diagnostic performance of radiomics score, ABUS model, and radiomics nomogram in the training and test sets was shown in Table 4. Furthermore, 200 times five-fold cross validation was performed to prove the robustness of radiomics nomogram in the training set, with a mean AUC of 0.863, a mean sensitivity of 0.861, a mean specificity of 0.831, and a mean accuracy of 0.839.

Fig. 5
figure 5

ROC Curves and Calibration curves of the model in the training and test sets. A ROC Curves of radiomics nomogram (blue, AUC:0.876), radiomics score (red, AUC:0.794), and ABUS model (black, AUC:0.772) in the training set. B ROC Curves of the radiomics nomogram (blue, AUC:0.851), radiomics score (red, AUC:0789), and ABUS model (black, AUC:0.736) in the test set. C, D Calibration curves of radiomics nomogram in the training (C) and test set (D). ROC, Receiver operating characteristic; US, ultrasound; ALN, axillary lymph node

Table 4 The diagnostic performance of radiomics score, ABUS model, and radiomics nomogram in the training and test sets

The Calibration curves (Fig. 5 C and D) of the radiomics nomogram in the training and test sets showed an accurate agreement between the prediction of ALN tumor burden and pathological verification. The Hosmer–Lemeshow test showed X-squared = 5.926, P = 0.655(P > 0.05) in the training set and X-squared = 11.856, P = 0.158(P > 0.05) in the test set. It meant there is no significant difference between the predicted result and the actual outcome.

Decision curves (Fig. 6) were plotted to assess the clinical value of the different models in the training set and test sets. If the threshold probabilities are between 0.1 to 1, the radiomics nomogram will receive maximized net benefit.

Fig. 6
figure 6

Decision curve of the radiomics nomogram (red line), radiomics score (purple line), ABUS model (blue line) and US-reported ALN status (yellow line) in the training set (A) and test (B) set. The vertical axis indicates the net benefit, the x-axis indicates the threshold probability. The black line indicates the presume that no patients showed ALN high burden, and the grey line indicates the presume that all patients showed ALN high burden

To evaluate the value of the ABUS radiomics nomogram in making the optimal treatment strategies. We compared the performance of US-reported ALN status with radiomics nomogram in predicting ALN tumor burden (confusion matrix shown in Additional file 3). In our study, the false negative rates of radiologist and radiomics nomogram were 7.3% (3/41), 14.6% (6/41) in the training set, and 11.7% (2/17),23.5% (4/17) in the test set. The false positive rates of radiologist and radiomics nomogram were 51.1% (88/172),17.4% (30/172) in the training set, and 43.8% (35/80),17.5% (14/80) in the test set.

Discussion

In the study, we constructed radiomics score based on ABUS image for predicting ALN tumor burden, which AUC values were 0.794 in the training set and 0.789 in the test set. The ABUS model, integrating diameter, hyperechoic halo, and retraction phenomenon, showed better performance, with AUC values of 0.772 in the training set and 0.736 in the test set. By multivariate logistic regression analysis, we developed a radiomics nomogram integrating the US-reported ALN status, retraction phenomenon and radiomics score, which showed best performance with AUC values 0.876 in the training set and 0.851 in the test set.

ABUS is a breakthrough in breast ultrasound that has been developed to meet the limitations of HHUS. ABUS provides three-dimensional imaging of the entire breast with multiplanar reconstructions, which has been demonstrated to improve diagnostic accuracy [24]. Several studies [11, 25,26,27,28,29,30,31] have focused on the ability of diagnosing breast lesions and assessment of response to neoadjuvant chemotherapy by ABUS. These studies indicated good diagnostic performance of ABUS. Recently, several studies [12, 32] showed some ABUS imaging features were significantly correlated with ALN status. In our study, ABUS model displayed adequate discriminative ability (AUC 0.772 in the training set and 0.736 in the test set). Three ABUS imaging features, diameter, hyperechoic halo, and retraction phenomenon, were obviously associated with ALN tumor burden (all P < 0.05), in this study. Some studies [33] reported that the size of primary breast cancer was significantly related to ALN status. Larger breast cancers have a wider range of glands invaded by cancer cells, and the probability of ALN metastasis through lymphatic drainage is also higher. In our study, retraction phenomenon was identified to be the strongest independent predictor with high diagnostic performance in differentiating high and low ALN tumor burden(P < 0.001). Desmoplastic reaction of breast malignancy, which can produce contraction of the surrounding tissues toward the mass and disrupt normal parallel tissue planes, might help explain the generation of retraction phenomenon [34]. The appearance of the retraction phenomenon indicates that breast cancer is highly aggressive and has poor response to treatment [35]. Jiang et al. [36] and Tang et al. [37]showed that the smaller and more superficial invasive carcinomas with lower histological grades tended to present with retraction phenomenon. A hyperechoic halo, also known as converging pattern, is caused by the compressed fibrous surrounding tissue or the infiltration between the tumor and the surrounding tissue [38, 39]. It reflects the degree of invasion of cancer cells and may be an important indicator of poor prognosis. Similarly, Tang et al. [37] reported that the malignant masses were associated with retraction phenomenon and discontinuous hyper-and hypoechoic rim (p < 0.001 for each).

However, visual assessment and qualitative descriptions for US features extremely rely on the personal clinical experience and subjective judgment of radiologists. Gillies et al. [40] have illustrated that tumor characteristics at the genetic and cellular levels can be captured from medical images by extracting and computing high-throughput features. Qiu et al. [41] extracted radiomics features from B-mode ultrasound and integrated the significant clinical characteristics of patients to construct radiomics model for predicting ALN metastasis. Therefore, this study comprised the radiomics signature and US-reported ALN status, with AUC values of 0.816 in the training cohort and 0.759 in the validation cohort. Similarly, Gao et al. [42] established a nomogram based on radiomics analysis of primary cancer B-mode ultrasound for predicting ALN tumor burden, with AUC values of 0.846 in the training cohort and 0.733 in the validation cohort. The results of above research were similar to our results. Jiang et al. [43] extracted radiomics features from shear-wave elastography and B-mode ultrasound and integrated Clinical characteristics of patients to construct radiomics model. The research showed in the performance of discriminating disease-free axillary (N0) and any axillary metastasis (N + (≥ 1)), it achieved a C-index of 0.845 for the training cohort and 0.817 for the validation cohort. The tool could also discriminate between low (N + (1–2)) and heavy metastatic ALN burden (N + (≥ 3)), with a C-index of 0.827 in the training cohort and 0.810 in the validation cohort.

ABUS provides standardized scanning protocols and uncouples detection from image acquisition, hence improving reproducibility, reducing operator dependency and radiologist workload [10]. To our knowledge, there is no published study that has showed whether the value of radiomics score from ABUS would benefit prediction of ALN tumor burden in early breast cancer. In our study, 13 radiomics features from axial images were screened out and used to established radiomics score. The ABUS radiomics score showed moderate discriminative ability (AUC 0.794 in the training set, 0.789 in the test set).

We constructed a radiomics nomogram that integrated radiomics score with retraction phenomenon and US-reported ALN status to improve its predictive accuracy for ALN tumor burden. Based on the proposed risk classifier, the ABUS radiomics nomogram was able to classify patients into low- and high-tumor burden groups. According to ACOSOG Z0011 trials, we should pay more attention to the ALN tumor burden, instead of metastasis or non-metastasis. The radiomics nomogram resulted in 17.4% (30/172) and 17.5% (14/80) false positives in the training and test sets, which was obviously below the rates in US-reported ALN status (51.1% (88/172), 43.8% (35/80) in the training and test sets). The results showed that the radiomics nomogram could evidently reduce the false positive rate, compared with US-reported ALN status. Accurate forecasting of ALN high tumor burden can assist in identifying what kind of initial axillary surgery can overlook the SLNB and undergo ALND specifically and assist in employing neoadjuvant chemotherapy or individualized adjuvant radiotherapy [44]. Besides, for patients with low ALN tumor burden, unnecessary treatment, and potential complications due to surgery can be avoided. Figure 7 showed an example, the prediction of ALN tumor burden in patient.

Fig. 7
figure 7

A case of radiomics nomogram. A 77-year-old woman who has a 3.1 cm diameter lesion, radiomics score = -1.301, retraction phenomenon positive and US-reported ALN positive, indicates a low burden, with a low probability of less than 45%. Pathology confirmed only one metastatic ALN in the patient. A coronal plane of ABUS examination. B axial plane of ABUS examination. C sagittal plane of ABUS examination. D Axillary US examination revealed suspiciously positive ALN (cortical thickening and lymphatic hilum disappeared). E The nomogram showed a low probability of high burden (< 45%), indicating ALN low tumor burden

Zhou et al. [45] firstly used deep learning methods in the evaluation of ALN status. Which included 974 primary tumor images of 756 patients with lymph node-negative early breast cancer, and validated the model in an independent validation set of 78 patients, with an AUC of 0.89. However, the logic behind the features and decisions based on deep learning is still a “black box”, which means the results based on deep learning are difficult to interpret [46]. On the contrary, in our study, “wavelet.LHH_glcm_JointEnergy”, “wavelet.HLL_firstorder_Kurtosis”, “wavelet.HHL_firstorder_Skewness”, and “wavelet.HHL_glcm_SumEntropy” were obviously associated with ALN tumor burden, indicating that these features strongly reflect ALN tumor burden. In our study, most of selected radiomics features are texture features after wavelet transform. Wavelet transform, which measures the resolution of image signals in different time, space, and frequency scale planes, is extremely useful for replaying even subtle but important texture information that is ignored by radiologists in low-contrast US images. Previous researches [19, 47] have demonstrated the texture features after wavelet transform are used to construct a prediction model.

There are several limitations to our study. First, the study was a retrospective study, inevitably there existed sample bias. Second, the sample originated from a single center and lacked external validation. Therefore, the number of positive patients was 58 (58/310,41/213 in the training set and 17/97 in the test set). The difference between the two sets may be caused by the small number of positive patients. In the future, more patients (particularly more positive patients), multi-center samples, and external validation are needed. Third, the three-dimensional ROI was manually delineated on the axis ABUS image by the radiologist. However, manually delineating the ROI was extremely time-consuming and inevitably involved inter-observer variability. In the future, we should research semi-automatic or automatic segmentation algorithms to overcome the above problems. Finally, we extracted radiomics features from intratumoral regions and failed to exploit peritumoral radiomics features in the study. The radiomics features of intratumoral and peritumoral regions should be integrated into further research.

Conclusions

The ABUS radiomics nomogram showed favorable ability for predicting ALN tumor burden, which may provide additional benefits in treatment strategies for patients with early breast cancer, especially for patients with low tumor burden. With a more individualized and precise assessment for ALN tumor burden, radiomics nomogram will assist clinicians to make the optimal treatment strategy and avoid overtreatment.

Availability of data and materials

The datasets used and analysed during the current study available from the corresponding author on reasonable request.

Abbreviations

ABUS:

Automated Breast Ultrasound

ALN:

Axillary lymph node

ALND:

Axillary lymph node dissection

US:

Ultrasound

SLN:

Sentinel lymph node

SLNB:

Sentinel lymph node biopsy

HHUS:

Handheld ultrasound

3D-ROI:

Three-dimensional region of interest

ICC:

Intraclass correlation coefficient

BI-RADS:

Breast Imaging-Reporting and Data System

LASSO:

Least absolute shrinkage and selection operator

ROC:

Receiver operating characteristic

AUC:

Areas under the receiver operating characteristic curve

PPV:

Positive predictive value

NPV:

Negative predictive value

OR:

Odds ratio

CI:

Confidence interval

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  2. Cardoso F, Kyriakides S, Ohno S, Penault-Llorca F, Poortmans P, Rubio IT, Zackrisson S, Senkus E, clinicalguidelines@esmo.org EGCEa. Early breast cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-updagger. Ann Oncol. 2019;30(8):1194–220.

    Article  CAS  PubMed  Google Scholar 

  3. Maxwell F, de Margerie MC, Bricout M, Cauderlier E, Chapelier M, Albiter M, Bourrier P, Espie M, de Kerviler E, de Bazelaire C. Diagnostic strategy for the assessment of axillary lymph node status in breast cancer. Diagn Interv Imaging. 2015;96(10):1089–101.

    Article  CAS  PubMed  Google Scholar 

  4. Giuliano AE, Ballman KV, McCall L, Beitsch PD, Brennan MB, Kelemen PR, Ollila DW, Hansen NM, Whitworth PW, Blumencranz PW, et al. Effect of Axillary Dissection vs No Axillary Dissection on 10-Year Overall Survival Among Women With Invasive Breast Cancer and Sentinel Node Metastasis: The ACOSOG Z0011 (Alliance) Randomized Clinical Trial. JAMA. 2017;318(10):918–26.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Lyman GH, Somerfield MR, Giuliano AE. Sentinel Lymph Node Biopsy for Patients With Early-Stage Breast Cancer: 2016 American Society of Clinical Oncology Clinical Practice Guideline Update Summary. J Oncol Pract. 2017;13(3):196–8.

    Article  PubMed  Google Scholar 

  6. Feng Y, Huang R, He Y, Lu A, Fan Z, Fan T, Qi M, Wang X, Cao W, Wang X, et al. Efficacy of physical examination, ultrasound, and ultrasound combined with fine-needle aspiration for axilla staging of primary breast cancer. Breast Cancer Res Treat. 2015;149(3):761–5.

    Article  PubMed  Google Scholar 

  7. Ahmed M, Jozsa F, Baker R, Rubio IT, Benson J, Douek M. Meta-analysis of tumour burden in pre-operative axillary ultrasound positive and negative breast cancer patients. Breast Cancer Res Treat. 2017;166(2):329–36.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Raber BM, Lin H, Shen Y, Shaitelman SF, Bedrosian I. Trends in Regional Nodal Management of Breast Cancer Patients with Low Nodal Burden. Ann Surg Oncol. 2019;26(13):4346–54.

    Article  PubMed  Google Scholar 

  9. Ahmed M, Douek M. Is axillary ultrasound imaging necessary for all patients with breast cancer? Br J Surg. 2018;105(8):930–2.

    Article  CAS  PubMed  Google Scholar 

  10. Vourtsis A. Three-dimensional automated breast ultrasound: Technical aspects and first results. Diagn Interv Imaging. 2019;100(10):579–92.

    Article  CAS  PubMed  Google Scholar 

  11. Zhang X, Lin X, Tan Y, Zhu Y, Wang H, Feng R, Tang G, Zhou X, Li A, Qiao Y. A multicenter hospital-based diagnosis study of automated breast ultrasound system in detecting breast cancer among Chinese women. Chin J Cancer Res. 2018;30(2):231–9.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Wang Q, Li B, Liu Z, Shang H, Jing H, Shao H, Chen K, Liang X, Cheng W. Prediction model of axillary lymph node status using automated breast ultrasound (ABUS) and ki-67 status in early-stage breast cancer. BMC Cancer. 2022;22(1):929.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Lin X, Wang J, Han F, Fu J, Li A. Analysis of eighty-one cases with breast lesions using automated breast volume scanner and comparison with handheld ultrasound. Eur J Radiol. 2012;81(5):873–8.

    Article  PubMed  Google Scholar 

  14. van Zelst JCM, Mann RM. Automated Three-dimensional Breast US for Screening: Technique, Artifacts, and Lesion Characterization. Radiographics. 2018;38(3):663–83.

    Article  PubMed  Google Scholar 

  15. Zheng FY, Yan LX, Huang BJ, Xia HS, Wang X, Lu Q, Li CX, Wang WP. Comparison of retraction phenomenon and BI-RADS-US descriptors in differentiating benign and malignant breast masses using an automated breast volume scanner. Eur J Radiol. 2015;84(11):2123–9.

    Article  PubMed  Google Scholar 

  16. Luo WQ, Huang QX, Huang XW, Hu HT, Zeng FQ, Wang W. Predicting Breast Cancer in Breast Imaging Reporting and Data System (BI-RADS) Ultrasound Category 4 or 5 Lesions: A Nomogram Combining Radiomics and BI-RADS. Sci Rep. 2019;9(1):11921.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Yu FH, Wang JX, Ye XH, Deng J, Hang J, Yang B. Ultrasound-based radiomics nomogram: A potential biomarker to predict axillary lymph node metastasis in early-stage invasive breast cancer. Eur J Radiol. 2019;119:108658.

    Article  PubMed  Google Scholar 

  18. Huang YQ, Liang CH, He L, Tian J, Liang CS, Chen X, Ma ZL, Liu ZY. Development and Validation of a Radiomics Nomogram for Preoperative Prediction of Lymph Node Metastasis in Colorectal Cancer. J Clin Oncol. 2016;34(18):2157–64.

    Article  PubMed  Google Scholar 

  19. Yang M, Liu H, Dai Q, Yao L, Zhang S, Wang Z, Li J, Duan Q. Treatment Response Prediction Using Ultrasound-Based Pre-, Post-Early, and Delta Radiomics in Neoadjuvant Chemotherapy in Breast Cancer. Front Oncol. 2022;12:748008.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong EEC, van Timmeren J, Sanduleanu S, Larue R, Even AJG, Jochems A, et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat Rev Clin Oncol. 2017;14(12):749–62.

    Article  PubMed  Google Scholar 

  21. Koelliker SL, Chung MA, Mainiero MB, Steinhoff MM, Cady B. Axillary lymph nodes: US-guided fine-needle aspiration for initial staging of breast cancer–correlation with primary tumor size. Radiology. 2008;246(1):81–9.

    Article  PubMed  Google Scholar 

  22. Zwanenburg A, Vallieres M, Abdalah MA, Aerts H, Andrearczyk V, Apte A, Ashrafinia S, Bakas S, Beukinga RJ, Boellaard R, et al. The Image Biomarker Standardization Initiative: Standardized Quantitative Radiomics for High-Throughput Image-based Phenotyping. Radiology. 2020;295(2):328–38.

    Article  PubMed  Google Scholar 

  23. Koo TK, Li MY. A Guideline of Selecting and Reporting Intraclass Correlation Coefficients for Reliability Research. J Chiropr Med. 2016;15(2):155–63.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Van Zelst JC, Platel B, Karssemeijer N, Mann RM. Multiplanar Reconstructions of 3D Automated Breast Ultrasound Improve Lesion Differentiation by Radiologists. Acad Radiol. 2015;22(12):1489–96.

    Article  PubMed  Google Scholar 

  25. Wang X, Huo L, He Y, Fan Z, Wang T, Xie Y, Li J, Ouyang T. Early prediction of pathological outcomes to neoadjuvant chemotherapy in breast cancer patients using automated breast ultrasound. Chin J Cancer Res. 2016;28(5):478–85.

    Article  PubMed  PubMed Central  Google Scholar 

  26. D’Angelo A, Rinaldi P, Belli P, D’Amico R, Carlino G, Grippo C, Giuliani M, Orlandi A, Infante A, Manfredi R. Usefulness of automated breast volume scanner (ABVS) for monitoring tumor response to neoadjuvant treatment in breast cancer patients: preliminary results. Eur Rev Med Pharmacol Sci. 2019;23(1):225–31.

  27. Guldogan N, Yilmaz E, Arslan A, Kucukkaya F, Atila N, Aribal E. Comparison of 3D-Automated Breast Ultrasound With Handheld Breast Ultrasound Regarding Detection and BI-RADS Characterization of Lesions in Dense Breasts: A Study of 592 Cases. Acad Radiol. 2022;29(8):1143–48.

  28. Zhang L, Bao LY, Tan YJ, Zhu LQ, Xu XJ, Zhu QQ, Shan YN, Zhao J, Xie LS, Liu J. Diagnostic Performance Using Automated Breast Ultrasound System for Breast Cancer in Chinese Women Aged 40 Years or Older: A Comparative Study. Ultrasound Med Biol. 2019;45(12):3137–44.

    Article  PubMed  Google Scholar 

  29. Vourtsis A, Kachulis A. The performance of 3D ABUS versus HHUS in the visualisation and BI-RADS characterisation of breast lesions in a large cohort of 1,886 women. Eur Radiol. 2018;28(2):592–601.

    Article  PubMed  Google Scholar 

  30. Choi WJ, Kim SH, Shin HJ, Bang M, Kang BJ, Lee SH, Chang JM, Moon WK, Bae K, Kim HH. Automated breast US as the primary screening test for breast cancer among East Asian women aged 40–49 years: a multicenter prospective study. Eur Radiol. 2021;31(10):7771–82.

    Article  PubMed  Google Scholar 

  31. Yang S, Gao X, Liu L, Shu R, Yan J, Zhang G, Xiao Y, Ju Y, Zhao N, Song H. Performance and Reading Time of Automated Breast US with or without Computer-aided Detection. Radiology. 2019;292(3):540–9.

    Article  PubMed  Google Scholar 

  32. Zhao F, Cai C, Liu M, Xiao J. Identification of the lymph node metastasis-related automated breast volume scanning features for predicting axillary lymph node tumor burden of invasive breast cancer via a clinical prediction model. Front Endocrinol (Lausanne). 2022;13:881761.

    Article  PubMed  Google Scholar 

  33. Schwartz RS, Erban JK. Timing of Metastasis in Breast Cancer. N Engl J Med. 2017;376(25):2486–8.

    Article  PubMed  Google Scholar 

  34. Kotsianos-Hermle D, Hiltawsky KM, Wirth S, Fischer T, Friese K, Reiser M. Analysis of 107 breast lesions with automated 3D ultrasound and comparison with mammography and manual ultrasound. Eur J Radiol. 2009;71(1):109–15.

    Article  CAS  PubMed  Google Scholar 

  35. Zheng FY, Lu Q, Huang BJ, Xia HS, Yan LX, Wang X, Yuan W, Wang WP. Imaging features of automated breast volume scanner: Correlation with molecular subtypes of breast cancer. Eur J Radiol. 2017;86:267–75.

    Article  PubMed  Google Scholar 

  36. Jiang J, Chen YQ, Xu YZ, Chen ML, Zhu YK, Guan WB, Wang XJ. Correlation between three-dimensional ultrasound features and pathological prognostic factors in breast cancer. Eur Radiol. 2014;24(6):1186–96.

    Article  PubMed  Google Scholar 

  37. Tang G, An X, Xiang H, Liu L, Li A, Lin X. Automated Breast Ultrasound: Interobserver Agreement, Diagnostic Value, and Associated Clinical Factors of Coronal-Plane Image Features. Korean J Radiol. 2020;21(5):550–60.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Chen L, Chen Y, Diao XH, Fang L, Pang Y, Cheng AQ, Li WP, Wang Y. Comparative study of automated breast 3-D ultrasound and handheld B-mode ultrasound for differentiation of benign and malignant breast masses. Ultrasound Med Biol. 2013;39(10):1735–42.

    Article  PubMed  Google Scholar 

  39. Kalmantis K, Dimitrakakis C, Koumpis C, Tsigginou A, Papantoniou N, Mesogitis S, Antsaklis A. The contribution of three-dimensional power Doppler imaging in the preoperative assessment of breast tumors: a preliminary report. Obstet Gynecol Int. 2009;2009:530579.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Gillies RJ, Kinahan PE, Hricak H. Radiomics: Images Are More than Pictures. They Are Data Radiology. 2016;278(2):563–77.

    PubMed  Google Scholar 

  41. Qiu X, Jiang Y, Zhao Q, Yan C, Huang M, Jiang T. Could Ultrasound-Based Radiomics Noninvasively Predict Axillary Lymph Node Metastasis in Breast Cancer? J Ultrasound Med. 2020;39(10):1897–905.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Gao Y, Luo Y, Zhao C, Xiao M, Ma L, Li W, Qin J, Zhu Q, Jiang Y. Nomogram based on radiomics analysis of primary breast cancer ultrasound images: prediction of axillary lymph node tumor burden in patients. Eur Radiol. 2021;31(2):928–37.

    Article  PubMed  Google Scholar 

  43. Jiang M, Li CL, Luo XM, Chuan ZR, Chen RX, Tang SC, Lv WZ, Cui XW, Dietrich CF. Radiomics model based on shear-wave elastography in the assessment of axillary lymph node status in early-stage breast cancer. Eur Radiol. 2022;32(4):2313–25.

  44. Li B, Zhao X, Wang Q, Jing H, Shao H, Zhang L, Cheng W. Prediction of high nodal burden in invasive breast cancer by quantitative shear wave elastography. Quant Imaging Med Surg. 2022;12(2):1336–47.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Zhou LQ, Wu XL, Huang SY, Wu GG, Ye HR, Wei Q, Bao LY, Deng YB, Li XR, Cui XW, et al. Lymph Node Metastasis Prediction from Primary Breast Cancer US Images Using Deep Learning. Radiology. 2020;294(1):19–28.

    Article  PubMed  Google Scholar 

  46. Afshar P, Mohammadi A, Plataniotis KN, Oikonomou A, Benali H. From handcrafted to deep-learning-based cancer radiomics: challenges and opportunities. ISPM. 2019;36(4):132–60.

    Google Scholar 

  47. Guo Y, Hu Y, Qiao M, Wang Y, Yu J, Li J, Chang C. Radiomics Analysis on Ultrasound for Prediction of Biologic Behavior in Breast Invasive Ductal Carcinoma. Clin Breast Cancer. 2018;18(3):e335–44.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (grant number: 82171953) and HAI YAN Science Foundation of Harbin Medical University Cancer Hospital (grant number: JJMS2020-01).

Author information

Authors and Affiliations

Authors

Contributions

All authors have made substantial contributions to this manuscript:Y.C. and H.J. conceived the study; Q.W. and Y.X. contributed to imaging data curation; B.L. and Y.C. performed data processing and analysis; Z.N. and H.S. wrote the original draft; H.J. and Q.W. reviewed and edited the manuscript and supervised the project. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Qiucheng Wang or Hui Jing.

Ethics declarations

Ethics approval and consent to participate

The retrospective study was approved by the Institutional Review Board of Harbin medical university cancer hospital and was granted a waiver of written informed consent for use of data. Institutional Review Board of Harbin medical university cancer hospital waived informed consent from all patients due to retrospective nature of the study. All methods were performed in accordance with the relevant guidelines and regulations. Research involving human participants, human material, or human data have been performed in accordance with the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

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: S1.

The details of R software used in this study.

Additional file 2: Table S1.

The details of radiomic features.

Additional file 3: Figure S1.

The confusion matrix of Radiomics nomogram and US-reported ALN status in the training and test sets.

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, Y., Xie, Y., Li, B. et al. Automated Breast Ultrasound (ABUS)-based radiomics nomogram: an individualized tool for predicting axillary lymph node tumor burden in patients with early breast cancer. BMC Cancer 23, 340 (2023). https://doi.org/10.1186/s12885-023-10743-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-10743-3

Keywords