Citations

Abstract BACKGROUND High-grade glioma (HGG) is the most common and aggressive form of primary brain cancer. Standard of care surgery employs Magnetic Resonance Imaging to identify the contrast enhancing (CE) core for resection, leaving a non-enhancing (NE) invasive rim of tumor which contributes to recurrence. Recent multi-omics analyses identified predominant neuronal and immune signatures within NE tumor. Here, we apply a graph network to integrate regionally heterogeneous data from a unique dataset of spatially-localized CE and NE biopsies from a large cohort of patients with HGG. METHODS RNA-sequencing of 1159 differentially expressed genes across 159 IDH-wt biopsies (74 patients) from NE/CE regions were represented as nodes in a Neo4j network. Candidate genes for network random walks occurred exclusively within NE communities detected by modularity optimization and localized by cell-type using single-cell RNA-sequencing (scRNAseq). CSF2RA and RBFOX1 were selected as sources (walk length = 5; 10,000 walks/source node). Frequently occurring genes (one-tailed test) in walks were ranked by correlation with CSF2RA/RBFOX1 expression in NE samples and analyzed by gene set enrichment analysis for biological processes (clusterProfiler). Random walk genes were visualized (DoHeatmap) against composite NE/CE scRNAseq profiles. RESULTS CSF2RA and RBFOX1 were upregulated in NE samples (logFC = 0.789, 1.85). Upregulation of CSF2RA, an immune receptor, in the Glycolytic/Plurimetabolic and RBFOX1, a neurodevelopmental regulator, in the Neuronal glioma subtypes were observed in NE samples (ANOVA p < 0.05). CSF2RA walk genes (n=46) enriched metabolic/proliferative signatures, while RBFOX1 walk genes (n=201) showed neurodevelopmental upregulation. scRNAseq localized CSF2RA to monocytes/macrophages while RBFOX1 was expressed in distinct astrocyte/neuroepithelial populations. CONCLUSION CSF2RA and RBFOX1’s gene networks modulate immune and neurodevelopmental processes influencing glioma progression in NE tumor. Understanding the region-specific roles of these networks may offer strategies for targeting cells in the invasive glioma rim.

Abstract BACKGROUND High-grade glioma (HGG) is the most common and aggressive form of primary brain cancer with poor prognosis. Standard of care surgery employs Magnetic Resonance Imaging (MRI) to identify the contrast enhancing (CE) core for resection, leaving behind an unresected non-enhancing (NE) invasive rim of tumor which borders healthy brain tissue and invariably contributes to recurrence. Recent multi-omic analyses identified predominant neuronal and immune signatures within NE tumor; understanding these signatures could provide key insights to combatting recurrence and improving patient outcomes. To this end, we utilized a novel graph database method to integrate multi-omic data from a unique dataset of multi-regional spatially-localized biopsies from the CE core and invasive NE rim from a large cohort of patients with HGG. MATERIAL AND METHODS 159 IDH-wildtype tumor biopsy samples from the NE and CE regions of 74 glioma patients were profiled by bulk RNA sequencing. Differentially expressed genes (DEG) (n=1159) across NE and CE samples were selected using EdgeR. Samples were hierarchically clustered (method = median) for each DEG and represented as nodes in a Neo4j graph network. Edges were generated between nodes sharing ≥ 1 sample. Modularity optimization was used to identify communities of densely connected nodes. Genes occurring exclusively within NE communities were localized by cell type using external single-cell RNA sequencing datasets and selected as candidate source nodes for NE network random walks. RBFOX1 was selected as a source (walk length = 5; 10,000 walks/source node). Random walks traversing samples with significantly elevated RBFOX1 expression (one-tailed test) were kept for analysis. Genes occurring at significantly high frequency (one-tailed test) within random walks were ranked by correlation with RBFOX1 expression in NE samples and analyzed by gene set enrichment analysis (GSEA) for biological processes using clusterProfiler. RESULTS RBFOX1, a neurodevelopmental regulator, was upregulated (logFC = 1.85) in NE samples and uniquely occurred in NE communities. We show RBFOX1 is differential expressed in the aggressive Neuronal glioma subtype and enriched in neuroepithelial/astrocyte subpopulations across single cell HGG datasets. NE network random walks identified 201 DEGs associated with RBFOX1. GSEA revealed upregulation of neurodevelopmental signatures and downregulation of proliferative pathways amongst RBFOX1’s network. CONCLUSION RBFOX1’s genetic network within NE tumor represents neurodevelopmental processes potentially underlying glioma invasiveness within the NE region. As such, understanding the mechanistic role of key genes within the RBFOX1 network may provide insight into novel strategies for targeting these invasive cells within an aggressive subtype of glioma in the NE rim.

