Skip to main content

Preoperative ultrasound radiomics analysis for expression of multiple molecular biomarkers in mass type of breast ductal carcinoma in situ

Abstract

Background

The molecular biomarkers of breast ductal carcinoma in situ (DCIS) have important guiding significance for individualized precision treatment. This study was intended to explore the significance of radiomics based on ultrasound images to predict the expression of molecular biomarkers of mass type of DCIS.

Methods

116 patients with mass type of DCIS were included in this retrospective study. The radiomics features were extracted based on ultrasound images. According to the ratio of 7:3, the data sets of molecular biomarkers were split into training set and test set. The radiomics models were developed to predict the expression of estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), Ki67, p16, and p53 by using combination of multiple feature selection and classifiers. The predictive performance of the models were evaluated using the area under the curve (AUC) of the receiver operating curve.

Results

The investigators extracted 5234 radiomics features from ultrasound images. 12, 23, 41, 51, 31 and 23 features were important for constructing the models. The radiomics scores were significantly (P < 0.05) in each molecular marker expression of mass type of DCIS. The radiomics models showed predictive performance with AUC greater than 0.7 in the training set and test set: ER (0.94 and 0.84), PR (0.90 and 0.78), HER2 (0.94 and 0.74), Ki67 (0.95 and 0.86), p16 (0.96 and 0.78), and p53 (0.95 and 0.74), respectively.

Conclusion

Ultrasonic-based radiomics analysis provided a noninvasive preoperative method for predicting the expression of molecular markers of mass type of DCIS with good accuracy.

Peer Review reports

Background

Breast ductal carcinoma in situ (DCIS) is a kind of malignant tumor originated in the ductal epithelial tissue, limited to the basement membrane [1]. DCIS is the second most common breast tumor, and accounts for approximately 20–30% [2]. Some DCIS had the potential to further develop into breast invasive cancer [3]. The clinical treatments of patients with DCIS include surgical resection, radiotherapy, chemotherapy and endocrine therapy, in which surgical resection includes simple focal resection and mastectomy, with different therapeutic effects [4]. Although the prognosis of DCIS is good, more than 14% of DCIS patients may develop invasive cancer without treatment within 10 years [5]. In the past 10 years, the incidence of DCIS has gradually increased, highlighting the understanding the importance of DCIS pathology [6]. However, the pathologic mechanism of the transition from DCIS to invasive carcinoma is still unclear, which produces clinical challenges of overdiagnosis and overtreatment in patients with DCIS [7]. Therefore, the investigators thought more studies were need to understand the potential of the pathological process of DCIS, in order to adapt to the current individualized, refined treatment.

Immunohistochemistry (IHC) can reflect the expression of molecular biomarkers in tumor tissue, which can further clarify the biological behaviors of tumors. The expression of different molecular biomarkers can lead to different biological behaviors and treatments. Some studies have shown that some molecular biomarkers were important indicators for predicting biological behavior and judging follow-up treatment in patients with DCIS, such as estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), Ki67, p16, and p53. ER and PR are the earliest molecular biomarkers of breast cancer. They are predictors of breast cancer prognosis and endocrine adjuvant therapy [8]. HER2 is a proto-oncogene, which is mainly involved in tumor signal transduction and cell proliferation. Its positive expression can lead to a high distant metastasis rate and poor prognosis of breast cancer. Ki67 is an antigenic nuclear protein that can be used as a proliferation marker. Its high expression is considered to be a biomarker of tumor invasion [9]. Ki67 has a good application prospect in predicting endocrine therapy response of breast cancer [10]. Defined as a tumor suppressor gene, p16 is considered to be an important cell cycle regulator [11]. P16 is closely related to abnormal methylation initiation. P53 is a common tumor suppressor gene. Impaired function of p53, such as p53 mutation, can lead to uncontrolled proliferation of damaged cells [12]. Therefore, accurate identification of the expression of molecular biomarkers can help stratify tumor risk and facilitate the development of personalized and accurate treatment plans.

Currently, the preoperative evaluation of the molecular biomarkers of DCIS mainly depends on IHC detection after biopsy. However, because the progression of tumors are dynamic process, there are differences in spatio-temporal evolution. In addition, the evaluation results of a few tissue biopsies do not necessarily represent the expression of the molecular biomarkers of the whole tumor [13]. Invasive procedures and potential risks limit its multiple applications in monitoring tumor progression and biological behavior. However, the preoperative monitoring of molecular biomarkers can dynamically identify the progression of tumors and the changes in biological behavior, which has great significance for the accurate formulation of treatment plans and the evaluation of curative effects. To avoid overdiagnosis and overtreatment of patients with DCIS, it is necessary to provide dynamic and accurate evaluation of biological behavior information for the clinic.

With the breakthrough of imaging technology, mammography is an important examination method for DCIS, which is sensitive to the detection of calcification [14]. Ultrasound (US) has become the main examination technology to detecting breast lesions [15], which is real-time, dynamic and non-invasive. There are two types of DCIS: mass and non-mass. Some studies suggested that the detection rate of US in 93 patients with mass type DCIS reached 77.4% [16]. The main characteristics of DCIS in ultrasound were: uneven low or slightly low echo, irregular shape, unclear borders, parallel skin, weakened posterior echo, calcification, and some blood flow signals [17, 18]. However, mammography or US-assisted screening could increase overdiagnosis because both tests primarily detect low-grade invasive cancers [19]. At the same time, radiologists are very time-consuming to accumulate experience and have strong personal subjectivity, which is another problem that needs to be solved. There is an urgent need for more advanced imaging evaluation methods to guide the diagnosis and treatment of DCIS.

