- 1Institute of Medical Sciences, School of Medicine, University of Aberdeen, Aberdeen, United Kingdom
- 2Department of Oncology, Aberdeen Royal Infirmary, Aberdeen, United Kingdom
- 3Department of Radiology, Royal Marsden Hospital, London, United Kingdom
- 4Department of Pathology, Aberdeen Royal Infirmary, Aberdeen, United Kingdom
- 5Breast Unit, Aberdeen Royal Infirmary, Aberdeen, United Kingdom
- 6Translational and Clinical Research Institute, Faculty of Medical Sciences, Newcastle University, Newcastle upon Tyne, United Kingdom
Introduction: The early identification of good responders to neoadjuvant chemotherapy (NACT) holds a significant potential in the optimal treatment of breast cancer. A recent Bayesian approach has been postulated to improve the accuracy of the intravoxel incoherent motion (IVIM) model for clinical translation. This study examined the prediction and early sensitivity of Bayesian IVIM to NACT response.
Materials and methods: Seventeen female patients with breast cancer were scanned at baseline and 16 patients were scanned after Cycle 1. Tissue diffusion and perfusion from Bayesian IVIM were calculated at baseline with percentage change at Cycle 1 computed with reference to baseline. Cellular proliferative activity marker Ki-67 was obtained semi-quantitatively with percentage change at excision computed with reference to core biopsy.
Results: The perfusion fraction showed a significant difference (p = 0.042) in percentage change between responder groups at Cycle 1, with a decrease in good responders [−7.98% (−19.47–1.73), n = 7] and an increase in poor responders [10.04% (5.09–28.93), n = 9]. There was a significant correlation between percentage change in perfusion fraction and percentage change in Ki-67 (p = 0.042). Tissue diffusion and pseudodiffusion showed no significant difference in percentage change between groups at Cycle 1, nor was there a significant correlation against percentage change in Ki-67. Perfusion fraction, tissue diffusion, and pseudodiffusion showed no significant difference between groups at baseline, nor was there a significant correlation against Ki-67 from core biopsy.
Conclusion: The alteration in tumour perfusion fraction from the Bayesian IVIM model, in association with cellular proliferation, showed early sensitivity to good responders in NACT.
Clinical trial registration: https://clinicaltrials.gov/ct2/show/NCT03501394, identifier NCT03501394.
1 Introduction
Neoadjuvant chemotherapy (NACT) is increasingly used in breast cancer, evolving from originally downstaging inoperable breast tumours to allow surgical excision (1) to facilitating potential breast and axillae conservation (2). However, NACT not only is costly at an estimated £6,000 per patient for a typical six-cycle regimen of 5-fluorouracil/epirubicin/cyclophosphamide (FEC 100) in the National Health Service (3) but also often leads to adverse side effects and subsequent severe physical and emotional distress (4, 5). Although NACT improves rates of pathological complete response (pCR) (6, 7) and disease-free survival (7, 8), poor responders to NACT might receive earlier and timely mastectomy or breast conservation (9). RECIST criterion, the current approach to estimate residual disease load based on tumour size (10) at the halfway point of NACT (11), has limited accuracy at a relatively late stage of treatment, demanding more precise radiological approaches.
The loss of tumour cellularity is the central histological marker of cellular damage in tumours responding to NACT (12). Diffusion-weighted imaging (DWI), although sensitive to cellularity (13, 14) with the potential of identifying responders after one cycle of NACT (15), is susceptible to biological noise and limited to large cohort studies (16), and is therefore inadequate for response-guided NACT (17). Apparent diffusion coefficient (ADC) from DWI (18) is effective in differentiation of tumour from healthy tissue and benign lesions (19, 20). An increase in ADC at the halfway point of 12 weeks of NACT-predicted pCR, however, may not reach clinical relevance with the receiver operating characteristics curve at an area under 0.6 (21). Diffusion tensor imaging yielded a significant increase in prime diffusion coefficient (λ1) and ADC in good responders compared to poor responders at the completion of NACT, although baseline diffusion metrics did not predict good response (22). Diffusion kurtosis imaging approximates the deviation from the tensor model using kurtosis, with a lower mean kurtosis at baseline associated with pCR at four cycles of NACT in patients with breast cancer (23). We have shown that q-space imaging was more effective in the evaluation of cellularity in breast cancer; however, the method was not suitable for routine clinical application due to the demand on high field gradient and long scan duration (24). Intravoxel incoherent motion (IVIM), incorporating tissue diffusion and blood microcirculation as two independent components (18), showed improved diagnostic sensitivity in breast cancer (25). However, IVIM is prone to misfitting as a result of the high susceptibility in the main algorithm to biological noise (26). A Bayesian probability (BP) approach has been suggested to improve fitting accuracy and reduce variability in the estimation of tissue diffusion and blood microcirculation (27).
We therefore hypothesise that the Bayesian IVIM model may differentiate good from poor responders at baseline and after Cycle 1 of NACT with association from tumour proliferative activity, providing a non-invasive biomarker sensitive to prediction and early response to NACT.
2 Materials and methods
We hence conducted a prospective, longitudinal study of NACT in 17 female patients with breast cancer using the Bayesian IVIM model (Figure 1). The study was approved by the London Research Ethics Committee (Identifier: 17/LO/1777) and registered as a clinical trial [NCT03501394]. The planned study incorporated four MRI scans across the entire NACT, but was interrupted and closed prematurely due to the COVID-19 pandemic. Therefore, analysis was conducted on MRI scans acquired at baseline and after Cycle 1 only.
Figure 1 Intravoxel incoherent motion (IVIM) images were acquired before neoadjuvant chemotherapy (NACT) at baseline and after Cycle 1. A Bayesian probability (BP) IVIM model was used to compute perfusion fraction (f), tissue diffusion (D), and pseudodiffusion (D*) for the assessment of prediction and early tumour response to NACT. The baseline and percentage change in f, D, and D* at Cycle 1 were examined between good and poor responders, with patients grouped according to the Miller–Payne system for pathological response (RQ1). Medians of baseline and percentage change in f, D, and D* were compared against tumour cellular proliferation marker Ki-67 at core biopsy and percentage change in Ki-67, respectively, from immunostaining in histopathology (RQ2).
Clinical Procedure: Seventeen female patients (age 37–71 years), with grade II or III invasive breast carcinoma from core biopsy and planned for NACT were recruited into the study. Patients with a previous history of breast cancer or receiving hormonal treatment were not eligible. All patients received 5-fluorouracil 500 mg/m2, epirubicin 100 mg/m2, and cyclophosphamide 500 mg/m2 (FEC) once every 21 days for the first three cycles, and docetaxel 100 mg/m2 once every 21 days for the remaining three cycles (28, 29). Two patients with HER2-positive breast cancer additionally received pertuzumab and trastuzumab for a year (30, 31). MRI scans were performed at 5–10 days (median: 7) before the start of the treatment and 10–14 days (median: 12) after Cycle 1. MRI was acquired from 17 patients at baseline and 16 patients at Cycle 1 due to complications in one patient. Standard clinical histopathological examination was performed for each patient to determine histological grade, and immunostaining of Ki-67, a nuclear marker of cellular proliferation associated with worse survival outcomes (32), was conducted in a single batch. The histology results were obtained from core biopsies before NACT and resected residual tumours after six cycles respectively, with appropriate positive controls (33). The pathological response was assessed on the resected tumours, and the good responders and poor responders were identified as above (grades 4 and 5) or below (grades 1, 2, and 3) 90% reduction in cellularity, respectively, according to the Miller–Payne system (12). The percentage change in Ki-67 was computed as the difference between biopsy and excision, normalised to biopsy: [Ki-67 in resected tumour – Ki-67 in core biopsy]/Ki-67 in core biopsy × 100%.
Magnetic Resonance Imaging: All images were acquired on a 3 T clinical whole-body MRI scanner (Achieva TX, Philips Healthcare, Best, The Netherlands), using body coil for uniform transmission and a 16-channel breast coil for signal detection. Patients were in prone position with the imaging volume centered on the breast affected by tumour. IVIM images were acquired in the sagittal orientation using pulsed gradient spin echo (PGSE) sequence with single-shot echo planar imaging (EPI) at 10 diffusion weightings (b-values at 0, 30, 60, 90, 120, 250, 400, 600, 800, and 1,000 s/mm2) (34). For each b-value, diffusion gradients were applied along three orthogonal directions, and the image was computed as the average across the three directions. Images were acquired with a diffusion time (δ/Δ) of 13.1/25.4 ms, a field of view (FOV) of 240 mm × 240 mm, an in-plane resolution of 2.5 mm × 2.5 mm, a slice thickness of 5 mm, an acceleration factor of 2, a repetition time (TR) of 2,400 ms, and an echo time (TE) of 50 ms.
Image Analysis: Bayesian IVIM was performed in MATLAB (R2020a, Mathworks, Natick, MA, USA). The tumour was delineated on dynamic contrast-enhanced MRI by a consultant radiologist in ImageJ (v1.58k, National Institute of Health, Bethesda, MD, USA), with adjustment of image resolution to match IVIM images and conservative definition of tumour boundary to avoid the necrotic, hemorrhagic, and cystic areas. The size of the tumours was evaluated based on the longest diameter from the high-resolution dynamic contrast-enhanced (DCE)-MRI (21, 35, 36). The Bayesian algorithm estimated the joint posterior distribution using the Rician noise likelihood function and uniform joint prior distribution, based on previous literature for Bayesian IVIM model fitting (37). The Bayesian fitting used a Markov chain Monte Carlo setup with Gibbs sampling and Metropolis-Hastings algorithm to derive a marginalised parameter distribution. The step-length parameters were updated every 2,000 iterations, with a total of 20,000 iterations. The conventional IVIM analysis algorithms, including nonlinear least squares full fitting, segmented-unconstrained, and segmented-constrained (38), were also deployed in supplementation to the study (Supplementary Data: Appendix A). The correction for the noise floor (39) was not undertaken since the data have a sufficiently high signal-to-noise ratio (SNR), and the same consistent approach was adopted for all the longitudinal data. The median perfusion fraction (f), tissue diffusion (D), and pseudodiffusion (D*) within the tumour, representing volume fraction between capillary blood and tissue water, mean diffusivity of the tissue, and vascular blood flow motion, respectively, were calculated for baseline and Cycle 1. The percentage change in perfusion fraction, diffusion, and pseudodiffusion at Cycle 1 was computed with reference to baseline: [Cycle1(f/D/D*) – Baseline(f/D/D*)]/Baseline(f/D/D*) × 100% (34).
Statistical Analysis: Statistical analysis was performed using the R software (v3.6.3, The R Foundation for Statistical Computing, Vienna, Austria). The normality of the distribution was assessed using the Shapiro–Wilk test. The measures at baseline and percentage change at Cycle 1 of perfusion fraction, diffusion, and pseudodiffusion were compared between good and poor responder groups using Wilcoxon rank sum test to determine the prediction and early sensitivity of the markers. The correlation of perfusion fraction, diffusion, and pseudodiffusion at baseline against Ki-67 from core biopsy for treatment-naïve prognosis was performed using Spearman’s rank correlation test. The percentage change in perfusion fraction, diffusion, and pseudodiffusion against percentage change in Ki-67 for treatment-altered prognosis was also performed using Spearman’s test. A p-value < 0.05 was considered statistically significant.
3 Results
The patient demographics is shown in Table 1. Among the 17 patients, there were 8 good responders and 9 poor responders at baseline, and due to complications, 1 patient did not complete an MR scan at Cycle 1. There was no significant difference in age and tumour size at baseline between good and poor responders. There was no significant difference in the change in tumour size between good and poor responders at Cycle 1 (Table 1).
There was no significant difference in Bayesian perfusion fraction (p = 0.481), tissue diffusion (p = 0.743), and pseudodiffusion (p = 0.673) at baseline between good and poor responders (Table 2; Supplementary Figure A1). There was also no significant difference in perfusion fraction, tissue diffusion, and pseudodiffusion at baseline between good and poor responders from full fitting and segmented analyses (Supplementary Table A1, Supplementary Figure A2). There was no significant correlation in Bayesian perfusion fraction (p = 0.480), tissue diffusion (p = 0.174), and pseudodiffusion (p = 0.474) at baseline against Ki-67 from core biopsy (Table 2, Supplementary Figure A1). There was also no significant correlation in perfusion fraction, tissue diffusion, and pseudodiffusion at baseline against Ki-67 from core biopsy from full fitting and segmented analyses (Supplementary Table A1, Supplementary Figure A3).
Table 2 Comparison of IVIM-derived parameters between responder groups before and after the first cycle of NACT and the association with Ki-67.
There was a significant difference (p = 0.042) in percentage change in Bayesian perfusion fraction between good and poor responders at Cycle 1, with a decrease in good responders [−7.98% (−19.47–1.73), n = 7] against an increase in poor responders [10.04% (5.09–28.93), n = 9] (Figure 2A, Table 2). There was no significant difference in percentage change in perfusion fraction between good and poor responders at Cycle 1 from full fitting and segmented analyses (Figure 3; Supplementary Table A1). There was a significant correlation in percentage change in Bayesian perfusion fraction (p = 0.042, Figure 4A, Table 2) against percentage change in Ki-67. There was no significant correlation in percentage change in perfusion fraction against percentage change in Ki-67 from full fitting and segmented analyses (Figure 5; Supplementary Table A1). There was no significant difference in percentage change in Bayesian tissue diffusion (p = 0.606, Figure 2B, Table 2) and pseudodiffusion (p = 0.918, Figure 2C, Table 2) between good and poor responders at Cycle 1. There was no significant correlation in percentage change in Bayesian tissue diffusion (p = 0.217, Figure 4B, Table 2) and pseudodiffusion (p = 0.075, Figure 4C, Table 2) against percentage change in Ki-67. There was also no significant difference in percentage change in tissue diffusion and pseudodiffusion between good and poor responders (Figure 3, Supplementary Table A1), nor was there correlation against percentage change in Ki-67 from full fitting and segmented analyses (Figure 5; Supplementary Table A1).
Figure 2 Percentage change in (A) perfusion fraction (f), (B) tissue diffusion (D), and (C) pseudodiffusion (D*) between good and poor responders at first treatment cycle (Cycle 1) from the Bayesian probability (BP) IVIM model. There was a significant difference in percentage change in perfusion fraction between good and poor responders, but not in tissue diffusion and pseudodiffusion. Each dot represents the percentage change in f, D, and D* from an individual patient. Error bar represents median (IQR). Statistically significant p-values (<0.05) are shown on the upper right corner with “*”.
Figure 3 Percentage change in perfusion fraction (f), tissue diffusion (D), and pseudodiffusion (D*) between good and poor responders from nonlinear least squares (Free), segmented-unconstrained (SU), and segmented-constrained (SC) IVIM models. The percentage change in f, D, and D* between good and poor responders from (A–C) Free, (D–F) SU, and (G–I) SC algorithms are shown in dot plots. Each dot represents the IVIM-derived parameter of an individual patient. Error bar represents median (IQR).
Figure 4 Percentage change in perfusion fraction (f), tissue diffusion (D), and pseudodiffusion (D*) against percentage change in the tumour cellular proliferation marker Ki-67. The correlation of percentage change in (A) f, (B) D, and (C) D* at Cycle 1 against percentage change in Ki-67 in resected tumour is shown in scatter plots. Spearman’s rank correlation coefficient (rho (ρ)) was used for correlation analysis and respective ρ score and p-value are shown on each plot. Statistically significant p-values (<0.05) are marked by “*”.
Figure 5 Percentage change in perfusion fraction (f), tissue diffusion (D), and pseudodiffusion (D*) against percentage change in Ki-67. The correlations of percentage change in f, D, and D* from (A–C) nonlinear least squares (Free), (D–F) segmented-unconstrained (SU), and (G–I) segmented-constrained (SC) algorithms at Cycle 1 against percentage change in tumour cellular proliferation marker Ki-67 in resected tumour are shown in scatter plots. Spearman’s rank correlation coefficient (rho (ρ)) was used for correlation analysis and respective ρ score and p-value are shown on each plot.
The parametric maps from IVIM analysis from a typical good and poor responder at baseline and Cycle 1 are shown in Figure 6. The Ki-67-stained microscopy slides from a typical good and poor responder at core biopsy and excision are shown in Figure 7.
Figure 6 Parametric maps from IVIM Bayesian analysis of f, D, and D* from a typical (A) good responder and (B) poor responder at baseline (V1) and Cycle 1 (V2) of NACT (overlaid on diffusion weighted images, b = 1000 s/mm2). Images were acquired with a field of view of 240 mm × 240 mm, an in-plane resolution of 2.5 mm × 2.5 mm, a repetition time of 2,400 ms, and an echo time of 50 ms.
Figure 7 Ki-67 stained microscopy slides from a typical good and poor responder of neoadjuvant chemotherapy (NACT). (A) In the good responder, the Ki-67 score was 17.5% in the core biopsy and 0.8% in the resected tumour. (B) In the poor responder, the Ki-67 score was 23.7% in the core biopsy and 12.4% in the resected tumour. Sections at the greatest dimension of the specimens are shown. Magnification, ×10.
4 Discussion
In this study, we investigated predictive and early response markers for NACT in breast cancer using perfusion fraction, diffusion, and pseudodiffusion derived from BP IVIM. We found that perfusion fraction showed a significant alteration between baseline and Cycle 1 in good responders compared to poor responders, and the alteration is correlated with the change in proliferative activity accumulated through the whole course of NACT across the cohort. However, we did not observe significant differences in alterations in diffusion or pseudodiffusion at Cycle 1 between groups or their correlation against change in proliferative activities. We further did not observe significant differences in imaging markers at baseline between groups, or any significant correlation against proliferative activities at baseline.
The imaging markers of tissue diffusion, perfusion fraction, and pseudodiffusion at baseline did not predict NACT response, indicating the absence of evidence to use tissue diffusion and perfusion at baseline to guide NACT. The results were in agreement with imaging markers of diffusion tensor imaging and ADC at baseline that did not have predictive value for pCR after eight cycles of NACT (22). The results also agreed with a recent study showing that pretreatment tissue diffusion, perfusion fraction, and pseudodiffusion from the segmented constrained model were not predictors of response in patients undergoing a comparable regimen of NACT (36). Diffusion and perfusion metrics estimate cellularity and angiogenesis, respectively, and the lack of a difference between responder groups indicated that a tumour with high cell density and vascular abnormality at initial presentation might not determine the effectiveness of NACT, despite an initial poorer prognosis. Diffusion and perfusion metrics showed no correlation with Ki-67 prior to NACT, indicating no direct correlation between imaging markers of cellularity and angiogenesis with treatment-naïve prognosis, although tissue sampling error could not be excluded.
There was an early significant decrease in perfusion fraction f in good responders, indicating that perfusion fraction might be a sensitive marker in the early identification of a successful NACT. The increase in stiffness of capillary vasculature obstructs microcirculation (40), leading to a modulation of perfusion in the tumour (41). Perfusion fraction has been shown to drastically decrease following a reduction in vascular blood flow motion, despite a subtle structural change in the functional capillary network (18). The susceptibility to systemic changes was lower in comparison to diffusion and pseudodiffusion as independent measures for the physiological response in cellularity and angiogenesis subsequent to cell apoptosis (42). Bayesian-derived perfusion fraction not only showed the potential of perfusion fraction as a marker to predict pathological complete response after one cycle of NACT, in agreement with a previous study (36), but also demonstrated a higher sensitivity since the full fitting and segmented analyses conducted in supplementation to Bayesian showed no group difference. The use of probability constraints on neighboring voxels in the Bayesian model led to less susceptibility of perfusion fraction and pseudodiffusion to the impact of noise, and improved the robustness of fitting (27). However, the higher demand on computing power may delay early adoption, whereas segmented analysis has an added advantage due to the faster processing time and initial validity in a recent study (36). Although perfusion volume ratio from DCE-MRI has been suggested as a marker of responders after one treatment cycle (43), DCE-MRI suffers from nonspecific contrast enhancement from post-treatment changes, including reactive inflammation, necrosis, and peritumoural oedema (18), requiring inputs from more than one radiologist (35). The sensitivity of DCE-MRI to angiogenesis (44) is also limited by the accuracy in the measurement of arterial input function in kinetic hemodynamic models (45) and specialist quantitative deconvolution analysis (16). The inconsistency in terminology between radiology and research practice Dickie et al. 20231 further hinders the wider clinical adoption of quantitative perfusion maps from DCE-MRI for early response in NACT. IVIM, incorporating tissue diffusion and perfusion, shows clinical relevance in the current and previous studies (34, 36, 41, 46), does not require contrast, and has a clearer set of terminology to aid clinical translation. However, the higher susceptibility to noise demands extended acquisition time to reach submillimeter resolution sufficient for accurate determination of tumour size.
The significant correlation between percentage change in perfusion fraction after one cycle and percentage change in Ki-67 indicates a strong association between capillary blood-to-tumour water volume ratio with proliferative activities. Although a causal relationship for the primary impact of NACT on proliferative activity or blood supply could not be established, a reduction in metabolic demand from stunned proliferation and limitation of blood supply from restricted perfusion are both central characteristics of a successful NACT (46). The association between proliferative activity and perfusion has been shown in cell and ex vivo studies as central to tumour development (47, 48). Ki-67 was positively correlated with median (35) and mean (49) tumour perfusion fraction respectively in cross-sectional studies. Thus, an increase in proliferative activity has a corresponding increase in volume fraction between capillary blood and tissue water. Bayesian-derived perfusion fraction showed that good responders with a greater decrease in Ki-67 across NACT also had a greater decrease in perfusion fraction at one cycle, therefore enhancing the critical evidence in the clinical population from a longitudinal study. However, simultaneous full fitting and segmented analyses showed no correlation between change in perfusion fraction and Ki-67. A decrease from high pre-NACT (>35%) to low post-NACT (<15%) Ki-67 showed a sustained low recurrence (<20%) at 3 years after diagnosis (32), and post-NACT Ki-67 proliferative index is an independent prognostic marker in addition to pCR (32). The results showed the potential of perfusion fraction in early response for treatment-altered prognosis and the clinical relevance of an imaging biomarker in the targeted evaluation of the impact of NACT on breast tumours.
There was no significant difference in alteration of tissue diffusion D between responder groups, indicating that cellularity might not be the correct biological target to reveal the effectiveness of NACT at Cycle 1. It has been shown that an increase in tissue diffusion at the second (34, 41) and third (46) cycle was associated with good response in NACT; however, the time points are at a later stage of NACT and metabolic change at an earlier time point might precede morphological change in cellularity and hence tissue diffusion (34, 41, 46). There was a limited number of cytological or histological studies on changes in cellularity and metabolism during the early phase of NACT, potentially due to the heterogeneity across tumour and the fact that biopsy suffers from partial sampling error. There was a decrease in cellularity in biopsy obtained from good responders after two cycles of NACT (50), although the authors in the current study did not find any study on the direct assessment of cellularity after one cycle of conventional NACT. However, a low tumour cellularity in biopsy at day 15 in patients treated with anti-HER-based chemotherapy (including lapatinib and trastuzumab) (51) and a decrease in cellular proliferative activity of Ki-67 after one cycle of conventional NACT (52) predicted good responders. There was no significant difference in alteration of pseudodiffusion D* between responder groups, in agreement with previous breast cancer treatment studies (34, 53). The results might be due to the higher variability in vascular blood flow motion within the capillary bed (54). The lack of association between alterations in tissue diffusion and pseudodiffusion against alterations in Ki-67 showed an absence of evidence for a direct link between early response in cellularity and vascular blood flow motion against change in proliferative activity in the course of NACT.
Bayesian algorithm offers a robust assessment, and an improved estimation of perfusion fraction in association with pathology. The results of the study suggest that perfusion fraction might be a sensitive biomarker of NACT to improve treatment planning, reduce side effects, and expedite precision medicine. Mammography and breast ultrasound have been proposed at the halfway point of NACT to measure the residual tumour size using the RECIST criteria (55); however, tumour regression is not an accurate predictor of response at the first (56) or second (34) cycle of NACT. There was no correlation in size reduction with tumour grade decrease after two cycles of NACT (50), and a reduction in size of the tumours was seen in both small and large tumours (57), potentially due to the formation of islands of nonviable tumour cells subsequent to NACT (50). Perfusion fraction has the potential for tumour perfusion rate characterisation and responder identification after the first cycle, and the correlation with the change in Ki-67 showed that perfusion fraction might have a unique prognostic value in response-guided NACT prior to surgical intervention.
This investigation was a prospective, registered clinical trial that recruited consecutive patients, and set timing for individual MRI scans ensured comparability between patients (46). This study on patient data provided important clinical evidence to a previous study that used simulated and volunteer data (38) and showed that the Bayesian model might ensure greater accuracy of perfusion fraction in association with pathology for differentiation between good and poor responders. A threshold might not be clear cut, and hence, IVIM will contribute to NACT early responder identification but not as a standalone test. Future large cohort studies that will give an accurate estimation of sensitivity and specificity are required to demonstrate the potential of the Bayesian IVIM model to support early response markers in breast cancer management. A three-direction acquisition scheme was utilised due to limited acquisition time (58) and potential risk of overfitting with DTI parameters very sensitive to noise (39); however, a six-direction scheme (or more) might be used to mitigate the impact of anisotropy in the breast (39, 59, 60) in a future study. The current analysis might also benefit from a multi-compartmental IVIM model to account for the exchange between the extracellular and intracellular compartment that affects the quantification of diffusion and pseudodiffusion, since there was a characteristic change in cellular fibrous tissue (stroma) after NACT and the stromal component of the tumour is critical in tumour biology (50, 61).
5 Conclusion
The alteration in perfusion fraction from the Bayesian IVIM model supported the differentiation of good responders from poor responders at the first treatment cycle, and warrants further investigation in comparison to full fitting and segmented analyses in large cohort studies. Early treatment-induced changes in perfusion fraction might serve as non-invasive biomarker to facilitate the delivery of response-guided NACT and the development of an optimal treatment plan.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by London Research Ethics Committee. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
SMC: Data curation, Formal Analysis, Investigation, Methodology, Project administration, Validation, Visualization, Writing – original draft, Writing – review & editing. W-SW: Formal Analysis, Investigation, Methodology, Writing – original draft. NS: Data curation, Investigation, Project administration, Writing – review & editing. RS: Conceptualization, Funding acquisition, Investigation, Supervision, Writing – review & editing. TM: Conceptualization, Funding acquisition, Investigation, Supervision, Writing – review & editing. TG: Conceptualization, Funding acquisition, Investigation, Supervision, Writing – review & editing. EH: Conceptualization, Formal Analysis, Funding acquisition, Investigation, Supervision, Writing – review & editing. YM: Conceptualization, Funding acquisition, Investigation, Supervision, Writing – review & editing. JH: Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The study was jointly funded by the National Health Service Grampian Endowment Research Fund (16/11/047), Friends of Aberdeen and North Centre for Haematology, Oncology and Radiotherapy (RS17 004) and Tenovus Scotland (G16.09). SMC PhD study was jointly supported by Elphinstone scholarship, Roland Sutton Academic Trust and John Mallard scholarship and is currently funded by Cancer Research UK (C68628/A28312). NS PhD study was supported by Biotechnology and Biological Sciences Research Council (1654748, BB/M010996/1). The funding sources were not involved in the study design, in the collection, analysis and interpretation of data, in the writing of the report nor in the decision to submit the article for publication.
Acknowledgments
The authors would like to thank Dr. Matthew Clemence (Philips Healthcare Clinical Science, UK) for clinical scientist support, Ms. Erica Banks and Ms. Alison McKay for patient recruitment support, Ms. Teresa Morris and Ms. Dawn Younie for logistics support, Ms. Beverly McLennan, Ms. Nichola Crouch, Ms. Laura Reid, Mr. Mike Hendry for radiographer support, and Dr. Gordon Urquhart for providing access to the patients.
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.1277556/full#supplementary-material
Footnotes
- ^ Dickie, B.R., Ahmed, Z., Arvidsson, J., Bell, L.C., Buckley, D.L., Debus, C., et al. (2023). A community-endorsed open-source lexicon for contrast agent-based perfusion MRI: A consensus guidelines report from the ISMRM Open Science Initiative for Perfusion Imaging (OSIPI). Magn Reson Med. Online ahead of print. doi: 10.1002/mrm.29840
References
1. Torrisi R, Marrazzo E, Agostinetto E, De Sanctis R, Losurdo A, Masci G, et al. Neoadjuvant chemotherapy in hormone receptor-positive/HER2-negative early breast cancer: When, why and what? Crit Rev Oncol Hematol (2021) 160:103280. doi: 10.1016/j.critrevonc.2021.103280
2. Banys-Paluchowski M, Gasparri ML, de Boniface J, Gentilini O, Stickeler E, Hartmann S, et al. Surgical management of the axilla in clinically node-positive breast cancer patients converting to clinical node negativity through neoadjuvant chemotherapy: current status, knowledge gaps, and rationale for the EUBREAST-03 AXSANA study. Cancers (Basel) (2021) 13(7):1565. doi: 10.3390/cancers13071565
3. NHS England, NHS Improvement. National Cost Collection 2019/20 Report - National Schedule of NHS Costs. London: NHS (2021).
4. Redana S, Sharp A, Lote H, Mohammed K, Papadimitraki E, Capelan M, et al. Rates of major complications during neoadjuvant and adjuvant chemotherapy for early breast cancer: An off study population. Breast (2016) 30:13–8. doi: 10.1016/j.breast.2016.07.019
5. Cain H, Macpherson IR, Beresford M, Pinder SE, Pong J, Dixon JM. Neoadjuvant therapy in early breast cancer: treatment considerations and common debates in practice. Clin Oncol (R Coll Radiol) (2017) 29(10):642–52. doi: 10.1016/j.clon.2017.06.003
6. von Minckwitz G, Untch M, Blohmer JU, Costa SD, Eidtmann H, Fasching PA, et al. Definition and impact of pathologic complete response on prognosis after neoadjuvant chemotherapy in various intrinsic breast cancer subtypes. J Clin Oncol (2012) 30(15):1796–804. doi: 10.1200/JCO.2011.38.8595
7. Potter DA, Herrera-Ponzanelli CA, Hinojosa D, Castillo R, Hernandez-Cruz I, Arrieta VA, et al. Recent advances in neoadjuvant therapy for breast cancer. Fac Rev (2021) 10:2. doi: 10.12703/r/10-2
8. Kim MM, Allen P, Gonzalez-Angulo AM, Woodward WA, Meric-Bernstam F, Buzdar AU, et al. Pathologic complete response to neoadjuvant chemotherapy with trastuzumab predicts for improved survival in women with HER2-overexpressing breast cancer. Ann Oncol (2013) 24(8):1999–2004. doi: 10.1093/annonc/mdt131
9. Golshan M, Loibl S, Wong SM, Houber JB, O'Shaughnessy J, Rugo HS, et al. Breast conservation after neoadjuvant chemotherapy for triple-negative breast cancer: surgical results from the brighTNess randomized clinical trial. JAMA Surg (2020) 155(3):e195410. doi: 10.1001/jamasurg.2019.5410
10. Bosch AM, Kessels AG, Beets GL, Rupa JD, Koster D, van Engelshoven JM, et al. Preoperative estimation of the pathological breast tumour size by physical examination, mammography and ultrasound: a prospective study on 105 invasive tumours. Eur J Radiol (2003) 48(3):285–92. doi: 10.1016/S0720-048X(03)00081-0
11. 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 Cancer (2009) 45(2):228–47. doi: 10.1016/j.ejca.2008.10.026
12. Ogston KN, Miller ID, Payne S, Hutcheon AW, Sarkar TK, Smith I, et al. A new histological grading system to assess response of breast cancers to primary chemotherapy: prognostic significance and survival. Breast (2003) 12(5):320–7. doi: 10.1016/S0960-9776(03)00106-1
13. Baliyan V, Das CJ, Sharma R, Gupta AK. Diffusion weighted imaging: Technique and applications. World J Radiol (2016) 8(9):785–98. doi: 10.4329/wjr.v8.i9.785
14. Partridge SC, Nissan N, Rahbar H, Kitsch AE, Sigmund EE. Diffusion-weighted breast MRI: Clinical applications and emerging techniques. J Magn Reson Imaging (2017) 45(2):337–55. doi: 10.1002/jmri.25479
15. Gao W, Guo N, Dong T. Diffusion-weighted imaging in monitoring the pathological response to neoadjuvant chemotherapy in patients with breast cancer: a meta-analysis. World J Surg Oncol (2018) 16(1):145–y. doi: 10.1186/s12957-018-1438-y
16. Prevos R, Smidt ML, Tjan-Heijnen VC, van Goethem M, Beets-Tan RG, Wildberger JE, et al. Pre-treatment differences and early response monitoring of neoadjuvant chemotherapy in breast cancer patients using magnetic resonance imaging: a systematic review. Eur Radiol (2012) 22(12):2607–16. doi: 10.1007/s00330-012-2653-5
17. von Minckwitz G, Blohmer JU, Costa SD, Denkert C, Eidtmann H, Eiermann W, et al. Response-guided neoadjuvant chemotherapy for breast cancer. JCO (2013) 31(29):3623–30. doi: 10.1200/JCO.2012.45.0940
18. Mendez AM, Fang LK, Meriwether CH, Batasin SJ, Loubrie S, Rodríguez-Soto AE, et al. Diffusion breast MRI: current standard and emerging techniques. Front Oncol (2022) 12:844790. doi: 10.3389/fonc.2022.844790
19. Chen X, Li WL, Zhang YL, Wu Q, Guo YM, Bai ZL. Meta-analysis of quantitative diffusion-weighted MR imaging in the differential diagnosis of breast lesions. BMC Cancer (2010) 10:693. doi: 10.1186/1471-2407-10-693
20. Chen L, Liu M, Bao J, Xia Y, Zhang J, Zhang L, et al. The correlation between apparent diffusion coefficient and tumor cellularity in patients: a meta-analysis. PloS One (2013) 8(11):e79008. doi: 10.1371/journal.pone.0079008
21. Partridge SC, Zhang Z, Newitt DC, Gibbs JE, Chenevert TL, Rosen MA, et al. Diffusion-weighted MRI findings predict pathologic response in neoadjuvant treatment of breast cancer: the ACRIN 6698 multicenter trial. Radiology (2018) 289(3):618–27. doi: 10.1148/radiol.2018180273
22. Furman-Haran E, Nissan N, Ricart-Selma V, Martinez-Rubio C, Degani H, Camps-Herrero J. Quantitative evaluation of breast cancer response to neoadjuvant chemotherapy by diffusion tensor imaging: Initial results. J Magn Reson Imaging (2018) 47(4):1080–90. doi: 10.1002/jmri.25855
23. Zhang D, Geng X, Suo S, Zhuang Z, Gu Y, Hua J. The predictive value of DKI in breast cancer: Does tumour subtype affect pathological response evaluations? Magn Reson Imaging (2022) 85:28–34. doi: 10.1016/j.mri.2021.10.013
24. Senn N, Masannat Y, Husain E, Siow B, Heys SD, He J. q-Space Imaging Yields a Higher Effect Gradient to Assess Cellularity than Conventional Diffusion-weighted Imaging Methods at 3.0 T: A Pilot Study with Freshly Excised Whole-Breast Tumors. Radiol Imaging Cancer (2019) 1(1):e190008. doi: 10.1148/rycan.2019190008
25. Liu C, Liang C, Liu Z, Zhang S, Huang B. Intravoxel incoherent motion (IVIM) in evaluation of breast lesions: comparison with conventional DWI. Eur J Radiol (2013) 82(12):782. doi: 10.1016/j.ejrad.2013.08.006
26. Orton MR, Collins DJ, Koh DM, Leach MO. Improved intravoxel incoherent motion analysis of diffusion weighted imaging by data driven Bayesian modeling. Magn Reson Med (2014) 71(1):411–20. doi: 10.1002/mrm.24649
27. Vidić I, Jerome NP, Bathen TF, Goa PE, While PT. Accuracy of breast cancer lesion classification using intravoxel incoherent motion diffusion-weighted imaging is improved by the inclusion of global or local prior knowledge with bayesian methods. J Magn Reson Imaging (2019) 50(5):1478–88. doi: 10.1002/jmri.26772
28. Roché H, Fumoleau P, Spielmann M, Canon JL, Delozier T, Serin D, et al. Sequential adjuvant epirubicin-based and docetaxel chemotherapy for node-positive breast cancer patients: the FNCLCC PACS 01 Trial. J Clin Oncol (2006) 24(36):5664–71. doi: 10.1200/JCO.2006.07.3916
29. Zaheed M, Wilcken N, Willson ML, O'Connell DL, Goodwin A. Sequencing of anthracyclines and taxanes in neoadjuvant and adjuvant therapy for early breast cancer. Cochrane Database Syst Rev (2019) 2(2):CD012873. doi: 10.1002/14651858.CD012873.pub2
30. Gianni L, Pienkowski T, Im YH, Tseng LM, Liu MC, Lluch A, et al. 5-year analysis of neoadjuvant pertuzumab and trastuzumab in patients with locally advanced, inflammatory, or early-stage HER2-positive breast cancer (NeoSphere): a multicentre, open-label, phase 2 randomised trial. Lancet Oncol (2016) 17(6):791–800. doi: 10.1016/S1470-2045(16)00163-7
31. Liedtke C, Thill M, Jackisch C, Thomssen C, Müller V, Janni W, et al. AGO recommendations for the diagnosis and treatment of patients with early breast cancer: update 2017. Breast Care (Basel) (2017) 12(3):172–83. doi: 10.1159/000477575
32. von Minckwitz G, Schmitt WD, Loibl S, Müller BM, Blohmer JU, Sinn BV, et al. Ki67 measured after neoadjuvant chemotherapy for primary breast cancer. Clin Cancer Res (2013) 19(16):4521–31. doi: 10.1158/1078-0432.CCR-12-3628
33. Cheung SM, Husain E, Mallikourti V, Masannat Y, Heys S, He J. Intra-tumoural lipid composition and lymphovascular invasion in breast cancer via non-invasive magnetic resonance spectroscopy. Eur Radiol (2021) 31(6):3703–11. doi: 10.1007/s00330-020-07502-4
34. Kim Y, Kim SH, Lee HW, Song BJ, Kang BJ, Lee A, et al. Intravoxel incoherent motion diffusion-weighted MRI for predicting response to neoadjuvant chemotherapy in breast cancer. Magn Reson Imaging (2018) 48(5):27–33. doi: 10.1016/j.mri.2017.12.018
35. You C, Li J, Zhi W, Chen Y, Yang W, Gu Y, et al. The volumetric-tumour histogram-based analysis of intravoxel incoherent motion and non-Gaussian diffusion MRI: association with prognostic factors in HER2-positive breast cancer. J Transl Med (2019) 17(1):182–6. doi: 10.1186/s12967-019-1911-6
36. Almutlaq ZM, Wilson DJ, Bacon SE, Sharma N, Stephens S, Dondo T, et al. Evaluation of monoexponential, stretched-exponential and intravoxel incoherent motion MRI diffusion models in early response monitoring to neoadjuvant chemotherapy in patients with breast cancer-A preliminary study. J Magn Reson Imaging (2022) 56(4):1079–88. doi: 10.1002/jmri.28113
37. Jalnefjord O, Andersson M, Montelius M, Starck G, Elf AK, Johanson V, et al. Comparison of methods for estimation of the intravoxel incoherent motion (IVIM) diffusion coefficient (D) and perfusion fraction (f). MAGMA (2018) 31(6):715–23. doi: 10.1007/s10334-018-0697-5
38. Barbieri S, Donati OF, Froehlich JM, Thoeny HC. Impact of the calculation algorithm on biexponential fitting of diffusion-weighted MRI in upper abdominal organs. Magn Reson Med (2016) 75(5):2175–84. doi: 10.1002/mrm.25765
39. Iima M, Partridge SC, Le Bihan D. Six DWI questions you always wanted to know but were afraid to ask: clinical relevance for breast diffusion MRI. Eur Radiol (2020) 30(5):2561–70. doi: 10.1007/s00330-019-06648-0
40. Suresh S. Biomechanics and biophysics of cancer cells. Acta Biomater (2007) 3(4):413–38. doi: 10.1016/j.actbio.2007.04.002
41. Che S, Zhao X, Ou Y, Li J, Wang M, Wu B, et al. Role of the intravoxel incoherent motion diffusion weighted imaging in the pre-treatment prediction and early response monitoring to neoadjuvant chemotherapy in locally advanced breast cancer. Med (Baltimore) (2016) 95(4):e2420. doi: 10.1097/MD.0000000000002420
42. Paran Y, Bendel P, Margalit R, Degani H. Water diffusion in the different microenvironments of breast cancer. NMR Biomed (2004) 17(4):170–80. doi: 10.1002/nbm.882
43. Tudorica A, Oh KY, Chui SY, Roy N, Troxell ML, Naik A, et al. Early prediction and evaluation of breast cancer response to neoadjuvant chemotherapy using quantitative DCE-MRI. Transl Oncol (2016) 9(1):8–17. doi: 10.1016/j.tranon.2015.11.016
44. Li L, Wang K, Sun X, Wang K, Sun Y, Zhang G, et al. Parameters of dynamic contrast-enhanced MRI as imaging markers for angiogenesis and proliferation in human breast cancer. Med Sci Monit (2015) 21:376–82. doi: 10.12659/MSM.892534
45. Petralia G, Summers PE, Agostini A, Ambrosini R, Cianci R, Cristel G, et al. Dynamic contrast-enhanced MRI in oncology: how we do it. Radiol Med (2020) 125(12):1288–300. doi: 10.1007/s11547-020-01220-z
46. Bedair R, Priest AN, Patterson AJ, McLean MA, Graves MJ, Manavaki R, et al. 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. Eur Radiol (2017) 27(7):2726–36. doi: 10.1007/s00330-016-4630-x
47. Baker E, Whiteoak N, Hall L, France J, Wilson D, Bhaskar P. Mammaglobin-A, VEGFR3, and ki67 in human breast cancer pathology and five year survival. Breast Cancer (Auckl) (2019) 13:1178223419858957. doi: 10.1177/1178223419858957
48. Ayoub NM, Jaradat SK, Al-Shami KM, Alkhalifa AE. Targeting angiogenesis in breast cancer: current evidence and future perspectives of novel anti-angiogenic approaches. Front Pharmacol (2022) 13:838133. doi: 10.3389/fphar.2022.838133
49. Feng W, Gao Y, Lu X, Xu Y, Guo Z, Lei J. Correlation between molecular prognostic factors and magnetic resonance imaging intravoxel incoherent motion histogram parameters in breast cancer. Magn Reson Imaging (2022) 85(1):262–70. doi: 10.1016/j.mri.2021.10.027
50. Sethi D, Sen R, Parshad S, Khetarpal S, Garg M, Sen J. Histopathologic changes following neoadjuvant chemotherapy in various Malignancies. Int J Appl Basic Med Res (2012) 2(2):111–6. doi: 10.4103/2229-516X.106353
51. Nuciforo P, Pascual T, Cortes J, Llombart-Cussac A, Fasani R, Pare L, et al. A predictive model of pathologic response based on tumor cellularity and tumor-infiltrating lymphocytes (CelTIL) in HER2-positive breast cancer treated with chemo-free dual HER2 blockade. Ann Oncol (2018) 29(1):170–7. doi: 10.1093/annonc/mdx647
52. Skoog L, Rutqvist LE, Wilking N. Analysis of hormone receptors and proliferation fraction in fine-needle aspirates from primary breast carcinomas during chemotherapy or tamoxifen treatment. Acta Oncol (1992) 31(2):139–41. doi: 10.3109/02841869209088893
53. Suo S, Yin Y, Geng X, Zhang D, Hua J, Cheng F, et al. Diffusion-weighted MRI for predicting pathologic response to neoadjuvant chemotherapy in breast cancer: evaluation with mono-, bi-, and stretched-exponential models. J Transl Med (2021) 19(1):236–3. doi: 10.1186/s12967-021-02886-3
54. Lemke A, Laun FB, Simon D, Stieltjes B, SChad LR. An in vivo verification of the intravoxel incoherent motion effect in diffusion-weighted imaging of the abdomen. Magn Reson Med (2010) 64(6):1580–5. doi: 10.1002/mrm.22565
55. Keune JD, Jeffe DB, Schootman M, Hoffman A, Gillanders WE, Aft RL. Accuracy of ultrasonography and mammography in predicting pathologic response after neoadjuvant chemotherapy for breast cancer. Am J Surg (2010) 199(4):477–84. doi: 10.1016/j.amjsurg.2009.03.012
56. Pavlov MV, Bavrina AP, Plekhanov VI, Golubyatnikov GY, Orlova AG, Subochev PV, et al. Changes in the tumor oxygenation but not in the tumor volume and tumor vascularization reflect early response of breast cancer to neoadjuvant chemotherapy. Breast Cancer Res (2023) 25(1):12. doi: 10.1186/s13058-023-01607-6
57. Sheereen S, Lobo FD, Kumar B, Manoj Kumar S, Santosh Reddy G, Patel W, et al. Histopathological changes in breast cancers following neoadjuvant chemotherapy: implications for assessment of therapy-induced cytological and stromal changes for better clinical outcome and effective patient care. Asian J Oncol (2018) 04(02):61. doi: 10.1055/s-0038-1676909
58. Baltzer P, Mann RM, Iima M, Sigmund EE, Clauser P, Gilbert FJ, et al. Diffusion-weighted imaging of the breast-a consensus and mission statement from the EUSOBI International Breast Diffusion-Weighted Imaging working group. Eur Radiol (2020) 30(3):1436–50. doi: 10.1007/s00330-019-06510-3
59. Partridge SC, Ziadloo A, Murthy R, White SW, Peacock S, Eby PR, et al. Diffusion tensor MRI: preliminary anisotropy measures and mapping of breast tumors. J Magn Reson Imaging (2010) 31(2):339–47. doi: 10.1002/jmri.22045
60. Teruel JR, Goa PE, Sjobakk TE, Ostlie A, Fjosne HE, Bathen TF. Diffusion weighted imaging for the differentiation of breast tumors: From apparent diffusion coefficient to high order diffusion tensor imaging. J Magn Reson Imaging (2016) 43(5):1111–21. doi: 10.1002/jmri.25067
Keywords: diffusion, cellularity, microcirculation, perfusion fraction, pathological response
Citation: Cheung SM, Wu W-S, Senn N, Sharma R, McGoldrick T, Gagliardi T, Husain E, Masannat Y and He J (2023) Towards detection of early response in neoadjuvant chemotherapy of breast cancer using Bayesian intravoxel incoherent motion. Front. Oncol. 13:1277556. doi: 10.3389/fonc.2023.1277556
Received: 14 August 2023; Accepted: 21 November 2023;
Published: 06 December 2023.
Edited by:
Junzhong Xu, Vanderbilt University Medical Center, United StatesReviewed by:
Yen-Peng Liao, University of Texas Southwestern Medical Center, United StatesXiaoyu Jiang, Vanderbilt University, United States
Copyright © 2023 Cheung, Wu, Senn, Sharma, McGoldrick, Gagliardi, Husain, Masannat and He. 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: Sai Man Cheung, Zy5jaGV1bmdAYWJkbi5hYy51aw==