Dataset Information

Integrated Analysis Reveals Prognostic Value and Immune Correlates of CD86 Expression in Lower Grade Glioma.

ABSTRACT:

Background

CD86 has great potential to be a new target of immunotherapy by regulating cancer immune response. However, it remains unclear whether CD86 is a friend or foe in lower-grade glioma (LGG).

Methods

The prognostic value of CD86 expression in pan-cancer was analyzed using Cox regression and Kaplan-Meier analysis with data from the cancer genome atlas (TCGA). Cancer types where CD86 showed prognostic value in overall survival and disease-specific survival were identified for further analyses. The Chinese Glioma Genome Atlas (CGGA) dataset were utilized for external validation. Quantitative real-time PCR (qRT-PCR), Western blot (WB), and Immunohistochemistry (IHC) were conducted for further validation using surgical samples from Jiangsu Province hospital. The correlations between CD86 expression and tumor immunity were analyzed using the Estimation of Stromal and Immune cells in Malignant Tumours using Expression data (ESTIMATE) algorithm, Tumor IMmune Estimation Resource (TIMER) database, and expressions of immune checkpoint molecules. Gene Set Enrichment Analysis (GSEA) was performed using clusterprofiler r package to reveal potential pathways.

Results

Pan-cancer survival analysis established CD86 expression as an unfavorable prognostic factor in tumor progression and survival for LGG. CD86 expression between Grade-II and Grade-III LGG was validated using qRT-PCR and WB. Additionally, CD86 expression in LGG with unmethylated O(6)-methylguanine-DNA-methyltransferase (MGMT) promoter was significantly higher than those with methylated MGMT (P<0.05), while in LGG with codeletion of 1p/19q it was significantly downregulated as opposed to those with non-codeletion (P<2.2*10-16). IHC staining validated that CD86 expression was correlated with MGMT status and X1p/19q subtypes, which was independent of tumor grade. Multivariate regression validated that CD86 expression acts as an unfavorable prognostic factor independent of clinicopathological factors in overall survival of LGG patients. Analysis of tumor immunity and GSEA revealed pivotal role of CD86 in immune response for LGG.

Conclusions

Integrated analysis shows that CD86 is an unfavorable prognostic biomarker in LGG patients. Targeting CD86 may become a novel approach for immunotherapy of LGG.

SUBMITTER: Qiu H

PROVIDER: S-EPMC8089378 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:ObjectiveAs a prevalent and infiltrative cancer type of the central nervous system, the prognosis of lower-grade glioma (LGG) in adults is highly heterogeneous. Recent evidence has demonstrated the prognostic value of the mRNA expression-based stemness index (mRNAsi) in LGG. Our aim was to develop a stemness index-based signature (SI-signature) for risk stratification and survival prediction.MethodsDifferentially expressed genes (DEGs) between LGG in the Cancer Genome Atlas (TCGA) and normal brain tissue samples from the Genotype-Tissue Expression (GTEx) project were screened out, and the weighted gene correlation network analysis (WGCNA) was employed to identify the mRNAsi-related gene sets. Meanwhile, the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed for the functional annotation of the key genes. ESTIMATE was used to calculate tumor purity for acquiring the correct mRNAsi. Differences in overall survival (OS) between the high and low mRNAsi (corrected mRNAsi) groups were compared using the Kaplan Meier analysis. By combining the Lasso regression with univariate and multivariate Cox regression, the SI-signature was constructed and validated using the Chinese Glioma Genome Atlas (CGGA).ResultsThere was a significant difference in OS between the high and low mRNAsi groups, which was also observed in the two corrected mRNAsi groups. Based on threshold limits, 86 DEGs were most significantly associated with mRNAsi via WGCNA. Seven genes (ADAP2, ALOX5AP, APOBEC3C, FCGRT, GNG5, LRRC25, and SP100) were selected to establish a risk signature for primary LGG. The ROC curves showed a fair performance in survival prediction in both the TCGA and the CGGA validation cohorts. Univariate and multivariate Cox regression revealed that the risk group was an independent prognostic factor in primary LGG. The nomogram was developed based on clinical parameters integrated with the risk signature, and its accuracy for predicting 3- and 5-years survival was assessed by the concordance index, the area under the curve of the time-dependent receiver operating characteristics curve, and calibration curves.ConclusionThe SI-signature with seven genes could serve as an independent predictor, and suggests the importance of stemness features in risk stratification and survival prediction in primary LGG.

