ABSTRACT Glioma intratumoral heterogeneity enables adaptation to challenging microenvironments and contributes to universal therapeutic resistance. Here, we integrated 914 single-cell DNA methylomes, 55,284 single-cell transcriptomes, and bulk multi-omic profiles across 11 adult IDH-mutant or IDH-wild-type gliomas to delineate sources of intratumoral heterogeneity. We found that local DNA methylation instability, or epimutation burden, was elevated in more aggressive tumors, reflected intratumoral variability, linked with transcriptional disruption, and associated with environmental stress response. We show that the activation of cell-state specific transcription factors is impacted by epimutations and that loosened epigenetic control may facilitate cellular plasticity. Our analyses support that somatic copy number alterations (SCNAs) promote epigenetic instability and that SCNAs largely precede epigenetic and transcriptomic diversification during glioma evolution. We confirmed the link between genetic and epigenetic instability by analyzing larger cohorts of bulk longitudinally collected and spatially separated DNA methylation data. Increased DNA methylation instability was associated with accelerated disease progression, and recurrently selected DNA methylation changes were enriched for environmental stress response pathways. Our work provides an integrative framework to better understand glioma evolution and highlights the importance of epigenetic heterogeneity in shaping therapeutic response.

Single-cell multimodal glioma analyses reveal epigenetic regulators of cellular plasticity and environmental stress response