Functional 4-D clustering for characterizing intratumor heterogeneity in dynamic imaging: evaluation in FDG PET as a prognostic biomarker for breast cancer

Rhea Chitalia1, Varsha Viswanath1, Austin R. Pantel2, Lanell M. Peterson3, Aimilia Gastounioti2, Éric A. Cohen2, Mark Muzi3, Joel S. Karp2, David A. Mankoff2, Despina Kontos2
1Department of Bioengineering, University of Pennsylvania, Philadelphia, USA
2Department of Radiology, University of Pennsylvania, Philadelphia, USA
3Department of Radiology, University of Washington, Seattle, USA

Tóm tắt

Abstract Purpose Probe-based dynamic (4-D) imaging modalities capture breast intratumor heterogeneity both spatially and kinetically. Characterizing heterogeneity through tumor sub-populations with distinct functional behavior may elucidate tumor biology to improve targeted therapy specificity and enable precision clinical decision making. Methods We propose an unsupervised clustering algorithm for 4-D imaging that integrates Markov-Random Field (MRF) image segmentation with time-series analysis to characterize kinetic intratumor heterogeneity. We applied this to dynamic FDG PET scans by identifying distinct time-activity curve (TAC) profiles with spatial proximity constraints. We first evaluated algorithm performance using simulated dynamic data. We then applied our algorithm to a dataset of 50 women with locally advanced breast cancer imaged by dynamic FDG PET prior to treatment and followed to monitor for disease recurrence. A functional tumor heterogeneity (FTH) signature was then extracted from functionally distinct sub-regions within each tumor. Cross-validated time-to-event analysis was performed to assess the prognostic value of FTH signatures compared to established histopathological and kinetic prognostic markers. Results Adding FTH signatures to a baseline model of known predictors of disease recurrence and established FDG PET uptake and kinetic markers improved the concordance statistic (C-statistic) from 0.59 to 0.74 (p = 0.005). Unsupervised hierarchical clustering of the FTH signatures identified two significant (p < 0.001) phenotypes of tumor heterogeneity corresponding to high and low FTH. Distributions of FDG flux, or Ki, were significantly different (p = 0.04) across the two phenotypes. Conclusions Our findings suggest that imaging markers of FTH add independent value beyond standard PET imaging metrics in predicting recurrence-free survival in breast cancer and thus merit further study.

Từ khóa


Tài liệu tham khảo

Polyak K. Breast cancer: origins and evolution. J Clin Invest. 2007;117:3155–63.

Polyak K. Heterogeneity in breast cancer. J Clin Invest. 2011;121:3786–8.

Marusyk A, Polyak K. Tumor heterogeneity: causes and consequences. Biochimica et Biophysica Acta (BBA)-reviews on. Cancer. 2010;1805:105–17.

Sala E, Mema E, Himoto Y, Veeraraghavan H, Brenton J, Snyder A, et al. Unravelling tumour heterogeneity using next-generation imaging: radiomics, radiogenomics, and habitat imaging. Clin Radiol. 2017;72:3–10.

O'Connor JP, Rose CJ, Waterton JC, Carano RA, Parker GJ, Jackson A. Imaging intratumor heterogeneity: role in therapy response, resistance, and clinical outcome. Clin Cancer Res. 2015;21:249–57.

Andor N, Graham TA, Jansen M, Xia LC, Aktipis CA, Petritsch C, et al. Pan-cancer analysis of the extent and consequences of intratumor heterogeneity. Nat Med. 2016;22:105.

Lebron L, Greenspan D, Pandit-Taskar N. PET imaging of breast cancer: role in patient management. PET Clin. 2015;10:159–95.

McDonald ES, Mankoff DA, Mach RH. Novel strategies for breast cancer imaging: new imaging agents to guide treatment. J Nucl Med. 2016;57:69S–74S.

Dunnwald LK, Doot RK, Specht JM, Gralow JR, Ellis GK, Livingston RB, et al. PET tumor metabolism in locally advanced breast cancer patients undergoing neoadjuvant chemotherapy: value of static versus kinetic measures of fluorodeoxyglucose uptake. Clin Cancer Res. 2011;17:2400–9.

Humbert O, Lasserre M, Bertaut A, Fumoleau P, Coutant C, Brunotte F, et al. Breast cancer blood flow and metabolism on dual-acquisition 18F-FDG PET: correlation with tumor phenotype and neoadjuvant chemotherapy response. J Nucl Med. 2018;59:1035–41.

Groheux D, Cochet A, Humbert O, Alberini J-L, Hindié E, Mankoff D. 18F-FDG PET/CT for staging and restaging of breast cancer. J Nucl Med. 2016;57:17S–26S.

Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. 2014;5:4006.

Li H, Zhu Y, Burnside ES, Drukker K, Hoadley KA, Fan C, et al. MR imaging radiomics signatures for predicting the risk of breast cancer recurrence as given by research versions of MammaPrint, Oncotype DX, and PAM50 gene assays. Radiology. 2016;281:382–91.

Chitalia R, Rowland J, McDonald ES, Pantalone L, Cohen EA, Gastounioti A, et al. Imaging phenotypes of breast cancer heterogeneity in pre-operative breast Dynamic Contrast Enhanced Magnetic Resonance Imaging (DCE-MRI) scans predict 10-year recurrence. Clin Cancer Res. 2019:clincanres.4067.2018. https://doi.org/10.1158/1078-0432.ccr-18-4067.