Breast lesions are diagnosed and screened by various imaging methods, such as mammography, US, and magnetic resonance imaging (MRI). All three examination have some limitations [20]. Radiomics is a hot subject of artificial intelligence that is applied in the medical imaging field, which is the cornerstone of precision science in the future. Radiomics is defined as the extraction of high-throughput features from single or multiple medical image patterns to select features that are closely associated with tumors,and the ultimate goal is to construct prediction models based on features to provide accurate tumor phenotypic analysis information and accurate treatment decision-making [21]. Radiomics highlights the image features that are not visible to the naked eye, thus significantly enhancing the predictive power of medical imaging [22]. Radiomics has been developed in a wide range of fields, such as disease diagnosis and biological behavior judgment. For example, the US-radiomics model developed by Luo WQ et al. had better performance in distinguishing breast lesions than breast imaging reporting and data system (BI-RADS) [23]. Lin F et al. found that the radiomics score was more effective than the clinical radiological model in benign and malignant breast lesions (< 1 cm) [24]. These series of studies showed that radiomics had better performance than traditional imaging features in the diagnosis of breast diseases to some extent.

Thus, this retrospective study intended to further clarify the relationship between US-radiomics and molecular markers of DCIS. Radiomics models had been developed to noninvasively evaluate the expression of molecular markers to help achieve accurate risk stratification and treatment for patients with DCIS.

Materials and methods

Study cohort

Clinical data of 400 patients with DCIS who were pathologically confirmed by surgery were retrospectively analyzed by the investigators. The data were based on the pathology reports of the first affiliated hospital of Guangxi medical university from January 2015 to July 2020. Further inclusion and exclusion criteria for this cohort study were as follows. Inclusion criteria: (1) primary breast DCIS, (2) IHC results of molecular biomarkers; and (3) preoperative US data within one month. Exclusion criteria: (1) non-mass DCIS, including manifestations of ductal dilation, diffuse calcification, and diffuse distribution of lesions; (2) unclear image of target lesions; (3) secondary DCIS or postoperative recurrence of DCIS; (4) preoperative treatment history of radiotherapy, chemotherapy and traditional Chinese medicine; and (5) lack of clinical data.

This study finally enrolled a total of 116 patients with DCIS. The IHC (ER, PR, HER2, Ki67, p16, and p53) conformed to the diagnositic criteria of the department pathology in this hospital, and were classified as positive or negative. The IHC results of Ki67 were positive or high expression (Ki67 >  = 14%) and negative or low expression (Ki67 < 14%) [25]. The number of patients enrolled for each molecular biomarker were listed as follows: 112 cases (ER), 109 cases (PR), 94 cases (HER2), 107 cases (Ki67), 74 cases (p16), and 116 cases (p53) (Fig. 1).

Fig. 1
figure 1

Study cohort. a Workflow of study cohort inclusion. b Up-set plot of the expression of molecular markers shared between different samples

Image collection and tumor segmentation

Each US radiologist involved in image collection had over 5 years of experience in the field of breast. Before collecting US data, all radiologists were strictly trained. GE Logiq E9 (GE Healthcare, United States), Aloka EZU-MT28-S1 (Aloka, Japan) and MYLAB CLASS C (MYLAB, Italy) medical ultrasound diagnositic instruments were utilized for image collection. The breast probe was selected, and the frequency was set to 7–14 MHz. The patients took the supine position, put their hands on the head, and fully exposed the breast area and armpits on both sides. The lesions were scanned from multiple angles, and the largest clear image of the lesions were selected. The following ultrasonic characteristics of the lesions were recorded: BI-RADS classification, location, size, shape, boundary, internal echo, calcification, posterior echo changes, ductal dilatation, blood flow signal distribution and axillary lymph nodes.

These images were imported into the ITKSNAP software (version 3.8.0). To avoid subjective compliance, two radiologists with five years of working experience manually delineated the region of interest (ROI) of the lesions. The radiologists disregarded the diagnosis and pathological results of the patients [26]. After the discussion, when there was a big difference between the two radiologists, the third radiologist with 10 years of experience re-examined and confirmed the final boundary. This process provided reliable DCIS area contours and ensured the accuracy of feature extraction.

Image pre-processing and feature extraction

The Intelligence Foundry software (version 1.3, GE Healthcare, Shanghai, China) was applied for radiomics analysis. Figure 2 summarized the main flow of radiomics analysis. The software relied on algorithms provided by the Pyradiomics package that comply with the image biomarker standardization initiative (IBSI, version 2016) [27]. Features were automatically calculated and extracted by the Pyradiomics extractor. The maximum number of features extraction of the software was 5234, including: 122 original, 48 intraperinodular textural transition (ipris), 468 co-occurrence of local anisotropic gradient orientations (CoLIAGe), 432 wavelets + local binary pattern (LBP), 2,944 shearlets, 1,080Gabors, 80 phased congruency-based local binary pattern (PLBP) and 60 wavelet-based improved local binary pattern (WILBP) features (Additional file 1). Before feature extraction, the images were pre-processed: the gray value of the image was discretized with a bin size of 256, and the original features were extracted. The features of wavelets + LBP, Shearlets, Gabors, PLBP and WILBP were extracted by wavelet transform, shearlet transform and garber operator transform on the gray value matrix of the original images, respectively [28] (Fig. 3).

