Chinese Glioma Genome Atlas Datasets Research Articles

Elevated HSPB1 Expression Is Associated with a Poor Prognosis in Glioblastoma Multiforme Patients.

Glioblastoma multiforme (GBM) is a highly aggressive form of brain cancer. This study investigated the clinical predictive value of heat shock protein β1 (HSPB1) in patients with GBM. A correlation was established between HSPB1 expression and GBM progression using data from The Cancer Genome Atlas (TCGA) dataset, Chinese Glioma Genome Atlas dataset, Gene Expression Omnibus dataset, and Human Protein Atlas database. A survival analysis was conducted and an HSPB1-based nomogram was constructed to evaluate the prognostic value of HSPB1 in patients with GBM. Based on TCGA data mining, we discovered that HSPB1 was significantly elevated in patients with GBM and may reflect their response to immunotherapy. In survival analysis, it appeared to have a predictive role in the prognosis of patients with GBM. Five signaling pathways were significantly enriched in the high HSPB1 expression phenotype according to the gene set enrichment analysis. In addition, a significant association was found between HSPB1 expression and immune checkpoints, tumor immune infiltration, tumor immune microenvironment, and immune cell markers in glioma. Overall, our results suggest that HSPB1 may regulate the function of immune cells, serve as a new immunotherapy target, and predict the response to immunotherapy in patients with GBM. HSPB1 appears to serve as a potential predictor of the clinical prognosis and response to immunotherapy in patients with GBM. It may be possible to identify patients who are likely to benefit from immunotherapy by assessing the expression level of HSPB1.

Construction and validation of a risk model based on the key SNARE proteins to predict the prognosis and immune microenvironment of gliomas.

Synaptic transmission between neurons and glioma cells can promote glioma progression. The soluble N-ethylmaleimide-sensitive fusion factor attachment protein receptors (SNARE) play a key role in synaptic functions. We aimed to construct and validate a novel model based on the SNARE proteins to predict the prognosis and immune microenvironment of glioma. Differential expression analysis and COX regression analysis were used to identify key SRGs in glioma datasets, and we constructed a prognostic risk model based on the key SRGs. The prognostic value and predictive performance of the model were assessed in The Cancer Genome Atlas (TCGA) and Chinese glioma Genome Atlas (CGGA) datasets. Functional enrichment analysis and immune-related evaluation were employed to reveal the association of risk scores with tumor progression and microenvironment. A prognostic nomogram containing the risk score was established and assessed by calibration curves and time-dependent receiver operating characteristic curves. We verified the changes of the key SRGs in glioma specimens and cells by real-time quantitative PCR and Western blot analyses. Vesicle-associated membrane protein 2 (VAMP2) and vesicle-associated membrane protein 5 (VAMP5) were identified as two SRGs affecting the prognoses of glioma patients. High-risk patients characterized by higher VAMP5 and lower VAMP2 expression had a worse prognosis. Higher risk scores were associated with older age, higher tumor grades, IDH wild-type, and 1p19q non-codeletion. The SRGs risk model showed an excellent predictive performance in predicting the prognosis in TCGA and CGGA datasets. Differentially expressed genes between low- and high-risk groups were mainly enriched in the pathways related to immune infiltration, tumor metastasis, and neuronal activity. Immune score, stromal score, estimate score, tumor mutational burden, and expression of checkpoint genes were positively correlated with risk scores. The nomogram containing the risk score showed good performance in predicting the prognosis of glioma. Low VAMP2 and high VAMP5 were found in different grades of glioma specimens and cell lines. We constructed and validated a novel risk model based on the expression of VAMP2 and VAMP5 by bioinformatics analysis and experimental confirmation. This model might be helpful for clinically predicting the prognosis and response to immunotherapy of glioma patients.

Integrative analysis of multi-omics data reveals a pseudouridine-related lncRNA signature for prediction of glioma prognosis and chemoradiotherapy sensitivity

