Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 10 March 2022
Sec. Cancer Imaging and Image-directed Interventions
This article is part of the Research Topic Advances and New Insights into Cancer Characterization: When Novel Imaging Meets Quantitative Imaging Biomarkers View all 24 articles

Comprehensive Risk System Based on Shear Wave Elastography and BI-RADS Categories in Assessing Axillary Lymph Node Metastasis of Invasive Breast Cancer—A Multicenter Study

Huiting Zhang&#x;Huiting Zhang1†Yijie Dong&#x;Yijie Dong1†Xiaohong JiaXiaohong Jia1Jingwen ZhangJingwen Zhang1Zhiyao LiZhiyao Li2Zhirui ChuanZhirui Chuan2Yanjun XuYanjun Xu3Bin HuBin Hu4Yunxia HuangYunxia Huang5Cai ChangCai Chang5Jinfeng XuJinfeng Xu6Fajin DongFajin Dong6Xiaona XiaXiaona Xia7Chengrong WuChengrong Wu7Wenjia HuWenjia Hu8Gang WuGang Wu8Qiaoying LiQiaoying Li9Qin ChenQin Chen10Wanyue DengWanyue Deng10Qiongchao JiangQiongchao Jiang11Yonglin MouYonglin Mou12Huannan YanHuannan Yan12Xiaojing XuXiaojing Xu13Hongju YanHongju Yan13Ping ZhouPing Zhou14Yang ShaoYang Shao14Ligang CuiLigang Cui15Ping HePing He15Linxue QianLinxue Qian16Jinping LiuJinping Liu16Liying ShiLiying Shi17Yanan ZhaoYanan Zhao18Yongyuan XuYongyuan Xu18Yanyan Song*Yanyan Song19*Weiwei Zhan*Weiwei Zhan1*Jianqiao Zhou*Jianqiao Zhou1*
  • 1Department of Ultrasound, Ruijin Hospital, Shanghai Jiaotong University School of Medicine, Shanghai, China
  • 2Department of Medical Ultrasound, Yunnan Cancer Hospital & The Third Affiliated Hospital of Kunming Medical University, Kunming, China
  • 3Department of Ultrasound in Medicine, Shanghai Jiao Tong University Affiliated Sixth People’s Hospital, Shanghai Institute of Ultrasound in Medicine, Shanghai, China
  • 4Department of Ultrasound, Minhang Hospital, Fudan University, Shanghai, China
  • 5Department of Ultrasonography, Fudan University Shanghai Cancer Center, Shanghai Medical College, Fudan University, Shanghai, China
  • 6Department of Ultrasound, Shenzhen People’s Hospital, The Second Clinical Medical College, Jinan University, and The First Affiliated Hospital, Southern University of Science and Technology, Shenzhen, China
  • 7Department of Ultrasound Medicine, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China
  • 8Department of Ultrasound, People’s Hospital of Henan Province, Zhengzhou, China
  • 9Department of Ultrasound Diseases, Tangdu Hospital, Four Military Medical University, Xi’an, China
  • 10Department of Ultrasound, Sichuan Provincial People’s Hospital, University of Electronic Science and Technology of China, Chengdu, China
  • 11Department of Ultrasound, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, China
  • 12Department of Ultrasound, General Hospital of Northern Theater Command, Shenyang, China
  • 13Department of Ultrasound, Affiliated Hangzhou First People’s Hospital, Zhejiang University School of Medicine, Hangzhou, China
  • 14Department of Ultrasound, The Third Xiangya Hospital of Central South University, Changsha, China
  • 15Department of Ultrasound, Peking University Third Hospital, Beijing, China
  • 16Department of Ultrasound, Beijing Friendship Hospital, Capital Medical University, Beijing, China
  • 17Department of Ultrasound, Affiliated Hospital of Guizhou Medical University, Guizhou, China
  • 18Department of Ultrasound, Second Affiliated Hospital of Zhejiang University, School of Medicine, Hangzhou, China
  • 19Department of Biostatistics, Institute of Medical Sciences, Shanghai Jiaotong University School of Medicine, Shanghai, China

Purpose: To develop a risk stratification system that can predict axillary lymph node (LN) metastasis in invasive breast cancer based on the combination of shear wave elastography (SWE) and conventional ultrasound.

Materials and Methods: A total of 619 participants pathologically diagnosed with invasive breast cancer underwent breast ultrasound examinations were recruited from a multicenter of 17 hospitals in China from August 2016 to August 2017. Conventional ultrasound and SWE features were compared between positive and negative LN metastasis groups. The regression equation, the weighting, and the counting methods were used to predict axillary LN metastasis. The sensitivity, specificity, and the areas under the receiver operating characteristic curve (AUC) were calculated.

Results: A significant difference was found in the Breast Imaging Reporting and Data System (BI-RADS) category, the “stiff rim” sign, minimum elastic modulusof the internal tumor and peritumor region of 3 mm between positive and negative LN groups (p < 0.05 for all). There was no significant difference in the diagnostic performance of the regression equation, the weighting, and the counting methods (p > 0.05 for all). Using the counting method, a 0–4 grade risk stratification system based on the four characteristics was established, which yielded an AUC of 0.656 (95% CI, 0.617–0.693, p < 0.001), a sensitivity of 54.60% (95% CI, 46.9%–62.1%), and a specificity of 68.99% (95% CI, 64.5%–73.3%) in predicting axillary LN metastasis.

