Gliomas are the most common and malignant intracranial tumors in adults. Recent studies have revealed the significance of functional genomics for glioma pathophysiological studies and treatments. However, access to comprehensive genomic data and analytical platforms is often limited. Here, we developed the Chinese Glioma Genome Atlas (CGGA), a user-friendly data portal for the storage and interactive exploration of cross-omics data, including nearly 2000 primary and recurrent glioma samples from Chinese cohort. Currently, open access is provided to whole-exome sequencing data (286 samples), mRNA sequencing (1018 samples) and microarray data (301 samples), DNA methylation microarray data (159 samples), and microRNA microarray data (198 samples), and to detailed clinical information (age, gender, chemoradiotherapy status, WHO grade, histological type, critical molecular pathological information, and survival data). In addition, we have developed several tools for users to analyze the mutation profiles, mRNA/microRNA expression, and DNA methylation profiles, and to perform survival and gene correlation analyses of specific glioma subtypes. This database removes the barriers for researchers, providing rapid and convenient access to high-quality functional genomic data resources for biological studies and clinical applications. CGGA is available at http://www.cgga.org.cn.

Patients with high ALDH1A3-expressing glioblastoma (ALDH1A3^hi GBM) show limited benefit from postoperative chemoradiotherapy. Understanding the mechanisms underlying such resistance in these patients is crucial for the development of new treatments. Here, we show that the interaction between ALDH1A3 and PKM2 enhances the latter's tetramerization and promotes lactate accumulation in glioblastoma stem cells (GSCs). By scanning the lactylated proteome in lactate-accumulating GSCs, we show that XRCC1 undergoes lactylation at lysine 247 (K247). Lactylated XRCC1 shows a stronger affinity for importin α, allowing for greater nuclear transposition of XRCC1 and enhanced DNA repair. Through high-throughput screening of a small-molecule library, we show that D34-919 potently disrupts the ALDH1A3-PKM2 interaction, preventing the ALDH1A3-mediated enhancement of PKM2 tetramerization. In vitro and in vivo treatment with D34-919 enhanced chemoradiotherapy-induced apoptosis of GBM cells. Together, our findings show that ALDH1A3-mediated PKM2 tetramerization is a potential therapeutic target to improve the response to chemoradiotherapy in ALDH1A3^hi GBM.

Preoperative MRI is one of the most important clinical results for the diagnosis and treatment of glioma patients. The objective of this study was to construct a stable and validatable preoperative T₂-weighted MRI-based radiomics model for predicting the survival of gliomas. A total of 652 glioma patients across three independent cohorts were covered in this study including their preoperative T₂-weighted MRI images, RNA-seq and clinical data. Radiomic features (1731) were extracted from preoperative T₂-weighted MRI images of 167 gliomas (discovery cohort) collected from Beijing Tiantan Hospital and then used to develop a radiomics prediction model through a machine learning-based method. The performance of the radiomics prediction model was validated in two independent cohorts including 261 gliomas from the The Cancer Genomae Atlas database (external validation cohort) and 224 gliomas collected in the prospective study from Beijing Tiantan Hospital (prospective validation cohort). RNA-seq data of gliomas from discovery and external validation cohorts were applied to establish the relationship between biological function and the key radiomics features, which were further validated by single-cell sequencing and immunohistochemical staining. The 14 radiomic features-based prediction model was constructed from preoperative T₂-weighted MRI images in the discovery cohort, and showed highly robust predictive power for overall survival of gliomas in external and prospective validation cohorts. The radiomic features in the prediction model were associated with immune response, especially tumour macrophage infiltration. The preoperative T₂-weighted MRI radiomics prediction model can stably predict the survival of glioma patients and assist in preoperatively assessing the extent of macrophage infiltration in glioma tumours.

