CANONICAL CORRELATION ANALYSIS IN GLIOBLASTOMA REVEALS ASSOCIATIONS BETWEEN EXPRESSION OF RADIOMIC SIGNATURES AND GENOMICS

Abstract PURPOSE Understanding the molecular underpinnings of imaging signatures of glioblastoma can provide insights into the biologic basis of tumor formation and progression as well as in vivo surrogate markers of molecular events driving the tumor’s phenotype. Through machine learning (ML), this study demonstrates that distinct imaging subtypes of glioblastoma are related to specific molecular alterations. METHODS Pre-operative multi-parametric MRI (T1, T2, T1CE, T2-FLAIR, DSC-MRI, DTI-MRI) of 669 IDH-wildtype subjects with glioblastoma were retrospectively collected and radiomic features, including descriptors of morphology, intensity, histogram, and texture, were extracted. Imaging subtypes were identified by a feature selection and clustering approach. Genomic markers, obtained using next generation sequencing (NGS) panel of 27 key glioblastoma genes, were available in 358/669 patients. Canonical correlation analysis (CCA) was conducted within each imaging subtype between the selected imaging features and genetic variables to seek maximum correlations between combinations of variables in imaging and genomic sets, and hence elucidate the molecular drivers of respective subtypes. RESULTS Three distinct imaging subtypes were identified by clustering on 50 selected features, representing characteristics of morphology, tumor neo-angiogenesis (DSC-derived features), and cellular density (DTI-derived features). These subtypes yielded differentiable overall survival based on Kaplan-Meier analysis. The canonical coefficients of each subtype revealed the distinction of the underlying genomic characteristics: one exhibited frequently mutated [ARID2, NTRK1], another subtype showed increased frequency of mutation in [ATRX, EGFR, PIK3R1], while the third was associated with all these genes and [NF1, PIK3CA, RB1], additionally. CONCLUSION We derived three distinct radiomic MRI subtypes for glioblastoma that highly correlate with the patients' survival and molecular genetic characteristics. Investigating the relationship between imaging and genomic information may enable identification of molecularly- and phenotypically-consistent tumor subtypes, which would offer non-invasive approaches for characterizing heterogeneity of glioblastoma, further facilitating patient stratification and treatment planning.