Skip to main content

Development and validation of a radiopathomic model for predicting pathologic complete response to neoadjuvant chemotherapy in breast cancer patients



Neoadjuvant chemotherapy (NAC) has become the standard therapeutic option for early high-risk and locally advanced breast cancer. However, response rates to NAC vary between patients, causing delays in treatment and affecting the prognosis for patients who do not sensitive to NAC.

Materials and methods

In total, 211 breast cancer patients who completed NAC (training set: 155, validation set: 56) were retrospectively enrolled. we developed a deep learning radiopathomics model(DLRPM) by Support Vector Machine (SVM) method based on clinicopathological features, radiomics features, and pathomics features. Furthermore, we comprehensively validated the DLRPM and compared it with three single-scale signatures.


DLRPM had favourable performance for the prediction of pathological complete response (pCR) in the training set (AUC 0.933[95% CI 0.895–0.971]), and in the validation set (AUC 0.927 [95% CI 0.858–0.996]). In the validation set, DLRPM also significantly outperformed the radiomics signature (AUC 0.821[0.700–0.942]), pathomics signature (AUC 0.766[0.629–0.903]), and deep learning pathomics signature (AUC 0.804[0.683–0.925]) (all p < 0.05). The calibration curves and decision curve analysis also indicated the clinical effectiveness of the DLRPM.


DLRPM can help clinicians accurately predict the efficacy of NAC before treatment, highlighting the potential of artificial intelligence to improve the personalized treatment of breast cancer patients.

Peer Review reports


Breast cancer is the most common cancer among women in the world, and it related incidence rate is continuously rising [1, 2]. Neoadjuvant chemotherapy (NAC) has become the standard therapeutic option for early high-risk and locally advanced breast cancer [3]. When breast cancer patients have a pathological complete response (pCR) to NAC, it can help patients lower the stage and shrink the tumor to receive more conservative treatment, and its event free survival (EFS) and overall survival (OS) are significantly improved [4, 5]. However, because of the heterogeneity and complexity of tumors, not all patients benefit from NAC. For patients who are not sensitive to treatment, although disease progression rarely occurs during NAC [6], the long-term treatment process will still have side effects [7, 8], which may also lead to missing the best time to change the treatment plan. Currently, there is an urgent requirement for accurate prediction of the response before the NAC, which is critical for breast cancer patients who are destined to have no response.

Radiomics could predict effectively pCR in patients with breast cancer. In the process of NAC, some studies used the images of different treatment nodes to fuse each other to predict pCR [9, 10], but the image data required for model construction was obtained after NAC, and the clinical practicability of the model was poor. Some scholars have made optimization on this basis [11, 12], it could benefit patients in some extent, but the radiomics features only provide tumor information from a macroscopic perspective in vitro, it cannot completely reflect information about pCR. As another source of medical images, histopathology combined with machine learning can help in risk stratification, prognosis prediction and adjuvant chemotherapy efficacy prediction [13,14,15,16]. Pathomics differ from radiomics in that they provide microstructural information about the tumor microenvironment, which can complement tumor heterogeneity and enhance the predictive power of existing models. So we hypothesized that a multi-scale model integrating the features of radiomics and pathomics could efficiently predict pCR.

In this study, we aimed to develop and validate a deep learning radiopathomics model(DLRPM) for the prediction of pCR to NAC in patients with breast cancer using Contrast-Enhanced Computed Tomography (CECT) images and whole slide images (WSIs). This proposed DLRPM can be used for early adjustment therapy in non-PCR patients to improve pCR rates and avoid toxic side effects. This might provide clinicians with treatment strategies to improve the effectiveness of individual therapy.

Materials and methods


This retrospective study was approved by the Institutional Review Board of the Southwestern Medical University Hospital (No. KY2022216), and the requirement for written informed consent was waived. This study collected 1532 patients with breast cancer who underwent a CECT examination between January 2020 and March 2022 from the Picture Archiving and Communications System (PACS).

The inclusion criteria were as follows: (a) Pathological biopsy confirmed non-specific invasive breast cancer with no distant metastasis; (b) The patient has undergone 6–8 cycles of NAC; (c) surgery was performed after NAC; (d) The available clinical data. A total of 245 patients fulfilling the inclusion criteria were enrolled.

Exclusion criteria were as follows: (a)No histopathological evaluation results. (b) Lack of images in venous phase or poor image quality. (c)Synchronous tumors or history of other malignancy.

A total of 211 breast cancer patients with non-specific invasive were enrolled in the study from February 2020 and March 2022, we divided patients into training set and independent validation set in chronological order. Patients who performed their first NAC treatment before September 2021 were used as the training set, the remaining patients were used as the validation set, the ratio of training set to validation set was about 7:3. A flowchart of the patients’ collection is shown in Fig. 1.

