Profiles of immune-related genes and immune cell infiltration in the tumor microenvironment of diffuse lower-grade gliomas

The tumor microenvironment is highly correlated with tumor occurrence, progress, and prognosis. We aimed to investigate the immune-related gene (IRG) expression and immune infiltration pattern in the tumor microenvironment of lower-grade glioma (LGG). We employed the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm to calculate immune and stromal scores and identify prognostic IRG based on The Cancer Genome Atlas data set. The potential molecular functions of these genes were explored with the help of functional enrichment analysis and the protein–protein interaction network. Remarkably, three cohorts that were downloaded from the Chinese Glioma Genome Atlas database were analyzed to further verify the prognostic values of these genes. Moreover, the Tumor IMmune Estimation Resource (TIMER) algorithm was used to estimate the abundance of infiltrating immune cells and explore the immune infiltration pattern in LGG. And unsupervised cluster analysis determined three clusters of the immune infiltration pattern and indicated that CD8⁺ T cells and macrophages were significantly associated with LGG outcomes. Altogether, our study identified a list of prognostic IRGs and provided a perspective to explore the immune infiltration pattern in LGG.     

Ferroptosis-related gene signature correlates with the tumor immune features and predicts the prognosis of glioma patients

Background: Glioma is a malignant intracranial tumor and the most fatal cancer. The role of ferroptosis in the clinical progression of gliomas is unclear.Method: Univariate and least absolute shrinkage and selection operator (Lasso) Cox regression methods were used to develop a ferroptosis-related signature (FRSig) using a cohort of glioma patients from the Chinese Glioma Genome Atlas (CGGA), and was validated using an independent cohort of glioma patients from The Cancer Genome Atlas (TCGA). A single-sample gene set enrichment analysis (ssGSEA) was used to calculate levels of the immune infiltration. Multivariate Cox regression was used to determine the independent prognostic role of clinicopathological factors and to establish a nomogram model for clinical application.Results: We analyzed the correlations between the clinicopathological features and ferroptosis-related gene (FRG) expression and established an FRSig to calculate the risk score for individual glioma patients. Patients were stratified into two subgroups with distinct clinical outcomes. Immune cell infiltration in the glioma microenvironment and immune-related indexes were identified that significantly correlated with the FRSig, the tumor mutation burden (TMB), copy number alteration (CNA), and immune checkpoint expression was also significantly positively correlated with the FRSig score. Ultimately, an FRSig-based nomogram model was constructed using the independent prognostic factors age, World Health Organization (WHO) grade, and FRSig score.Conclusion: We established the FRSig to assess the prognosis of glioma patients. The FRSig also represented the glioma microenvironment status. Our FRSig will contribute to improve patient management and individualized therapy by offering a molecular biomarker signature for precise treatment.     

Mutational profiling of low-grade gliomas identifies prognosis and immunotherapy-related biomarkers and tumour immune microenvironment characteristics

Low-grade glioma (LGG) is a heterogeneous tumour with the median survival rate less than 10 years. Therefore, it is urgent to develop efficient immunotherapy strategies of LGG. In this study, we analysed mutation profiles based on the data of 510 LGG patients from the Cancer Genome Atlas (TCGA) database and investigated the prognostic value of mutated genes and evaluate their immune infiltration. Tumor Immune Dysfunction and Exclusion (TIDE) algorithm was used to indicate the characteristics of gliomas that respond to immune checkpoint blockade (ICB) therapy. Univariate and multivariate cox regression analysis was performed to identify indicators to construct the nomogram model. 485 (95.47%) of 508 LGG samples showed gene mutation, and 9 mutated genes were significantly related to overall survival (OS), among which 6 mutated genes were significantly correlated with OS between mutation and wildtypes. Immune infiltration and immune score analyses revealed that these six mutated genes were significantly associated with tumour immune microenvironment in LGG. The response of LGG with different characteristics to ICB was evaluated by TIDE algorithm. Finally, CIC gene was screened through both univariate and multivariate Cox regression analyses, and the nomogram model was established to determine the potential prognostic value of CIC in LGG. Our study provides comprehensive analysis of mutated genes in LGG, supporting modulation of mutated genes in the management of LGG.     

Overexpression of oncostatin M receptor regulates local immune response in glioblastoma