BackgroundGlioblastoma is the most common type of glioma with a high incidence and poor prognosis, and effective medical treatment remains challenging. Pseudouridine (Ψ) is the first post-transcriptional modification discovered and one of the most abundant modifications to RNA. However, the prognostic value of Ψ-related lncRNAs (ΨrLs) for glioma patients has never been systematically evaluated. This study aims to construct a risk model based on ΨrLs signature and to validate the predictive efficiency of the model. MethodTranscriptomic data, genomic data, and relevant clinical data of glioma patients were extracted from the Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). ΨrLs with significant correlation with Ψ-related genes were identified, and univariate Cox regression, least absolute shrinkage and selection operator (LASSO) regression, and multivariate Cox regression were used to further select biomarkers and construct a ΨrLs signature risk model. Then, the expression of lncRNAs of ΨrLs signature in multiple glioma cell lines was detected by qPCR. Further, ROC analysis, stratification analysis, correlation analysis, survival analysis, nomogram, enrichment analysis, immune infiltration analysis, chemoradiotherapy sensitivity analysis, somatic mutation, and recurrent copy number variation (CNV) analysis were used to validate the predictive efficiency of ΨrLs signature in TCGA and CGGA datasets. ResultsA four-lncRNA ΨrLs signature (DNAJC27-AS1, GDNF-AS1, ZBTB20-AS4, and DNMBP-AS1) risk model was constructed. By ROC analysis, stratified analysis, correlation analysis, survival analysis, and nomogram, the signature showed satisfactory predictive efficiency. Functional enrichment analysis revealed the differences in immune-related biological processes between high- and low-risk groups. Immune infiltration analysis showed that the high-risk group had lower tumor purity and higher stromal, immune and ESTIMATE scores. Mitoxantrone was identified as effective drug for low-risk group of glioma patients. Key genes in glioma development, including IDH1, EGFR, PTEN, etc., were differentially mutated between risk groups. The main recurrent CNVs in low-risk groups were 19q13.42 deletion and 7q34 amplification; 10q23.31 deletion and 12q14.1 in the high-risk group. ConclusionsOur study identified a four-lncRNA ΨrLs signature that effectively predicts the prognosis of glioma patients and may serve as a diagnostic tool. Risk scores of glioma patients generated by the signature is associated with immune-related biological processes and chemoradiotherapy sensitivity. These findings may inform the development of more targeted and effective therapies for glioma patients.

The evaluation of six genes combined value in glioma diagnosis and prognosis.

Glioma is the most common and fatal type of brain tumour. Owing to its aggressiveness and lethality, early diagnosis and prediction of patient survival are very important. This study aimed to identify key genes and biomarkers for glioma that can guide clinicians in making rapid diagnosis and prognostication. Data mining of The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), Repository of Molecular Brain Neoplasia Data, and Genotype-Tissue Expression Project brain expression data revealed significantly differentially expressed genes (DEGs), and the risk scores of individual patients were calculated. WGCNA was utilized to screen for genes most related to clinical diagnosis. Prognostic genes associated with glioma were selected via combining the LASSO regression with univariate and multivariate Cox regression and protein-protein interaction network analyses. Then, a nomogram was constructed. And CGGA dataset was utilized to validated. The protein expression levels of the signature were detected using the human protein atlas. Drug response prediction was carried out using the package "pRRophetic". A six-gene signature (KLF6, CHI3L1, SERPINE1, ANGPT2, TGFBR1, and PTX3) was identified and used to stratify patients into low- and high-risk groups. Survival, ROC curve, and Cox analyses clarified that the six hub genes were a favourable independent prognostic factor for patients with glioma. A nomogram was set up by integrating clinical parameters with risk signatures, showing high precision for predicting 2-, 3-, 4-, 5-years survival. In addition, the expression of most genes was consistent with protein expression. Furthermore, the sensitivity to the top ten drugs in the GDSC database of the high-risk group was significantly higher than the low-risk group. Based on genetic profiles and clinicopathological features, including age, grade, isocitrate dehydrogenase mutation status, we constructed a comprehensive prognostic model for patients with glioma. These signatures can be regarded as biomarkers to predict the prognosis of gliomas, possibly providing more therapeutic strategies for future clinical research.

Development and validation of a predictive model in diagnosis and prognosis of primary glioblastoma patients based on Homeobox A family

BackgroundHomeobox A (HOXA) family is involved in the development of malignancies as either tumor suppressors or oncogenes. However, their roles in glioblastoma (GBM) and clinical significance have not been fully elucidated.MethodsHOXA mutation and expressions in pan-cancers were investigated using GSCA and Oncomine, which in GBM were validated by cBioPortal, Chinese Glioma Genome Atlas (CGGA), and The Cancer Genome Atlas (TCGA) datasets. Kaplan–Meier analyses were conducted to determine prognostic values of HOXAs at genetic and mRNA levels. Diagnostic roles of HOXAs in tumor classification were explored by GlioVis and R software. Independent prognostic HOXAs were identified using Cox survival analyses, the least absolute shrinkage and selection operator (LASSO) regression, quantitative real-time PCR, and immunohistochemical staining. A HOXAs-based nomogram survival prediction model was developed and evaluated using Kaplan–Meier analysis, time-dependent Area Under Curve, calibration plots, and Decision Curve Analysis in training and validation cohorts.ResultsHOXAs were highly mutated and overexpressed in pan-cancers, especially in CGGA and TCGA GBM datasets. Genetic alteration and mRNA expression of HOXAs were both found to be prognostic. Specific HOXAs could distinguish IDH mutation (HOXA1-7, HOXA9, HOXA13) and molecular GBM subtypes (HOXA1-2, HOXA9-11, HOXA13). HOXA1/2/3/10 were confirmed to be independent prognostic members, with high expressions validated in clinical GBM tissues. The HOXAs-based nomogram model exhibited good prediction performance and net benefits for patients in training and validation cohorts.ConclusionHOXA family has diagnostic values, and the HOXAs-based nomogram model is effective in survival prediction, providing a novel approach to support the treatment of GBM patients.