Fig. 1
figure 1

Flow diagram of patient cohort selection

Workflow of study

The workflow of this study is shown in Fig. 2, including (1) image acquisition, (2) feature extraction, (3) feature selection, (4) model construction, and (5) model validation.

Fig. 2
figure 2

Workflow of Study. The images were preprocessed for feature extraction. After feature evaluation and model construction, four sets of features [radiomic signature (RS), pathologic signature (PS), deep learning pathologic signature (DLPS) and clinical features] were generated and further used to construct DLRPM. The performance of DLRPM in predicting pCR before NAC was validated in validation set

The CECT images and WSIs acquisition

All patients received contrast-enhanced CT chest examination (Netherlands, Philips Medical Systems) before NAC treatment. The scanning procedure was as follows: The contrast agent (iodohexol, 320 mg/mL) was injected into the median cubital vein with a double-barrel high-pressure syringe (dose 1.0 mL/kg, flow rate 3.0 mL/s). The CT value of blood vessels at the level of the aortic arch was monitored after injection of contrast agent. The Enhanced CT scans are automatically triggered when the CT value reaches around 250 HU. And venous phase scans were performed after a delay of 30s.

The pathologists collected all biopsy samples of breast cancer patients using crude needle puncture before NAC. Firstly, biopsy tissue was soaked in 10% formalin for 4 h and buried in immunohistochemical paraffin wax. Subsequently, the biopsied tissue was severed at 4-μm intervals and stained with hematoxylin and eosin (H&E) for pathological evaluation. A pathologist with 8 years of experience in pathological diagnosis scanned all H&E-stained histopathological slides using a digital slide scanner (KFBio KF-PRO-020) at 10 × magnification to obtain WSIs of the breast cancer patients, and images were digitized as kbf. format files, which were managed with the KF-Viewer software (version

Pathological complete response assessment

In accordance with the National Comprehensive Cancer Network (NCCN) guideline [17], all patients received six or eight cycles of NAC. The NAC regimens were based on taxane or taxane and anthracycline; all human epidermal growth factor receptor 2 (HER2) positive patients also received trastuzumab. At the end of treatment, we performed an initial imaging assessment of the efficacy of NAC according to the Response Evaluation Criteria in Solid Tumors (RECIST) 1.1 [18]. Subsequently, the final pCR status of each patient was determined by the pathological findings after surgery (Fig. 3). pCR was defined as the complete absence of invasive tumor cells in the breast and axillary lymph nodes, regardless of the presence of residual ductal carcinoma in situ (ypT0/isypN0) (Fig. 3A).

Fig. 3
figure 3

CECT and histology images from complete responder (A) and partial responder (B) before NAC and after 8 courses of NAC. In the CECT image, it is seen that the tumor in the complete responder have completely dissipated in the post-NAC image, and stromal tissue with no visible tumor cells was presented in the pathological images. But CECT and histology images from and partial responder shows residual tumor cells but reduced compared to baseline

Radiomics feature extraction

The volume of interest (VOI) segmentation was performed using 3D-Slicer software. All manual segmentation of the CECT images were performed by 2 practicing experienced radiologists. Both radiologists were blinded to the patient's clinical data when they evaluated the CECT images. Firstly, the VOIs covering the whole tumor (VOI 1) were segmented. After manual tumor segmentation, we automatically segment the peritumoral regions (VOI 2) (Figure S1). The regions (2-mm radius) surrounding the tumor were defined as the peritumoral regions. If the peritumor regions were beyond the parenchyma of the breast after the spread, the portion beyond the parenchyma was removed manually.

According to the instructions of the Image Biomarker Standardization Initiative [19], radiomics features were extracted from the VOI 1 and VOI 2 using PyRadiomics, before feature extraction, PyRadiomics was also used for image preprocessing. A total of 3814 radiomics features were extracted from two VOIs per patient. These features included First Order, Shape-based (2D and 3D), Gray Level Cooccurence Matrix (GLCM), Gray Level Run Length Matrix (GLRLM), Gray Level Size Zone Matrix (GLSZM), Gray Level Dependece Matrix (GLDM); and Filters included Laplacian of Gaussian (LoG), Wavelet, Square, Square Root, Logarithm, Exponential, Gradient, Local Binary Pattern (LBP) 3D. To make features reproducible, an interclass correlation coefficient (ICC) higher than 0.75 was considered credible [20].

Pathomics feature extraction

In KF Viewer, WSIs were magnified 10 × , the pathologist (W.Q) selected the sample area containing nuclear pleomorphism, mitosis, carcinoma infiltration, cancer invasion, tumor cell differentiation and pathological grading, and obtained five typical non-overlapping screenshots with a field of vision of 1534 × 832 pixels, and then confirmed by the other pathologist (J.M).They had 3 years and 8 years of experience in breast cancer pathological diagnosis respectively. We saved the selected screenshots as format files (. jpg, 300dpi). If the two pathologists have different opinions, they will consult the third pathologist to make a decision. All screenshots were cut into small frame tiles (512 × 512 pixels) by sampling without overlap for subsequent analysis (Figure S2).

We used CellProfiler (version 4.0.7) [21], an open-source image analysis software developed by Broad Institute (Cambridge, Massachusetts), to extract quantitative pathomics features of selected pathological screenshots. Based on the “Unmix Colors” module to separate H&E-stained images and convert them into haematoxylin-stained and eosin-stained greyscale images, The H&E-stained images were also converted to greyscale images using the “ColorToGray” module (Figure S3). We measured images twice, in the first measurement, we obtained 136 original features, which summarize the three types of images in general. For the second measurement, we made a careful exploration of hematoxylin images. First, we identified the primary and secondary objects, and then measured them. After measurement, we took their mean, median and standard deviation as our research characteristics and 1054 pathomics features were obtained, the extracted features were aggregated by mean of the values for every 10 tiles in each WSI. Detailed method of feature extraction in Figure S4.

Resnet50 was employed to extract deep learning features of pathomics. Before extracting features, All the small tiles went through color normalization with the Vahadane method based on staintools (Figure S2), which was an open-source package based on python for stain normalization and augmentation. The input area was 512 × 512 pixels, and the transfer learning took the pretrained weights of Resnet50 on the ImageNet dataset as the initial weights of the model, The model was fine-tuned using data from our data. Resnet50 was adjusted from the original multi-classification task to a binary classification task, we extracted the deep learning features from the last layer of resnet50, and the principal component analysis (PCA) algorithm further compressed the deep learning features, the extracted features were aggregated by mean of the values for every 10 tiles in each WSI. We obtained a total of 200 deep learning features.

Feature selection and signature construction

The radiomics, pathomics features and pathomics deep learning features can reveal tumour information from macroenvironment and microenvironment perspectives, repectively. However, these features were high-dimensional data, which had an adverse impact on predicting the pCR to NAC. Therefore, we should obtain the features which were most closely related to pCR in the training set. Firstly, all variables were normalized, and a U-test was performed on each feature as a preliminary selection to remove redundant features. To sufficiently extract discriminative features in this process, the threshold of p value was determined with 0.05. Subsequently, considering the dependence between features, we perform correlation analysis on the features, if the correlation coefficient between the two features was greater than 0.9, one of them was excluded. Then, the least absolute shrinkage and selection operator (LASSO) algorithm was utilized to select the extracted features [22], and tenfold cross-validation was used to select the value of Lambda to determine the optimal features.

Based on the above three types of optimal features, we constructed three distinct single-scale prediction models by Support Vector Machine (SVM) method [23], The best regularization parameter C and Gamma (γ) for Gaussian Radial Basis Function (RBF) kernel were determined by fivefold cross validation and grid search. Then, each model prediction value was used to construct signatures, named radiomics signature (RS), pathomics signature (PS), and deep learning pathomics signature (DLPS), respectively.

DLRPM development and validation

Independent clinical predictor and three single-scale signatures were used to construct the DLRPM for the integrated prediction of pCR in breast cancer patients by a similar non-linear SVM method. We used the following methods to comprehensively evaluate the model. Receiver operating characteristic (ROC) curve analysis, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were employed to evaluate the discrimination performance [24]. The calibration of models were assessed using the calibration curve and the Hosmer–Lemeshow test was used to assess the goodness-of-fit of the models.The Decision curve analysis (DCA) of all models was performed to quantify the net benefit of patients under different threshold probabilities in the sets to assess the clinical value of the predictive models in our study [25].

In addition, the net reclassification improvement (NRI) test and integrated discrimination improvement (IDI) test were calculated to compare the performance of DLRPM and single-scale signatures.

Statistical analysis

All statistical analysis was performed using R studio (version 4.1.1; R Studio, and Jupyter Notebook (version 6.4.11). Differences of categorical variables were calculated with the chi-square test or Fisher’s exact test. The differences of continuous variables were analyzed using independent t-test or Mann–Whitney U test. All tests were two-sided, and two-tailed p < 0.05 were considered statistically significant.


Demographic and clinicopathological characteristics

The baseline characteristics of all patients are summarized in Table 1. A total of 211 patients with breast cancer were enrolled in this study, patients with pCR accounted for 39.35% (61/155) and 37.50% (21/56) of the training and validation sets, respectively, the findings on HER2 was significantly correlated with the pCR status in both the training and validation set (p < 0.05). In addition, there were no statistically significant differences between the two sets (p > 0.05) (Table S1).

Table 1 Clinical characteristics of patients

Radio-pathomics feature selection and signature construction

The inter-observer reproducibility of the feature extraction was excellent, with inter-observer ICCs ranging from 0.758 to 0.953 for CECT. U-test and Spearman correlation coefficient analysis was performed exclude redundant features, which resulted in 154 radiomic features, 21 pathomics Feature and 8 deep learning Pathomics Feature per patient.

Then, LASSO was adopted to deeply select the pre-existing features. For LASSO, the subset of the best features depended on the choice of lambda value, and we used fivefold cross validation to find the best lambda value. To further simplify the model, step-forward feature selection was then conducted to reduce the optimal features, we pick the lambda value with one standard error to select features, which resulted in 30 radiomic features and 12 pathomics Feature, for deep learning pathomics Feature, Lasso did not further reduce features (Fig. 4) (Table S2, 3).

Fig. 4
figure 4

Feature selection process. Radiomics features (A, B) and pathomics features (C, D) were selected by the LASSO model with tuning parameter (λ) using fivefold cross-validation via minimum and 1se criteria

We constructed three predictive signatures by non-linear SVM method respectively. In training set, the grid search with fivefold cross validation found the optimal parameters of the three models, RS (C = 4.67, gamma = 0.0015), PS (C = 1, gamma = 0.0028) and DLPS (C = 2.15, gamma = 0.031). Raincloud plot (Fig. 5) visualized the different distributions of the samples in training and validation sets, it indicated Three single-scale signatures already had a certain discriminant ability (p < 0·05, Table 2).

Fig. 5
figure 5

The raincloud plot visualizes prediction probability of RS, PS and DLPS, it shows the sample distribution locations and interval sample densities for the training (A)and validation set(B) of signatures

Table 2 Predicted probability of signatures and in training and validation set

DLRPM development and validation

By integrating HER2, RS, PS and DLPS in the training set, we developed the DLRPM comprehensive prediction model by using the nonlinear support vector machine method. The same method was used to find the optimal parameter of DLRPM (C = 5, gamma = 0.01). ROC curves were used to assess discrimination performance, DLRPM accurately predicted pCR in Training set (AUC 0·933[95% CI 0.895–0.971]) and validation set (0.927 [95% CI 0.858–0.996]). The sensitivity of DLRPM was markedly high in validation set (94.28%), whereas the specificity remained moderate (76.19%). The NPV of DLRPM exceeded 90% in validation set, whereas the PPV was around 71.25% (Table 3, Fig. 6A).

Table 3 Discrimination performance of predict models for predicting pCR status in breast cancer patients
Fig. 6
figure 6

ROC analysis of predict models for predicting pCR in the training set (A) and validation set (B), respectively. C Calibration curves of models in training set on discriminating Non-pCR versus pCR. D Decision curve analysis in training set using RS, PS, DLPS and DLRPM

In validation set, RS and DLPS yielded marginally AUC values of 0.821(0.700–0.942) and 0.804(0.683–0.925) (Fig. 6B), whereas PS had a lower AUC value of 0.766(0.629–0.903). the DeLong test, NRI and IDI showed that performance of the three single scale models had no significant difference(all p > 0.05).

Compared with single-scale prediction models, DLRPM showed superior to evaluate the discrimination performance. The improvements in discriminative ability were confirmed by the NRI tests (all p < 0·05) and IDI (all p < 0·05). Calibration plots demonstrated good agreement between all model prediction and the actual observation for detecting pCR (Fig. 6C). The Hosmer–Lemeshow test showed non-significant statistics in both groups (p > 0.05). The DCA plots showed that the DLRPM provided better net benefit compared with single-scale prediction models, it indicated DLRPM had a better clinical benefit (Fig. 6D).


In this study, we developed a multi-scale integrated model for prediction of pCR to NAC in breast cancer patients before treatment based on SVM algorithm, combining CECT images with WISs. In the independent validation set, DLRPM had a better performance in terms of discriminative ability, calibration, and clinical utility relative to other single-scale models. In clinical practice, patients who predicted for pCR using DLRPM should be given aggressive NAC therapy, and intensive follow-up strategies were used to improve survival and quality of life. DLRPM assisted clinicians predict accurately the efficacy of NAC treatment before it is administered, which is critical for developing patient treatment plans and optimizing overall patient management.

Predictive biomarkers of response to NAC for breast cancer patients have been the focus topic of research [26], but not all biomarkers were applicable to clinical practice. Previous studies have explored the predictive effectiveness of genetic biomarkers, but they have not been applied in clinical practice because they are costly and and not widespread [27, 28]. Secondly, we should consider the predictive power of the model. Clinicopathological factors were currently used to estimate the potential benefit of NAC, however, satisfactory performance cannot be achieved based on clinical characteristics alone [29,30,31]. Full digitalization of the stained tissue sections has become feasible because of advances in slide scanning technology and reductions in the cost of digital storage. In our study, DLRPM not only had good predictive power (AUCs > 0.9), but DLRPM had a stable source of modeling and validation data in clinical practice, providing assurance for subsequent large-scale studies.

Previous studies had demonstrated that radiomic biomarker can predict pCR to NAC in breast cancer patients. Many studies applied pre-treatment radiomic features combined with clinicopathological factors for efficacy prediction, but their predictive efficacy was unstable [32,33,34]; others predicted pCR using pretreatment and post-treatment US images, and their model had better predictive efficacy in the external validation set [9], but promising results were mainly attributed to post-treatment US images, which provided direct information on tumor regression, and the model was unable to provide early estimates of treatment response to guide the implementation of NAC because the images required for the construction of model were obtained late and clinical utility was poor. In contrast, the US images from the second and fourth courses of treatment were matched with US images before treatment, respectively, and This staged prediction pipeline benefited patients to a certain extent [11], and treatment regimens were adjusted based on the prediction results. However, for patients with breast cancer, accurate prediction of neoadjuvant chemotherapy efficacy before treatment can help maximise patient’s benefit.

In the pathological examination, the pathologist used light microscope to determine the benign and malignant tumor, the growth mode and differentiation degree of tumor cells under the light microscope. Because of the limitation of microscope magnification, pathologists were not able to describe the microscopic information of each slide in detail. In recent years,with the application of cell analysis software and the development of deep learning algorithms [14, 15], several scholars researchers have extracted image features from digital pathology slices for quantitative analysis. Breast cancer was a highly heterogeneous tumor, and tumor microenvironment will change when tumor responds to NAC [35]. This subtle change was not detectable by the naked eye. Pathomics can capture the microstructure of tumors and provide the characteristics of cells and microenvironment in tumor lesions. Previous studies had shown pathomics features were used to predict the efficacy and prognosis of adjuvant chemotherapy for gastric cancer [16], and have achieved good performance. In our study, pathomics features also showed a stable predictive ability in predicting the pCR to NAC.

DLRPM integrates macroscopic radiomics and microscopic pathomics features for integrated prediction of pCR to NAC in breast cancer patients. In the radiomics workflow, we added a 2 mm peritumoral region for feature extraction, and the great predictive potential of the peritumoral region features had been demonstrated in the Ning’s study [36]. In our study, the peritumoral region features accounted for 50% of the total radiomics features selected by LASSO and contributed significantly to the predictive power. In the pathomics workflow, we learned from the experience of the current study in CellProfile extraction and performed two extractions of our pathology images. In our analysis, we found that the predictive power of the two types of features can complement each other. DLRPM has high sensitivity and NPV in the independent validation set, indicating that the model can reliably identify individuals without pathological complete response. This study is consistent with the findings of Li’s study [37]. and they could avoid subsequent ineffective treatments, breast cancer patients who were destined not to respond will benefit from the predicted results of DLRPM.

Although our research is innovative, there were also several limitations in our current study. Firstly, this study is a retrospective study, and all of breast cancer patients were obtained from a single medical institution. Considering the limited number of study samples, we will further obtain a large number of sample data from multiple medical institutions and perform prospective studies to validate the generalization and accuracy of DLRPM. Secondly, All the features of radiomics are derived from the venous phase images of CECT and lack of diversity. In the future, multiphase CT data can be collected for feature enrichment.Finally,VOI segmentation of tumor is not automatic, and the probability of error in artificial semi-automatic segmentation is large and difficult to find. This may be overcome by automatic segmentation artificial intelligence system in the future.

In conclusion, we established DLRPM based on the characteristics of radiomics and pathomics to predict the complete pathological response of breast cancer patients to neoadjuvant chemotherapy. This model can help clinicians accurately predict the efficacy of neoadjuvant chemotherapy before treatment, highlighting the potential of artificial intelligence to improve the personalized treatment of breast cancer patients.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.


  1. Giaquinto AN, Sung H, Miller KD, Kramer JL, Newman LA, Minihan A, Jemal A, Siegel RL. Breast cancer statistics, 2022. CA Cancer J Clin. 2022;72(6):524–41.

    Article  PubMed  Google Scholar 

  2. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2022. CA Cancer J Clin. 2022;72(1):7–33.

    Article  PubMed  Google Scholar 

  3. Loi S. The ESMO clinical practise guidelines for early breast cancer: diagnosis, treatment and follow-up: on the winding road to personalized medicine. Ann Oncol. 2019;30(8):1183–4.

    Article  CAS  PubMed  Google Scholar 

  4. Guarneri V, Griguolo G, Miglietta F, Conte PF, Dieci MV, Girardi F. Survival after neoadjuvant therapy with trastuzumab-lapatinib and chemotherapy in patients with HER2-positive early breast cancer: a meta-analysis of randomized trials. ESMO Open. 2022;7(2): 100433.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Cortazar P, Zhang L, Untch M, Mehta K, Costantino JP, Wolmark N, Bonnefoi H, Cameron D, Gianni L, Valagussa P, et al. Pathological complete response and long-term clinical benefit in breast cancer: the CTNeoBC pooled analysis. Lancet. 2014;384(9938):164–72.

    Article  PubMed  Google Scholar 

  6. Wang H, Mao X. Evaluation of the Efficacy of Neoadjuvant Chemotherapy for Breast Cancer. Drug Des Devel Ther. 2020;14:2423–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Geyer CE, Sikov WM, Huober J, Rugo HS, Wolmark N, O’Shaughnessy J, Maag D, Untch M, Golshan M, Lorenzo JP, et al. Long-term efficacy and safety of addition of carboplatin with or without veliparib to standard neoadjuvant chemotherapy in triple-negative breast cancer: 4-year follow-up data from BrighTNess, a randomized phase III trial. Ann Oncol. 2022;33(4):384–94.

    Article  CAS  PubMed  Google Scholar 

  8. Jang MK, Park S, Park C, Doorenbos AZ, Go J, Kim S. Body composition change during neoadjuvant chemotherapy for breast cancer. Front Oncol. 2022;12: 941496.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Jiang M, Li CL, Luo XM, Chuan ZR, Lv WZ, Li X, Cui XW, Dietrich CF. Ultrasound-based deep learning radiomics in the assessment of pathological complete response to neoadjuvant chemotherapy in locally advanced breast cancer. Eur J Cancer. 2021;147:95–105.

    Article  CAS  PubMed  Google Scholar 

  10. Wu L, Ye W, Liu Y, Chen D, Wang Y, Cui Y, Li Z, Li P, Li Z, Liu Z, et al. An integrated deep learning model for the prediction of pathological complete response to neoadjuvant chemotherapy with serial ultrasonography in breast cancer patients: a multicentre, retrospective study. Breast Cancer Res. 2022;24(1):81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Gu J, Tong T, He C, Xu M, Yang X, Tian J, Jiang T, Wang K. Deep learning radiomics of ultrasonography can predict response to neoadjuvant chemotherapy in breast cancer at an early stage of treatment: a prospective study. Eur Radiol. 2022;32(3):2099–109.

    Article  PubMed  Google Scholar 

  12. Liu Y, Wang Y, Wang Y, Xie Y, Cui Y, Feng S, Yao M, Qiu B, Shen W, Chen D, et al. Early prediction of treatment response to neoadjuvant chemotherapy based on longitudinal ultrasound images of HER2-positive breast cancer patients by Siamese multi-task network: A multicentre, retrospective cohort study. EClinicalMedicine. 2022;52: 101562.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Cao R, Yang F, Ma SC, Liu L, Zhao Y, Li Y, Wu DH, Wang T, Lu WJ, Cai WJ, et al. Development and interpretation of a pathomics-based model for the prediction of microsatellite instability in Colorectal Cancer. Theranostics. 2020;10(24):11080–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Chen S, Jiang L, Zheng X, Shao J, Wang T, Zhang E, Gao F, Wang X, Zheng J. Clinical use of machine learning-based pathomics signature for diagnosis and survival prediction of bladder cancer. Cancer Sci. 2021;112(7):2905–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Chen S, Jiang L, Gao F, Zhang E, Wang T, Zhang N, Wang X, Zheng J. Machine learning-based pathomics signature could act as a novel prognostic marker for patients with clear cell renal cell carcinoma. Br J Cancer. 2022;126(5):771–7.

    Article  CAS  PubMed  Google Scholar 

  16. Chen D, Fu M, Chi L, Lin L, Cheng J, Xue W, Long C, Jiang W, Dong X, Sui J, et al. Prognostic and predictive value of a pathomics signature in gastric cancer. Nat Commun. 2022;13(1):6903.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Goetz MP, Gradishar WJ, Anderson BO, Abraham J, Aft R, Allison KH, Blair SL, Burstein HJ, Dang C, Elias AD, et al. NCCN Guidelines Insights: Breast Cancer, Version 3.2018. J Natl Compr Canc Netw. 2019;17(2):118–26.

    Article  CAS  PubMed  Google Scholar 

  18. Eisenhauer EA, Therasse P, Bogaerts J, Schwartz LH, Sargent D, Ford R, Dancey J, Arbuck S, Gwyther S, Mooney M, et al. New response evaluation criteria in solid tumours: revised RECIST guideline (version 1.1). Eur J Cancer. 2009;45(2):228–47.

    Article  CAS  PubMed  Google Scholar 

  19. Zwanenburg A, Vallières M, Abdalah MA, Aerts H, Andrearczyk V, Apte A, Ashrafinia S, Bakas S, Beukinga RJ, Boellaard R, et al. The image biomarker standardization initiative: standardized quantitative radiomics for high-throughput image-based phenotyping. Radiology. 2020;295(2):328–38.

    Article  PubMed  Google Scholar 

  20. Koo TK, Li MY. A guideline of selecting and reporting intraclass correlation coefficients for reliability research. J Chiropr Med. 2016;15(2):155–63.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Carpenter AE, Jones TR, Lamprecht MR, Clarke C, Kang IH, Friman O, Guertin DA, Chang JH, Lindquist RA, Moffat J, et al. Cell Profiler: image analysis software for identifying and quantifying cell phenotypes. Genome Biol. 2006;7(10):R100.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Hu JY, Wang Y, Tong XM, Yang T. When to consider logistic LASSO regression in multivariate analysis? Eur J Surg Oncol. 2021;47(8):2206.

    Article  PubMed  Google Scholar 

  23. Huang MW, Chen CW, Lin WC, Ke SW, Tsai CF. SVM and SVM ensembles in breast cancer prediction. PLoS One. 2017;12(1): e0161501.

    Article  PubMed  PubMed Central  Google Scholar 

  24. de Hond AAH, Steyerberg EW, van Calster B. Interpreting area under the receiver operating characteristic curve. Lancet Digit Health. 2022;4(12):e853–5.

    Article  PubMed  Google Scholar 

  25. Vickers AJ, Cronin AM, Elkin EB, Gonen M. Extensions to decision curve analysis, a novel method for evaluating diagnostic tests, prediction models and molecular markers. BMC Med Inform Decis Mak. 2008;8:53.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Derouane F, van Marcke C, Berlière M, Gerday A, Fellah L, Leconte I, et al. Predictive biomarkers of response to neoadjuvant chemotherapy in breast cancer: current and future perspectives for precision medicine. Cancers (Basel). 2022;14(16):3876.

  27. Gong C, Cheng Z, Yang Y, Shen J, Zhu Y, Ling L, Lin W, Yu Z, Li Z, Tan W, et al. A 10-miRNA risk score-based prediction model for pathological complete response to neoadjuvant chemotherapy in hormone receptor-positive breast cancer. Sci China Life Sci. 2022;65(11):2205–17.

    Article  CAS  PubMed  Google Scholar 

  28. Chen L, Huang S, Liu Q, Kong X, Su Z, Zhu M, Fang Y, Zhang L, Li X, Wang J. PD-L1 protein expression is associated with good clinical outcomes and nomogram for prediction of disease free survival and overall survival in breast cancer patients received neoadjuvant chemotherapy. Front Immunol. 2022;13: 849468.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Ouldamer L, Bendifallah S, Pilloy J, Arbion F, Body G, Brisson C, Lavoué V, Lévêque J, Daraï E. Risk scoring system for predicting breast conservation after neoadjuvant chemotherapy. Breast J. 2019;25(4):696–701.

    Article  CAS  PubMed  Google Scholar 

  30. Arici S, Sengiz Erhan S, Geredeli C, Cekin R, Sakin A, Cihan S. The clinical importance of androgen receptor status in response to neoadjuvant chemotherapy in turkish patients with local and locally advanced breast cancer. Oncol Res Treat. 2020;43(9):435–40.

    Article  CAS  PubMed  Google Scholar 

  31. Haque W, Verma V, Hatch S, Suzanne Klimberg V, Brian Butler E, Teh BS. Response rates and pathologic complete response by breast cancer molecular subtype following neoadjuvant chemotherapy. Breast Cancer Res Treat. 2018;170(3):559–67.

    Article  CAS  PubMed  Google Scholar 

  32. Pesapane F, Rotili A, Botta F, Raimondi S, Bianchini L, Corso F, et al. Radiomics of MRI for the prediction of the pathological response to neoadjuvant chemotherapy in breast cancer patients: a single referral centre analysis. Cancers (Basel). 2021;13(17):4271.

  33. Guo L, Du S, Gao S, Zhao R, Huang G, Jin F, et al. Delta-radiomics based on dynamic contrast-enhanced mri predicts pathologic complete response in breast cancer patients treated with neoadjuvant chemotherapy. Cancers (Basel). 2022;14(14):3515.

  34. Wang Z, Lin F, Ma H, Shi Y, Dong J, Yang P, Zhang K, Guo N, Zhang R, Cui J, et al. Contrast-enhanced spectral mammography-based radiomics nomogram for the prediction of neoadjuvant chemotherapy-insensitive breast cancers. Front Oncol. 2021;11: 605230.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Urueña C, Lasso P, Bernal-Estevez D, Rubio D, Salazar AJ, Olaya M, Barreto A, Tawil M, Torregrosa L, Fiorentino S. The breast cancer immune microenvironment is modified by neoadjuvant chemotherapy. Sci Rep. 2022;12(1):7981.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Mao N, Shi Y, Lian C, Wang Z, Zhang K, Xie H, Zhang H, Chen Q, Cheng G, Xu C, et al. Intratumoral and peritumoral radiomics for preoperative prediction of neoadjuvant chemotherapy effect in breast cancer based on contrast-enhanced spectral mammography. Eur Radiol. 2022;32(5):3207–19.

    Article  CAS  PubMed  Google Scholar 

  37. Feng L, Liu Z, Li C, Li Z, Lou X, Shao L, Wang Y, Huang Y, Chen H, Pang X, et al. Development and validation of a radiopathomics model to predict pathological complete response to neoadjuvant chemoradiotherapy in locally advanced rectal cancer: a multicentre observational study. Lancet Digit Health. 2022;4(1):e8–17.

    Article  CAS  PubMed  Google Scholar 

Download references


Not applicable.


This research was supported by Sichuan Science and Technology Program (Grant No. 2022YFS0616, 2021YFQ0002).

Author information

Authors and Affiliations



Jieqiu Zhang: Visualization,Writing- Original draft preparation, Software. Qi Wu: Conceptualization, Methodology. Wei Yin: Investigation, Data curation. Lu Yang: Supervision. Bo Xiao: Software, Validation. Jianmei Wang: Writing- Reviewing and Editing, Data curation. Xiaopeng Yao: Writing- Reviewing and Editing. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Jianmei Wang or Xiaopeng Yao.

Ethics declarations

Ethics approval and consent to participate

This retrospective study was approved by the ethics committee Review Board of the Affiliated Hospital of Southwest Medical University (protocol code KY2022216, 20 June 2022) and in conformity to the Declaration of Helsinki.The need for informed consent was waived by the ethics committee Review Board of the Affiliated Hospital of Southwest Medical University (protocol code KY2022216), because of the retrospective nature of the study.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest.

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:

 Figure S1. VOI segmentation. (A) The process of lesion delineation in a patient with breast cancer who achieved pCR in the evaluation of the efficacy of NAC. (B) The process of lesion delineation in a patient with breast cancer who achieved Non-pCR in the evaluation of the efficacy of NAC. Figure S2. (A)Segmentation Process of WSI, (B)color normalization. Figure S3. Image preprocessing in CellProBased on the “Unmix Colors” module to separate H&E-stained images and convert them into haematoxylin-stained and eosin-stained greyscale images, The H&E-stained images were also converted to greyscale images using the “ColorToGray” module. Figure S4. (A)Pipeline 1. First, greyscale H&E, haematoxylin and eosin images were assessed by using the “MeasureImageQuality” module with three types of features, including blur features, intensity features and threshold features. Subsequently, ‘MeasureColocalization’ module measured the colocalization and correlation between intensities in haematoxylin images and eosin images on a pixel-by-pixel basis. Next, ‘MeasureGranularity’ module outputed spectra of size measurements of the textures in three types of images. Finally, ‘MeasureTexture’ module measured the degree and nature of textures within three types of images to quantify their roughness and smoothness. (B)Pipeline 2. Haematoxylin-stained images were segmented via ‘IdentifyPrimaryObjects’module and‘IdentifySecondaryObjects’ module ,Quantitative image features of object shape, size, texture, and pixel intensity distribution were further extracted via multiple modules, including measure models of ‘Object Intensity Distribution’, ‘Object Intensity’, ‘Texture’, and ‘Object Size Shape’. Table S1. Baseline characteristics in the training and validation Sets. Table S2. The selected radiomics features. Table S3. The selected pathomics features. Table S4. The selected deep learning pathomics features.

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 The Creative Commons Public Domain Dedication waiver ( 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

Zhang, J., Wu, Q., Yin, W. et al. Development and validation of a radiopathomic model for predicting pathologic complete response to neoadjuvant chemotherapy in breast cancer patients. BMC Cancer 23, 431 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Breast cancer
  • Radiomics
  • Pathomics
  • Radiopathomics
  • CECT