Smith, Kris A.

Abstract High-grade glioma continues to have dismal survival owing in part to its intra- and inter-patient heterogeneity. Standard clinical protocol collects tumor samples with the aim of providing or confirming a diagnosis and determining the status of a few key genes (e.g. IDH1, MGMT). However, this protocol is unable to capture the diversity within tumor regions, immune expression or normal cell abundances that play key roles in the development of the disease. To overcome this, during surgery we collect image-localized multi-regional biopsies to characterize this disease heterogeneity. Data collection is ongoing, and we currently have 202 samples from 58 patients with available bulk RNA-Seq. With a single-cell reference dataset from Columbia University, we used CIBERSORTx, a deconvolution method, to predict relative abundances of 7 normal, 6 glioma, and 5 immune cell subpopulations for each sample. We used Cox proportional hazard models with first-order statistics (mean, minimum, maximum) of abundances within patients to determine whether these are prognostic, then used TCGA RNA-Seq data to validate these findings. We found that one glioma proneural subtype was significantly beneficial for patient survival relative to other glioma subtypes across all statistics (and showed significance in TCGA). Proliferative and mesenchymal glioma subtypes also showed significance for one or two of the calculated statistics. Oligodendrocyte progenitor cell abundances were consistently significantly beneficial for patient survival in our cohort, while an increased abundance of abnormal (reactive) astrocytes is associated with poor prognosis. We also ran these analyses within invasive margin and core tumor regions to determine location-specific population abundances associated with patient survival. In conclusion, understanding the in vivo diversity of cellular subpopulations within high grade glioma is important for treatment stratification and patient benefit.

Abstract AIMS High-grade glioma continues to have dismal survival with current standard-of-care treatment, owing in part to its intra- and inter-patient heterogeneity. Typical diagnostic biopsies are taken from the dense tumor core to determine the presence of abnormal cells and the status of a few key genes (e.g. IDH1, MGMT). However, the tumor core is typically resected, leaving behind possibly genetically, transcriptomically and/or phenotypically distinct invasive margins that repopulate the disease. As these remaining populations are the ones ultimately being treated, it is important to know their compositional differences from the tumor core. We aim to identify the phenotypic niches defined by the relative composition of key cellular populations and understand their variation amongst patients. METHOD We have established an image-localized research biopsy study, that samples from both the invasive margin and tumor core. From this protocol, we currently have 202 samples from 58 patients with available bulk RNA-Seq, collected between Mayo Clinic and Barrow Neurological Institute. Using a single-cell reference dataset from our collaborators at Columbia University, we used CIBERSORTx, a deconvolution method, to predict relative abundances of 7 normal, 6 glioma, and 5 immune cell states for each sample. RESULTS We find that these cell state abundances connect to patient survival and show regional differences. For example, proneural glioma states were higher in invasive regions, whereas proliferative and mesenchymal states were higher in the tumor core. CONCLUSIONS Our analysis demonstrates a need to characterize the residual tissue following glioma resection to better understand the recurrent disease.

