Abstract The intra- and inter-patient heterogeneity in high-grade glioma (HGG) continues to contribute to its poor prognosis. Clinical biopsies are often harvested from limited regions and typically are not image localized. Thus, they fail to capture the diversity within tumor regions, immune expression or normal cell abundances that play key roles in tumor development. It is important to gain an understanding of these subpopulation ecologies, their spatial resolution, and interactions between them that may then be exploited for future therapeutic benefit. Further, these may differ by patient characteristics such as sex, age at diagnosis and treatment status. Using an ongoing image-localized biopsy collection protocol, we have so far evaluated the bulk transcriptomics of 202 multi-regional biopsies from 58 patients to characterize HGG heterogeneity. These samples were processed through Monocle, a reverse graph embedding algorithm that groups samples into states and orders them along developmental trajectories. Deconvolution methods were previously used to predict relative abundances of 7 normal, 6 glioma, and 5 immune cell subpopulations for each sample, which we have now overlaid on the Monocle graph. Monocle classified HGG into 4 main states along a three-pronged trajectory. These states reveal distinct population ecologies with associated enriched gene pathways. We also note significant immune pathway enrichments that differ between state and patient-reported sex. Together, these algorithms reveal a simple transcriptomic trajectory that helps us understand the development and evolution of HGG. Characterizing the in vivo diversity within and between high grade gliomas is important for understanding prognosis, stratifying future treatments and ultimately improving patient outcome.
Read full abstract