Identification of a Prognostic Hypoxia-Associated Gene Set in IDH-Mutant Glioma

WITHDRAWN: LonP1 Drives Proneural Mesenchymal Transition in IDH1-R132H Diffuse Glioma

The authors have withdrawn their manuscript owing to massive revision and data validation. Therefore, the authors do not wish this work to be cited as reference for the project. If you have any questions, please contact the corresponding author.

WITHDRAWN: Dual targeting of mitochondrial Lon peptidase 1 and chymotrypsin-like protease by small molecule BT317, as potential therapeutics in malignant astrocytomas

The authors have withdrawn their manuscript owing to massive revision and data validation. Therefore, the authors do not wish this work to be cited as reference for the project. If you have any questions, please contact the corresponding author.

Overexpression of CD99 is associated with tumor adaptiveness and indicates the tumor recurrence and therapeutic responses in gliomas

Glioma undergoes adaptive changes, leading to poor prognosis and resistance to treatment. CD99 influences the migration and invasion of glioma cells and plays an oncogene role. However, whether CD99 can affect the adaptiveness of gliomas is still lacking in research, making its clinical value underestimated. Here, we enrolled our in-house and public multiomics datasets for bioinformatic analysis and conducted immunohistochemistry staining to investigate the role of CD99 in glioma adaptive response and its clinical implications. CD99 is expressed in more adaptative glioma subtypes and cell states. Under hypoxic conditions, CD99 is upregulated in glioma cells and is associated with angiogenesis and metabolic adaptations. Gliomas with over-expressed CD99 also increased the immunosuppressive tumor-associated macrophages. The relevance with tumor adaptiveness of CD99 presented clinical significance. We discovered that CD99 overexpression is associated with short-time recurrence and validated its prognostic value. Additionally, Glioma patients with high expression of CD99 were resistant to chemotherapy and radiotherapy. The CD99 expression was also related to anti-angiogenic and immune checkpoint inhibitor therapy response. Inhibitors of the PI3K-AKT pathway have therapeutic potential against CD99-overexpressing gliomas. Our study identified CD99 as a biomarker characterizing the adaptive response in glioma. Gliomas with high CD99 expression are highly tolerant to stress conditions such as hypoxia and antitumor immunity, making treatment responses dimmer and tumor progression. Therefore, for patients with CD99-overexpressing gliomas, tumor adaptiveness should be fully considered during treatment to avoid drug resistance, and closer clinical monitoring should be carried out to improve the prognosis.

Hypoxia-induced ROS aggravate tumor progression through HIF-1α-SERPINE1 signaling in glioblastoma

Hypoxia, as an important hallmark of the tumor microenvironment, is a major cause of oxidative stress and plays a central role in various malignant tumors, including glioblastoma. Elevated reactive oxygen species (ROS) in a hypoxic microenvironment promote glioblastoma progression; however, the underlying mechanism has not been clarified. Herein, we found that hypoxia promoted ROS production, and the proliferation, migration, and invasion of glioblastoma cells, while this promotion was restrained by ROS scavengers N-acetyl-L-cysteine (NAC) and diphenyleneiodonium chloride (DPI). Hypoxia-induced ROS activated hypoxia-inducible factor-1α (HIF-1α) signaling, which enhanced cell migration and invasion by epithelial-mesenchymal transition (EMT). Furthermore, the induction of serine protease inhibitor family E member 1 (SERPINE1) was ROS-dependent under hypoxia, and HIF-1α mediated SERPINE1 increase induced by ROS via binding to the SERPINE1 promoter region, thereby facilitating glioblastoma migration and invasion. Taken together, our data revealed that hypoxia-induced ROS reinforce the hypoxic adaptation of glioblastoma by driving the HIF-1α-SERPINE1 signaling pathway, and that targeting ROS may be a promising therapeutic strategy for glioblastoma.

The integrative analysis based on super-enhancer related genes for predicting different subtypes and prognosis of patient with lower-grade glioma

Objective: Emerging evidence revealed that super-enhancer plays a crucial role in the transcriptional reprogramming for many cancers. The purpose aimed to explored how the super-enhancer related genes affects the prognosis and tumor immune microenvironment (TIME) of patients with low-grade glioma (LGG). Methods: In this study, the differentially expressed genes (DEGs) between LGG cohorts and normal brain tissue cohort were identified by the comprehensive analysis of the super-enhancer (SE) related genes. Then non-negative matrix factorization was performed to seek the optimal classification based on the DEGs, while investigating prognostic and clinical differences between different subtypes. Subsequently, a prognostic related signature (SERS) was constructed for the comprehensive evaluation in term of individualized prognosis, clinical characteristics, cancer markers, genomic alterations, and immune microenvironment of patients with LGG. Results: Based on the expression profiles of 170 DEGs, we identified three SE subtypes, and the three subtypes showed significant differences in prognostic, clinicopathological features. Then, nine optimal SE-related genes were selected to construct the SERS through the least absolute shrinkage and selection operator Cox regression analysis. Survival analysis showed that SERS had strong and stable predictive ability for the prognosis of LGG patients in the The Cancer Genome Atlas, China Glioma Genome Atlas, and Remdrandt cohorts, respectively. We also found that SERS was highly correlated with clinicopathological features, tumor immune microenvironment, cancer hallmarks, and genomic alterations in LGG patients. In addition, the predictive power of SERS for immune checkpoint inhibitor treatment is also superior. The qRT-PCR results and immunohistochemical results also confirmed the difference in the expression of four key genes in normal cells and tumors, as well as in normal tissues and tumor tissues. Conclusion: The SERS could be suitable to utilize individualized prognosis prediction and immunotherapy options for LGG patients in clinical application.