Fig. 2
figure 2

Workflow of radiomics analysis

Fig. 3
figure 3

The process of quantifying features. a Delineation of the ROIs. b Gray level co-occurrence matrix (GLCM), run length matrices (RLM), and histogram feature extraction. c The classification of 5234 features

Data grouping and data cleaning

To balance the initial distribution of data, each sub-data set was randomly split into training set and test set in a ratio of 7:3. Based on the difference in the image extraction feature quantization caused by different medical ultrasound diagnositic instruments and parameters, the combat method was employed to solve this problem. The combat method could be used to coordinate and correct the differences between different machines and different center images. Some studies had applied this method to the MRI images [29]. In addition, the median value of the feature quantization value was applied to fill the missing sample. The min–max normalization method was employed to normalize the feature data to improve the comparability between features. It converted the original data to the range of [0, 1] by linearization, which realized the proportional scaling of the original data.

Feature importance

The purpose of the study is difficult to explain with thousands of radiomics features of high-dimensional data. Feature importance analysis helps to explain the importance of features for subsequent model constructing. Multiple combination techniques were applied to explain the importance of features: First, Spearman correlation coefficient test was used to eliminate high correlation features with threshold values (0.75, 0.85, 0.95). This test was a statistical index to measure the correlation between two variables. Three dimensionality reduction methods (least absolute shrinkage and selection operator (LASSO) [30], random forests (RF) [31], and support vector machine-recursive feature elimination (SVM-RFE) [32]) separated or jointed statistical tests for selecting the important features. In the statistical test, if the data accorded with the normal distribution, the t-test was adopted; otherwise, the Mann–Whitney U test was adopted.

Predictive radiomics models

Machine learning algorithms were developed based on Python environment. Five machine-learning-based classifiers (decision tree (DT), k-nearest neighbors (KNN), logistics regression (LR), naive Bayes (NB), and support vector machine (SVM)) were employed to predict the expression levels of the molecular biomarkers of DCIS [33, 34], and the score of each model was calculated. In addition, the fivefold cross-validation method was explored to improve the accuracy of the models. The test set was used to evaluate the reliability of the models.

To accurately evaluate the predictive ability of radiomics models, the receiver operating curve (ROC), the area under the curve (AUC), accuracy (ACC), precision (PREC), sensitivity (Sn) and specificity (Sp) were adopted for the evaluation. The closer the AUC was to 1, the higher the diagnostic efficiency was. In this study, only the best classification results of the classifier were shown.

Results

Patient characteristics and molecular biomarkers of interest

The mean age of all patients was 48.8 ± 11.1 years, and the age range was 29–84 years. The characteristics parameters were shown in Table 1. The ultrasonographic features of the patients were similar to those reported in the literatures. The expression of ER, PR, HER2, Ki67, p16, and p53 were as follows: 49 patients with ER-negatives and 63 patients with ER-positives; 53 patients with PR-negatives and 59 patients with PR-positives; 36 patients with HER2-negatives and 58 patients with HER2-positives; 45 patients with Ki67-negatives and 62 patients with Ki67-positives; 29 patients with p16-negatives and 45 patients with p16-positives; 34 patients with p53-negatives and 82 patients with p53-positives.

Table 1 Patient characteristics and molecular biomarkers of interest

Radiomics analysis

The correlation clustering heatmaps among 5234 features of each molecular biomarkers (ER, PR, HER2, Ki67, p16, and p53) were shown in Fig. 4. A list of 18 feature importance methods were obtained, and the combination feature selection methods for optimal modeling results of each molecular biomarkers were as follows: Spearman0.75 + Statistical Test + RF, Spearman0.75 + Statistical Test + RF, Spearman0.75 + LASSO, Spearman0.75 + Statistical Test + RF, Spearman0.75 + Statistical Test + SVM-RFE, and Spearman0.85 + SVM-RFE. 12 features, 23 features, 41 features, 20 features, 31 features and 23 features were important for constructing prediction models. The heatmaps of the model features were presented in Fig. 5.

Fig. 4
figure 4

Correlation cluster analysis of 5234 radiomics features. The Pearson correlation test was used to analyze the correlation between features, and the "pheatmap" R software package was applied to draw heat maps. a ER; b PR; c HER2; d Ki67; e p16; f p53

Fig. 5
figure 5

Important features for each molecular biomarkers. a ER; b PR; c HER2; d Ki67; e p16; f p53