Project description:BackgroundAbnormal expression of the eukaryotic initiation factor 3 (eIF3) subunits plays critical roles in tumorigenesis and progression, and also has potential prognostic value in cancers. However, the expression and clinical implications of eIF3 subunits in glioma remain unknown.MethodsExpression data of eIF3 for patients with gliomas were obtained from the Chinese Glioma Genome Atlas (CGGA) (n = 272) and The Cancer Genome Atlas (TCGA) (n = 595). Cox regression, the receiver operating characteristic (ROC) curves and Kaplan-Meier analysis were used to study the prognostic value. Gene oncology (GO) and gene set enrichment analysis (GSEA) were utilized for functional prediction.ResultsIn both the CGGA and TCGA datasets, the expression levels of eIF3d, eIF3e, eIF3f, eIF3h and eIF3l highly were associated with the IDH mutant status of gliomas. The expression of eIF3b, eIF3i, eIF3k and eIF3m was increased with the tumor grade, and was associated with poorer overall survival [All Hazard ratio (HR) > 1 and P < 0.05]. By contrast, the expression of eIF3a and eIF3l was decreased in higher grade gliomas and was associated with better overall survival (Both HR < 1 and P < 0.05). Importantly, the expression of eIF3i (located on chromosome 1p) and eIF3k (Located on chromosome 19q) were the two highest risk factors in both the CGGA [eIF3i HR = 2.068 (1.425-3.000); eIF3k HR = 1.737 (1.166-2.588)] and TCGA [eIF3i HR = 1.841 (1.642-2.064); eIF3k HR = 1.521 (1.340-1.726)] databases. Among eIF3i, eIF3k alone or in combination, the expression of eIF3i was the more robust in stratifying the survival of glioma in various pathological subgroups. The expression of eIF3i was an independent prognostic factor in IDH-mutant lower grade glioma (LGG) and could also predict the 1p/19q codeletion status of IDH-mutant LGG. Finally, GO and GSEA analysis showed that the elevated expression of eIF3i was significantly correlated with the biological processes of cell proliferation, mRNA processing, translation, T cell receptor signaling, NF-κB signaling and others.ConclusionsOur study reveals the expression alterations during glioma progression, and highlights the prognostic value of eIF3i in IDH-mutant LGG.

Project description:Previous studies have shown that the prognosis of patients with lower-grade glioma (LGG) is closely related to the infiltration of immune cells and the expression of long non-coding RNAs (lncRNAs). In this paper, we applied single-sample gene set enrichment analysis (ssGSEA) algorithm to evaluate the expression level of immune genes from tumor tissues in The Cancer Genome Atlas (TCGA) database, and divided patients into the high immune group and the low immune group, which were separately analyzed for differential expression. Venn analysis was taken to select 36 immune-related lncRNAs. To construct a prognostic model of LGG based on immune-related lncRNAs, we divided patients into a training set and a verification set at a ratio of 2:1. Univariate Cox regression and the Least Absolute Shrinkage and Selection Operator (LASSO) regression were performed to select 11 immune-related lncRNAs associated with the prognosis of LGG, and based on these selected lncRNAs, the risk scoring model was constructed. Through Kaplan-Meier analysis, the overall survival (OS) of patients in the high-risk group was significantly lower than that of the low-risk group. Then, established a nomogram including age, gender, neoplasm histologic grade, and risk score. Meanwhile, the predictive performance of the model was evaluated by calculating the C-index, drawing the calibration chart, the clinical decision curve as well as the Receiver Operating Characteristic (ROC) curve. Similar results were obtained by utilizing the validation data to verify the above consequences. Based on the TIMER database, the correlation analysis showed that the 11 immune-related lncRNAs risk score of LGG were in connection with the infiltration of the subtypes of immune cells. Subsequently, we performed enrichment analysis, whose results showed that these immune-related lncRNAs played important roles in the progress of LGG. In conclusion, these 11 immune-related lncRNAs have the potential to predict the prognosis of patients with LGG, which may play a key role in the development of LGG.