Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 17 October 2023
Sec. Cancer Imaging and Image-directed Interventions
This article is part of the Research Topic Breast Cancer Imaging: Clinical Translation of Novel Methods View all 6 articles

A clinical-radiomics nomogram based on multimodal ultrasound for predicting the malignancy risk in solid hypoechoic breast lesions

Guo ShiyanGuo ShiyanJiang LiqingJiang LiqingYan YueqiongYan YueqiongZhang Yan*Zhang Yan*
  • Department of Ultrasound, Third Xiangya Hospital, Central South University, Changsha, Hunan, China

Background: In routine clinical examinations, solid hypoechoic breast lesions are frequently encountered, but accurately distinguishing them poses a challenge. This study proposed a clinical-radiomics nomogram based on multimodal ultrasound that enhances the diagnostic accuracy for solid hypoechoic breast lesions.

Method: This retrospective study analyzed ultrasound strain elastography (SE) and automated breast volume scanner images (ABVS) of 423 solid hypoechoic breast lesions from 423 female patients in our hospital between August 2019 and May 2022. They were assigned to the training (n=296) and validation (n=127) groups in a 7:3 ratio by generating random numbers. Radiomics features were extracted and screened from ABVS and SE images, followed by the calculation of the radiomics score (Radscore) based on these features. Subsequently, a nomogram was constructed through multivariate logistic regression to assess the malignancy risk in breast lesions by combining Radscore with Breast Imaging Reporting and Data System (BI-RADS) scores and clinical risk factors associated with breast malignant lesions. The diagnostic performance, calibration performance, and clinical usefulness of the nomogram were assessed by the area under the curve (AUC) of the receiver operating characteristic curve, the calibration curve, and the decision analysis curve, respectively.

Results: The diagnostic performance of the nomogram is significantly superior to that of both the clinical diagnostic model (BI-RADS model) and the multimodal radiomics model (SE+ABVS radiomics model) in training (AUC: 0.972 vs 0.930 vs 0.941) and validation group (AUC:0.964 vs 0.916 vs 0.933). In addition, the nomogram also exhibited a favorable goodness-of-fit and could lead to greater net benefits for patients.

Conclusion: The nomogram enables a more effective assessment of the malignancy risk of solid hypoechoic breast lesions; therefore, it can serve as a new and efficient diagnostic tool for clinical diagnosis.

1 Introduction

As the most prevalent cancer in the world, breast cancer poses a grave threat to people’s health and survival (1). Given its high metastatic tendency and high mortality rate (2, 3), coupled with the significant differences in treatment modalities for benign and malignant breast tumors, early definitive diagnosis is a critical first step in the therapeutic management of breast lesions, which plays a crucial role in improving patient outcomes and survival (35).

With the recent advancements in ultrasound imaging technology, ultrasound plays an increasingly important role in the detection of breast lesions. Strain elastography (SE) allows for a quick and intuitive display of differences in elasticity coefficients within the lesion through color-coded imaging, therefore, it serves as a powerful diagnostic aid to offer valuable reference values for lesion diagnosis (6, 7). Automated breast volume scanner (ABVS) provides good reproducibility of diagnostic results due to its standardized operating procedures (8), it can acquire the whole breast volume information and perform multiplanar imaging on the acquired information. Studies have shown that ABVS exhibits comparable diagnostic accuracy to handheld ultrasound scanners for detecting breast lesions, while also providing additional information (9, 10). In routine ultrasound examinations, it is frequent to encounter patients with solid hypoechoic breast lesions, physicians can make an initial assessment of the malignancy risk of breast lesions based on their morphological appearance on ABVS image and elastic performance on SE images. The combination of ABVS and SE imaging techniques demonstrates significant diagnostic efficacy in evaluating breast lesions (11). However, the dependability of diagnostic outcomes generated by conventional imaging techniques is largely contingent on the proficiency of the examining physician and is markedly susceptible to interobserver variability (12).

Radiomics is in line with the current trend toward precision medicine, as it transforms ordinary visual images into high-throughput data through deep mining of medical images, allowing for the capturing of the internal heterogeneity of the entire tumor in a non-invasive manner (1315). Therefore, it may provide novel biomarkers to facilitate diagnosis for better clinical decision-making. There are already several radiomics studies on ultrasound (US), mammography, and magnetic resonance (MR) in breast cancer diagnosis that have yielded promising results (1629). However, there have been no studies on the combination of ABVS and UE radiomics features with clinical ultrasound factors for the diagnosis of breast cancer. Therefore, we conducted a radiomics analysis on SE and ABVS images, then combined these features with traditional imaging risk assessments and other clinical risk factors, resulting in a novel nomogram to help physicians accurately diagnose solid hypoechoic breast lesions.

2 Materials and methods

2.1 Patients