Glioblastoma (GBM) is the most lethal cancer in central nervous system. It is urgently needed to look for novel therapeutics for GBM. Oncostatin M receptor (OSMR) is a cytokine receptor gene of IL-6 family and has been reported to be involved in regulating GBM tumorigenesis. However, the role of OSMR regulating the disrupted immune response in GBM need to be further investigated. Three gene expression profiles, Chinese Glioma Genome Atlas (CGGA), The Cancer Genome Atlas (TCGA), and Gene Expression Omnibus (GEO) data set (GSE16011), were enrolled in our study and used for OSMR expression and survival analysis. The expression of OSMR was further verified with immunohistochemistry and western blot analysis in glioma tissues. Microenvironment cell populations-counter (MCP-counter) was applied for analyzing the relationship between OSMR expression and nontumor cells. The functions of OSMR in GBM was investigated by Gene Ontology, Gene set enrichment analysis (GSEA), gene set variation analysis and so on. The analysis of cytokine receptor activity-related genes in glioma identifies OSMR as a gene with an independent predictive factor for progressive malignancy in GBM. Furthermore, OSMR expression is a prognostic marker in the response prediction to radiotherapy and chemotherapy. OSMR contributes to the regulation of local immune response and extracellular matrix process in GBM. Our findings define an important role of OSMR in the regulation of local immune response in GBM, which may suggest OSMR as a possible biomarker in developing new therapeutic immune strategies in GBM.     

Comprehensive analysis reveals a prognostic and therapeutic biomarker CD3D in the breast carcinoma microenvironment

Breast carcinoma (BRCA) is the most common carcinoma among women worldwide. Despite the great progress achieved in early detection and treatment, morbidity and mortality rates remain high. In the present study, we make a systematic analysis of BRCA using TCGA database by applying CIBERSORT and ESTIMATE computational methods, uncovered CD3D as a prognostic biomarker by intersection analysis of univariate COX and protein–protein interaction (PPI). It revealed that high CD3D expression was strongly associated with poor survival of BRCA, based on The Cancer Genome Atlas (TCGA) database and online websites. Gene Set Enrichment Analysis (GSEA) revealed that the high CD3D expression group was mainly enriched for the immune-related pathways and the low CD3D expression group was mainly enriched for metabolic-related activities. Based on CIBERSORT analysis, the difference test and correlation test suggested that CD3D had a strong correlation with T cells, particularly CD8 + T cells, which indicated that CD3D up-regulation may increase T cell immune infiltration in the TME and induce antitumor immunity by activating T lymphocytes. Furthermore, the correlation analysis showed that CD3D expression had a strongly positive correlation with immune checkpoints, which indicating that the underlying mechanism involves CD3D mediated regulation of T cell functions in BRCA, and single cell RNA-seq analysis revealed that CD3D correlate with CD8 + T cells and it is itself highly expressed in CD8 + T cells. In summary, we identified a prognostic biomarker CD3D in BRCA, which was associated with lymphocyte infiltration, immune checkpoints and could be developed for innovative therapeutics of BRCA.     

Characterization of an endoplasmic reticulum stress-related signature to evaluate immune features and predict prognosis in glioma

Endoplasmic reticulum (ER) stress has considerable impact on cell growth, proliferation, metastasis, invasion, angiogenesis and chemoradiotherapy resistance in various cancers. However, the effect of ER stress on the outcomes of glioma patients remains unclear. In this study, we established an ER stress risk model based on The Cancer Genome Atlas (TCGA) glioma data set to reflect immune characteristics and predict the prognosis of glioma patients. Survival analysis indicated that there were significant differences in the overall survival (OS) of glioma patients with different ER stress-related risk scores. Moreover, the ER stress-related risk signature, which was markedly associated with the clinicopathological properties of glioma patients, could serve as an independent prognostic indicator. Functional enrichment analysis revealed that the risk model correlated with immune and inflammation responses, as well as biosynthesis and degradation. In addition, the ER stress-related risk model indicated an immunosuppressive microenvironment. In conclusion, we present an ER stress risk model that is an independent prognostic factor and indicates the general immune characteristics in the glioma microenvironment.     

Prognostic power of a lipid metabolism gene panel for diffuse gliomas

