ADAMTSL4, a Secreted Glycoprotein, Is a Novel Immune-Related Biomarker for Primary Glioblastoma Multiforme

The Proteomic Landscape of Monocytes in Response to Colorectal Cancer Cells

Colorectal cancer (CRC) involves a complex interaction between tumor cells and immune cells, notably monocytes, leading to immunosuppression. This study explored these interactions using in vitro coculture systems of THP-1 cells and CRC cell lines, employing quantitative proteomics to analyze protein changes in monocytes. Multiple analytical methods were utilized to delineate the altered proteomic landscape, identify key proteins, and their associated functional pathways for comprehensive data analysis. Differentially expressed proteins (DEPs) were selected and validated by cross-referencing them with publicly available TCGA and GEO data sets to explore their potential clinical significance. Our analysis identified 161 up-regulated and 130 down-regulated DEPs. The enrichment results revealed impairments in adhesion and innate immune functions in monocytes, potentially facilitating cancer progression. The down-regulation of FN1, THSB1, and JUN may contribute to these impairments. Furthermore, the overexpression of ADAMTSL4, PRAM1, GPNMB, and NPC2 on monocytes was associated with unfavorable prognostic outcomes in CRC patients, suggesting potential biomarkers or therapeutic targets. This study illustrated the proteomic landscape of monocytes in response to CRC cells, providing clues for future investigations of the crosstalk between cancer cells and monocytes within the tumor microenvironment.

Human serum proteomics reveals a molecular signature after one night of sleep deprivation

Study Objectives

Sleep deprivation is highly prevalent and caused by conditions such as night shift work or illnesses like obstructive sleep apnea. Compromised sleep affects cardiovascular-, immune-, and neuronal systems. Recently, we published human serum proteome changes after a simulated night shift. This pilot proteomic study aimed to further explore changes in human blood serum after 6 hours of sleep deprivation at night.

Methods

Human blood serum samples from eight self-declared healthy females were analyzed using Orbitrap Eclipse mass spectrometry (MS-MS) and high-pressure liquid chromatography. We used a within-participant design, in which the samples were taken after 6 hours of sleep at night and after 6 hours of sleep deprivation the following night. Systems biological databases and bioinformatic software were used to analyze the data and comparative analysis were done with other published sleep-related proteomic datasets.

Results

Out of 494 proteins, 66 were found to be differentially expressed proteins (DEPs) after 6 hours of sleep deprivation. Functional enrichment analysis revealed the associations of these DEPs with several biological functions related to the altered regulation of cellular processes such as platelet degranulation and blood coagulation, as well as associations with different curated gene sets.

Conclusions

This study presents serum proteomic changes after 6 hours of sleep deprivation, supports previous findings showing that short sleep deprivation affects several biological processes, and reveals a molecular signature of proteins related to pathological conditions such as altered coagulation and platelet function, impaired lipid and immune function, and cell proliferation. Data are available via ProteomeXchange with identifier PXD045729. This paper is part of the Genetic and other molecular underpinnings of sleep, sleep disorders, and circadian rhythms including translational approaches Collection.

Loss of p53 epigenetically modulates epithelial to mesenchymal transition in colorectal cancer

Epithelial to Mesenchymal transition (EMT) drives cancer metastasis and is governed by genetic and epigenetic alterations at multiple levels of regulation. It is well established that loss/mutation of p53 confers oncogenic function to cancer cells and promotes metastasis. Though transcription factors like ZEB1, SLUG, SNAIL and TWIST have been implied in EMT signalling, p53 mediated alterations in the epigenetic machinery accompanying EMT are not clearly understood. This work attempts to explore epigenetic signalling during EMT in colorectal cancer (CRC) cells with varying status of p53. Towards this, we have induced EMT using TGFβ on CRC cell lines with wild type, null and mutant p53 and have assayed epigenetic alterations after EMT induction. Transcriptomic profiling of the four CRC cell lines revealed that the loss of p53 confers more mesenchymal phenotype with EMT induction than its mutant counterparts. This was also accompanied by upregulation of epigenetic writer and eraser machinery suggesting an epigenetic signalling cascade triggered by TGFβ signalling in CRC. Significant agonist and antagonistic relationships observed between EMT factor SNAI1 and SNAI2 with epigenetic enzymes KDM6A/6B and the chromatin organiser SATB1 in p53 null CRC cells suggest a crosstalk between epigenetic and EMT factors. The observed epigenetic regulation of EMT factor SNAI1 correlates with poor clinical outcomes in 270 colorectal cancer patients taken from TCGA-COAD. This unique p53 dependent interplay between epigenetic enzymes and EMT factors in CRC cells may be exploited for development of synergistic therapies for CRC patients presenting to the clinic with loss of p53.

