A new prognostic model for glioblastoma multiforme based on coagulation-related genes.

Glioblastoma multiforme (GBM) is the most aggressive, common, and lethal type of primary brain tumor. Multiple cancers have been associated with abnormalities in the coagulation system that facilitate tumor invasion and metastasis. In GBM, the prognostic value and underlying mechanism of coagulation-related genes (CRGs) have not been explored. RNA sequencing (RNA-seq) and clinical information on GBM were obtained from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA), respectively. Following the identification of differentially expressed CRGs (DECRGs) between GBM and control samples, the survival-related DECRGs were selected via univariate and multivariate Cox regression analyses to establish a prognostic signature. The prognostic performance and clinical utility of the prognostic signature were assessed by the Kaplan-Meier (KM) analysis and receiver operating characteristic (ROC) curve analysis, and a nomogram was constructed. The signature genes-related underlying mechanisms were analyzed according to gene set enrichment analysis (GSEA), Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and single-cell analysis. Finally, the difference in immune cell infiltration, stromal score, immune score, and Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression data (ESTIMATE) score were compared between different risk groups. A 5-gene prognostic signature (PLAUR, GP6, C5AR1, SERPINA5, F2RL2) was established for overall survival (OS) prediction of GBM patients. The predicted efficiency of the prognostic signature was confirmed in TGGA-GBM dataset and validated in the CGGA-GBM dataset, revealing that it could differentiate GBM patients from controls well, and high risk score was accompanied with poor prognosis. Moreover, biological process (BP) and signaling pathway analyses showed that signature genes were mainly enriched in the functions of blood coagulation and tumor invasion and metastasis. Moreover, high-risk patients exhibited higher levels of immune cell infiltration, stromal score, immune score, and ESTIMATE score than that of low-risk patients. An analysis of coagulation-related prognostic signatures was conducted in this study, as well as how signature genes may affect GBM progress, providing information that might provide new ideas for the development of GBM-related molecular targeted therapies.

Zhou M ，Deng Y ，Fu Y ，Liang R ，Liu Y ，Liao Q ... -  《-》

Construction of an immune-related gene signature for the prognosis and diagnosis of glioblastoma multiforme.

Increasing evidence has suggested that inflammation is related to tumorigenesis and tumor progression. However, the roles of immune-related genes in the occurrence, development, and prognosis of glioblastoma multiforme (GBM) remain to be studied. The GBM-related RNA sequencing (RNA-seq), survival, and clinical data were acquired from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), Chinese Glioma Genome Atlas (CGGA), and Gene Expression Omnibus (GEO) databases. Immune-related genes were obtained from the Molecular Signatures Database (MSigDB). Differently expressed immune-related genes (DE-IRGs) between GBM and normal samples were identified. Prognostic genes associated with GBM were selected by Kaplan-Meier survival analysis, Least Absolute Shrinkage and Selection Operator (LASSO)-penalized Cox regression analysis, and multivariate Cox analysis. An immune-related gene signature was developed and validated in TCGA and CGGA databases separately. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to explore biological functions of the signature. The correlation between immune cell infiltration and the signature was analyzed by single-sample gene set enrichment analysis (ssGSEA), and the diagnostic value was investigated. The gene set enrichment analysis (GSEA) was performed to explore the potential function of the signature genes in GBM, and the protein-protein interaction (PPI) network was constructed. Three DE-IRGs [Pentraxin 3 (PTX3), TNFSF9, and bone morphogenetic protein 2 (BMP2)] were used to construct an immune-related gene signature. Receiver operating characteristic (ROC) curves and Cox analyses confirmed that the 3-gene-based prognostic signature was a good independent prognostic factor for GBM patients. We found that the signature was mainly involved in immune-related biological processes and pathways, and multiple immune cells were disordered between the high- and low-risk groups. GSEA suggested that PTX3 and TNFSF9 were mainly correlated with interleukin (IL)-17 signaling pathway, nuclear factor kappa B (NF-κB) signaling pathway, tumor necrosis factor (TNF) signaling pathway, and Toll-like receptor signaling pathway, and the PPI network indicated that they could interact directly or indirectly with inflammatory pathway proteins. Quantitative real-time PCR (qRT-PCR) indicated that the three genes were significantly different between target tissues. The signature with three immune-related genes might be an independent prognostic factor for GBM patients and could be associated with the immune cell infiltration of GBM patients.

Yu Z ，Yang H ，Song K ，Fu P ，Shen J ，Xu M ，Xu H ... -  《Frontiers in Oncology》

A Prognostic Microenvironment-Related Immune Signature via ESTIMATE (PROMISE Model) Predicts Overall Survival of Patients With Glioma.