Development of an Independent Prognostic Signature Based on Three Hypoxia-Related Genes for Breast Cancer

Background. Hypoxia was considered to be a prognostic indicator in a variety of solid tumors. This study aims at identifying the hypoxia-related genes (HRGs) in breast cancer (BC) and the feasibility of HRGs as a prognostic indicator. Methods. We downloaded the mRNA expression data of BC patients from TCGA and GEO databases. The LASSO Cox regression analysis was applied to screen the hub HRGs to establish a prognostic Risk Score. The independence of Risk Score was assessed by multivariate Cox regression analysis. And the immune checkpoint analysis was also performed. In addition, we also detected the expression level of hub HRGs in MCF-10A cells, MCF-7 cells, and SK-BR-3 cells by RT-qPCR. Results. Three HRGs were identified as hub genes with prognostic value in BC, including CA9, PGK1, and SDC1. The Risk Score constructed by these three genes could efficiently distinguish the prognosis of different BC patients and has been shown to be an independent prognostic indicator. In the high-risk group, patients had lower overall survival and poorer prognosis. In addition, the expression levels of five immune checkpoints (PD1, CTLA4, TIGIT, LAG3, and TIM3) in the high-risk group were significantly higher than those in the low-risk group. Moreover, the expression levels of PGK1 and SDC1 in BC cells were significantly increased. Conclusion. In this study, we established an efficiently model based on three optimal HRGs (CA9, PGK1, and SDC1) could clearly distinguish the prognosis of different BC patients.

Molecular markers related to patient outcome in patients with IDH-mutant astrocytomas grade 2 to 4: A systematic review

BACKGROUND

Grading and classification of IDH-mutant astrocytomas has shifted from solely histology towards histology combined with molecular diagnostics. In this systematic review, we give an overview of all currently known clinically relevant molecular markers within IDH-mutant astrocytomas grade 2 to 4.

METHODS

A literature search was performed in five electronic databases for English original papers on patient outcome with respect to a molecular marker as determined by DNA/RNA sequencing, micro-arrays, or DNA methylation profiling in IDH-mutant astrocytomas grade 2 to 4. Papers were included if molecular diagnostics were performed on tumour tissue of at least 15 IDH-mutant astrocytoma patients, and if the investigated molecular markers were not limited to the diagnostic markers MGMT, ATRX, TERT, and/or TP53.

RESULTS

The literature search identified 4508 unique articles, published between August 2012 and December 2021, of which ultimately 44 articles were included. Numerous molecular markers from these papers were significantly correlated to patient outcome. The associations between patient outcome and non-canonical IDH mutations, PI3K mutations, high expression of MSH2, high expression of RAD18, homozygous deletion of CDKN2A/B, amplification of PDGFRA, copy number neutral loss of chromosomal arm 17p, loss of chromosomal arm 19q, the G-CIMP-low DNA methylation cluster, high total CNV, and high tumour mutation burden were confirmed in multiple studies.

CONCLUSIONS

Multiple genetic and epigenetic markers are associated with survival in IDH-mutant astrocytoma patients. Commonly affected are the RB signalling pathway, the RTK-PI3K-mTOR signalling pathway, genomic stability markers, and (epigenetic) gene regulation.

LncRNA-FAM66C Was Identified as a Key Regulator for Modulating Tumor Microenvironment and Hypoxia-Related Pathways in Glioblastoma