Characterization of tumor-associated reactive astrocytes in gliomas by single-cell and bulk tumor sequencing.

Astrocytes constitute approximately 30% of cells in gliomas and play important roles in synapse construction and survival. Recently, JAK/STAT pathway activation associated with a new type of astrocyte was reported. However, the implications of these tumor-associated reactive astrocytes (TARAs) in glioma are not known. We comprehensively assessed TARAs in gliomas, both in single cells and at the bulk tumor level, by analyzing five independent datasets. First, we analyzed two single-cell RNA sequencing datasets of 35,563 cells from 23 patients to estimate the infiltration level of TARAs in gliomas. Second, we collected clinical information and genomic and transcriptomic data of 1,379 diffuse astrocytoma and glioblastoma samples from the Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas datasets to evaluate the genomic, transcriptomic and clinical characteristics of TARA infiltration. Third, we downloaded expression profiles of recurrent glioblastoma samples from patients receiving PD-1 inhibitors to analyze the predictive value of TARAs for immune checkpoint inhibition. Single-cell RNA sequencing data showed TARAs were abundant in the glioma micro-environment (15.7% in the CGGA dataset and 9.1% in the Gene Expression Omnibus GSE141383 dataset, respectively). Bulk tumor sequencing data showed that the extent of TARA infiltration was highly associated with major clinical and molecular features of astrocytic gliomas. Patients with more TARA infiltration were more likely to have MUC16, FLG, and PICK3A mutations, chromosome 9p21.3, 10q23.3, and 13q14.2 deletions and 7p11.2 amplification. Gene Ontology analysis revealed that the high level of astrocyte infiltration was characterized by immune and oncogenic pathways, such as the inflammatory response, positive regulation of the JAK-STAT cascade, positive regulation of NIK/NF-kappa B signaling and the tumor necrosis factor biosynthetic process. Patients with greater TARA infiltration showed inferior prognosis. Meanwhile, the extent of reactive astrocyte infiltration exhibited a predictive value for recurrent glioblastoma patients undergoing anti-PD-1 immune therapy. TARA infiltration might promote glioma tumor progression and can be used as a diagnostic, predictive and prognostic marker in gliomas. Prevention of TARA infiltration might be a new therapeutic strategy for glioma.

Development and validation of an immune infiltration/tumor proliferation-related Notch3 nomogram for predicting survival in patients with primary glioblastoma.

Background: Notch receptors (Notch 1/2/3/4), the critical effectors of the Notch pathway, participate in the tumorigenesis and progression of many malignancies. However, the clinical roles of Notch receptors in primary glioblastoma (GBM) have not been fully elucidated. Methods: The genetic alteration-related prognostic values of Notch receptors were determined in the GBM dataset from The Cancer Genome Atlas (TCGA). Two GBM datasets from TCGA and Chinese Glioma Genome Atlas (CGGA) were used to explore the differential expression between Notch receptors and IDH mutation status, and GBM subtypes. The biological functions of Notch Receptors were explored by Gene Ontology and KEGG analysis. The expression and prognostic significance of Notch receptors were determined in the TCGA and CGGA datasets and further validated in a clinical GBM cohort by immunostaining. A Notch3-based nomogram/predictive risk model was constructed in the TCGA dataset and validated in the CGGA dataset. The model performance was evaluated by receiver operating curves, calibration curves, and decision curve analyses. The Notch3-related phenotypes were analyzed via CancerSEA and TIMER. The proliferative role of Notch3 in GBM was validated in U251/U87 glioma cells by Western blot and immunostaining. Results: Notch receptors with genetic alterations were associated with poor survival of GBM patients. Notch receptors were all upregulated in GBM of TCGA and CGGA databases and closely related to the regulation of transcription, protein-lysine N-methyltransferase activity, lysine N-methyltransferase activity, and focal adhesion. Notch receptors were associated with Classical, Mesenchymal, and Proneural subtypes. Notch1 and Notch3 were closely correlated with IDH mutation status and G-CIMP subtype. Notch receptors displayed the differential expression at the protein level and Notch3 showed a prognostic significance in a clinical GBM cohort. Notch3 presented an independent prognostic role for primary GBM (IDH1 mutant/wildtype). A Notch3-based predictive risk model presented favorable accuracy, reliability, and net benefits for predicting the survival of GBM patients (IDH1 mutant/wildtype and IDH1 wildtype). Notch3 was closely related to immune infiltration (macrophages, CD4+ T cells, and dendritic cells) and tumor proliferation. Conclusion: Notch3-based nomogram served as a practical tool for anticipating the survival of GBM patients, which was related to immune-cell infiltration and tumor proliferation.

