Abstract BACKGROUND Standard of care surgery for high grade glioma (HGG) employs Magnetic Resonance Imaging to identify the contrast enhancing (CE) core for resection, leaving a non-enhancing (NE) invasive rim of tumor which contributes to recurrence. Multi-parametric MRI (mpMRI) has recently been used to map regional biological heterogeneity of neuronal and immune signatures in HGG, and continues to offer a promising strategy for non-invasive surveillance of tumor progression. Here, we present a graph network analysis relating mpMRI to biological signatures in the NE tumor rim and show the utility of communities of imaging features for discriminating tumor phenotypes. METHODS 53 mpMRI features across 101 spatially-localized RNA-sequenced biopsy samples from 64 patients with IDH-wt HGG were represented as nodes within a Neo4j graph network. Modularity optimization was used to identify communities of densely connected nodes. Communities of nodes representing high proportions of NE samples (one-tailed test) were kept for further analysis. Gene set enrichment analysis (GSEA) on samples within selected imaging communities was analyzed for biological process enrichment using clusterProfiler. Single sample GSEA (ssGSEA) using ssGSEA2 was used to annotate neuronal and immune cell states across samples within selected communities. mpMRI feature selection within each community identified the most relevant imaging features for isolating samples classified into neuronal and immune cell states. RESULTS 11 imaging communities gathered a high proportion of NE samples. Three NE communities showed significant enrichment of neuronal signatures, while one showed significant enrichment of immune pathways. One neuronal-predominant community showed unique upregulation of 14 pathways related to neuronal migration/astrocyte projection. CONCLUSION Network imaging communities can identify clusters of NE glioma samples showing contrasting neuronal and immune signatures, including distinct upregulation of migration and projection pathways. Mapping these imaging feature communities to specific biological processes implicated in glioma progression may aid treatment planning and prognostication for glioma patients.
Read full abstract