The retrospective study was approved by the institutional review board at our hospital. The inclusion criteria were: (1) Patients who underwent both ABVS and SE examinations at our hospital between August 2019 and May 2022 and subsequently underwent biopsy or surgical resection within two weeks with a pathologically confirmed diagnosis. (2) Patients with complete imaging data. (3) Patients’ breast lesions were hypoechoic solid lesions. The exclusion criteria were: (1) Patients whose images were of poor quality; (2) Patients who underwent aspiration or clinical treatment before examining target lesions. Eventually, we included 423 solid hypoechoic breast lesions from 423 female patients. By generating random numbers, they were allocated into training and validation groups in the ratio of 7:3. The flow is shown in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1 The grouping process of this study.

2.2 Image acquisition and assessment

In this study, all images used were obtained using the ACUSON S2000 US machine and its accompanying ABVS system.

2.2.1 SE image

Patients were instructed to breathe normally while lying supine on the examination bed with their breasts fully exposed. A 9L4 probe in two-dimensional ultrasound mode was used to examine the breast in all planes. The imaging mode was then switched to the elastic mode when scanning the largest two-dimensional section of the lesions, and the patient was required to cooperate by holding her breath. The probe is placed perpendicularly over the breast without applying any pressure during the capture of SE images, with the lesion positioned at the center of an elasticity sampling window at least twice the size of the area of interest.

2.2.2 ABVS image

Instruct the patient to raise both arms over the head and remain in the supine position. A sufficient amount of coupling agent was applied uniformly to the breast. Before scanning, parameters such as depth and overall gain were adjusted to achieve optimal image quality. During the scanning procedure, patients were instructed to breathe normally. Each breast was routinely scanned in two positions using a 14L5BV high-frequency linear array automatic scanning US probe. The nipples were marked after the scanning, and the acquired images were saved and transferred to a workstation for processing and analysis. Images with the maximum section image of the target lesion in coronal, transverse, and sagittal planes were selected for subsequent region of interest (ROI) segmentation and feature extraction.

Refer to the BI-RADS criteria which were defined by the American College of Radiology in 2013 (30), we evaluated the morphology, margin, border, orientation, posterior echogenicity of the lesion, microcalcification within the lesion, and conditions of retraction in its coronal plane on the saved images. This is followed by a combination with ultrasound elastic strain performance of the lesions (elasticity score using a 5-point scale (31)), enabling an accurate classification of the risk of malignancy of the lesions.

All of the above steps were performed by a physician who has over ten years of expertise in ultrasound breast disease diagnosis.

2.3 Extraction and selection of radiomics features

The ABVS and SE images were sequentially imported into 3D Slicer 5.2.1 for image processing, manual segmentation of the ROI and extraction of radiomics features. We delineated ROI for lesions in both SE and ABVS images. Notably, in the ABVS images, we delineated the ROI on the coronal, sagittal, and transverse planes of the lesions, respectively. Further details can be found in Figure 2. This procedure was performed with the participation of two physicians. Physician A, with five years of experience in ultrasound-based breast disease diagnosis, performed outlining for all the lesions. Physician B, with eight years of experience in ultrasound-based breast disease diagnosis, conducted lesion outlining on the training group to validate ROI outlining reproducibility. Then, we utilized the Pyradiomics package within 3D Slicer to extract radiomics features from the SE and ABVS images, respectively. The features extracted encompassed first-order statistics features, texture features (including the gray level co-occurrence matrix (glcm), gray level dependence matrix (gldm), gray level run length matrix (glrlm), gray level size zone matrix (glszm), and the neighbouring gray tone difference matrix (ngtdm)), as well as post-wavelet transformed features. Subsequently, we subjected the extracted features to screening. By utilizing the intra-class correlation coefficient (ICC) analysis, we can identify features that exhibit high levels of reproducibility (ICC > 0.75) (32), subsequently, radiomics features extracted from the region of interest segmented by physician A were utilized for further analysis. All feature values were normalized using Zscore. The radiomics features of both modalities were subjected to dimensionality reduction through the Mann-Whitney U test and least absolute shrinkage and selection operator (LASSO) regression and to identify features with strong qualitative diagnostic ability for solid hypoechoic breast lesions.

FIGURE 2
www.frontiersin.org

Figure 2 An instance of manually delineating a region of interest (ROI). The strain elastography (SE) and automated breast volume scanner (ABVS) images of a 41-year-old female with a solid hypoechoic lesion measuring approximately 16x11x12mm on her left breast. The lesion was irregular in shape, parallel in position, with still well-defined borders, sharp margins, and scattered microcalcifications visible internally, and exhibited no significant posterior echogenicity change or retraction in the coronal plane, and the ultrasound elasticity score was 4, finally, the lesion was classified as BI-RADS category 4a. Pathological examination confirmed it as invasive ductal carcinoma. ROI segmentation was performed on both the SE image (A) and ABVS coronal image (B), with delineation along the boundary of the lesion followed by uniform outward expansion of its edges by 3 mm to encompass some surrounding tissue.ROI segmentation was performed on ABVS transverse (C) and sagittal (D) images, respectively, and meticulous delineation was performed along the lesion’s contour and borders on these two views.

2.4 Development of models

2.4.1 Radiomics models