M2 macrophage marker CHI3L2 could serve as a potential prognostic and immunological biomarker in glioma by integrated single-cell and bulk RNA-Seq analysis.

CHI3L2 plays a crucial role in multiple cancers, but its importance in glioma remains unclear. Hence, we comprehensively integrated bulk RNA-sequencing (RNA-seq), proteomics and single-cell RNA-seq (scRNA-seq) to determine the roles of CHI3L2 in gliomas. Bulk RNA-seq, proteomics and scRNA-seq data of CHI3L2 in glioma were obtained from online databases. The quantitative real-time polymerase chain reaction (qRT-PCR) and immunohistochemistry (IHC) were conducted to verify the CHI3L2 expression. Then, univariate and multivariate Cox regression analyses, Norman charts and gene set enrichment analysis (GSEA) were performed. Finally, the associations between CHI3L2 and tumor immunity were explored. The expression of CHI3L2 was markedly higher in glioma cancers compared with normal tissues from analysis of the data of the Cancer Genome Atlas and Chinese Glioma Genome Atlas datasets and as verified by GSE4290, GSE50161, qRT-PCR and IHC results (p < 0.05). High expression of CHI3L2 suggested poor overall survival (OS) prognosis in gliomas (p < 0.05). CHI3L2 might also serve as an independent predictor of OS for gliomas (p < 0.05) and we also constructed a Norman chart to predict these patients' survival prognosis with good performance. GSEA analysis showed that CHI3L2 might be involved with eight pathways in gliomas. Regarding tumor immunity, CHI3L2 was found to be significantly involved with immune cell infiltration levels of low-grade glioma, the tumor immune microenvironment, immune checkpoints and immune cells in both low-grade glioma and glioblastoma (p < 0.05). Additionally, scRNA-seq data for CHI3L2 in glioma from the TISCH2 website showed that CHI3L2 is mainly expressed in astrocytes, endothelial cells, CD8+ T cells, mono/macrophage cells, etc. CONCLUSIONS: CHI3L2 presents prognostic and immunological values in glioma, providing novel therapeutic targets for glioma patients.

Dynamic Regulation Genes at Microtubule Plus Ends: A Novel Class of Glioma Biomarkers

Glioma is the most prevalent and aggressive primary nervous system tumor with an unfavorable prognosis. Microtubule plus-end-related genes (MPERGs) play critical biological roles in the cell cycle, cell movement, ciliogenesis, and neuronal development by coordinating microtubule assembly and dynamics. This research seeks to systematically explore the oncological characteristics of these genes in microtubule-enriched glioma, focusing on developing a novel MPERG-based prognostic signature to improve the prognosis and provide more treatment options for glioma patients. First, we thoroughly analyzed and identified 45 differentially expressed MPERGs in glioma. Based on these genes, glioma patients were well distinguished into two subgroups with survival and tumor microenvironment infiltration differences. Next, we further screened the independent prognostic genes (CTTNBP2, KIF18A, NAV1, SLAIN2, SRCIN1, TRIO, and TTBK2) using 36 prognostic-related differentially expressed MPERGs to construct a signature with risk stratification and prognostic prediction ability. An increased risk score was related to the malignant progression of glioma. Therefore, we also designed a nomogram model containing clinical factors to facilitate the clinical use of the risk signature. The prediction accuracy of the signature and nomogram model was verified using The Cancer Genome Atlas and Chinese Glioma Genome Atlas datasets. Finally, we examined the connection between the signature and tumor microenvironment. The signature positively correlated with tumor microenvironment infiltration, especially immunoinhibitors and the tumor mutation load, and negatively correlated with microsatellite instability and cancer stemness. More importantly, immune checkpoint blockade treatment and drug sensitivity analyses confirmed that this prognostic signature was helpful in anticipating the effect of immunotherapy and chemotherapy. In conclusion, this research is the first study to define and validate an MPERG-based signature closely associated with the tumor microenvironment as a reliable and independent prognostic biomarker to guide personalized choices of immunotherapy and chemotherapy for glioma patients.

Differences of macrophages in the tumor microenvironment as an underlying key factor in glioma patients.