In the development of immunotherapies in gliomas, the tumor microenvironment (TME) needs to be investigated. We aimed to construct a prognostic microenvironment-related immune signature via ESTIMATE (PROMISE model) for glioma. Stromal score (SS) and immune score (IS) were calculated via ESTIMATE for each glioma sample in the cancer genome atlas (TCGA), and differentially expressed genes (DEGs) were identified between high-score and low-score groups. Prognostic DEGs were selected via univariate Cox regression analysis. Using the lower-grcade glioma (LGG) data set in TCGA, we performed LASSO regression based on the prognostic DEGs and constructed a PROMISE model for glioma. The model was validated with survival analysis and the receiver operating characteristic (ROC) in TCGA glioma data sets (LGG, glioblastoma multiforme [GBM] and LGG+GBM) and Chinese glioma genome atlas (CGGA). A nomogram was developed to predict individual survival chances. Further, we explored the underlying mechanisms using gene set enrichment analysis (GSEA) and Cibersort analysis of tumor-infiltrating immune cells between risk groups as defined by the PROMISE model. We obtained 220 upregulated DEGs and 42 downregulated DEGs in both high-IS and high-SS groups. The Cox regression highlighted 155 prognostic DEGs, out of which we selected 4 genes (CD86, ANXA1, C5AR1, and CD5) to construct a PROMISE model. The model stratifies glioma patients in TCGA as well as in CGGA with distinct survival outcome (P<0.05, Hazard ratio [HR]>1) and acceptable predictive accuracy (AUCs>0.6). With the nomogram, an individualized survival chance could be predicted intuitively with specific age, tumor grade, Isocitrate dehydrogenase (IDH) status, and the PROMISE risk score. ROC showed significant discrimination with the area under curves (AUCs) of 0.917 and 0.817 in TCGA and CGGA, respectively. GSEA between risk groups in both data sets were significantly enriched in multiple immune-related pathways. The Cibersort analysis highlighted four immune cells, i.e., CD 8 T cells, neutrophils, follicular helper T (Tfh) cells, and Natural killer (NK) cells. The PROMISE model can further stratify both LGG and GBM patients with distinct survival outcomes.These findings may help further our understanding of TME in gliomas and shed light on immunotherapies.

Qiu H ，Li Y ，Cheng S ，Li J ，He C ，Li J ... -  《Frontiers in Oncology》

Identification of cuproptosis-related subtypes, construction of a prognosis model, and tumor microenvironment landscape in gastric cancer.

Cuproptosis is a novel identified regulated cell death (RCD), which is correlated with the development, treatment response and prognosis of cancer. However, the potential role of cuproptosis-related genes (CRGs) in the tumor microenvironment (TME) of gastric cancer (GC) remains unknown. Transcriptome profiling, somatic mutation, somatic copy number alteration and clinical data of GC samples were downloaded from the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database to describe the alterations of CRGs from genetic and transcriptional fields. Differential, survival and univariate cox regression analyses of CRGs were carried out to investigate the role of CRGs in GC. Cuproptosis molecular subtypes were identified by using consensus unsupervised clustering analysis based on the expression profiles of CRGs, and further analyzed by GO and KEGG gene set variation analyses (GSVA). Genes in distinct molecular subtypes were also analyzed by GO and KEGG gene enrichment analyses (GSEA). Differentially expressed genes (DEGs) were screened out from distinct molecular subtypes and further analyzed by GO enrichment analysis and univariate cox regression analysis. Consensus clustering analysis of prognostic DEGs was performed to identify genomic subtypes. Next, patients were randomly categorized into the training and testing group at a ratio of 1:1. CRG Risk scoring system was constructed through logistic least absolute shrinkage and selection operator (LASSO) cox regression analysis, univariate and multivariate cox analyses in the training group and validated in the testing and combined groups. Real-time quantitative polymerase chain reaction (RT-qPCR) was used to evaluate the expression of key Risk scoring genes. Sensitivity and specificity of Risk scoring system were examined by using receiver operating characteristic (ROC) curves. pRRophetic package in R was used to investigate the therapeutic effects of drugs in high- and low- risk score group. Finally, the nomogram scoring system was developed to predict patients' survival through incorporating the clinicopathological features and CRG Risk score. Most CRGs were up-regulated in tumor tissues and showed a relatively high mutation frequency. Survival and univariate cox regression analysis revealed that LIAS and FDX1 were significantly associated with GC patients' survival. After consensus unsupervised clustering analysis, GC patients were classified into two cuproptosis molecular subtypes, which were significantly associated with clinical features (gender, age, grade and TNM stage), prognosis, metabolic related pathways and immune cell infiltration in TME of GC. GO enrichment analyses of 84 DEGs, obtained from distinct molecular subtypes, revealed that DEGs primarily enriched in the regulation of metabolism and intracellular/extracellular structure in GC. Univariate cox regression analysis of 84 DEGs further screened out 32 prognostic DEGs. According to the expression profiles of 32 prognostic DEGs, patients were re-classified into two gene subtypes, which were significantly associated with patients' age, grade, T and N stage, and survival of patients. Nest, the Risk score system was constructed with moderate sensitivity and specificity. A high CRG Risk score, characterized by decreased microsatellite instability-high (MSI-H), tumor mutation burden (TMB) and cancer stem cell (CSC) index, and high stromal and immune score in TME, indicated poor survival. Four of five key Risk scoring genes expression were dysregulated in tumor compared with normal samples. Moreover, CRG Risk score was greatly related with sensitivity of multiple drugs. Finally, we established a highly accurate nomogram for promoting the clinical applicability of the CRG Risk scoring system. Our comprehensive analysis of CRGs in GC demonstrated their potential roles in TME, clinicopathological features, and prognosis. These findings may improve our understanding of CRGs in GC and provide new perceptions for doctors to predict prognosis and develop more effective and personalized therapy strategies.

Wang J ，Qin D ，Tao Z ，Wang B ，Xie Y ，Wang Y ，Li B ，Cao J ，Qiao X ，Zhong S ，Hu X ... -  《Frontiers in Immunology》

Identification of Hypoxia Prognostic Signature in Glioblastoma Multiforme Based on Bulk and Single-Cell RNA-Seq.

Ahmed YB ，Ababneh OE ，Al-Khalili AA ，Serhan A ，Hatamleh Z ，Ghammaz O ，Alkhaldi M ，Alomari S ... -  《Cancers》