Grinband, Jack

PURPOSE Age ≥40 years is regarded as a high-risk feature and an indication for adjuvant chemoradiotherapy for patients with lower-grade glioma in clinical practice guidelines. It is unclear whether age remains a relevant prognostic factor for contemporary definitions of lower-grade gliomas in the molecular era. METHODS The Prospective Gliomas Research (PROGRES) database contains individual patient-level data from 11 prospective clinical trials or observational registries of histologically defined lower-grade 2-3 oligodendroglioma or astrocytoma. We determined the association of age (18-39 years v ≥40 years) with progression-free survival (PFS) stratified by isocitrate dehydrogenase 1 or 2 ( IDH1/2 ) status, using log-rank tests and Cox regression models. We validated our findings in a separate multi-institutional retrospective cohort (Retrospective Glioma Research [REGRES] database). RESULTS We identified 1,619 and 1,292 eligible patients in the PROGRES and REGRES cohorts, respectively. IDH -wildtype tumors were more common in patients 40 years and older (38% v 5%, odds ratio: 11.3 [95% CI, 6.5 to 19.7]). Age was associated with PFS in IDH -wildtype (5-year PFS for ≥40 v 18-39 years: 6% v 24%, hazard ratio [HR], 1.74 [95% CI, 1.21 to 2.50]) but not in IDH -mutant glioma (60% v 59%, HR, 0.89 [95% CI, 0.76 to 1.05], P interaction < .001). In IDH -wildtype tumors, older age predicted aggressive molecular features, including TERT promoter mutation (65% v 28%), EGFR amplification (41% v 15%), and chromosome +7/–10 alteration (57% v 25%). In a pooled analysis of four clinical trials, age was not predictive of a benefit from chemoradiotherapy versus radiotherapy alone for IDH -mutant glioma. CONCLUSION In the absence of additional clinical or molecular risk factors, age alone should not be considered an indication for administration or deferral of adjuvant treatment. Practice guidelines should be revised to reflect contemporary prognostic factors in the molecular era.

Abstract PURPOSE Age over 40 years is commonly regarded as a “high risk” feature and potential indication for adjuvant chemoradiotherapy for patients with lower grade glioma in most clinical practice guidelines. However, these recommendations are based on retrospective observations and their subsequent incorporation into legacy clinical trials. It is unclear whether age remains a relevant prognostic factor in the molecular era. METHODS The Prospective Gliomas Research [PROGRES] Database contains individual patient-level data from 11 prospective clinical trials or observational registries, including seven international trials from the RTOG, EORTC, or Alliance cooperative groups. We determined the association of age with progression-free survival (PFS) in low grade or anaplastic glioma, stratified by isocitrate dehydrogenase (IDH1/2) status. We validated our findings in a separate multi-institutional retrospective cohort (Retrospective Glioma Research [REGRES] Database). RESULTS We identified 1,619 and 1,271 eligible patients in the PROGRES and REGRES cohorts, respectively. IDH-wildtype tumors were more common in patients &gt;40 years (38% vs, 5%, odds ratio: 11.3 [95%CI 6.5 – 19.7]). Age was associated with PFS in IDH-wildtype (5-year PFS: 6% vs, 24%, P =.003) but not in IDH-mutant glioma (60% vs, 59%, P = 0.22, P value of interaction &lt;.001). In IDH-wildtype tumors, older age predicted aggressive molecular features, including TERT promoter mutation (65% vs, 28%), EGFR amplification (41% vs, 15%), and chromosome +7/-10 alteration (57% vs, 25%). Age was not predictive of a benefit from chemoradiotherapy versus radiotherapy alone in IDH-mutant glioma. CONCLUSION Older age is a putative surrogate for IDH-wildtype status, but was not an independent predictor of progression in IDH-mutant lower grade glioma and did not predict benefit from chemoradiotherapy. In the absence of additional clinical or molecular risk factors, age alone should not be considered an indication for administration or deferral adjuvant treatment. Practice guidelines should be revised to reflect contemporary prognostic factors in the molecular era.