Glioblastoma (GBM) is the most common primary malignant brain tumor in adults, with few effective treatments. EGFR alterations, including expression of the truncated variant EGFRvIII, are among the most frequent genomic changes in these tumors. EGFRvIII is known to preferentially signal through STAT5 for oncogenic activation in GBM, yet targeting EGFRvIII has yielded limited clinical success to date. In this study, we employed patient-derived xenograft (PDX) models expressing EGFRvIII to determine the key points of therapeutic vulnerability within the EGFRvIII-STAT5 signaling axis in GBM. Our findings reveal that exogenous expression of paralogs STAT5A and STAT5B augments cell proliferation and that inhibition of STAT5 phosphorylation in vivo improves overall survival in combination with temozolomide (TMZ). STAT5 phosphorylation is independent of JAK1 and JAK2 signaling, instead requiring Src family kinase (SFK) activity. Saracatinib, an SFK inhibitor, attenuates phosphorylation of STAT5 and preferentially sensitizes EGFRvIII+ GBM cells to undergo apoptotic cell death relative to wild-type EGFR. Constitutively active STAT5A or STAT5B mitigates saracatinib sensitivity in EGFRvIII+ cells. In vivo, saracatinib treatment decreased survival in mice bearing EGFR WT tumors compared to the control, yet in EGFRvIII+ tumors, treatment with saracatinib in combination with TMZ preferentially improves survival.