Abstract High grade glioma (HGG) represents a group of devastating diseases with dismal prognosis. Surgical resection of the contrast enhancing (CE) region of HGG remains the mainstay of treatment, but recurrence inevitably arises from the unresected non-contrast enhancing (NE) region, surgically inaccessible due to cancer cell invasion into healthy brain tissue. Due to its critical role in recurrence, understanding of the NE region is central to the improvement of clinical outcomes. We reveal the biological characteristics of this region through image localized multi-regional sampling. We linked microenvironmental characteristics measured by multi-parametric MRI to genomic mutations and transcriptional phenotypes using mixed effect modeling which allowed us to control for individualized patient effects. We first confirmed that T2 is a significant indicator of IDH mutation status in the NE region, being the first description of such a relationship in a HGG cohort. We found the combination of EGFR amplification and CDKN2A homozygous loss was associated with a significantly lower mean diffusivity (MD) compared to double wild type tumors in the NE region, indicating the presence of greater cellular packing and proliferation in EGFR amplification/CDKN2A loss regions. Finally, using single cell pathway based tumor classifications, we showed that nK2, a DSC-MRI metric representing cell size heterogeneity, correlated positively with neuronal signature and negatively with glycolytic/plurimetabolic signature within the NE tumor, indicating that glycolytic/plurimetobolic tumors possessed a high amount of cell size heterogeneity compared to neuronal samples. This hypothesis was supported using digital reference object (DRO) modeling which confirmed that cell size and heterogeneity drove the differential nK2 signal between neuronal and glycolytic/plurimetabolic samples. We identified immune cell infiltrate as one possible mechanism of increased cell size heterogeneity using transcriptomic signature analysis which found more immune cell signatures within glycolytic/plurimetobolic tumors compared to neuronal. Collectively this study demonstrates the central role of multi-parametric MRI as a non-invasive measure of tumor biology and a tool for understanding the clinically critical NE region which can then inform new therapies targeting this region of HGG recurrence. Citation Format: Matthew Flick, Taylor Weiskittel, Kevin Meng-Lin, Fulvio D'Angelo, Francesca Caruso, Shannon Ensign, Mylan Blomquist, Luija Wang, Christopher Sereduk, Gustavo De Leon, Ashley Nespodzany, Javier Urcuyo, Ashlynn Gonzalez, Lee Curtin, Kyle Singleton, Aliya Anil, Natenael Simmineh, Erika Lewis, Teresa Noviello, Reyna Patel, Panwen Wang, Junwen Wang, Jennifer Eschbacher, Andrea Hawkins-Daarud, Pamela Jackson, Kris Smith, Peter Nakaji, Bernard Bendok, Richard Zimmerman, Chandan Krishna, Devi Patra, Naresh Patel, Mark Lyons, Matthew Neal, Kliment Donev, Maciej Mrugala, Alyx Porter, Scott Beeman, Yuxiang Zhou, Leslie Baxter, Christopher Plaisier, Jing Li, Hu Li, Anna Lasorella, Chad Quarles, Kristin Swanson, Michele Ceccarelli, Antonio Iavarone, Nhan Tran, Leland Hu. Multi-parametric MRI maps regional heterogeneity of high grade glioma phenotypes. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5621.

Abstract High grade gliomas (HGG) are aggressive primary brain malignancies typified by diffuse invasion, genetic heterogeneity, and a universally fatal outcome. MRI-defined contrast-enhancing (CE) tumor burden serves as the clinical standard that guides maximal surgical resection and post-therapy response assessment. However, HGGs also comprise an invasive non-enhancing (NE) tumor margin that extends beyond the CE core and harbors the cells that contribute to recurrence. Sampling restrictions have hindered the comprehensive study of these NE HGG cell populations driving tumor progression. Herein, we present an integrated multi-omic analysis of 313 spatially matched multi-regional CE and NE tumor biopsies from 68 HGG patients, performing whole exome and RNA sequencing of both IDH wild-type and IDH mutant HGGs. We report spatially restricted molecular profiles in IDH-mutant HGG, highlighting a concern for sampling bias given the importance of molecular diagnosis and prognostication in IDH-mutant HGG. Regardless of IDH status, we found that NE tumor regions harbored the highest proportion of private mutations, which reflects an increased development of regional genomic complexity in infiltrative tumor. The multiregional genomic profiling of our IDH wild-type HGG cohort reveals that EGFR and NF1 somatic alterations occur as mutually exclusive events in 98.7% of tumors. However, we resolved rare low allele frequency co-alterations of EGFR and NF1 within the NE region. We find this co-occurrence enriched in recurrent tumors, pointing to the early emergence of NF1 inactivation in the NE regions. We constructed genomic models predictive of recurrent disease from both NE and CE regions, which highlight the occurrence of clonal EGFR copy number alterations and NF1 loss as clonal or subclonal events, respectively, emphasizing the regional and temporal complexity of well-studied canonical driver alterations. We detailed the spatially unique acquisition of multiple distinct EGFR alterations giving rise to intra-tumoral EGFR mosaicism, a challenge in the implementation of EGFR directed therapies. Our study also identified two transcriptomic clusters delineated by the significant overrepresentation of neuronal (NEU) and glycolytic/plurimetabolic (GPM) pathway-based functional states in the NE region. NE regions of the NEU subtype harbor the greatest proportion of private mutations, suggesting these infiltrative tumor cells accumulate alterations without clonal expansion. GPM populations conversely displayed a less branched phylogeny and were transcriptionally enriched in immune cell signatures. This phenotypic dichotomy between GPM and NEU populations supports the growing body of evidence that invasive GBM cells either take on a neuronal phenotype for active invasion or a more metabolic phenotype involving interaction with astrocytes, other glial cells, and infiltrating immune cells. Citation Format: Mylan R. Blomquist, Leland S. Hu, Fulvio D'Angelo, Taylor M. Weiskittel, Francesca P. Caruso, Shannon P. Fortin Ensign, Christopher Sereduk, Gustavo De Leon, Lee Curtin, Javier Urcuyo, Ashlynn Gonzalez, Ashley Nespodzany, Teresa Noviello, Jennifer M. Eschbacher, Kris A. Smith, Peter Nakaji, Bernard R. Bendok, Richard S. Zimmerman, Chandan Krishna, Devi Patra, Naresh Patel, Mark Lyons, Kliment Donev, Maciej Mrugala, Alyx Porter, Anna Lasorella, Kristin R. Swanson, Michele Ceccarelli, Antonio Iavarone, Nhan L. Tran. Multiregional sampling of high grade glioma identifies regional biologic signatures [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 1507.

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 Treatment options for glioblastoma remain limited, particularly for those who are not eligible for traditional resection, whether due to lesion location or inability to tolerate open craniotomy. Maximal-safe resection followed by radiation with concurrent and adjuvant temozolomide offers the best outcomes for patients. Unfortunately, not all tumors are amenable to conventional surgical resection at the time of diagnosis with only about 1/3 of patients able to receive a gross-total resection and 15-25% of patients receiving biopsy only, thus reducing their projected overall survival to 9 months. Laser interstitial thermal therapy (LITT) is a minimally invasive, cytoreductive tool, that has demonstrated safety as a surgical approach to treat primary brain tumors. METHODS Data from LAANTERN prospective multicenter registry (NCT02392078) was analyzed to determine clinical outcomes for patients with new and recurrent IDH wild-type glioblastoma (N=89). Demographics, intraprocedural data, adverse events, KPS, health-economics, and survival data were prospectively collected then analyzed separately for newly diagnosed GBM (N=29) and recurrent GBM (N=60). RESULTS Median overall-survival was 9.73 months (95% CI: 5.16, 15.91) for newly diagnosed patients and median post-procedure survival was 8.97 (6.94, 12.36) months for recurrent patients. Median overall-survival for newly diagnosed patients receiving post-LITT chemoradiotherapy was 16.14 months (6.11, not reached). The median length of hospital stay was 50 hours and 80% of patients were discharged to home. CONCLUSIONS LITT offers an effective cytoreductive approach for patients with newly diagnosed and recurrent IDH wild-type glioblastoma. Importantly, its use in newly diagnosed patients who receive post-LITT chemoradiotherapy leads to a median OS similar to that of patients treated with conventional surgical resection. LITT remains an important alternative for patients with inoperable tumors or those not amenable to resection. Enrollment in LAANTERN is ongoing and these cohorts will be revisited as data continues to mature. Benefits beyond cytoreduction are also being actively explored.