Although the role of hypoxia has been greatly explored and unveiled in glioblastoma (GBM), the mechanism of hypoxia-related long non-coding (lnc) RNAs has not been clearly understood. This study aims to reveal the crosstalk among hypoxia-related lncRNAs, tumor microenvironment (TME), and tumorigenesis for GBM. Gene expression profiles of GBM patients were used as a basis for identifying hypoxia-related lncRNAs. Unsupervised consensus clustering was conducted for classifying samples into different molecular subtypes. Gene set enrichment analysis (GSEA) was performed to analyze the enrichment of a series of genes or gene signatures. Three molecular subtypes were constructed based on eight identified hypoxia-related lncRNAs. Oncogenic pathways, such as epithelial mesenchymal transition (EMT), tumor necrosis factor-α (TNF-α) signaling, angiogenesis, hypoxia, P53 signaling, and glycolysis pathways, were significantly enriched in C1 subtype with poor overall survival. C1 subtype showed high immune infiltration and high expression of immune checkpoints. Furthermore, we identified 10 transcription factors (TFs) that were highly correlated with lncRNA-FAM66C. Three key lncRNAs (ADAMTS9-AS2, LINC00968, and LUCAT1) were screened as prognostic biomarkers for GBM. This study shed light on the important role of hypoxia-related lncRNAs for TME modulation and tumorigenesis in GBM. The eight identified hypoxia-related lncRNAs, especially FAM66C may serve as key regulators involving in hypoxia-related pathways.

Regulation of Glucose Transport in Adipocytes by Interleukin-4

Metabolic abnormalities such as obesity, insulin resistance, and type 2 diabetes mellitus are known to be associated with adipose tissue inflammation and impaired secretion of cytokines. Anti-inflammatory cytokine interleukin-4 (IL-4) was found to promote insulin sensitivity, glucose tolerance, and reduce lipid accumulation in vivo through multiple mechanisms, including direct regulation of lipolysis in adipocytes. However, little is known about its role in adipocyte glucose metabolism. This study reveals that IL-4 upregulates glucose uptake in adipocytes without additional activation of the insulin-dependent IRS1 (insulin receptor substrate 1)-Akt (protein kinase B) pathway. Moreover, the main transcription factor STAT6 (signal transducer and activator of transcription 6), regulated by IL-4, was not involved in adipocyte glucose uptake. The proteomic results showed that IL-4 upregulates expression of proteins involved in mitochondrial biogenesis, renewal, and glucose oxidation. Our study provides a new hypothesis, explaining protective effects of IL-4 against metabolic abnormalities through activation of adipocytes glucose utilization and maintenance of mitochondrial function under metabolic overload conditions.

PLOD1 acts as a tumor promoter in glioma via activation of the HSF1 signaling pathway

Procollagen-lysine, 2-oxoglutarate 5-dioxygenase 1 (PLOD1) is a collagen-related lysyl hydroxylase and its prognostic value in glioma patients was verified. However, its biological function in glioma has yet to be fully investigated. The PLOD1 mRNA status and clinical significance in gliomas were assessed via the GEPIA database. Overexpression or targeted depletion of PLOD1 was carried out in the human glioma cell line U87 and verified by western blotting. CCK8 and colony formation assays were implemented to examine the impact of PLOD1 on the proliferative and colony-forming phenotypes of U87 cells. Luciferase reporter assays and HSF1-specific pharmacologic inhibitors (KRIBB11) were employed to determine the regulatory relationship between PLOD1 and heat shock factor 1 (HSF1). High expression of PLOD1 was observed in tissue samples of glioblastoma multiforme (GBM) and brain lower-grade glioma (LGG). GEPIA overall survival further demonstrated that both GBM and LGG patients with high PLOD1 displayed worse clinical outcomes compared with those with low PLOD1. Overexpression and targeted depletion of PLOD1 enhanced and suppressed U87 cell proliferation and colony formation, respectively. Luciferase reporter assays showed that PLOD1 significantly enhanced the transcriptional activity of HSF1 in HEK293T cells. PLOD1 deficiency in U87 cells inhibited HSF1-induced survivin accumulation, whereas KRIBB11 also blocked the PLOD1-overexpressing induced survivin expression. An inhibitor of HSF1 signaling events abolished the increased clonogenic potential caused by PLOD1 overexpression in U87 cells. High expression of PLOD1 can increase the proliferation and colony formation of U87 cells by activating the HSF1 signaling pathway. This study suggested PLOD1/HSF1 as an effective therapeutic target for gliomas.

Identifying a hypoxia related score to predict the prognosis of bladder cancer: a study with The Cancer Genome Atlas (TCGA) database

Background

Recurrence is common in bladder cancer, with a hypoxic tumor microenvironment (TME) playing a role in genetic instability and prognosis of bladder cancer. However, we still lack practical hypoxia related model for predicting the prognosis of bladder cancer. In this study, we identified new prognosis-related hypoxia genes and established a new hypoxia score related signature.

Methods

The Gene Set Variation Analysis (GSVA) algorithm was utilized to calculate the hypoxia score of bladder cancer cases found on the The Cancer Genome Atlas (TCGA) database on the gene expression profiles. The cases were first divided into low- and high-hypoxia score groups and then differentially expressed genes (DEGs) expression analysis was conducted. Hypoxia-related genes were identified using weighted gene co-expression network analysis (WGCNA). We then conducted a protein-protein interaction (PPI) network and carried out functional enrichment analysis of the genes that overlapped between DEGs and hypoxia-related genes. LASSO Cox regression analysis was used to establish a hypoxia-related prognostic signature, which was validated using the GSE69795 dataset downloaded from GEO database.

