- 1Institute of Biomedical Engineering and Instrumentation, Hangzhou Dianzi University, Hangzhou, China
- 2Department of Mathematics, The Chinese University of Hong Kong, Shatin, Hong Kong SAR, China
- 3Zhejiang Cancer Hospital, Institute of Basic Medicine and Cancer (IBMC), Chinese Academy of Sciences, Hangzhou, Zhejiang, China
Purpose: During neoadjuvant chemotherapy (NACT), breast tumor morphological and vascular characteristics are usually changed. This study aimed to evaluate the tumor shrinkage pattern and response to NACT by preoperative multiparametric magnetic resonance imaging (MRI), including dynamic contrast-enhanced MRI (DCE-MRI), diffuse weighted imaging (DWI) and T2 weighted imaging (T2WI).
Method: In this retrospective analysis, female patients with unilateral unifocal primary breast cancer were included for predicting tumor pathologic/clinical response to NACT (n=216, development set, n=151 and validation set, n=65) and for discriminating the tumor concentric shrinkage (CS) pattern from the others (n=193; development set, n=135 and validation set, n=58). Radiomic features (n=102) of first-order statistical, morphological and textural features were calculated on tumors from the multiparametric MRI. Single- and multiparametric image-based features were assessed separately and were further combined to feed into a random forest-based predictive model. The predictive model was trained in the testing set and assessed on the testing dataset with an area under the curve (AUC). Molecular subtype information and radiomic features were fused to enhance the predictive performance.
Results: The DCE-MRI-based model showed higher performance (AUCs of 0.919, 0.830 and 0.825 for tumor pathologic response, clinical response and tumor shrinkage patterns, respectively) than either the T2WI or the ADC image-based model. An increased prediction performance was achieved by a model with multiparametric MRI radiomic feature fusion.
Conclusions: All these results demonstrated that multiparametric MRI features and their information fusion could be of important clinical value for the preoperative prediction of treatment response and shrinkage pattern.
Introduction
Neoadjuvant chemotherapy (NACT) is routinely used to treat locally advanced tumors to be operable, allowing for breast-preserving surgery (1, 2). Patients who achieve pathologic complete response (pCR) to NACT tend to have improved disease-free and overall survival compared with patients with residual invasive disease (3–5). Nevertheless, only a subset (varying from 10% to 25%) of the patients achieved pCR after NACT, as reported by a larger meta-analysis with 3776 patients (6). Therefore, predicting the response to NACT before treatment is important for the accurate management of breast cancer.
During NACT, breast cancers can also show different regression/shrinkage patterns that are associated with breast cancer treatment outcomes (7). Tumor concentric shrinkage (CS) is associated with good survival in low-grade breast cancer (8) and is regarded as a more suitable marker than pCR to choose candidates for breast-conserving surgery (1). Therefore, the preoperative identification of patients who have a sufficient response to NACT and how the tumor changes during treatment are indispensable in guiding chemotherapy and surgical treatment (9).
Existing studies have attempted to predict the tumor response to NACT before treatment using a magnetic resonance imaging (MRI)-based predictive model (10–12). Radiomics extracted from DCE-MRI was used to predict tumor response in breast cancer (13–15). In addition to DCE-MRI, T2-weighted imaging (T2WI) and apparent diffusion coefficient (ADC) derived from diffuse weighted imaging (DWI) are also used as predictors for evaluating NACT responses (16, 17).
To increase the prediction accuracy, multiparametric MRI (DCE-MRI, T2WI and ADC) features were analyzed and/or combined (18–20). Specifically, DCE-MRI and T2WI radiomics were implemented, which showed that tumor morphologic and quantitative enhancement kinetics can increase the specificity in predicting breast cancer treatment responses (21). Texture analysis of tumors from T2WI and ADC has been conducted for predicting recurrence in breast cancer patients treated with NACT and surgery (22). The addition of histological information of hormone receptor status, Ki-67 index, and MRI variables with radiomics were also used, which showed enhanced discrimination power in predicting tumor response (23).
Despite these advances, how multiparametric images correlate with the tumor response to NACT and the shrinkage pattern still need to be explored. In this study, we sought to enhance the prediction of tumor shrinkage patterns and responses after full cycles of treatment by using radiomics from multiparametric MR image maps and to comprehensively investigate the effectiveness of various parametric images in prediction.
Materials and methods
Patients
This retrospective study, which was approved by the Institutional Review Board, initially included 656 female unilateral breast cancer patients. A total of 216 and 193 unilateral unifocal primary breast cancer patients remained for evaluation of response to NACT and shrinkage pattern (Figure 1) after the selection with an exclusion criterion as follows: 1) surgery or any other treatment before chemotherapy (n=295); 2) at least one of the images among DCE-MRI, T2WI and DWI was missing (n=36); 3) no preoperative image before the initial chemotherapy (n=66); and 4) histopathologic information for residual tumor size was not available (n=43). Additionally, 23 samples were further excluded from the tumor shrinkage pattern analysis dataset because the follow-up data during treatment were not available for determining the type of shrinkage.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1057841/fonc-13-1057841-HTML/image_m/fonc-13-1057841-g001.jpg)
Figure 1 Flow chart of patients included and excluded in tumor pathologic/clinical response and shrinkage pattern datasets.
Treatment response evaluation
All patients underwent standard NACT treatment with six to eight cycles of taxotere–epirubicin–cyclophosphamide (TEC) regimen (24). The duration of the NACT is about 4 or 5 months, and the time interval between the treatments is approximately 20 days. The pCR to NACT was defined as the absence of residual cancer in the surgical specimen. After NACT, 25 patients achieved pCR and 191 patients were non-pCR.
Tumor clinical response was evaluated by monitoring the changes in the longest diameter (LD) of the tumor and kinetic features proposed in the Response Evaluation Criteria in Solid Tumors (RECIST) 1.1 guidelines, which are widely used as recommended tumor measurements (25). According to these guidelines, patients were categorized as having a complete response (CR), partial response (PR), stable disease (SD), or progressive disease (PD). A PR was defined as at least a 30% decrease in the longest tumor diameter (LD). PD indicated an increase of less than 20% or a decrease of less than 30% in the LD. After the completion of the whole treatment procedure, 34 patients achieved CR, 156 achieved PR, one achieved PD and 25 achieved SD. Similar to a previous study, clinical responders (n=190) were defined as PR and CR, and non-clinical responders were defined as PD and SR (n=26) (13).
For tumor shrinkage evaluation, the pattern of shrinkage by NACT can be categorized as CS and non-CS (8). The CS pattern was defined as 1) only reduced tumor size and no small mass surrounding the main tumor or only small foci (less than 5 mm) surrounding the main tumor (26) and 2) CR after treatment. The non-CS pattern was defined as a tumor with 1) diffuse decrease; 2) decrease in intensity only; and 3) no change or enlargement (i.e., PD or SD). Patients in the tumor shrinkage dataset (n=193) were grouped as CS (n=166) and non-CS (n=27).
MRI acquisition
All patients received preoperative MRI examinations in a prone position using a 3.0 T scan system (Siemens Medical Systems, Erlangen, Germany) with a dedicated eight-channel breast array coil. Fat-suppressed bilateral sagittal three-dimensional T2WI, DWI and DCE-MRI were acquired for each patient. For DCE-MRI, one precontract image (S0) followed by five postcontrast images (from S1 to S5) was obtained after injection of contrast agent with a bolus of 0.1 mmol/kg gadobutrol at a rate of 2 mL/second. The time interval between the precontrast and postcontrast image series was 60 seconds. ADC maps were generated from DWIs at b values of 50 and 1000 s/mm2. Detailed MR imaging parameters for DCE-MRI, DWI and T2WI are illustrated in Supplementary Table 1.
Image preprocessing
Image normalization was performed and the pixel value was mapped to a fixed gray area ranging from 0 to 800. A fuzzy C-mean clustering method was used to segment the tumor region of interest (ROI) on DCE-MRI at the intermediate postcontrast series that usually has the highest enhancement to facilitate visual inspection (27). To minimize segmentation errors, the obtained tumor ROI was examined and manually corrected by our investigator. Thereafter, the generated tumor ROI on the DCE-MRI was aligned to the DWI with a b value of 50 and the T2WI map using the Statistical Parametric Modeling toolbox version 12 (SPM12) platform (http://www.fil.ion.ucl.ac.uk/spm/software/spm12). The aligned tumor ROI in DWI was mapped to the ADC image. The precontrast (S0) and first to fifth postcontrast images (termed S1 to S5) were obtained. The subtraction maps, namely, the subtraction between the first, second, third and fifth postcontrast images and S0 (termed S10, S20, S30 and S50, respectively) and that generated by subtraction between the fifth and second postcontrast images (termed S52), were analyzed.
Radiomic feature analysis
Radiomic features (n=102) were extracted from each image series using the publicly available software package pyradiomics (28). These features include morphologic (n=14), first-order statistics (n=18), and texture features (n=70). Specifically, texture features include the gray level cooccurrence matrix (GLCM) (n=24), the gray-level run-length matrix (GLRLM) (n=16), the gray-level size-zone matrix (GLSZM) (n=16) and the gray-level spatial dependence matrix (GLDM) (n=14)-based features. A detailed list of the features is illustrated in Supplementary Table 2. Tumor radiomics was extracted on the DCE-MRI series (from S1 to S5), the subtraction maps (S10, S20, S30, S50 and S52) and the ADC and T2W images.
Pathological examination
The status of the estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor (HER2), and Ki-67 expression level was acquired by using streptavidin-peroxidase (SP) immunohistochemistry (29). Hormone receptor (HR) positivity was determined as ER- or PR-positive, and negativity was defined as both ER- and PR-negative. Ki-67 expression levels greater than 14% were determined to be positive. Breast cancer subtypes were determined as follows: luminal A, HR positive and HER2 negative and low Ki-67 expression level; luminal B, HR positive and either HER2 positive or HER2 negative with high levels of Ki-67; basal like, HR negative and HER2 negative; and HER2-enriched, HR negative and HER2 positive.
Statistical analysis
The differences in the categorized variables, i.e., the histological characteristics between the groups (pCR and non-pCR; CS and non-CS groups; clinical response and non-clinical response), were evaluated using a test. The difference between the continuous values, i.e., age, was evaluated using analysis of variance (ANOVA). Each molecular subtype information was binarized, which was one if it was a particular subtype and zero if it was any other subtype.
A recursive feature elimination (RFE) method was used to rank the features that were most relevant to the target with a random forest model used as a base classifier. The ranked features were then sequentially added into the predictive model, in which the parameters were tuned using a grid search method with a 10-fold cross-validation (CV) framework. The importance of the image features in the predictive model was assessed by using the mean decrease in the Gini score of the random forest over all the CV loops.
Thereafter, the tuned model parameters with the optimal feature subset were used to establish a predictive model using all the samples in the training set. The model performance was assessed on the testing set with an area under the receiver operating characteristic (ROC) curve (AUC).
We have also provided the model performance evaluation by using a nested five-fold CV. Specifically, the inner loop was used to optimize model parameters of the random forest using grid search under ten-fold CV on the training data, while the outer loop produced a prediction score for the testing data. The DeLong test was used to compare the AUCs. The sensitivity and specificity were calculated for the ROC curve by using the Youden index to maximize their summation. P values less than 0.05 were considered significant. Statistical analysis and machine learning methods were performed using R (version 4.0) and MATLAB (MathWorks, Natick, Massachusetts, version 2018 b).
Results
Patient characteristics
Patient characteristics regarding the pathologic response, clinical response and tumor shrinkage pattern are illustrated in Table 1 and Supplementary Tables 3 and 4, respectively. For the pathologic or clinical response dataset (n=216), samples were separated into a development (n=151) and a validation (n=65) set. The tumor shrinkage pattern dataset (n=193) was also separated into development (n=135) and validation (n=58) datasets.
From Table 1 and Supplementary Table 3, there was no significant association between tumor NACT responses (pathologic and clinical responses) and the patient histological characteristics, including age, family history, menopausal status, maximum tumor diameter, PR status, ER status, HER2 status, Ki-67 status, molecular subtypes, and lymph node status. Moreover, no significant differences in these features were observed between the CS and non-CS groups. A significant association between the CS pattern and the tumor response was observed with a p value of 0.02.
Individual parametric MRI radiomic analysis
Individual features associated with the tumor response and shrinkage pattern to NACT were assessed for DCE-MRI series/subtraction images and T2WI and ADC images. The higher voxel volume and the maximum 2D diameter were associated with the non-pCR and non-CS patterns of the tumor during treatment (Figure 2; Supplementary Figure 1, respectively). For DCE-MR image series or maps, individual radiomic features obtained from the subtraction images, particularly the S50 and S30 maps, are better than those obtained using the precontrast (S0) image. For parametric images other than DCE-MRI, radiomics from T2WI is relatively better than that from ADC images.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1057841/fonc-13-1057841-HTML/image_m/fonc-13-1057841-g002.jpg)
Figure 2 Individual feature performance for predicting pathologic response to neoadjuvant chemotherapy. The bar plot is drawn on the right side if a higher value was observed in pCR than in non-pCR, while the bar plot is drawn on the right side if the opposite was true. S10, S20, S30 and S50 represent the subtraction map between the first, second, third and fifth postcontrast images and precontrast image. S52 represents subtraction between the fifth and second postcontrast images.
Figure 3 shows an example of a texture feature (dependent variance) calculated from the subtraction image between the fifth and second postcontrast images (S52). A high level of this feature was significantly (p=0.002) correlated with the CS pattern. Moreover, a similar pattern of higher values in the clinical responders than in the non-clinical responders (p=0.004) was observed for this feature. In other words, a higher dependent variance value in terms of tumor heterogeneity is associated with a good response to NACT or CS. The T2WI-based dependence variance feature was also significantly higher in responders than in non-clinical responders (p=0.008) or higher in the CS pattern and non-CS pattern (p=0.003) (Figure 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1057841/fonc-13-1057841-HTML/image_m/fonc-13-1057841-g003.jpg)
Figure 3 Texture feature of dependence values for discriminating concentric shrinkage (CS) pattern/response and non-CS shrinkage pattern/nonresponse tumors. This texture feature was calculated from the image subtraction between the fifth and the second postcontrast image series (S52). A 64-year-old woman with cancer in the left breast at (A) pretreatment and (B) posttreatment S52 image who responded well to chemotherapy with a CS pattern. A 65-year-old woman with cancer in the right breast at (C) pretreatment and (D) posttreatment S52 who had a poor response to chemotherapy with a non-CS pattern. The boxplot shows a significantly higher dependence variance feature value in the (E) clinical responder group than in the non-clinical responder group, and the same higher value was observed (F) in patients who were CS than in non-CS.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1057841/fonc-13-1057841-HTML/image_m/fonc-13-1057841-g004.jpg)
Figure 4 Texture feature of dependence variance feature at T2 weighting imaging (T2WI) for discriminating concentric shrinkage (CS) pattern/response and non-CS shrinkage pattern/nonresponse tumors. A 40-year-old woman with cancer in the right breast at (A) pretreatment and (B) posttreatment T2WI who has a good response to chemotherapy with a CS pattern. A 48-year-old woman with cancer in the left breast at (C) pretreatment and (D) posttreatment T2WI who had a poor response to chemotherapy with a non-CS pattern. The boxplot shows a significantly higher dependence variance feature value in the (E) clinical responder group than in the non-clinical responder group, and the same higher value was observed (F) in patients who were CS than in non-CS.
Examples of the distributions of multiparametric image features in shrinkage patterns and treatment responses before treatment and after NACT are illustrated (Figure 5). Figure 5A shows an example of a patient with PR to NACT and a CS pattern, while Figure 5B illustrates a patient with nonresponse to treatment and diffuse decrease during NACT. A higher level of the difference variance value was observed in the clinical response/CS patient (37.3 and 41.6 for DCE-MRI and ADC images, respectively) than in the diffuse decrease/non-CS patient (31.2 and 35 for DCE-MRI and ADC images, respectively). Additionally, the high value of the 10th percentile tumor ADC in this patient was associated with a good response/CS pattern to NACT.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1057841/fonc-13-1057841-HTML/image_m/fonc-13-1057841-g005.jpg)
Figure 5 Illustration of the imaging features at pretreatment MRI and post NAC images. From the left to the right columns, the maximum intensity projection (MIP) at the first and third precontrast image series of the DCE-MRI, the third postcontrast image, the T2WI and ADC images are shown. (A) An example of a 47-year-old woman who had a good response to NACT and a concentric shrinkage (CS) pattern at preoperative and postoperative MRI examination. (B) An example of a 42-year-old woman who had a poor response or nonconcentric shrinkage pattern during NACT. The gray level dependence matrix (GLDM)-based dependence variance values from DCE and T2WI were 37.3 and 41.6, respectively, for the response and CS patients and 31.2 and 35.0, respectively, for the nonresponse and non-CS patients. The ADC values of the 10th percentile were 717 and 474 for the response/CS and the nonresponse/non-CS patients, respectively. FV, feature value.
Single parametric image-based predictive model
Radiomics from single parametric images was combined to establish predictive models. DCE-MRI series along with their subtraction maps and the T2WI and ADC images were assessed in the validation set separately (Table 2). Regarding DCE-MRI, the image subtraction-based predictive model had overall better prediction performance than that based on the original image series. Among these, the predictive model using the S50 image had the highest predictive performance for separating patients with pCR from those with non-pCR. The predictive model using the S52 image had the highest prediction performance for tumor clinical response (AUC=0.811, sensitivity of 0.654 at specificity of 1.0), and for tumor shrinkage prediction (AUC = 0.805, sensitivity of 0.820 at specificity of 0.750).
The T2WI-based model achieved an AUC of 0.741 for predicting tumor shrinkage pattern and AUCs of 0.732 and 0.752 for the prediction of pathological and clinical responses to NACT, respectively (Table 3). The predictive model using radiomics from ADC had a relatively lower prediction performance, with AUCs of 0.625, 0.719 and 0.715 for the prediction of tumor pathologic response, clinical response and shrinkage pattern to NACT, respectively (Table 3).
Multiparametric image feature fusion-based prediction
Radiomics from the DCE-MRI series and multiparametric images were used and feature selection was performed. The models were evaluated in the validation set, which generated AUCs of 0.919, 0.830 and 0.825 for tumor pathologic response, clinical response and tumor shrinkage patterns, respectively (Table 3). When the features from the multiparametric images were combined with feature selection, the model generated a better performance in terms of AUCs of 0.943, 0.919 and 0.905 for these three tasks. The multiparametric image-based predictive model was significantly better than that based on T2WI in the prediction of tumor pathological and clinical response to NACT (p=0.009 and p=0.039, respectively). Thereafter, molecular subtype information was included as an additional feature along with features from ADC, T2WI and DCE. For tumor pathologic response prediction, the addition of luminal A, luminal B, HER2 and basal-like tumors generated AUCs of 0.954, 0.934, 0.952 and 0.950, while for tumor clinical response prediction, the AUCs were 0.948, 0.926, 0.913 and 0.914, respectively. For the tumor shrinkage pattern, the performances for the three tasks were improved while separately adding the four molecular subtypes with AUCs of 0.912, 0.905, 0.893 and 0.900. Among the four molecular subtypes, the inclusion of the luminal A subtype showed an optimal increment in the predictive model performance (Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1057841/fonc-13-1057841-HTML/image_m/fonc-13-1057841-g006.jpg)
Figure 6 Receiver operating characteristic (ROC) curves for the models in predicting (A) tumor pathologic (B) clinical response and (C) shrinkage pattern.
Feature importance in the multivariate predictive models for the two tasks was obtained (Table 4). Among them, 10 features were common in the multivariate predictive models, including the informational measure of correlation 1, cluster shade, run length nonuniformity normalized, short-run emphasis from S52, dependence variance and large dependence emphasis from T2WI, the dependence variance from S30, the 90th percentile from S10, and the maximum 2D diameter row and autocorrelation from the ADC. In the multivariate predictive model, the greatest number of features (5 out of 16) was obtained from the subtraction image map (i.e., S52) in DCE-MRI.
The predictive model performance was also evaluated by using nested 5-fold CV (Supplementary Table 5). From this table, similar results to those based on the testing dataset were observed. For example, T2WI based predictive model showed lower prediction performance in pathologic response, clinical response and shrinkage pattern with AUCs of 0.758, 0.708 and 0.731, respectively. The multiparametric image-based predictive model showed improved prediction performance with AUCs of 0.908, 0.868 and 0.843 for these three tasks.
Discussion
A preoperative prediction of the tumor shrinkage pattern and the response to NACT in breast cancer was conducted using multiparametric radiomics from DCE-MRI, T2WI and ADC. Higher tumor heterogeneity evaluated by radiomic features, including the dependent variance value within the tumor, was associated with a good tumor response or CS pattern. The DCE-MRI subtraction image series-based model showed more discriminative power than either T2WI or ADC. Better performance was achieved by combining the multiparametric image features than by using a single series-based model. The addition of luminal A subtype information into the predictive model yields the highest performance.
A previous study used multiparametric images for the prediction of response to NACT and survival outcomes with 38 breast cancer samples (18). Tumor morphology, such as lesion size and volume distribution, is the most discriminative feature. A related multiparametric MRI study identified the tumor sphericity and mean absolute deviation feature value from T2WI and ADC, respectively, for predicting the response to treatment and shrinkage pattern (30). In our study, various types of features, including statistical, morphological and texture features, were identified from the multiparametric images. Among them, the dependence variance value was identified as the most discriminative feature. This feature measures the variance in gray level dependence size, which is defined as the number of connected voxels within a given distance that are dependent on the center voxel in the image. A high value of this feature reflects the high spread of the dependences of the tumor pixel value, which is correlated with a good response to the tumor or CS pattern during NACT. In a previous study, the addition of ER and node status into the radiomic model enhanced the prediction results for pCR or tumor shrinkage size (31). Our results demonstrated that the inclusion of the luminal A subtype data enhanced the predictive power of the model, which indicates a complementary information to imaging features. Women with luminal A breast cancer are less likely to achieve pCR after NACT than those with other molecular subtypes (32).
Comprehensive analysis of multiparametric images in our study demonstrated that the DCE-MRI-based model had higher performance (in terms of AUC) in tumor shrinkage and response prediction than T2WI and ADC images. A related study using multiparametric image radiomics reported a relatively lower prediction performance for a pretreatment ADC-based model (33), which is similar to our study. Compared with a lower sensitivity, the T2WI and ADC images showed a higher specificity in either the shrinkage pattern or response prediction. The combination of ADC, T2WI and DCE-MRI radiomics enhanced the overall model performance.
The predictive effectiveness for DCE-MRI was evaluated in different image series and maps, including the precontrast and subtraction images from different time points. It is interesting to note that the subtraction image had better predictive performance than either the precontrast or the postcontrast images. Among these, the subtraction images (e.g., S50 and S52) have better prediction performance than the other image series. This may be explained by the fact that the subtraction operation enhanced imaging evaluation of tumor angiogenesis and tumor heterogeneity by analyzing dynamic patterns of enhancement (34). The results may suggest that the washout pattern (35) calculated from the subtraction image between the last postcontrast images and the intermediate image (usually showing the highest enhancement) or precontrast image can be better than that using the original images.
Limitations
Our study has several limitations. First, the number of patients in different groups (e.g., CS and non-CS) was imbalanced, which could have induced biased outcomes. Additional studies with larger datasets should be conducted to refine the results. Second, the imaging data were acquired from a single hospital/cohort with unified imaging protocols. The predictive model is limited to these specific imaging protocols.
Conclusions
In summary, radiomics features of multiparametric images were used to predict tumor response and shrinkage patterns during NACT. In addition to molecular subtypes, the high imaging heterogeneity defined by radiomic features is associated with a good response or CS pattern. Further investigation in a larger and external validation dataset with different imaging protocols and longitudinal images (36) is warranted before our method can be used for breast cancer treatment management and 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 authors.
Ethics statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
MF, TZ and LL designed the study. XW, JY, YL, KW, TX performed the image processing, machine learning methods, and statistical analyses. MF and LL wrote the manuscript. TZ and SC analyzed the data. All authors reviewed the manuscript.
Funding
This work was supported by the National Key R&D program of China (2021YFE0203700), the National Natural Science Foundation of China (62271178, U21A20521), the Natural Science Foundation of Zhejiang Province of China (LR23F010002, LGF22H180037) and the Zhejiang Provincial Medicine and Health Science and Technology Plan General Project (2023KY577).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1057841/full#supplementary-material
Abbreviations
CS, Concentric shrinkage; NACT, Neoadjuvant chemotherapy.
References
1. Chen AM, Meric-Bernstam F, Hunt KK, Thames HD, Oswald MJ, Outlaw ED, et al. Breast conservation after neoadjuvant chemotherapy: the MD Anderson cancer center experience. J Clin Oncol (2004) 22(12):2303–12. doi: 10.1200/JCO.2004.09.062
2. Killelea BK, Yang VQ, Mougalian S, Horowitz NR, Pusztai L, Chagpar AB, et al. Neoadjuvant chemotherapy for breast cancer increases the rate of breast conservation: results from the national cancer database. J Am Coll Surg (2015) 220(6):1063–9. doi: 10.1016/j.jamcollsurg.2015.02.011
3. Wang H, Mao X. Evaluation of the efficacy of neoadjuvant chemotherapy for breast cancer. Drug Des Devel Ther (2020) 14:2423–33. doi: 10.2147/DDDT.S253961
4. Cortazar P, Zhang L, Untch M, Mehta K, Costantino JP, Wolmark N, et al. Pathological complete response and long-term clinical benefit in breast cancer: the CTNeoBC pooled analysis. Lancet (2014) 384(9938):164–72. doi: 10.1016/S0140-6736(13)62422-8
5. Guarneri V, Broglio K, Kau SW, Cristofanilli M, Buzdar AU, Valero V, et al. Prognostic value of pathologic complete response after primary chemotherapy in relation to hormone receptor status and other factors. J Clin Oncol (2006) 24(7):1037–44. doi: 10.1200/JCO.2005.02.6914
6. Kong X, Moran MS, Zhang N, Haffty B, Yang Q. Meta-analysis confirms achieving pathological complete response after neoadjuvant chemotherapy predicts favourable prognosis for breast cancer patients. Eur J Cancer (2011) 47(14):2084–90. doi: 10.1016/j.ejca.2011.06.014
7. Goorts B, Dreuning KMA, Houwers JB, Kooreman LFS, Boerma EG, Mann RM, et al. MRI-Based response patterns during neoadjuvant chemotherapy can predict pathological (complete) response in patients with breast cancer. Breast Cancer Res: BCR (2018) 20(1):34. doi: 10.1186/s13058-018-0950-x
8. Fukada I, Araki K, Kobayashi K, Shibayama T, Takahashi S, Gomi N, et al. Pattern of tumor shrinkage during neoadjuvant chemotherapy is associated with prognosis in low-grade luminal early breast cancer. Radiology (2018) 286(1):49–57. doi: 10.1148/radiol.2017161548
9. Sutton EJ, Onishi N, Fehr DA, Dashevsky BZ, Sadinski M, Pinker K, et al. A machine learning model that classifies breast cancer pathologic complete response on MRI post-neoadjuvant chemotherapy. Breast Cancer Res: BCR (2020) 22(1):57. doi: 10.1186/s13058-020-01291-w
10. Reig B, Heacock L, Lewin A, Cho N, Moy L. Role of MRI to assess response to neoadjuvant therapy for breast cancer. J Magnetic Resonance Imag JMRI (2020) 52(6):1587–606. doi: 10.1002/jmri.27145
11. Li Y, Chen Y, Zhao R, Ji Y, Li J, Zhang Y, et al. Development and validation of a nomogram based on pretreatment dynamic contrast-enhanced MRI for the prediction of pathologic response after neoadjuvant chemotherapy for triple-negative breast cancer. Eur Radiology (2021) 32(3):1676–87. doi: 10.1007/s00330-021-08291-0
12. Reig B, Lewin AA, Du L, Heacock L, Toth HK, Heller SL, et al. Breast MRI for evaluation of response to neoadjuvant therapy. Radiographics (2021) 41(3):665–79. doi: 10.1148/rg.2021200134
13. Fan M, Wu G, Cheng H, Zhang J, Shao G, Li L. Radiomic analysis of DCE-MRI for prediction of response to neoadjuvant chemotherapy in breast cancer patients. Eur J Radiology (2017) 94:140–7. doi: 10.1016/j.ejrad.2017.06.019
14. Braman NM, Etesami M, Prasanna P, Dubchuk C, Gilmore H, Tiwari P, et al. Intratumoral and peritumoral radiomics for the pretreatment prediction of pathological complete response to neoadjuvant chemotherapy based on breast DCE-MRI. Breast Cancer Res: BCR (2017) 19(1):57. doi: 10.1186/s13058-017-0846-1
15. Fan M, Cui Y, You C, Liu L, Gu Y, Peng W, et al. Radiogenomic signatures of oncotype DX recurrence score enable prediction of survival in estrogen receptor-positive breast cancer: a multicohort study. Radiology (2022) 302(3):516–24. doi: 10.1148/radiol.2021210738
16. Pereira NP, Curi C, Osorio C, Marques EF, Makdissi FB, Pinker K, et al. Diffusion-weighted magnetic resonance imaging of patients with breast cancer following neoadjuvant chemotherapy provides early prediction of pathological response - a prospective study. Sci Rep (2019) 9(1):16372. doi: 10.1038/s41598-019-52785-3
17. Chamming's F, Ueno Y, Ferre R, Kao E, Jannot AS, Chong J, et al. Features from computerized texture analysis of breast cancers at pretreatment MR imaging are associated with response to neoadjuvant chemotherapy. Radiology (2018) 286(2):412–20. doi: 10.1148/radiol.2017170143
18. Tahmassebi A, Wengert GJ, Helbich TH, Bago-Horvath Z, Alaei S, Bartsch R, et al. Impact of machine learning with multiparametric magnetic resonance imaging of the breast for early prediction of response to neoadjuvant chemotherapy and survival outcomes in breast cancer patients. Invest Radiol (2019) 54(2):110–7. doi: 10.1097/RLI.0000000000000518
19. Liu Z, Li Z, Qu J, Zhang R, Zhou X, Li L, et al. Radiomics of multiparametric MRI for pretreatment prediction of pathologic complete response to neoadjuvant chemotherapy in breast cancer: a multicenter study. Clin Cancer Res an Off J Am Assoc Cancer Res (2019) 25(12):3538–47. doi: 10.1158/1078-0432.CCR-18-3190
20. Huang Y, Chen W, Zhang X, He S, Shao N, Shi H, et al. Prediction of tumor shrinkage pattern to neoadjuvant chemotherapy using a multiparametric MRI-based machine learning model in patients with breast cancer. Front Bioeng Biotechnol (2021) 9:662749. doi: 10.3389/fbioe.2021.662749
21. Taourel P, Pages E, Millet I, Bourgier C, Rouanet P, Jacot W, et al. Magnetic resonance imaging in breast cancer management in the context of neo-adjuvant chemotherapy. Crit Rev Oncol Hematol (2018) 132:51–65. doi: 10.1016/j.critrevonc.2018.09.012
22. Eun NL, Kang D, Son EJ, Youk JH, Kim JA, Gweon HM. Texture analysis using machine learning-based 3-T magnetic resonance imaging for predicting recurrence in breast cancer patients treated with neoadjuvant chemotherapy. Eur Radiology (2021) 31(9):6916–28. doi: 10.1007/s00330-021-07816-x
23. Kim SY, Cho N, Choi Y, Lee SH, Ha SM, Kim ES, et al. Factors affecting pathologic complete response following neoadjuvant chemotherapy in breast cancer: development and validation of a predictive nomogram. Radiology (2021) 299(2):290–300. doi: 10.1148/radiol.2021203871
24. Gu X, Zhang Y, Chen L, Guo J, Zhang WH. Efficacy of neo-adjuvant chemotherapy with TEC regimen on breast cancer. Cancer Chemother Pharmacol (2015) 75(2):301–8. doi: 10.1007/s00280-014-2646-y
25. Eisenhauer EA, Therasse P, Bogaerts J, Schwartz LH, Sargent D, Ford R, et al. New response evaluation criteria in solid tumours: revised RECIST guideline (version 1.1). Eur J Canc (2009) 45(2):228–47. doi: 10.1016/j.ejca.2008.10.026
26. Green MC, Buzdar AU, Smith T, Ibrahim NK, Valero V, Rosales MF, et al. Weekly paclitaxel improves pathologic complete remission in operable breast cancer when compared with paclitaxel once every 3 weeks. J Clin Oncol (2005) 23(25):5983–92. doi: 10.1200/JCO.2005.06.232
27. Fan M, Zhang P, Wang Y, Peng W, Wang S, Gao X, et al. Radiomic analysis of imaging heterogeneity in tumours and the surrounding parenchyma based on unsupervised decomposition of DCE-MRI for predicting molecular subtypes of breast cancer. Eur Radiology (2019) 29(8):4456–67. doi: 10.1007/s00330-018-5891-3
28. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational radiomics system to decode the radiographic phenotype. Cancer Res (2017) 77(21):e104–7. doi: 10.1158/0008-5472.CAN-17-0339
29. Hammond ME, Hayes DF, Wolff AC, Mangu PB, Temin S. American Society of clinical oncology/college of american pathologists guideline recommendations for immunohistochemical testing of estrogen and progesterone receptors in breast cancer. J Oncol Pract (2010) 6(4):195–7. doi: 10.1200/JOP.777003
30. Zhuang X, Chen C, Liu Z, Zhang L, Zhou X, Cheng M, et al. Multiparametric MRI-based radiomics analysis for the prediction of breast tumor regression patterns after neoadjuvant chemotherapy. Transl Oncol (2020) 13(11):100831. doi: 10.1016/j.tranon.2020.100831
31. Yan S, Wang W, Zhu B, Pan X, Wu X, Tao W. Construction of nomograms for predicting pathological complete response and tumor shrinkage size in breast cancer. Cancer Manag Res (2020) 12:8313–23. doi: 10.2147/CMAR.S270687
32. 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. doi: 10.1007/s10549-018-4801-3
33. Eun NL, Kang D, Son EJ, Park JS, Youk JH, Kim JA, et al. Texture analysis with 3.0-T MRI for association of response to neoadjuvant chemotherapy in breast cancer. Radiology (2020) 294(1):31–41. doi: 10.1148/radiol.2019182718
34. Karahaliou A, Vassiou K, Arikidis NS, Skiadopoulos S, Kanavou T, Costaridou L. Assessing heterogeneity of lesion enhancement kinetics in dynamic contrast-enhanced MRI for breast cancer diagnosis. Br J Radiology (2010) 83(988):296–309. doi: 10.1259/bjr/50743919
35. Wu J, Gong G, Cui Y, Li R. Intratumor partitioning and texture analysis of dynamic contrast-enhanced (DCE)-MRI identifies relevant tumor subregions to predict pathological response of breast cancer to neoadjuvant chemotherapy. J Magnetic Resonance Imag JMRI (2016) 44(5):1107–15. doi: 10.1002/jmri.25279
Keywords: breast cancer, radiomics, multiparametric MRI, neoadjuvant chemotherapy, concentric shrinkage
Citation: Fan M, Wu X, Yu J, Liu Y, Wang K, Xue T, Zeng T, Chen S and Li L (2023) Multiparametric MRI radiomics fusion for predicting the response and shrinkage pattern to neoadjuvant chemotherapy in breast cancer. Front. Oncol. 13:1057841. doi: 10.3389/fonc.2023.1057841
Received: 30 September 2022; Accepted: 19 April 2023;
Published: 03 May 2023.
Edited by:
Valeria Romeo, University of Naples Federico II, ItalyReviewed by:
Bingsheng Huang, Shenzhen University, ChinaFrédérique Frouin, Centre de Recherche de l’Institut Curie, France
Copyright © 2023 Fan, Wu, Yu, Liu, Wang, Xue, Zeng, Chen and Li. 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: Shujun Chen, chensj@zjcc.org.cn; Lihua Li, lilh@hdu.edu.cn