Based on radiomics features of SE and ABVS images, logistic regression analysis was utilized to build the radiomics model, including the SE radiomics model, the ABVS radiomics model, the SE+ABVS radiomics model that was developed by combining the two. The radiomics score (Radscore) for each lesion was computed by weighting the coefficients of features in the SE+ABVS radiomics model.

2.4.2 Clinical diagnostic model

The BI-RADS model was constructed through logistic regression analysis of lesion’s BI-RADS categories in the training group.

2.5 Development and performance validation of nomogram

This study performed a univariate analysis in order to determine the risk predictive variables associated with breast cancer(P<0.05), which were then combined with the results of conventional imaging assessment and radiomics analysis. Based on these findings, we integrated relevant clinical risk factors, the Radscore, and the BI-RADS category of lesions to develop a nomogram for assessing the malignancy risk in such breast lesions by multivariate logistic regression analysis. Subsequently, the nomogram’s diagnostic performance was compared to that of the BI-RADS model and SE+ABVS radiomics model. To evaluate the diagnostic performance of the models, we calculated the area under the receiver operating characteristic curve (AUC) for each model in the training group, validation group, and in the BI-RADS category 4 lesions within both groups. Furthermore, the DeLong test was used to examine differences in AUC values between different models. The nomogram’s goodness of fit was investigated graphically and by calculating significance by plotting the calibration curve and conducting the Hosmer-Lemeshow test. Lastly, clinical decision analysis curves were drawn for quantifying the net benefits of the BI-RADS model, SE+ABVS radiomics models, and nomogram at various threshold probabilities.

2.6 Statistical analysis

SPSS 23.0, R 4.2.2, and MedCalc 19.6.0 were utilized for statistical analysis and graph plotting. The ‘psych’, ‘survival’, ‘glmnet’, ‘rms’, ‘ResourceSelection’, and ‘rmda’ packages were used in R. We performed normality tests on each group of data and selected the appropriate hypothesis test based on the results to compare the distribution of data between the training and validation groups. The study has chosen a significance level of 0.05 as the threshold for detecting statistical differences.

3 Results

3.1 Comparison of clinical basis information and sonographic features

The study included 423 breast lesions that were pathologically confirmed to include 215 benign lesions and 208 malignant lesions. Table 1 demonstrates that both the clinical basis data and sonographic features of lesions were evenly distributed in the training and validation groups, indicating no statistically significant differences between the two groups (P>0.05). Furthermore, the univariate risk analysis revealed that patients with malignant lesions had significantly higher age and lesion’s maximum diameter compared to those with benign lesions in both groups (P<0.05). However, no correlation was seen between the location of the lesion and the malignant risk of the lesion (P>0.05). Hence, we regarded age and lesion size as predictor variables in the context of breast cancer. Regarding the sonographic features of the lesions, there were statistically significant differences (P<0.01) observed in morphology, borders, margins, orientation, microcalcifications, retraction condition of the coronal plane, and elasticity scores between benign and malignant lesions within both groups. while no statistical differences were found in posterior echogenicity (P>0.05). This study assessed the malignancy risk of lesions by these sonographic features of them, and the BI-RADS categories obtained were also significantly different in benign and malignant lesions (P<0.01).

TABLE 1
www.frontiersin.org

Table 1 Clinical basis information and sonographic features of patients with breast lesions.

3.2.Screening of radiomics features