Results

Results from Kaplan-Meier analysis showed that patients with a high hypoxia score had significantly poor overall survival compared to patients with low hypoxia score. We selected 270 DEGs between low- and high-hypoxia score groups, while WGCNA analysis identified 1,313 genes as hypoxia-related genes. A total of 170 genes overlapped between DEGs and hypoxia-related genes. LASSO algorithms identified 29 genes associated with bladder cancer prognosis, which were used to construct a novel 29-gene signature model. The prognostic risk model performed well, since the receiver operating characteristic (ROC) curve showed an accuracy of 0.802 (95% CI: 0.759–0.844), and Cox proportional hazards regression analysis proved the model an independent predictor with hazard ratio (HR) =1.789 (95% CI: 1.585–2.019) (P<0.001). The low-risk score patients had remarkably longer overall survival than patients with a higher score (survival rate 71.06% vs. 23.66%) in the The Cancer Genome Atlas (TCGA) cohort (P<0.0001) and in the dataset GSE69795 (P=0.0079).

Conclusions

We established a novel 29-gene hypoxia-related signature model to predict the prognosis of bladder cancer cases. This model and identified hypoxia-related genes may further been used as biomarkers, assisting the evaluation of prognosis of bladder cancer cases and decision making in clinical practice.

Metabolic adaptations in cancers expressing isocitrate dehydrogenase mutations

The most frequently mutated metabolic genes in human cancer are those encoding the enzymes isocitrate dehydrogenase 1 (IDH1) and IDH2; these mutations have so far been identified in more than 20 tumor types. Since IDH mutations were first reported in glioma over a decade ago, extensive research has revealed their association with altered cellular processes. Mutations in IDH lead to a change in enzyme function, enabling efficient conversion of 2-oxoglutarate to R-2-hydroxyglutarate (R-2-HG). It is proposed that elevated cellular R-2-HG inhibits enzymes that regulate transcription and metabolism, subsequently affecting nuclear, cytoplasmic, and mitochondrial biochemistry. The significance of these biochemical changes for tumorigenesis and potential for therapeutic exploitation remains unclear. Here we comprehensively review reported direct and indirect metabolic changes linked to IDH mutations and discuss their clinical significance. We also review the metabolic effects of first-generation mutant IDH inhibitors and highlight the potential for combination treatment strategies and new metabolic targets.

Integrated Analysis of Weighted Gene Coexpression Network Analysis Identifying Six Genes as Novel Biomarkers for Alzheimer's Disease

BACKGROUND

Alzheimer's disease (AD) is a chronic progressive neurodegenerative disease; however, there are no comprehensive therapeutic interventions. Therefore, this study is aimed at identifying novel molecular targets that may improve the diagnosis and treatment of patients with AD.

METHODS

In our study, GSE5281 microarray dataset from the GEO database was collected and screened for differential expression analysis. Genes with a P value of <0.05 and ∣log2FoldChange | >0.5 were considered differentially expressed genes (DEGs). We further profiled and identified AD-related coexpression genes using weighted gene coexpression network analysis (WGCNA). Functional enrichment analysis was performed to determine the characteristics and pathways of the key modules. We constructed an AD-related model based on hub genes by logistic regression and least absolute shrinkage and selection operator (LASSO) analyses, which was also verified by the receiver operating characteristic (ROC) curve.

RESULTS

In total, 4674 DEGs were identified. Nine distinct coexpression modules were identified via WGCNA; among these modules, the blue module showed the highest positive correlation with AD (r = 0.64, P = 3e - 20), and it was visualized by establishing a protein-protein interaction network. Moreover, this module was particularly enriched in "pathways of neurodegeneration-multiple diseases," "Alzheimer disease," "oxidative phosphorylation," and "proteasome." Sixteen genes were identified as hub genes and further submitted to a LASSO regression model, and six genes (EIF3H, RAD51C, FAM162A, BLVRA, ATP6V1H, and BRAF) were identified based on the model index. Additionally, we assessed the accuracy of the LASSO model by plotting an ROC curve (AUC = 0.940).

CONCLUSIONS

Using the WGCNA and LASSO models, our findings provide a better understanding of the role of biomarkers EIF3H, RAD51C, FAM162A, BLVRA, ATP6V1H, and BRAF and provide a basis for further studies on AD progression.

Integrated genomic and transcriptomic analysis suggests KRT18 mutation and MTAP are key genetic alterations related to the prognosis between astrocytoma and glioblastoma

Background

