Citations

Abstract Characterizing intra- and inter-tumoral heterogeneity of glioblastoma has historically relied on discrete classifications of malignant cell populations leaving immune and other cell populations, known to exist admixed with the malignant tumor cells, relatively neglected. Manifold learning algorithms can manage deconvolving multiple cell populations and are often used to track cell state transitions in single cell transcriptomics. We applied a manifold learning approach to TCGA microarray data (Nf525) and bulk transcriptomics of 134 image localized biopsies across 30 patients with primary and 9 with recurrent glioblastoma to further elucidate how to organize biopsies across a spectrum of possible tissue states. The algorithm revealed a low-dimensional manifold graph for which each biopsy lives across 3 polarizing tissue states - one that is associated with diffusely invaded brain, one that is enriched in mesenchymal genes, and one that is enriched in classical proliferative tumor signatures. We deconvolved the bulk transcriptomics of the image-localized biopsies to reveal the relative abundance of 18 malignant, immune, and other cell subpopulations in each biopsy. Overlaying the cellular decomposition onto the manifold graph visualizing the tissue state distributions revealed that transitions between states correlate with changes in cellular cohabitation composition. The tumor cellular cohabitation ecologies have the lowest diversity, as inferred by ecological measures such as Shannon entropy and evenness, at the distal poles of the graph when compared to the transitional arms. Further, we found that the relationship between imaging appearance of contrast enhancement on T1-weighted MRI and the biopsy cellular composition varies with sex and primary vs recurrent biopsy status. The limited spectrum of possible tissue states revealed by the manifold learning is suggestive of a limited continuum along which tumor and non-tumoral cell populations can cohabitate. Such a limited low-dimensional biological space may constrain the dynamics of tumor biology in a predictable manner.

Abstract Brain tumor associated edema is a significant cause of patient morbidity and mortality often requiring continual treatment with differing levels of success across patients. Edema is a broad term that can indicate increased local water as well as diffusely infiltrating tumor cells. Edema is visualized as hyperintense regions on T2-weighted (T2W) magnetic resonance imaging (MRI) and is typically associated with blood-brain barrier (BBB) breakdown and tumor aggressiveness. Understanding the molecular mechanisms driving imaging patterns of “edema” could provide insights into clinical imaging interpretation. We have an ongoing image-guided biopsy study that allows us to link biopsy molecular markers with locoregional MRI patterns of edema. Further, we previously developed a physics-based method to estimate edema abundance (i.e., edema map) from T2W MRIs. Our goal was to identify connections between BBB-associated molecular factors and edema abundance in brain tumors. Our cohort included 38 patients (female: 15, male: 23) with 129 image-guided biopsies (female: 62, male: 67). We correlated image-localized edema map values with the mean transcriptional frequency for 57 genes related to BBB function. Additionally, we examined correlations separately according to patient reported sex (i.e., male and female) and imaging phenotype (i.e., ENH: enhancing and NE: non-enhancing). We utilized multiple comparisons corrections with a 5% false discovery rate to determine significance. For the overall cohort, we observed significant positive correlations for the HIF1A (p< 0.001) and SOX2 markers with edema. For NE samples, significant correlations included APOE (p=0.001), HIF1A (p< 0.001), PIK3CA (p< 0.001), PTCH1 (p< 0.001), and SOX2 (p< 0.001). Amongst female samples, a significant correlation with PTCH1 (p=0.002) was observed. There were no significant correlations noted for male and enhancing sub-cohorts. Significant correlations between molecular markers of BBB and edema map values could lead to clinical biomarkers for edema or tumor aggressiveness.

Abstract Background Magnetic Resonance Imaging (MRI) is the mainstay for neurosurgical oncology but not for informing us about glioma biology. An obstacle to developing MR-based glioma biomarkers is the absence of rigorous correlation between MRI features and glioma biology as assessed in multi-regional biopsies, within and across patients. Methods We directly addressed this obstacle by collating a unique cohort of 202 MRI-localized biopsies from 58 patients. We define a low-dimensional transcriptional pseudotime continuum along which heterogeneous high-grade glioma (HGG) samples organize both within and across patients. Results We observe three polarized transcriptional tissue states: infiltrated brain, immune/inflammatory, and proliferative associated with patterns of cohabitation of cellular subpopulations. The states and deconvolved populations show correlation with enhancement status on T1Gd MRI. Moreover, discrete MRI habitats, regions sharing common imaging features, defined as combinations of high or low signal intensity across multiparametric MRI revealed 14 MRI habitats. We order the MRI habitats according to the average pseudotime on the transcriptional continuum. We find that MRI habitats with low pseudotime (associated with early tumor development and diffusely invaded brain tissue) localized at the periphery of the tumor whilst high pseudotime either proliferative or immune/inflammatory states were towards the core of the lesion. We find that composition of MRI habitats are impacted by MGMT status and patient sex. Conclusion This suggests that ongoing aggregation of MRI-localized biopsies may augment our projection of biology onto MRI habitats to support tracking of the evolution of cellular ecologies within and across each patient’s tumor over time. Key Points Trajectory inference of HGG samples reveals a continuum connecting three polarized tissue states reflecting transitions in cellular population ecologies and gene set enrichment Combining multiparametric MRI features into regional MRI habitats reveal clinically-distinct patterns of cell co-habitation and tissue state MGMT methylation status shifts the cellular patterns of cohabitation within MRI habitats Importance of the Study Regional heterogeneity is a defining characteristic of high-grade glioma (HGG). Mapping the transcriptomic biology of multiregionally sampled MRI-localized biopsies from high-grade gliomas onto a low-dimensional space reveals a continuum between three primary tissue states: infiltrated brain, immune/inflammatory, and proliferative. Each state is associated with distinct cell cohabitation ecologies. Transitions between the tissue states coordinate with changes in tumor composition, revealing potential direct connections between tumor biology shifts and their manifestation on clinical imaging.

