We analyzed DCE-MR images from 132 women with locally advanced breast cancer from the I-SPY1 trial to evaluate changes of intra-tumor heterogeneity for augmenting early prediction of pathologic complete response (pCR) and recurrence-free survival (RFS) after neoadjuvant chemotherapy (NAC). Utilizing image registration, voxel-wise changes including tumor deformations and changes in DCE-MRI kinetic features were computed to characterize heterogeneous changes within the tumor. Using five-fold cross-validation, logistic regression and Cox regression were performed to model pCR and RFS, respectively. The extracted imaging features were evaluated in augmenting established predictors, including functional tumor volume (FTV) and histopathologic and demographic factors, using the area under the curve (AUC) and the C-statistic as performance measures. The extracted voxel-wise features were also compared to analogous conventional aggregated features to evaluate the potential advantage of voxel-wise analysis. Voxel-wise features improved prediction of pCR (AUC = 0.78 (±0.03) vs 0.71 (±0.04), p < 0.05 and RFS (C-statistic = 0.76 ( ± 0.05), vs 0.63 ( ± 0.01)), p < 0.05, while models based on analogous aggregate imaging features did not show appreciable performance changes (p > 0.05). Furthermore, all selected voxel-wise features demonstrated significant association with outcome (p < 0.05). Thus, precise measures of voxel-wise changes in tumor heterogeneity extracted from registered DCE-MRI scans can improve early prediction of neoadjuvant treatment outcomes in locally advanced breast cancer.

Purpose

Image registration plays a vital role in spatially aligning multiple MRI scans for better longitudinal assessment of tumor morphological features. The objective was to evaluate the effect of registration accuracy of six established deformable registration methods(ANTs, DRAMMS, ART, NiftyReg, SSD-FFD, and NMI-FFD) on the predictive value of extracted radiomic features when modeling recurrence-free-survival(RFS) for women after neoadjuvant chemotherapy(NAC) for locally advanced breast cancer.

Methods

130 women had DCE-MRI scans available from the first two visits in the ISPY1/ACRIN-6657 cohort. We calculated the transformation field from each of the different deformable registration methods, and used it to compute voxel-wise parametric-response-maps(PRM) for established four kinetic features.104-radiomic features were computed from each PRM map to characterize intra-tumor heterogeneity. We evaluated performance for RFS using Cox-regression, C-statistic, and Kaplan-Meier(KM) plots.

Results

A baseline model(F1:Age, Race, and Hormone-receptor-status) had a 0.54 C-statistic, and model F2(baseline + functional-tumor-volume at early treatment visit(FTV₂)) had 0.63. The F2+ANTs had the highest C-statistic(0.72) with the smallest landmark differences(5.40±4.40mm) as compared to other models. The KM curve for model F2 gave p=0.004 for separation between women above and below the median hazard compared to the model F1(p=0.31). A models augmented with radiomic features, also achieved significant KM curve separation(p<0.001) except the F2+ART model.

Conclusion

Incorporating image registration in quantifying changes in tumor heterogeneity during NAC can improve prediction of RFS. Radiomic features of PRM maps derived from warping the DCE-MRI kinetic maps using ANTs registration method further improved the early prediction of RFS as compared to other methods.