Interim heterogeneity changes measured using entropy texture features on T2-weighted MRI at 3.0 T are associated with pathological response to neoadjuvant chemotherapy in primary breast cancer

European Radiology, May 2017

Objectives To investigate whether interim changes in hetereogeneity (measured using entropy features) on MRI were associated with pathological residual cancer burden (RCB) at final surgery in patients receiving neoadjuvant chemotherapy (NAC) for primary breast cancer. Methods This was a retrospective study of 88 consenting women (age: 30–79 years). Scanning was performed on a 3.0 T MRI scanner prior to NAC (baseline) and after 2–3 cycles of treatment (interim). Entropy was derived from the grey-level co-occurrence matrix, on slice-matched baseline/interim T2-weighted images. Response, assessed using RCB score on surgically resected specimens, was compared statistically with entropy/heterogeneity changes and ROC analysis performed. Association of pCR within each tumour immunophenotype was evaluated. Results Mean entropy percent differences between examinations, by response category, were: pCR: 32.8%, RCB-I: 10.5%, RCB-II: 9.7% and RCB-III: 3.0%. Association of ultimate pCR with coarse entropy changes between baseline/interim MRI across all lesions yielded 85.2% accuracy (area under ROC curve: 0.845). Excellent sensitivity/specificity was obtained for pCR prediction within each immunophenotype: ER+: 100%/100%; HER2+: 83.3%/95.7%, TNBC: 87.5%/80.0%. Conclusions Lesion T2 heterogeneity changes are associated with response to NAC using RCB scores, particularly for pCR, and can be useful across all immunophenotypes with good diagnostic accuracy. Key Points • Texture analysis provides a means of measuring lesion heterogeneity on MRI images. • Heterogeneity changes between baseline/interim MRI can be linked with ultimate pathological response. • Heterogeneity changes give good diagnostic accuracy of pCR response across all immunophenotypes. • Percentage reduction in heterogeneity is associated with pCR with good accuracy and NPV.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1007%2Fs00330-017-4850-8.pdf

Interim heterogeneity changes measured using entropy texture features on T2-weighted MRI at 3.0 T are associated with pathological response to neoadjuvant chemotherapy in primary breast cancer