Secreted Glycoproteins That Regulate Synaptic Function: the Dispatchers in the Central Nervous System

Glycoproteins are proteins that contain oligosaccharide chains. As widely distributed functional proteins in the body, glycoproteins are essential for cellular development, cellular function maintenance, and intercellular communication. Glycoproteins not only play a role in the cell and the membrane, but they are also secreted in the intercell. These secreted glycoproteins are critical to the central nervous system for neurodevelopment and synaptic transmission. More specifically, secreted glycoproteins play indispensable roles in neurite growth mediation, axon guiding, synaptogenesis, neuronal differentiation, the release of synaptic vesicles, subunit composition of neurotransmitter receptors, and neurotransmitter receptor trafficking among other things. Abnormal expressions of secreted glycoproteins in the central nervous system are associated with abnormal neuron development, impaired synaptic organization/transmission, and neuropsychiatric disorders. This article reviews the secreted glycoproteins that regulate neuronal development and synaptic function in the central nervous system, and the molecular mechanism of these regulations, providing reference for research about synaptic function regulation and related central nervous system diseases.

Assessment tool based on fatty acid metabolic signatures for predicting the prognosis and treatment response in bladder cancer

Background

Fatty acid metabolism (FAM) is closely connected with tumorigenesis as well as disease progression and affects the efficacy of platinum-based drugs. Exploring biomarkers related to FAM in bladder cancer (BLCA) is essential to improve cancer prognosis.

Methods

High-throughput sequencing data from The Cancer Genome Atlas (TCGA) were bioinformatically resolved to identify molecular subtypes of fatty acid metabolic profiles in BLCA using coherent clustering analysis. Based on fatty acid metabolic profile, a prognostic model was created using COX and LASSO COX models. CIBERSORT, Estimation of STromal and Immune cells in MAlignant Tumours using Expression (ESTIMATE), MCP-Count, and single sample gene set enrichment analysis (ssGSEA) were used to assess the differences in tumor microenvironment (TME) among different molecular subtypes, prognostic groups. Kaplan-Meier (K-M) survival curve was plotted to assess patients' prognosis. Receiver operating characteristic curve (ROC) and the clinical prognostic value of prognostic models was evaluated by the Nomogram.

Results

Three molecular subtypes (FAMC1, FAMC2, FAMC3) of fatty acid metabolic patterns were determined. FAMC1 showed significant prognostic advantage with immunoreactivity. Five key prognostic FAMGs were identified and RiskScore was developed. We found that patients with low RiskScore showed significantly better immune microenvironment status, survival and response to immunotherapy. Similarly, both Nomogram and RiskScore demonstrated excellent prognostic value.

Conclusions

In conclusion, our study showed that the RiskScore was closely related to the clinical traits of BLCA patients. The RiskScore may provide essential clinical guidance for predicting prognosis and treatment response in bladder cancer.

Identification and validation of a novel senescence-related biomarker for thyroid cancer to predict the prognosis and immunotherapy

Introduction

Cellular senescence is a hallmark of tumors and has potential for cancer therapy. Cellular senescence of tumor cells plays a role in tumor progression, and patient prognosis is related to the tumor microenvironment (TME). This study aimed to explore the predictive value of senescence-related genes in thyroid cancer (THCA) and their relationship with the TME.

Methods

