A novel analytical model of MGMT methylation pyrosequencing offers improved predictive performance in patients with gliomas

Unveiling hierarchy and spatial distribution of O6-methylguanine-DNA methyltransferase promoter methylation in World Health Organization grade 2-3 gliomas

This study investigated the role of O6-methylguanine-DNA methyltransferase promoter (MGMTp) methylation hierarchy and heterogeneity in grade 2-3 gliomas, focusing on variations in chemotherapy benefits and resection dependency. A cohort of 668 newly diagnosed grade 2-3 gliomas, with comprehensive clinical, radiological, and molecular data, formed the basis of this analysis. The extent of resection was categorized into gross total resection (GTR ≥100%), subtotal resection (STR >90%), and partial resection (PR ≤90%). MGMTp methylation levels were examined using quantitative pyrosequencing. Our findings highlighted the critical role of GTR in improving the prognosis for astrocytomas (IDH1/2-mutant and 1p/19q non-codeleted), contrasting with its lesser significance for oligodendrogliomas (IDH1/2 mutation and 1p/19q codeletion). Oligodendrogliomas demonstrated the highest average MGMTp methylation levels (median: 28%), with a predominant percentage of methylated cases (average methylation levels >20%). Astrocytomas were more common in the low-methylated group (10%-20%), while IDH wild-type gliomas were mostly unmethylated (<10%). Spatial distribution analysis revealed a decrement in frontal lobe involvement from methylated, low-methylated to unmethylated cases (72.8%, 59.3%, and 47.8%, respectively). In contrast, low-methylated and unmethylated cases were more likely to invade the temporal-insular region (19.7%, 34.3%, and 40.4%, respectively). Astrocytomas with intermediate MGMTp methylation were notably associated with temporal-insular involvement, potentially indicating a moderate response to temozolomide and underscoring the importance of aggressive resection strategies. In conclusion, our study elucidates the complex interplay of MGMTp methylation hierarchy and heterogeneity among grade 2-3 gliomas, providing insights into why astrocytomas and IDH wild-type lower-grade glioma might derive less benefit from chemotherapy.

Accurate and rapid molecular subgrouping of high-grade glioma via deep learning-assisted label-free fiber-optic Raman spectroscopy

Molecular genetics is highly related with prognosis of high-grade glioma. Accordingly, the latest WHO guideline recommends that molecular subgroups of the genes, including IDH, 1p/19q, MGMT, TERT, EGFR, Chromosome 7/10, CDKN2A/B, need to be detected to better classify glioma and guide surgery and treatment. Unfortunately, there is no preoperative or intraoperative technology available for accurate and comprehensive molecular subgrouping of glioma. Here, we develop a deep learning-assisted fiber-optic Raman diagnostic platform for accurate and rapid molecular subgrouping of high-grade glioma. Specifically, a total of 2,354 fingerprint Raman spectra was obtained from 743 tissue sites (astrocytoma: 151; oligodendroglioma: 150; glioblastoma (GBM): 442) of 44 high-grade glioma patients. The convolutional neural networks (ResNet) model was then established and optimized for molecular subgrouping. The mean area under receiver operating characteristic curves (AUC) for identifying the molecular subgroups of high-grade glioma reached 0.904, with mean sensitivity of 83.3%, mean specificity of 85.0%, mean accuracy of 83.3%, and mean time expense of 10.6 s. The diagnosis performance using ResNet model was shown to be superior to PCA-SVM and UMAP models, suggesting that high dimensional information from Raman spectra would be helpful. In addition, for the molecular subgroups of GBM, the mean AUC reached 0.932, with mean sensitivity of 87.8%, mean specificity of 83.6%, and mean accuracy of 84.1%. Furthermore, according to saliency maps, the specific Raman features corresponding to tumor-associated biomolecules (e.g. nucleic acid, tyrosine, tryptophan, cholesteryl ester, fatty acid, and collagen) were found to contribute to the accurate molecular subgrouping. Collectively, this study opens up new opportunities for accurate and rapid molecular subgrouping of high-grade glioma, which would assist optimal surgical resection and instant post-operative decision-making.

Genome-wide DNA methylation analysis identifies potent CpG signature for temzolomide response in non-G-CIMP glioblastomas with unmethylated MGMT promoter: MGMT-dependent roles of GPR81

PURPOSES

To identify potent DNA methylation candidates that could predict response to temozolomide (TMZ) in glioblastomas (GBMs) that do not have glioma-CpGs island methylator phenotype (G-CIMP) but have an unmethylated promoter of O-6-methylguanine-DNA methyltransferase (unMGMT).

METHODS

The discovery-validation approach was planned incorporating a series of G-CIMP-/unMGMT GBM cohorts with DNA methylation microarray data and clinical information, to construct multi-CpG prediction models. Different bioinformatic and experimental analyses were performed for biological exploration.

RESULTS

By analyzing discovery sets with radiotherapy (RT) plus TMZ versus RT alone, we identified a panel of 64 TMZ efficacy-related CpGs, from which a 10-CpG risk signature was further constructed. Both the 64-CpG panel and the 10-CpG risk signature were validated showing significant correlations with overall survival of G-CIMP-/unMGMT GBMs when treated with RT/TMZ, rather than RT alone. The 10-CpG risk signature was further observed for aiding TMZ choice by distinguishing differential outcomes to RT/TMZ versus RT within each risk subgroup. Functional studies on GPR81, the gene harboring one of the 10 CpGs, indicated its distinct impacts on TMZ resistance in GBM cells, which may be dependent on the status of MGMT expression.

