Zadeh, Gelareh

Abstract Brain metastases (BMs) are a common and often fatal progression of systemic cancers, affecting up to 25% of patients. Despite advances in surgical resection and radiotherapy, the median survival remains limited to 10–16 months. While genomic profiling has enabled the development of targeted therapies, there remains a paucity of prognostic tools for clinical use. To address this, we explored the utility of DNA methylation profiling—an epigenetic marker increasingly recognized for its diagnostic and prognostic value. We profiled fresh-frozen tissue from 327 surgically resected BM samples of lung, breast, melanoma, and gastrointestinal (GI) origin using the Illumina Infinium MethylationEPIC BeadChip array. Unsupervised analysis via partitioning around medoids (PAM) clustering identified five robust groups. Four of these correlated closely with primary tumor origin. However, a distinct “Poly-Origin Cluster” emerged, comprising tumors from multiple primary sites with a shared epigenetic signature. Poly-Origin tumors showed improved survival outcomes compared to origin-aligned clusters (20.2 vs. 10.1 months, p=0.0018), independent of clinical covariates. Notably, these tumors were enriched for immune cell infiltration—including CD14+ macrophages, CD19+ B-cells, and CD56+ NK cells—based on deconvolution analysis, and exhibited reduced genomic instability as measured by copy number variation. These tumors were enriched for pathways related to cell signaling. We validated the prognostic significance of the Poly-Origin methylation signature in a publicly available independent cohort (n=96), and found that high signature concordance correlated with extended survival. Finally, we demonstrated the feasibility of detecting this signature in plasma-derived cell-free DNA using cfMeDIP-seq, with high classification accuracy (AUC = 0.98). These findings reveal a clinically meaningful subtypes of BM not captured by tissue origin or mutation status alone. Our work suggests that DNA methylation profiling may be a powerful tissue- and liquid biopsy–based tool that holds promise as a non-invasive prognostic biomarker in the management of brain metastasis.

Abstract A major barrier in understanding adult diffuse gliomas lies in the inability to track tumor evolution in real time. To address this, we employed a novel method known as continuous lineage tracing, which integrates CRISPR/Cas9-based expressed DNA barcoding with single-cell RNA sequencing, enabling a phylogenetic approach to studying tumor development. Patient-derived glioma-initiating cell lines were engineered with expressed barcodes targeted by CRISPR/Cas9 and implanted into mice to create intracranial xenografts. Tumors underwent single-cell RNA sequencing; expressed barcodes were used to infer clonal relationships, and transcriptomic profiles enabled cell state classification. Phylogenetic lineage trees were reconstructed using lineage inference algorithms to characterize cell fitness, expansion, and plasticity. Our analysis uncovered extensive intra-clonal cell state heterogeneity, indicating active phenotype switching prior to therapy. We identified consistent transcriptional programs associated with tumor engraftment and in vivo clonal advantage. Lineage tracing revealed gene expression signatures linked to key phenotypes: fitness, enriched for neural-mesenchymal and injury-response pathways; expansion, associated with RNA splicing; and plasticity, correlated with cell cycle and DNA repair programs. Glioma stem cells appeared to span a transcriptional continuum from undifferentiated, high-fitness states to more differentiated, low-fitness states, with high-fitness cells potentially representing transitional phenotypes. We then validated these findings in a cohort of 185 surgically resected gliomas with matched bulk RNA and proteomic data. Phylogenetic gene signatures differed markedly between IDH-mutant and IDH-wildtype tumors. When stratified by tumor type, both fitness and expansion signatures were significantly associated with overall survival in GBM (median 12.6 months, p=0.041) and oligodendrogliomas (median 66 months, p=0.027). This study demonstrates the utility of continuous lineage tracing to reconstruct tumor evolution and identify transcriptional programs linked to tumor growth and prognosis. Our approach provides a powerful framework for dissecting glioma biology and identifying potential therapeutic vulnerabilities.

