Singleton, Kyle

Glioblastoma is the most common form of brain cancer in adults and remains one of the most lethal forms of cancer. Treating glioblastoma is uniquely challenging due to its aggressive nature and its critical location, the brain, which must remain largely undamaged. Due to these difficulties, there is a lack of effective treatment options. Computer models simulating tumors and their response to therapy, known as digital twins, provide a promising new strategy. These model systems are mathematical equations that capture the biological knowledge of how cells interact during tumor growth and in response to treatment. By solving these equations, we can create virtual patients that provide accurate predictions of real patients’ tumor progression. This approach has the potential to improve outcomes by designing therapies for the specific needs of each patient, moving beyond the limitations of the “one-size-fits-all” approach that is currently used.

Abstract High-grade glioma (HGG) continues to have a dismal prognosis owing in part to its heterogeneous and invasive nature both within and between patients. The most aggressive form of high-grade glioma, called glioblastoma, portends a median survival of around 15 months. Like many cancers, high-grade glioma rewires metabolism towards its own benefit, shifting away from the typical reliance of brain tissue on glucose towards other sources such as fatty acids and acetate. Using a cohort of image-localized high-grade glioma biopsies collected from diverse imageable regions with associated transcriptomics (58 patients, 202 biopsies), we are able to assess the spatial transcriptomic landscape of HGG. We have previously shown that the transcriptome of HGG can be characterized into a continuous transition between states, and we have shown that the cellular ecosystem follows a similar pattern, using a combination of trajectory inference and cellular deconvolution. These patterns are strongly associated with the non-enhancing and contrast-enhancing regions of the tumors, which are typically considered for treatment planning and progression assessment. Here, we applied differential expression of hallmark pathways to determine how metabolism factors into this paradigm. As expected, we find overall higher levels of metabolism-associated pathway expression in the contrast enhancing regions compared to the non-enhancing regions. Interestingly, we found metabolic expression differs between the ecologies associated with transcriptional subtypes of high-grade glioma. For example, fatty acid (p=0.003), xenobiotic (p<0.0001) and heme (p<0.0001) metabolism pathways were significantly upregulated for mesenchymal tumor samples, while significantly downregulated for proliferative tumor samples (p=0.02, p<0.0001, p<0.0001, respectively). As we continue to further understand the complex spatial and temporally heterogeneous landscape of these aggressive tumors, it is important to consider the metabolites on which they depend, to better inform prognosis and develop new treatments.

Abstract High-grade gliomas (HGGs) are aggressive, infiltrative brain tumors with poor survival outcomes and limited therapeutic options. Their spatial and molecular heterogeneity necessitates large, well-curated datasets integrating imaging, clinical, and genomic information to enable precision medicine and drive novel treatment development. Using the Medtronic StealthStation neuronavigation system, we collected over 1,000 spatially localized intraoperative biopsy samples from more than 250 patients undergoing resection for glioma at Mayo Clinic and partnering sites. Each biopsy site was recorded with 3D coordinates, labeled in real time, flash-frozen, and annotated with metadata including collection time, anatomical location, and surgical plan. Samples were processed and entered into a secure database along with the corresponding sampling coordinates, frozen pathology, whole-exome and RNA sequencing data, 5-ALA fluorescence status, patient demographics, and T1, T1Gd, T2, and T2-FLAIR MRI. This multimodal dataset supports radiogenomic research by enabling precise mapping of imaging features to molecular characteristics. A subset of this data, including 268 biopsies from 52 patients (25 male, 27 female) with reported intraoperative 5-ALA fluorescence, was used to develop a radiomics model predicting locoregional patterns of 5-ALA positivity based on preoperative MRI. Image features were extracted from the biopsy sample coordinates across multiple MRI sequences and paired with corresponding 5-ALA fluorescence status to train the model for regional fluorescence prediction. The final model achieved 83% accuracy on the training set and 86% on the validation set using a 70-15-15 data split; this included 208 samples for training and 42 for validation, stratified by patient to ensure no overlap of samples from the same individual across sets. These findings demonstrate the feasibility of predicting 5-ALA fluorescence using radiomic features from standard preoperative MRI, offering a promising noninvasive strategy for improving intraoperative decision-making and enabling more precise, personalized resections in patients with high-grade gliomas.

