Assessment of early treatment response to neoadjuvant chemotherapy in breast cancer using non-mono-exponential diffusion models: a feasibility study comparing the baseline and mid-treatment MRI examinations

European Radiology, Nov 2016

Objectives To assess the feasibility of the mono-exponential, bi-exponential and stretched-exponential models in evaluating response of breast tumours to neoadjuvant chemotherapy (NACT) at 3 T. Methods Thirty-six female patients (median age 53, range 32–75 years) with invasive breast cancer undergoing NACT were enrolled for diffusion-weighted MRI (DW-MRI) prior to the start of treatment. For assessment of early response, changes in parameters were evaluated on mid-treatment MRI in 22 patients. DW-MRI was performed using eight b values (0, 30, 60, 90, 120, 300, 600, 900 s/mm2). Apparent diffusion coefficient (ADC), tissue diffusion coefficient (D t), vascular fraction (ƒ), distributed diffusion coefficient (DDC) and alpha (α) parameters were derived. Then t tests compared the baseline and changes in parameters between response groups. Repeatability was assessed at inter- and intraobserver levels. Results All patients underwent baseline MRI whereas 22 lesions were available at mid-treatment. At pretreatment, mean diffusion coefficients demonstrated significant differences between groups (p < 0.05). At mid-treatment, percentage increase in ADC and DDC showed significant differences between responders (49 % and 43 %) and non-responders (21 % and 32 %) (p = 0.03, p = 0.04). Overall, stretched-exponential parameters showed excellent repeatability. Conclusion DW-MRI is sensitive to baseline and early treatment changes in breast cancer using non-mono-exponential models, and the stretched-exponential model can potentially monitor such changes. Key points • Baseline diffusion coefficients demonstrated significant differences between complete pathological responders and non-responders. • Increase in ADC and DDC at mid-treatment can discriminate responders and non-responders. • The ƒ fraction at mid-treatment decreased in responders whereas increased in non-responders. • The mono- and stretched-exponential models showed excellent inter- and intrarater repeatability. • Treatment effects can potentially be assessed by non-mono-exponential diffusion models.

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-016-4630-x.pdf

Assessment of early treatment response to neoadjuvant chemotherapy in breast cancer using non-mono-exponential diffusion models: a feasibility study comparing the baseline and mid-treatment MRI examinations

