Graphical AbstractHighlights d Characterization of the mutational landscape of secondary glioblastoma d Clonal and subclonal METex14 promote glioma progression and mark worse prognosis d PLB-1001 is a highly selective, efficient, and BBB-permeable MET kinase inhibitor d PLB-1001 provides a safe and efficacious therapeutic approach for glioma treatment SUMMARY Low-grade gliomas almost invariably progress into secondary glioblastoma (sGBM) with limited therapeutic option and poorly understood mechanism. By studying the mutational landscape of 188 sGBMs, we find significant enrichment of TP53 mutations, somatic hypermutation, MET-exon-14-skipping (METex14), PTPRZ1-MET (ZM) fusions, and MET amplification. Strikingly, METex14 frequently co-occurs with ZM fusion and is present in $14% of cases with significantly worse prognosis. Subsequent studies show that METex14 promotes glioma progression by prolonging MET activity. Furthermore, we describe a MET kinase inhibitor, PLB-1001, that demonstrates remarkable potency in selectively inhibiting MET-altered tumor cells in preclinical models. Importantly, this compound also shows blood-brain barrier permeability and is subsequently applied in a phase I clinical trial that enrolls MET-altered chemo-resistant glioma patients. Encouragingly, PLB-1001 achieves partial response in at least two advanced sGBM patients with rarely significant side effects, underscoring the clinical potential for precisely treating gliomas using this therapy.

No abstract

Brain tumors are among the most challenging human tumors for which the mechanisms driving progression and heterogeneity remain poorly understood. We combined single-cell RNA-seq with multisector biopsies to sample and analyze single-cell expression profiles of gliomas from 13 Chinese patients. After classifying individual cells, we generated a spatial and temporal landscape of glioma that revealed the patterns of invasion between the different sub-regions of gliomas. We also used single-cell inferred CNVs and pseudotime trajectories to inform on the crucial branches that dominate tumor progression. The dynamic cell components of the multi-region biopsy analysis allowed us to spatially deconvolute with unprecedented accuracy the transcriptomic features of the core and those of the periphery of glioma at single cell level. Through this rich and geographically detailed dataset, we were also able to characterize and construct the chemokine and chemokine receptor interactions that exist among different tumor and non-tumor cells. This study provides the first spatial-level analysis of the cellular states that characterize human gliomas. It also presents an initial molecular map of the cross-talks between glioma cells and the surrounding microenvironment with single cell resolution.

Brain tumors are among the most challenging human tumors for which the mechanisms driving progression and heterogeneity remain poorly understood. We combined single-cell RNA-seq with multisector biopsies to sample and analyze single-cell expression profiles of gliomas from 13 Chinese patients. After classifying individual cells, we generated a spatial and temporal landscape of glioma that revealed the patterns of invasion between the different sub-regions of gliomas. We also used single-cell inferred CNVs and pseudotime trajectories to inform on the crucial branches that dominate tumor progression. The dynamic cell components of the multi-region biopsy analysis allowed us to spatially deconvolute with unprecedented accuracy the transcriptomic features of the core and those of the periphery of glioma at single cell level. Through this rich and geographically detailed dataset, we were also able to characterize and construct the chemokine and chemokine receptor interactions that exist among different tumor and non-tumor cells. This study provides the first spatial-level analysis of the cellular states that characterize human gliomas. It also presents an initial molecular map of the crosstalks between glioma cells and the surrounding microenvironment with single cell resolution.

Although O(6)-methylguanine DNA methyltransferase (MGMT) promoter methylation status is an important marker for glioblastoma multiforme (GBM), there is considerable variability in the clinical outcome of patients with similar methylation profiles. The present study aimed to refine the prognostic and predictive value of MGMT promoter status in GBM by identifying a micro (mi)RNA risk signature. Data from The Cancer Genome Atlas was used for this study, with MGMT promoter-methylated samples randomly divided into training and internal validation sets. Data from The Chinese Glioma Genome Atlas was used for independent validation. A five miRNA-based risk signature was established for MGMT promoter-methylated GBM to distinguish cases as high- or low-risk with distinct prognoses, which was confirmed using internal and external validation sets. Importantly, the prognostic value of the signature was significant in different cohorts stratified by clinicopathologic factors and alkylating chemotherapy, and a multivariate Cox analysis found it to be an independent prognostic marker along with age and chemotherapy. Based on these three factors, we developed a quantitative model with greater accuracy for predicting the 1-year survival of patients with MGMT promoter-methylated GBM. These results indicate that the five-miRNA signature is an independent risk predictor for GBM with MGMT promoter methylation and can be used to identify patients at high risk of unfavorable outcome and resistant to alkylating chemotherapy, underscoring its potential for personalized GBM management.

