Malignant brain tumor-glioblastoma is not only difficult to treat but also hard to study and model. One of the reasons for these is their heterogeneity, i.e., individual tumors consisting of cancer cells that are unlike each other. Such diverse cells can thrive due to the simultaneous co-evolution of anatomic niches and adaption into zones with distorted homeostasis of oxygen. It dampens cytotoxic and immune therapies as the response depends on the cellular composition and its adaptation to hypoxia. We explored what transcriptome reposition strategies are used by cells in the different areas of the tumor. We created the hypoxic map by differential expression analysis between hypoxic and cellular features using RNA sequencing data cross-referenced with the tumor's anatomic features (Ivy Glioblastoma Atlas Project). The molecular functions of genes differentially expressed in the hypoxic regions were analyzed by a systematic review of the gene ontology analysis. To put a hypoxic niche signature into a clinical context, we associated the model with patients' survival datasets (The Cancer Genome Atlas). The most unique class of genes in the hypoxic area of the tumor was associated with the process of autophagy. Both hypoxic and cellular anatomic features were enriched in immune response genes whose, along with autophagy cluster genes, had the power to predict glioblastoma patient survival. Our analysis revealed that transcriptome responsive to hypoxia predicted worse patients' outcomes by driving tumor cell adaptation to metabolic stress and immune escape.
Keywords: autophagy; glioblastoma; hypoxia; immune response; transcriptome.