Lipid metabolism reprogramming plays important role in cell growth, proliferation, angiogenesis and invasion in cancers. However, the diverse lipid metabolism programmes and prognostic value during glioma progression remain unclear. Here, the lipid metabolism‐related genes were profiled using RNA sequencing data from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) database. Gene ontology (GO) and gene set enrichment analysis (GSEA) found that glioblastoma (GBM) mainly exhibited enrichment of glycosphingolipid metabolic progress, whereas lower grade gliomas (LGGs) showed enrichment of phosphatidylinositol metabolic progress. According to the differential genes of lipid metabolism between LGG and GBM, we developed a nine‐gene set using Cox proportional hazards model with elastic net penalty, and the CGGA cohort was used for validation data set. Survival analysis revealed that the obtained gene set could differentiate the outcome of low‐ and high‐risk patients in both cohorts. Meanwhile, multivariate Cox regression analysis indicated that this signature was a significantly independent prognostic factor in diffuse gliomas. Gene ontology and GSEA showed that high‐risk cases were associated with phenotypes of cell division and immune response. Collectively, our findings provided a new sight on lipid metabolism in diffuse gliomas.     

A five‐microRNA signature for individualized prognosis evaluation and radiotherapy guidance in patients with diffuse lower‐grade glioma

Despite the prognostic value of IDH and other gene mutations found in diffuse glioma, markers that judge individual prognosis of patients with diffuse lower‐grade glioma (LGG) are still lacking. This study aims to develop an expression‐based microRNA signature to provide survival and radiotherapeutic response prediction for LGG patients. MicroRNA expression profiles and relevant clinical information of LGG patients were downloaded from The Cancer Genome Atlas (TCGA; the training group) and the Chinese Glioma Genome Atlas (CGGA; the test group). Cox regression analysis, random survival forests‐variable hunting (RSFVH) screening and receiver operating characteristic (ROC) were used to identify the prognostic microRNA signature. ROC and TimeROC curves were plotted to compare the predictive ability of IDH mutation and the signature. Stratification analysis was conducted in patients with radiotherapy information. Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to explore the biological function of the signature. We identified a five‐microRNA signature that can classify patients into low‐risk or high‐risk group with significantly different survival in the training and test datasets (P < 0.001). The five‐microRNA signature was proved to be superior to IDH mutation in survival prediction (AUCtraining = 0.688 vs 0.607). Stratification analysis found the signature could further divide patients after radiotherapy into two risk groups. GO and KEGG analyses revealed that microRNAs from the prognostic signature were mainly enriched in cancer‐associated pathways. The newly discovered five‐microRNA signature could predict survival and radiotherapeutic response of LGG patients based on individual microRNA expression.     

Identification of immune-enhanced molecular subtype associated with BRCA1 mutations,immune checkpoints and clinical outcome in ovarian carcinoma

Ovarian carcinoma has the highest mortality among the malignant tumours in gynaecology, and new treatment strategies are urgently needed to improve the clinical status of ovarian carcinoma patients. The Cancer Genome Atlas (TCGA) cohort were performed to explore the immune function of the internal environment of tumours and its clinical correlation with ovarian carcinoma. Finally, four molecular subtypes were obtained based on the global immune-related genes. The correlation analysis and clinical characteristics showed that four subtypes were all significantly related to clinical stage; the immune scoring results indicated that most immune signatures were upregulated in C3 subtype, and the majority of tumour-infiltrating immune cells were upregulated in both C3 and C4 subtypes. Compared with other subtypes, C3 subtype had a higher BRCA1 mutation, higher expression of immune checkpoints, and optimal survival prognosis. These findings of the immunological microenvironment in tumours may provide new ideas for developing immunotherapeutic strategies for ovarian carcinoma.     

Identification of UGP2 as a progression marker that promotes cell growth and motility in human glioma

