Tran, Nhan L.

Abstract Cellular heterogeneity within glioblastoma multiforme (GBM) is a key driver of resistance to standard-of-care therapeutic approaches such as DNA damage (temozolomide, TMZ), angiogenesis inhibition (Bevacizumab, BM), radiation necrosis, and surgical resection. As such, understanding and targeting heterogeneous tumor-associated populations is a critical aim of ongoing GBM research. Central nervous system (CNS) resident microglia constitute significant percentages of primary GBM samples, often up to 50% of bulk tumor. Additionally, glioblastoma-associated microglia have been implicated in intra-tumoral immunosuppression, and are thus a key barrier to successful implementation of immunotherapies that seek to enhance T-cell mediated anti-tumoral cytotoxicity. As such, insight into targetable pathways that upregulate immunosuppressive cytokines within GBM-associated microglia is necessary for overcoming immuno-therapeutic failure. Here, we identify expression of the cytokine TWEAK (TNFSF12) and its receptor Fn14 (TNFRSF12A) within the myeloid compartment of primary GBM samples via scRNA-seq. Furthermore, we demonstrate TWEAK-mediated reprogramming of GBM-associated microglia towards an immunosuppressive phenotype capable of inhibiting CD8+ T cell activation and proliferation. These findings unveil a potent, targetable axis of immunosuppression within GBM and are thus an exciting step towards enhancing efficacy of existing T-cell based anti-cancer strategies for GBM patients. Citation Format: Angad Beniwal, Ryan Eghlimi, Kenneth D. Alexander, Mylan R. Blomquist, Nhan Tran. Role of TWEAK-Fn14 axis in driving immunosuppression in glioblastoma multiforme [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 6442.

Abstract Background Although the epidermal growth factor receptor (EGFR) is a frequent oncogenic driver in glioblastoma (GBM), efforts to therapeutically target this protein have been largely unsuccessful. The present preclinical study evaluated the novel EGFR inhibitor WSD-0922. Methods We employed flank and orthotopic patient-derived xenograft models to characterize WSD-0922 and compare its efficacy to erlotinib, a potent EGFR inhibitor that failed to provide benefit for GBM patients. We performed long-term survival studies and collected short-term tumor, plasma, and whole-brain samples from mice treated with each drug. We utilized mass spectrometry to measure drug concentrations and spatial distribution and to assess the impact of each drug on receptor activity and cellular signaling networks. Results WSD-0922 inhibited EGFR signaling as effectively as erlotinib in in vitro and in vivo models. While WSD-0922 was more CNS penetrant than erlotinib in terms of total concentration, comparable concentrations of both drugs were measured at the tumor site in orthotopic models, and the concentration of free WSD-0922 in the brain was significantly less than the concentration of free erlotinib. WSD-0922 treatment provided a clear survival advantage compared to erlotinib in the GBM39 model, with marked suppression of tumor growth and most mice surviving until the end of the study. WSD-0922 treatment preferentially inhibited phosphorylation of several proteins, including those associated with EGFR inhibitor resistance and cell metabolism. Conclusions WSD-0922 is a highly potent inhibitor of EGFR in GBM, and warrants further evaluation in clinical studies.

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.

Genome-wide fragmentation patterns in cell-free DNA (cfDNA) in plasma are strongly influenced by cellular origin due to variation in chromatin accessibility across cell types. Such differences between healthy and cancer cells provide the opportunity for development of novel cancer diagnostics. Here, we investigated whether analysis of cfDNA fragment end positions and their surrounding DNA sequences reveals the presence of tumor-derived DNA in blood. We performed genome-wide analysis of cfDNA from 521 samples and analyzed sequencing data from an additional 2147 samples, including healthy individuals and patients with 11 different cancer types. We developed a metric based on genome-wide differences in fragment positioning, weighted by fragment length and GC content [information-weighted fraction of aberrant fragments (iwFAF)]. We observed that iwFAF strongly correlated with tumor fraction, was higher for DNA fragments carrying somatic mutations, and was higher within genomic regions affected by copy number amplifications. We also calculated sample-level means of nucleotide frequencies observed at genomic positions spanning fragment ends. Using a combination of iwFAF and nine nucleotide frequencies from three positions surrounding fragment ends, we developed a machine learning model to differentiate healthy individuals from patients with cancer. We observed an area under the receiver operative characteristic curve (AUC) of 0.91 for detection of cancer at any stage and an AUC of 0.87 for detection of stage I cancer. Our findings remained robust with as few as 1 million fragments analyzed per sample, demonstrating that analysis of fragment ends can become a cost-effective and accessible approach for cancer detection and monitoring.

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 The brain is notoriously difficult to access because of its location and protective mechanisms, both internal and external. Imaging-based approaches, such as multiparametric magnetic resonance imaging (mpMRI), provide a method to bypass these obstacles. mpMRI can capture a wide range of data from tumors embedded in brain tissue through various imaging features. However, the biological function captured by each imaging feature remains largely unknown, particularly within the aggressive disease glioblastoma multiforme (GBM). We integrated our GBM SYstems Genetics Network AnaLysis (gbmSYGNAL) gene regulatory network (Plaisier, et al. 2016), based on the multi-omic patient data from The Cancer Genome Atlas (TCGA), with mpMRI-localized patient tumor biopsies. The gbmSYGNAL network was reconstructed within the new image-localized biopsies, and associations revealed a strong connection between the spoiled gradient recalled echo + contrast (SPGR+C) feature and the immune hallmarks of cancer. The SPGR+C method generates an image by employing a contrast agent that bypasses the blood brain barrier (BBB) through leaky blood vessels that have weakened from the force of the rapidly growing GBM tumor. Furthermore, a significant association between SPGR+C signal and the presence of both infiltrating macrophages and T cells indicated that peripheral immune cells may also cross the BBB at these leaky sites. Separately, the mean diffusivity (MD) & T2 features were linked to the sustained angiogenesis hallmark. These features displayed high association with angiogenic genes, such as FLT1, KDR, and PDGFRB, and correlated with the presence of histologically identified vessels. Detection of these disease-relevant functions via mpMRI features has powerful clinical applications, providing a non-invasive approach to GBM diagnosis, prognosis, and treatment direction.

Abstract Brain tumor associated edema is a significant cause of patient morbidity and mortality often requiring continual treatment with differing levels of success across patients. Edema is a broad term that can indicate increased local water as well as diffusely infiltrating tumor cells. Edema is visualized as hyperintense regions on T2-weighted (T2W) magnetic resonance imaging (MRI) and is typically associated with blood-brain barrier (BBB) breakdown and tumor aggressiveness. Understanding the molecular mechanisms driving imaging patterns of “edema” could provide insights into clinical imaging interpretation. We have an ongoing image-guided biopsy study that allows us to link biopsy molecular markers with locoregional MRI patterns of edema. Further, we previously developed a physics-based method to estimate edema abundance (i.e., edema map) from T2W MRIs. Our goal was to identify connections between BBB-associated molecular factors and edema abundance in brain tumors. Our cohort included 38 patients (female: 15, male: 23) with 129 image-guided biopsies (female: 62, male: 67). We correlated image-localized edema map values with the mean transcriptional frequency for 57 genes related to BBB function. Additionally, we examined correlations separately according to patient reported sex (i.e., male and female) and imaging phenotype (i.e., ENH: enhancing and NE: non-enhancing). We utilized multiple comparisons corrections with a 5% false discovery rate to determine significance. For the overall cohort, we observed significant positive correlations for the HIF1A (p< 0.001) and SOX2 markers with edema. For NE samples, significant correlations included APOE (p=0.001), HIF1A (p< 0.001), PIK3CA (p< 0.001), PTCH1 (p< 0.001), and SOX2 (p< 0.001). Amongst female samples, a significant correlation with PTCH1 (p=0.002) was observed. There were no significant correlations noted for male and enhancing sub-cohorts. Significant correlations between molecular markers of BBB and edema map values could lead to clinical biomarkers for edema or tumor aggressiveness.

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 Glioblastoma (GBM) is the most common and aggressive primary brain tumor with a median survival time of 12-15 months and 5-year survival rate of 5%. Despite the application of advanced genetics and biological breakthroughs in GBM, the outcome for GBM remains dismal. Previous genome-wide association studies (GWAS) for GBM shows that the incidence rate for GBM in males is 1.6 times higher in men than women, with men displaying overall worse outcomes. In-vivo studies in gender-matched isogenic murine astrocytic cell lines show that the rate of astrocytic transformation is higher in males than in females, corroborating the GWAS findings in human GBM. Thus, understanding how sex effects GBM biology may provide insight into the mechanism of gliomagenesis and response to therapy. To explore the molecular mechanisms by which transformation of GBM increases in males, we used DESeq2 to assess differentially expressed genes between male and female patients. We found that a long non-coding RNA WNT5A-AS1 is overexpressed in males [log2(fold change) of 1.40; FDR = 0.036], regardless of cofounding factors such as IDH1 status, age, and race. Next, we examined if expression of WNT5A-AS1 is correlated with survival within male patients versus female patients in the TCGA-GBM dataset. In males, WNT5A-AS1’s expression is significantly correlated with patient overall survival (P = 0.042, N=109), but not in female (P = 0.585, N=60). Male patients with high expression of WNT5A-AS1 display poor survival outcomes compared to males with low expression. WNT5A-AS1 is located in the upstream, promoter region of WNT5A. Q-PCR quantification of WNT5A expression increases in male GBM specimens as compared to GBM specimens from females. The role that WNT5A-AS1 plays in regulating WNT5A expression is unclear to date. Interestingly, WNT5A expression shows a positive correlation with patient survival in males, with high expression displaying better survival outcomes (p=0.0272). We have data obtained from a subset of patient samples showing an inverse relationship between WNT5A and WNT5A-AS1 expression. This suggests that there may be an inverse correlation between WNT5A-AS1 and WNT5A expression, which elucidates that WNT5A-AS1 may be regulating WNT5A expression. Current ongoing studies will evaluate the role of WNT5A-AS1/WNT5A in sex-differences of GBM progression. Citation Format: Aunay Miller, Junwen Wang, Joseph Loftus, Nhan L. Tran. Long non-coding RNA WNT5A-AS1 shows sex-dimorphic differences in survival for males with glioblastoma [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 1556.