Velásquez, Carlos

Abstract Background Glioblastoma (GB) exhibits marked tumor microenvironmental heterogeneity, contributing to both therapy resistance and poor survival outcomes. Intercellular adhesion molecule 1 (ICAM1) is an inflammation-associated adhesion molecule implicated in immune-stromal interactions, but its clinical and biological significance in GB remains incompletely defined. Methods We performed ICAM1 immunohistochemistry, bulk RNA sequencing, DNA methylation analysis, single-nucleus RNA sequencing, and preoperative MRI-derived morphologic feature analysis across multi-institutional cohorts of treatment-naïve IDH-wildtype GB. Results High ICAM1 expression at both transcript and protein levels is associated with shorter overall survival, particularly within the mesenchymal transcriptional subtype. ICAM1 expression is inversely associated with promoter methylation. Single-nucleus RNA-sequencing analysis shows that ICAM1-high tumors are enriched for mesenchymal, hypoxia, stress-associated, and immunosuppressive myeloid programs, suggesting that ICAM1 reflects a multicompartment inflammatory tumor ecosystem. ICAM1 expression, tumor surface irregularity, and patient age are independently prognostic and minimally correlated, capturing complementary clinical, imaging, and molecular features of GB biology. When combined, our Clinical-Imaging-Molecular (CIM) framework refines outcome prediction across independent GB cohorts using routinely available clinical data, conventional imaging, and standard immunohistochemistry. Conclusions ICAM1 characterizes an immune-remodeled, stress-adapted GB ecosystem associated with poor clinical outcomes. Integrating patient age, surface irregularity, and ICAM1 expression yields a scalable and clinically accessible prognostic framework for GB, supporting tailored outcome prediction, particularly in low- and middle-resource settings without routine access to advanced molecular profiling.

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.