Background: PD-L1 has been widely reported as immune check points in a range of malignancies as well as some immune-originated diseases. In glioma, the role of PD-L1 remains unclear. We aimed at investigating its role at transcriptome level and relationship with clinical practice. Method and patients: In total, 976 glioma samples with transcriptome data, including 301 microarray data from Chinese Glioma Genome Atlas (CGGA project) and 675 RNAseq data from TCGA project, were enrolled into our study. Clinical and IDH mutation data were also available. R language was used as the main tool for statistical analysis and graphical work. Results: PD-L1 expression was found to be positively correlated with WHO grade of glioma. PD-L1 seemed to express more in mesenchymal subtype according to TCGA transcriptional classification scheme and may contribute as a potential marker for mesenchymal subtype in glioblastoma. Pearson correlation test indicated that PD-L1 showed robust correlation with PD1, PD-L2, and CD80 in CGGA dataset. Subsequent gene ontology analysis based on significantly correlated genes of PD-L1 revealed that PD-L1 seemed to be profoundly associated with T cell activation. To further investigate the relationship between PD-L1 expression and immune response, we selected a series of immune signatures, which were then transformed into metagenes, and found that PD-L1 expression was particularly paralleled with T-cells and macrophages-related immune response instead of B cell linage-related immune response. In line with the corresponding biological process, PD-L1 exhibited predictive value for glioma patients: Higher PD-L1 indicated significantly shorter survival, especially in glioblastoma. Conclusion: PD-L1 is upregulated in glioblastoma, and is synergistic with other check point members. Moreover, PD-L1 is significantly associated with T-cell activation and macrophage-related immune response and predicts much worse survival for patients, warranting clinical trials of PD1/PD-L1 checkpoint inhibitors for potential glioma treatment.

Isocitrate dehydrogenase (IDH) wild-type glioblastoma (GBM) has a dismal prognosis. A better understanding of tumor evolution holds the key to developing more effective treatment. Here we study GBM's natural evolutionary trajectory by using rare multifocal samples. We sequenced 61,062 single cells from eight multifocal IDH wild-type primary GBMs and defined a natural evolution signature (NES) of the tumor. We show that the NES significantly associates with the activation of transcription factors that regulate brain development, including MYBL2 and FOSL2. Hypoxia is involved in inducing NES transition potentially via activation of the HIF1A–FOSL2 axis. High-NES tumor cells could recruit and polarize bone marrow–derived macrophages through activation of the FOSL2–ANXA1–FPR1/3 axis. These polarized macrophages can efficiently suppress T-cell activity and accelerate NES transition in tumor cells. Moreover, the polarized macrophages could upregulate CCL2 to induce tumor cell migration. Significance: GBM progression could be induced by hypoxia via the HIF1A–FOSL2 axis. Tumor-derived ANXA1 is associated with recruitment and polarization of bone marrow–derived macrophages to suppress the immunoenvironment. The polarized macrophages promote tumor cell NES transition and migration. This article is highlighted in the In This Issue feature, p. 2711

Background: Researches on immunotherapy of glioma has been increasing exponentially in recent years. However, autoimmune-like side effects of current immune checkpoint blockade hindered the clinical application of immunotherapy in glioma. The discovery of the TIM-3, a tumor-specific immune checkpoint, has shed a new light on solution of this dilemma. We aimed at investigating the role of TIM-3 at transcriptome level and its relationship with clinical practice in glioma. Methods: A cohort of 325 glioma patients with RNA-seq data from Chinese Glioma Genome Atlas (CGGA project) was analyzed, and the results were well validated in TCGA RNA-seq data of 699 gliomas. R language was used as the main tool for statistical analysis and graphical work. Results: TIM-3 was enriched in glioblastoma (the most malignant glioma) and IDH-wildtype glioma. TIM-3 can act as a potential marker for mesenchymal molecular subtype according to TCGA transcriptional classification scheme in glioma. TIM-3 was closely related to immune functions in glioma, especially T cell mediated immune response to tumor cell and T cell mediated cytotoxicity directed against tumor cell target. Moreover, TIM-3 and PD-L1 played almost exactly the same inflammatory activation functions in glioma. Clinically, high expression of TIM-3 was an independent indicator of poor prognosis. Conclusion: The expression of TIM-3 is closely related to the pathology and molecular pathology of glioma. Meanwhile, in glioma TIM-3 plays a specific role in T cell tumor immune response. Therefore, TIM-3 is a promising target for immunotherapeutic strategies, providing an alternative treatment when glioma gains resistance to antibodies of PD-1/PD-L1.