The methylation status of the promoter of MGMT gene is a crucial factor influencing clinical decision-making in patients with gliomas. MGMT pyrosequencing results are often dichotomized by a cut-off value based on an average of several tested CpGs. However, this method frequently results in a "gray zone", representing a dilemma for physicians. We therefore propose a novel analytical model for MGMT methylation pyrosequencing. MGMT CpG heterogeneity was investigated in 213 glioma patients in two tested cohorts: cohort A in which CpGs 75-82 were tested and cohort B in which CpGs 72-78 were tested. The predictive performances of the novel and traditional averaging models were compared in 135 patients who received temozolomide using receiver operating characteristic curves and Kaplan-Meier curves, and in patients stratified according to isocitrate dehydrogenase gene mutation status. The results were validated in an independent cohort of 65 consecutive patients with high-grade gliomas from the Chinese Glioma Genome Atlas database. Heterogeneity of MGMT promoter CpG methylation level was observed in most gliomas. The optimal cut-off value for each individual CpG varied from 4-16%. The current analysis defined MGMT promoter methylation as occurring when at least three CpGs exceeded their respective cut-off values. This novel analysis could accurately predict the prognosis of patients in the methylation "gray zone" according to the standard averaging method, and improved the area under the curves from 0.67, 0.76, and 0.67 to 0.70, 0.84, and 0.72 in cohorts A, B, and the validation cohort, respectively, demonstrating superiority of this analytical method in all three cohorts. Furthermore, the advantages of the novel analysis were retained regardless of WHO grade and isocitrate dehydrogenase gene mutation status. In conclusion, this novel analytical model offers an improved clinical predictive performance for MGMT pyrosequencing results and is suitable for clinical use in patients with gliomas.

BackgroundIncreasing evidence has shown that long non-coding RNAs (lncRNAs) are important prognostic biomarkers and epigenetic regulators with critical roles in cancer initiation and progression. However, the expression and clinical prognostic value of antisense lncRNAs in diffuse glioma patients remain unknown.MethodsHere, we profiled differentially expressed antisense lncRNAs in glioma using RNA sequencing data from Chinese Glioma Genome Atlas database. Cox regression was performed to evaluate the prognostic value. Gene oncology (GO) and gene set enrichment analysis (GSEA) were used for functional analysis of antisense LncRNAs.ResultsExpression profiling identified 169 aberrantly expressed antisense lncRNAs between lower grade glioma (LGG) (grade II and III) and glioblastoma multiforme (GBM), 113 antisense lncRNAs between LGG IDH-wt and IDH-mut samples, and 70 antisense lncRNAs between GBM IDH-wt and IDH-mut samples, respectively. Among them, three antisense lncRNAs (WDFY3-AS2, MCM3AP-AS1 and LBX2-AS1) were significantly associated with prognosis and malignant progression of patients. WDFY3-AS2, the top one of downregulated antisense lncRNAs in GBM with fold change of 0.441 (P < 0.001), showed specific decreased expression in classical, mesenchymal, LGG IDH-wt, GBM IDH-wt or MGMT promoter unmethylated stratified patients. Chi square test found that WDFY3-AS2 was significantly associated with the clinical and molecular features of glioma. Univariate and multivariate Cox regression analysis indicated that WDFY3-AS2 was independently correlated with overall survival (OS) of patients. Kaplan–Meier analysis found that patients with high WDFY3-AS2 expression had longer OS than the low expression ones in the stratified cohorts. Additionally, GO and GSEA showed that gene sets correlated with WDFY3-AS2 expression were involved in regulation of synaptic transmission, glutamate receptor and TNF signaling pathway.ConclusionOur findings provided convincing evidence that WDFY3-AS2 is a novel valuable prognostic biomarker for diffuse glioma.Electronic supplementary materialThe online version of this article (10.1186/s12935-018-0603-2) contains supplementary material, which is available to authorized users.

1p/19q codeletion, which leads to the abnormal expression of 1p19q genes in oligodendroglioma, is associated with chemosensitivity and favorable prognosis. Here, we aimed to explore the clinical implications of 1p19q gene expression in 1p/19q non-codel gliomas. We analyzed expression of 1p19q genes in 668 1p/19q non-codel gliomas obtained from The Cancer Genome Atlas (n = 447) and the Chinese Glioma Genome Atlas (n = 221) for training and validation, respectively. The expression of 1p19q genes was significantly correlated with the clinicopathological features and overall survival of 1p/19q non-codel gliomas. Then, we derived a risk signature of 25 selected 1p19q genes that not only had prognosis value in total 1p/19q non-codel gliomas but also had prognosis value in stratified gliomas. The prognosis value of the risk signature was superior than known clinicopathological features in 1p/19q non-codel gliomas and was also highly associated with the following features: loss of CDKN2A/B copy number in mutant-IDH-astrocytoma; telomerase reverse transcriptase (TERT) promoter mutation, combined chromosome 7 gain/chromosome 10 loss and epidermal growth factor receptor amplification in wild-type-IDH-astrocytoma; classical and mesenchymal subtypes in glioblastoma. Furthermore, genes enriched in the biological processes of cell division, extracellular matrix, angiogenesis significantly correlated to the signature risk score, and this is also supported by the immunohistochemistry and cell biology experiments. In conclusion, the expression profile of 1p19q genes is highly associated with the malignancy and prognosis of 1p/19q non-codel gliomas. A 25-1p19q-gene signature has powerfully predictive value for both malignant molecular pathological features and prognosis across distinct subgroups of 1p/19q non-codel gliomas.