Ninety models were obtained by constructing prediction models with five classifiers, and the performance of the models were presented in Fig. 6 (Additional file 2). The optimal radiomics models were constructed by DT, SVM, KNN, SVM, KNN and KNN classifiers, respectively, and showed above moderate predictive performance in predicting the expression of molecular markers of DCIS (Table 2). Radiomics scores of training set and test set were significantly different in each molecular marker expression (training set, P < 0.001, test set, P < 0.05). The predictive performance of the radiomics models of each molecular biomarker in the training set: ER (AUC, 0.94, 95% confidence interval (CI) 0.89–0.99), PR (AUC, 0.90, 95% CI 0.83–0.97), HER2 (AUC, 0.94, 95% CI 0.89–0.99), Ki67 (AUC, 0.95, 95% CI 0.90–0.99), p16 (AUC, 0.96, 95% CI 0.91–1.00), p53 (AUC, 0.95, 95% CI 0.90–0.99), respectively (Fig. 7). The calibration curve of the prediction models in the training set confirmed the better consistency of the models (Fig. 8). The radiomics models showed predictive performance with AUC greater than 0.7 in the test set: ER (AUC, 0.84, 95% CI 0.68–0.99), PR (AUC, 0.78, 95% CI, 0.60–0.96), HER2 (AUC, 0.74, 95% CI 0.74–0.99), Ki67 (AUC, 0.86, 95% CI 0.67–0.97), p16 (AUC, 0.78, 95% CI 0.59–0.97), p53 (AUC, 0.74, 95% CI 0.55–0.93), respectively (Fig. 9).

Fig. 6
figure 6

Heat maps of evaluation indicators for ninety radiomics prediction models. a ER; b PR; c HER2; d Ki67; e p16; f p53

Table 2 Evaluation of radiomics models in each DCIS molecular biomarkers
Fig. 7
figure 7

Performance of the radiomics models in the training set. a ER; b PR; c HER2; d Ki67; e p16; f p53

Fig. 8
figure 8

Calibration curves of the radiomics models in the training set. The oblique dashed line represents the perfect prediction of the ideal model. The solid line represents the performance of the radiomics model, and the dotted line near the diagonal indicates a better prediction. a ER; b PR; c HER2; d Ki67; e p16; f p53

Fig. 9
figure 9

Performance of the radiomics models in the test set. a ER; b PR; c HER2; d Ki67; e p16; f p53

Discussion

This study was the first non-invasive comprehensive analysis based on US-radiomics to predict the expression of molecular markers of DCIS. The investigators recruited only 116 patients with DCIS for this study, but it was exciting to see that the radiomics models showed more than moderate predictive performance in predicting molecular biomarker expression of DCIS.

DCIS is a malignant tumor with good prognosis, but it is heterogeneous in morphology and genetics. Before the imaging examination was performed, the diagnosis of DCIS was only due to the appearance of nipple discharge and/or palpable mass symptoms, which accounted for only 2% of DCIS detected. It showed that DCIS with hidden symptoms were easily missed [35]. With the screening of imaging technology (mammography, US and MRI), the detection rate of DCIS had gradually increased. This detection rate included symptomatic DCIS, and whether there was overdiagnosis in the detection of insidious DCIS was also a hot topic of controversy [36], Unfortunately, the diagnosis of DCIS marked women as at risk of invasive breast cancer, so women diagnosed with DCIS may suffer serious psychological distress, leading to the progression of DCIS [37]. In addition, the current treatment methods were also facing the controversy over the treatment of some patients [38]. Therefore, the main clinical challenge in DCIS has been to distinguish between patients who have a better chance of developing invasive cancer and require more treatment and those who are less likely to develop DCIS and need less or no treatment [39]. Immunohistochemical markers can explain the changes in the biological behaviors of tumors on the molecular level. More and more studies have pointed out the changes in molecular markers associated with the progression of DCIS to invasive cancer [40]. For example, Zhang GJ et al. [41] found that 79% of DCIS patients were positive for P53 when studying the occurrence and development of breast cancer. Davis et al. [42] demonstrated that high Ki67 expression was an independent predictor of postoperative recurrence in patients with DCIS. Cornfield DB had found a higher recurrence rate with PR > 3.5% using tree structure survival [43]. The results showed that the changes in the biological behavior of DCIS were closely related to the expression of molecular biomarkers.

About various imaging technologies, they also have application limitations [44]. Mammography is the main method of early breast cancer detection, but it is closely related to the density of the lesion and the possibility of covering the lesion [45]. However, Chinese women have dense breasts, so they had certain limitations in finding suspicious lesions in the dense tissues of breasts through mammography [46, 47]. The traditional mammography diagnosis method will cause trauma to the patient to a certain extent and reduce the patient's treatment compliance. Due to its high sensitivity to soft tissues, US can better show lesions in dense glands and has become the primary imaging method for Chinese women to screen and diagnose breast diseases. For non-mass DCIS, US is difficult to recognize [48]. Therefore, this retrospective study only examined mass DCIS, which is a limitation of the study. MRI has considerable advantages in detecting breast lesions, but its specificity is limited by several factors that affect image quality, such as magnetic field and gradient strength, coil performance, contrast agent efficacy and menstrual cycle [49].

Radiomics mainly studies the quantitative features that are related to biology in medical images. Radiomics features are considered the invisible tissue infrastructure components of the object to be imaged, which can serve as a valuable method for studying cancer by imaging, such as MRI. Radiomics can provide in vivo visualization and quantitative analysis of the imaging features of the whole imaging mass. Therefore, radiomics is a precision medical method for non-invasive diagnosis, evaluation of efficacy, biological behavior [50]. Currently, radiomics mainly relies machine learning algorithms to identify meaningful features of image training data set, and for further interpretation of the information and the optimization, so as to accurately predict the content of the research. An independent data set is applied to test the universality of the model and provides feedback for further optimization of the model [51]. To a certain extent, it improves the utilization of image information and enables differential diagnosis of diseases on more subtle levels that cannot be recognized by the naked eye.