Astrocytoma and glioblastoma (GBM) are the two main subtypes of glioma, with the 2016 World Health Organization Classification of Tumors of the Central Nervous System (CNS WHO) classifying them into different grades. GBM is the most malignant among all CNS tumors with a 5-year survival rate of less than 5%. Although the prognosis of patients with astrocytoma is better than that of GBM in general, patients with anaplastic astrocytoma (AA) and isocitrate dehydrogenase (IDH) wild type have a similar prognosis as GBM and entail a high risk of progression. Exploring the molecular driving force behind the malignant phenotype of astrocytoma and GBM will help explain the diversity of glioma and discover new drug targets.

Methods

We enrolled 12 patients with astrocytoma and 12 patients with GBM and performed whole-exome sequencing (WES) and RNA sequencing analysis on tumor samples from the patients.

Results

We found that the somatic mutation of KRT18, which is associated with cell apoptosis and adhesion by interacting with receptor 1-associated protein (TRADD) and pinin, was significantly enriched in astrocytoma, but rare in GBM. Copy number loss of MTAP, which is closely related to a poor prognosis of glioma, was found to be significantly enriched in GBM. In addition, different somatic copy number alteration (SCNA), gene expression, and immune cell infiltration patterns between astrocytoma and GBM were found.

Conclusions

This study revealed the distinct characteristics of astrocytoma and GBM at the DNA and RNA level. Somatic mutation of KRT18 and copy number loss of MTAP, two key genetic alterative genes in astrocytoma and GBM, have the potential to become therapeutic targets in glioma.

Integrated Gene Expression and Methylation Analyses Identify DLL3 as a Biomarker for Prognosis of Malignant Glioma

Glioma is one of the most common neurological malignancies worldwide. Delta-like ligand 3 (DLL3), an inhibitory ligand-driven activation of the Notch pathway, has been shown to be significantly associated with overall survival in patients with glioma. Therefore, the purpose of this study was to determine whether DLL3 as a biomarker in glioma is associated with patients' clinicopathological features and prognosis. We identified differences in transcriptome and promoter methylation in the Chinese Glioma Genome Atlas (CGGA) in patients with malignant glioma with shorter (less than 1 year) and longer (greater than 3 years) survival time. Further analysis of The Cancer Genome Atlas (TCGA) revealed that four genes (DLL3, TSPAN15, RTN1, PAK7) are highly associated with patient prognosis and play an indispensable role in evolution. We chose the expression level of DLL3 in glioma patients for our study. Patients were divided into groups with low and high expression of DLL3 according to the cutoff values obtained, and Kaplan-Meier and Cox analysis were used to examine the correlation between DLL3 gene expression and patient survival. We then performed a gene set enrichment analysis (GSEA) to identify significantly enriched signaling pathways. Our results confirmed that the overall survival of patients with low DLL3 expression was significantly shorter than that of patients with high DLL3 expression. GSEA showed that the signaling pathways of the immune process and immune response, among others, were enhanced with the DLL3 low-expression phenotype. Collectively, our findings signify that DLL3 is a potent prognostic factor for glioma, which can provide a viable approach for glioma prognostic assessment and valuable insights for anti-tumor immune-targeted therapies.

Hypoxia-induced PLOD1 overexpression contributes to the malignant phenotype of glioblastoma via NF-κB signaling

Procollagen lysyl hydroxylase 1 (PLOD1) is highly expressed in malignant tumors such as esophageal squamous cell carcinoma, gastric cancer, and colorectal cancer. Bioinformatics analysis revealed that PLOD1 is associated with the progression of GBM, particularly the most malignant mesenchymal subtype (MES). Moreover, in the TCGA and CGGA datasets, the mean survival time of patients with high PLOD1 expression was significantly shorter than that of patients with low expression. The clinical samples confirmed this result. Therefore, we aimed to investigate the effect of PLOD1 on the development of mesenchymal GBM in vitro and in vivo and its possible mechanisms. Molecular experiments were conducted on the patient-derived glioma stem cells and found that PLOD1 expressed higher in tumor tissues and cancer cell lines of patients with GBM, especially in the MES. PLOD1 also enhanced tumor viability, proliferation, migration, and promoted MES transition while inhibited apoptosis. Tumor xenograft results also indicated that PLOD1 overexpression significantly promotes malignant behavior of tumors. Mechanistically, bioinformatics analysis further revealed that PLOD1 expression was closely associated with the NF-κB signaling pathway. Besides, we also found that hypoxic environments also enhanced the tumor-promoting effects of PLOD1. In conclusion, overexpression of PLOD1 may be an important factor in the enhanced invasiveness and MES transition of GBM. Thus, PLOD1 is a potential treatment target for mesenchymal GBM or even all GBM.

Comprehensive analysis of PLOD family members in low-grade gliomas using bioinformatics methods