Abstract INTRODUCTION Dendritic cells (DCs) are potent antigen presenting cells that can be exploited to initiate an adaptive anti-tumoral immune response. DC vaccine clinical trials for primary glioblastoma (GBM) have reported prolonged progression-free survival without any impact on overall survival (OS). We report a radiomics approach that identifies a subpopulation of patients with prolonged OS in clinical trial MC1272. METHODS Twenty adults with primary GBM undergoing standard-of-care therapy were enrolled in MC1272. Autologous DCs were pulsed with allogenic GBM cell lines to generate vaccines that were administered for up to twelve cycles. Standard brain imaging was obtained at the initiation of treatment and two months afterwards. An independent cohort of image-localized biopsies underwent RNA sequencing followed by cellular deconvolution to estimate T cell abundance. A machine learning model was trained to predict intratumoral T cell abundance from imaging features, and the model was applied to MC1272 patient imaging. RESULTS Voxelwise predictions of T cell abundance were generated for each patient’s pre- and post-treatment images. The difference in total intratumoral T cell abundance between imaging timepoints classified patients into increasing or decreasing T cell groups. Patients whose T cells increased had longer OS (median, 21 months) than those whose T cells decreased (median, 10 months; p=0.0035). The significance remained in a Cox proportional hazards analysis that adjusted for patient age and sex (p=0.011). CONCLUSIONS A spatially-resolved radiomics model detected that an intratumoral T cell influx after DC vaccine therapy was associated with prolonged OS. The “post-treatment” imaging in this study was obtained a mere two months after DC vaccine initiation, suggesting that our radiomics model can provide an early indication of treatment responsiveness and prognosis in these patients.

Abstract Background : Sex is recognized as a significant determinant of outcome among glioblastoma patients, but the relative prognostic importance of glioblastoma features has not been thoroughly explored for sex differences. Methods : Combining multi-modal MR images, biomathematical models, and patient clinical information, this investigation assesses which pretreatment variables have a sex-specific impact on the survival of glioblastoma patients (299 males and 195 females). Results : Among males, tumor (T1Gd) radius was a predictor of overall survival (HR=1.027, p=0.044). Among females, higher tumor cell net invasion rate was a significant detriment to overall survival (HR=1.011, p<0.001). Female extreme survivors had significantly smaller tumors (T1Gd) (p=0.010 t-test), but tumor size was not correlated with female overall survival (p=0.955 CPH). Both male and female extreme survivors had significantly lower tumor cell net proliferation rates than other patients (M p=0.004, F p=0.001, t-test). Conclusion : Despite similar distributions of the MR imaging parameters between males and females, there was a sex-specific difference in how these parameters related to outcomes.

Abstract Background : Sex is recognized as a significant determinant of outcome among glioblastoma patients, but the relative prognostic importance of glioblastoma features has not been thoroughly explored for sex differences. Methods : Combining multi-modal MR images, biomathematical models, and patient clinical information, this investigation assesses which pretreatment variables have a sex-specific impact on the survival of glioblastoma patients (299 males and 195 females). Results : Among males, tumor (T1Gd) radius was a predictor of overall survival (HR=1.027, p=0.044). Among females, higher tumor cell net invasion rate was a significant detriment to overall survival (HR=1.011, p<0.001). Female extreme survivors had significantly smaller tumors (T1Gd) (p=0.010 t-test), but tumor size was not correlated with female overall survival (p=0.955 CPH). Both male and female extreme survivors had significantly lower tumor cell net proliferation rates than other patients (M p=0.004, F p=0.001, t-test). Conclusion : Despite similar distributions of the MR imaging parameters between males and females, there was a sex-specific difference in how these parameters related to outcomes.