Glioblastoma (GBM) is one of the most common highly malignant primary brain tumor with poor prognosis. This study aimed to explore the possible mechanism by bioinformatics method and detect potential function of UGP2 of GBM. Gene expression microarray data of GSE4412 and messenger RNA-sequencing data of GBM with samples clinical information were downloaded from the Gene Expression Omnibus database and The Cancer Genome Atlas database, respectively. Differentially expressed genes (DEGs) analysis using the Kyoto Encyclopedia of Genes and Genomes and Gene Ontology based on R language. A total of 1000 common DEGs were identified in GBM samples, including 353 upregulated and 647 downregulated genes. Based on the random survival forest model, we identified UDP-glucose pyrophosphorylase 2 (UGP2) (upregulated gene) had a significant effect on GBM prognosis. Functional enrichment showed that UGP2 was enriched in the biological progresses of cell proliferation, migration, and invasion. Furthermore, UGP2 expression is aberrantly overexpressed in human glioma and positively correlated with pathologic grade. A loss-of-function study showed that knockdown of UGP2 decreases U251 cell growth, migration, and invasion in vivo and vitro. We proposed the development and progression of human glioma were associated with survival based on bioinformatics analysis. We also found that UGP2 might function as prognostic markers in the pathogenesis of GBM.     

Metabolic reprogramming-associated genes predict overall survival for rectal cancer

Metabolic reprogramming has become a hot topic recently in the regulation of tumour biology. Although hundreds of altered metabolic genes have been reported to be associated with tumour development and progression, the important prognostic role of these metabolic genes remains unknown. We downloaded messenger RNA expression profiles and clinicopathological data from The Cancer Genome Atlas and the Gene Expression Omnibus database to uncover the prognostic role of these metabolic genes. Univariate Cox regression analysis and lasso Cox regression model were utilized in this study to screen prognostic associated metabolic genes. Patients with high-risk demonstrated significantly poorer survival outcomes than patients with low-risk in the TCGA database. Also, patients with high-risk still showed significantly poorer survival outcomes than patients with low-risk in the GEO database. What is more, gene set enrichment analyses were performed in this study to uncover significantly enriched GO terms and pathways in order to help identify potential underlying mechanisms. Our study identified some survival-related metabolic genes for rectal cancer prognosis prediction. These genes might play essential roles in the regulation of metabolic microenvironment and in providing significant potential biomarkers in metabolic treatment.     

基于单细胞转录组的多级别胶质瘤异质性及免疫微环境分析揭示了潜在的预后生物标志物

脑胶质瘤（glioma）是中枢神经系统最常见的内在肿瘤,具有发病率高、预后较差等特点。本研究旨在鉴定多形性胶质母细胞瘤（glioblastoma multiforme,GBM）和低级别胶质瘤（lower-grade gliomas,LGG）之间的差异表达基因（differentially expressed genes,DEGs）,以探讨不同级别胶质瘤的预后影响因素。从NCBI基因表达综合数据库中收集了胶质瘤的单细胞转录组测序数据,其中包括来自3个数据集的共29 097个细胞样本。对于不同分级的人脑胶质瘤进行分析,经过滤得到21 071个细胞,通过基因本体分析、京都基因与基因组百科全书途径分析,从差异表达基因中筛选出70个基因,我们通过查阅文献,聚焦到delta样典型Notch配体3（delta like canonical Notch ligand 3,DLL3）这个基因。基于TCGA的基因表达谱交互分析（gene expression profiling interactive analysis,GEPIA）数据库用于探索LGG和GBM中DLL3基因的表达差异,采用基因表达谱交互式分析和肿瘤免疫学估计资源（tumor immune estimation resource,TIMER）数据库,研究关键基因在不同分级的脑胶质瘤中的表达,预测了与免疫治疗密切相关的生物标志物。cBioPortal数据库用于探索DLL3表达与25个免疫检查点之间的关系。基因集富集分析（gene set enrichment analysis,GSEA）进一步确定了与中心基因相关的途径。最后,在中国胶质瘤基因组图谱（Chinese glioma genome atlas,CGGA）中验证了生物标志物在预后和预测中的疗效。这些结果发现,预后基因与肿瘤增殖和进展有关,通过生物学信息和生存分析,表明这些基因可能作为一种有前途的预后生物标志物,并作为选择治疗策略的新靶点。     

TGFβ2 is a prognostic‐related biomarker and correlated with immune infiltrates in gastric cancer