Abstract Intratumoral heterogeneity poses a significant challenge to the diagnosis and treatment of glioblastoma (GBM). This heterogeneity is further exacerbated during GBM recurrence, as treatment-induced reactive changes produce additional intratumoral heterogeneity that is ambiguous to differentiate on clinical imaging. There is an urgent need to develop non-invasive approaches to map the heterogeneous landscape of histopathological alterations throughout the entire lesion for each patient. We propose to predictively fuse Magnetic Resonance Imaging (MRI) with the underlying intratumoral heterogeneity in recurrent GBM using machine learning (ML) by leveraging image-localized biopsies with their associated locoregional MRI features. To this end, we develop BioNet, a biologically-informed neural network model, to predict regional distributions of three tissue-specific gene modules: proliferating tumor, reactive/inflammatory cells, and infiltrated brain tissue. BioNet offers valuable insights into the integration of multiple implicit and qualitative biological domain knowledge, which are challenging to describe in mathematical formulations. BioNet performs significantly better than a range of existing methods on cross-validation and blind test datasets. Voxel-level prediction maps of the gene modules by BioNet help reveal intratumoral heterogeneity, which can improve surgical targeting of confirmatory biopsies and evaluation of neuro-oncological treatment effectiveness. The non-invasive nature of the approach can potentially facilitate regular monitoring of the gene modules over time, and making timely therapeutic adjustment. These results also highlight the emerging role of ML in precision medicine.

Abstract Glioblastoma (GBM) is the most aggressive primary brain tumor with a median survival of 14 months. GBMs are challenging to treat due to their heterogeneous nature. It has also been seen that these tumors have sex differences in their cellular subtypes as well as imaging. Radiomics has the potential to provide a non-invasive, spatial understanding of genetic and epigenetic diversity in these complex tumors and to aid in treatment planning. We have an ongoing study to obtain image-localized biopsies from GBM patients, allowing us to complete radiomic analysis and make connections between immunohistochemistry (IHC) and magnetic resonance imaging (MRI) features. We sought to determine if the patterns on imaging were correlated with underlying tumor biology. We focused on immunohistochemistry (IHC) markers of key features of tumor biology including SOX2 for stem-like tumor cells, CD68 for immune response and Ki67 for proliferation kinetics. Our study included 38 patients with a total of 99 biopsies (bxs): 27 males with 77 bxs and 11 females with 22 bxs. Biopsies were sectioned and stained for the SOX2, CD68, and KI67 markers. We computed 18 first-order radiomic features at each biopsy location for patients’ multimodal MRIs: T1W, T1Gd, T2W, FLAIR, apparent diffusion coefficient, diffusion weighted imaging (DWI) and susceptibility weighted imaging. We then performed correlation analysis between each radiomic feature and marker abundance for each IHC stain. Overall, we found sex-distinct patterns connecting imaging with these IHC markers. For example, amongst female patients, DWI held more prominent correlations with SOX2 than in males. Whereas there were more correlations between CD68 IHC abundance and T1Gd imaging features in males compared to females. Taken together, the overall patterns connecting locoregional imaging features to these IHC markers showed sex-distinct patterns suggesting the potential for sex to be an important biological variable when interpreting the biology underlying imaging changes.

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 Glioblastoma (GBM) is a devastating primary brain tumor known for its heterogeneity with a median survival of 15 months. Clinical imaging remains the primary modality to assess brain tumor response, but it is nearly impossible to distinguish between tumor growth and treatment response. Ki67 is a marker of active cell proliferation that shows inter- and intra-patient heterogeneity and should change under many therapies. In this work, we assessed the utility of a semi-supervised deep learning approach for regionally predicting high-vs-low Ki67 in GBM patients based on MRI. We used both labeled and unlabeled datasets to train the model. Labeled data included 114 MRI-localized biopsies from 43 unique GBM patients with available immunohistochemistry Ki67 labels. Unlabeled data included nine repeat routine pretreatment paired scans of newly-diagnosed GBM patients acquired within three days. Data augmentation techniques were utilized to enhance the size of our data and increase generalizability. Data was split between training, validation, and testing sets using 65-15-20 percent ratios. Model inputs were 16x16x3 patches around biopsies on T1Gd and T2 MRIs for labeled data, and around randomly selected patches inside the T2 abnormal region for unlabeled data. The network was a 4-conv layered VGG-inspired architecture. Training objective was accurate prediction of Ki67 in labeled patches and consistency in predictions across repeat unlabeled patches. We measured final model accuracy on held-out test samples. Our promising preliminary results suggest potential for deep learning in deconvolving the spatial heterogeneity of proliferative GBM subpopulations. If successful, this model can provide a non-invasive readout of cell proliferation and reveal the effectiveness of a given cytotoxic therapy dynamically during the patient's routine follow up. Further, the spatial resolution of our approach provides insights into the intra-tumoral heterogeneity of response which can be related to heterogeneity in localization of therapies (e.g. radiation therapy, drug dose heterogeneity).

