Analysis of transcription profiles for the identification of master regulators as the key players in glioblastoma

Glioblastoma (GBM) is the most common malignant brain tumor with poor overall survival. Current treatment management for GBM has low efficacy, mainly due to high inter-patient heterogeneity. The transcription profiles in GBM define cell properties essential for tumor progression. We have developed an approach for the identification of master regulators (MRs) that are responsible for the gene expression changes in GBM. The approach is based on transcription factor enrichment analysis with subsequent "upstream" analysis in the signaling network. The main feature of the approach is that all calculations are performed for transcription profiles from individual samples, which allows taking into account GBM transcription heterogeneity. We identified 451 MRs that were up-regulated or down-regulated and, thus, were important parts of positive feedback loops. The number of MRs in the samples correlated with the degree of tumor immune infiltration, while the differences in MR profiles were generally consistent with the known GBM subtypes: mesenchymal, classical, and proneural. MRs densely interact with each other in the signaling network that may be associated with the robustness to pharmacological intervention. We identified 102 receptors among MRs, which is coherent with the importance of cell-cell interactions for GBM progression. The role of some of them in GBM is not currently investigated: lysophosphatidic acid receptors 5 and 6, sphingosine-1-phosphate receptor 4, lysophosphatidylserine receptors GPR34 and GPR174, and G protein-coupled receptors 84 and 132 for fatty acids. Information on the revealed MRs can be used to search for novel therapeutic strategies to treat GBM.