RESUMEN
Analysis of tumor gene expression is an important approach for the classification and identification of therapeutic vulnerabilities. However, targeting glioblastoma (GBM) based on molecular subtyping has not yet translated into successful therapies. Here, we present an integrative approach based on molecular pathways to expose new potentially actionable targets. We used gene set enrichment analysis (GSEA) to conduct an unsupervised clustering analysis to condense the gene expression data from bulk patient samples and patient-derived gliomasphere lines into new gene signatures. We identified key targets that are predicted to be differentially activated between tumors and were functionally validated in a library of gliomasphere cultures. Resultant cluster-specific gene signatures associated not only with hallmarks of cell cycle and stemness gene expression, but also with cell-type specific markers and different cellular states of GBM. Several upstream regulators, such as PIK3R1 and EBF1 were differentially enriched in cells bearing stem cell like signatures and bear further investigation. We identified the transcription factor E2F1 as a key regulator of tumor cell proliferation and self-renewal in only a subset of gliomasphere cultures predicted to be E2F1 signaling dependent. Our in vivo work also validated the functional significance of E2F1 in tumor formation capacity in the predicted samples. E2F1 inhibition also differentially sensitized E2F1-dependent gliomasphere cultures to radiation treatment. Our findings indicate that this novel approach exploring cancer pathways highlights key therapeutic vulnerabilities for targeting GBM.