Conclusion: A 0–4 grade risk stratification system was developed based on SWE characteristics and BI-RADS categories, and this system has the potential to predict axillary LN metastases in invasive breast cancer.

Introduction

Female breast cancer, the first cause of the death in malignant tumors (1), is the most common neoplasm in 20–59-year-old women; early diagnosis can reduce 40% of deaths (25). According to the American Joint Committee on Cancer, tumor size (T), node (N), and distant metastasis (M) have served as the global standard for conveying disease status among clinicians (6). Sentinel lymph node biopsy (SLNB) is regarded as the standard method for identifying the staging and determining the clinical arrangement. However, SLNB was recommended only for T1 or T2 tumors (7, 8), and multiple factors contribute to its limited pathological findings, suggesting that tedious and time-consuming SLNB is not optimal (9). Besides, SLNB may cause complications like arm numbness or upper limb edema in 3.5%–10.9% patients because of the increasing anesthesia time (10). Thus, an accurate and non-invasive examination for evaluating axillary LN status is expected before surgery or clinical treatment. However, showing a great difference among various studies, the sensitivity and specificity of axillary ultrasound in detection of LN metastasis ranged from 45.2% to 92.7% and from 40.5% to 93.9%, respectively (1114).

For breast ultrasound examinations, as the suspicious imaging features increase, the Breast Imaging Reporting and Data System (BI-RADS) categories of the breast tumor increase accordingly. Tumors in BI-RADS category 4 or 5 indicate a risk of malignancy, and a detailed assessment of the SLN status is required (15). Shear wave elastography (SWE) is a novel ultrasonography based on the speed of shear wave velocity in different tissues, which can provide quantitative (that is kPa or m/s) and qualitative information (such as “stiff rim” sign) in differentiation of breast tumors (10, 12, 16). It has been reported that LN metastasis of breast cancer was associated with breast lesion elastic modulus (17). Besides, the “stiff rim” sign, usually appearing in malignant breast lesions with increased peritumoral stiffness, was regarded as the infiltration of cancer cells into the interstitial tissues or a desmoplastic reaction (16), which is an independent prognostic factor predicting tumor recurrence and patient death (18, 19). Thus, it is reasonable to hypothesize that the “stiff rim” sign may be potentially associated with LN metastasis. However, to our best knowledge, no study has detected the relationship between SWE information of peritumor tissues and axillary LN metastasis.

Hence, our study aims at assessing the risk of breast cancer LN metastasis in ultrasound examination by combining quantitative and qualitative SWE features, as well as BI-RADS categories.

Materials and Methods

Patients

This prospective study was approved by the institutional ethics committee, and written informed consent to participate was acquired before examinations. From August 2016 to August 2017, there were a total of 689 patients with malignant breast tumors from a multicenter consisting of 17 hospitals in China diagnosed by surgical pathology enrollment, with conventional ultrasound and SWE images acquired before treatment. The pathology of LNs from SLNB or axillary LN dissection was acquired. All macrometastases, micrometastases, and isolated tumor cells were counted as node positive. The exclusion criteria were set as follows: (1) patients with non-invasive malignant tumors according to the pathology; (2) patients who accepted chemotherapy before surgery; and (3) patients who accepted chemotherapy or biopsy before ultrasound examination. In total, we have excluded 70 patients, including 34 ductal carcinoma in situ, 11 solid papillary carcinoma, 6 intraductal papillary carcinoma, 2 cases with neoadjuvant chemotherapy before surgery, 6 cases of secondary tumor, and 11 cases without LN pathology. According to the rules established by the study, only one lesion was evaluated per patient, i.e., the most suspicious lesion in ultrasound examination or the largest one among the same BI-RADS category was selected in patients with multiple masses. Finally, our study included 619 invasive breast cancer participants (age range, 22–91 y; mean age, 52.16 + 11.21 y), including 1 male and 618 female patients.

SWE Image Acquisition and Analysis

All the ultrasound examinations of patients from 17 hospitals used the Resona 7 ultrasound system (Mindray Medical International, Shenzhen, China) equipped with an L11-3 high-frequency probe. Prior to collecting data, all participating radiologists received systematic training in conventional ultrasound and SWE breast examination. Moreover, all radiologists from the multicenter had more than 3 years of experience in breast ultrasound elastography. Following standard conventional ultrasound examination, SWE was performed according to the already well-established method (20). The SWE display scale was set at 140 kPa according to the protocol provided by the manufacturer.