Low-grade gliomas (LGGs) is a primary invasive brain tumor that grows slowly but is incurable and eventually develops into high malignant glioma. Novel biomarkers for the tumorigenesis and lifetime of LGG are critically demanded to be investigated. In this study, the expression levels of procollagen-lysine, 2-oxoglutarate 5-dioxygenases (PLODs) were analyzed by ONCOMINE, HPA and GEPIA. The GEPIA online platform was applied to evaluate the interrelation between PLODs and survival index in LGG. Furthermore, functions of PLODs and co-expression genes were inspected by the DAVID. Moreover, we used TIMER, cBioportal, GeneMINIA and NetworkAnalyst analysis to reveal the mechanism of PLODs in LGG. We found that expression levels of each PLOD family members were up-regulated in patients with LGG. Higher expression of PLODs was closely related to shorter disease-free survival (DFS) and overall survival (OS). The findings showed that LGG cases with or without alterations were significantly correlated with the OS and DFS. The mechanism of PLODs in LGG may be involved in response to hypoxia, oxidoreductase activity, Lysine degradation and immune cell infiltration. In general, this research has investigated the values of PLODs in LGG, which could serve as biomarkers for diagnosis, prognosis and potential therapeutic targets of LGG patients.

Expression and Prognostic Role of PLOD1 in Malignant Glioma

Background

Malignant glioma is rarely curable, and factors that influence the prognosis of glioma patients are not fully understood. Lysyl hydroxylases such as PLOD1 promote the cross-linking in extracellular matrix (ECM) molecules, which contribute to ECM structural stability and maturation. However, the expression and prognostic role of PLOD1 in malignant glioma remained to be determined.

Methods

The expression of PLOD1 was evaluated by immunohistochemistry in 72 malignant glioma patients from Shenzhen People's hospital. The mRNA expression profiles and clinical information of malignant glioma patients were obtained from public databases, including TCGA, CGGA, Rembrandt, and Gravendeel. The correlation between gene expression and tumor grade, and IDH1/2 status and 1p19q status were evaluated. The association between gene expression and overall survival of malignant glioma patients was examined using Kaplan-Meier survival analysis. GO and KEGG pathways were analyzed by Metascape. Transwell invasion assays were performed to determine the effect of PLOD1 on migration and invasion of glioma cells in vitro.

Results

PLOD1 expression was significantly elevated in malignant glioma tissues compared with non-tumor brain tissues. Besides, elevated levels of PLOD1 were significantly correlated with high tumor grade, wildtype IDH1/2 status, and 1p19q non-codel in all the four public datasets and in-house cohort. Malignant glioma patients with high PLOD1 expression had better overall survival compared to those with low PLOD1 expression. More importantly, patients with IDH1/2 mutations, 1p19q codeletions, and PLOD1 overexpression had the best overall survival. GO enrichment pathway analysis indicated that PLOD1 participates in regulating the extracellular matrix. Transwell invasion assay, which revealed that inhibiting PLOD1 reduced cell invasion in both U87 and U251 cells.

Conclusion

PLOD1 serves as a potential prognostic marker and therapeutic target for malignant glioma.

Cystathionine-γ-lyase drives antioxidant defense in cysteine-restricted IDH1-mutant astrocytomas

BACKGROUND

Mutations in isocitrate dehydrogenase 1 or 2 (IDH1/2) define glioma subtypes and are considered primary events in gliomagenesis, impacting tumor epigenetics and metabolism. IDH enzyme activity is crucial for the generation of reducing potential in normal cells, yet the impact of the mutation on the cellular antioxidant system in glioma is not understood. The aim of this study was to determine how glutathione (GSH), the main antioxidant in the brain, is maintained in IDH1-mutant gliomas, despite an altered NADPH/NADP balance.

METHODS

Proteomics, metabolomics, metabolic tracer studies, genetic silencing, and drug targeting approaches in vitro and in vivo were applied. Analyses were done in clinical specimen of different glioma subtypes, in glioma patient-derived cell lines carrying the endogenous IDH1 mutation and corresponding orthotopic xenografts in mice.

RESULTS

We find that cystathionine-γ-lyase (CSE), the enzyme responsible for cysteine production upstream of GSH biosynthesis, is specifically upregulated in IDH1-mutant astrocytomas. CSE inhibition sensitized these cells to cysteine depletion, an effect not observed in IDH1 wild-type gliomas. This correlated with an increase in reactive oxygen species and reduced GSH synthesis. Propargylglycine (PAG), a brain-penetrant drug specifically targeting CSE, led to delayed tumor growth in mice.

CONCLUSIONS

We show that IDH1-mutant astrocytic gliomas critically rely on NADPH-independent de novo GSH synthesis via CSE to maintain the antioxidant defense, which highlights a novel metabolic vulnerability that may be therapeutically exploited.

Retracted: Identification and validation of hypoxic markers of glioblastoma multiforme in methylation, somatic CNV characteristics based on large-scale samples