Senescence-related genes were identified from the Molecular Signatures Database and used to conduct consensus clustering across TCGA-THCA. Differentially expressed genes (DEGs) were identified between the clusters used to perform multivariate Cox regression and least absolute shrinkage and selection operator regression (LASSO) analyses to construct a senescence-related signature. TCGA dataset was randomly divided into training and test datasets to verify the prognostic ability of the signature. Subsequently, the immune cell infiltration pattern, immunotherapy response, and drug sensitivity of the two subtypes were analyzed. Finally, the expression of signature genes was detected across TCGA-THCA and GSE33630 datasets, and further validated by RT-qPCR.

Results

Three senescence clusters were identified based on the expression of 432 senescence-related genes. Then, 23 prognostic DEGs were identified in TCGA dataset. The signature, composed of six genes, showed a significant relationship with survival, immune cell infiltration, clinical characteristics, immune checkpoints, immunotherapy response, and drug sensitivity. Low-risk THCA shows a better prognosis and higher immunotherapy response than high-risk THCA. A nomogram with perfect stability constructed using signature and clinical characteristics can predict the survival of each patient. The validation part demonstrated that ADAMTSL4, DOCK6, FAM111B, and SEMA6B were expressed at higher levels in the tumor tissue, whereas lower expression of MRPS10 and PSMB7 was observed.

Discussion

In conclusion, the senescence-related signature is a promising biomarker for predicting the outcome of THCA and has the potential to guide immunotherapy.

Genes associated with diagnosis and prognosis of Burkitt lymphoma

Burkitt lymphoma (BL) is one of the most aggressive forms of non-Hodgkin's lymphomas that affect children and young adults. The expression of genes and other molecular markers during carcinogenesis can be the basis for diagnosis, prognosis and the design of new and effective drugs for the management of cancers. The aim of this study was to identify genes that can serve as prognostic and therapeutic targets for BL. We analysed RNA-seq data of BL transcriptome sequencing projects in Africa using standard RNA-seq analyses pipeline. We performed pathway enrichment analyses, protein-protein interaction networks, gene co-expression and survival analyses. Gene and pathway enrichment analyses showed that the differentially expressed genes are involved in tube development, signalling receptor binding, viral protein interaction, cell migration, external stimuli response, serine hydrolase activity and PI3K-Akt signalling pathway. Protein-protein interaction network analyses revealed the genes to be highly interconnected, whereas module analyses revealed 25 genes to possess the highest interaction score. Overall survival analyses delineated six genes (ADAMTSL4, SEMA5B, ADAMTS15, THBS2, SPON1 and THBS1) that can serve as biomarkers for prognosis for BL management.

APOL4, a Novel Immune-Related Prognostic Biomarker for Glioma

Glioma is the common, most aggressive and poorest prognostic tumor type in the brain. More and more biomarkers associated with glioma treatment, prognosis, and immunity are being discovered. Here, we aimed to explore the underlying biological functions and prognostic predictive value of Apolipoprotein L4 (APOL4) in glioma. We downloaded the expression data of APOL4 and clinical information from several databases and used R software for preprocessing. The clinical significance of APOL4 in a glioma outcome was explored by the Cox regression analysis and Kaplan–Meier survival analysis. In addition, immune infiltrates and microenvironmental indicators were assessed by CIBERSORT and TIMER. GO and KEGG analyses were used to analyze the potential functions of APOL4 in gliomas. APOL4 expression was increased in glioma specimens compared to normal tissues and correlated dramatically with the WHO grade. A survival analysis showed a shorter overall survival (OS) in glioma patients with APOL4 overexpression, and a Cox regression analysis showed that APOL4 was an independent prognostic factor for the OS of glioma patients. GSEA, GO, and KEGG enrichment analyses showed remarkable enrichment in immune-related pathways. APOL4 expression was positively correlated with immune infiltration (including DC cells, neutrophils, CD8+ T cells, B cells, macrophages, CD4+ T cells, etc.) and microenvironmental parameters (including immune, stromal, and ESTIMATE scores) in gliomas. Glioma patients with a higher expression of APOL4 may be more sensitive to immune checkpoint inhibitors (ICI). In conclusion, these findings suggest that APOL4 is associated with the tumor grade and immune infiltrates; APOL4 may be a new and potential biomarker for therapeutic and prognostic evaluations that may further suggest the therapeutic efficacy of immunotherapy.