Macrophages, the major immune cells in glioma microenvironment, are closely related to tumor prognosis. Further studies are needed to investigate macrophages, which will be helpful to fully understand the role of it and early achieve clinical translation. A total of 1334 glioma cases were enrolled in this study from 3 databases. In our works, the single cell cohorts from GSE89567, GSE84465, and the Chinese Glioma Genome Atlas (CGGA) datasets were used to analyze the key genes of macrophage. The bulk sequencing data from the Cancer Genome Atlas (TCGA) and CGGA datasets were respectively divided into the training set and validation set to test prognostic value of the key genes from single cell analysis. Quantitative and functional differences significantly emerge in macrophage clusters between LGG and GBM. Firstly, we used the Seurat R package to identify 281 genes differentially expressed genes in macrophage clusters between LGG and GBM. Furthermore, based on these genes, we developed a predictive risk model to predict prognosis and reflect the immune microenvironment in glioma. The risk score calculation formula was yielded as follows: Risk score = (0.11 × EXPMACC1) + (-0.31 × EXPOTUD1) + (-0.09 × EXPTCHH) + (0.26 × EXPADPRH) + (-0.40× EXPABCG2) + (0.21 × EXPPLBD1) + (0.12 × EXPANG) + (0.29 × EXPQPCT). The risk score was independently related to prognosis. Further, significant differences existed in immunological characteristics between the low- and high-risk score groups. What is more, mutation analysis found different genomic patterns associated with the risk score. This study further confirms that the proportion of macrophage infiltration is not only significantly different, but the function of them is also different. The signature, identified from the differentially expressed macrophage-related genes impacts poor prognosis and short overall survival and may act as therapeutic targets in the future.

Identification of tumor-associated antigens and immune subtypes of lower-grade glioma and glioblastoma for mRNA vaccine development

BackgroundmRNA became a promising therapeutic approach in many diseases. This study aimed to identify the tumor antigens specifically expressed in tumor cells for lower-grade glioma (LGG) and glioblastoma (GBM) patients.MethodsIn this work, the mRNA microarray expression profile and clinical data were obtained from 301 samples in the Chinese Glioma Genome Atlas (CGGA) database, the mRNA sequencing data and clinical data of 701 samples were downloaded from The Cancer Genome Atlas (TCGA) database. Genetic alterations profiles were extracted from CGGA and cBioPortal datasets. R language and GraphPad Prism software were applied for the statistical analysis and graph work.ResultsPTBP1 and SLC39A1, which were overexpressed and indicated poor prognosis in LGG patients, were selected as tumor-specific antigens for LGG patients. Meanwhile, MMP9 and SLC16A3, the negative prognostic factors overexpressed in GBM, were identified as tumor-specific antigens for GBM patients. Besides, three immune subtypes (LGG1-LGG3) and eight WGCNA modules were identified in LGG patients. Meanwhile, two immune subtypes (GBM1–GBM2) and 10 WGCNA modules were selected in GBM. The immune characteristics and potential functions between different subtypes were diversity. LGG2 and GBM1 immune subtype were associated with longer overall survival than other subtypes.ConclusionIn this study, PTBP1 and SLC39A1 are promising antigens for mRNA vaccines development in LGG, and MMP9 and SLC16A3 were potential antigens in GBM. Our analyses indicated that mRNA vaccine immunotherapy was more suitable for LGG2 and GBM1 subtypes. This study was helpful for the development of glioma immunotherapies.

Integrin subunit alpha V is a potent prognostic biomarker associated with immune infiltration in lower-grade glioma

As a member of integrin receptor family, ITGAV (integrin subunit α V) is involved in a variety of cell biological processes and overexpressed in various cancers, which may be a potential prognostic factor. However, its prognostic value and potential function in lower-grade glioma (LGG) are still unclear, and in terms of immune infiltration, it has not been fully elucidated. Here, the expression preference, prognostic value, and clinical traits of ITGAV were investigated using The Cancer Genome Atlas database (n = 528) and the Chinese Glioma Genome Atlas dataset (n = 458). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses and gene set enrichment analysis (GSEA) were used to explore the biological function of ITGAV. Using R package "ssGSEA" analysis, it was found thatthe ITGAV mRNA expression level showed intense correlation with tumor immunity, such as tumor-infiltrating immune cells and multiple immune-related genes. In addition, ITGAV is associated with some immune checkpoints and immune checkpoint blockade (ICB) and response to chemotherapy. and the expression of ITGAV protein in LGG patients was verified via immunohistochemistry (IHC). ITGAV expression was higher in LGG tissues than in normal tissues (P < 0.001) and multifactor analysis showed that ITGAV mRNA expression was an independent prognostic factor for LGG overall survival (OS; hazard ratio = 2.113, 95% confidence interval = 1.393-3.204, P < 0.001). GSEA showed that ITGAV expression was correlated with Inflammatory response, complement response, KRAS signal, and interferon response. ssGSEA results showed a positive correlation between ITGAV expression and Th2 cell infiltration level. ITGAV mRNA was overexpressed in LGG, and high ITGAV mRNA levels were found to be associated with poor protein expression and poor OS. ITGAV is therefore a potential biomarker for the diagnosis and prognosis of LGG and may be a potential immunotherapy target.