Retraction: Wang J, Xu S, Lv W, Shan A, Shi F, Xu J, Mei S, and Yang Y. Identification and validation of hypoxic markers of glioblastoma multiforme in methylation, somatic CNV characteristics based on large-scale samples. J Biochem Mol Toxicol. 2020; e22585. https://doi.org/10.1002/jbt.22585 The above article, published online on 24 July 2020 in Wiley Online Library (wileyonlinelibrary.com) has been retracted by agreement between the authors, the journal's Editor-in-Chief, Hari K. Bhat, and Wiley Periodicals, LLC. The retraction has been agreed due to the authors raising concerns about the accuracy of their data following their article's acceptance. The authors had agreed with the Editor-in-Chief to withdraw their manuscript before publication. However, unfortunately, the Version of Record was mistakenly published online on the journal's Early View service before it could be withdrawn. The Publisher apologizes to the authors for this error.

A novel hypoxic tumor microenvironment signature for predicting the survival, progression, immune responsiveness and chemoresistance of glioblastoma: a multi-omic study

The hypoxic tumor microenvironment (TME) was reported to promote the aggressive phenotype, progression, recurrence, and chemoresistance of glioblastoma (GBM). We developed and validated a hypoxia gene signature for individualized prognostic prediction in GBM patients. In total, 259 GBM-specific hypoxia-related genes (HRGs) were obtained in hypoxic cultured GBM cells compared with normoxic cells. By applying the k-means algorithm, TCGA GBM patients were divided into two subgroups, and the patients in Cluster 1 exhibited high HRG expression patterns, older age, and poor prognosis, which was validated in the CGGA cohort. Cox regression analyses were performed to generate an HRG-based risk score model consisting of five HRGs, which could reliably discriminate the overall survival (OS) and progression-free survival (PFS) of high- and low-risk patients in both the TCGA training and CGGA validation cohorts. Then, nomograms with the hypoxia signature for OS and PFS prediction were constructed for individualized survival prediction, better treatment decision-making, and follow-up scheduling. Finally, functional enrichment, immune infiltration, immunotherapy response prediction and chemotherapy resistance analyses demonstrated the vital roles of the hypoxic TME in the development, progression, multitherpy resistance of GBM. The hypoxia gene signature could serve as a promising prognostic predictor and potential therapeutic target to combat chemoresistant GBM.

Large-Scale Drug Screening in Patient-Derived IDHmut Glioma Stem Cells Identifies Several Efficient Drugs among FDA-Approved Antineoplastic Agents

The discovery of the isocitrate dehydrogenase (IDH) mutation in glioma led to a paradigm shift on how we see glioma biology. Difficulties in cultivating IDH mutant glioma stem cells (IDH^mut GSCs) resulted in a paucity of preclinical models in IDH^mut glioma, limiting the discovery of new effective chemotherapeutic agents. To fill this gap, we used six recently developed patient-derived IDH^mut GSC lines and performed a large-scale drug screening with 147 Food and Drug Administration (FDA)-approved anticancer drugs. GSCs were subjected to the test compounds for 72 h in concentrations ranging from 0.0001 to 1 µM. Cell viability was assessed by CellTiterGlo and the induction of apoptosis by flow cytometry with Annexin V/propidium iodide staining. The initial screen was performed with two IDH^mut GSC lines and identified seven drugs (bortezomib, carfilzomib, daunorubicin, doxorubicin, epirubicin, omacetaxine, plicamycin) with a substantial antiproliferative activity, as reflected by half maximal inhibitory concentrations (IC₅₀) below 1 µM and maximum inhibitory effects (E_max) below 25%. These findings were validated in an additional four IDH^mut GSC lines. The candidate drugs, of which plicamycin and omacetaxine are known to cross the blood brain barrier, were used for subsequent cell death analyses. A significant induction of apoptosis was observed at IC₅₀ values of the respective drugs. In summary, we were able to identify seven FDA-approved drugs that should be further taken into clinical investigations for the treatment of IDH^mut gliomas.

A novel molecular link between HOXA9 and WNT6 in glioblastoma identifies a subgroup of patients with particular poor prognosis

Despite much effort to improve treatments, patients with malignant glioma still present a very poor prognosis that has not changed significantly in the last decades. In this context, it is crucial to better understand glioma pathogenesis to identify new molecular prognostic subgroups and therapeutic targets. WNT6 was recently identified as a new oncogenic molecule in glioblastoma (GBM), with prognostic value in patients, but the mechanisms underlying WNT6 aberrant expression in glioma are still unknown. WNT6 was overexpressed in a subset of gliomas independently of IDH mutations, 1p/19q codeletion status, and WNT6 gene copy number. Interestingly, WNT6 expression is associated with the DNA methylation levels of particular CpG regions at both the WNT6 promoter and the gene body in glioma patient samples. HOXA9, a transcription factor previously associated with poorer clinical outcome in GBM, was identified as a novel transcriptional regulator of WNT6, activating the WNT/β‐catenin pathway in vitro and in vivo. In various cohorts of glioma patients, mRNA levels of WNT6 and HOXA9 were significantly correlated, extending our in vitro and in vivo findings into the clinical setting. Interestingly, this novel molecular link between WNT6 and HOXA9 was not limited to glioma, as they were co‐expressed also in patients with other tumor types. Clinically, WNT6 was a prognostic biomarker of shorter survival in GBM, independently of HOXA9 expression. Concomitant high expression of both WNT6 and HOXA9 identified a subgroup of patients with particularly dismal survival. These findings describe novel WNT6 regulatory mechanisms in GBM, establishing particular DNA methylation patterns and HOXA9 as critical regulators of WNT6 expression in glioma. This HOXA9‐WNT6 molecular link supports WNT signaling in GBM cells and is a powerful prognostic biomarker, highlighting the clinical relevance of this axis in patients. Novel therapies targeting WNT6‐HOXA9 signaling may thus be useful for this deadly disease.