Chitalia RD, Kontos D. Role of texture analysis in breast MRI as a cancer biomarker: a review. J Magn Reson Imaging. 2019;49:927–38.

Eary JF, O'Sullivan F, O'Sullivan J, Conrad EU. Spatial heterogeneity in sarcoma 18F-FDG uptake as a predictor of patient outcome. J Nucl Med. 2008;49:1973–9.

Stoyanova R, Huang K, Sandler K, Cho H, Carlin S, Zanzonico PB, et al. Mapping tumor hypoxia in vivo using pattern recognition of dynamic contrast-enhanced MRI data. Transl Oncol. 2012;5:437–IN2.

Cherezov D, Goldgof D, Hall L, Gillies R, Schabath M, Müller H, et al. Revealing tumor habitats from texture heterogeneity analysis for classification of lung cancer malignancy and aggressiveness. Sci Rep. 2019;9:1–9.

Chang Y-CC, Ackerstaff E, Tschudi Y, Jimenez B, Foltz W, Fisher C, et al. Delineation of tumor habitats based on dynamic contrast enhanced MRI. Sci Rep. 2017;7:9746.

Jardim-Perassi BV, Martinez G, Gillies R. Habitat imaging of tumor evolution by magnetic resonance imaging (MRI). Radiomics and Radiogenomics: Technical Basis and Clinical Applications 2019:115.

Mankoff DA, Dunnwald LK, Gralow JR, Ellis GK, Schubert EK, Tseng J, et al. Changes in blood flow and metabolism in locally advanced breast cancer treated with neoadjuvant chemotherapy. J Nucl Med. 2003;44:1806–14.

Dunnwald LK, Gralow JR, Ellis GK, Livingston RB, Linden HM, Specht JM, et al. Tumor metabolism and blood flow changes by positron emission tomography: relation to survival in patients treated with neoadjuvant chemotherapy for locally advanced breast cancer. J Clin Oncol. 2008;26:4449.

Muzi M, Vesselle H, Grierson JR, Mankoff DA, Schmidt RA, Peterson L, et al. Kinetic analysis of 3′-deoxy-3′-fluorothymidine PET studies: validation studies in patients with lung cancer. J Nucl Med. 2005;46:274–82.

Strydhorst J, Buvat I. Redesign of the GATE PET coincidence sorter. Phys Med Biol. 2016;61:N522.

Karp JS, Viswanath V, Geagan MJ, Muehllehner G, Pantel AR, Parma MJ, Perkins AE, Schmall JP, Werner ME, Daube-Witherspoon ME. PennPET Explorer: design and preliminary performance of a whole-body imager. J Nucl Med. 2020;61(1):136–43.

Viswanath V, Pantel AR, Daube-Witherspoon ME, Doot R, Muzi M, Mankoff DA, Karp JS. Quantifying bias and precision of kinetic parameter estimation on the PennPET Explorer, a long axial field-of-view scanner. IEEE Transactions on Radiation and Plasma Medical Sciences. 2020;4(6):35–749.

Li SZ. Markov random field models in computer vision. European conference on computer vision: Springer, Berlin, Heidelberg. 1994;361–70.

Dempster AP, Laird NM, Rubin DB. Maximum likelihood from incomplete data via the EM algorithm. J R Stat Soc Ser B Methodol. 1977;39:1–22.

Von Luxburg U. A tutorial on spectral clustering. Stat Comput. 2007;17:395–416.

Dice LR. Measures of the amount of ecologic association between species. Ecology. 1945;26:297–302.

Chen JL, Gunn SR, Nixon MS, Gunn RN. Markov random field models for segmentation of PET images. Biennial International Conference on Information Processing in Medical Imaging: Springer. 2001;468–74.

Bhattacharyya A. On a measure of divergence between two statistical populations defined by their probability distributions. Bull Calcutta Math Soc. 1943;35:99–109.

Jahani N, Cohen E, Hsieh M-K, Weinstein SP, Pantalone L, Hylton N, et al. Prediction of treatment response to neoadjuvant chemotherapy for breast cancer via early changes in tumor heterogeneity captured by Dce-MRi registration. Sci Rep. 2019;9:1–12.

Monti S, Tamayo P, Mesirov J, Golub T. Consensus clustering: a resampling-based method for class discovery and visualization of gene expression microarray data. Mach Learn. 2003;52:91–118.

Liu Y, Hayes DN, Nobel A, Marron J. Statistical significance of clustering for high-dimension, low–sample size data. J Am Stat Assoc. 2008;103:1281–93.

Dagogo-Jack I, Shaw AT. Tumour heterogeneity and resistance to cancer therapies. Nature reviews Clinical oncology. 2018;15(2):81.

Soret M, Bacharach SL, Buvat I. Partial-volume effect in PET tumor imaging. J Nucl Med. 2007;48:932–45.

Chan SR, Fowler AM, Allen JA, Zhou D, Dence CS, Sharp TL, et al. Longitudinal noninvasive imaging of progesterone receptor as a predictive biomarker of tumor responsiveness to estrogen deprivation therapy. Clin Cancer Res. 2015;21:1063–70.

O'Sullivan F. Locally constrained mixture representation of dynamic imaging data from PET and MR studies. Biostatistics. 2006;7:318–38.

O’Sullivan F, Muzi M, Mankoff DA, Eary JF, Spence AM, Krohn KA. Voxel-level mapping of tracer kinetics in PET studies: a statistical approach emphasizing tissue life tables. Ann Appl Stat. 2014;8:1065.