MELK is a prognostic biomarker and correlated with immune infiltration in glioma

Glioma accounts for the vast majority of primary brain tumors with inevitable recurrence and poor prognosis. Maternal embryonic leucine zipper kinase (MELK) is overexpressed in multiple human tumors and could activate a variety of oncogenic-associated signal pathways. However, its role in the glioma microenvironment is still largely unknown. We collected the RNA sequence data and clinical information of gliomas from the Chinese Glioma Genome Atlas (CGGA), The Cancer Genome Atlas (TCGA), and the Gene Expression Omnibus (GEO) databases, and investigated MELK expression and its correlation with clinicopathologic features and prognosis in glioma. Moreover, the relationship between MELK expression and immune cell infiltration in the tumor microenvironment of gliomas was explored through single-sample gene set enrichment analysis (ssGSEA) and CIBERSORT. In addition, gene set enrichment analysis (GSEA) and Metascape online analysis were performed to find out signaling pathways enriched by differentially expressed genes (DEGs) between high- and low-MELK expression groups. Finally, immunohistochemistry was performed to validate our findings. Data analysis of CGGA and GEO datasets showed that MELK was significantly upregulated in gliomas than in normal brain tissues, and MELK expression was obviously correlated with clinicopathologic features, including age, WHO grade, histological subtype, IDH mutant status, 1p19q codeletion status, and PRS type. Stratified analysis, Cox regression analysis, and nomogram model revealed that high expression of MELK predicted poor survival; hence, MELK could serve as an independent prognostic biomarker for glioma. Moreover, results from enrichment pathway analysis indicated that the immune system process, angiogenesis, apoptosis, cell cycle, and other oncogenic-related signal pathways were significantly enriched between high- and low-MELK expression groups. Immune infiltration analysis demonstrated that increased MELK expression was significantly correlated with higher immune scores, higher fractions of immunocytes (T cells, NK cells resting, macrophages, resting mast cells, and neutrophils), and higher expression levels of immune checkpoints (B7-H3, CTLA4, LAG3, PD-1, PD-L1, and TIM3). Finally, immunohistochemistry analysis validated our findings that high expression of MELK relates to increased malignancy and poor prognosis of glioma. Our findings identified that MELK could act as an independent prognostic indicator and potential immunotherapy target for glioma. In conclusion, these findings suggested that DDOST mediated the immunosuppressive microenvironment of gliomas and could be an important biomarker in diagnosing and treating gliomas.

Establishment and validation of a novel prognostic model for lower-grade glioma based on senescence-related genes.

Increasing studies have indicated that senescence was associated with tumorigenesis and progression. Lower-grade glioma (LGG) presented a less invasive nature, however, its treatment efficacy and prognosis prediction remained challenging due to the intrinsic heterogeneity. Therefore, we established a senescence-related signature and investigated its prognostic role in LGGs. The gene expression data and clinicopathologic features were from The Cancer Genome Atlas (TCGA) database. The experimentally validated senescence genes (SnGs) from the CellAge database were obtained. Then LASSO regression has been performed to build a prognostic model. Cox regression and Kaplan-Meier survival curves were performed to investigate the prognostic value of the SnG-risk score. A nomogram model has been constructed for outcome prediction. Immunological analyses were further performed. Data from the Chinese Glioma Genome Atlas (CGGA), Repository of Molecular Brain Neoplasia Data (REMBRANDT), and GSE16011 were used for validation. The 6-SnG signature has been established. The results showed SnG-risk score could be considered as an independent predictor for LGG patients (HR=2.763, 95%CI=1.660-4.599, P&lt;0.001). The high SnG-risk score indicated a worse outcome in LGG (P&lt;0.001). Immune analysis showed a positive correlation between the SnG-risk score and immune infiltration level, and the expression of immune checkpoints. The CGGA datasets confirmed the prognostic role of the SnG-risk score. And Kaplan-Meier analyses in theadditional datasets (CGGA, REMBRANDT, and GSE16011) validated the prognostic role of the SnG-signature (P&lt;0.001 for all). The SnG-related prognostic model could predict the survival of LGG accurately. This study proposed a novel indicator for predicting the prognosis of LGG and provided potential therapeutic targets.

Deep learning algorithm reveals two prognostic subtypes in patients with gliomas