TGFβ2 is an essential regulator of immune cell functionality, but the mechanisms whereby it drives immune infiltration in gastric cancer remain uncertain. The Oncomine and Tumor Immunoassay Resource (TIMER) databases were used for assessing the expression of TGFβ2, after which TIMER was used to explore the relationship between TGFβ2 and tumour immune infiltration. Finally, we assessed how TGFβ2 expression correlated with the expression of a set of marker genes associated with immune infiltration using TIMER and GEPIA. We determined TGFβ2 expression to be significantly correlated with outcome in multiple types of cancer in the Cancer Genome Atlas (TCGA), with the effect being particularly pronounced in gastric cancer. Furthermore, elevated TGFβ2 expression was found to be significantly correlated with gastric cancer N staging, and with the expression of a variety of immune markers associated with particular immune cell subsets. These results indicate that TGFΒ2 is associated with patient outcome and tumour immune cell infiltration in multiple cancer types. This suggests that TGFβ2 is a key factor which governs immune cell recruitment to gastric cancer tumours, potentially playing a vital role in governing immune cell infiltration and thus representing a valuable prognostic biomarker in gastric cancer patients.     

Molecular prognostic of nine parthanatos death-related genes in glioma,particularly in COL8A1 identification

Post-operative progression and chemotherapy resistance are the main causes of treatment failure in glioma patients. There is a lack of ideal prediction models for post-operative glioma patient progression and drug sensitivity. We aimed to develop a prognostic model of parthanatos mRNA biomarkers for glioma outcomes. A total of 11 parthanatos genes were obtained from ParthanatosCluster database. ConsensusClusterPlus and R “Limma” package were used to cluster The Cancer Genome Atlas (TCGA)-glioma cohort and analyze the differential mRNAs. Univariate Cox regression analysis, random survival forest model, and least absolute shrinkage and selection operator (LASSO) regression analysis were used to determine the nine ParthanatosScore prognostic genes combination. ParthanatosScore was verified by 656 patients and 979 patients in TCGA and CGCA-LGG/GBM datasets. Differences in genomic mutations, tumor microenvironments, and functional pathways were assessed. Drug response prediction was performed using pRRophetic. Kaplan–Meier survival analysis was analyzed. Finally, COL8A1 was selected to evaluate its potential biological function and drug sensitivity of temozolomide and AZD3759 in glioma cells. ParthanatosScore obtained a combination of nine glioma prognostic genes, including CD58, H19, TNFAIP6, FTLP3, TNFRSF11B, SFRP2, LOXL1, COL8A1, and FABP5P7. In the TCGA-LGG/GBM dataset, glioma prognosis was poor in high ParthanatosScore. Low-score glioma patients were sensitive to AZD3759_1915, AZD5582_1617, AZD8186_1918, Dasatinib_1079, and Temozolomide_1375, while high-score patients were less sensitive to these drugs. Compared with HA cells, COL8A1 was significantly over-expressed in LN229 and U251 cells. Silencing COL8A1 inhibited the malignant characterization of LN229 and U251 cells. Temozolomide and AZD3759 also promoted parthanatos gene expression in glioma cells. Temozolomide and AZD3759 inhibited COL8A1 expression and cell viability and promoted apoptosis in glioma cells and PGM cells. ParthanatosScore can accurately predict clinical prognosis and drug sensitivity after glioma surgery. Silencing COL8A1 inhibited the malignant characterization. Temozolomide and AZD3759 inhibited COL8A1 expression and cell viability and promoted apoptosis and parthanatos gene expression, which is a target to improve glioma.     

Protein disulphide isomerase can predict the clinical prognostic value and contribute to malignant progression in gliomas

Increasing evidence from structural and functional studies has indicated that protein disulphide isomerase (PDI) has a critical role in the proliferation, survival and metastasis of several types of cancer. However, the molecular mechanisms through which PDI contributes to glioma remain unclear. Here, we aimed to investigate whether the differential expression of 17 PDI family members was closely related to the different clinicopathological features in gliomas from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas data sets. Additionally, four subgroups of gliomas (cluster 1/2/3/4) were identified based on consensus clustering of the PDI gene family. These findings not only demonstrated that a poorer prognosis, higher WHO grade, lower frequency of isocitrate dehydrogenase mutation and higher 1p/19q non-codeletion status were significantly correlated with cluster 4 compared with the other clusters, but also indicated that the malignant progression of glioma was closely correlated with the expression of PDI family members. Moreover, we also constructed an independent prognostic marker that can predict the clinicopathological features of gliomas. Overall, the results indicated that PDI family members may serve as possible diagnostic markers in gliomas.     