CONCLUSIONS

The 64 TMZ efficacy-related CpGs and in particular the 10-CpG risk signature may serve as promising predictive biomarker candidates for guiding optimal usage of TMZ in G-CIMP-/unMGMT GBMs.

Biomedical applications of terahertz spectra in clinical and molecular pathology of human glioma

Gliomas are the most common type of primary tumor originating in the central nervous system of adults. Tumor histological type, pathological grade, and molecular pathology are significant prognosis and predictive factors. In this study, we were aiming to predict histological type and molecular pathological features based on terahertz time-domain spectroscopy technology. Nine gliomas with different grades, one meningioma, and one lymphoma were enrolled. There were significant differences in terahertz absorption coefficient between normal brain tissue, tumoral-periphery, and tumoral-center tissue in specific frequency bands (0.2-1.4 THz). Histological type, pathological grade, and glioma-specific biomarkers were closely related to the terahertz absorption coefficient in both tumoral-periphery and tumoral-center tissues. Interestingly, tumoral-periphery showed more obvious differences than tumoral-center tissues in almost all aspects. All the results show that the terahertz technology has potential application value in the intraoperative real-time glioma recognition and diagnosis of glioma histological and molecular pathological features.

Molecular pathology and clinical implications of diffuse glioma

ABSTRACT

The prognosis for diffusely infiltrating gliomas at World Health Organization (WHO) grade 2-4 remains dismal due to their heterogeneity. The rapid development of genome-wide molecular-profiling-associated studies has greatly promoted the accuracy of glioma classification. Thus, the latest version of the WHO classification of the central nervous system tumors published in 2021 has incorporated more molecular biomarkers together with histological features for the diagnosis of gliomas. Advanced usage of molecular pathology in clinical diagnostic practice provides also new opportunities for the therapy of patients with glioma, including surgery, radiotherapy and chemotherapy, targeted therapy, immunotherapy, and more precision clinical trials. Herein, we highlight the updates in the classification of gliomas according to the latest WHO guidelines and summarize the clinically relevant molecular markers by focusing on their applications in clinical practice. We also review the advances in molecular features of gliomas, which can facilitate the development of glioma therapies, thereby discussing the challenges and future directions of molecular pathology toward precision medicine for patients with glioma.

Volumetric Analysis of Glioblastoma for Determining Which CpG Sites Should Be Tested by Pyrosequencing to Predict Temozolomide Efficacy

The aim of the present study was to determine which individual or combined CpG sites among O6-methylguanine DNA methyltransferase CpG 74–89 in glioblastoma mainly affects the response to temozolomide resulting from CpG methylation using statistical analyses focused on the tumor volume ratio (TVR). We retrospectively examined 44 patients who had postoperative volumetrically measurable residual tumor tissue and received adjuvant temozolomide therapy for at least 6 months after initial chemoradiotherapy. TVR was defined as the tumor volume 6 months after the initial chemoradiotherapy divided by that before the start of chemoradiotherapy. Predictive values for TVR as a response to adjuvant therapy were compared among the averaged methylation percentages of individual or combined CpGs using the receiver operating characteristic curve. Our data revealed that combined CpG 78 and 79 showed a high area under the curve (AUC) and a positive likelihood ratio and that combined CpG 76–79 showed the highest AUC among all combinations. AUCs of consecutive CpG combinations tended to be higher for CpG 74–82 in exon 1 than for CpG 83–89 in intron 1. In conclusion, the methylation status at CpG sites in exon 1 was strongly associated with TVR reduction in glioblastoma.

MGMT and Whole-Genome DNA Methylation Impacts on Diagnosis, Prognosis and Therapy of Glioblastoma Multiforme

Epigenetic changes in DNA methylation contribute to the development of many diseases, including cancer. In glioblastoma multiforme, the most prevalent primary brain cancer and an incurable tumor with a median survival time of 15 months, a single epigenetic modification, the methylation of the O⁶-Methylguanine-DNA Methyltransferase (MGMT) gene, is a valid biomarker for predicting response to therapy with alkylating agents and also, independently, prognosis. More recently, the progress from single gene to whole-genome analysis of DNA methylation has allowed a better subclassification of glioblastomas. Here, we review the clinically relevant information that can be obtained by studying MGMT gene and whole-genome DNA methylation changes in glioblastomas, also highlighting benefits, including those of liquid biopsy, and pitfalls of the different detection methods. Finally, we discuss how changes in DNA methylation, especially in glioblastomas bearing mutations in the Isocitrate Dehydrogenase (IDH) 1 and 2 genes, can be exploited as targets for tailoring therapy.

Expression and Potential Biomarkers of Regulators for M7G RNA Modification in Gliomas

Gliomas are the most frequent primary malignant brain tumors of the central nervous system, causing significant impairment and death. There is mounting evidence that N7 methylguanosine (m7G) RNA dysmethylation plays a significant role in the development and progression of cancer. However, the expression patterns and function of the m7G RNA methylation regulator in gliomas are yet unknown. The goal of this study was to examine the expression patterns of 31 critical regulators linked with m7G RNA methylation and their prognostic significance in gliomas. To begin, we systematically analyzed patient clinical and prognostic data and mRNA gene expression data from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. We found that 17 key regulators of m7G RNA methylation showed significantly higher expression levels in gliomas. We then divided the sample into two subgroups by consensus clustering. Cluster 2 had a poorer prognosis than cluster 1 and was associated with a higher histological grade. In addition, cluster 2 was significantly enriched for cancer-related pathways. Based on this discovery, we developed a risk model involving three m7G methylation regulators. Patients were divided into high-risk and low-risk groups based on risk scores. Overall survival (OS) was significantly lower in the high-risk group than in the low-risk group. Further analysis showed that the risk score was an independent prognostic factor for gliomas.

