Molecular landscapes of glioblastoma cell lines revealed a group of patients that do not benefit from WWOX tumor suppressor expression

INTRODUCTION: Glioblastoma (GBM) is notorious for its clinical and molecular heterogeneity, contributing to therapeutic failure and a grim prognosis. WWOX is one of the tumor suppressor genes important in nervous tissue or related pathologies, which was scarcely investigated in GBM for reliable associations with prognosis or disease progression despite known alterations. Recently, we observed a phenotypic heterogeneity between GBM cell lines (U87MG, T98G, U251MG, DBTRG-05MG), among which the anti-GBM activity of WWOX was generally corresponding, but colony growth and formation were inconsistent in DBTRG-05MG. This prompted us to investigate the molecular landscapes of these cell lines, intending to translate them into the clinical context. METHODS: U87MG/T98G/U251MG/DBTRG-05MG were subjected to high-throughput sequencing, and obtained data were explored via weighted gene co-expression network analysis, differential expression analysis, functional annotation, and network building. Following the identification of the most relevant DBTRG-distinguishing driver genes, data from GBM patients were employed for, e.g., differential expression analysis, survival analysis, and principal component analysis. RESULTS: Although most driver genes were unique for each cell line, some were inversely regulated in DBTRG-05MG. Alongside driver genes, the differentially-expressed genes were used to build a WWOX-related network depicting protein–protein interactions in U87MG/T98G/U251MG/DBTRG-05MG. This network revealed processes distinctly regulated in DBTRG-05MG, e.g., microglia proliferation or neurofibrillary tangle assembly. POLE4 and HSF2BP were selected as DBTRG-discriminating driver genes based on the gene significance, module membership, and fold-change. Alongside WWOX, POLE4 and HSF2BP expression was used to stratify patients into cell lines-resembling groups that differed in, e.g., prognosis and treatment response. Some differences from a WWOX-related network were certified in patients, revealing genes that clarify clinical outcomes. Presumably, WWOX overexpression in DBTRG-05MG resulted in expression profile change resembling that of patients with inferior prognosis and drug response. Among these patients, WWOX may be inaccessible for its partners and does not manifest its anti-cancer activity, which was proposed in the literature but not regarding glioblastoma or concerning POLE4 and HSF2BP. CONCLUSION: Cell lines data enabled the identification of patients among which, despite high expression of WWOX tumor suppressor, no advantageous outcomes were noted due to the cancer-promoting profile ensured by other genes.