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.

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

NoncoRNA (http://www.ncdtcdb.cn:8080/NoncoRNA/) is a manually curated database of experimentally supported non-coding RNAs (ncRNAs) and drug target associations that aim to potentially provide a high-quality data resource for exploring drug sensitivity/resistance-related ncRNAs in various human cancers. ncRNA are RNA molecular that do not encode proteins, but are involved in gene regulation and cellular functions in variety of human diseases, including neurodegenerative diseases and cancers. Here, we developed NoncoRNA which contained 8233 entries between 5568 ncRNAs and 154 drugs in 134 cancers. Each entry in the NoncoRNA contains detailed information on the ncRNAs, drugs, and cancers, the ncRNA expression pattern and experimental detection techniques, drug response and other targets, literature references, and other information. NoncoRNA offers a user-friendly, open access web interface to easily browse, search, and download data. NoncoRNA also provides a submission page for researchers to submit newly validated ncRNA-drug-cancer associations. NoncoRNA might serve as an immeasurable resource for understanding the roles of ncRNAs in cancer therapy.

Transcriptomic data derived from bulk sequencing have been applied to delineate the tumor microenvironment (TME) and define immune subtypes in various cancers, which may facilitate the design of immunotherapy treatment strategies. We herein gathered published gene expression data from diffuse lower‐grade glioma (LGG) patients to identify immune subtypes. Based on the immune gene profiles of 402 LGG patients from The Cancer Genome Atlas, we performed consensus clustering to determine robust clusters of patients, and evaluated their reproducibility in three Chinese Glioma Genome Atlas cohorts. We further integrated immunogenomics methods to characterize the immune environment of each subtype. Our analysis identified and validated three immune subtypes—Im1, Im2, and Im3—characterized by differences in lymphocyte signatures, somatic DNA alterations, and clinical outcomes. Im1 had a higher infiltration of CD8+ T cells, Th17, and mast cells. Im2 was defined by elevated cytolytic activity, exhausted CD8+ T cells, macrophages, higher levels of aneuploidy, and tumor mutation burden, and these patients had worst outcome. Im3 displayed more prominent T helper cell and APC coinhibition signatures, with elevated pDCs and macrophages. Each subtype was associated with distinct somatic alterations. Moreover, we applied graph structure learning‐based dimensionality reduction to the immune landscape and revealed significant intracluster heterogeneity with Im2 subtype. Finally, we developed and validated an immune signature with better performance of prognosis prediction. Our results demonstrated the immunological heterogeneity within diffuse LGG and provided valuable stratification for the design of future immunotherapy.

CircRNADisease v2.0 is an enhanced and reliable database that offers experimentally verified relationships between circular RNAs (circRNAs) and various diseases. It is accessible at http://cgga.org.cn/circRNADisease/ or http://cgga.org.cn:9091/circRNADisease/. The database currently includes 6998 circRNA-disease entries across multiple species, representing a remarkable 19.77-fold increase compared to the previous version. This expansion consists of a substantial rise in the number of circRNAs (from 330 to 4246), types of diseases (from 48 to 330) and covered species (from human only to 12 species). Furthermore, a new section has been introduced in the database, which collects information on circRNA-associated factors (genes, proteins and microRNAs), molecular mechanisms (molecular pathways), biological functions (proliferation, migration, invasion, etc.), tumor and/or cell line and/or patient-derived xenograft (PDX) details, and prognostic evidence in diseases. In addition, we identified 7 159 865 relationships between mutations and circRNAs among 30 TCGA cancer types. Due to notable enhancements and extensive data expansions, the circRNADisease 2.0 database has become an invaluable asset for both clinical practice and fundamental research. It enables researchers to develop a more comprehensive understanding of how circRNAs impact complex diseases.

Clonal evolution drives cancer progression and therapeutic resistance. Recent studies have revealed divergent longitudinal trajectories in gliomas, but early molecular features steering posttreatment cancer evolution remain unclear. Here, we collected sequencing and clinical data of initial-recurrent tumor pairs from 544 adult diffuse gliomas and performed multivariate analysis to identify early molecular predictors of tumor evolution in three diffuse glioma subtypes. We found that CDKN2A deletion at initial diagnosis preceded tumor necrosis and microvascular proliferation that occur at later stages of IDH-mutant glioma. Ki67 expression at diagnosis was positively correlated with acquiring hypermutation at recurrence in the IDH–wild-type glioma. In all glioma subtypes, MYC gain or MYC- target activation at diagnosis was associated with treatment-induced hypermutation at recurrence. To predict glioma evolution, we constructed CELLO2 (Cancer EvoLution for LOngitudinal data version 2), a machine learning model integrating features at diagnosis to forecast hypermutation and progression after treatment. CELLO2 successfully stratified patients into subgroups with distinct prognoses and identified a high-risk patient group featured by MYC gain with worse post-progression survival, from the low-grade IDH-mutant-noncodel subtype. We then performed chronic temozolomide-induction experiments in glioma cell lines and isogenic patient-derived gliomaspheres and demonstrated that MYC drives temozolomide resistance by promoting hypermutation. Mechanistically, we demonstrated that, by binding to open chromatin and transcriptionally active genomic regions, c-MYC increases the vulnerability of key mismatch repair genes to treatment-induced mutagenesis, thus triggering hypermutation. This study reveals early predictors of cancer evolution under therapy and provides a resource for precision oncology targeting cancer dynamics in diffuse gliomas.