Tumor mutation burden (TMB) is a useful biomarker to predict prognosis and the efficacy of immune checkpoint inhibitors (ICIs). In this study, we aimed to explore the prognostic value of TMB and the potential association between TMB and immune infiltration in lower-grade gliomas (LGGs). Somatic mutation and RNA-sequencing (RNA-seq) data were downloaded from the Cancer Genome Atlas (TCGA) database. TMB was calculated and patients were divided into high- and low-TMB groups. After performing differential analysis between high- and low-risk groups, we identified six hub TMB and immune-related genes that were correlated with overall survival in LGGs. Then, Gene Set Enrichment Analysis was performed to screen significantly enriched GO terms between the two groups. Moreover, an immune-related risk score system was developed by LASSO Cox analysis based on the six hub genes and was validated with the Chinese Glioma Genome Atlas dataset. Using the TIMER database, we further systematically analyzed the relationships between mutants of the six hub genes and immune infiltration levels, as well as the relationships between the immune-related risk score system and the immune microenvironment in LGGs. The results showed that TMB was negatively correlated with OS and high TMB might inhibit immune infiltration in LGGs. Furthermore, the risk score system could effectively stratify patients into low- and high-risk groups in both the training and validation datasets. Multivariate Cox analysis demonstrated that TMB was not an independent prognostic factor, but the risk score was. Higher infiltration of immune cells (B cells, CD4 + T cells, CD8 + T cells, neutrophils, macrophages, and dendritic cells) and higher levels of immune checkpoints (PD-1, CTLA-4, LAG-3, and TIM-3) were found in patients in the high-risk group. Finally, a novel nomogram model was constructed and evaluated to estimate the overall survival of LGG patients. In summary, our study provided new insights into immune infiltration in the tumor microenvironment and immunotherapies for LGGs.

Rationale: Diffuse glioma patients have high mortality and recurrence despite multimodal therapies. This study aims to identify the potential tumor antigens for mRNA vaccines and subtypes suitable for the immunotherapy of patients with diffuse glioma. Methods: Gene expression profiles and corresponding clinical information were obtained from the Chinese Glioma Genome Atlas (CGGA) and the Cancer Genome Atlas (TCGA) databases. Genetic alterations were extracted from cBioPortal. Differential gene analysis, survival analysis, correlation analysis, consensus clustering analysis, and immune cell infiltration analysis were conducted based on the various databases. Finally, the hub genes, the modules related to tumor antigens, and the immune subtypes were identified using WGCNA method. Results: Three over-expressed, amplified, and mutated tumor antigens, including KDR, COL1A2, and SAMD9, were associated with clinical outcomes. The expression of the three genes had a positive correlation with the abundance of antigen-presenting cells (APCs) and APC marker expression. Subsequently, three immune subtypes (Ims1, Ims2, and Ims3) were distinguished in the TCGA cohort, which exhibited distinct molecular, cellular, and clinical characteristics consistent with the CGGA cohort. Diffuse gliomas with subtype Ims1 were more malignant with immunosuppressive phenotypes and more associated with poor prognosis than the other two subtypes. The three antigens and the immune checkpoints were differentially expressed among the three immune subtypes. Finally, functional enrichment analysis of the genes related to tumor antigens and immune subtypes suggested that they are enriched in many immune-associated processes. Conclusions: KDR, COL1A2, and SAMD9 are potential antigens for developing mRNA vaccines against diffuse glioma. The results suggest that immunotherapy targeting these three antigens is more suitable for patients with subtype Ims1. This study provides insights into immunotherapy for diffuse glioma.

Methyltransferase-like 7B (METTL7B) is a member of the methyltransferase-like protein family that plays an important role in the development and progression of tumors. However, its prognostic value and the correlation of METTL7B expression and tumor immunity in some cancers remain unclear. By analyzing online data, we found that METTL7B is abnormally overexpressed in multiple human tumors and plays an important role in the overall survival (OS) of patients with 8 cancer types and disease-free survival (DFS) of patients with 5 cancer types. Remarkably, METTL7B expression was positively correlated with the OS and DFS of patients with lower-grade glioma (LGG). In addition, a positive correlation between METTL7B expression and immune cell infiltration in LGG was observed. Moreover, we identified a strong correlation between METTL7B expression and immune checkpoint gene expression in kidney chromophobe (KICH), LGG and pheochromocytoma and paraganglioma (PCPG). Furthermore, METTL7B was involved in the extracellular matrix (ECM) and immune-related pathways in LGGs. Finally, in vitro experiments showed that knockdown of METTL7B inhibited the growth, migration, invasion and the epithelial–mesenchymal transition (EMT) of LGG cells. METTL7B expression potentially represents a novel prognostic biomarker due to its significant association with immune cell infiltration in LGG.

IntroductionGlioma is the most common primary cancer of the central nervous system with dismal prognosis. Long noncoding RNAs (lncRNAs) have been discovered to play key roles in tumorigenesis in various cancers, including glioma. Because of the relevance between immune infiltrating and clinical outcome of glioma, identifying immune-related lncRNAs is urgent for better personalized management.Materials and methodsSingle-sample gene set enrichment analysis (ssGSEA) was applied to estimate immune infiltration, and glioma samples were divided into high immune cell infiltration group and low immune cell infiltration group. After screening differentially expressed lncRNAs in two immune groups, least absolute shrinkage and selection operator (LASSO) Cox regression analysis was performed to construct an immune-related prognostic signature. Additionally, we explored the correlation between immune infiltration and the prognostic signature.ResultsA total of 653 samples were appropriate for further analyses, and 10 lncRNAs were identified as immune-related lncRNAs in glioma. After univariate Cox regression and LASSO Cox regression analysis, six lncRNAs were identified to construct a prognostic signature for glioma, which could be taken as independent prognostic factors in both univariate and multivariate Cox regression analyses. Moreover, risk score was significantly correlated with all the 29 immune-related checkpoint expression (p < 0.05) in ssGSEA except neutrophils (p = 0.43).ConclusionThe study constructed an immune-related prognostic signature for glioma, which contributed to improve clinical outcome prediction and guide immunotherapy.

Glioma is the common histological subtype of malignancy in the central nervous system, with high morbidity and mortality. Glioma cancer stem cells (CSCs) play essential roles in tumor recurrence and treatment resistance. Thus, exploring the stem cell-related genes and subtypes in glioma is important. In this study, we collected the RNA-sequencing (RNA-seq) data and clinical information of glioma patients from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) databases. With the differentially expressed genes (DEGs) and weighted gene correlation network analysis (WGCNA), we identified 86 mRNA expression-based stemness index (mRNAsi)-related genes in 583 samples from TCGA RNA-seq dataset. Furthermore, these samples from TCGA database could be divided into two significantly different subtypes with different prognoses based on the mRNAsi corresponding gene, which could also be validated in the CGGA database. The clinical characteristics and immune cell infiltrate distribution of the two stemness subtypes are different. Then, functional enrichment analyses were performed to identify the different gene ontology (GO) terms and pathways in the two different subtypes. Moreover, we constructed a stemness subtype-related risk score model and nomogram to predict the prognosis of glioma patients. Finally, we selected one gene (ETV2) from the risk score model for experimental validation. The results showed that ETV2 can contribute to the invasion, migration, and epithelial-mesenchymal transition (EMT) process of glioma. In conclusion, we identified two distinct molecular subtypes and potential therapeutic targets of glioma, which could provide new insights for the development of precision diagnosis and prognostic prediction for glioma patients.