Breast radiomics studies are mostly applied to the prediction of the molecular classification, lymph node metastasis and molecular markers of invasive ductal carcinoma. For example, Demircioglu A et al. [52] constructed radiomics models for predicting Ki67 expression in invasive breast cancer based on eight features extracted from MRI images, with an AUC of 0.81. Zhou et al. [53] explored the significance of MRI-radiomics models for predicting the expression of HER2 in patients with invasive breast cancer before surgery; the validation set AUC reached 0.81. There are few reports on DICS with radiomics. However, there are clinical challenges in the diagnosis and treatment of patients with DCIS. Tumor progression and treatment decisions are affected by multiple tumor molecular biomarkers, which require to comprehensively analyze and evaluate the molecular biomarkers of DCIS. To expand the application of radiomics in DCIS, the investigators carried out this study to assess the feasibility of molecular biomarkers of DCIS. The investigators believe that information obtained from multiple molecular biomarkers can help explain the underlying pathological process of DCIS.

In this study, the first highlight, the first comprehensive analysis of molecular markers of DCIS was conducted based on radiomics. Second highlight, there were thousands of radiomics features, including eight classifications: original feature can reflect the number of voxels in the images, intensity distribution, pixel pair frequency, image average gray value, size and shape of ROI (https://pyradiomics.readthedocs.io/en/latest/features.html); Ipris features capture nodular heterogeneity and differential growth patterns; CoLIAGe features can distinguish disease phenotypes that have similar morphologic appearances [54]; wavelets features represent most of the edge information in images; Shearlet features are better for processing high-dimensional signals; Gabors features extract the edge and gradient information of image and reflect the spatial frequency feature; PLBP, and WILBP: PLBP features are an oriented local texture descriptor that combines the phase congruency approach with the LBP. The third highlight of this study was to construct dozens of prediction models by combining multiple classifiers with multiple feature selection to select the optimal prediction results. RF and SVM-RFE had significant performance in feature selection of multiple molecular markers, KNN and SVM classification performed well too. Finally, through the verification of the test set, the prediction models all showed moderate performance.

There were also some shortcomings in our study. First, this retrospective study had the problem of small sample size. It was necessary to increase the sample size or multi-center cooperation to construct universal models. Second, when the radiologists manually delineated the ROIs, there were a certain degree of subjectivity to the contours of the lesions, which may lead to poor robustness of the models. In addition, the delineation process was done by only one radiologist. Third, the investigators only investigated the features extracted from the largest section, which could not represent the whole tumor. Due to the limitations of US, it was not possible to conduct three-dimensional studies similar to other imaging studies.

Conclusion

The application of machine learning-based radiomics analysis provided a non-invasive method for predicting the expression of multiple molecular biomarkers in DCIS, with good prediction performance. This study also demonstrated the potential of radiomics in pathologic assessment and individualized precision therapy.

Availability of data and materials

The datasets supporting the conclusions of this article were included within the article and its additional files.

Abbreviations

DCIS:

Ductal carcinoma in situ

ER:

Estrogen receptor

PR:

Progesterone receptor

HER2:

Human epidermal growth factor receptor-2

AUC:

Area under the curve

CI:

Confidence interval

IHC:

Immunohistochemistry

US:

Ultrasound

MRI:

Magnetic resonance imaging

BI-RADS:

Breast imaging reporting and data system

ROI:

Region of interest

CoLIAGe:

Co-occurrence of local anisotropic gradient orientations

LBP:

Local binary pattern

PLBP:

Phased congruency-based local binary pattern

WILBP:

Wavelet-based improved local binary pattern

GLCM:

Gray level co-occurrence matrix

RLM:

Run length matrices

LASSO:

Least absolute shrinkage and selection operator

RF:

Random forests

SVM-RFE:

Support vector machine-recursive feature elimination

DT:

Decision tree

KNN:

K-nearest neighbors

LR:

Logistics regression

NB:

Naive Bayes

SVM:

Support vector machine

ROC:

Receiver operating curve

ACC:

Accuracy

PREC:

Precision

Sn:

Sensitivity

Sp:

Specificity

References

  1. Wellings S, Jensen HJJI. On the origin and progression of ductal carcinoma in the human breast. Origin Breast Carcinoma. 1973;50(5):1111–8.

    CAS  Google Scholar 

  2. Liu Y, West R, Weber JD, Colditz GA. Race and risk of subsequent aggressive breast cancer following ductal carcinoma in situ. Cancer. 2019;125(18):3225–33.

    Article  CAS  PubMed  Google Scholar 

  3. Villanueva H, Grimm S, Dhamne S, Rajapakshe K, Visbal A, Davis CM, Ehli EA, Hartig SM, Coarfa C, Edwards DP. The emerging roles of steroid hormone receptors in ductal carcinoma in situ (DCIS) of the breast. J Mammary Gland Biol Neoplasia. 2018;23(4):237–48.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Kuerer HM, Albarracin CT, Yang WT, Cardiff RD, Brewster AM, Symmans WF, Hylton NM, Middleton LP, Krishnamurthy S, Perkins GH, et al. Ductal carcinoma in situ: state of the science and roadmap to advance the field. J Clin Oncol. 2009;27(2):279–88.

    Article  PubMed  Google Scholar 

  5. Schnitt SJ. Diagnosis of ductal carcinoma in situ in an era of de-escalation of therapy. Modern. 2020;34:1–7.

    Google Scholar 

  6. Martínez-Pérez C, Turnbull AK, Ekatah GE, Arthur LM, Sims AH, Thomas JS, Dixon JM. Current treatment trends and the need for better predictive tools in the management of ductal carcinoma in situ of the breast. Cancer Treat Rev. 2017;55:163–72.

    Article  PubMed  Google Scholar 

  7. Shah C, Wobb J, Manyam B, Kundu N, Arthur D, Wazer D, Fernandez E, Vicini F. Management of ductal carcinoma in situ of the breast: a review. JAMA Oncol. 2016;2(8):1083–8.

    Article  PubMed  Google Scholar 

  8. Giuliano M, Schettini F, Rognoni C, Milani M, Jerusalem G, Bachelot T, De Laurentiis M, Thomas G, De Placido P, Arpino G, et al. Endocrine treatment versus chemotherapy in postmenopausal women with hormone receptor-positive, HER2-negative, metastatic breast cancer: a systematic review and network meta-analysis. Lancet Oncol. 2019;20(10):1360–9.

    Article  CAS  PubMed  Google Scholar 

  9. Hida AI, Omanovic D, Pedersen L, Oshiro Y, Ogura T, Nomura T, Kurebayashi J, Kanomata N, Moriya T. Automated assessment of Ki-67 in breast cancer: the utility of digital image analysis using virtual triple staining and whole slide imaging. Histopathology. 2020;77(3):471–80.

    Article  PubMed  Google Scholar 

  10. Guarneri V, Dieci MV, Bisagni G, Frassoldati A, Bianchi GV, De Salvo GL, Orvieto E, Urso L, Pascual T, Paré L, et al. De-escalated therapy for HR+/HER2+ breast cancer patients with Ki67 response after 2-week letrozole: results of the PerELISA neoadjuvant study. Ann Oncol. 2019;30(6):921–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Kim M, Katayose Y, Rojanala L, Shah S, Sgagias M, Jang L, Jung YJ, Lee SH, Hwang SG, Cowan KH. Induction of apoptosis in p16INK4A mutant cell lines by adenovirus-mediated overexpression of p16INK4A protein. Cell Death Differ. 2000;7(8):706–11.

    Article  CAS  PubMed  Google Scholar 

  12. Shan M, Zhang X, Liu X, Qin Y, Liu T, Liu Y, Wang J, Zhong Z, Zhang Y, Geng J, et al. P16 and p53 play distinct roles in different subtypes of breast cancer. PLoS ONE. 2013;8(10):e76408.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Juan MW, Yu J, Peng GX, Jun LJ, Feng SP, Fang LP. Correlation between DCE-MRI radiomics features and Ki-67 expression in invasive breast cancer. Oncol Lett. 2018;16(4):5084–90.

    PubMed  PubMed Central  Google Scholar 

  14. Holmberg L, Wong YN, Tabár L, Ringberg A, Karlsson P, Arnesson LG, Sandelin K, Anderson H, Garmo H, Emdin S. Mammography casting-type calcification and risk of local recurrence in DCIS: analyses from a randomised study. Br J Cancer. 2013;108(4):812–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Won SY, Park HS, Kim EK, Kim SI, Moon HJ, Yoon JH, Park VY, Park S, Kim MJ, Cho YU, et al. Survival rates of breast cancer patients aged 40 to 49 years according to detection modality in Korea: screening ultrasound versus mammography. Korean J Radiol. 2020;22:159–67.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Watanabe T, Yamaguchi T, Tsunoda H, Kaoku S, Tohno E, Yasuda H, Ban K, Hirokaga K, Tanaka K, Umemoto T, et al. Ultrasound image classification of ductal carcinoma in situ (DCIS) of the breast: analysis of 705 DCIS lesions. Ultrasound Med Biol. 2017;43(5):918–25.

    Article  PubMed  Google Scholar 

  17. Li W, Zhou Q, Xia S, Wu Y, Fei X, Wang Y, Tao L, Fan J, Zhou W. Application of contrast-enhanced ultrasound in the diagnosis of ductal carcinoma in situ: analysis of 127 cases. J Ultrasound Med. 2020;39(1):39–50.

    Article  PubMed  Google Scholar 

  18. Moon HJ, Kim EK, Kim MJ, Yoon JH, Park VY. Comparison of clinical and pathologic characteristics of ductal carcinoma in situ detected on mammography versus ultrasound only in asymptomatic patients. Ultrasound Med Biol. 2019;45(1):68–77.

    Article  PubMed  Google Scholar 

  19. Evans A, Vinnicombe S. Overdiagnosis in breast imaging. Breast. 2017;31:270–3.

    Article  PubMed  Google Scholar 

  20. Conti A, Duggento A, Indovina I, Guerrisi M, Toschi N. Radiomics in breast cancer classification and prediction. Semin Cancer Biol. 2020;1:323–4.

    Google Scholar 

  21. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology. 2016;278(2):563–77.

    Article  PubMed  Google Scholar 

  22. Pinker K, Chin J, Melsaether AN, Morris EA, Moy L. Precision medicine and radiogenomics in breast cancer: new approaches toward diagnosis and treatment. Radiology. 2018;287(3):732–47.

    Article  PubMed  Google Scholar 

  23. 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  CAS  Google Scholar 

  24. Lin F, Wang Z, Zhang K, Yang P, Ma H, Shi Y, Liu M, Wang Q, Cui J, Mao N, 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.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Cardoso F, Bartlett JMS, Slaets L, van Deurzen CHM, van Leeuwen-Stok E, Porter P, Linderholm B, Hedenfalk I, Schröder C, Martens J, et al. Characterization of male breast cancer: results of the EORTC 10085/TBCRC/BIG/NABCG International Male Breast Cancer Program. Ann Oncol. 2018;29(2):405–17.

    Article  CAS  PubMed  Google Scholar 

  26. Hu HT, Wang Z, Huang XW, Chen SL, Zheng X, Ruan SM, Xie XY, Lu MD, Yu J, Tian J, et al. Ultrasound-based radiomics score: a potential biomarker for the prediction of microvascular invasion in hepatocellular carcinoma. Eur Radiol. 2019;29(6):2890–901.

    Article  PubMed  Google Scholar 

  27. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, Beets-Tan RGH, Fillion-Robin JC, Pieper S, Aerts H. Computational radiomics system to decode the radiographic phenotype. Cancer Res. 2017;77(21):e104–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Peng Y, Lin P, Wu L, Wan D, Zhao Y, Liang L, Ma X, Qin H, Liu Y, Li X, et al. Ultrasound-based radiomics analysis for preoperatively predicting different histopathological subtypes of primary liver cancer. Front Oncol. 2020;10:1646.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Lucia F, Visvikis D, Vallières M, Desseroit MC, Miranda O, Robin P, Bonaffini PA, Alfieri J, Masson I, Mervoyer A, et al. External validation of a combined PET and MRI radiomics model for prediction of recurrence in cervical cancer patients treated with chemoradiotherapy. Eur J Nucl Med Mol Imaging. 2019;46(4):864–77.

    Article  PubMed  Google Scholar 

  30. Atabaki-Pasdar N, Ohlsson M, Viñuela A, Frau F, Pomares-Millan H, Haid M, Jones AG, Thomas EL, Koivula RW, Kurbasic A, et al. Predicting and elucidating the etiology of fatty liver disease: a machine learning modeling and validation study in the IMI DIRECT cohorts. PLoS Med. 2020;17(6):e1003149.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Kapwata T, Gebreslasie MT. Random forest variable selection in spatial malaria transmission modelling in Mpumalanga Province, South Africa. Geospat Health. 2016;11(3):434.

    Article  PubMed  Google Scholar 

  32. Tian XP, Su N, Wang L, Huang WJ, Liu YH, Zhang X, Huang HQ, Lin TY, Ma SY, Rao HL, et al. A CpG Methylation Classifier to Predict Relapse in Adults with T-Cell Lymphoblastic Lymphoma. Clin Cancer Res. 2020;26(14):3760–70.

    Article  CAS  PubMed  Google Scholar 

  33. Naeem SM, Mabrouk MS, Marzouk SY, Eldosoky MA. A diagnostic genomic signal processing (GSP)-based system for automatic feature analysis and detection of COVID-19. Brief Bioinform. 2020;22:1197–206.

    Article  Google Scholar 

  34. Dey N, Rajinikanth V, Fong SJ, Kaiser MS, Mahmud M. Social group optimization-assisted Kapur’s entropy and morphological segmentation for automated detection of COVID-19 infection from computed tomography images. Cognit Comput. 2020;54:1–13.

    Google Scholar 

  35. Rosner D, Bedwani RN, Vana J, Baker HW, Murphy GP. Noninvasive breast carcinoma: results of a national survey by the American College of Surgeons. Ann Surg. 1980;192(2):139–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Groen EJ, Elshof LE, Visser LL, Rutgers EJT, Winter-Warnars HAO, Lips EH, Wesseling J. Finding the balance between over- and under-treatment of ductal carcinoma in situ (DCIS). Breast. 2017;31:274–83.

    Article  PubMed  Google Scholar 

  37. Davey C, White V, Warne C, Kitchen P, Villanueva E, Erbas B. Understanding a ductal carcinoma in situ diagnosis: patient views and surgeon descriptions. Eur J Cancer Care. 2011;20(6):776–84.

    Article  CAS  Google Scholar 

  38. Hunter NB, Kilgore MR, Davidson NE. The long and winding road for breast cancer biomarkers to reach clinical utility. Clin Cancer Res. 2020;26:5543–5.

    Article  PubMed  PubMed Central  Google Scholar 

  39. van Seijen M, Lips EH, Thompson AM, Nik-Zainal S, Futreal A, Hwang ES, Verschuur E, Lane J, Jonkers J, Rea DW, et al. Ductal carcinoma in situ: to treat or not to treat, that is the question. Br J Cancer. 2019;121(4):285–92.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Visser LL, Elshof LE, Van de Vijver K, Groen EJ, Almekinders MM, Sanders J, Bierman C, Peters D, Hofland I, Broeks A, et al. Discordant marker expression between invasive breast carcinoma and corresponding synchronous and preceding DCIS. Am J Surg Pathol. 2019;43(11):1574–82.

    Article  PubMed  Google Scholar 

  41. Zhang GJ, Kimijima I, Abe R, Kanno M, Katagata N, Hara K, Watanabe T, Tsuchiya A. Correlation between the expression of apoptosis-related bcl-2 and p53 oncoproteins and the carcinogenesis and progression of breast carcinomas. Clin Cancer Res. 1997;3(12 Pt 1):2329–35.

    CAS  PubMed  Google Scholar 

  42. Davis JE, Nemesure B, Mehmood S, Nayi V, Burke S, Brzostek SR, Singh M. Her2 and Ki67 biomarkers predict recurrence of ductal carinoma in situ. Appl Immunohistochem Mol Morphol. 2016;24(1):20–5.

    Article  CAS  PubMed  Google Scholar 

  43. Cornfield DB, Palazzo JP, Schwartz GF, Goonewardene SA, Kovatich AJ, Chervoneva I, Hyslop T, Schwarting R. The prognostic significance of multiple morphologic features and biologic markers in ductal carcinoma in situ of the breast: a study of a large cohort of patients treated with surgery alone. Cancer. 2004;100(11):2317–27.

    Article  PubMed  Google Scholar 

  44. Lee CH, Dershaw DD, Kopans D, Evans P, Monsees B, Monticciolo D, Brenner RJ, Bassett L, Berg W, Feig S, et al. Breast cancer screening with imaging: recommendations from the Society of Breast Imaging and the ACR on the use of mammography, breast MRI, breast ultrasound, and other technologies for the detection of clinically occult breast cancer. J Am Coll Radiol. 2010;7(1):18–27.

    Article  PubMed  Google Scholar 

  45. Vourtsis A, Berg WA. Breast density implications and supplemental screening. Eur Radiol. 2019;29(4):1762–77.

    Article  PubMed  Google Scholar 

  46. Gartlehner G, Thaler K, Chapman A, Kaminski-Hartenthaler A, Berzaczy D, Van Noord MG, Helbich TH. Mammography in combination with breast ultrasonography versus mammography for breast cancer screening in women at average risk. Cochrane Database Syst Rev. 2013;2013(4):Cd009632.

    PubMed Central  Google Scholar 

  47. Pu H, Peng J, Xu F, Liu N, Wang F, Huang X, Jia Y. Ultrasound and clinical characteristics of false-negative results in mammography screening of dense breasts. Clin Breast Cancer. 2020;20(4):317–25.

    Article  PubMed  Google Scholar 

  48. Gunawardena DS, Burrows S, Taylor DB. Non-mass versus mass-like ultrasound patterns in ductal carcinoma in situ: is there an association with high-risk histology? Clin Radiol. 2020;75(2):140–7.

    Article  CAS  PubMed  Google Scholar 

  49. Partridge SC, Nissan N, Rahbar H, Kitsch AE, Sigmund EE. Diffusion-weighted breast MRI: clinical applications and emerging techniques. J Magn Reson Imaging. 2017;45(2):337–55.

    Article  PubMed  Google Scholar 

  50. Alderson PO, Summers RM. The evolving status of radiomics. J Natl Cancer Inst. 2020;112:869–70.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Lundervold AS, Lundervold A. An overview of deep learning in medical imaging focusing on MRI. Z Med Phys. 2019;29(2):102–27.

    Article  PubMed  Google Scholar 

  52. Demircioglu A, Grueneisen J, Ingenwerth M, Hoffmann O, Pinker-Domenig K, Morris E, Haubold J, Forsting M, Nensa F, Umutlu L. A rapid volume of interest-based approach of radiomics analysis of breast MRI for tumor decoding and phenotyping of breast cancer. PLoS ONE. 2020;15(6):e0234871.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Zhou J, Tan H, Li W, Liu Z, Wu Y, Bai Y, Fu F, Jia X, Feng A, Liu H, et al. Radiomics signatures based on multiparametric MRI for the preoperative prediction of the her2 status of patients with breast cancer. Acad Radiol. 2020;23:568–98.

    Google Scholar 

  54. Prasanna P, Tiwari P, Madabhushi A. Co-occurrence of local anisotropic gradient orientations (CoLlAGe): a new radiomics descriptor. Sci Rep. 2016;6:37241.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: HQ and YL; Methodology: PL, DW and XL; Formal analysis and investigation: LW, YZ, PL, HQ, YL, DW, XL, YH, HY; Writing—original draft preparation: LW and YZ; Writing—review and editing: LW and YZ; Resources: HY and YH; Supervision: HY and YH. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yun He or Hong Yang.

Ethics declarations

Ethics approval and consent to participate

This retrospective breast DCIS study was approved by the ethics committee of the First Affiliated Hospital of Guangxi Medical University. Informed consent was waived. This study on the implementation of all procedures are in line with the National Research Council of moral standards.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have 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

. 5234 radiomics features matrix file.

Additional file 2

. Modeling matrix file for molecular biomarkers.

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

Wu, L., Zhao, Y., Lin, P. et al. Preoperative ultrasound radiomics analysis for expression of multiple molecular biomarkers in mass type of breast ductal carcinoma in situ. BMC Med Imaging 21, 84 (2021). https://doi.org/10.1186/s12880-021-00610-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12880-021-00610-7

Keywords