Tumor antigens and immune subtypes of glioblastoma: the fundamentals of mRNA vaccine and individualized immunotherapy development

PURPOSE

Glioblastoma (GBM) is the most common primary brain tumor in adults and is notorious for its lethality. Given its limited therapeutic measures and high heterogeneity, the development of new individualized therapies is important. mRNA vaccines have exhibited promising performance in a variety of solid tumors, those designed for glioblastoma (GBM) need further development. The aim of this study is to explore tumor antigens for the development of mRNA vaccines against GBM and to identify potential immune subtypes of GBM to identify the patients suitable for different immunotherapies.

METHODS

RNA-seq data and the clinical information of 143 GBM patients was extracted from the TCGA database; microarray data and the clinical information of 181 GBM patients was obtained from the REMBRANDT cohort. A GBM immunotherapy cohort of 17 patients was obtained from a previous literature. GEPIA2, cBioPortal, and TIMER2 were used to identify the potential tumor antigens. Immune subtypes and gene modules were identified using consensus clustering; immune landscape was constructed using graph-learning-based dimensionality reduction analysis.

RESULTS

Nine potential tumor antigens associated with poor prognosis and infiltration of antigen-presenting cells were identified in GBM: ADAMTSL4, COL6A1, CTSL, CYTH4, EGFLAM, LILRB2, MPZL2, SAA2, and LSP1. Four robust immune subtypes and seven functional gene modules were identified and validated in an independent cohort. Immune subtypes had different cellular and molecular characteristics, with IS1, an immune cold phenotype; IS2, an immune hot and immunosuppressive phenotype; IS3, a relatively immune cold phenotype, second only to IS1; IS4, having a moderate tumor immune microenvironment. Immune landscape revealed the immune distribution of the GBM patients. Additionally, the potential value of immune subtypes for individualized immunotherapy was demonstrated in a GBM immunotherapy cohort.

CONCLUSIONS

ADAMTSL4, COL6A1, CTSL, CYTH4, EGFLAM, LILRB2, MPZL2, SAA2, and LSP1 are the candidate tumor antigens for mRNA vaccine development in GBM, and IS1 GBM patients are best suited for mRNA vaccination, IS2 patients are best suited for immune checkpoint inhibitor. This study provides a theoretical framework for GBM mRNA vaccine development and individualized immunotherapy strategies.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1186/s40537-022-00643-x.

The potential prognostic values of the ADAMTS-like protein family: an integrative pan-cancer analysis

Background

A disintegrin-like and metalloproteinase domain with thrombospondin type 1 motifs (ADAMTS)-like proteins, including ADAMTSL1-6 and papilin, which are part of the mammalian ADAMTS superfamily, appear to be relevant to extracellular matrix function and the regulation of ADAMTS protease activity. Their roles in tumor initiation and progression and regulating the tumor microenvironment (TME) are now recognized.

Methods

In the present study, a comprehensive investigation of the pan-cancer effects of ADAMTSLs and their associations with patient survival, drug responses, and the TME was performed by integrating The Cancer Genome Atlas (TCGA) data and annotated data resources.

Results

The expression of ADAMTSL family members was found to be dysregulated in many cancer types. More importantly, their expression was frequently associated with patients’ overall survival (OS), drug responses, and the TME. ADAMTSL1, ADAMTSL4, and ADAMTSL5 were primarily associated with aggressive phenotypes, while PAPLN was more frequently associated with a favorable prognosis. In a non-small cell lung cancer (NSCLC) cohort, Thrombospondin Type 1 Domain Containing 4 (THSD4) (ADAMTSL6) and Papilin (PAPLN) were associated with immune checkpoint inhibitor (ICI) sensitivity in samples from the Gene Expression Omnibus repository (GSE135222). Twenty and 30 proteins related to THSD4 and PAPLN, respectively, were identified through a proteomic analysis of 18 Chinese lung adenocarcinoma patients.

Conclusions