Abstract Background : Sex is recognized as a significant determinant of outcome among glioblastoma patients, but the relative prognostic importance of glioblastoma features has not been thoroughly explored for sex differences. Methods : Combining multi-modal MR images, biomathematical models, and patient clinical information, this investigation assesses which pretreatment variables have a sex-specific impact on the survival of glioblastoma patients (299 males and 195 females). Results : Among males, tumor (T1Gd) radius was a predictor of overall survival (HR=1.027, p=0.044). Among females, higher tumor cell net invasion rate was a significant detriment to overall survival (HR=1.011, p<0.001). Female extreme survivors had significantly smaller tumors (T1Gd) (p=0.010 t-test), but tumor size was not correlated with female overall survival (p=0.955 CPH). Both male and female extreme survivors had significantly lower tumor cell net proliferation rates than other patients (M p=0.004, F p=0.001, t-test). Conclusion : Despite similar distributions of the MR imaging parameters between males and females, there was a sex-specific difference in how these parameters related to outcomes.

Abstract Glioblastoma (GBM) is the most aggressive primary brain tumor with a median overall survival of 15 months with standard-of-care treatment. GBM patients sometimes present with a cystic component, which can be identified through magnetic resonance imaging (MRI). Previous studies suggest that cysts occur in 7–22% of GBM patients and have reported mixed results regarding whether cystic GBM have a survival benefit compared to noncystic GBM. Using our large retrospective cohort of 493 first-diagnosis GBM patients, we aim to elucidate this link between cystic GBM and survival. Within this cohort, 88 patients had a significant cystic component at presentation as identified on MRI. Compared to noncystic GBM (n=405), cystic GBM patients had significantly better overall survival (15 vs 22 months median, log-rank, p=0.001) and were significantly younger at the time of presentation (t-test, p=0.002). However, within patients that received current standard-of-care treatment (n=184), cystic GBM (n=40) was not as beneficial for outcome (22 vs 25 months, log-rank, p=0.3). We also did not observe a significant survival benefit when comparing this standard-of-care cystic cohort to cystic GBM patients diagnosed before the standard was established (n=19, 25 vs 23 months, log-rank, p=0.3), but the analogous result for noncystic GBM patients gives a sizeable benefit, as expected (n=144, n=111, respectively, 22 vs 12 months, log-rank p < 0.0001). Together, these results on current standard-of-care may explain later studies that note no significant survival benefit for cystic GBM patients receiving current standard-of-care. We also report differences in the absolute and relative sizes of imaging abnormalities on MRI and in prognostic impact of cysts based on sex. We discuss current hypotheses for these observed differences, including the possibility that the presence of a cyst could be indicative of a less aggressive tumor.

Abstract Glioblastomas (GBMs) are known to be complex tumors comprising multiple subpopulations of genetically-distinct cancer cells; it is thought that this genetic variation is a major factor in the lack of observed survival benefit of treatment regimes that target one of these subpopulations. The field of radiogenomics seeks to study correlations between MRI patterns and genetic features of GBM tumors. Spatial radiogenomic maps produced using machine-learning (ML) methods that are trained against information from image-localized patient biopsies identify regions where particular cancer sub-populations are predicted to occur within a GBM, thus non-invasively characterizing the regional genetic variability of these tumors. These tumor subpopulations may also interact with one another, in ways which may be of a competitive or cooperative nature to varying degrees. It is important to ascertain the nature of these interactions, as they may have implications for treatment response to targeted therapies, and characterization of the spatio-temporal dynamics of these co-evolving sub-populations will shed light on why some therapies fail. Here we combine mathematical modeling techniques and spatially-resolved radiogenomic maps to study the nature of these interactions between molecularly-distinct GBM subpopulations. We model the interactions between cell populations using a partial differential equation based formalism. The model is parameterized using radiogenomic ML maps from which we infer the nature of interactions between subpopulations. Furthermore, using maps as inputs, the model turns static maps into dynamic information, thus providing insight into how these subpopulations composing the tumor change over time and the effect this has on observed treatment response for individual patients.

Abstract Glioblastoma (GBM) is the most aggressive primary brain tumor with a median survival of only 15 months with standard of care treatment. Lacunarity, a quantitative morphological measure of how shapes fill space, and fractal dimension, another morphological measure of the complexity of pixel arrangement, of segmented necrotic regions on gadolinium-enhanced T1 weighted (T1gd) MRI have previously been shown to distinguish both overall survival (OS) and progression free survival (PFS) in GBM (n = 95). In our larger patient cohort (n = 389), we sought to validate or refute previously published results connecting morphological metrics and patient survival. We identified pretreatment necrotic regions of our retrospective first-diagnosis GBM patient cohort using segmented T1gd MRI enhancing regions. We calculated lacunarity and fractal dimension across all T1gd MRI slices with enhancing tumor, and used the median lacunarity and fractal dimension values for our analysis. We find that a lacunarity threshold can significantly distinguish OS (14 months vs 19 months median, log-rank p = 0.015, n = 389) and a fractal dimension threshold can significantly distinguish PFS (8 months vs 11 months median, log-rank p = 0.015, n = 123). We believe that morphological metrics such as lacunarity and fractal dimension could play a role in standard-of-care prognostic considerations at tumor presentation. This link between morphological and survival metrics could be driven by underlying biological phenomena or microenvironmental factors that should be further explored.