Assessment of early treatment response to neoadjuvant chemotherapy in breast cancer using non-mono-exponential diffusion models: a feasibility study comparing the baseline and mid-treatment MRI examinations Reem Bedair 0 1 2 Andrew N. Priest 0 1 2 Andrew J. Patterson 0 1 2 Mary A. McLean 0 1 2 Martin J. Graves 0 1 2 Roido Manavaki 0 1 2 Andrew B. Gill 0 1 2 Oshaani Abeyakoon 0 1 2 John R. Griffiths 0 1 2 Fiona J. Gilbert 0 1 2 0 Cancer Research UK Cambridge Institute, University of Cambridge, Li Ka Shing Centre , Cambridge CB2 0RE , UK 1 Department of Radiology, Addenbrookes Hospital, Cambridge University Hospitals NHS Foundation Trust , Hills Road, Cambridge CB2 0QQ , UK 2 Department of Radiology, School of Clinical Medicine, University of Cambridge , Box 218, Cambridge Biomedical Campus, Hills Road, Cambridge CB2 0QQ , UK Objectives To assess the feasibility of the mono-exponential, bi-exponential and stretched-exponential models in evaluating response of breast tumours to neoadjuvant chemotherapy (NACT) at 3 T. Methods Thirty-six female patients (median age 53, range 32-75 years) with invasive breast cancer undergoing NACT were enrolled for diffusion-weighted MRI (DW-MRI) prior to the start of treatment. For assessment of early response, changes in parameters were evaluated on mid-treatment MRI in 22 patients. DW-MRI was performed using eight b values (0, 30, 60, 90, 120, 300, 600, 900 s/mm2). Apparent diffusion coefficient (ADC), tissue diffusion coefficient (Dt), vascular fraction (ƒ), distributed diffusion coefficient (DDC) and alpha (α) parameters were derived. Then t tests compared the baseline and changes in parameters between response groups. Repeatability was assessed at inter- and intraobserver levels. Breast carcinoma; Diffusion-weighted MRI; Neoadjuvant treatment; Quantitative evaluation; Tumour biomarkers - Results All patients underwent baseline MRI whereas 22 lesions were available at mid-treatment. At pretreatment, mean diffusion coefficients demonstrated significant differences between groups (p < 0.05). At mid-treatment, percentage increase in ADC and DDC showed significant differences between responders (49 % and 43 %) and non-responders (21 % and 32 %) (p = 0.03, p = 0.04). Overall, stretched-exponential parameters showed excellent repeatability. Conclusion DW-MRI is sensitive to baseline and early treatment changes in breast cancer using non-mono-exponential models, and the stretched-exponential model can potentially monitor such changes. Key points Baseline diffusion coefficients demonstrated significant differences between complete pathological responders and non-responders. Increase in ADC and DDC at mid-treatment can discriminate responders and non-responders. The ƒ fraction at mid-treatment decreased in responders whereas increased in non-responders. The mono- and stretched-exponential models showed excellent inter- and intrarater repeatability. Treatment effects can potentially be assessed by non-monoexponential diffusion models. Abbreviations α Alpha, Intravoxel heterogeneity index ADC Apparent diffusion coefficient Neoadjuvant chemotherapy (NACT) has increasingly been utilized in the treatment of breast cancer to enable breastc o n s e r v i n g s u r g e r y a n d i m p r o v e r e s e c t a b i l i t y. Approximately 80 % of patients have been found to respond to NACT, but only 6–25 % of patients show complete pathological response (pCR) [1–3]. Therefore, functional imaging techniques have been investigated for the prediction of response early after initiating therapy. Diffusion-weighted magnetic resonance imaging (DWMRI) derives image contrast from the differences in the mobility of water protons between tissues. DW imaging in breast protocols generally involves the acquisition of images at two b values to quantitatively determine the mono-exponential relationship between signal attenuation and b value. The apparent diffusion coefficient (ADC) enables characterisation of lesions based on differences in water diffusivity which in turn reflect tissue cellularity and integrity of membranes [4]. The association with cellular density makes ADC ideal for monitoring cytotoxic responses [5]. However, it has been shown that the degree of signal attenuation in breast tissue becomes non-linear with increasing b value. Bi-exponential signal decay has been observed over a 2 range of b values, where a small increase (0 − 200 s/mm ) results in a steep reduction in the measured signal intensity. This has been related to perfusion in the microcapillary circulation (expressed as pseudo-diffusion coefficient Dp and vascular fraction ƒ). The signal then attenuates more gradually over the range of higher b values (>200 s/m2) enabling the measurement of true tissue diffusivity (Dt). This phenomenon, known as intravoxel incoherent motion (IVIM), enables the separation of molecular diffusion from perfusion, provided that a wide range of low and high b values are used [6]. Whilst the IVIM model has the advantage that its components can be associated with distinct physical phenomena [7, 8], estimates of perfusion have not been fully investigated in practice for the assessment of therapeutic response in breast cancer. To overcome the assumptions associated with the biexponential model, Bennett et al. introduced the stretchedexponential approach, which models the continuous distribution of diffusion compartments attenuating at different rates (termed DDC, distributed diffusion coefficient). The plot of the signal intensity vs. b value becomes characteristically stretched, indicating deviation from the single exponential decay. This is denoted by the parameter α, which provides a new type of image contrast that relates to the degree of intravoxel water diffusion heterogeneity ranging from 0 to 1. A numerically high α index (i.e. α approaching 1) represents low intravoxel heterogeneity indicative of mono-exponential diffusion-weighted signal decay, whereas a numerically low α index (i.e. approaching 0) represents a high degree of diffusion heterogeneity exhibited as multi-exponential signal decay [9–11]. In general, these non-Gaussian diffusion models offer more parameters, which provide a better fit to the diffusion data. Preliminary studies suggest that multi-exponential models implemented in tumours of the brain [12, 13], head and neck [14, 15], abdomen [16] and prostate [17–19] can offer additional information on tissue heterogeneity, vascularity and cellularity beyond ADC. The purpose of this work was to assess the feasibility of diffusion parameters obtained from the mono-exponential, biexponential and stretched-exponential diffusion models in evaluating response of breast tumours to NACT at 3 T. Materials and methods Patient population and study design The local institutional review boards and ethics committees approved this prospective study and written informed consent was obtained from all patients. Patients were eligible if they were at least 18 years of age, had pathologically confirmed invasive breast cancer and were undergoing NACT as a first line of treatment. Patients were ineligible if they had poor renal or liver functions, allergy to gadolinium-based contrast agents, metals implants or a pacemaker. Assuming a moderate effect size between the diagnostic performance of diffusion parameters and response (effect size = 0.6), a sample size of 19 patients would be needed to yield a power of 80 % with 95 % confidence levels. Between February 2014 and September 2015, 40 female patients (median age 55, range 32 − 75 years) presented with a palpable lump and underwent core biopsies under ultrasound guidance by an experienced radiologist in the outpatient clinic prior to the start of treatment. According to our local protocol, patients received six cycles of NACT. The regimen consisted of docetaxel 100 mg/m2 once every 21 days for three cycles, followed by fluorouracil 500 mg/m2, epirubicin 100 mg/m2, with cyclophosphamide 500 mg/m2 (FEC) once every 21 days for three cycles if the tumour was negative for human epidermal growth factor 2 (HER2−) on biopsy. Two patients, however, were started on weekly taxols for 12 weeks at a lower dose because of their age, as studies show that weekly doses may offer the same benefit with fewer side effects [20, 21]. They subsequently received the FEC part of their treatment as described. This combination was reversed with the HER2+ cancers. In addition, docetaxel was combined with the HER2targeted agent trastuzumab (Herceptin®; Genentech Inc., CA) in these patients for the last three cycles. This therapeutic combination is based on the recommendation of the international consensus conference for neoadjuvant systemic therapy in primary breast cancer [1]. Patient and tumour characteristics are summarized in Table 1. All MR examinations were performed on a 3.0-T system (MR750, GE Healthcare) with a dedicated eight-channel breast array coil. The MRI protocol included T1- and T2weighted axial images, DWI and dynamic contrast-enhanced (DCE) series (Table 2). DW-MRI was performed utilizing a single-shot spin-echo echo planar imaging sequence at eight b values (0, 30, 60, 90, 120, 300, 600 and 900 s/mm2). Subsequently, DCE-MRI data was acquired using a three-dimensional segmented k-space spoiled gradientecho technique. Five acquisitions were obtained before contrast agent injection and then once every 10 s after bolus injection of 0.1 mmol/kg of gadopentetate dimeglumine (Magnevist, Bayer Schering, Berlin) for 8 min 7 s. The total examination time was approximately 25 min. Patients underwent MR imaging before the start of chemotherapy, after completion of three cycles and at the end of therapy. For the purposes of assessing early response, data from the baseline and midway scans only were analysed. The baseline MRI was performed approximately 10 days after histological confirmation of malignancy (range 7–11 days). The median time interval between pretreatment MRI and the start of therapy was 1 day (range 0–2). The median interval between the third cycle of chemotherapy and the second MRI was 21 days (range 15–22). Tumour characteristics and patient outcome Mixed carcinomaa Unless otherwise indicated, data are number of patients a Histology showed invasive carcinoma of mixed types; one of which was mucinous and ductal type, and the other showing lobular growth pattern with tubule formation classified as mixed lobular and ductal carcinoma b Tumours were classified as oestrogen receptor positive (and progesterone receptor positive) if more than 10 % of the cells were stained positively c Tumours were classified as human epidermal growth factor 2 (HER2) positive when they scored 3+ at immunohistochemistry or when gene amplification was observed with fluorescence in situ hybridization (FISH) After completion of NACT, all patients received breast and axillary surgery, radiotherapy and endocrine therapy as appropriate. Tumours were identified on the post-contrast T1-weighted images by two breast radiologists in consensus (R.B, F.J.G, Cambridge, UK) with 4 and 20 years of experience in breast MRI. Lesion site, size, morphology (mass vs. non-mass-like lesion), enhancement pattern (heterogeneous vs. homogeneous enhancement) and kinetic features (curve type; type I (progressive), type II (plateau) or type III (washout)) were recorded. For both mass and non-mass-like lesion demarcation, the early subtracted Contrast-enhanced T1-weighted SE spin-echo, 2D SS-EPI 2-dimensional single-shot echo planar imaging, 3D SPGR 3-dimensional spoiled gradient-recalled echo, FOV field of view, ASSET array coil spatial sensitivity encoding technique, TR repetition time, TE echo time, RF radiofrequency DCE images were used as reference where freehand regions of interest (ROIs) were manually drawn on the imaging slice with the largest tumour dimension on the b = 900 s/mm2 image by one radiologist (R.B). Care was taken to avoid tumour borders and areas of necrosis. Intraobserver variability was evaluated by redrawing the ROIs 3 weeks after the initial measurements. A third radiologist (O.A) with 3 years’ experience in breast MRI manually reoutlined the tumours to assess the inter-rater variability of measurements. All readers were blinded to the pathological findings and therapeutic responses. Diffusion analyses were performed using in-house software developed in MATLAB (The Mathworks, Natick, MA). To enable comparison with previous DWI breast studies, the ADC was calculated from two b values (0, 900 s/ mm2). All b values were used for the stretched-exponential and IVIM models. For the non-Gaussian models, data was fitted using a non-linear least-squares approach. ROIs were analysed on a voxel-wise basis and parameters were expressed as means over the volumes measured. Parametric maps of diffusion coefficients were generated. Details of the quantitative diffusion models are provided in the supplementary material. Percentage change in parameters was calculated as: (Parametermid − Parameterpre)/Parameterpre × 100, where Parameterpre and Parametermid are the baseline and midtreatment measurements respectively. Histological analysis and response assessment Histopathological assessment was performed after surgical excision following the last cycle of chemotherapy. Tumour type, grade, hormone-receptor status and HER2 expression were obtained from reports of the core biopsies or surgical specimens. Although molecular testing is prognostic, it is expensive and not widely available. Therefore according to Onitilo et al. [22, 23], four immunohistochemical (IHC) categories have been identified on the basis of the hormonalreceptor status (oestrogen and progesterone receptor ER, PR) and HER2 overexpression. These groups are ER+PR+/ HER2−, ER+PR+/HER+, ER−PR−/HER2+ and ER−PR −/HER2− (i.e. triple-negative breast cancers, TNBC). Tumour response was assessed in the excision specimens by expert pathologists. Three categories were defined: (i) pathological complete response (pCR) with or without the presence of ductal carcinoma in situ (DCIS); (ii) partial response to therapy; ranging from minimal residual disease up to greater than 50 % of tumour cellularity; and (iii) no evidence of response. The last two categories were considered pathological non-complete response (pNCR). All statistics were calculated from the logarithm of the parameters, assessed for normality using the Shapiro–Wilk test and back-transformed where appropriate [24]. The mean value ± standard deviation of ADC and the non-Gaussian parameters were reported. The pretreatment characteristics between pCR and pNCR were compared using the unpaired two-tailed t test. The IHC subtypes were compared between response groups at baseline using the Chi-squared test. The imaging parameters across the four subtypes were also compared using the one-way analysis of variance (ANOVA). The correlation between the diffusion coefficients was evaluated using Pearson’s correlation. Receiver operating characteristics (ROC) curves were constructed to assess the performance of parameters in differentiating between pCR and pNCR, and areas under the curve (AUCs) were compared. Repeatability of the baseline measurements was assessed at inter- and intraobserver levels using the intraclass correlation coefficient (ICC). Model comparisons were made using the corrected Akaike information criterion (AICc), which imposes a penalty for additional parameters in the model [25]. The AICc was computed for the three models and averaged over all lesion voxels, which was subsequently averaged over all subjects. Statistical analyses were performed using the software SPSS (v. 21.0, Chicago, IL). As this study is primarily descriptive, p values are presented as raw values and not corrected for multiple comparisons. For purposes of discussion and similar to Orton et al. [5], multiple comparisons were accounted for by adjusting the p value significance threshold by a correction factor of 5. Thus statistically significant comparisons were set at p < 0.05 (no correction) and highly significant comparisons at p < 0.01. Four patients opted out of the study prior to imaging, therefore baseline analysis included 36 patients with unifocal malignant lesions. Twenty-five lesions (69 %) were classified as mass enhancement while 11 were non-mass-like lesions (31 %). All lesions showed heterogeneous enhancement on DCE-MRI; 27 lesions (75 %) demonstrated a predominately type III enhancement curve while the rest (9/36, 25 %) showed type II pattern. After surgery, 14 patients (39 %) showed pCR and 22 (61 %) were pNCR. Prediction of response: pretreatment tumour characteristics The mean tumour diameter was 4.7 cm (median 4.0 cm, range 1.2 − 12 cm). No significant difference was found between responders (4.1 ± 0.4 cm) and non-responders (5.1 ± 0.5 cm, p = 0.194). According to IHC subtype, 33 % (12/36) were ER+PR+/ HER−, 33 % ER+PR+/HER2+, 3 % (1/36) ER−PR−/HER2+, and 31 % (11/36) were categorised as TNBCs. All but one patient in the triple-negative group showed pCR, whereas 17 % (2/12) of both ER+PR+/HER2− and ER+PR+/HER2+ tumours showed complete response. Significant differences were observed between response groups with respect to histological subtype (p < 0.001). Prediction of response: pretreatment diffusion parameters Pretreatment mean ADC, DDC and Dt values across the cohort were 1.11 × 10−3, 1.09 × 10−3 and 0.93 × 10−3 mm2/s respectively. Figures 1 and 2 show representative images and parametric maps of a complete responder. An example of the mono-, bi- and stretched-exponential curves fitted to the data from one pixel is illustrated in Fig. 3. Table 3 shows the mean diffusion coefficients of the response groups from all three models. Responders showed significantly lower ADC, DDC and Dt values (0.92 ± 0.03, 0.93 ± 0.04 and 0.85 ± 0.05 × 10−3mm2/s respectively) relative to non-responders (1.20 ± 0.02, 1.25 ± 0.03 and 1.02 ± 0.05 × 10−3mm2/s respectively) (p < 0.01, p < 0.01, p = 0.02 respectively). The ƒ fraction showed higher fractions in responders; however, this was not significant (p = 0.09). On subgroup analysis, ƒ was found to be significantly higher in responders of the TNBC subtype (12.4 ± 4.1 % vs. 10.9 ± 1.2 %, p = 0.01). Although lower α values were found in non-responders compared to eventual responders, this was not significant (0.81 ± 0.02 vs. 0.84 ± 0.02, p = 0.07). Figure 4 shows the ROC curves comparing the diagnostic efficacy of diffusion coefficients relative to pathological response. DDC demonstrated a larger AUC (0.756, p = 0.01) compared to ADC and Dt (0.749, p = 0.01 and 0.641, p = 0.15 respectively). The DDC cut-off to differentiate response groups on pretreatment MRI (1.141 × 10−3 mm2/s) yielded the highest sensitivity (81 %) and specificity (72 %) (Table 4). Mid-treatment assessment All patients attended a second MRI midway through treatment; however, 14 patients were excluded for the following reasons: patients opted out of the research examination (n = 3), images were unusable as a result of suboptimal fat suppression (n = 2), subjects not receiving the standard six cycles of chemotherapy (n = 2), three patients did not proceed with the full diffusion imaging sequence prior to surgery after NACT as they found it intolerable. Additionally, four patients were excluded from analysis as their follow-up MRI was not Fig. 1 Representative images from pretreatment MRI of a 49year-old female patient with cancer of the right breast: a axial DW image showing hyperintense tumour with restricted diffusion on the b900 s/mm2 image. b ADC map was generated from two b values (0, 900 s/mm2). ROI was drawn on the primary lesion and copied to the ADC map (ADC = 0.92 ± 0.094 × 10−3 mm2/s) performed in a timely manner at mid-therapy (post cycle 3). These patients developed liver and kidney lesions that were later found to be benign. Finally, 22 malignant lesions were available for analysis on the midway examination. Eight patients (36 %) showed pCR and 14 (64 %) were pNCRs. Prediction of response: change in tumour size A significant difference was found in tumour size as responders showed a smaller mean tumour size (1.5 ± 0.2 cm vs. 2.9 ± 0.5 cm, p < 0.05). Fig. 2 Parametric maps of the biexponential and stretchedexponential models in the same female patient as in Fig. 1: a Dt = 0.80 ± 0.28 × 10−3 mm2/s, b ƒ = 11.8 ± 1.3 %, c DDC = 0.98 ± 0.12 × 10−3 mm2/s, d α = 0.84 ± 0.18. It should be noted that the α values in tumours express the intravoxel heterogeneity, whereas the other maps show intervoxel heterogeneity between tissues. At histopathology, the lesion was identified as a grade 3 invasive ductal carcinoma of TNBC subtype. The patient underwent wide local excision of the lesion. On the excision specimen, no invasive components were seen and the patient was considered a complete pathological responder to NACT Fig. 3 Mono-exponential (green), bi-exponential (blue) and stretchedexponential curves (red) fitted to one pixel in a breast lesion in a pretreatment scan Prediction of response: change in diffusion parameters An increase in the mean values of the diffusion coefficients was observed after three cycles of chemotherapy, with the percentage increase in ADC and DDC showing a statistically significant difference between responders (49 % and 43 %) and non-responders (21 %, 32 %, p = 0.03 and p = 0.04 respectively). However the increase in Dt did not show a significant difference between response groups (36 % vs. 23 %, p = 0.14). Moreover, the decrease in ƒ fraction found in responders (29 %) was substantially different from the increase in ƒ Fig. 4 ROCs for the response prediction of the pretreatment diffusion coefficients from the mono-exponential, bi-exponential and stretchedexponential models. The DDC demonstrated the largest AUC (0.75, p = 0.01) compared with ADC and Dt (0.74, p = 0.01 and 0.641, p = 0.15 respectively). The cut-off to differentiate between response groups on pretreatment MRI for DDC (1.141 × 10−3 mm2/s) yielded the highest measures of accuracy (sensitivity 81 %, specificity 72 %) observed in pNCR (5 %, p = 0.05). Responders also showed a larger increase in α compared to pNCR (7 % vs. 5 %). This, however, was not significant (p = 0.68) (Table 3). Figure 5 Baseline (n = 36) ADC (×10−3mm2/s) DDC (×10−3mm2/s) Dt (×10−3mm2/s) ADC (×10−3mm2/s) DDC (×10−3mm2/s) Dt (×10−3mm2/s) 1.52 ± 0.32 (↑49 %) 1.51 ± 0.15 (↑43 %) 0.91 ± 0.07 (↑7 %) 1.30 ± 0.14 (↑36 %) 8.48 ± 1.54 (↓29 %) 1.27 ± 0.18 (↑21 %) 1.40 ± 0.12 (↑32 %) 0.86 ± 0.11 (↑5 %) 1.28 ± 0.15 (↑23 %) 10.53 ± 2.51 (↑5 %) Responders (% change) Non-responders (% change) p value of %change Baseline and mid-treatment values of the various parameters (units × 10−3 mm2 /s except where shown). Values in parenthesis are the percentage change. P values using a two-tailed independent t test compared between response groups. Significant statistics have p < 0.05, while highly significant statistics have p < 0.01 which includes a correction for multiple comparisons). Unless otherwise indicated, data represent mean values ± standard deviation (percentage change between pretreatment and mid-treatment values) of the mono-exponential (ADC), biexponential (Dt, ƒ) and stretched-exponential (DDC, α) parameters *Values show statistically significant differences **Values show highly significant differences Table 4 ROC analysis of the various parameters (units × 10−3 mm2/s except where shown) indicating the sensitivity and specificity measures at the respective cut-off values and their significance levels Cut-off value (×10−3mm2/s) Values show statistical significance of the pretreatment mono-exponential (ADC), bi-exponential (Dt, ƒ) and stretched-exponential (DDC, α) parameters AUC area under the curve shows the mean change in the diffusion coefficients between response groups. Repeatability of measurements Overall parameters of the mono- and stretched-exponential models showed excellent inter- and intrarater repeatability. Bi-exponential parameters ranged from excellent to fair (Table 5). A significant correlation was found between ADC and the non-Gaussian diffusion coefficients (p < 0.001 for both). A higher correlation, however, was observed between ADC and DDC (r = 0.89) compared to ADC and Dt (r = 0.76). The AICc from the ROI data showed that stretchedexponential was the preferred model at baseline and midtreatment measurements, showing lower values (62.6 ± 8.2, 66.6 ± 9.2) compared to the mono-exponential and biexponential models (baseline, 81.7 ± 9.4 and 82.9 ± 9.7; mid-treatment, 85.6 ± 10.2 and 88.6 ± 11.6 respectively). It is generally accepted that highly proliferating malignant breast lesions result in a packed cellular microstructure, showing more restricted diffusion and decreased ADC [26]. However the choice of b value may affect the calculated measurement, which is influenced by multiple pools diffusing at different rates, confounding the assessment of tissue diffusivity [27]. Alternative diffusion models have been devised to account for the more complex non-Gaussian diffusion behaviour of biological tissues in vivo. Le Bihan et al. have shown that signal attenuation is not only a result of random microscopic motion of water molecules influenced by cell density but is also dependent on microperfusion within the voxel [6]. Fig. 5 Mean distribution of the diffusion coefficients of the three models before the start of treatment and after three cycles of chemotherapy in complete and non-complete responders. There is an overall increase in parameters in both groups; however, a larger increase is noted in the ADC and DDC of complete responders. When the increase in mean values of ADC and DDC were compared between response groups at midtreatment, a significant difference was observed; p = 0.03 and p = 0.04 respectively. However, the increase in Dt did not show statistical significance between response groups (p = 0.14). Error bars represent the 95 % confidence interval Table 5 Inter- and intraobserver repeatability measures of the diffusion parameters using the intraclass correlation (ICC) metrics Interobserver agreement Intraobserver agreement ADC (×10−3mm2/s) 0.816 (0.657, 0.912) DDC (×10−3mm2/s) 0.789 (0.601, 0.898) Dt (×10−3mm2/s) 0.699 (0.641, 0.859) α (unit-less) 0.808 (0.760, 0.938) ƒ (%) 0.605 (0.569, 0.868) 0.910 (0.823, 0.954) 0.860 (0.725, 0.929) 0.778 (0.756, 0.889) 0.822 (0.648, 0.910) 0.695 (0.531, 0.893) Data represents intraclass correlation with 95 % confidence intervals This has become of particular interest, as tumour angiogenesis is seen an important determinant in the outcome of patients. Numerically, however, it has been shown that fitting three parameters for the bi-exponential model can be unstable [28, 29]. In this study we assessed two non-Gaussian models acquired at eight b values to capture the various diffusion properties of tissues and compared them with the monoexponential model using the standard b values commonly implemented in DW-MRI of the breast. Our baseline results were consistent with Sigmund et al. showing lower Dt values compared to ADC [30]. This was expected on exclusion of the perfusion effect, as the lower b values were employed to capture the much higher pseudodiffusion. Furthermore, a strong correlation was found between the mono- and stretched-exponential diffusion coefficients. This suggests that DDC can be interpreted in the same way as ADC with the observation of a continuous distribution of diffusion components within the microenvironment. Similar to Park et al. [31], our study also shows an inverse correlation between tumour diffusion coefficients and therapeutic response, where substantially higher pretreatment values of ADC, DDC, Dt and lower ƒ fractions were more suggestive of necrotic, less viable lesions. These lesions have often been found to be hypoxic, more aggressive and less sensitive to chemotherapy [32]. The mean ƒ fraction of tumours was also able to differentiate responders in the most biologically aggressive subtype (TNBC), showing complete response in about one-third of tumours, consistent with previous reports [33]. However, a full statistical evaluation was not possible because of the small sample size of the other subtypes. Although the attribution of the diffusion index α to a biological correlate is still under investigation, lower values seen in non-responsive lesions indicate a more heterogeneous microenvironment within the imaging voxels [19]. This parameter could in turn be viewed as a reflection of the microstructural complexity of the tissue, suggesting changes in the degree of cellular pleomorphism, vascular heterogeneity and presence of microscopic necrosis [13]. Although non-significant (p = 0.68), responders in our study showed an increase in the structural homogeneity (α closer to 1) at mid-treatment compared to non-responders. These preliminary results are similar to those found by Orton et al. in abdominal and pelvic tumours where changes in α were substantially different between response groups when assessed before the start of therapy, 7 days and 28 following treatment [5]. Neoadjuvant chemotherapy regimens used in the treatment of breast cancer commonly consist of an anthracycline in combination or sequentially administered with taxanes. The mechanism of action of anthracyclines is mainly attributed to the inhibition of DNA synthesis, preventing the replication of rapidly dividing cells [34], whereas taxanes have been reported to have an anti-angiogenic effect with selective shutdown of microvessel formation [35, 36]. This is supported by the significant increase in DDC seen on mid-treatment (p = 0.04), suggesting a change in the distribution of diffusion compartments resulting from the early breakdown of the vascular endothelial cells. On IVIM analysis, we also showed a large decrease in the perfusion fraction in lesions showing complete response, which is in contrast to its increase in lesions failing to respond to treatment. This result is supported by previous findings [37, 38]. When repeatability of measurements was assessed, the mono- and stretched-exponential showed high ICC measures, suggesting equivalent robustness of the derived parameters to ADC. This is consistent with recent studies that have found the stretched exponential equivalent to or outperforming other models [5, 12, 16, 39]. This work suffered a few limitations. First, the chemotherapy regimen in our population varied, which made it difficult to draw specific conclusions on the biological effect of each drug as detected by the diffusion models. However, this should not have impacted the final response assessment as the population received the same combination by the end of NACT. Second, although we evaluated the repeatability of the ROI measurements, we did not perform validation studies. Orton et al., however, showed that the parameters derived from the stretched-exponential model are highly reproducible and could potentially serve as quantitative biomarkers for response assessment in abdominal and pelvic tumours [5]. In this study, we implemented the stretched-exponential model using eight b values to characterise the deviation from the Gaussian curve and compared it with the bi-exponential model. Past studies have used 4−5 b values and produced similar results [13, 14, 19], which affords the potential of increased DW-MRI capability for characterising tissue properties over an extended b value range at reasonable scan times. Finally, our study was limited by the small sample size with differing rates of responders within the tumour subtypes. However, we showed that the DW-MRI could be sensitive to baseline and early treatment effects caused by NACT using the DDC parameter of the stretched-exponential model, as demonstrated by the large AUC with satisfactory sensitivity and specificity measures, and the ƒ fraction of the biexponential model. Unlike in abdominal and pelvic tumours [5], α did not differentiate between response groups in the breast. This may relate to the timing of the follow-up MRI, which was performed after the third cycle of NACT (i.e. 6 weeks of treatment). Therefore follow-up imaging at an earlier time point may elucidate the change in α between responders and non-responders. In conclusion, this feasibility study showed that DW-MRI is sensitive to baseline and early treatment changes in breast cancer using non-mono-exponential models, which offer additional imaging biomarkers that can potentially provide insights into the cellular compartments and membranes and may become more sensitive to treatment-induced tissue changes. Our results show that the stretched-exponential model can potentially monitor such changes. This data supports the wider use of these models in assessing treatment effects beyond that routinely measured with ADC. Acknowledgments The project was supported by the Addenbrooke’s Charitable Trust and the NIHR comprehensive Biomedical Research Centre (BRC) and the Experimental Cancer Medicine Centre (ECMC) awards to Cambridge University Hospitals NHS Foundation Trust in partnership with the University of Cambridge. The scientific guarantor of this publication is Professor Fiona J. Gilbert. The authors of this manuscript declare no relationships to any companies whose product or services may be related to the subject matter of the article. The study has received funding from the Addenbrookes Charitable Trust and the NIHR comprehensive Biomedical Research Centre (BRC) and the Experimental Cancer Medicine Centre (ECMC) awards to Cambridge University Hospitals NHS Foundation Trust in partnership with the University of Cambridge. No complex statistical methods were necessary for this paper. Institutional review board approval was obtained. Written informed consent was obtained from all subjects (patients) in this study. Methodology: Prospective, diagnostic or prognostic 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. Kaufmann M , von Minckwitz G , Mamounas E , Cameron D , Carey L , Cristofanilli M et al ( 2012 ) Recommendations from an international consensus conference on the current status and future of neoadjuvant systemic therapy in primary breast cancer . Ann Surg Oncol 19 : 1508 - 1516 2. Bear HD ( 2006 ) Sequential preoperative or postoperative docetaxel added to preoperative doxorubicin plus cyclophosphamide for operable breast cancer: National Surgical Adjuvant Breast and Bowel Project Protocol B-27. J Clin Oncol 24 : 2019 - 2027 3. Cho N , Im S-A , Park I-A , Lee K-H , Li M , Han W et al ( 2014 ) Breast cancer: early prediction of response to neoadjuvant chemotherapy using parametric response maps for MR imaging . Radiology 272 : 385 - 396 4. Woodhams R , Ramadan S , Inoue Y ( 2011 ) Diffusion-weighted imaging of the breast: principles and clinical. Radiographics 31 : 1059 - 1084 5. Orton MR , Messiou C , Collins D , Morgan VA , Tessier J , Young H et al ( 2016 ) Diffusion-weighted MR imaging of metastatic abdominal and pelvic tumours is sensitive to early changes induced by a VEGF inhibitor using alternative diffusion attenuation models . Eur Radiol 26 : 1412 - 1419 6. Le Bihan D , Breton E , Lallemand D , Grenier P , Cabanis E , LavalJeantet M ( 1986 ) MR imaging of intravoxel incoherent motions: application to diffusion and perfusion in neurologic disorders . Radiology 161 : 401 - 407 7. Iacconi C , Giannelli M , Marini C , Cilotti A , Moretti M , Viacava P et al ( 2010 ) The role of mean diffusivity (MD) as a predictive index of the response to chemotherapy in locally advanced breast cancer: a preliminary study . Eur Radiol 20 : 303 - 308 8. Li X-R , Cheng L-Q , Liu M , Zhang Y-J , Wang J-D , Zhang A-L et al ( 2012 ) DW-MRI ADC values can predict treatment response in patients with locally advanced breast cancer undergoing neoadjuvant chemotherapy . Med Oncol 29 : 425 - 431 9. Bennett KM , Schmainda KM , Bennett RT , Rowe DB , Lu H , Hyde JS ( 2003 ) Characterization of continuously distributed cortical water diffusion rates with a stretched-exponential model . Magn Reson Med 50 : 727 - 734 10. Bennett KM , Hyde JS , Rand SD , Bennett R , Krouwer HGJ , Rebro KJ et al ( 2004 ) Intravoxel distribution of DWI decay rates reveals C6 glioma invasion in rat brain . Magn Reson Med 52 : 994 - 1004 11. Bennett KM , Hyde JS , Schmainda KM ( 2006 ) Water diffusion heterogeneity index in the human brain is insensitive to the orientation of applied magnetic field gradients . Magn Reson Med 56 : 235 - 239 12. Bai Y , Lin Y , Tian J , Shi D , Cheng J , Em H et al ( 2016 ) Grading of gliomas by using monoexponential, biexponential, and stretched exponential diffusion-weighted MR imaging and diffusion kurtosis MR imaging . Radiology 278 : 496 - 504 13. Kwee TC , Galbán CJ , Tsien C , Junck L , Sundgren PC , Ivancevic MK et al ( 2010 ) Comparison of apparent diffusion coefficients and distributed diffusion coefficients in high-grade gliomas . J Magn Reson Imaging 31 : 531 - 537 14. Lai V , Lee VHF , Lam KO , Sze HCK , Chan Q , Khong PL ( 2015 ) Intravoxel water diffusion heterogeneity MR imaging of nasopharyngeal carcinoma using stretched exponential diffusion model . Eur Radiol 25 : 1708 - 1713 15. Yuan J , Ka D , Yeung W , Mok GSP , Bhatia KS , Wang YJ et al ( 2014 ) Non-Gaussian analysis of diffusion weighted imaging in head and neck at 3T : a pilot study in patients with nasopharyngeal carcinoma . PLoS One 9 : e87024 16. Winfield JM , DeSouza NM , Priest AN , Wakefield JC , Hodgkin C , Freeman S et al ( 2015 ) Modelling DW-MRI data from primary and metastatic ovarian tumours . Eur Radiol 25 : 2033 - 2040 17. Mazaheri Y , Afaq A , Rowe DB , Lu Y , Shukla-Dave A , Grover J ( 2012 ) Diffusion-weighted magnetic resonance imaging of the prostate: improved robustness with stretched exponential modeling . J Comput Assist Tomogr 36 : 695 - 703 18. Jambor I , Merisaari H , Taimen P , Boström P , Minn H , Pesola M et al ( 2014 ) Evaluation of different mathematical models for diffusion-weighted imaging of normal prostate and prostate cancer using high b-values: a repeatability study . Magn Reson Med 1998 : 1988 - 1998 19. Liu X , Zhou L , Peng W , Wang H , Zhang Y ( 2015 ) Comparison of stretched-Exponential and monoexponential model diffusionweighted imaging in prostate cancer and normal tissues . J Magn Reson Imaging 42 : 1078 - 1085 20. Eniu A , Palmieri FM , Perez EA ( 2005 ) Weekly administration of docetaxel and paclitaxel in metastatic or advanced breast cancer . Oncologist 10 : 665 - 685 21. Budd G , Barlow W , Moore H , Hobday T , Stewart J , Isaacs C et al ( 2013 ) Comparison of two schedules of paclitaxel as adjuvant therapy for breast cancer . J Clin Oncol 31(suppl; abstr CRA1008) 22. Parise C , Caggiano V ( 2014 ) Breast cancer survival defined by the ER/PR/HER2 subtypes and a surrogate classification according to tumor grade and immunohistochemical biomarkers . J Cancer Epidemiol 2014 : 1 - 11 23. Onitilo AA , Engel JM , Greenlee RT , Mukesh BN ( 2009 ) Breast cancer subtypes based on ER/PR and Her2 expression: comparison of clinicopathologic features and survival . Clin Med Res 7 : 4 - 13 24. Keene ON ( 1995 ) The log transformation is special . Stat Med 14 : 811 - 819 25. Glatting G , Kletting P , Reske SN , Hohl K , Ring C ( 2007 ) Choosing the optimal fit function: comparison of the Akaike information criterion and the F-test . Med Phys 34 : 4285 - 4292 26. Pickles MD , Gibbs P , Lowry M , Turnbull LW ( 2006 ) Diffusion changes precede size reduction in neoadjuvant treatment of breast cancer . Magn Reson Imaging 24 : 843 - 847 27. Taouli B , Koh D-M ( 2010 ) Diffusion-weighted MR imaging of the liver . Radiology 254 : 47 - 66 28. Hoff B , Chenevert TL , Bhojani MS , Kwee TC , Rehemtulla A , Le Bihan D et al ( 2010 ) Assessment of multiexponential diffusion features as MRI cancer therapy response metrics . Magn Reson Med 64 : 1499 - 1509 29. Koh D-M , Collins DJ , Orton MR ( 2011 ) Intravoxel incoherent motion in body diffusion-weighted MRI: reality and challenges . Am J Roentgenol 196 : 1351 - 1361 30. Sigmund EE , Cho GY , Kim S , Finn M , Moccaldi M , Jensen JH et al ( 2011 ) Intravoxel incoherent motion imaging of tumor microenvironment in locally advanced breast cancer . Magn Reson Med 65 : 1437 - 1447 31. Park SH , Moon WK , Cho N , Song IC , Chang JM , Park I-A et al ( 2010 ) Diffusion-weighted MR imaging: pretreatment prediction of response to neoadjuvant chemotherapy in patients with breast cancer . Radiology 257 : 56 - 63 32. Li SP , Padhani AR , Makris A ( 2011 ) Dynamic contrast-enhanced magnetic resonance imaging and blood oxygenation leveldependent magnetic resonance imaging for the assessment of changes in tumor biology with treatment . JNCI Monogr 2011 : 103 - 107 33. Goldhirsch A , Wood WC , Coates AS , Gelber RD , Thurlimann B , Senn H-J ( 2011 ) Strategies for subtypes-dealing with the diversity of breast cancer: highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2011 . Ann Oncol 22 : 1736 - 1747 34. Nitiss J ( 2009 ) Targeting DNA topoisomerase II in cancer chemotherapy . Nat Rev Cancer 9 : 338 - 350 35. Sweeney CJ , Miller KD , Sissons SE , Nozaki S , Heilman DK , Shen J et al ( 2001 ) The antiangiogenic property of docetaxel is synergistic with a recombinant humanized monoclonal antibody against vascular endothelial growth factor or 2-methoxyestradiol but antagonized by endothelial growth factors . Cancer Res 61 : 3369 - 3372 36. Li S , Makris A , Taylor N , Beresford M , Ah-See M , Stirling J et al ( 2009 ) Greater anti-angiogenic effects with docetaxel compared with anthracyclines in the neoadjuvant treatment of breast cancer . Cancer Res 69 : 5013 37. Mankoff D , Dunnwald LK , Gralow JR , Ellis GK , Schubert EK , Tseng J et al ( 2003 ) Changes in blood flow and metabolism in locally advanced breast cancer treated with neoadjuvant chemotherapy . J Nucl Med 44 : 1806 - 1814 38. Semple S , Staff RT , Heys SD , Redpath TW , Welch AE , Ahearn TS et al ( 2006 ) Baseline MRI delivery characteristics predict change in invasive ductal breast carcinoma PET metabolism as a result of primary chemotherapy administration . Ann Oncol 17 : 1393 - 1398 39. Panek R , Borri M , Orton M , Flynn EO , Morgan V , Giles SL et al ( 2015 ) Evaluation of diffusion models in breast cancer . Med Phys 42 : 4833 - 4839


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1007%2Fs00330-016-4630-x.pdf

Reem Bedair, Andrew N. Priest, Andrew J. Patterson, Mary A. McLean, Martin J. Graves, Roido Manavaki, Andrew B. Gill, Oshaani Abeyakoon, John R. Griffiths, Fiona J. Gilbert. Assessment of early treatment response to neoadjuvant chemotherapy in breast cancer using non-mono-exponential diffusion models: a feasibility study comparing the baseline and mid-treatment MRI examinations, European Radiology, 2017, 2726-2736, DOI: 10.1007/s00330-016-4630-x