Sosa, Julio

Abstract Glioblastoma (GBM) progression and recurrence are often driven by therapy-resistant cellular states emerging within hypoxic tumor niches. While hypoxia is a hallmark of GBM, the spatial transcriptomic heterogeneity that underpins resistance in these regions remains poorly characterized. Here, we employ pimonidazole labeling in conjunction with high-resolution spatial and single cell transcriptomics across 29 primary, treatment-naive human GBM samples, syngeneic mouse models, and patient-derived glioma-initiating cells exposed to hypoxia to define the hypoxia-specific tumor microenvironment and its associated cellular programs. We reveal spatially organized sub-niches within hypoxic regions, each supporting diverse cellular programs yet coexisting tumor lineage states with minimal transitions between cell states. Phylogenetic reconstruction based on inferred copy number variations demonstrates consistent evolutionary trajectories from non-hypoxic to hypoxic tumor loci. Specific copy number alterations were enriched in hypoxic regions and associated with immune infiltration and poor prognosis in MGMT-methylated tumors. Cross-species analyses confirmed the conservation of hypoxia-associated gene programs and spatial architectures. Together, our study delineates the cellular, evolutionary, and spatial complexity of hypoxia in GBM, identifies hypoxia-induced genomic alterations that may promote tumor fitness, and offers insights into new therapeutic vulnerabilities within hypoxic niches.

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.