BackgroundGlioblastoma (GBM) is the most prominent and aggressive primary brain tumor in adults. Anoikis is a specific form of programmed cell death that plays a key role in tumor invasion and metastasis. The presence of anti-anoikis factors is associated with tumor aggressiveness and drug resistance.MethodsThe non-negative matrix factorization algorithm was used for effective dimension reduction for integrated datasets. Differences in the tumor microenvironment (TME), stemness indices, and clinical characteristics between the two clusters were analyzed. Difference analysis, weighted gene coexpression network analysis (WGCNA), univariate Cox regression, and least absolute shrinkage and selection operator regression were leveraged to screen prognosis-related genes and construct a risk score model. Immunohistochemistry was performed to evaluate the expression of representative genes in clinical specimens. The relationship between the risk score and the TME, stemness, clinical traits, and immunotherapy response was assessed in GBM and pancancer.ResultsTwo definite clusters were identified on the basis of anoikis-related gene expression. Patients with GBM assigned to C1 were characterized by shortened overall survival, higher suppressive immune infiltration levels, and lower stemness indices. We further constructed a risk scoring model to quantify the regulatory patterns of anoikis-related genes. The higher risk score group was characterized by a poor prognosis, the infiltration of suppressive immune cells and a differentiated phenotype, whereas the lower risk score group exhibited the opposite effects. In addition, patients in the lower risk score group exhibited a higher frequency of isocitrate dehydrogenase (IDH) mutations and a more sensitive response to immunotherapy. Drug sensitivity analysis was performed, revealing that the higher risk group may benefit more from drugs targeting the PI3K/mTOR signaling pathway.ConclusionWe revealed potential relationships between anoikis-related genes and clinical features, TME, stemness, IDH mutation, and immunotherapy and elucidated their therapeutic value.

Background Glioblastoma (GBM) is the most malignant, aggressive and recurrent primary brain tumor. Cell senescence can cause irreversible cessation of cell division in normally proliferating cells. According to studies, senescence is a primary anti-tumor mechanism that may be seen in a variety of tumor types. It halts the growth and spread of tumors. Tumor suppressive functions held by cellular senescence provide new directions and pathways to promote cancer therapy. Methods We comprehensively analyzed the cell senescence-associated genes expression patterns. The potential molecular subtypes were acquired based on unsupervised cluster analysis. The tumor immune microenvironment (TME) variations, immune cell infiltration, and stemness index between 3 subtypes were analyzed. To identify genes linked with GBM prognosis and build a risk score model, we used weighted gene co-expression network analysis (WGCNA), univariate Cox regression, Least absolute shrinkage and selection operator regression (LASSO), and multivariate Cox regression analysis. And the correlation between risk scores and clinical traits, TME, GBM subtypes, as well as immunotherapy responses were estimated. Immunohistochemistry (IHC) and cellular experiments were performed to evaluate the expression and function of representative genes. Then the 2 risk scoring models were constructed based on the same method of calculation whose samples were acquired from the CGGA dataset and TCGA datasets to verify the rationality and the reliability of the risk scoring model. Finally, we conducted a pan-cancer analysis of the risk score, assessed drug sensitivity based on risk scores, and analyzed the pathways of sensitive drug action. Results The 3 potential molecular subtypes were acquired based on cell senescence-associated genes expression. The Log-rank test showed the difference in GBM patient survival between 3 potential molecular subtypes (P = 0.0027). Then, 11 cell senescence-associated genes were obtained to construct a risk-scoring model, which was systematically randomized to distinguish the train set (n = 293) and the test set (n = 292). The Kaplan-Meier (K-M) analyses indicated that the high-risk score in the train set (P < 0.0001), as well as the test set (P = 0.0053), corresponded with poorer survival. In addition, the high-risk score group showed a poor response to immunotherapy. The reliability and credibility of the risk scoring model were confirmed according to the CGGA dataset, TCGA datasets, and Pan-cancer analysis. According to drug sensitivity analysis, it was discovered that LJI308, a potent selective inhibitor of RSK pathways, has the highest drug sensitivity. Moreover, the GBM patients with higher risk scores may potentially be more beneficial from drugs that target cell cycle, mitosis, microtubule, DNA replication and apoptosis regulation signaling. Conclusion We identified potential associations between clinical characteristics, TME, stemness, subtypes, and immunotherapy, and we clarified the therapeutic usefulness of cell senescence-associated genes.