Abstract High-grade gliomas represent the most common type of primary adult malignant brain tumor historically diagnosed and graded from histologic criteria alone. Gliomas harboring isocitrate dehydrogenase (IDH) 1 or 2 mutations, which are present in more than 80% of World Health Organization (WHO) grade 2 or 3 tumors, portend a favorable prognosis as compared to IDH-wildtype tumors. However, recent molecular profiling has identified the presence of additional alterations in these lower grade tumors that transform the biology toward an aggressive high-grade phenotype clinically that are independent of histologic grading and confer a worse clinical prognosis. Of these, the loss of CDKN2A/B is the strongest implicated alteration and is sufficient to define a grade 4 tumor in IDH-mutant astrocytomas regardless of histologic appearance. However, there remain no effective therapies targeted at molecular subgroups in aggressive gliomas to date. Here we have profiled the intratumoral heterogeneity via whole exome and RNA sequencing and copy number alteration analysis within a cohort of de-identified IDH-mutant astrocytoma patient tumors collected through the Mayo Clinic and Barrow Neurological Institute (BNI) biorepositories to determine the spatial evolution of aggressive molecular features. Within this cohort we have collected 42 image-localized intratumoral biopsy specimens representing 13 IDH-mutant patient tumors. We identified 11 biopsies from 3 tumors that harbored loss of CDKN2A/B which included 1 tumor that displayed intratumoral heterogeneity of CDKN2A/B expression across subclones. We performed an unsupervised cluster analysis to profile the transcriptional signatures in IDH-mutant astrocytomas across both differential CDKN2A/B status and tumor geographic location comparing contrast enhancing versus non-enhancing regions, and we analyzed the MRI features in IDH-mutant tumors across CDKN2A/B status. We found that tumors with CDKN2A/B alteration demonstrated lower signal intensity on delayed post-contrast T2*W images, which could represent lower pre-contrast T2*, unique contrast agent kinetics, or dissimilar diffusion properties related to increased cellular density and/or greater heterogeneity of cell size. Interestingly, we uncovered spatial heterogeneity of additional molecular drivers of high-grade glioma including KRAS, MYCN, and PDGFRA. While EGFR alterations are rarely reported in IDH-mutant tumors, we identified one treatment naïve patient tumor that harbored EGFR alteration only in the contrast non-enhancing region. Based on this work we propose there is early spatial evolution of aggressive tumor subclones prior to primary tumor diagnosis, and multiregional molecular profiling may better inform tumor biology and predict therapeutic strategies in IDH-mutant astrocytomas. Citation Format: Shannon P. Fortin Ensign, Fulvio D'Angelo, Francesca P. Caruso, Christopher Sereduk, Kyle Singleton, Jennifer Eschbacher, Gina Mazza, Maciej Mrugala, Alyx Porter, Bernard Bendok, Richard Zimmerman, Kris Smith, Peter Nakaji, Kamala Swanson, Junwen Wang, Nathenael Semmineh, Christopher Quarles, Anna Lasorella, Kristin Swanson, Michele Ceccarelli, Antonio Iavarone, Leland Hu, Nhan L. Tran. Multiregional imaging and genomic analysis in high-grade IDH-mutant astrocytoma identifies intratumor subclones with aggressive features [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 770.

