The purpose of this study was to identify perfusional subregions sharing similar kinetic characteristics from dynamic contrast-enhanced magnetic resonance imaging (MRI) using data-driven clustering, and to evaluate the effect of perfusional heterogeneity based on those subregions on patients' survival outcomes in various risk models. From two hospitals, 308 and 147 women with invasive breast cancer who underwent preoperative MRI between October 2011 and July 2012 were retrospectively enrolled as development and validation cohorts, respectively. Using the Cox-least absolute shrinkage and selection operator model, a habitat risk score (HRS) was constructed from the radiomics features from the derived habitat map. An HRS-only, clinical, combined habitat, and two conventional radiomics risk models to predict patients' disease-free survival (DFS) were built. Patients were classified into low-risk or high-risk groups using the median cutoff values of each risk score. Five habitats with distinct perfusion patterns were identified. An HRS was an independent risk factor for predicting worse DFS outcomes in the HRS-only risk model (hazard ratio = 3.274 [95% CI = 1.378-7.782]; p = 0.014) and combined habitat risk model (hazard ratio = 4.128 [95% CI = 1.744-9.769]; p = 0.003) in the validation cohort. In the validation cohort, the combined habitat risk model (hazard ratio = 4.128, p = 0.003, C-index = 0.760) showed the best performance among five different risk models. The quantification of perfusion heterogeneity is a potential approach for predicting prognosis and may facilitate personalized, tailored treatment strategies for breast cancer.
Keywords: breast cancer; heterogeneity; kinetics; magnetic resonance imaging; perfusion; prognosis; radiomics.