MGMT gene promoter methylation by pyrosequencing method correlates volumetric response and neurological status in IDH wild-type glioblastomas

Purpose

Although the usefulness of O⁶-methylguanine DNA methyltransferase (MGMT) promoter methylation analysis for predicting response to chemoradiotherapy and the prognosis of patients with glioblastoma has been widely reported, there is still no consensus regarding how to define MGMT promoter methylation percentage (MGMTpm%) cutoffs by pyrosequencing method. The aim of this study was to determine the optimal cutoff value of MGMT promoter methylation status using volumetric analysis focused on the tumor volume ratio (TVR) measured by MRI.

Methods

This retrospective study included newly diagnosed IDH wild-type glioblastoma patients with residual tumor after surgery, followed by local radiotherapy with temozolomide. TVR was defined as the tumor volume at 6 months after the initial chemoradiotherapy administration divided by the tumor volume before the start of therapy. The mean MGMTpm% of 16 CpG islands (74–89) was analyzed using pyrosequencing. We statistically analyzed the correlation between MGMTpm%, TVR, and change in Karnofsky performance status.

Results

The study included 44 patients with residual tumors. Thirteen (92.9%) of 14 patients with MGMTpm% ≥ 23.9% showed 50% or more volumetric response, leading to prolonged survival, and 17 (70.8%) of 24 patients with MGMTpm% < 8.2% had progressive disease after initial chemoradiotherapy administration. Three (50.0%) of six patients with MGMTpm% 8.2% to < 23.9% had stable disease or partial response.

Conclusion

Evaluation of MGMTpm% by pyrosequencing is important in predicting the volumetric response and prognosis of glioblastoma patients with residual tumors.

Should temozolomide be used on the basis of O6-methylguanine DNA methyltransferase status in patients with advanced neuroendocrine tumors? A systematic review and meta-analysis

BACKGROUND

Temozolomide (TEM) is an active treatment in metastatic neuroendocrine tumors (NETs). Patients affected by glioblastoma multiforme or advanced melanoma treated with TEM who have deficiency of O⁶-methylguanine DNA methyltransferase (MGMT) have a better responses and survival. However, the predictive role of MGMT in patients with NETs treated with TEM is still debated.

METHODS

We conducted a systematic review of the literature and meta-analysis, based on PRISMA methodology, searching in the main databases (PubMed, Embase, Scopus, Web of Science, Cochrane Library and clinical trial.gov) and the proceedings of the main international congresses, until April 26, 2021.

RESULTS

Twelve out of 616 articles were selected for our analysis, regarding a total of 858 NET patients treated with TEM-based chemotherapy. The status of MGMT had been tested in 513 (60%) patients, using various methods. The pooled overall response rate (ORR) was higher in MGMT-deficient compared with MGMT-proficient NETs, with a risk difference of 0.31 (95% confidence interval, CI: 0.13-0.50; p < 0.001; I²: 73%) and risk ratio of 2.29 (95% CI: 1.34-3.91; p < 0.001; I²: 55%). The pooled progression free survival (PFS) (hazard ratio, HR = 0.56; 95% CI: 0.43-0.74; p < 0.001) and overall survival (OS) (HR = 0.41; 95% CI: 0.20-0.62; p = 0.011) were longer in MGMT-deficient versus MGMT-proficient NETs.

CONCLUSIONS

Our meta-analysis suggested that MGMT status may be predictive of TEM efficacy. However, due to the high heterogeneity of the evaluated studies the risk of biases should be considered. On this hypothesis future homogeneous prospective studies are warranted.

Combining MGMT promoter pyrosequencing and protein expression to optimize prognosis stratification in glioblastoma

Pyrosequencing (PSQ) represents the golden standard for MGMT promoter status determination. Binary interpretation of results based on the threshold from the average of several CpGs tested would neglect the existence of the “gray zone”. How to define the gray zone and reclassify patients in this subgroup remains to be elucidated. A consecutive cohort of 312 primary glioblastoma patients were enrolled. CpGs 74‐81 in the promoter region of MGMT were tested by PSQ and the protein expression was assessed by immunohistochemistry (IHC). Receiver operating characteristic curves were constructed to calculate the area under the curves (AUC). Kaplan‐Meier plots were used to estimate the survival rate of patients compared by the log‐rank test. The optimal threshold of each individual CpG differed from 5% to 11%. Patients could be separated into the hypomethylated subgroup (all CpGs tested below the corresponding optimal thresholds, n = 126, 40.4%), hypermethylated subgroup (all CpGs tested above the corresponding optimal thresholds, n = 108, 34.6%), and the gray zone subgroup (remaining patients, n = 78, 25.0%). Patients in the gray zone harbored an intermediate prognosis. The IHC score instead of the average methylation levels could successfully predict the prognosis for the gray zone (AUC for overall survival, 0.653 and 0.519, respectively). Combining PSQ and IHC significantly improved the efficiency of survival prediction (AUC: 0.662, 0.648, and 0.720 for PSQ, IHC, and combined, respectively). Immunohistochemistry is a robust method to predict prognosis for patients in the gray zone defined by PSQ. Combining PSQ and IHC could significantly improve the predictive ability for clinical outcomes.

