Curtin, Lee

Abstract Magnetic resonance imaging (MRI) is key to clinically managing brain tumor patients, however connecting the biology to imaging remains challenging. We previously developed a two-compartment model of MRI signal intensity to quantitatively estimate relative edema abundance from T2-weighted MRIs. Using this model, we previously identified the fatty acid metabolism (FAM) and oxidative phosphorylation (OxPhos) pathways as sex-distinct, with both pathways amplified for high edema in males and low edema in females. The purpose of this project was to further delineate sex-distinct biology associated with MRI-estimated brain tumor edema abundance. We analyzed 179 multiregional samples (Female: 75; Male: 104) from 55 high grade glioma patients (Female: 21; Male: 34) for bulk RNA-Seq. Patients’ pre-surgical multiparametric MRIs were preprocessed and segmented for abnormal regions and normal tissue. Utilizing the segmentations and preprocessed images we estimated the relative fractions of extracellular and intracellular space based on the edema mathematical model. Samples were characterized by their edema scores, which were analyzed using differential expression for high and low edema, gene set enrichment analysis (GSEA) using MSigDB hallmarks, and leading edge interpretation. Based on transcriptomic leading edge analyses of high and low edema samples from sex-separated cohorts, we identified common and sex-distinct enrichment of the FAM and OxPhos pathways underlying edema patterns. Of the OxPhos pathway leading edge genes, 45 were common, 36 were unique to females, and 66 were unique to males. Of the FAM pathway leading edge genes, 8 were common, 39 were unique to females, and 29 were unique to males. Notably, expression of both IDH1 and IDH2 were increased for males in regions of high edema in the OxPhos pathway. IDH3a was decreased for females in regions of low edema in the OxPhos pathway. These data suggest that there may be sex-distinct metabolism underlying MRI measurable edema formation.