Comprehensive analysis of prognostic immune-related genes in the tumor microenvironment of cutaneous melanoma

Cutaneous malignant melanoma (hereafter called melanoma) is one of the most aggressive cancers with increasing incidence and mortality rates worldwide. In this study, we performed a systematic investigation of the tumor microenvironmental and genetic factors associated with melanoma to identify prognostic biomarkers for melanoma. We calculated the immune and stromal scores of melanoma patients from the Cancer Genome Atlas (TCGA) using the ESTIMATE algorithm and found that they were closely associated with patients’ prognosis. Then the differentially expressed genes were obtained based on the immune and stromal scores, and prognostic immune-related genes further identified. Functional analysis and the protein–protein interaction network further revealed that these genes enriched in many immune-related biological processes. In addition, the abundance of six infiltrating immune cells was analyzed using prognostic immune-related genes by TIMER algorithm. The unsupervised clustering analysis using immune-cell proportions revealed eight clusters with distinct survival patterns, suggesting that dendritic cells were most abundant in the microenvironment and CD8⁺ T cells and neutrophils were significantly related to patients’ prognosis. Finally, we validated these genes in three independent cohorts from the Gene Expression Omnibus database. In conclusion, this study comprehensively analyzed the tumor microenvironment and identified prognostic immune-related biomarkers for melanoma.     

Exploration of a novel prognostic risk signatures and immune checkpoint molecules in endometrial carcinoma microenvironment

In this study, we devoted to investigate immune-related genes and tumor microenvironment (TME) in EC based on The Cancer Genome Atlas (TCGA) database. In total 799 up-regulated and 139 down-regulated immune-related and differentially expressed genes in EC were investigated for functional annotations and prognosis. By a conjoint Cox regression analysis, we built two risk models for OS and DFS, as well as the consistent nomograms. Immune-related pathways were revealed mostly enriched in the low-risk group. By further analyzing TME based on the risk signatures, the higher immune cell infiltration and activation, lower tumor purity and higher tumor mutational burden were found in low-risk group, which presented a better prognosis. Both the expression and immunophenoscore of immune checkpoints PD-1, CTLA4, PD-L1 and PD-L2 increased significantly in low-risk group. These findings may provide new ideas for novel biomarkers and immunotherapy targets in EC.     

SAMSN1 Is Highly Expressed and Associated with a Poor Survival in Glioblastoma Multiforme

Objectives

To study the expression pattern and prognostic significance of SAMSN1 in glioma.

Methods

Affymetrix and Arrystar gene microarray data in the setting of glioma was analyzed to preliminarily study the expression pattern of SAMSN1 in glioma tissues, and Hieratical clustering of gene microarray data was performed to filter out genes that have prognostic value in malignant glioma. Survival analysis by Kaplan-Meier estimates stratified by SAMSN1 expression was then made based on the data of more than 500 GBM cases provided by The Cancer Genome Atlas (TCGA) project. At last, we detected the expression of SAMSN1 in large numbers of glioma and normal brain tissue samples using Tissue Microarray (TMA). Survival analysis by Kaplan-Meier estimates in each grade of glioma was stratified by SAMSN1 expression. Multivariate survival analysis was made by Cox proportional hazards regression models in corresponding groups of glioma.

Results

With the expression data of SAMSN1 and 68 other genes, high-grade glioma could be classified into two groups with clearly different prognoses. Gene and large sample tissue microarrays showed high expression of SAMSN1 in glioma particularly in GBM. Survival analysis based on the TCGA GBM data matrix and TMA multi-grade glioma dataset found that SAMSN1 expression was closely related to the prognosis of GBM, either PFS or OS (P<0.05). Multivariate survival analysis with Cox proportional hazards regression models confirmed that high expression of SAMSN1 was a strong risk factor for PFS and OS of GBM patients.

Conclusion

SAMSN1 is over-expressed in glioma as compared with that found in normal brains, especially in GBM. High expression of SAMSN1 is a significant risk factor for the progression free and overall survival of GBM.