BackgroundGliomas are highly complex and heterogeneous tumors, rendering prognosis prediction challenging. The advent of deep learning algorithms and the accessibility of multi-omic data represent a new approach for the identification of survival-sensitive subtypes. Herein, an autoencoder-based approach was used to identify two survival-sensitive subtypes using RNA sequencing (RNA-seq) and DNA methylation (DNAm) data from The Cancer Genome Atlas (TCGA) dataset. The subtypes were used as labels to build a support vector machine model with cross-validation. We validated the robustness of the model on Chinese Glioma Genome Atlas (CGGA) dataset. DNAm-driven genes were identified by integrating DNAm and gene expression profiling analyses using the R MethylMix package and carried out for further enrichment analysis.ResultsFor TCGA dataset, the model produced a high C-index (0.92 ± 0.02), low brier score (0.16 ± 0.02), and significant log-rank p value (p < 0.0001). The model also had a decent performance for CGGA dataset (CGGA DNAm: C-index of 0.70, brier score of 0.21; CGGA RNA-seq: C-index of 0.79, brier score of 0.18). Moreover, we identified 389 DNAm-driven genes of survival-sensitive subtypes, which were significantly enriched in the glutathione metabolism pathway.ConclusionsOur study identified two survival-sensitive subtypes of glioma and provided insights into the molecular mechanisms underlying glioma development; thus, potentially providing a new target for the prognostic prediction of gliomas and supporting personalized treatment strategies.

Large-scale bulk RNA-seq analysis defines immune evasion mechanism related to mast cell in gliomas.

Accumulating evidence has demonstrated that the immune cells have an emerging role in controlling anti-tumor immune responses and tumor progression. The comprehensive role of mast cell in glioma has not been illustrated yet. In this study, 1,991 diffuse glioma samples were collected from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). xCell algorithm was employed to define the mast cell-related genes. Based on mast cell-related genes, gliomas were divided into two clusters with distinct clinical and immunological characteristics. The survival probability of cluster 1 was significantly lower than that of cluster 2 in the TCGA dataset, three CGGA datasets, and the Xiangya cohort. Meanwhile, the hypoxic and metabolic pathways were active in cluster 1, which were beneficial to the proliferation of tumor cells. A potent prognostic model based on mast cell was constructed. Via machine learning, DRG2 was screened out as a characteristic gene, which was demonstrated to predict treatment response and predict survival outcome in the Xiangya cohort. In conclusion, mast cells could be used as a potential effective prognostic factor for gliomas.

CDH6 as a prognostic indicator and marker for chemotherapy in gliomas.

Glioma is the most malignant cancer of the central nervous system. There are various therapies for treating gliomas, but their outcomes are not satisfactory. Therefore, new targets for glioma treatment are needed. This study examined the cadherin-6 (CDH6) expression in gliomas using The Cancer Genome Atlas and Chinese Glioma Genome Atlas datasets. CDH6 expression positively correlated with the World Health Organization (WHO) tumor grade and negatively correlated with patient prognosis. A significant decrease in CDH6 promoter methylation was identified with an increase in the WHO grade severity. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses suggested that CDH6 might be involved in cell–cell interactions and immune processes in the glioma microenvironment. Weighted gene co-expression network analysis revealed a correlation between CDH6 and cell adhesion molecules, focal adhesions, phosphatidylinositol 3-kinase-protein kinase B signaling pathways, nuclear division, chromosome segregation, mitotic nuclear division, and immune-related pathways. CDH6 strongly correlated with immunosuppressive cells, including regulatory T cells, monocytes, macrophages, tumor-associated macrophages, and myeloid-derived suppressor cells. It also showed correlations with immune-active cells such as B cells, CD8+ T cells, and dendritic cells. Single-cell analysis showed that CDH6 was expressed mainly in astrocyte (AC)-like malignant cells. Differentially expressed genes of AC-like malignant cells were found to be associated with stress response, membranous processes, viral infections, and several types of cancers. Potential drugs associated with high CDH6 expression were also predicted, including AMG-22, rutin, CCT128930, deforolimus, bis(maltolato)oxovanadium, anagrelide, vemurafenib, CHIR-98014, and AZD5582. Thus, this study showed that CDH6 correlates with glioma immune infiltration, it is expressed mainly in AC-like malignant cells, and it may act as a new target for glioma therapy.

Brain radiotoxicity-related 15CAcBRT gene expression signature predicts survival prognosis of glioblastoma patients.