A Signature of Nine lncRNA Methylated Genes Predicts Survival in Patients With Glioma

Glioma is one of the most common malignant tumors of the central nervous system, and its prognosis is extremely poor. Aberrant methylation of lncRNA promoter region is significantly associated with the prognosis of glioma patients. In this study, we investigated the potential impact of methylation of lncRNA promoter region in glioma patients to establish a signature of nine lncRNA methylated genes for determining glioma patient prognosis. Methylation data and clinical follow-up data were obtained from The Cancer Genome Atlas (TCGA). The multistep screening strategy identified nine lncRNA methylated genes that were significantly associated with the overall survival (OS) of glioma patients. Subsequently, we constructed a risk signature that containing nine lncRNA methylated genes. The risk signature successfully divided the glioma patients into high-risk and low-risk groups. Compared with the low-risk group, the high-risk group had a worse prognosis, higher glioma grade, and older age. Furthermore, we identified two lncRNAs termed PCBP1-AS1 and LINC02875 that may be involved in the malignant progression of glioma cells by using the TCGA database. Loss-of-function assays confirmed that knockdown of PCBP1-AS1 and LINC02875 inhibited the proliferation, migration, and invasion of glioma cells. Therefore, the nine lncRNA methylated genes signature may provide a novel predictor and therapeutic target for glioma patients.

Clinical validation of a novel quantitative assay for the detection of MGMT methylation in glioblastoma patients

Background

The promoter hypermethylation of the methylguanine-DNA methyltransferase gene is a frequently used biomarker in daily clinical practice as it is associated with a favorable prognosis in glioblastoma patients treated with temozolamide. Due to the absence of adequately standardized techniques, international harmonization of the MGMT methylation biomarker is still an unmet clinical need for the diagnosis and treatment of glioblastoma patients.

Results

In this study we carried out a clinical validation of a quantitative assay for MGMT methylation detection by comparing a novel quantitative MSP using double-probe (dp_qMSP) with the conventional MSP in 100 FFPE glioblastoma samples. We performed both technologies and established the best cutoff for the identification of positive-methylated samples using the quantitative data obtained from dp_qMSP. Kaplan–Meier curves and ROC time dependent curves were employed for the comparison of both methodologies.

Conclusions

We obtained similar results using both assays in the same cohort of patients, in terms of progression free survival and overall survival according to Kaplan–Meier curves. In addition, the results of ROC(t) curves showed that dp_qMSP increases the area under curve time-dependent in comparison with MSP for predicting progression free survival and overall survival over time. We concluded that dp_qMSP is an alternative methodology compatible with the results obtained with the conventional MSP. Our assay will improve the therapeutic management of glioblastoma patients, being a more sensitive and competitive alternative methodology that ensures the standardization of the MGMT-biomarker making it reliable and suitable for clinical use.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01044-2.

Predictive value of MGMT promoter methylation on the survival of TMZ treated IDH-mutant glioblastoma

Objective

O6methylguanine-DNA methyltransferase (MGMT) promoter methylation is a biomarker widely used to predict the sensitivity of IDH-wildtype glioblastoma to temozolomide therapy. Given that the IDH status has critical effects on the survival and epigenetic features of glioblastoma, we aimed to assess the role of MGMT promoter methylation in IDH-mutant glioblastoma.

Methods

This study included 187 IDH-mutant glioblastomas and used 173 IDH-wildtype glioblastomas for comparison. Kaplan-Meier curves and multivariate Cox regression were used to study the predictive effects.

Results

Compared with IDH-wildtype glioblastomas, IDH-mutant glioblastomas showed significantly higher (P < 0.0001) MGMT promoter methylation. We demonstrated that MGMT promoter methylation status, as determined by a high cutoff value (≥30%) in pyrosequencing, could be used to significantly stratify the survival of 50 IDH-mutant glioblastomas receiving temozolomide therapy (cohort A); this result was validated in another cohort of 25 IDH-mutant glioblastomas (cohort B). The median progression-free survival and median overall survival in cohort A were 9.33 and 13.76 months for unmethylated cases, and 18.37 and 41.61 months for methylated cases, and in cohort B were 6.97 and 9.10 months for unmethylated cases, and 23.40 and 26.40 months for methylated cases. In addition, we confirmed that the MGMT promoter methylation was significantly (P = 0.0001) correlated with longer OS in IDH-mutant patients with GBM, independently of age, gender distribution, tumor type (primary or recurrent/secondary), and the extent of resection.

Conclusions

MGMT promoter methylation has predictive value in IDH-mutant glioblastoma, but its cutoff value should be higher than that for IDH-wildtype glioblastoma.

[ALKBH5 suppresses migration and invasion of human trophoblast cells by inhibiting epithelial-mesenchymal transition]

OBJECTIVE

To investigate the effects of ALKBH5 on migration, invasion and epithelial-mesenchymal transition (EMT) of human trophoblast cells.

METHODS

The expression plasmid of ALKBH5 or a negative control plasmid (ALKBH5-NC) was transfected in human trophoblast HTR-8 /SVneo cells, and the expressions of ALKBH5 mRNA and protein were detected by qRT-PCR and Western blotting. Transwell assay was used to assess the changes in migration and invasion abilities of the trophoblast cells after the transfection. Western blotting was performed to detect the expressions of EMT-related proteins in the cells including vimentin, fibronectin, E-cadherin, N-cadherin, MMP9 and MMP2.

