Erik Toorens, Sanjay Saxena, Donald M. O'Rourke, Hamed Akbari, Arati Desai, Christos Davatzikos, MacLean Nasrallah, Danni Tu, Spyridon Bakas, Elizabeth Mamourian, Anahita Fathi Kazerooni, Chiharu Sako, Costas Koumenis, Stephen J Bagley, Russell T. Shinohara, Vishnu Bashyam, Tapan Ganguly, and Robert A. Lustig
PURPOSE Multi-omics data integration captures tumor characteristics at multiple scales [i.e., microscopic (genomics and epigenetics), macroscopic (radiomics), clinical manifestation], provides a more comprehensive assessment of patient’s risk, and facilitates personalized therapies. In this work, we investigated the synergistic value of such multiple data sources for risk stratification and prediction of overall survival in IDH-wildtype glioblastoma tumors. METHODS Quantitative conventional and deep radiomics were extracted from pre-operative multi-parametric structural MRI (T1, T1Gd, T2, T2-FLAIR) of 501 patients with newly diagnosed glioblastoma. 389/501 and 112/501 patients formed our discovery and replication cohorts, respectively. Conventional radiomics were extracted from CaPTk, and deep radiomics from a pre-trained VGG-19 model. Multivariate SVM classification was performed on the discovery cohort to stratify patients in high, medium, and low-risk groups, using recursive feature elimination and 5-fold cross-validation. This model was independently tested on the replication cohort, and a radiomic-based survival prediction index (SPIradiomics) was calculated for each patient. Multi-stage integration of omics data, i.e., clinical (age, gender, extent of resection (EOR)), SPIradiomics, epigenetics (MGMT promoter methylation), and genomics (27 clinically relevant gene mutations via next-generation sequencing (NGS)), was performed using multivariate Cox proportional hazards (Cox-PH) model for stratification of the risk in the replication cohort. RESULTS Cox-PH modeling resulted in a concordance index (c-index) of 0.65 (95% CI:0.6–0.7) for clinical data, 0.67 (95% CI:0.62–0.72) for clinical and epigenetics, 0.70 (95% CI:0.65–0.75) for clinical and radiomics, 0.72 (95% CI:0.68–0.77) for clinical, epigenetics, and radiomics, and 0.75 (95% CI:0.71 – 0.78) for the multi-omics combination of all data; highlighting the added value of each layer of information in prediction of the patient’s risk. CONCLUSION Our results reinforce the synergistic value of integrated diagnostic methods for improving risk assessment of patients with glioblastoma that may pave the path towards a more personalized treatment planning.