An online measurement of the tumor size and SWE parameters was done instantly by the examining radiologist. In measuring SWE parameters, the outline of the lesion was drawn by tracing, and the function of the “shell” was equipped to acquire the elastic modulus of the peritumoral tissues in regions of 1, 2, and 3 mm outside the boundary of internal tumors. Values of mean elasticity (Emean), maximum elasticity (Emax), minimum elastic modulus (Emin), standard deviation (Esd) evaluating both internal tumors and peritumoral tissues as well as intratumor ratio between mean elastic modulus of breast lesions, and normal fatty tissue (Eratio) were acquired and recorded. In assessing the location of the breast tumor, it is based on the location of the center of the lump in the 4 quadrants of the breast (upper/lower outer/inner quadrants). Upon completion of the evaluation, both B-mode images and SWE images of each patient were stored in a hard disk of the ultrasound system and subsequently sent to the study center (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1 B-mode and SWE images of a patient with invasive breast tumor. The “stiff rim” sign presents on the SWE images. (A) SWE images and parameters of outlined tumor. The Emax, Emean, Emin, and Esd of internal tumor were 34.63, 118.90, 3.60, and 16.82 kPa. (B) SWE images and parameters of shell 1 mm. The Emax, Emean, Emin, and Esd of shell 1 mm were 45.35, 202.76, 0.84, and 31.62 kPa. (C) SWE images and parameters of shell 2 mm. The Emax, Emean, Emin, and Esd of shell 2 mm were 48.00, 287.41, 0.84, and 36.93 kPa. (D) SWE images and parameters of shell 3 mm. The Emax, Emean, Emin, and Esd of shell 3 mm were 49.17, 296.74, 0.84, and 37.21 kPa.

At the study center, three radiologists with more than 8 years of experience in breast ultrasound evaluated ultrasound images according to consensus principles and proposed BI-RADS classification based on conventional ultrasound. At the same time, the presence of the “stiff rim” sign of the tumor was determined based on the SWE images. SWE images are colored in blue, green, orange, and red to sequentially show the progressive increase in tissue stiffness. When compared with the stiffness of the surrounding normal breast tissues and the interior tumor tissues, the peritumoral region, up to approximately 3 mm outside the boundary of internal tumors, showed increased stiffness (coded in orange or red), then it is judged as positive for the “stiff rim” sign (16).

Statistical Analysis

Statistical analyses were performed using SPSS Version 23.0 (IBM, Armonk, NY, USA) and SAS 9.4(Statistical Analysis System, v.9.4). A p value less than 0.05 was considered for significant differences. The Mann–Whitney U test was used to compare the statistical difference in continuous variables between positive and negative LN groups, and the chi-square test was used to verify whether there was a statistical difference in categorical variables between the two groups. Significant variables were included in the final logistic regression analysis. The selection criteria alpha was set to 0.05. The odds ratio (OR) and 95% confidence interval (CI) for each feature were obtained from the logistic regression model. To predict axillary LN metastasis, the regression equation, the weighting, and the counting methods were used. The OR for the BI-RADS category and each SWE feature from the final logistic regression model was used in the weighting and counting methods. In the weighting method, the ORs were standardized using rounding to ensure appropriate application. In the counting method, the number of features with OR >1 was counted (21). The area under the receiver operating characteristic curve (AUC), sensitivity, specificity, and 95% CI of different variables were calculated, with histologic diagnosis as the gold standard. Finally, the risk system was established by balancing the diagnostic performance and ease of use. Multivariate logistic analysis was used to evaluate the OR of each risk grade. The Cochran–Armitage test was used to determine the change in the probability of axillary LN metastasis with the risk grade increased.

Results

Pathologic Results

A total of 619 invasive breast cancers included 508 cases of invasive ductal carcinoma and 111 cases of other types (Table 1). Among them, 174 (28.1%) patients were diagnosed with positive LN metastasis and 445 (71.9%) were negative. The mean size of the tumor was 20.99 ± 9.04 mm (range 4.7–89 mm).

TABLE 1
www.frontiersin.org

Table 1 Pathologic diagnosis in 619 patients.

Clinical and BI-RADS Characteristics

Features of categorical variables including staging of tumor, tumor location, and tumor BI-RADS categories were compared between negative and positive LN groups, as shown in Table 2. There was a significant difference in constituent ratio between 5 BI-RADS categories (p < 0.001). In the 619 masses, 481 lesions were defined as BI-RADS 3-4C, with about 25% (ranging from 22.4% to 27.3%) cases in positive metastasis LN, while in 138 masses of BI-RADS 5, patients with positive LN metastasis were up to 43.5%. There was no significant difference in tumor staging (p = 0.340) and location (p = 0.390) in the two groups.

TABLE 2
www.frontiersin.org

Table 2 Clinical and ultrasound characteristics between negative and positive LN groups.

Qualitative SWE Characteristics

For the qualitative SWE feature, i.e., the “stiff rim” sign, there was a significant difference between the positive and negative groups for lymph node metastasis (p < 0.001). In patients with negative LN metastasis, 68.3% (304/445) cases presented a “stiff rim” sign in SWE images, while in positive LN groups, up to 82.2% (143/174) cases presented a “stiff rim” sign, significantly higher than that of negative LN groups.

Quantitative SWE Characteristics

The characteristics of continuous SWE parameters in internal tumors and peritumor tissues between the two groups are shown in Table 3. The Emean, Emax, Emin, Esd, and Eratio of internal tumors in positive LN groups were similar to that of negative LN groups (p > 0.05 for all). The Emean, Emax, and Esd of shell 1, 2, and 3 mm as well as Emin of shell 1 and 2 mm showed no significant difference between the two groups (p > 0.05 for all). However, the Emin of internal tumor (p = 0.005) and shell 3 mm (p = 0.007) was significantly lower in positive LN groups than that in the negative group. The AUC, sensitivity, specificity, 95% CI, and cutoff values are shown in Table 4. Only Emin of internal tumor (p <0.01) and shell 3 mm (p < 0.01) could predict axillary LN metastasis, with cutoff values of 3.31 and 3.52, respectively.

TABLE 3
www.frontiersin.org

Table 3 Quantitative evaluation of SWE in the internal and peritumor tissues between positive and negative LN metastasis groups.

TABLE 4
www.frontiersin.org

Table 4 Quantitative SWE evaluation of the tumor and peritumor shell for the predicting of axillary LN metastasis.

Risk Grade of LN Metastasis and Its Diagnostic Performance

The multiple logistic regression analysis indicated that the following four features were significantly associated with LN metastasis: (1) the breast tumor was assessed as BI-RADS 5; (2) the Emin of intrinsic tumor was no more than 3.31 kPa; (3) the Emin of shell 3 mm was no more than 3.52 kPa; and (4) the tumor showed a positive “stiff rim” sign in SWE. The OR values of the four features were 2.155, 1.654, 1.564, and 1.900, respectively (p < 0.05 for all). In the weighting method, the weight values of all four features were 2 based on the OR values of the above four features (Table 5). In the prediction of LN metastasis, the logistic regression equation yielded the highest AUC of 0.659 (95% CI, 0.621–0.697), while the weighting method and the counting method had the same AUC of 0.656 (95% CI, 0.617–0.693). There was no statistical difference among the three (p > 0.05 for all) (Table 6 and Figure 2).

TABLE 5
www.frontiersin.org

Table 5 Odds ratios for the suspicious features between positive and negative LN groups and the corresponding weighting values.

TABLE 6
www.frontiersin.org

Table 6 Comparison of the logistic regression equation, the weighting, and the counting methods in predicting the axillary LN metastasis.

FIGURE 2
www.frontiersin.org

Figure 2 The area under the receiver operating characteristic (ROC) curves of the logistic regression equation, the weighting method, and the counting method in predicting axillary lymph node metastasis of invasive breast cancer.

Following the principle of balancing the diagnostic performance and ease of use, the counting method was used to establish the risk system for predicting axillary LN metastases. Thus, a 0–4 grade risk stratification system was established based on the four suspicious features mentioned above: The LN metastasis risk grade of the breast cancer is defined from 0 to 4 according to the number of suspicious features, i.e., the mass was judged to be risk grade 4 when all suspicious features were present and grade 0 when no suspicious features were present. The risk grade of LN metastasis for the 619 patients was determined according to the above rule.

Subsequently, we compared the difference of risk grade between patients with positive and negative metastasis LNs (Table 7), and a significant difference was shown between the two groups (p < 0.001). The percentage of patients with positive LNs from grades 0 to 4 was constantly increased, with 7.1% (3/42) for grade 0, 19.2% (30/156) for grade 1, 24.5% (46/188) for grade 2, 35.6% (64/180) for grade 3, and 58.5% (31/53) for grade 4, respectively, which is shown in Figure 3. The Cochran–Armitage trend test indicated that the risk of axillary LN metastasis was increased with the increase in grade level (p < 0.0001). Results of the logistic analysis of the comprehensive grade showed that the OR of axillary LN metastasis was increased progressively with rising grade. The risk of axillary LN metastasis was 3.095 times higher in patients with risk grade 1 than in patients with risk grade 0 (p = 0.074), and 4.211–18.318 times higher in patients with risk grades 2–4 (p < 0.05 for all). The risk grade model predicted axillary LN metastasis with an AUC of 0.656 (95% CI, 0.617–0.693), a sensitivity of 54.60% (95% CI, 46.9%–62.1%), a specificity of 68.99% (95% CI, 64.5%–73.3%), and a cutoff value of grade 2.

TABLE 7
www.frontiersin.org

Table 7 Proportion and odds ratios for each risk grade in negative and positive LN metastasis groups.

FIGURE 3
www.frontiersin.org

Figure 3 Proportion distribution of positive and negative LN metastasis in each grade of the risk system.

Since risk category 1 was not different from category 0, we also compared the difference of OR values between grade 1 vs. grades 2, 3, and 4 and grades 0–1 vs. grades 2, 3, and 4 in the risk system. Results showed that the risk of axillary LN metastasis was higher in patients with risk grades 3 and 4 than in patients with risk grades 0–1 and grade 1 (p < 0.05 for all) (Table 8).

TABLE 8
www.frontiersin.org

Table 8 Differences of odds ratios between grade 1 and grades 2, 3, 4, as well as between grades 0–1 and grades 2, 3, 4.

Discussion

The risk of LN metastasis constantly exists in breast cancer patients, and it is expected to be provided in an early stage in order to guide the clinical decision and improve the tumor survival. SWE has been used to predict axillary LN metastasis in invasive breast cancer. Previous studies have reported that higher elasticity of the tumor was related to lymph node metastasis of breast cancer (22, 23), and further study has revealed that the mean elastic modulus of internal breast tumor could independently predict axillary LN metastasis (17). Wen et al. reported that the Emax and Emean of primary tumors were higher in LN metastasis cases (24). In our study, however, quantitative SWE parameters including Emean, Emax, Esd, and Eratio did not show a significant difference between patients with positive and negative LN metastasis groups (p > 0.05 for all), indicating that these SWE parameters of intrinsic tumor and peritumor tissues might not provide sufficient information in predicting LN metastasis for our research group. The difference between the results of this study and previous studies may be due to the difference of the study sample and the difference of the ultrasound instruments used.

In our study, SWE parameters including the Emin of internal tumor and shell 3 mm and the “stiff rim” sign showed a significant difference between patients with positive and negative LN metastasis groups (p < 0.05 for all). Zhou et al. have reported that the “stiff rim” sign could select the malignant breast lesions, and possible reasons may contribute to cancer cell infiltration into interstitial tissues or a desmoplastic reaction (16), indicating that a qualitative assessment of peritumoral tissue stiffness may reflect the progress of invasive malignant tumors to some extent. Peritumoral stroma and peritumoral invasion played a critical role in the spreading and metastasis of breast cancer (18, 19, 25, 26) and were also responsible for the formation of the “stiff rim” sign, as well as the attenuation of shear wave energy (16). The low-quality shear waves caused by attenuation may be interpreted as slow-speed shear waves (16). Therefore, the “stiff rim” sign was mostly found in the positive LN metastasis group, and the Emin of the positive LN metastasis group was lower than that of the negative LN metastatic group.

There was a significant difference in the constituent ratio between 5 BI-RADS categories (p < 0.001). In the 619 masses, 481 lesions were defined as BI-RADS 3-4C, with about 25% (ranging from 22.4% to 27.3%) cases in positive metastasis LN, while in 138 masses of BI-RADS 5, patients with positive LN metastasis were up to 43.5%. Guo et al. revealed that irregular shape and higher color Doppler flow imaging grades were associated with axillary LN metastasis of breast cancer (27). Bae et al. showed that architectural distortion of breast tumors was associated with LN metastasis of breast cancer (28). Previous studies also indicated that a non-circumscribed margin and heterogeneous internal echo were potential predictors of axillary LN metastasis independently (28, 29). Possible reasons might be that the non-circumscribed margins reflect the invasiveness of tumors and aggressive of metastatic lesions to lymphocytes. Besides, tumor growth rate might influence the internal echo of the masses, and rapid tumor growth could lead to an increased likelihood of LN metastasis (29). These ultrasound features, such as irregular shape, non-circumscribed margin, and architectural distortion, were considered as determining characteristics of ultrasound BI-RADS, which may explain the higher incidence of LN metastasis in BI-RADS 5 tumors in this study.

A previous study revealed that tumor staging was an independent predictor of breast cancer LN metastasis (30). A tumor size of 20 mm or greater on ultrasound was an independent predictor of LN metastasis (31). In the present study, however, there was no significant difference in the number of tumors evaluated as T1 and T2/T3 in LN-positive and -negative groups (p = 0.340), which was similar to the results of the previous study; that is, tumor size on ultrasound was not an independent predictor of axillary LN metastasis (28). Breast tumor location has been acknowledged as an independent factor of tumor prognosis, and malignant tumors in the lower inner quadrant were related to a worse prognosis because of the internal mammary lymphatic pathway without sufficient studies (32). However, there was no significant difference in tumor location (p = 0.390) in LN-positive and -negative groups in our study.

The grade system has been extensively used in categorical variables of clinical and imaging data analysis (1, 33). In the prediction of LN metastasis, there was no statistical difference among the logistic regression equation, the weighting method, and the counting method (p > 0.05 for all). For ease of use, a 0–4 grade risk stratification system for LN metastasis was established based on the counting method in this study. The risk system was based on quantitative and qualitative SWE features as well as BI-RADS category. To the best of our knowledge, no studies have combined quantitative and qualitative SWE characteristics of tumors, especially peritumor tissues, and the BI-RADS category to predict axillary LN metastasis in breast invasive cancer. Our results showed that the risk of axillary LN metastasis of breast cancer increased as the risk grade rises from 0 to 4, with the increase of the four significant factors (BI-RADS, “stiff rim” sign, Emin of internal tumor, and shell 3 mm). A patient with a BI-RADS 5 breast tumor exhibiting the “stiff rim” sign (value of Emin in internal tumor and shell 3 mm ≤3.31 and ≤3.52, respectively) was assessed as grade 4 in our risk system. The OR of each grade in the risk system gradually increased as the increasing levels were 3.095 (p = 0.074, 95% CI = 0.896–10.696) for grade 1, 4.211 (p < 0.05, 95% CI = 1.243–14.271) for grade 2, 7.172 (p < 0.05, 95% CI = 2.132–24.132) for grade 3, and 18.318 (p < 0.05, 95% CI = 5.016–66.892) for grade 4, respectively, compared with grade 0, indicating that patients with a higher risk grade were more likely to have LN metastasis. The diagnostic performance of this risk grade system is significant in assessing the risk grade of LN metastasis, although the AUC is not excellent (0.656, 95% CI, 0.617–0.693).

The merit of our research is that this is a multicenter study covering 17 hospitals, which effectively improves the reliability of the research data. However, the study still has some limitations. Firstly, the evaluation of the BI-RADS category and the SWE features of breast lesions were influenced by the experience of the radiologist. However, previous studies have confirmed that the reproducibility ranged from moderate to substantial for the BI-RADS category and ranged from substantial to almost perfect for both the quantitative SWE features and the “stiff rim” sign (16, 34, 35). Secondly, all SWE information were acquired from the same ultrasound system (Resona 7, Mindray), which would limit the generalizability of the study results. Further studies using different ultrasound SWE systems are warranted in the future. Third, the AUC of the risk system is not very high. In the future, better predictive models may be obtained if the information provided by ultrasound, MRI, mammography, and clinical data is utilized jointly.

Conclusion

A 0–4 grade risk stratification system, based on quantitative and qualitative SWE features as well as BI-RADS category, for predicting axillary LN metastasis in invasive breast cancer was established in this study. Although the diagnostic performance of this risk system was not excellent, there is no doubt that it has the potential to predict axillary LN metastasis by the combination of SWE and conventional ultrasound.

Data Availability Statement

The datasets presented in this article are not readily available because the generated datasets belong to department of ultrasound, Ruijin Hospital, Shanghai Jiaotong University School of Medicine. Requests to access the datasets should be directed to JQZ, emhvdXN1MzBAMTI2LmNvbQ==.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of Ruijin Hospital, Shanghai Jiaotong University School of Medicine. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

JWZ drafted and edited the manuscript. HTZ, YJD, XHJ, JWZ, ZYL, ZRC, YJX, BH, YXH, CC, JFX, FJD, XNX, CRW, WJH, GW, QYL, QC, WYD, QCJ, YLM, HNY, XJX, HJY, PZ, YS, LGC, PH, LXQ, JPL, LYS, YNZ, and YYX collected the clinical and ultrasound data. YJD and WWZ analyzed the images. XHJ collected the ultrasound images. XHJ and JWZ collected the autopsy and pathology data. HTZ prepared for the literature. YYS and HTZ completed the statistical analysis. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the National Natural Science Foundation of China (grant no. 82071928).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We would like to express our sincere gratitude to all doctors from the 17 hospitals participating in the multicenter study who provided data for this study.

References

1. T. I. A. f. R. o. Cancer. Latest Global Cancer Data: Cancer Burden Rises to 19.3 Million New Cases and 10.0 Million Cancer Deaths in 2020. Geneva: World Health Organization (2020).

Google Scholar

2. Araujo T, Aresta G, Castro E, Rouco J, Aguiar P, Eloy C, et al. Classification of Breast Cancer Histology Images Using Convolutional Neural Networks. PloS One (2017) 12(6):e0177544. doi: 10.1371/journal.pone.0177544

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Cao Z, Duan L, Yang G, Yue T, Chen Q. An Experimental Study on Breast Lesion Detection and Classification From Ultrasound Images Using Deep Learning Architectures. BMC Med Imaging (2019) 19(1):51. doi: 10.1186/s12880-019-0349-x

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Zeng H, Zheng R, Zhang S, Zou X, Chen W. Female Breast Cancer Statistics of 2010 in China: Estimates Based on Data From 145 Population-Based Cancer Registries. J Thorac Dis (2014) 6(5):466–70. doi: 10.3978/j.issn.2072-1439.2014.03.03

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Tao C, Chen K, Han L, Peng Y, Li C, Hua Z, et al. New One-Step Model of Breast Tumor Locating Based on Deep Learning. J Xray Sci Technol (2019) 27(5):839–56. doi: 10.3233/XST-190548

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Giuliano AE, Edge SB, Hortobagyi GN. Eighth Edition of the AJCC Cancer Staging Manual: Breast Cancer. Ann Surg Oncol (2018) 25(7):1783–5. doi: 10.1245/s10434-018-6486-6

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Lyman GH, Giuliano AE, Somerfield MR, Benson AB, Bodurka 3DC, Burstein HJ, et al. American Society of Clinical Oncology Guideline Recommendations for Sentinel Lymph Node Biopsy in Early-Stage Breast Cancer. J Clin Oncol (2005) 23(30):7703–20. doi: 10.1200/JCO.2005.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Mercado CL. BI-RADS Update. Radiologic Clinics of North America (2014) 52(3):481-7. doi: 10.1016/j.rcl.2014.02.008

CrossRef Full Text | Google Scholar

9. Ehteshami Bejnordi B, Veta M, Johannes van Diest P, van Ginneken B, Karssemeijer N, Litjens G, et al. Diagnostic Assessment of Deep Learning Algorithms for Detection of Lymph Node Metastases in Women With Breast Cancer. JAMA (2017) 318(22):2199–210. doi: 10.1001/jama.2017.14585

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Zheng X, Yao Z, Huang Y, Yu Y, Wang Y, Liu Y, et al. Deep Learning Radiomics can Predict Axillary Lymph Node Status in Early-Stage Breast Cancer. Nat Commun (2020) 11(1):1236. doi: 10.1038/s41467-020-15027-z

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Valente SA, Levine GM, Silverstein MJ, Rayhanabad JA, Weng-Grumley JG, Ji L, et al. Accuracy of Predicting Axillary Lymph Node Positivity by Physical Examination, Mammography, Ultrasonography, and Magnetic Resonance Imaging. Ann Surg Oncol (2012) 19(6):1825–30. doi: 10.1245/s10434-011-2200-7

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Chou Y-H, Tiu C-M, Lai Y-C, Lin Y-H. Elastography of Breast Tumors: Pitfalls. Ultrasound Med Biol (2019) 45:S22. doi: 10.1016/j.ultrasmedbio.2019.07.484

CrossRef Full Text | Google Scholar

13. Nori J, Vanzi E, Bazzocchi M, Bufalini FN, Distante V, Branconi F, et al. Role of Axillary Ultrasound Examination in the Selection of Breast Cancer Patients for Sentinel Node Biopsy. Am J Surg (2007) 193(1):16–20. doi: 10.1016/j.amjsurg.2006.02.021

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Lee MC, Eatrides J, Chau A, Han G, Kiluk JV, Khakpour N, et al. Consequences of Axillary Ultrasound in Patients With T2 or Greater Invasive Breast Cancers. Ann Surg Oncol (2011) 18(1):72–7. doi: 10.1245/s10434-010-1171-4

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Mercado CL. BI-RADS Update. Radiol Clin North Am (2014) 52(3):481–7. doi: 10.1016/j.rcl.2014.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Zhou J, Zhan W, Chang C, Zhang X, Jia Y, Dong Y, et al. Breast Lesions: Evaluation With Shear Wave Elastography, With Special Emphasis on the "Stiff Rim" Sign. Radiology (2014) 272(1):63–72. doi: 10.1148/radiol.14130818

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Evans A, Rauchhaus P, Whelehan P, Thomson K, Purdie CA, Jordan LB, et al. Does Shear Wave Ultrasound Independently Predict Axillary Lymph Node Metastasis in Women With Invasive Breast Cancer? Breast Cancer Res Treat (2014) 143(1):153–7. doi: 10.1007/s10549-013-2747-z

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Colleoni M, Rotmensz N, Maisonneuve P, Sonzogni A, Pruneri G, Casadio C, et al. Prognostic Role of the Extent of Peritumoral Vascular Invasion in Operable Breast Cancer. Ann Oncol (2007) 18(10):1632–40. doi: 10.1093/annonc/mdm268

PubMed Abstract | CrossRef Full Text | Google Scholar

19. de Mascarel I, Bonichon F, Durand M, Mauriac L, MacGrogan G, Soubeyran I, et al. Obvious Peritumoral Emboli: An Elusive Prognostic Factor Reappraised. Multivariate Analysis of 1320 Node-Negative Breast Cancers. Eur J Cancer (1998) 34(1):58–65. doi: 10.1016/s0959-8049(97)00344-4

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Barr RG, Nakashima K, Amy D, Cosgrove D, Farrokh A, Schafer F, et al. WFUMB Guidelines and Recommendations for Clinical Use of Ultrasound Elastography: Part 2: Breast. Ultrasound Med Biol (2015) 41(5):1148–60. doi: 10.1016/j.ultrasmedbio.2015.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Zhou J, Song Y, Zhan W, Wei X, Zhang S, Zhang R, et al. Thyroid Imaging Reporting and Data System (TIRADS) for Ultrasound Features of Nodules: Multicentric Retrospective Study in China. Endocrine (2021) 72(1):157–70. doi: 10.1007/s12020-020-02442-x

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Evans A, Whelehan P, Thomson K, McLean D, Brauer K, Purdie C, et al. Invasive Breast Cancer: Relationship Between Shear-Wave Elastographic Findings and Histologic Prognostic Factors. Radiology (2012) 263(3):673–7. doi: 10.1148/radiol.12111317

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Youk JH, Gweon HM, Son EJ, Kim JA, Jeong J. Shear-Wave Elastography of Invasive Breast Cancer: Correlation Between Quantitative Mean Elasticity Value and Immunohistochemical Profile. Breast Cancer Res Treat (2013) 138(1):119–26. doi: 10.1007/s10549-013-2407-3

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Wen X, Yu X, Tian Y, Liu Z, Cheng W, Li H, et al. Quantitative Shear Wave Elastography in Primary Invasive Breast Cancers, Based on Collagen-S100A4 Pathology, Indicates Axillary Lymph Node Metastasis. Quant Imaging Med Surg (2020) 10(3):624–33. doi: 10.21037/qims.2020.02.18

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Wernicke M, Roitman P, Manfre D, Stern R. Breast Cancer and the Stromal Factor. The "Prometastatic Healing Process" Hypothesis. Med (B Aires) (2011) 71(1):15–21.

Google Scholar

26. Auvinen P, Tammi R, Parkkinen J, Tammi M, Agren U, Johansson R, et al. Hyaluronan in Peritumoral Stroma and Malignant Cells Associates With Breast Cancer Spreading and Predicts Survival. Am J Pathol (2000) 156(2):529–36. doi: 10.1016/S0002-9440(10)64757-8

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Guo Q, Dong Z, Zhang L, Ning C, Li Z, Wang D, et al. Ultrasound Features of Breast Cancer for Predicting Axillary Lymph Node Metastasis. J Ultrasound Med (2018) 37(6):1354–3. doi: 10.1002/jum.14469

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Bae MS, Shin SU, Song SE, Ryu HS, Han W, Moon WK. Association Between US Features of Primary Tumor and Axillary Lymph Node Metastasis in Patients With Clinical T1-T2N0 Breast Cancer. Acta Radiol (2018) 59(4):402–8. doi: 10.1177/0284185117723039

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Zhang H, Sui X, Zhou S, Hu L, Huang X. Correlation of Conventional Ultrasound Characteristics of Breast Tumors With Axillary Lymph Node Metastasis and Ki-67 Expression in Patients With Breast Cancer. J Ultrasound Med (2019) 38(7):1833–40. doi: 10.1002/jum.14879

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Stachs A, Thi AT, Dieterich M, Stubert J, Hartmann S, Glass A, et al. Assessment of Ultrasound Features Predicting Axillary Nodal Metastasis in Breast Cancer: The Impact of Cortical Thickness. Ultrasound Int Open (2015) 1(1):E19–24. doi: 10.1055/s-0035-1555872

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Zhu AQ, Li XL, An LW, Guo LH, Fu HJ, Sun LP, et al. Predicting Axillary Lymph Node Metastasis in Patients With Breast Invasive Ductal Carcinoma With Negative Axillary Ultrasound Results Using Conventional Ultrasound and Contrast-Enhanced Ultrasound. J Ultrasound Med (2020) 39(10):2059–70. doi: 10.1002/jum.15314

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Gou ZC, Liu XY, Xiao Y, Zhao S, Jiang YZ, Shao ZM. Decreased Survival in Patients With Carcinoma of Axillary Tail Versus Upper Outer Quadrant Breast Cancers: A SEER Population-Based Study. Cancer Manag Res (2018) 10:1133–41. doi: 10.2147/CMAR.S165291

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Liu Z, Wang R, Zhou J, Zheng Y, Dong Y, Luo T, et al. Ultrasound Lymphatic Imaging for the Diagnosis of Metastatic Central Lymph Nodes in Papillary Thyroid Cancer. Eur Radiol (2021) 31(11):8458–67. doi: 10.1007/s00330-021-07958-y

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Schwab F, Redling K, Siebert M, Schotzau A, Schoenenberger CA, Zanetti-Dallenbach R. Inter- and Intra-Observer Agreement in Ultrasound BI-RADS Classification and Real-Time Elastography Tsukuba Score Assessment of Breast Lesions. Ultrasound Med Biol (2016) 42(11):2622–9. doi: 10.1016/j.ultrasmedbio.2016.06.017

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Cosgrove DO, Berg WA, Dore CJ, Skyba DM, Henry JP, Gay J, et al. Shear Wave Elastography for Breast Masses is Highly Reproducible. Eur Radiol (2012) 22(5):1023–32. doi: 10.1007/s00330-011-2340-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: breast neoplasms, lymphatic metastasis, ultrasonography, elasticity imaging techniques, risk assessment

Citation: Zhang H, Dong Y, Jia X, Zhang J, Li Z, Chuan Z, Xu Y, Hu B, Huang Y, Chang C, Xu J, Dong F, Xia X, Wu C, Hu W, Wu G, Li Q, Chen Q, Deng W, Jiang Q, Mou Y, Yan H, Xu X, Yan H, Zhou P, Shao Y, Cui L, He P, Qian L, Liu J, Shi L, Zhao Y, Xu Y, Song Y, Zhan W and Zhou J (2022) Comprehensive Risk System Based on Shear Wave Elastography and BI-RADS Categories in Assessing Axillary Lymph Node Metastasis of Invasive Breast Cancer—A Multicenter Study. Front. Oncol. 12:830910. doi: 10.3389/fonc.2022.830910

Received: 07 December 2021; Accepted: 14 February 2022;
Published: 10 March 2022.

Edited by:

Kyung Hyun Sung, UCLA Health System, United States

Reviewed by:

Kaustav Bera, Maimonides Medical Center, United States
Wenwu Ling, Sichuan University, China

Copyright © 2022 Zhang, Dong, Jia, Zhang, Li, Chuan, Xu, Hu, Huang, Chang, Xu, Dong, Xia, Wu, Hu, Wu, Li, Chen, Deng, Jiang, Mou, Yan, Xu, Yan, Zhou, Shao, Cui, He, Qian, Liu, Shi, Zhao, Xu, Song, Zhan and Zhou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jianqiao Zhou, emhvdXN1MzBAMTI2LmNvbQ==; Yanyan Song, c3l5c3NtdUAxMjYuY29t; Weiwei Zhan, c2hhbmdoYWlydWlqaW51c0AxNjMuY29t

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.