Abstract INTRODUCTION Brain metastases (BM) are the most common intracranial tumor, with a poor prognosis of 10-16 months, and lung adenocarcinoma (LUAD) is the most common source. Currently, we cannot reliably predict which patients will develop BMs clinically, so most are detected after they develop and growing to cause deficits. Here, we use tumor DNA methylomes to predict which patients are likely to develop BM through the course of disease and we use both cerebrospinal fluid (CSF) and plasma methylomes for BM liquid biopsy, enabling their early detection and treatment. METHODS DNA methylomes were acquired from 419 tumor and plasma samples from 363 LUAD and BM patients. We built DNA methylation-based models that predict BM development and that also identify BMs non-invasively from CSF and plasma in discovery data subsets. A nomogram was built that predicts BM risk using methylation scores and clinical stage together. Models were evaluated in independent validation sets. RESULTS Our model predicting BM development was accurate (Univariable Cox: HR=5.7, 95% CI=1.9-17.2), p=0.0023), independent of clinical factors (Multivariable Cox: HR=8.9, 95% CI=2.0-40.5, p=0.0046), and more predictive than cancer staging used currently clinically (AUC=0.81 versus 0.65). Our combinatorial nomogram using methylation and staging data together showed enhanced BM prediction (Univariable Cox: HR=17.2, 95% CI=4.1-71.3, p<0.0001; AUROC=0.82). Our models identifying BM through biofluid methylomes accurately distinguished BM from differential diagnoses of gliomas and lymphomas in plasma (AUROC=0.80, 95% CI=0.68-0.93) and with greater accuracy in CSF (AUROC=0.93, 95% CI=0.71-1.0). CONCLUSION Overall, we show the first approaches for robust molecular prediction of BM development and for reliable liquid biopsy of BM. This work is expected to transform cancer management by identifying high-risk lung cancer patients after initial treatment for inclusion in screening protocols using neuroimaging and liquid biopsy, to identify BM development early with less advanced disease and more optimal therapeutic outcomes.

Abstract INTRODUCTION Recurrent IDH-mutant gliomas pose a significant therapeutic challenge, with limited treatment options following progression after standard therapy. Combining PARP inhibitors with immune checkpoint blockade has been proposed as a synergistic strategy in IDH-mutant high-grade gliomas, leveraging vulnerabilities in homologous recombination repair and increased PD-L1 expression following PARP inhibition. OBJECTIVE This phase II trial (NCT03991832) evaluated the combination of the PARP inhibitor olaparib and the PD-L1 inhibitor durvalumab in patients with recurrent IDH-mutant glioma. We also investigated the potential of the plasma tumor methylome as a non-invasive biomarker of treatment response. METHODS Twenty-nine patients (median age 40.5 years; 41% female) were enrolled between January 2020 and February 2023. All patients received olaparib (300 mg twice daily) and durvalumab (1,500 mg IV every four weeks) until radiographic or clinical progression. Plasma samples were collected at baseline and monthly, and cell-free methylated DNA immunoprecipitation sequencing (cfMeDIP-seq) was performed. RESULTS The objective response rate was 10%, and median overall survival was 9.3 months. Longitudinal cfMeDIP-seq profiling enabled development of a circulating methylome classifier that accurately distinguished responders from non-responders. Integration with matched tumor transcriptomic and methylation data revealed enrichment of immune and DNA repair pathways in responders. Whole-exome sequencing identified no consistent mutational correlates. Spatial transcriptomic analysis demonstrated a more interactive, immune-rich tumor microenvironment and reduced malignant cell state diversity in responders. CONCLUSION This study supports the safety of combined PARP and PD-L1 blockade in recurrent IDH-mutant glioma and highlights the plasma methylome as a promising non-invasive biomarker of therapeutic response.

Abstract Meningiomas are the most common primary intracranial tumors and exhibit clinical heterogeneity. Radiotherapy (RT) remains the only adjuvant therapy, but response is variable and biomarkers are limited. RTOG-0539 is the first prospective phase 2 trial to stratify meningioma patients for adjuvant RT based on clinical risk. Here, we apply modern molecular tools to this cohort and uncover correlates of RT response. To do so, tumor tissue from 100 RTOG-0539 patients was profiled using DNA methylation arrays, RNA sequencing, and whole-exome sequencing. Recurrence scores, Molecular Groups, gene expression, and copy number alterations were compared across clinical groups and between RT responders and non-responders. Using this approach, modern grading criteria, including brain invasion, TERT mutation, CDKN2A/B deletion and 1p/1q status, would reclassify 10% of tumors and alter treatment group assignment in 7%. Notably, non-responders to RT had more frequent 1p and 14q loss, and greater copy number burden. In addition, transcriptomic and epigenetic profiling revealed immune-related signatures in RT responders and upregulation of cell cycle–related pathways in RT non-responders, several of which overlapped with targets of vorinostat, a histone deacetylase inhibitor validated in aggressive meningioma models. Of particular importance, the Proliferative Molecular Group was an independent predictor of post-RT recurrence in multivariable analysis, outperforming WHO grade. In conclusion, multi-omic analysis of the RTOG-0539 cohort shows that updated WHO grading criteria, incorporating molecular and cytogenetic features, improve risk stratification. However, molecular classification, particularly the Proliferative group, remains an independent and stronger predictor of RT response. These findings support integrating molecular biomarkers alongside modern grading frameworks to guide treatment and trial design in meningioma.