RESULTS

ALKBH5 mRNA and protein expressions were significantly higher in ALKBH5 group than in the control group (P < 0.05). Over-expression of ALKBH5 significantly attenuated migration and invasion abilities of HTR-8/Svneo cells (P < 0.05). Compared with the control cells, the cells overexpressing ALKBH5 showed an up-regulated expression of E-cadherin and down-regulated expressions of vimentin, fibronectin, N-cadherin, MMP9 and MMP2 (P < 0.05).

CONCLUSIONS

ALKBH5 is involved in the pathogenesis of preeclampsia by inhibiting EMT of trophoblast cells and hence reducing their migration and invasion abilities.

Transcriptional Characteristics of IDH-Wild Type Glioma Subgroups Highlight the Biological Processes Underlying Heterogeneity of IDH-Wild Type WHO Grade IV Gliomas

Isocitric dehydrogenase (IDH)-wild type diffuse gliomas, which have a poorer prognosis than their IDH-mutant counterparts, are also accompanied with high heterogeneity. Here, we aimed to identify the key biological processes associated with the three groups of IDH-wild type diffuse gliomas in 323 patients. By The Consortium to Inform Molecular and Practical Approaches to CNS Tumor Taxonomy (cIMPACT-NOW) update 3 recommendation, the three groups are Group A, diffuse astrocytic glioma, World Health Organization (WHO) grade II/III; Group B, diffuse astrocytic glioma, with one (or more) of the three genetic alterations: TERT promoter mutation, EGFR gene amplification, gain of chromosome 7 combined with loss of chromosome 10, WHO grade IV; and Group C, glioblastoma, WHO grade IV. Consistent with their histologic and genetic molecular features, we successfully identified that biological activities associated with “cell cycle” and “cell mitosis” are significantly elevated in Group B compared with Group A; microenvironment-related hallmarks “angiogenesis” and “hypoxia,” and biological processes of “extracellular matrix,” “immune response,” and “positive regulation of transcriptional activities” were more enriched in Group C than Group B. We also constructed a nine-gene signature from differentially expressed genes among the three groups to further stratify the WHO grade IV gliomas (Groups B and C) whose survival cannot be clearly stratified by current classification systems. This signature was an independent prognosis factor for WHO grade IV gliomas and had better prognostic value than other known factors in both training and validation dataset. In addition, the signature risk score was positively correlated with the amount of infiltrated immune cells, expression of immune checkpoints, and the genes enriched in biological processes of “immune response,” “cell cycle,” and “extracellular matrix.” The bioinformatic analysis results were also validated by immunohistochemistry and patient-derived cell proliferation assay. Overall, our findings revealed the key biological processes underlying the new classifications of IDH-wild type diffuse glioma. Meanwhile, we constructed a signature, which could properly stratify the prognosis, cell proliferation activates, extracellular matrix-mediated biological activities, and immune-microenvironment of IDH-wild type WHO grade IV gliomas.

Influence of Intratumor Heterogeneity on the Predictivity of MGMT Gene Promoter Methylation Status in Glioblastoma

Glioblastoma is the most aggressive tumor of the central nervous system. Prognosis is poor, even in the presence of a methylated state of MGMT gene promoter, which represents the biomarker with the highest prognostic/predictive value for the standard treatment of patients. Among patients with a methylated MGMT status, we identified an intermediate range of methylation above the standard 9% cut-off (gray zone) in which the predictive strength of the marker was lost. In an effort to improve the evaluation of the biomarker in clinical decision-making, we are carrying out a retrospective study, performing an in-depth analysis of samples used for diagnosis to understand how molecular heterogeneity, a hallmark of glioblastoma, impacts the evaluation of MGMT gene promoter methylation. Preliminary data from samples belonging to the "gray zone" tend to confirm the hypothesis of a mismatch between methylation values used for clinical decision-making and those included in our in-depth analysis. Confirmation of these data would help to better define the predictive power of MGMT promoter methylation status and greatly facilitate clinical decision-making.

Long-term efficacy of surgical resection with or without adjuvant therapy for treatment of secondary glioblastoma in adults

Background

There are limited studies on treatment strategies and associated clinical outcomes in patients with secondary glioblastoma (sGBM). We sought to investigate the prognostic factors and treatment decisions in a retrospective cohort of patients with sGBM.

Methods

One hundred and seventy-one patients with sGBM who met the screening criteria were included in this study. Kaplan–Meier survival analysis and Cox survival analysis were used to detect prognostic factors. R (v3.5.0) and SPSS software (v25.0, IBM) were used to perform statistical analyses.

Results

The median overall survival was 303 days (range 23–2237 days) and the median progression-free survival was 229 days (range 33–1964 days) in patients with sGBM. When assessing the relationship between adjuvant treatment outcome and extent of resection (EOR), the results showed that patients underwent gross total resection can benefit from postoperative radiotherapy and chemotherapy, but not in patients underwent subtotal resection. In addition, we also found that aggressive adjuvant therapy can significantly improve clinical outcomes of IDH1-mutated patients but no significant prognostic value for IDH1-wildtyped patients. The univariate Cox regression analyses demonstrated that EOR, adjuvant therapy, and postoperative Karnofsky Performance Scores were prognostic factors for patients with sGBM, and multivariate COX analysis confirmed that adjuvant therapy and EOR were independent prognostic factors.

Conclusions