Abstract Introduction Distinguishing between tumor and treatment effect in post-treatment glioma, although crucial for clinical management, is difficult because contrast-enhancing regions are mixtures of recurrent tumor and reactive tissue, and definitive histopathological criteria do not exist. This study disentangles the marked intra-tumoral heterogeneity in the treatment-recurrent setting by developing an unsupervised framework to algorithmically categorize intraoperative MRI-localized biopsies into three clinically-relevant tissue clusters based on joint analysis of RNA sequencing and histopathological data. Methods A retrospective cohort of 84 MRI-localized biopsies from 37 patients with post-treatment recurrent glioblastoma underwent mRNA extraction and quantification via PLATEseq protocol. For 48 of 84 biopsies, a neighboring piece of tissue underwent quantitative histopathology based on labeling index (LI) for SOX2, CD68, NeuN, Ki67, and H&E. Correlation between LIs and gene expression for these 48 samples was performed. Genes significantly correlated (p&lt;0.05) with ≥1 marker were used for hierarchical clustering of correlation matrix, identifying three mutually-exclusive tissue-specific gene sets. These sets were then used to perform ssGSEA to categorize each of 84 biopsies into one of three tissue types. Results Correlation analysis identified 7779 genes significantly correlated with ≥1 histopathological marker. Clustering revealed three gene sets associated with specific markers: SetA-3688 genes associated with SOX2/Ki67/H&E; SetB-2418 genes associated with CD68; SetC-1673 genes associated with NeuN. ssGSEA using these sets categorized each biopsy into one of three tissue types: 27 biopsies enriched in SetA, 28 in SetB, and 29 in SetC. Conclusions Using MRI-localized biopsies with both RNAseq and histopathological data, this algorithmic approach allows development of three orthogonal tissue-specific gene sets that can be applied to characterize the heterogeneity in post-treatment recurrent glioma: SetA: genes correlated with SOX2/Ki67/H&E, representing recurrent tumor; SetB: genes correlated with CD68 (microglia) representing reactive tissue consistent with treatment effect; SetC: genes correlated with NeuN (neurons), representing infiltrated brain.

Abstract High-grade gliomas (HGGs) nearly always recur after standard initial treatment, and the resulting mixture of recurrent tumor and treatment-induced reactive changes presents major diagnostic challenges. Anatomical imaging, such as MRI, cannot adequately distinguish progressive disease from treatment effect (pseudo-progression). Furthermore, there is marked intra-tumoral heterogeneity, such that some areas of a tumor may demonstrate necrotic treatment effect and others frank recurrence. Due to this difficulty reliably differentiating between these two clinical findings, analytic methods using multiple modalities are necessary to further our understanding of this disease process. To this end, we sought to correlate radiographic, histopathologic and molecular features of surgically sampled post-treatment suspected recurrence to identify markers distinguishing tumor growth from treatment effect. We performed Stimulated Raman Histology (SRH) imaging and highly multiplexed RNA-sequencing (PLATE-seq) on 84 MRI-localized biopsies from 39 patients with clinically suspected recurrent HGG. The SRH images were classified as recurrent tumor or gliotic/reactive tissue using a convolutional neural network trained on an independent cohort including a large set of recurrent HGG, and an automated cell-counting algorithm was used to quantify cellularity from the SRH image of each sample. Differential gene expression analysis of the PLATE-seq data was used to identify gene sets that distinguish recurrent tumor from treatment effect, and single sample gene set variation analysis (GSVA) was used to further assess the molecular and cellular composition of each MRI-localized sample. The histopathologic and molecular features of each sample were also correlated with the MRI features of the corresponding biopsy sites, and this data is currently being used to train machine learning models that predict the distribution of recurrent tumor and treatment-induced reactive changes within a patient’s radiographic lesion. These predictive radiomic models will help to guide neurosurgical sampling, and improve our ability to monitor glioma progression and response to therapy.