Abstract High grade gliomas (HGGs) are aggressive brain tumors that are difficult to remove completely, leading to poor survival rates. Contrast enhancement on T1Gd magnetic resonance imaging (MRI) has traditionally been used to guide resection, however it often underestimates true tumor burden, adding a layer of uncertainty for surgeons during pre-operative planning. Fluorescence guided surgery using 5-aminolevulinic acid (5-ALA) has emerged as an intraoperative solution and understanding it in conjunction with MRI features can address the limitations of routine T1Gd imaging during surgical planning. Here we investigate sex-specific patterns in HGGs using 5-ALA fluorescence and imaging features. We analyzed a patient dataset of 202 multiregional MRI-localized biopsies from 42 HGG patients (20 male and 22 female). Data included patient demographics, 5-ALA status (positive or negative), MRI annotations (CE = contrast enhancement, NE = non-enhancement, BAT = brain around tumor non-enhancement), and survival outcomes. Statistical analyses included Wilcoxon rank-sum tests, Kruskal-Wallis tests, and chi-squared tests to assess differences in imaging annotations and 5-ALA positive vs negative status among sex cohorts. Among all samples, CE and BAT were associated with higher 5-ALA positivity, while NE was associated with higher 5-ALA negativity (p = 0.02119). Female patients were less likely to have detectable 5-ALA in their multiregional biopsy samples than males (p < 0.0001). CE at the biopsy location was more likely to be indicative of 5-ALA positivity in males than females (p = 0.001927); whereas NE was associated with higher 5-ALA negativity in females than males (p = 0.03426). BAT was associated with higher 5-ALA positivity in males than females (p = 0.03895). Our study reveals sex-specific patterns in 5-ALA fluorescence and image features among HGG patients. These unique features may serve as potential pre-surgical markers for tumor regions that are more likely to benefit from 5-ALA-guided resection.

Abstract The intra- and inter-patient heterogeneity in high-grade glioma (HGG) continues to contribute to its poor prognosis. Clinical biopsies are often harvested from limited regions and typically are not image localized. Thus, they fail to capture the diversity within tumor regions, immune expression or normal cell abundances that play key roles in tumor development. It is important to gain an understanding of these subpopulation ecologies, their spatial resolution, and interactions between them that may then be exploited for future therapeutic benefit. Further, these may differ by patient characteristics such as sex, age at diagnosis and treatment status. Using an ongoing image-localized biopsy collection protocol, we have so far evaluated the bulk transcriptomics of 202 multi-regional biopsies from 58 patients to characterize HGG heterogeneity. These samples were processed through Monocle, a reverse graph embedding algorithm that groups samples into states and orders them along developmental trajectories. Deconvolution methods were previously used to predict relative abundances of 7 normal, 6 glioma, and 5 immune cell subpopulations for each sample, which we have now overlaid on the Monocle graph. Monocle classified HGG into 4 main states along a three-pronged trajectory. These states reveal distinct population ecologies with associated enriched gene pathways. We also note significant immune pathway enrichments that differ between state and patient-reported sex. Together, these algorithms reveal a simple transcriptomic trajectory that helps us understand the development and evolution of HGG. Characterizing the in vivo diversity within and between high grade gliomas is important for understanding prognosis, stratifying future treatments and ultimately improving patient outcome.

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 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 Image-based mathematical modeling is emerging as an advantageous tool to integrate into the care routine for patients with brain tumors. These tools can be used for a variety of tasks including predicting tumor growth, indicating spatial genomic alterations, and guiding radiation therapy. Often, the methods to produce these maps are developed in silos, using cohorts of retrospective patients and switching to a per patient map calculation can be an afterthought. Automatic pipelines are needed to enable creating mathematical maps in a clinically-relevant workflow. To this end, we developed a pipeline to create an image-based mathematical map of edema abundance based on T2-weighted (T2W) magnetic resonance imaging that includes preprocessing steps, abnormality and tissue type segmentation, brain extraction, map generation, and conversion between image file formats. The python-based pipeline was designed to take in T1-weighted (T1W), T1W with gadolinium contrast (T1Gd), T2W, and fluid attenuated inversion recovery (FLAIR) DICOM images for a single patient as input and convert them to NIfTIs for downstream steps. The pipeline preprocessed images by correcting bias field fluctuations using the N4 algorithm and normalizing intensities. The T1W, T1Gd, and FLAIR are then registered to the T2W. T1Gd and FLAIR are used in a deep learning algorithm to delineate the brain and abnormality including enhancing tumor, necrosis, and peripheral edema. Normal tissue segmentations are performed using the statistical parametric mapping (SPM) segmentation routine. The segmentations and preprocessed images are all utilized in the edema mathematical model. The resulting edema map is then reformatted as a DICOM and output for further use by clinicians and clinical imaging systems. Our future goals include integrating our pipeline in a clinical workstation and sending final maps to our hospital’s picture archiving and communication system (PACS) to allow for deployment of these models for patient care in a clinically efficient timeline.