For patients with sGBM, aggressive postoperative adjuvant therapy after gross total resection was recommended. However, we did not detect a benefit in IDH1-wildtype patients in our cohort.

Clinicopathological characteristics and survival of spinal cord astrocytomas

Background

Due to their rarity, the clinicopathological characteristics and prognostic factors of spinal cord gliomas are still unclear. Here, we aimed to clarify these issues in a cohort of 108 spinal cord astrocytomas.

Methods

We characterized the clinicopathological characteristics, including 2016 World Health Organization (WHO) grade, age, sex, location, segment length, resection, pre‐ and postsurgery, Modified McCormick Scale (MMS), radio‐ and chemotherapy, and Ki‐67 and H3 K27M mutations, in 108 spinal cord astrocytomas through heatmaps. The Cox regression analysis and Kaplan‐Meier curves were used to study the prognostic value of these clinicopathological features.

Results

There are a total 38 H3 K27M‐mutant tumors, including 31 cases with histological grade II/III tumors. The age of low‐grade astrocytoma patients (WHO grade I/II, n = 54) was significantly younger (27.0 vs 35.5 years, P = .001) than those with high‐grade tumors (WHO grade III/IV, n = 54). All patients underwent surgical resection with neurophysiological monitoring, and the surgery did not result in significant changes in MMS. The presurgery MMS was associated with overall survival in the high‐grade subgroup (P = .008) but not in the low‐grade subgroup (P = .312). While, the high content of resection improved the survival of only patients with low‐grade astrocytomas (P = .016) but not those with high‐grade astrocytomas (P = .475). Both the low‐grade and high‐grade astrocytomas had no obvious benefit from neither adjuvant chemotherapy nor radiotherapy (all P > .05).

Conclusions

We characterized the clinicopathological characteristics and their prognostic values in 108 spinal cord astrocytomas, which could help with evidence‐based management of spinal cord astrocytomas.

A novel Cas9-targeted long-read assay for simultaneous detection of IDH1/2 mutations and clinically relevant MGMT methylation in fresh biopsies of diffuse glioma

Molecular biomarkers provide both diagnostic and prognostic results for patients with diffuse glioma, the most common primary brain tumor in adults. Here, we used a long-read nanopore-based sequencing technique to simultaneously assess IDH mutation status and MGMT methylation level in 4 human cell lines and 8 fresh human brain tumor biopsies. Currently, these biomarkers are assayed separately, and results can take days to weeks. We demonstrated the use of nanopore Cas9-targeted sequencing (nCATS) to identify IDH1 and IDH2 mutations within 36 h and compared this approach against currently used clinical methods. nCATS was also able to simultaneously provide high-resolution evaluation of MGMT methylation levels not only at the promoter region, as with currently used methods, but also at CpGs across the proximal promoter region, the entirety of exon 1, and a portion of intron 1. We compared the methylation levels of all CpGs to MGMT expression in all cell lines and tumors and observed a positive correlation between intron 1 methylation and MGMT expression. Finally, we identified single nucleotide variants in 3 target loci. This pilot study demonstrates the feasibility of using nCATS as a clinical tool for cancer precision medicine.

Systematically characterize the clinical and biological significances of 1p19q genes in 1p/19q non-codeletion glioma

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.

m6A RNA methylation regulators contribute to malignant progression and have clinical prognostic impact in gliomas

N6-methyladenosine (m⁶A) RNA methylation, associated with cancer initiation and progression, is dynamically regulated by the m⁶A RNA methylation regulators (“writers”, “erasers” and “readers”). Here, we demonstrate that most of the thirteen main m⁶A RNA methylation regulators are differentially expressed among gliomas stratified by different clinicopathological features in 904 gliomas. We identified two subgroups of gliomas (RM1/2) by applying consensus clustering to m⁶A RNA methylation regulators. Compared with the RM1 subgroup, the RM2 subgroup correlates with a poorer prognosis, higher WHO grade, and lower frequency of IDH mutation. Moreover, the hallmarks of epithelial-mesenchymal transition and TNFα signaling via NF-κB are also significantly enriched in the RM2 subgroup. This finding indicates that m⁶A RNA methylation regulators are closely associated with glioma malignancy. Based on this finding, we derived a risk signature, using seven m⁶A RNA methylation regulators, that is not only an independent prognostic marker but can also predict the clinicopathological features of gliomas. Moreover, m⁶A regulators are associated with the mesenchymal subtype and TMZ sensitivity in GBM. In conclusion, m⁶A RNA methylation regulators are crucial participants in the malignant progression of gliomas and are potentially useful for prognostic stratification and treatment strategy development.

The molecular characteristics of spinal cord gliomas with or without H3 K27M mutation