Abstract Magnetic resonance imaging (MRI) is a centerpiece of clinical management of brain tumor patients, yet the biology underlying even commonly interpreted imaging changes is unclear. The purpose of this project was to understand what biology is associated with brain tumor edema abundance and if that association was sex distinct. We have developed a two-compartment model of MRI signal intensity to quantitatively estimate relative edema abundance from T2-weighted MRIs. Multiparametric MRIs were preprocessed with bias field correction and intensity normalization. The images were co-registered and segmented for abnormal regions and normal tissue. The segmentations and preprocessed images were utilized in the edema mathematical model to estimate the fraction of extracellular space (Fecs) and fraction of intracellular space (Fics). Through an ongoing image-localized biopsy collection protocol, we analyzed 179 samples (Female: 75; Male: 104) from 55 patients (Female: 21; Male: 34) for bulk RNA-Seq. Samples were characterized by their edema Fecs scores and underwent differential expression, gene set enrichment analysis using MSigDB hallmarks and cellular deconvolution using CIBERSORTx. Through differential expression and gene set enrichment analyses of these spatial biopsies, we have identified sex-distinct gene expression and pathways corresponding to regions of high vs low edema. Within highly edematous regions, male samples were amplified for respiratory pathways, while female samples were suppressed. Connecting changes on MRIs with molecular markers from image-localized biopsies provides an opportunity to identify the biological drivers of brain tumor associated imageable edema. Specifically, elucidating the sex-distinct patterns connecting changes on MRI with cellular populations and molecular pathways could be used to better interpret imaging of patients and interpreting treatment response. Citation Format: Pamela R. Jackson, Lee Curtin, Sara Ranjbar, Kyle W. Singleton, Maciej M. Mrugala, Richard S. Zimmerman, Bernard R. Bendok, Peter Canoll, Kristin R. Swanson. Sex-distinct patterns of molecular pathways associated with brain tumor edema [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6168.

Abstract It has been widely acknowledged that limited access to the human brain, along with the importance of its structures, presents significant challenges in obtaining brain tumor tissue. This poses difficulties for clinical decision-making, as healthcare professionals often have to rely solely on imaging to address complex problems. There is currently no definitive method to make these decisions other than relying on informed yet subjective human judgment. That said, imaging abnormalities are a result of a presence of particular tumor-related biology. Because of this, and the heterogeneity frequently observed in brain tumors, we have developed an image-localized biopsy protocol that has been employed for research purposes at Mayo Clinic Arizona for over five years. This multidisciplinary approach involves the collaborative efforts from neurosurgeons, neuro-oncologists, scrub-techs, nurses, radiologists, as well as multiple research labs, all working within a well-established pipeline. The implementation of this protocol has allowed us to obtain image-localized biopsies during surgical procedures, enabling the correlation of local imaging findings with corresponding tissue features through the utilization of machine learning models. The primary objective of this endeavor is to enhance our understanding of the underlying biological mechanisms of brain tumors based on patients' imaging data. This protocol has been successful for our research purposes, and with minor adjustments, it can be readily adopted by any hospital in order to collect more data about these incredibly heterogeneous diseases. While the use of neuronavigation during surgery has become commonplace in many medical institutions, localization of tissue samples is not. Utilizing pre-existing infrastructure in the operating room is a small burden with big potential positive outcomes. Sample locations can be useful to clinicians during tumor board discussions, to better assess disease heterogeneity. These minor refinements in the surgical workflow have the capacity to provide clinicians with additional clinically-relevant data without disrupting patient care.

Abstract Glioblastoma (GBM) is the most aggressive and most common primary malignant brain tumor in adults, with a poor median survival time of 15 months. One of the key challenges in treating GBM is its highly heterogeneous nature, with multiple distinct subtypes that have been shown to occur on both inter- and intra-patient levels. Three main classifications, known as classical/proliferative, mesenchymal and proneural have become commonly demonstrated phenotypes. The cell cycle is a fundamental and highly conserved process that controls faithful division of cells; dysregulation of the cell cycle is known to be a key driver in many cancers. However, how the cell cycle is differently regulated between these subtypes has not been well classified in vivo. We investigate these three GBM subtypes using a recently published single nucleus RNAseq (snRNAseq) data set. We compare cell cycle regulation/dysregulation among these three subtypes using Tricycle, an R/Bioconductor package that utilises dimension reduction via principal component analysis and transfer learning to infer cell cycle position from any snRNAseq data set. We find that the classical GBM subtype has the highest proportion of actively dividing cells (cells in: S/G2/M phases), while the mesenchymal and proneural subtypes have a very low proportion of actively dividing cells. This supports the idea of a proliferation-migration dichotomy between GBM subtypes. We use this proportion of actively proliferating cells to calibrate a minimal spatiotemporal mathematical model for GBM tumor growth that accounts for the differences in cell cycle regulation between these three GBM subtypes.

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 Identification of key phenotypic regions such as necrosis, contrast enhancement, and edema on magnetic resonance imaging (MRI) is important for understanding disease evolution and treatment response in patients with glioma. Manual delineation is time intensive and not feasible for a clinical workflow. Automating phenotypic region segmentation overcomes many issues with manual segmentation, however, current glioma segmentation datasets focus on pre-treatment, diagnostic scans, where treat-ment effects and surgical cavities are not present. Thus, existing automatic segmenta-tion models are not applicable to post-treatment imaging that is used for longitudinal evaluation of care. Here, we present a comparison of three-dimensional convolutional neural networks (nnU-Net architecture) trained on large temporally defined pre-treatment, post-treatment, and mixed cohorts. We used a total of 1563 imaging timepoints from 854 patients curated from 13 different institutions as well as diverse public data sets to understand the capabilities and limitations of automatic segmenta-tion on glioma images with different phenotypic and treatment appearance. We as-sessed the performance of models using Dice coefficients on test cases from each group comparing predictions with manual segmentations generated by trained techni-cians. We demonstrate that training a combined model can be as effective as models trained on just one temporal group. The results highlight the importance of a diverse training set, that includes images from the course of disease and with effects from treatment, in the creation of a model that can accurately segment glioma MRIs at multiple treatment time points.

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.