Interim heterogeneity changes measured using entropy texture features on T2-weighted MRI at 3.0 T are associated with pathological response to neoadjuvant chemotherapy in primary breast cancer Shelley Henderson 0 1 2 3 4 5 Colin Purdie 0 1 2 3 4 5 Caroline Michie 0 1 2 3 4 5 Andrew Evans 0 1 2 3 4 5 Richard Lerski 0 1 2 3 4 5 Marilyn Johnston 0 1 2 3 4 5 Sarah Vinnicombe 0 1 2 3 4 5 Alastair M. Thompson 0 1 2 3 4 5 0 Department of Oncology, Ninewells Hospital and Medical School , Dundee, UK DD1 9SY 1 Department of Pathology, Ninewells Hospital and Medical School , Dundee, UK DD1 9SY 2 Department of Medical Physics, Ninewells Hospital and Medical School , Dundee, UK DD1 9SY 3 Department of Breast Surgical Oncology, University of Texas MD Anderson Cancer Centre , Houston, TX 77030 , USA 4 Department of Clinical Radiology, Ninewells Hospital and Medical School , Dundee, UK DD1 9SY 5 Division of Imaging and Technology, Ninewells Hospital and Medical School, University of Dundee , Dundee, UK DD1 9SY Objectives To investigate whether interim changes in hetereogeneity (measured using entropy features) on MRI were associated with pathological residual cancer burden (RCB) at final surgery in patients receiving neoadjuvant chemotherapy (NAC) for primary breast cancer. Methods This was a retrospective study of 88 consenting women (age: 30-79 years). Scanning was performed on a 3.0 T MRI scanner prior to NAC (baseline) and after 2-3 cycles of treatment (interim). Entropy was derived from the grey-level co-occurrence matrix, on slice-matched baseline/ interim T2-weighted images. Response, assessed using RCB score on surgically resected specimens, was compared statistically with entropy/heterogeneity changes and ROC analysis performed. Association of pCR within each tumour immunophenotype was evaluated. Magnetic resonance imaging; Breast cancer; Neoadjuvant chemotherapy; Lesion hetereogeneity; Pathogical response - * Shelley Henderson Neoadjuvant chemotherapy (NAC) for primary breast cancer is clinically useful to downstage locally advanced breast cancer and/or reduce the extent of surgery from mastectomy to breast conservation. Response to NAC varies by subtype of breast cancer and chemotherapy regimen [1, 2] but prediction of response remains clinically challenging [3, 4]. A pathological complete response (pCR) to NAC is intrinsically linked with survival outcomes, particularly in triple negative and hormone receptor positive disease [5–7] and, in the future, where confidence is high that a pCR has been achieved, post-treatment surgery might not be required [8, 9]. Therefore, early identification of such patients becomes crucial for improving and personalising patient management. Conversely, identification of patients unlikely to respond well early in treatment may allow modification to the cytotoxic regimen or dose, although the benefits of such a strategy would need to be proven in prospective trials [10]. Breast magnetic resonance imaging (MRI) is increasingly used to monitor patients during NAC therapy, allowing repeated examinations, with dynamic contrast-enhanced imaging (DCE) providing a measure of tumour vascularity. However, breast MRI can under- or overestimate the amount of residual tumour in around 30% of women [11], and the accuracy of response prediction is dependent on immunohistochemical tumour subtype [12]. Quantitative imaging biomarkers play an important role in repeat imaging studies [13], particularly when comparative changes in response to a therapy require quantification. To this end, a number of MRI metrics have been explored, including MRI DCE parameters [14], volume changes [3, 15] and diffusion changes [16]. The application of textural analysis (TA) to medical imaging studies shows promise. TA statistically models spatial distributions of pixel grey-levels in order to classify or segment the data [17] and has been used successfully in neuroimaging for tissue differentiation and classification [18, 19]. For breast MRI, TA has been used to differentiate between malignant and benign lesions [20, 21] and in discrimination between breast cancer subtypes [22–24]. Recently, heterogeneity, as assessed using TA, has found a role in the diagnosis, characterization [25] and treatment response assessment of various cancer types [26, 27], with initial reports suggesting an association with patient survival [28]. Preliminary work utilising TA [29] indicates that heterogeneity changes in response to NAC may be seen using T2weighted imaging, as well as on DCE-MRI. T2-weighted imaging has the advantage that signal intensity is directly related to underlying morphology (via intracellular and extracellular/ extravascular space) and can provide useful information in cases where contrast material cannot be administered. T2weighted images are also less susceptible to degradation (e.g. movement) and are likely to be easier to standardise across multiple sites. As previous studies using TA have been relatively small, we explored application of these principles in a larger patient cohort, correlating with a recognised, validated measure of pathological response to NAC for breast cancerthe residual cancer burden score (RCB) [30]. Patients with a pathologic complete response (pCR) have favourable survival rates compared with those having minimal, moderate or extensive residual disease (RCB-I/RCB-II/RCB-III) [30]. This study sought to determine if changes in lesion heterogeneity could be linked with response to NAC, as measured using ultimate pathological outcome and assess whether the r e s p o n s e a s s o c i a t i o n w a s d e p e n d e n t o n c a n c e r immunophenotype. Measurement repeatability was also considered, as is essential when evaluating a potential quantitative tool. Patients were scanned prior to NAC, after NAC cycle 2 (n = 16) or 3 (n = 72) and prior to surgery to assess final imaging response. At interim MRI examination, all patients had received only FEC treatment. The median time between final Assessed for eligibility (n =117 patients) NAC Study Cohort (n = 88 patients) Excluded (n = 29 patients) Received first cycle of NAC prior to MRI examination (n = 14) Interim MRI examination performed after 4 cycles of NAC (n = 6) Excessive movement during MRI examination (n=4) Proceeded to surgery after interim MRI examination (n=2) No final pathologic information available (n=2) Non-standard treatment (n=1) 6xFEC (n=18) 3xFEC; 3 DOC (n=41) 3xFEC; 3xDOC+TRA (n=26) Fig. 1 Diagram of recruitment and study population including neoadjuvant chemotherapy (NAC) regime and response as measured using residual cancer burden (RCB) score. FEC fluorouracil, epirubicin (75) and cyclophosphamide, DOC docetaxel, TRA trastuzumab, TDM-1 trastuzumab emtansine MRI examination and surgical resection was 28 days (range: 6–57 days). All MR examinations were carried out using a 32channel 3.0 Tesla (T) Siemens Magnetom Trio scanner (Erlangen, Germany) and a 7-channel breast coil. A T2 Patient and tumour characteristics Patients (total n = 88) weighted sequence was acquired using a turbo-spin echo sequence (TR/TE/α = 8,000 ms/86 ms/150°, field of view = 340 × 340 mm, matrix = 320 × 320, voxel size 1.1 mm3), with parallel imaging factor × 2 and turbo factor × 21. The rest of the examination consisted of a T1 non-fat-saturated acquisition followed by the DCE acquisition. DICOM images were anonymised and exported off-line to perform TA. Texture analysis was performed using MaZda version 4.7 [32] on baseline and interim T2-weighted images. Postcontrast subtracted images from both examinations of each patient were considered side-by-side to identify maximum lesion diameter on each; subsequently these slices were matched to the appropriate slice on the T2 acquisition for region of interest (ROI) positioning on comparable slices (see Fig. 2). All image analysis was performed blinded to patient outcome, clinical and pathology information. T2-weighted images were magnified to visualise the tumour border and an ROI drawn around the entire tumour on the imaging slice demonstrating maximum lesion diameter, with no marker clip present. Histograms were visualised to ensure no fat was erroneously present within Fig. 2 Region of interest (ROIs) (red) drawn for texture analysis on two patients, with slices matched between baseline (left column) and interim (right column) examinations. Insert images highlight lesion locations without ROI overlay. Image (a) and (b) are from slice-matched T2-weighted images from a woman who had a RCB-III at final pathology, while (c) and (d) are from slice-matched images in a patient that ultimately achieved a pCR at end of treatment. Images (e) and (f) highlight the importance of correlation with the subtracted dynamic contrastenhanced imaging (DCE) images, particularly in the case of non-mass enhancement, as shown the ROI (see Fig. 2). A 10 × 10 pixel ROI was placed in healthy, normal contralateral breast parenchyma to exclude systematic inter-scan differences (as shown in Fig. 2(c,d)). These were placed in comparable regions of breast parenchyma at baseline and interim to ensure consistency and care was taken to exclude presence of vessels within the ROI. Data was rescaled to 6 bits/pixel and the image histogram normalized to within 3 standard deviations of the mean to minimise brightness and contrast variations. Entropy features, representing heterogeneity, were derived from the grey-level co-occurrence matrix [33], as these are reported in the literature as most appropriate for cancer imaging [29]. Inter-pixel distances of n = 2 and n = 5 were calculated to depict fine and coarse texture. Raw feature values were exported for statistical analysis. Intra- and interobserver error was calculated based on a subset of 20 baseline images that were analysed twice by SAH (8 years TA experience), with a 1-month interval, and once by RAL (30+ years TA experience). Assessment of response RCB score [30], was calculated by a specialist breast pathologist (CP) using final pathology of resected specimen, based on tumour bed dimensions, cellularity and axillary node burden and reported in terms of the RCB index [30]. Patients were categorised in terms of whether they achieved a pCR post treatment, or had residual disease (RCB-I, RCB-II, RCB-III). All statistics were calculated using SPSS (v21; IBM Corporation, New York, NY, USA) with p < 0.05 considered significant, and p < 0.001 highly significant. Intra- and interobserver measurement error of entropy features were calculated using coefficients of variation (CoVs) and Pearson’s intraclass correclation coefficient (ICC). Data was not normally distributed therefore non-parametric tests were used throughout. Absolute baseline and interim entropy values were compared using Wilcoxon tests in both normal and lesion tissue, and percentage changes between examinations assessed using Kruskal-Wallis tests (across all response categories) and Mann-Whitney U-tests (between response categories). Receiver operating characteristic (ROC) curves were generated and areas under the curve (AUCs) calculated. Optimal thresholds were derived using Youden’s index, and resultant sensitivity, specificity, accuracy, and positive and negative predictive values (PPVs and NPVs) calculated. Of 88 patients analysed, 107 lesions were identified; however, TA was restricted to the largest, index lesion in each case. At final pathological examination, 16 patients were categorized as pCR, 12 as RCB-I, 43 as RCB-II and 17 as RCB-III (Fig. 1). Intra-observer ICC values demonstrated a high correlation for entropy features (ICC = 0.813, p < 0.001) with corresponding CoV value of 5.4%. Inter-observer ICC values demonstrated similar correlations, with entropy ICC = 0.790 (p < 0.001) and a slightly poorer CoV value of 8.0%. Mean baseline ROI size was 795 pixels (range: 202–1,995), corresponding to a physical size of 962 mm2 (244– 2,414 mm2); interim ROI size was 454 pixels (77–1,561) with physical size 549 mm2 (93–1,889 mm2). There were no significant differences in baseline TA features between response categories (p > 0.866; Kruskal-Wallis). No significant differences were found within normal parenchyma between baseline and interim examinations (p > 0.419; Wilcoxon test) or in pair-wise comparisons of difference in entropy features between response groups (p > 0.166; MannWhitney U-test) (Fig. 3), suggesting no random or systematic inter-scan temporal differences. Patients with pCR, RCB-I or RCB-II scores had decreases in fine and coarse entropy features between baseline and interim examinations (p < 0.006; Wilcoxon test), indicating lesion image heterogeneity reduction. No significant differences were found for RCB-III patients (p > 0.397; Wilcoxon test) (Table 2, Fig. 3). Median decreases in coarse entropy between baseline and interim MRI were pCR: 32.8%, RCB-I: 10.5%, RCB-II 9.7% and RCB-III: 3.0% (Fig. 4). For both entropy features, patients with pCR had a significantly greater reduction than those patients with residual disease (p < 0.014; Mann-Whitney U-test). The ROC curve for fine and coarse entropy is shown in Fig. 5, with calculated AUC for pCR classification of 0.834 using fine features and 0.845 using coarse features. Optimal thresholds for pCR identification were derived as 11% (fine entropy) and 20% (coarse entropy) using Youden’s index. Resulting sensitivities, specificities, accuracies, NPVs and PPVs are shown in Table 3. T u m o u r s w e r e t h e n c o n s i d e r e d i n t e r m s o f immunophenotype. There were significant differences between pCR and residual disease groups when considered across all immunophenotypes (ER+: p < 0.001; HER2+: p < 0.007; TNBC: p < 0.001; Mann-Whitney U-test). Median percentage changes in coarse entropy categories for ER+, HER2+ and TNBC (Table 4) confirmed reduction in lesion image heterogeneity was associated with pathological response as assessed by RCB with changes in heterogeneity being greatest in the pCR group with ER+ tumours. Patients with TNBC who ultimately achieved a pCR had a lower reduction in heterogeneity with a median 23% reduction indicating a likely pCR post-NAC. Optimal thresholds were derived using Youden’s index as ER: 28%, HER2+: 30% and TNBC: 18%. These resulted in good diagnostic characteristics, as presented in Table 5. Neoadjuvant chemotherapy is increasingly used in breast cancer management, particularly for TNBC and HER2+ disease, and a pCR is directly related to overall survival [5–7, 34]. Image analysis carries potential to assess response during therapy, prior to surgical resection, but most studies using breast MRI have focused on changes in DCE or apparent diffusion coefficient [16, 35]. This study considered changes in T2Fig. 3 Coarse entropy features at baseline and interim examinations in both normal tissue and lesion in each response category. Within normal tissue there were no significant differences between visits for any of the residual cancer burden (RCB) response categories (p > 0.419; Wilcoxon test), while within the lesion entropy features demonstrate a reduction in lesion heterogeneity at the interim examination weighted MRI entropy features (representing heterogeneity) during NAC and correlated with pathological RCB score posttreatment. We demonstrate a significant association between Table 2 Median differences between baseline and interim examinations for coarse and fine entropy in each residual cancer burden (RCB) response category, and Wilcoxon statistical tests (*indicates significance at p < 0.05 level) pCR pathological complete response RCB-I RCB-II RCB-III Response Baseline Interim Lesion Normal heterogeneity changes and RCB, and thus, indirectly, with clinical outcome [30, 34]. Importantly, the technique has performed well for each immunophenotype and thus further investigation is of real clinical importance. Absolute differences in fine entropy features between baseline and interim examinations dichotomised data into two statistically distinct and clinically important groups; pCR and RCB-I, II and III (residual disease). ROC analysis demonstrates that a threshold reduction of 20% in coarse entropy resulted in an association of patients who would ultimately have a pCR, with an accuracy of 85.2% (AUC 0.845). Importantly, these thresholds result in a very good NPV for pCR (Tables 3 and 5), therefore clinicians could potentially identify at interim MRI patients unlikely to achieve a pCR on therapy completion, with over 95% accuracy. Based on our results, coarse entropy features perform better than fine features, with greater differences between baseline and interim MRI. This larger threshold in heterogeneity reduction for pCR leads to improved diagnostic accuracy (Table 3, Fig. 5). Fine and coarse entropy use has been reported elsewhere [29] and relates to the scale at which structural features are enhanced within textural findings. Previous hypotheses Fig. 4 Percentage reduction in coarse entropy features between baseline and interim examinations for each response category as measured using the residual cancer burden (RCB) score. Statistical comparisons are indicated (#p < 0.05, ##p < 0.001; Mann-Whitney U-test) suggest fine features relate to parenchyma while coarser features represent underlying vasculature; however these are based on CT images and have not been validated in MR imaging. No significant differences were found in heterogeneity of the contralateral breast, lending weight to findings being attributable to therapeutic effects. While chemotherapy-induced changes have been reported within normal breast tissue [36], most likely relating to vascular damage [37], these are less likely to affect T2-weighted images. There is scant literature concerning repeatability of TA in clinical cohorts; however, our intra-observer measurement error of CoV = 5.4% is in line with Carballido-Gamio et al. who report 2.8–6.6% intra-observer repeatability in entropy Fig. 5 Receiver operating characteristic (ROC) curve for fine and coarse entropy features in the classification of pathological complete response (pCR). Areas under the curve are 0.834 and 0.845, respectively Table 3 Test performance results in the identification of pathological complete responders (pCRs) when using reduction in entropy features – both coarse and fine Table 5 Test performance results in the identification of pathological complete responders (pCRs) when using reduction in coarse entropy features for ER+, HER2+ and triple negative breast cancers (TNBCs) Entropy coarse (20% reduction) Entropy fine (11% reduction) AUROC area under the receiver operating characteristic curve, PPV positive-predictive value, NPV negative-predictive value measures of knee cartilage [38] and Mathias et al. who quote 1.1% within the spinal cord [39]. We can find no studies that consider inter-observer repeatability measures, however, estimate our CoV at 8%. With estimated error below 10%, this could potentially mean the technique could be comparable across multicentre studies. The reduction in entropy has previously been reported [29]; however, we have extended this work to consider the association of pCR versus residual disease at the end of treatment. While this was a pilot study, our patient numbers were larger than comparable studies [29, 40, 41], an accepted measurement of pathological response was used [30] and the important consideration of immunophenotype was included. Although baseline prediction of ultimate response to NAC would have real clinical importance, we found no link in our cohort of patients. Reports in the literature on baseline heterogeneity response prediction are contradictory; with Ahmed et al. reporting increased pre-treatment entropy in association with poor NAC response [42], whereas studies by Teruel [40] and Michoux [41] suggest the opposite. It should be noted, h o w e v e r , t h a t n o n e o f t h e s e s t u d i e s c o n s i d e r immunophenotype, which is particularly relevant in the NAC response setting. Prediction of pCR at interim MRI has been widely investigated, and two recent papers have considered various metrics for this purpose. O’Flynn et al. evaluated a range of Table 4 Median percentage changes in lesion image heterogeneity as assessed using coarse entropy features for ER+, HER2+ and TNBC in each response category (n = number of cancers per group) ER+ (n = 26) HER2+ (n = 29) TNBC (n = 33) TNBC triple negative breast cancer, ER oestrogen receptor, RCB residual cancer burden ER+ (n = 26) HER2+ (n = 29) TNBC (n = 33) a threshold derived using Youdens index PPV positive predictive value, NPV negative predictive value parameters including volumes, pharmacokinetics and enhancement ratios. They found the best pCR predictor was volumetric changes, which resulted in an area under the ROC curve of 0.773 and sensitivity/specificity of 71.4%/ 76.9% [43]. Li et al. focused on optimising enhancement thresholds for calculating tumour volumes and obtained a value of AUC = 0.73 [44]. This group also considered immunophenotype with a good AUC of 0.77 (ER+/HER2-), 0.75 (HER2+) and 0.85 (TNBC). Our study has shown higher diagnostic accuracies, and it would therefore be of great interest to evaluate the two combined techniques in future studies. Given the association with RCB score, TA has the potential to support trial development where, after an imaging assessment of complete response to NAC, no resection would be performed, only percutaneous sampling of the tumour bed and radiotherapy treatment to the breast. Such a concept, while not new [9, 45], may have a resurgence of interest given the increasing efficacy of NAC, if a sufficiently accurate noninvasive modality for monitoring, such as MRI TA, were available. Texture analysis has wide application in cancer imaging for differentiation, classification and treatment monitoring [22–29, 40, 42]. Spatial resolution is key in generalisability across scanners [46], and while the technique is computationally demanding, both in analysis and interpretation, incorporation into commercial, validated post-processing software should be possible. Presently, different software applications result in entropy values with different orders of magnitude and therefore utility of absolute entropy features in diagnosis is unlikely to come to fruition until standardisation occurs. However, the most useful, and perhaps clinically relevant, utility is in comparative studies – be that comparison of tissue types or in follow-up studies, such as treatment monitoring. Currently, the precise correlate of TA with underlying anatomical or pathological structures is unknown and other groups are working on linking imaging heterogeneity with underlying histopathological and genetic heterogeneity [24]. While most TA of breast images have been performed on post-contrast images, these are subject to movement, and registration algorithms can introduce artificial texture into the images. In a recent comparative study, TA of T2 images was found to be more discriminatory between clinical response groups in patients undergoing NAC than DCE-MRI [29]. The underlying mechanism for signal intensity distributions in DCE-MRI and T2-weighted imaging is entirely different. Signal intensity in DCE-MRI is derived from contrast agent distribution and vascular leakage, while T2-weighted imaging is predominantly influenced by intracellular and extracellular/ extravascular space and is therefore more likely to directly reflect tumour morphology. The heterogeneous nature of cancer subtypes included in this study led to the use different chemotherapy regimes, which reflects clinical reality. However, interim analysis was performed on patients who had only undergone FEC treatment. Another potential weakness of the study is that interim scans were performed after two or three cycles of NAC, yet despite this variability we were still able to show statistically significant differences between response categories. Finally, there is a potential for a slight mismatch in analysed imaging slice between baseline and interim MRI examination. This mismatch is likely to be small, especially in relation to tumour size changes, imaging slice thickness and the use of anatomical landmarks and marker clips for localisation of ROI positioning. In conclusion, lesion heterogeneity changes are associated with response to NAC using RCB scores, and could be useful in identification of patients likely to achieve a pCR, across all immunophenotypes. These data suggest changes in MRI entropy features may provide a clinically useful early indication of response to NAC and warrants further investigation and consideration in multivariate analyses for prediction of pathological results. Acknowledgements The authors would like to thank Breast Cancer Now for partial funding of this research, as well as the patients who took part in this study. Compliance with ethical standards Conflict of interest The authors of this manuscript declare no relationships with any companies whose products or services may be related to the subject matter of the article. This study has received partial funding by Breast Cancer Statistics and biometry sary for this paper. No complex statistical methods were necesEthical approval Institutional Review Board approval was not required because of the retrospective nature of the study of anonymised images. Informed consent Written informed consent was obtained from all subjects (patients) in this study. performed at one institution Open Access This article is distributed under the terms of the Creative C o m m o n s A t t r i b u t i o n 4 . 0 I n t e r n a t i o n a l L i c e n s e ( h t t p : / / creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 1. Esserman LJ , Berry DA , Cheang MCU et al ( 2012 ) Chemotherapy response and recurrence-free survival in neoadjuvant breast cancer depends on biomarker profiles: results from the I-SPY 1 TRIAL (CALGB 150007/150012;ACRIN 6657) . Breast Cancer Res Treat 132 : 1049 - 1062 2. Houssami N , Macaskill P , von Minckwitz G , Marinovich ML , Mamounas E ( 2012 ) Meta-analysis of the association of breast cancer subtype and pathologic complete response to neoadjuvant chemotherapy . Eur J Cancer 48 : 3342 - 3354 3. Hylton NM , Blume JD , Bernreuter WK et al ( 2012 ) Locally advanced breast cancer: MR imaging for prediction of response to neoadjuvant chemotherapy- results from the ACRIN 6657/I- SPY TRIAL. Radiology 263 : 663 - 672 4. Loo CE , Straver ME , Rodenhuis S et al ( 2011 ) Magnetic resonance imaging response monitoring of breast cancer during neoadjuvant chemotherapy: relevance of breast cancer subtype . J Clin Oncol 29 : 660 - 666 5. Spring L , Greenup R , Reynolds K , Smith BL , Moy B , Bardia A ( 2016 ) Abstract 1439: Pathological complete response after neoadjuvant chemotherapy predicts improved survival in all major subtypes of breast cancer: systematic review and meta-analyses of over 18,000 patients . Cancer Res 76 ( 14 Suppl): 1439 6. Cortazar P , Zhang L , Untch M et al ( 2014 ) Pathological complete response and long term clinical benefit analysis in breast cancer: the CTNeoBC pooled analysis . Lancet 384 : 164 - 172 7. Minckwitz G , Untch M , Blohmer JU et al ( 2012 ) Definition and impact of pathologic complete response on prognosis after neoadjuvant chemotherapy in various intrinsic breast cancer subtypes . J Clin Oncol 30 : 1796 - 1804 8. Schaefgen B , Mati M , Sinn HP et al ( 2016 ) Can routine imaging after neoadjuvant chemotherapy in breast cancer predict pathologic complete response ? Ann Surg Oncol 26 : 789 - 795 9. la Parra V , Kuerer HM ( 2016 ) Selective elimination of breast cancer surgery in exceptional responders: historical perspective and current trials . Breast Cancer Res 18 : 28 - 35 10. Barker AD , Sigman CC , Kelloff GJ , Hylton NM , Berry DA , Esserman LJ ( 2009 ) I-SPY 2: An adaptive breast cancer trial design in the setting of neoadjuvant chemotherapy . Clin Pharmacol Ther 86 : 97 - 100 11. Yeh R , Slanetz P , Kopans DB et al ( 2005 ) Prospective comparison of mammography, sonography, and MRI in patients undergoing neoadjuvant chemotherapy for palpable breast cancer . AJR Am J Roentgenol 184 ( 3 ): 868 - 877 12. De Los SJ , Cantor A , Amos KD et al ( 2013 ) Magnetic resonance imaging as a predictor of pathologic response in patients treated with neoadjuvant systemic treatment for operable breast cancer . Cancer 119 ( 10 ): 1776 - 1783 13. Quantitative Imaging Biomarkers Alliance. Radiological Society of North America Website. https://www.rsna.org/QIBA.aspx. Accessed October 21st , 2015 . 14. Pickles MD , Lowry M , Manton DJ , Gibbs P , Turnbull LW ( 2005 ) Role of dynamic contrast enhanced MRI in monitoring early response of locally advanced breast cancer to neoadjuvant chemotherapy . Breast Cancer Res Treat 91 : 1 - 10 15. Partridge SC , Gibbs JE , Lu Y et al ( 2005 ) MRI measurement of breast tumour volume predict response to neoadjuvant chemotherapy and recurrence-free survival . AJR Am J Roentgenol 184 ( 6 ): 1774 - 1781 16. Pickles MD , Gibbs P , Lowry M , Turnbull LW ( 2006 ) Diffusion changes precede size reduction in neoadjuvant treatment of breast cancer . Magn Reson Imaging 24 ( 7 ): 843 - 847 17. Castellano G , Bonilha L , Li LM , Cendes F ( 2004 ) Texture analysis of medical images . Clin Radiol 59 : 1061 - 1069 18. Nedelec J-F , Yu O , Chambron J , Macher J-P ( 2004 ) Texture analysis of the brain: from animal models to human applications . Dialogues Clin Neurosci 6 ( 2 ): 227 - 233 19. Yu O , Mauss Y , Namer IJ , Chambron J ( 2001 ) Existence of contralateral abnormalities revealed by texture analysis in unilateral intractable hippocampal epilepsy . Magn Reson Imaging 19 : 1305 - 1310 20. Gibbs P , Turnbull LW ( 2003 ) Textural analysis of contrastenhanced MR images of the breast . Magn Reson Med 50 : 92 - 98 21. Sinha S , Lucas-Quesada FA , DeBruhl ND et al ( 1997 ) Multifeature analysis of Gd-enhanced MR images of breast lesions . J Magn Reson Imaging 7 ( 6 ): 1016 - 1026 22. Holli K , Laaperi AL , Harrison L et al ( 2010 ) Characterisation of breast cancer types by texture analysis of magnetic resonance images . Acad Radiol 17 ( 2 ): 135 - 141 23. Waugh SA , Purdie CA , Jordan LB et al ( 2016 ) Magnetic resonance imaging texture analysis classification of primary breast cancer . Eur Radiol 26 ( 2 ): 322 - 330 24. Grimm LJ ( 2016 ) Breast MRI , radiogenomics: current status and research implications . JMRI 43 : 1269 - 1278 25. Aerts HJWL , Velazquez RE , Leijenaar RTH et al ( 2014 ) Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach . Nat Commun 5:4006 26. Davnall F , Yip CS , Ljungqvist G et al ( 2012 ) Assessment of tumour heterogeneity: an emerging imaging tool for clinical practice ? Insights Imaging 3 : 573 - 589 27. O'Connor JP , Rose CJ , Waterton JC , Carano RA , Parker GJ , Jackson A ( 2015 ) Imaging intratumour heterogeneity: role in therapy response, resistance and clinical outcome . Clin Cancer Res 21 : 249 - 257 28. Ganeshan B , Panayiotou E , Burnand K , Dizdarevic S , Miles K ( 2012 ) Tumour heterogeneity in non-small cell lung carcinoma assessed by CT texture analysis: a potential marker of survival . Eur Radiol 22 : 796 - 802 29. Parikh J , Selmi M , Charles-Edwards G et al ( 2014 ) Changes in primary breast cancer heterogeneity may augment midtreatment MR imaging assessment of response to neoadjuvant chemotherapy . Radiology 272 : 100 - 112 30. Symmans WF , Peintinger F , Hatzis C et al ( 2007 ) Measurement of residual cancer burden to predict survival after neoadjuvant chemotherapy . J Clin Oncol 25 ( 28 ): 4414 - 4422 31. Scottish Intercollegiate Guidelines Network. Treatment of primary breast cancer . Edinburgh: Health Improvement Scotland ; 2013 . Available from: http://www.sign. ac.uk/pdf/SIGN134 .pdf. Accessed September 30th , 2015 . 32. Szczpinski PM , Strzelecki M , Materka A , Klepaczko A ( 2009 ) MaZda- a software package for image texture analysis . Comput Methods Programs Biomed 94 ( 1 ): 66 - 76 33. Haralick RM , Shanmugam K , Dinstein I ( 1973 ) Textural features for image classification . IEEE Trans on Syst Man Cybern 3 : 610 - 621 34. Mougalian SS , Hernandez M , Lei X et al ( 2016 ) Ten-year outcomes of patients with breast cancer with cytologically confirmed axillary lymph node metastases and pathological complete response after primary systemic therapy . JAMA Oncol 2 : 508 - 516 35. Pickles MD , Lowry M , Manton DJ , Turnbull LW ( 2015 ) Prognostic value of DCE-MRI in breast cancer patients undergoing neoadjuvant chemotherapy: a comparison with traditional survival indicators . Eur Radiol 25 : 1097 - 1106 36. Chen J-H , Nie K , Bahri S et al ( 2010 ) Decrease in breast density in the contralateral normal breast of patients receiving neoadjuvant chemotherapy: MR imaging evaluation . Radiology 255 : 44 - 52 37. O'Sullivan TD , Leproux A , Chen J-H et al ( 2013 ) Optical imaging correlates with magnetic resonance imaging breast density and reveals composition changes during neoadjuvant chemotherapy . Breast Cancer Res 15 : R14 38. Carballido-Gamio J , Link TM , Majumder S ( 2008 ) New Techniques for cartilage magnetic resonance imaging relaxation time analysis: texture analysis of flattened cartilage and localised intra- and inter- subject comparisons . Magn Reson Med 59 : 1472 - 1477 39. Mathias JM , Tofts PS , Losseff NA ( 1999 ) Texture analysis of spinal cord pathology in multiple sclerosis . Magn Reson Med 42 : 929 - 935 40. Teruel JR , Heldahl MG , Goa PE et al ( 2014 ) Dynamic contrastenhanced MRI texture analysis for pre-treatment prediction of clinical and pathological response to neoadjuvant chemotherapy in patients with locally advanced breast cancer . NMR Biomed 27 : 887 - 896 41. Michoux N , Van den Broeck S , Lacoste L et al ( 2015 ) Texture analysis on MR images help predicting non-response to NAC in breast cancer . BMC Cancer 15 : 574 - 586 42. Ahmed A , Gibbs P , Pickles M , Turnbull L ( 2013 ) Texture analysis in assessment and prediction of chemotherapy response in breast cancer . J Magn Reson Imaging 38 : 89 - 101 43. Flynn E , Collins D, D'Arcy J , Schmidt M , de Souza N ( 2016 ) Multi-parametric MRI in the early prediction of response to neoadjuvant chemotherapy in breast cancer: value of non-modelled parameters . Eur J Radiol 85 : 837 - 842 44. Li W , Arasu V , Newitt D et al ( 2016 ) Effect of MR imaging contrast thresholds on prediction of neoadjuvant chemotherapy response in breast cancer subtypes: a subgroup analysis of the ACRIN 6657/ISPY 1 trial . Tomography 2: 378 - 387 45. Ring A , Webb A , Ashely S et al ( 2003 ) Is surgery necessary after complete clinical remission following neoadjuvant chemotherapy for early breast cancer ? J Clin Oncol 21 : 4540 - 4545 46. Waugh SA , Lerski RA , Bidaut L , Thompson AM ( 2011 ) The influence of field strengths and different clinical breast MRI protocols on the outcome of texture analysis using foam phantoms . Med Phys 38 : 5058 - 5066


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1007%2Fs00330-017-4850-8.pdf

Shelley Henderson, Colin Purdie, Caroline Michie, Andrew Evans, Richard Lerski, Marilyn Johnston, Sarah Vinnicombe, Alastair M. Thompson. Interim heterogeneity changes measured using entropy texture features on T2-weighted MRI at 3.0 T are associated with pathological response to neoadjuvant chemotherapy in primary breast cancer, European Radiology, 2017, 1-10, DOI: 10.1007/s00330-017-4850-8