Due to the rare incidence of spinal cord astrocytomas, their molecular features remain unclear. Here, we characterized the landscapes of mutations in H3 K27M, isocitrate dehydrogenase 1 (IDH1) R132H, BRAF V600E, and the TERT promoter in 83 diffuse spinal cord astrocytic tumors. Among these samples, thirty-five patients had the H3 K27M mutation; this mutant could be observed in histological grade II (40%), III (40%), and IV (20%) astrocytomas. IDH1 mutations were absent in 58 of 58 cases tested. The BRAF V600E mutation (7/57) was only observed in H3-wildtype astrocytomas, and was associated with a better prognosis in all histological grade II/III astrocytomas. TERT promoter mutations were observed in both H3 K27M-mutant (4/25) and -wildtype (9/33) astrocytomas, and were associated with a poor prognosis in H3-wildtype histological grade II/III astrocytomas. In the 2016 WHO classification of CNS tumors, H3 K27M-mutant diffuse midline gliomas, including spinal cord astrocytomas, are categorized as WHO grade IV. Here, we noticed that the median overall survival of histological grade II/III H3 K27M-mutant cases (n = 28) was significantly longer than that of either the total histological grade IV cases (n = 12) or the H3 K27M-mutant histological grade IV cases (n = 7). We also directly compared H3 K27M-mutant astrocytomas to H3-wildtype astrocytomas of the same histological grade. In histological grade II astrocytomas, compared to H3-wildtype cases (n = 37), H3 K27M-mutant patients (n = 14) had showed a significantly higher Ki-67-positive rate and poorer survival rate. However, no significant differences in these parameters were observed in histological grade III and IV astrocytoma patients. In conclusion, these findings indicate that spinal cord astrocytomas are considerably different from hemispheric and brainstem astrocytomas in terms of their molecular profiles, and that the histological grade cannot be ignored when assessing the prognosis of H3 K27M-mutant spinal cord astrocytomas.

Gene Expression Profiling Stratifies IDH-Wildtype Glioblastoma With Distinct Prognoses

Objectives: In the present study, we aimed to determine the candidate genes that may function as biomarkers to further distinguish patients with isocitrate dehydrogenase (IDH)-wildtype glioblastoma (GBM), which are heterogeneous with respect to clinical outcomes.

Materials and Methods: We selected 41 candidate genes associated with overall survival (OS) using univariate Cox regression from IDH-wildtype GBM patients based on RNA sequencing (RNAseq) expression data from the Chinese Glioma Genome Atlas (CGGA, n = 105) and The Cancer Genome Atlas (TCGA, n = 139) cohorts. Next, a seven-gene-based risk signature was formulated according to Least Absolute Shrinkage and Selection Operator (LASSO) regression algorithm in the CGGA RNAseq database as a training set, while another 525 IDH-wildtype GBM patient TCGA datasets, consisting of RNA sequencing and microarray data, were used for validation. Patient survival in the low- and high-risk groups was calculated using Kaplan-Meier survival curve analysis and the log-rank test. Uni-and multivariate Cox regression analysis was used to assess the prognosis value. Gene oncology (GO) and gene set enrichment analysis (GSEA) were performed for the functional analysis of the seven-gene-based risk signature.

Results: We developed a seven-gene-based signature, which allocated each patient to a risk group (low or high). Patients in the high-risk group had dramatically shorter overall survival than their low-risk counterparts in three independent cohorts. Univariate and multivariate analysis showed that the seven-gene signature remained an independent prognostic factor. Moreover, the seven-gene risk signature exhibited a striking prognostic validity, with AUC of 78.4 and 73.9%, which was higher than for traditional “age” (53.7%, 62.4%) and “GBM sub-type” (57.7%, 52.9%) in the CGGA- and TCGA-RNAseq databases, respectively. Subsequent bioinformatics analysis predicted that the seven-gene signature was involved in the inflammatory response, immune response, cell adhesion, and apoptotic process.

Conclusions: Our findings indicate that the seven-gene signature could be a potential prognostic biomarker. This study refined the current classification system of IDH-wildtype GBM and may provide a novel perspective for the research and individual therapy of IDH-wildtype GBM.

A Novel DNA Methylation-Based Signature Can Predict the Responses of MGMT Promoter Unmethylated Glioblastomas to Temozolomide

Glioblastoma (GBM) is the most malignant glioma, with a median overall survival (OS) of 14–16 months. Temozolomide (TMZ) is the first-line chemotherapy drug for glioma, but whether TMZ should be withheld from patients with GBMs that lack O6-methylguanine-DNA methyltransferase (MGMT) promoter methylation is still under debate. DNA methylation profiling holds great promise for further stratifying the responses of MGMT promoter unmethylated GBMs to TMZ. In this study, we studied 147 TMZ-treated MGMT promoter unmethylated GBM, whose methylation information was obtained from the HumanMethylation27 (HM-27K) BeadChips (n = 107) and the HumanMethylation450 (HM-450K) BeadChips (n = 40) for training and validation, respectively. In the training set, we performed univariate Cox regression and identified that 3,565 CpGs were significantly associated with the OS of the TMZ-treated MGMT promoter unmethylated GBMs. Functional analysis indicated that the genes corresponding to these CpGs were enriched in the biological processes or pathways of mitochondrial translation, cell cycle, and DNA repair. Based on these CpGs, we developed a 31-CpGs methylation signature utilizing the least absolute shrinkage and selection operator (LASSO) Cox regression algorithm. In both training and validation datasets, the signature identified the TMZ-sensitive GBMs in the MGMT promoter unmethylated GBMs, and only the patients in the low-risk group appear to benefit from the TMZ treatment. Furthermore, these identified TMZ-sensitive MGMT promoter unmethylated GBMs have a similar OS when compared with the MGMT promoter methylated GBMs after TMZ treatment in both two datasets. Multivariate Cox regression demonstrated the independent prognostic value of the signature in TMZ-treated MGMT promoter unmethylated GBMs. Moreover, we also noticed that the hallmark of epithelial–mesenchymal transition, ECM related biological processes and pathways were highly enriched in the MGMT unmethylated GBMs with the high-risk score, indicating that enhanced ECM activities could be involved in the TMZ-resistance of GBM. In conclusion, our findings promote our understanding of the roles of DNA methylation in MGMT umethylated GBMs and offer a very promising TMZ-sensitivity predictive signature for these GBMs that could be tested prospectively.