Abstract Background Treatment options for unresectable new and recurrent glioblastoma remain limited. Laser ablation has demonstrated safety as a surgical approach to treating primary brain tumors. The LAANTERN prospective multicenter registry (NCT02392078) data were analyzed to determine clinical outcomes for patients with new and recurrent IDH wild-type glioblastoma. Methods Demographics, intraprocedural data, adverse events, KPS, health economics, and survival data were prospectively collected and then analyzed on IDH wild-type newly diagnosed and recurrent glioblastoma patients who were treated with laser ablation at 14 US centers between January 2016 and May 2019. Data were monitored for accuracy. Statistical analysis included individual variable summaries, multivariable differences in survival, and median survival numbers. Results A total of 29 new and 60 recurrent IDH wild-type WHO grade 4 glioblastoma patients were treated. Positive MGMT promoter methylation status was present in 5/29 of new and 23/60 of recurrent patients. Median physician-estimated extent of ablation was 91%-99%. Median overall survival (OS) was 9.73 months (95% confidence interval: 5.16, 15.91) for newly diagnosed patients and median post-procedure survival was 8.97 months (6.94, 12.36) for recurrent patients. Median OS for newly diagnosed patients receiving post-LITT chemo/radiation was 16.14 months (6.11, not reached). Factors associated with improved survival were MGMT promoter methylation, adjuvant chemotherapy within 12 weeks, and tumor volume <3 cc. Conclusions Laser ablation is a viable option for patients with new and recurrent glioblastoma. Median OS for IDH wild-type newly diagnosed glioblastoma is comparable to outcomes observed in other tumor resection studies when those patients undergo radiation and chemotherapy following LITT.

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 Neuroimaging dogma for glioblastoma asserts that hyperintensity on T1Gd MRI reveals the bulk of the tumor, while T2/FLAIR signal indicates edema. However, it is unclear whether this edema results from immune response or increased tumor cells. Further, one significant driver of the known sex differences in glioblastoma may be differences in immune response, due to the X-linkage of many immune genes. Based on this, we hypothesized that assumptions regarding tumor cellularity in T2/FLAIR images should be tailored to the biological sex of the patient. METHODS Using a retrospective cohort of 18 primary glioblastoma patients receiving multiple image-localized biopsies (82 total) and standard MRI, we assessed: distance of biopsy from T1Gd and T2 areas; a pathologist’s score of percent tumor cell density; and an imaging-based invasion metric, D/ρ. This metric is derived from the biomathematical Proliferation-Invasion model of glioma growth, which features two parameters, net growth rate (ρ) and net invasion rate (D). Their ratio D/ρ is related to degree of invasion, and is estimated from volumetric measurements of MRI abnormalities. Additionally, 25 patient-derived xenograft models implanted in females were grown until moribund, at which point brains were excised and stained for DAPI (to show all cells) and Lamin (to highlight tumor cells). Image processing of lamin-stained sections defines contours of intensity correlating with cell density. RESULTS Outside both the T1Gd and T2 region, male patient biopsies had higher tumor cell densities than females. Males also tended to have higher invasion metrics. Although each set derived from different patients, preclinical metrics of invasion were positively correlated with clinical invasion in females but negatively correlated in males. CONCLUSION Our preliminary finding that cell distribution patterns correlate with imaging metrics differently between the sexes supports the hypothesis that the degree of tumor cell density represented on certain MRI sequences may be sex-specific.