Our findings extend understandings of the role of the ADAMTSL family in cancers and are a valuable resource on their clinical utility. This article provides insight into the clinical importance of next-generation sequencing technology to identify novel biomarkers for prognosis and investigate therapeutic strategy for clinical benefit.

MXRA5 Is a Novel Immune-Related Biomarker That Predicts Poor Prognosis in Glioma

Background

Glioma is the most common primary intracranial tumor and is associated with poor prognosis. Identifying effective biomarkers for glioma is particularly important. MXRA5, a secreted glycoprotein, is involved in cell adhesion and extracellular matrix remodeling and has been reported to be expressed in many cancers. However, the role and mechanism of action of MXRA5 in gliomas remain unclear. This study was aimed at investigating the role of MXRA5 at the transcriptome level and its clinical prognostic value.

Methods

In this study, RNA microarray data of 301 glioma patients from the Chinese Glioma Genome Atlas (CGGA) were collected as a training cohort and RNA-seq data of 702 glioma samples from The Cancer Genome Atlas (TCGA) were used for validation. We analyzed the clinical and molecular characteristics as well as the prognostic value of MXRA5 in glioma. In addition, the expression level of MXRA was evaluated in 28 glioma tissue samples.

Results

We found that MXRA5 expression was significantly upregulated in high-grade gliomas and IDH wild-type gliomas compared to controls. Receiver operating characteristic (ROC) analysis showed that MXRA5 is a potential marker of the mesenchymal subtype of glioblastoma multiforme (GBM). We found that MXRA5 expression is highly correlated with immune checkpoint molecule expression levels and tumor-associated macrophage infiltration. High MXRA5 expression could be used as an independent indicator of poor prognosis in glioma patients.

Conclusion

Our study suggests that MXRA5 expression is associated with the clinicopathologic features and poor prognosis of gliomas. MXRA5 may play an important role in the immunosuppressive microenvironment of glioma. As a secreted glycoprotein, MXRA5 is a potential circulating biomarker for glioma, deserving further investigation.

Chinese Glioma Genome Atlas (CGGA): A Comprehensive Resource with Functional Genomic Data from Chinese Glioma Patients

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.

DZIP3 is a key factor to stratify IDH1 wild-type lower-grade gliomas

BACKGROUND

Malignant glioma is the most common form of primary malignant brain cancer. Heterogeneity is the hallmark of glioma. DAZ-interacting zinc finger 3 (DZIP3), acts as an RNA-binding RING-type ubiquitin ligase; however, its function in glioma is yet unclear.

RESULTS

The DZIP3 expression was related to the World Health Organization (WHO) grade and isocitrate dehydrogenase 1(IDH1) status, as well as the clinical outcome. Malignant cases exhibit lower DZIP3 expression. DZIP3 was an independent predictive factor of good prognosis in all grade and lower grade gliomas (p < 0.0001). Gene enrichment analysis and immunohistochemistry indicated that DZIP3 affected the biological behavior of glioma through the angiogenesis pathway. Moreover, based on DZIP3 expression, IDH1 wild-type lower-grade gliomas could be divided into two groups with different survival time.

CONCLUSION

In conclusion, the loss of DZIP3 may be involved in the mechanism of angiogenesis in the invasive biological process of glioma. These findings laid an understanding of DZIP3-specific clinical features in glioma.

METHODS

A total of 325 glioma patients from the Chinese Glioma Genome Atlas (CGGA) RNA-seq cohort comprised the training cohort, while 265 patients from the GSE 16011 array cohort formed the validation cohort. The mRNA expression of DZIP3 and clinical characteristics was assessed. DZIP3 protein expression and microvessel density (MVD) were evaluated by immunohistochemistry (IHC).

Analysis of Long Noncoding RNA ZNF667-AS1 as a Potential Biomarker for Diagnosis and Prognosis of Glioma Patients

Objective