Abstract CIC-rearranged sarcoma (CRS) is a rare disease driven by a specific fusion protein involving the CIC gene. The occurrence in the brain is 3% in all CRS patients, and these tumors frequently metastasize to the brain. The most common rearrangement is with the double homeobox 4 (DUX4) transcription factor (CIC-DUX4), and others, such as CIC-NUTM1 fusions, have been identified in a subset of pediatric primitive neuroectodermal tumors. However, the molecular mechanisms by which CIC-fusions drive CRS remain unknown. Preliminary data show that CIC-DUX4/NUTM1 fusions activate JAK and its downstream effector STAT1/3. We hypothesize that JAK/STAT1/3 signaling cooperates with CIC-fusions to drive CIC-sarcomas by inducing ETV1/4/5 expression. Patient-derived CRS cell lines showed elevated levels of JAK1/STAT1/3 activation compared to fusion-negative sarcoma lines. Inhibition of JAK1 using Ruxolitinib and Solicitinib reduced STAT1/3 phosphorylation, downregulated ETV1/4/5 expression at both mRNA and protein levels, and diminished ETV5 promoter activity, cell proliferation, and tumorigenicity. Although the mechanism by which CIC-fusions activate oncogenic targets is still under investigation, histone acetylation appears to play a central role. STAT1/3 interacts with p300/CBP to enhance transcription, and STAT1 is necessary for p300 acetyltransferase activity. Ruxolitinib significantly reduced histone acetylation at ETV1/4/5 promoters in hMSC cells expressing CIC-DUX4/NUTM1, as well as in CRS cell lines. Unlike the p300 inhibitor C646, which causes global hypoacetylation, Ruxolitinib’s effects were promoter-specific, indicating its potential as a more targeted and less toxic therapeutic option. Importantly, we found that STAT1/3 binds to ETV1/4/5 promoters only in the presence of CIC-fusions. Luciferase assays confirmed that STAT1/3 alone cannot activate ETV5 transcription without CIC-fusions, revealing a novel cooperative mechanism. In vivo, Ruxolitinib treatment of CRS xenografts led to significant reductions in tumor volume, STAT1/3 activation, and ETV1/4/5 expression. These findings support JAK1/STAT1/3 inhibition as a promising therapeutic strategy for CRS and warrant further preclinical investigation.

Abstract INTRODUCTION While TERT promoter mutation (TPM) has been established as a marker of clinically aggressive meningiomas, this alteration is rare and found in less than 5% of all cases. However, a larger subset of meningiomas may exhibit aberrant TERT expression in the absence of TPMs. This study investigated the effect of TERT gene expression on clinical outcome in meningioma patients. METHODS Clinical and molecular data were retrospectively collected on 1241 meningiomas, split into a Toronto discovery cohort and a multi-institutional validation cohort. Sanger sequencing and bulk RNA sequencing were used to determine TPM status and TERT gene expression. The effect of TERT expression on progression-free survival (PFS) was assessed using Kaplan-Meier and Cox regression analysis. RESULTS While meningiomas with TPM showed expectedly higher TERT gene expression compared to wildtype (TP-WT) cases (p<0.0001), TERT expression was still detected in 28.7% (157/547) of TP-WT meningiomas. Meningiomas with TERT expression showed significantly worse PFS compared to meningiomas without any TERT expression. In fact, WHO grade 1 meningiomas with TERT expression had PFS outcomes resembling WHO grade 2 meningiomas lacking TERT expression (p=0.59). In turn, WHO grade 2 meningiomas with TERT expression had clinical outcomes similar to WHO grade 3 meningiomas without TERT expression (p=0.42). Furthermore, the proportion of meningiomas expressing TERT as well as overall TERT expression levels increased with increasing WHO grade. Multivariable analysis showed that TERT expression was significantly associated with worse PFS even when controlling for other known predictors of clinical outcome including TPM, CDKN2A/B loss, 1p/22q status and WHO grade (HR 1.85 [95% CI 1.33-2.57], p=0.00024). CONCLUSION TERT expression is a novel independent biomarker of outcome for meningiomas identifiable in up to one-third of cases that may be utilized to reclassify tumours to a higher WHO grade.