The SE and ABVS images generated 837 and 2511 radiomics features (ABVS cross plane, sagittal plane, and coronal plane, each generated 837 features), respectively. The training group’s SE radiomics features as well as the ABVS coronal plane, transverse plane, and sagittal plane radiomics features underwent sequential ICC analysis, Mann-Whitney U test, LASSO regression analysis with tenfold cross-validation for dimensionality reduction. Finally, a total of 14 features were selected, comprising four SE radiomics features and ten ABVS radiomics features (two from the coronal plane, three from the transverse plane, and five from the sagittal plane). All of these radiomics features are texture features, one of which was from the original image and thirteen were obtained after wavelet transform (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3 Screening of radiomics features. Selection of strain elastography (SE) radiomics features (A), automated breast volume scanner (ABVS) coronal plane radiomics features (B), ABVS transverse plane radiomics features (C), and ABVS sagittal plane radiomics features (D) using the least absolute shrinkage and selection operator (LASSO) regression model. The coefficient profiles of LASSO for each modal radiomics feature are presented on the left. The right shows that the tuning parameter λ (lambda) in the LASSO model was selected using tenfold cross-validation, and the binomial deviance was plotted as a function of log(λ), with vertical dashed lines drawn at the minimum deviation (log(λ.min)) and the 1 standard error of the minimum deviation (log(λ.1se)). Selected the non-zero coefficient features in the model when the horizontal coordinate was log(λ.1se).

3.3 Comparison of radiomics models

By comparing and validating the diagnostic efficacy of the radiomics models (Figure 4, Table 2), the AUC values of the selected ABVS and SE features for distinguishing between benign and malignant solid hypoechoic breast lesions were consistently above 0.8 in both the training and validation groups. Moreover, compared to any single-modality radiomics models, the SE+ABVS radiomics model, which integrated the radiomics features of two imaging modalities, demonstrated significantly higher AUC values in both training (All P<0.01) and validation groups (compared to the ABVS radiomics model: P<0.01, compared to SE radiomics model: P<0.05). These outcomes suggest that combining radiomics features from both SE and ABVS could enhance the accuracy of diagnostic models. Thus, the Radscore for each patient was obtained by weighting the corresponding coefficients for each feature in the SE+ABVS radiomics model., the formula is shown below, the Radscore for malignant lesions was found to be significantly higher than that for benign lesions within both groups. (Training group: 2.86 + 2.66, -2.34 + 1.80, P<0.01; Validation group: 2.60 + 2.31, -2.17 + 1.85, P<0.01).

FIGURE 4
www.frontiersin.org

Figure 4 The receiver operator characteristic curves for various radiomics models in the training (A) and validation groups (B).

TABLE 2
www.frontiersin.org

Table 2 The AUC values of radiomics models in the training and validation groups.

Radscore=1.944469*original glszm GrayLevelVariance_SE+0.03492*waveletLHL glszm GrayLevelNonUniformity_SE+0.634491*waveletLHH gldm LargeDependenceLowGrayLevelEmphasis_SE+0.465968*waveletLLL glcm SumSquares_SE+0.133391*waveletHHL glszm GrayLevelNonUniformity_ABVS_Coronal1.515639*waveletLLL glszm ZoneEntropy_ABVS_Coronal+0.010666*waveletLHH glszm ZoneEntropy_ ABVS_Transverse+1.292836*waveletHLL glszm ZoneEntropy_ ABVS_Transverse+0.043941*waveletHHL glszm GrayLevelNonUniformity_ ABVS_Transverse1.828919*waveletLHH glcm MCC_ABVS_Sagittal+0.813141*waveletLHH glszm ZoneEntropy_ABVS_Sagittal0.008279*waveletHLL glszm GrayLevelNonUniformity_ABVS_Sagittal+1.520859*waveletHLH glszm ZoneEntropy_ABVS_Sagittal+0.271946*waveletHHL glszm GrayLevelNonUniformity_ABVS_Sagittal25.66532

3.4 Evaluation of nomogram performance

Based on the clinical risk factors identified through univariate analysis, BI-RADS categories determined from imaging assessments, and Radscore obtained from radiomics analysis, we constructed a nomogram using multivariate logistic regression to visually assess the risk of malignancy in solid hypoechoic breast lesions, the nomogram incorporated the patient’s age, lesion’s maximum diameter, Radscore, and BI-RADS category. As illustrated in Figure 5, Radscore had the highest weightage followed by BI-RADS score while age and maximum diameter of the lesion exerted less influence on assessment results.

FIGURE 5
www.frontiersin.org

Figure 5 The Nomogram for predicting the malignant risk of solid hypoechoic breast lesions.

Figure 6 and Table 3 present that the BI-RADS model, SE+ABVS radiomics model, and nomogram are effective in predicting the malignancy risk in solid hypoechoic breast lesions, Notably, the nomogram exhibits superior diagnostic performance with higher AUC values (0.972, 0.964) in training and validation group compared to both the BI-RADS model (AUC: 0.930, 0.916) and SE+ABVS radiomics models (AUC: 0.941, 0.933). Furthermore, its difference with BI-RADS model and SE+ABVS radiomics model was statistically significant in both groups (P<0.05). Besides, we further compared the diagnostic efficacy of the three models for BI-RADS category 4 lesions within the two groups. The results revealed that the nomogram (AUC: 0.952, 0.930) consistently exhibited higher AUC values than both the BI-RADS model (AUC:0.844, 0.839) and SE+ABVS radiomics model (AUC:0.915, 0.899). Moreover, there were consistently statistically significant differences between the nomogram and BI-RADS model (All P<0.01). However, in comparison to the SE+ABVS model, the nomogram was only statistically different from it in the training group (P<0.05), but not in the validation group (P>0.05). Other than that, in terms of diagnostic sensitivity, specificity, and accuracy, Although the specificity of the nomogram was slightly inferior to that of the BI-RADS model in the training group, it significantly improved diagnostic sensitivity. Furthermore, its diagnostic parameters were at the highest level across all validation groups. These results suggest that the nomogram exhibited the best overall diagnostic performance. Finally, we observed that the AUC values of the SE+ABVS Radiomics model consistently outperformed those of the BI-RADS model, and a statistically significant difference was found between them when diagnosing BI-RADS category 4 lesions of the training group (P<0.05). This finding highlights the ability of radiomics analysis to detect deep-seated features within the images, ultimately leading to improved diagnostic efficiency.

FIGURE 6
www.frontiersin.org

Figure 6 The receiver operator characteristic curves of the BI-RADS model, SE+ABVS radiomics model, and Nomogram in the training group (A), the validation group (B), the BI-RADS category 4 lesions in the training group (C), and the BI-RADS category 4 lesions in the validation group (D).

TABLE 3
www.frontiersin.org

Table 3 The diagnostic parameters of the BI-RADS model, SE+ABVS radiomics model, and Nomogram in each group.

The calibration curve exhibits a favorable fit of the nomogram (Figures 7A, B). indicating that the predicted risk by the nomogram was close to the observed risks. The results from the Hosmer-Lemeshow test further proved that the differences between them did not present statistical significance in either the training group (P=0.70) or validation group (P=0.95).

FIGURE 7
www.frontiersin.org

Figure 7 The calibration curves for the Nomogram in the training group (A) and the validation group (B).The decision analysis curves of the BI-RADS mode, SE+ABVS radiomics model, and Nomogram in the validation group (C).

The clinical decision analysis curve (Figure 7C) indicates that utilizing the BI-RADS model, SE+ABVS radiomics model, and nomogram for decision-making significantly improved the net benefit for patients compared to the assumption of intervention for all lesions or no intervention at all. Furthermore, the nomogram provided a greater net benefit to patients compared to both the BI-RADS model and SE+ABVS radiomics model.

4 Discussion

The study combined radiomics features of ABVS and SE images with conventional imaging diagnosis criteria along with clinical risk factors for developing a clinical-radiomics nomogram that demonstrated excellent diagnostic efficacy, as well as good calibration capabilities, and significant clinical usefulness.

Although ABVS and SE examination techniques offer significant advantages in breast screening, the examiner’s naked eye remains incapable of capturing deep image information. Radiomics provides a pathway to capture internal tumor information at a more profound level. Wang et al. derived radiomics features from ABVS images and constructed multiple machine learning models for breast cancer diagnosis, the best of which was the support vector machine model with an AUC of 0.857 (17). Additionally, Liu et al. employed radiomics features extracted from SE images for breast cancer prediction, yielding a Radscore with an AUC of 0.866 in the test set (18). Besides, Ma et al. developed a multivariate logistic model by combining SE, B-mode, and ABVS coronal radiomics features, with an AUC value of 0.946 in the internal validation group (19). In this study, we performed radiomics analysis on both ABVS and SE images. For ABVS images, we delineated the ROI across sagittal, transverse, and coronal planes. While outlining the ROI in the ABVS coronal planes and the SE image, we incorporated a portion of the lesion’s peripheral tissues to capture additional information. As a result, the ABVS and SE features that we acquired demonstrated good predictive capabilities for breast cancer, and the combination of the two yielded a higher diagnostic efficacy than the BI-RADS model that obtained by a highly experienced physician based on visual assessment alone (AUC: 0.933 vs. 0.916). Additionally, this study analyzed clinical risk factors related to breast cancer and revealed that age and lesion size exhibited significantly higher values in the malignant group compared to the benign group, which is consistent with previous research findings (3335). Therefore, we developed a nomogram by integrating Radscore, patient’s age, maximum diameter of the lesion, and BI-RADS scores using multivariate logistic regression analysis. The AUC of this nomogram in the internal validation group was 0.964, which surpassed that of both the SE+ABVS radiomics model and the clinical model. Furthermore, we conducted an analysis on the clinical utility of this nomogram, and the decision analysis curves revealed that it could offer superior net benefit to patients across a broad range of threshold intervals. Consequently, the nomogram holds significant value as a point of reference for clinicians, particularly novice practitioners lacking diagnostic expertise in identifying suspicious lesions.

In addition, this nomogram has demonstrated significant advantages in the diagnosis of BI-RADS category 4 lesions. The appearance of these lesions on imaging can be highly deceptive, so they span a wide range of malignancy risks (36, 37), which makes clinical diagnosis extremely challenging, often necessitating biopsies to definitively determine the nature of such lesions (30). However, routine biopsy results are often influenced by the spatial heterogeneity of the lesion and operator expertise (38), while also being an invasive procedure with potential complications such as bleeding (39). The majority of radiomics studies for this category of lesions have predominantly utilized MR images, and these studies have yielded favorable outcomes (2629). Nevertheless, MR examinations are expensive, time-consuming, and not suitable for common screenings (40, 41). Based on ABVS images, Wang et al. integrated clinical ultrasound factors and Radscore to develop a nomogram for the diagnosis of BI-RADS category 4 lesions, which achieved an AUC value of 0.925 in the internal validation group and effectively minimized unnecessary biopsies (20). During this study, we constructed nomogram that also achieved an AUC value of 0.930 for the diagnosis of BI-RADS category 4 lesions in the validation group, surpassing the performance of the clinical model (AUC: 0.839), thereby further validating its good diagnostic efficacy. This may be attributed to the fact that the radiomics features selected for this study are all texture features with the majority derived from wavelet transform. Previous studies have demonstrated the value of wavelet transform-based texture features for the diagnosis of tumor lesions (42). The primary advantage of wavelet transform in image analysis lies in its multi-scale analysis capability, allowing it to capture the texture information of an image at various granularities. It possesses directional sensitivity, enabling it to accurately identify texture changes in multiple directions, while its time-frequency localization property allows it to keenly detect local variations in images. Additionally, wavelet transform can enhance image contrast, exhibit certain resistance to noise, and effectively compress image information, making feature extraction more robust and efficient (43). By quantifying the textural variances of breast lesions, we successfully captured the subtle heterogeneity within these lesions, thereby effectively distinguished between benign and malignant breast lesions.

Lambin et al. introduced the radiomics quality score (RQS) to provide a framework for clinical researchers to evaluate and guide their radiomics studies (44). This study has given comparatively detailed elaboration on image acquisition, feature extraction and screening, and model construction in order to ensure the reproducibility of the study. Two physicians independently delineated the lesions, effectively achieving multiple segmentations. The features extracted from both segmentations were then subjected to ICC analysis. Consequently, only the features demonstrating excellent repeatability and robustness were selected. To prevent model overfitting, we standardized the feature values. Features with strong discriminative ability were obtained through the U-test, LASSO regression with tenfold cross-validation from ABVS and SE images, respectively. Subsequently, we evaluated the constructed nomogram using calibration curves and the Hosmer-Lemeshow test, which demonstrated excellent calibration performance. Based on these analyses, the nomogram appears to be a robust and generalizable tool, offering accurate risk prediction with potential for practical clinical implementation. Although the development of this nomogram necessitates a combination of diverse factors, these data can be retrospectively obtained without imposing an additional examination burden on patients.

Admittedly, This study was subject to certain limitations: It was conducted as a single-center retrospectively study thus selection bias may have occurred, and lacked external validation, which necessitates further multicenter large-sample studies and prospective trials for the validation of our developed nomogram.

5 Conclusion

The nomogram developed in this study, which combined SE and ABVS radiomics features, with traditional imaging assessment criteria and clinical risk factors, it can serve as a reliable and non-invasive analytical tool to assist physicians in accurately assessing the malignancy risk in solid hypoechoic breast lesions, leading to better clinical decision-making.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by Institutional Review Board of the Third Xiangya Hospital, CSU. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

GS: Conceptualization, Formal Analysis, Methodology, Writing – original draft, Writing – review & editing. JL: Formal Analysis, Methodology, Writing – review & editing. YY: Data curation, Writing – review & editing. ZY: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Supervision, Writing – review & editing.

Funding

This work was supported by grants from the Natural Science Foundation of Hunan Province, China (No. 2022JJ40760) and Project of Hunan Provincial Health Commission (No.202209024823).

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.

References

1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. 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. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

2. AlSendi M, O'Reilly D, Zeidan YH, Kelly CM. Oligometastatic breast cancer: Are we there yet? Int J Cancer (2021) 149(8):1520–8. doi: 10.1002/ijc.33693

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Liang Y, Zhang H, Song X, Yang Q. Metastatic heterogeneity of breast cancer: Molecular mechanism and potential therapeutic targets. Semin Cancer Biol (2020) 60:14–27. doi: 10.1016/j.semcancer.2019.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Cedolini C, Bertozzi S, Londero AP, Bernardi S, Seriau L, Concina S, et al. Type of breast cancer diagnosis, screening, and survival. Clin Breast Cancer (2014) 14(4):235–40. doi: 10.1016/j.clbc.2014.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Yang R, Han Y, Yi W, Long Q. Autoantibodies as biomarkers for breast cancer diagnosis and prognosis. Front Immunol (2022) 13:1035402. doi: 10.3389/fimmu.2022.1035402

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Zhi H, Ou B, Luo BM, Feng X, Wen YL, Yang HY. Comparison of ultrasound elastography, mammography, and sonography in the diagnosis of solid breast lesions. J ultrasound Med (2007) 26(6):807–15. doi: 10.7863/jum.2007.26.6.807

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Balleyguier C, Ciolovan L, Ammari S, et al. Breast elastography: the technical process and its applications. Diagn interventional Imaging (2013) 94(5):503–13. doi: 10.1016/j.diii.2013.02.006

CrossRef Full Text | Google Scholar

8. Golatta M, Franz D, Harcos A, Junkermann H, Rauch G, Scharf A, et al. Interobserver reliability of automated breast volume scanner (ABVS) interpretation and agreement of ABVS findings with hand held breast ultrasound (HHUS), mammography and pathology results. Eur J Radiol (2013) 82(8):e332–6. doi: 10.1016/j.ejrad.2013.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Chen L, Chen Y, Diao XH, Fang L, Pang Y, Cheng AQ, Fang L, Pang Y, Cheng AQ, et al. 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. doi: 10.1016/j.ultrasmedbio.2013.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Zheng FY, Yan LX, Huang BJ, Xia HS, Wang X, Lu Q, et al. 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. doi: 10.1016/j.ejrad.2015.07.028

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Xu C, Wei S, Xie Y, Guan X, Fu N, Huang P, et al. Combined use of the automated breast volume scanner and the US elastography for the differentiation of benign from Malignant lesions of the breast. BMC Cancer (2014) 14:798. doi: 10.1186/1471-2407-14-798

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Guo R, Lu G, Qin B, Fei B. Ultrasound imaging technologies for breast cancer detection and management: A review. Ultrasound Med Biol (2018) 44(1):37–70. doi: 10.1016/j.ultrasmedbio.2017.09.012

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology (2016) 278(2):563–77. doi: 10.1148/radiol.2015151169

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Lambin P, Rios-Velazquez E, Leijenaar R, Carvalho S, van Stiphout RG, Granton P, et al. Radiomics: extracting more information from medical images using advanced feature analysis. Eur J Cancer (Oxford England: 1990) (2012) 48(4):441–6. doi: 10.1016/j.ejca.2011.11.036

CrossRef Full Text | Google Scholar

15. Sala E, Mema E, Himoto Y, Veeraraghavan H, Brenton JD, Snyder A, et al. Unravelling tumour heterogeneity using next-generation imaging: radiomics, radiogenomics, and habitat imaging. Clin Radiol (2017) 72(1):3–10. doi: 10.1016/j.crad.2016.09.013

PubMed Abstract | CrossRef Full Text | 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. doi: 10.1038/s41598-019-48488-4

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wang H, Yang X, Ma S, Zhu K, Guo S. An optimized radiomics model based on automated breast volume scan images to identify breast lesions: comparison of machine learning methods: comparison of machine learning methods. J Ultrasound Med (2022) 41(7):1643–55. doi: 10.1002/jum.15845

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Liu X, Zhang J, Zhou J, He Y, Xu Y, Zhang Z, et al. Multi-modality radiomics nomogram based on DCE-MRI and ultrasound images for benign and Malignant breast lesion classification. Front Oncol (2022) 12:992509. doi: 10.3389/fonc.2022.992509

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Ma Q, Shen C, Gao Y, Duan Y, Li W, Lu G, et al. Radiomics analysis of breast lesions in combination with coronal plane of ABVS and strain elastography. Breast Cancer (Dove Med Press) (2023) 15:381–90. doi: 10.2147/BCTT.S410356

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Wang SJ, HQ L, Yang T, Huang MQ, Zheng BW, Wu T, et al. Automated breast volume scanner (ABVS)-based radiomic nomogram: A potential tool for reducing unnecessary biopsies of BI-RADS 4 lesions. Diagnostics (Basel) (2022) 12(1):172. doi: 10.3390/diagnostics12010172

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Lin F, Wang Z, Zhang K, Yang P, Ma H, Shi Y, et al. Contrast-enhanced spectral mammography-based radiomics nomogram for identifying benign and Malignant breast lesions of sub-1 cm. Front Oncol (2020) 10:573630. doi: 10.3389/fonc.2020.573630

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Parekh VS, Jacobs MA. Integrated radiomic framework for breast cancer and tumor biology using advanced machine learning and multiparametric MRI. NPJ Breast Cancer (2017) 3:43. doi: 10.1038/s41523-017-0045-3

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zhang Q, Xiao Y, Suo J, Shi J, Yu J, Guo Y, et al. Sonoelastomics for breast tumor classification: A radiomics approach with clustering-based feature selection on sonoelastography. Ultrasound Med Biol (2017) 43(5):1058–69. doi: 10.1016/j.ultrasmedbio.2016.12.01

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Xu Z, Wang Y, Chen M, Zhang Q. Multi-region radiomics for artificially intelligent diagnosis of breast cancer using multimodal ultrasound. Comput Biol Med (2022) 149:105920. doi: 10.1016/j.compbiomed.2022.105920

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Ou X, Wang J, Zhou R, Zhu S, Pang F, Zhou Y, et al. Ability of F-FDG PET/CT radiomic features to distinguish breast carcinoma from breast lymphoma. Contrast Media Mol Imaging (2019) 2019:4507694. doi: 10.1155/2019/4507694

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Zhang R, Wei W, Li R, Li J, Zhou Z, Ma M, et al. An MRI-based radiomics model for predicting the benignity and Malignancy of BI-RADS 4 breast lesions. Front Oncol (2022) 11:733260. doi: 10.3389/fonc.2021.733260

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Cui Q, Sun L, Zhang Y, Zhao Z, Li S, Liu Y, et al. Value of breast MRI omics features and clinical characteristics in Breast Imaging Reporting and Data System (BI-RADS) category 4 breast lesions: an analysis of radiomics-based diagnosis. Ann Transl Med (2021) 9(22):1677. doi: 10.21037/atm-21-5441

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lyu Y, Chen Y, Meng L, Guo J, Zhan X, Chen Z, et al. Combination of ultrafast dynamic contrast-enhanced MRI-based radiomics and artificial neural network in assessing BI-RADS 4 breast lesions: Potential to avoid unnecessary biopsies. Front Oncol (2023) 13:1074060. doi: 10.3389/fonc.2023.1074060

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Hao W, Gong J, Wang S, Zhu H, Zhao B, Peng W. Application of MRI radiomics-based machine learning model to improve contralateral BI-RADS 4 lesion assessment. Front Oncol (2020) 10:531476. doi: 10.3389/fonc.2020.531476

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Mendelson EB, Böhm-Vélez M, Berg WA, et al. ACR BIRADS® Ultrasound, in: ACR BI-RADS® Atlas, Breast Imaging Reporting and Data System. American College of Radiology, Reston, 2013 VA. Available at: https://www.acr.org/Clinical-Resources/Reporting-and-Data-Systems/Bi-rads (Accessed 11 Jan 2023).

Google Scholar

31. Zhi H, Xiao XY, Ou B, Zhong WJ, Zhao ZZ, Zhao XB, et al. Could ultrasonic elastography help the diagnosis of small (≤2 cm) breast cancer with the usage of sonographic BI-RADS classification? Eur J Radiol (2012) 81(11):3216–21. doi: 10.1016/j.ejrad.2012.04.016

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Büsing KA, Kilian AK, Schaible T, Debus A, Weiss C, Neff KW. Reliability and validity of MR image lung volume measurement in fetuses with congenital diaphragmatic hernia and in vitro lung models. Radiology (2008) 246(2):553–61. doi: 10.1148/radiol.2462062166

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Marzbani B, Nazari J, Najafi F, Marzbani B, Shahabadi S, Amini M, et al. Dietary patterns, nutrition, and risk of breast cancer: a case-control study in the west of Iran. Epidemiol Health (2019) 41:e2019003. doi: 10.4178/epih.e2019003

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Sun YS, Zhao Z, Yang ZN, Xu F, Lu HJ, Zhu ZY, et al. Risk factors and preventions of breast cancer. Int J Biol Sci (2017) 13(11):1387–97. doi: 10.7150/ijbs.21635

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Gity M, Arabkheradmand A, Taheri E, Shakiba M. Diagnostic investigation of breast magnetic resonance imaging in Malignant and benign mass lesions. Arch Med Sci (2018) 14(5):1061–9. doi: 10.5114/aoms.2016.62281

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Zou X, Wang J, Lan X, Lin Q, Han F, Liu L, Lin Q, Han F, Liu L, et al. Assessment of diagnostic accuracy and efficiency of categories 4 and 5 of the second edition of the BI-RADS ultrasound lexicon in diagnosing breast lesions. Ultrasound Med Biol (2016) 42(9):2065–71. doi: 10.1016/j.ultrasmedbio.2016.04.020

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Liu G, Zhang MK, He Y, Liu Y, Li XR, Wang ZL. BI-RADS 4 breast lesions: could multi-mode ultrasound be helpful for their diagnosis? Gland Surg (2019) 8(3):258–70. doi: 10.21037/gs.2019.05.01

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Teng R, Wei Q, Zhou J, Dong M, Jin L, Hu W, et al. The influence of preoperative biopsy on the surgical method in breast cancer patients: a single-center experience of 3,966 cases in China. Gland Surg (2021) 10(3):1038–45. doi: 10.21037/gs-21-7

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Nakano S, Imawari Y, Mibu A, Otsuka M, Oinuma T. Differentiating vacuum-assisted breast biopsy from core needle biopsy: Is it necessary? Br J Radiol (2018) 91(1092):20180250. doi: 10.1259/bjr.20180250

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Wang L. Early diagnosis of breast cancer. Sensors (Basel) (2017) 17(7):1572. doi: 10.3390/s17071572

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Raber B, VJ B, Bedrosian I. How does MR imaging help care for my breast cancer patient? Perspective of a surgical oncologist. Magn Reson Imaging Clin N Am (2018) 26(2):281–8. doi: 10.1016/j.mric.2017.12.010

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Jiang Z, Yin J, Han P, Chen N, Kang Q, Qiu Y, et al. Wavelet transformation can enhance computed tomography texture features: a multicenter radiomics study for grade assessment of COVID-19 pulmonary lesions. Quant Imaging Med Surg (2022) 12(10):4758–70. doi: 10.21037/qims-22-252

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Sudarshan VK, Mookiah MR, Acharya UR, Chandran V, Molinari F, Fujita H, et al. Application of wavelet techniques for cancer diagnosis using ultrasound images: A Review. Comput Biol Med (2016) 69:97–111. doi: 10.1016/j.compbiomed.2015.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong EEC, van Timmeren J, et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat Rev Clin Oncol (2017) 14(12):749–62. doi: 10.1038/nrclinonc.2017.141

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: nomogram, breast, radiomics, automated breast volume scanner, strain elastography

Citation: Shiyan G, Liqing J, Yueqiong Y and Yan Z (2023) A clinical-radiomics nomogram based on multimodal ultrasound for predicting the malignancy risk in solid hypoechoic breast lesions. Front. Oncol. 13:1256146. doi: 10.3389/fonc.2023.1256146

Received: 10 July 2023; Accepted: 27 September 2023;
Published: 17 October 2023.

Edited by:

Luca Nicosia, European Institute of Oncology (IEO), Italy

Reviewed by:

Dania Cioni, University of Pisa, Italy
Zhiyong Zhou, Chinese Academy of Sciences (CAS), China

Copyright © 2023 Shiyan, Liqing, Yueqiong and Yan. 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: Zhang Yan, 546940642@qq.com

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.