Long noncoding RNAs (lncRNAs) have been strongly associated with various types of cancer. The present study aimed at exploring the diagnostic and prognostic value of lncRNA Zinc finger protein 667-antisense RNA 1 (ZNF667-AS1) in glioma patients. Patients and Methods. The expressions of ZNF667-AS1 were detected in 155 glioma tissues and matched normal brain tissue samples by qRT-PCR. The receiver operating characteristic (ROC) curve was performed to estimate the diagnostic value of ZNF667-AS1. The association between the ZNF667-AS1 expression and clinicopathological characteristics was analyzed by the chi-square test. The Kaplan-Meier method was performed to determine the influence of the ZNF667-AS1 expression on the overall survival and disease-free survival of glioma patients. The Cox regression analysis was used to evaluate the effect of independent prognostic factors on survival outcome. Cell proliferation was measured by the respective cell counting Kit-8 (CCK-8) assays.

Results

We observed that ZNF667-AS1 was significantly upregulated in glioma tissues compared to normal tissue samples (p < 0.01). Higher levels of ZNF667-AS1 were positively associated with the WHO grade (p = 0.018) and KPS score (p = 0.008). ROC assays revealed that the high ZNF667-AS1 expression had an AUC value of 0.8541 (95% CI: 0.8148 to 0.8934) for glioma. Survival data revealed that glioma patients in the high ZNF667-AS1 expression group had significantly shorter 5-year overall survival (p = 0.0026) and disease-free survival (p = 0.0005) time than those in the low ZNF667-AS1 expression group. Moreover, multivariate analyses confirmed that the ZNF667-AS1 expression was an independent predictor of the overall survival and disease-free survival for glioma patients. Functionally, we found that knockdown of ZNF667-AS1 suppressed the proliferation of glioma cells.

Conclusions

Our results suggest that ZNF667-AS1 could be used as a potential diagnostic and prognostic biomarker in glioma.

Cancer Cell Acid Adaptation Gene Expression Response Is Correlated to Tumor-Specific Tissue Expression Profiles and Patient Survival

The acidic pH of the tumor microenvironment plays a critical role in driving cancer development toward a more aggressive phenotype, but the underlying mechanisms are unclear. To this end, phenotypic and genotypic changes induced by adaptation of cancer cells to chronic acidosis have been studied. However, the generality of acid adaptation patterns across cell models and their correlation to the molecular phenotypes and aggressiveness of human cancers are essentially unknown. Here, we define an acid adaptation expression response shared across three cancer cell models, dominated by metabolic rewiring, extracellular matrix remodeling, and altered cell cycle regulation and DNA damage response. We find that many genes which are upregulated by acid adaptation are significantly correlated to patient survival, and more generally, that there are clear correlations between acid adaptation expression response and gene expression change between normal and tumor tissues, for a large subset of cancer patients. Our data support the notion that tumor microenvironment acidity is one of the key factors driving the selection of aggressive cancer cells in human patient tumors, yet it also induces a growth-limiting genotype that likely limits cancer cell growth until the cells are released from acidosis, for instance during invasion.

IGFBP3 as an indicator of lymph node metastasis and unfavorable prognosis for papillary thyroid carcinoma

Lymph node metastasis (LNM) is a usual event in papillary thyroid carcinoma (PTC) patients, which usually leads to poor prognosis. However, the molecular mechanisms of LNM remain unclear. Thus, we aimed to screen the possible key genes in the progression of LNM in PTC patients and further validate their roles. The study involved two phases: a discovery phase and a validation one. In the former phase, the candidate genes were screened by using bioinformatics methods. In the latter one, the genes were firstly assessed in a cohort from the cancer genome atlas (TCGA) to evaluate the associations of their expressions with clinical features and the prognostic values, and then, they were assessed at protein levels by using an immunohistochemical assay. Consequently, IGHBP3 was selected as the candidate gene, which might be enriched in several metabolism-related pathways and cancer progression-related pathways. High expressions of IGHBP3 have an association with gender, advanced clinical stages, high T stages, and the presence of LNM. Survival analysis indicated that IGHBP3 may affect the prognosis of PTC patients. The use of a tissue chip confirmed the view that IGHBP3 might play a crucial role in the LNM of PTC. In conclusion, IGHBP3 might be involved in the development of LNM in PTC patients. IGHBP3 over-expression might be a novel indicator and a potential target for PTC therapy.

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.