Abstract Glioblastoma (GBM) is the most common and lethal primary brain tumor in adults, despite current treatment options. One of the major contributing factors to this dismal outcome is a subset of cells called glioma stem cells (GSCs) that are highly tumorigenic, invasive, and resistant to conventional treatments. Thus, studies that focus on the properties of GSC biology are needed to develop novel therapies to improve overall survival. Here we report a novel role of STAT5 in GSC stemness and survival. In GBM tumors, STAT5 has been shown to be active in the infiltrative, residual population of cells remaining post-surgery. Knockdown of STAT5 results in the decrease in GSC stemness and therapeutic resistance. Interestingly, we find that expression of OLIG2, a key cell fate factor important for the tumorigenic potential of GSCs, affects STAT5 activation with a positive correlation between OLIG2 and STAT5 mRNA expression in GBM tumors. Furthermore, we find that in GSCs expressing EGFRvIII, STAT5 interacts with OLIG2 and pharmacological inhibition of STAT5 activation or knockdown of STAT5 expression results in decrease of OLIG2-mediated GSC invasion. Together, these data suggest a novel transcriptional complex, STAT5-OLIG2, in the potential role for GSC survival, maintenance, and invasion. Citation Format: Aunay Miller, Costanza Lo Cascio, Mylan Blomquist, James B. McNamara, Keely Orndorff, Matthew Dufault, Christopher Sereduk, Joseph Loftus, Shwetal Mehta, Nhan L. Tran. A novel transcriptional complex of OLIG2-STAT5 mediates glioma stem cell self renewal and invasion [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 256.

Abstract Glioblastoma (GBM) is the most common and lethal primary brain tumor in adults, despite current treatment options. One of the major contributing factors to this dismal outcome is a subset of cells called glioma stem cells (GSCs) that are highly tumorigenic, invasive, and resistant to conventional treatments. Thus, studies that focus on the properties of GSC biology are needed to develop novel therapies to improve overall survival. Here we report a novel role of STAT5 in GSC stemness and survival. In GBM tumors, STAT5 has been shown to be active in the infiltrative, residual population of cells remaining post-surgery. Knockdown of STAT5 results in the decrease in GSC stemness and therapeutic resistance. Interestingly, we find that expression of OLIG2, a key cell fate factor important for the tumorigenic potential of GSCs, affects STAT5 activation with a positive correlation between OLIG2 and STAT5 mRNA expression in GBM tumors. Furthermore, we find that in GSCs expressing EGFRvIII, STAT5 interacts with OLIG2 and pharmacological inhibition of STAT5 activation or knockdown of STAT5 expression results in decrease of OLIG2-mediated GSC invasion. Together, these data suggest a novel transcriptional complex, STAT5-OLIG2, in the potential role for GSC survival, maintenance, and invasion.

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 Glioblastoma (GBM) tumors contain a subpopulation of glioma stem cells (GSC) which can escape therapeutic effects through self-renewal and invasion. Epidermal growth factor receptor (EGFR) mutations are the most common alteration in GBM. EGFRvIII (deletion of exons 2-7) is observed in GBM (20-25%) resulting in a constitutively activated EGFR receptor. In addition, extracellular domain missense point mutations, (e.g. G598V, A289V, R108K) are observed in GBM (14.4%). These missense mutations in the EGFR extracellular domain also induces a highly active receptor by enhancing ligand binding. Earlier reports have shown that these EGFR missense point mutations have exhibited therapeutic resistance; however, the molecular mechanism of these missense point mutations in GBM remains unknown. Previous studies in our lab showed that EGFRvIII promotes GSC maintenance and therapeutic resistance via oligodendrocyte transcription factor 2 (OLIG2) signaling. OLIG2 is universally expressed in glioma and essential for reprogramming GBM cells into stem-like tumor cells. Increased OLIG2 expression promotes invasion, growth, proliferation, and survival. Here we report an increase in OLIG2 protein expression in GSCs harboring EGFR-G598V or EGFR-A289V mutations, suggesting these mutations may support stemness properties in GBM tumors. GSC expressing these point mutations display increased self-renewal, growth, and survival. In addition, EGFR-G598V and EGFR-A289V trigger distinct downstream signaling pathways in GSCs to regulate OLIG2 function. Hence, we hypothesize that GBM tumors harboring EGFR missense mutations (G598V, A289V) influence stem cell properties via multiple pathways to evade therapeutic effects and inhibiting these different signaling pathways will decrease therapeutic resistance in GBM. Understanding the mechanism of these EGFR mutations may lead to improved therapeutic strategies for GBM patients.

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 are associated with lower grade tumors and portend a favorable prognosis relative to IDH-wildtype glioma. However, molecular profiling has identified additional alterations in lower grade gliomas that confer an aggressive high-grade phenotype and a worse clinical prognosis. Of these, the loss of CDKN2A/B in IDH-mutant astrocytoma is the strongest implicated alteration and is sufficient to define a grade 4 tumor regardless of histologic grade. However, there remain no effective therapies targeted at molecular subgroups in aggressive gliomas to date. Here, we sought to elucidate the biologic pathways and functional outcomes associated with the acquisition of high-grade behavior under loss of CDKN2A/B in IDH-mutant astrocytoma to inform future therapeutic strategies. We analyzed a cohort of patient IDH-mutant astrocytomas with whole exome and RNA sequencing data and found several candidate genes/pathways regulated concurrent with CDKN2A/B loss, including regulators of embryonic development, extracellular matrix organization and regulation of vascular processes. Moreover, we identified 1 tumor that displayed intratumoral heterogeneity of CDKN2A/B expression across biopsies. Using grade 4 IDH-mutant astrocytoma patient-derived (PDX) and established cell-lines harboring CDKN2A/B homozygous deletion, we then stably re-expressed cell cycle regulators p14, p15, and p16 using a Tet-inducible lentiviral system and confirmed their decreased proliferation as compared to parental lines. IDH-mutant cell lines with re-expression of p14, p15 or p16 were characterized by differential transcription factor activation and distinct vascular formation patterns as compared to parental lines. These lines also displayed an increased cell migration and invasive potential and were associated with enhanced integrin signaling through integrin beta-3. Together, our work supports a hypothesis of regional tumor evolution with distinct biologic phenotypes across IDH-mutant gliomas.

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.