Glioblastoma is the most common and devastating primary brain cancer. Radiotherapy is standard of care; however, it is associated with brain radiation toxicity (BRT). This study used a multi-omics approach to determine whether BRT-related genes (RGs) harbor survival prognostic value and whether their encoded proteins represent novel therapeutic targets for glioblastoma. RGs were identified through analysis of single-nucleotide variants associated with BRT (R-SNVs). Functional relationships between RGs were established using Protein-Protein Interaction networks. The influence of RGs and their functional groups on glioblastoma prognosis was evaluated using clinical samples from the Glioblastoma Bio-Discovery Portal database and validated using the Chinese Glioma Genome Atlas dataset. The identification of clusters of radiotoxic and putative pathogenic variants in proteins encoded by RGs was achieved by computational 3D structural analysis. We identified the BRT-related 15CAcBRT molecular signature with prognostic value in glioblastoma, by analysis of the COMT and APOE protein functional groups. Its external validation confirmed clinical relevance independent of age, MGMT promoter methylation status, and IDH mutation status. Interestingly, the genes IL6, APOE, and MAOB documented significant gene expression levels alteration, useful for drug repositioning. Biological networks associated with 15CAcBRT signature involved pathways relevant to cancer and neurodegenerative diseases. Analysis of 3D clusters of radiotoxic and putative pathogenic variants in proteins coded by RGs unveiled potential novel therapeutic targets in neuro-oncology. 15CAcBRT is a BRT-related molecular signature with prognostic significance for glioblastoma patients and represents a hub for drug repositioning and development of novel therapies.

CCR5 as a Prognostic Factor in Lower-Grade Glioma is Involved in the Remodeling of the Tumor Microenvironment.

Background: Lower-grade gliomas (LGGs) carry a high risk of malignant transformation, leading to severe neurologic deterioration and ultimately, death. The tumor microenvironment (TME) plays an essential role in tumor maintenance, progression, and immunotherapy resistance. Therefore, the LGG TME deserves comprehensive exploration for a novel therapeutic target. Methods: The ESTIMATE algorithm was used to estimate infiltrating stromal and immune cells of LGG patients obtained from the Cancer Genome Atlas (TCGA) database. Kaplan–Meier analysis was performed to classify survival differences. TME-related differentially expressed genes were identified between the low- and high-immune/stromal groups. Hub genes were screened by constructing protein–protein interaction networks and performing the Cox regression analysis. Differential analysis, survival analysis, gene set enrichment analysis, and clinical relevance analysis specific to hub genes were evaluated by using the TCGA and the Chinese Glioma Genome Atlas datasets, and the results were validated by qRT-PCR, Western blotting, and immunohistochemistry in tissues from LGG patients. Results: The immune and stromal components in TME were negatively related to patient prognosis. Differentially expressed genes sharing immune score and stromal score were mainly involved in the immune response. C-C chemokine receptor type 5 (CCR5), as only a hub gene, was significantly higher in LGG patients than normal patients and negatively correlated with the prognosis of patients. High-expression CCR5 was positively related to immune-related and tumor progression pathways. CCR5 protein expression was higher in LGG with isocitrate dehydrogenase wildtype. Validated results showed that CCR5 was upregulated in LGG tissues at mRNA and protein levels and could affect immune cell infiltration. These results suggested that CCR5 was a potential indicator for the status of TME. Conclusion: Glioma cells remodel the immune microenvironment through the high expression of CCR5 and lead to a poor prognosis in patients with LGG. The inhibition of CCR5 may contribute to the efficacy of LGG immunotherapy.

Prognostic Value and Biological Function of Galectins in Malignant Glioma.

Malignant glioma is the most common solid tumor of the adult brain, with high lethality and poor prognosis. Hence, identifying novel and reliable biomarkers can be advantageous for diagnosing and treating glioma. Several galectins encoded by LGALS genes have recently been reported to participate in the development and progression of various tumors; however, their detailed role in glioma progression remains unclear. Herein, we analyzed the expression and survival curves of all LGALS across 2,217 patients with glioma using The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), and Rembrandt databases. By performing multivariate Cox analysis, we built a survival model containing LGALS1, LGALS3, LGALS3BP, LGALS8, and LGALS9 using TCGA database. The prognostic power of this panel was assessed using CGGA and Rembrandt datasets. ESTIMATE and CIBERSORT algorithms confirmed that patients in high-risk groups exhibited significant stromal and immune cell infiltration, immunosuppression, mesenchymal subtype, and isocitrate dehydrogenase 1 (IDH1) wild type. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), CancerSEA, and Gene Set Enrichment Analysis (GSEA) showed that pathways related to hypoxia, epithelial-to-mesenchymal transition (EMT), stemness, and inflammation were enriched in the high-risk group. To further elucidate the function of LGALS in glioma, we performed immunohistochemical staining of tissue microarrays (TMAs), Western blotting, and cell viability, sphere formation, and limiting dilution assays following lentiviral short hairpin RNA (shRNA)-mediated LGALS knockdown. We observed that LGALS expression was upregulated in gliomas at both protein and mRNA levels. LGALS could promote the stemness maintenance of glioma stem cells (GSCs) and positively correlate with M2-tumor-associated macrophages (TAMs) infiltration. In conclusion, we established a reliable survival model for patients with glioma based on LGALS expression and revealed the essential roles of LGALS genes in tumor growth, immunosuppression, stemness maintenance, pro-neural to mesenchymal transition, and hypoxia in glioma.