RNA processing genes characterize RNA splicing and further stratify lower-grade glioma

BACKGROUND

Aberrant expression of RNA processing genes may drive the alterative RNA profile in lower-grade gliomas (LGGs). Thus, we aimed to further stratify LGGs based on the expression of RNA processing genes.

METHODS

This study included 446 LGGs from The Cancer Genome Atlas (TCGA, training set) and 171 LGGs from the Chinese Glioma Genome Atlas (CGGA, validation set). The least absolute shrinkage and selection operator (LASSO) Cox regression algorithm was conducted to develop a risk-signature. The receiver operating characteristic (ROC) curves and Kaplan-Meier curves were used to study the prognosis value of the risk-signature.

RESULTS

Among the tested 784 RNA processing genes, 276 were significantly correlated with the OS of LGGs. Further LASSO Cox regression identified a 19-gene risk-signature, whose risk score was also an independently prognosis factor (P<0.0001, multiplex Cox regression) in the validation dataset. The signature had better prognostic value than the traditional factors "age", "grade" and "WHO 2016 classification" for 3- and 5-year survival both two datasets (AUCs > 85%). Importantly, the risk-signature could further stratify the survival of LGGs in specific subgroups of WHO 2016 classification. Furthermore, alternative splicing events for genes such as EGFR and FGFR were found to be associated with the risk score. mRNA expression levels for genes, which participated in cell proliferation and other processes, were significantly correlated to the risk score.

CONCLUSIONS

Our results highlight the role of RNA processing genes for further stratifying the survival of patients with LGGs and provide insight into the alternative splicing events underlying this role.

Pyrosequencing versus methylation-specific PCR for assessment of MGMT methylation in tumor and blood samples of glioblastoma patients

Circulating biomarkers in blood may provide an interesting alternative to risky tissue biopsies in the diagnosis and follow-up of glioblastoma patients. We have assessed MGMT methylation status in blood and tissue samples from unresected glioblastoma patients who had been included in the randomized GENOM-009 trial. Paired blood and tissue samples were assessed by methylation-specific PCR (MSP) and pyrosequencing (PYR). After establishing the minimum PYR cut-off that could yield a significant difference in overall survival, we assessed the sensitivity, specificity, positive predictive value and negative predictive value (NPV) of the analyses. Methylation could be detected in cfDNA by both MSP and PYR but with low concordance with results in tissue. Sensitivity was low for both methods (31% and 38%, respectively), while specificity was higher for MSP in blood than for PYR in plasma (96% vs 76%) and NPV was similar (56 vs 57%). Concordance of results in tissue by MSP and PYR was 84.3% (P < 0.001) and correlated with outcome. We conclude that detection of cfDNA in the blood of glioblastoma patients can be an alternative when tumor tissue is not available but methods for the detection of cfDNA in blood must improve before it can replace analysis in tumor tissue.

Systematically profiling the expression of eIF3 subunits in glioma reveals the expression of eIF3i has prognostic value in IDH-mutant lower grade glioma

Background

Abnormal expression of the eukaryotic initiation factor 3 (eIF3) subunits plays critical roles in tumorigenesis and progression, and also has potential prognostic value in cancers. However, the expression and clinical implications of eIF3 subunits in glioma remain unknown.

Methods

Expression data of eIF3 for patients with gliomas were obtained from the Chinese Glioma Genome Atlas (CGGA) (n = 272) and The Cancer Genome Atlas (TCGA) (n = 595). Cox regression, the receiver operating characteristic (ROC) curves and Kaplan–Meier analysis were used to study the prognostic value. Gene oncology (GO) and gene set enrichment analysis (GSEA) were utilized for functional prediction.

Results

In both the CGGA and TCGA datasets, the expression levels of eIF3d, eIF3e, eIF3f, eIF3h and eIF3l highly were associated with the IDH mutant status of gliomas. The expression of eIF3b, eIF3i, eIF3k and eIF3m was increased with the tumor grade, and was associated with poorer overall survival [All Hazard ratio (HR) > 1 and P < 0.05]. By contrast, the expression of eIF3a and eIF3l was decreased in higher grade gliomas and was associated with better overall survival (Both HR < 1 and P < 0.05). Importantly, the expression of eIF3i (located on chromosome 1p) and eIF3k (Located on chromosome 19q) were the two highest risk factors in both the CGGA [eIF3i HR = 2.068 (1.425–3.000); eIF3k HR = 1.737 (1.166–2.588)] and TCGA [eIF3i HR = 1.841 (1.642–2.064); eIF3k HR = 1.521 (1.340–1.726)] databases. Among eIF3i, eIF3k alone or in combination, the expression of eIF3i was the more robust in stratifying the survival of glioma in various pathological subgroups. The expression of eIF3i was an independent prognostic factor in IDH-mutant lower grade glioma (LGG) and could also predict the 1p/19q codeletion status of IDH-mutant LGG. Finally, GO and GSEA analysis showed that the elevated expression of eIF3i was significantly correlated with the biological processes of cell proliferation, mRNA processing, translation, T cell receptor signaling, NF-κB signaling and others.

Conclusions

Our study reveals the expression alterations during glioma progression, and highlights the prognostic value of eIF3i in IDH-mutant LGG.

Electronic supplementary material

The online version of this article (10.1186/s12935-019-0867-1) contains supplementary material, which is available to authorized users.