Formation of new lymphatic vessels in glioma: An immunohistochemical analysis

We investigated the distribution and formation of new lymphatic vessels in gliomas. Specimens from seven glioma cases were analyzed by immunohistochemical staining for CD34, lymphatic endothelial hyaluronic acid receptor 1 (LYVE‐1), prospero‐related homeobox 1 (Prox1), nestin, and hypoxia‐inducible factor 1α (HIF‐1α). Three types of vessels were observed in glioma specimens: LYVE‐1⁺ lymphatic vessels, CD34⁺ blood vessels, and LYVE‐1⁺/CD34⁺ blood vessels. Prox1⁺/LYVE‐1⁺ cells were distributed in some lymphatic vessels as well as among vascular endothelial cells and glioma cells. Nestin⁺ cells were scattered throughout the gliomas, and some lymphatic cells also expressed nestin. HIF‐1α⁺ Prox1⁺ cells were widely distributed within the glioma specimens. The present immunohistochemical analysis revealed upregulation of Prox1 and HIF‐1α in some glioma tissues as well as the differentiation of nestin⁺ tumor stem cells into LYVE‐1⁺ lymphatic vessels.

Integrating tumor hypoxic stress in novel and more adaptable strategies for cancer immunotherapy

Immunotherapy is poised to become an increasingly utilized therapy in the treatment of cancer. However, several abnormalities in the tumor microenvironment (TME) that can thwart the efficacy of immunotherapies have been established. Microenvironmental hypoxia is a determining factor in shaping aggressiveness, metastatic potential and treatment resistance of solid tumors. The characterization of this phenomenon could prove beneficial for determining a patient's treatment path and for the introduction of novel targetable factors that can enhance therapeutic outcome. Indeed, the ablation of hypoxia has the potential to sensitize tumors to immunotherapy by metabolically remodeling their microenvironment. In this review, we discuss the intrinsic contributions of hypoxia to cellular plasticity, heterogeneity, stemness and genetic instability in the context of immune escape. In addition, we will shed light on how managing hypoxia can ameliorate response to immunotherapy and how integrating hypoxia gene signatures could play a role in this pursuit.

Increased RLIP76 expression in IDH1 wild‑type glioblastoma multiforme is associated with worse prognosis

Mutation of the isocitrate dehydrogenase (IDH) gene is regarded a novel indicator for the prognosis of patients with glioma. However, the role of the IDH1 gene mutations in carcinogenesis and the mechanisms underlying their function in glioblastoma multiforme (GBM) remain unknown. The present study aimed to determine whether the association of RLIP76 with the different IDH1 mutational status could serve as a putative biomarker for improving disease prognosis. Quantitative PCR, western blotting and immunohistochemical staining assays were used to investigate the expression levels of RLIP76 in 124 patients with GBM with different IDH1 mutational status. In addition, the association between RLIP76 expression, IDH1 mutational status and clinicopathological characteristics was investigated. The effects of RLIP76 expression and IDH1 mutational status on cell proliferation, cell apoptosis, and cell signaling were examined by Cell Counting Kit-8, flow cytometry and western blot assays, respectively. The data demonstrated that IDH1 wild-type (IDH1^Wt) patients with low RLIP76 expression exhibited improved overall and progression-free survival. This effect was not observed in patients with IDH1 mutant (IDH1^Mut) GBM. In vitro assays demonstrated that knockdown of IDH1 or overexpression of the IDH1 R132H mutation suppressed cell proliferation and promoted cell apoptosis in U87 glioma cells. Mechanistic studies further indicated that although the IDH1 R132H mutant phenotype exhibited similar antitumor effects on GBM cells as those observed with the IDH1 knockdown, it acted via a different mechanism with regard to the regulation of the apoptosis signaling pathway. IDH1 R132H mutant cells promoted p53-induced apoptosis, while the IDH1 knockdown inhibited the RLIP76-dependent apoptotic pathway in glioma cells. The findings of the present study provided insight to the contribution of IDH1 mutation in the development of GBM and indicated that RLIP76 may be considered as a prognostic biomarker of IDH1^Wt GBM.

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.