Abstract BACKGROUND IDH-mutant gliomas are diffusely infiltrating, slow-growing brain tumors that inevitably undergo malignant progression to grade IV astrocytoma - a highly aggressive and uniformly fatal disease with dismal prognosis. Despite extensive research, the molecular mechanisms underlying this transformation remain poorly understood. While alterations in canonical driver genes such as CDKN2A, PTEN, and PDGFRA have been implicated, the broader landscape of genetic drivers involved in IDH-mutant glioma progression is largely uncharacterized. In contrast to many other malignancies, IDH-mutant gliomas acquire relatively few somatic point mutations during progression - typically in the range of 30 to 40 - suggesting that alternative mechanisms may be driving tumor evolution. Notably, these tumors exhibit widespread chromosomal instability, resulting in a heterogeneous and highly prevalent pattern of copy number alterations (CNAs). This observation raises the possibility that gene dosage changes due to CNAs play a critical role in tumor progression. MATERIAL AND METHODS To investigate this hypothesis, we performed integrative analysis of IDH-mutant glioma cohorts from the Mayo Clinic and TCGA datasets (~1,500 patients), identifying twelve highly recurrent CNAs significantly associated with patient outcomes. We refer to these as malignant IDH-associated CNAs (mIDH CNAs). We hypothesized that these recurrent CNAs harbor driver genes that confer selective growth advantages, thereby contributing to their recurrence and prognostic relevance. To systematically assess the functional impact of genes within these regions, we developed a novel in vivo CRISPR-based screening platform, termed CRISPR-Knock Out and Activation Linked Assay (CRISPR-KOALA). This approach utilizes stereotaxic delivery of lentiviral gRNA libraries - targeting mouse orthologs of ~1,800 genes within recurrently deleted regions and ~1,200 genes within amplified regions - into the brains of a genetically engineered mouse model of IDH-mutant glioma (rs557mut;IDHmut;P53KO;AtrxKO;Cas9GFP: Yanchus et. al Science 2022). RESULTS Mice injected with mIDH CNA-targeting libraries developed tumors significantly more rapidly than control cohorts. Deconvolution of the CRISPR-KOALA screen identified several putative driver genes, including regulators of the NOTCH signaling pathway, lipid metabolism modulators such as THEM6, and the NMDA receptor modulator GRINA - implicating dysregulation of metabolism, immune response, and glutamatergic signaling in the malignant transformation of IDH-mutant gliomas. CONCLUSION Together, these findings provide critical insights into the genetic mechanisms driving the progression of IDH-mutant gliomas and identify new candidate genes and pathways for potential therapeutic targeting.

Abstract BACKGROUND Glioblastoma (GBM) is a highly aggressive brain tumour with limited prognostic tools that are broadly implementable in routine clinical settings. Current models often depend on molecular platforms not universally accessible. We aimed to develop and validate a clinically practical prognostic model for GBM that integrates ICAM1 protein expression, imaging-derived tumor morphometrics, and clinical features. MATERIAL AND METHODS We analyzed ICAM1 expression in GBM at the DNA methylation, RNA, and protein levels using public datasets and institutional cohorts, confirming its association with tumor grade and poor survival. ICAM1 protein levels were then quantified via digital pathology in a multi-institutional cohort of 197 newly diagnosed, treatment-naïve GBM patients. Tumor surface irregularity (SIR) was calculated from pre-operative contrast-enhanced MRI using open-source imaging tools. Age at diagnosis was also recorded. These three variables—ICAM1 expression, SIR, and age—were integrated into a Clinical-Imaging-Molecular (CIM) model using multivariable Cox regression. The model was validated in an independent external GBM cohort (n = 20). Survival prediction was assessed using Kaplan-Meier curves, ROC analysis, and concordance index (c-index). RESULTS ICAM1 expression was associated with DNA hypomethylation and significantly correlated with reduced survival at transcriptomic and proteomic levels. Pathway analysis identified associations with immune activation and angiogenesis. In the clinical cohort, high ICAM1 expression, elevated tumor SIR, and older age independently predicted worse survival. The CIM model combining these factors stratified patients into distinct prognostic groups, achieving strong performance (c-index = 0.779, AUC = 0.778). External validation confirmed the model’s predictive utility (median survival difference: 9.1 months, P = 0.009). CONCLUSION We present a validated and clinically applicable prognostic model for GBM that integrates ICAM1 expression, imaging features, and age. The CIM model is scalable, interpretable, and feasible for implementation across academic and community settings, providing a robust tool for patient stratification, treatment planning, and trial design.