Dataset Information

Identification of Prognostic Genes in the Tumor Microenvironment of Hepatocellular Carcinoma.

ABSTRACT: Background: Hepatocellular carcinoma (HCC) is one of the most common malignant tumors in the world. The efficacy of immunotherapy usually depends on the interaction of immunomodulation in the tumor microenvironment (TME). This study aimed to explore the potential stromal-immune score-based prognostic genes related to immunotherapy in HCC through bioinformatics analysis. Methods: ESTIMATE algorithm was applied to calculate the immune/stromal/Estimate scores and tumor purity of HCC using the Cancer Genome Atlas (TCGA) transcriptome data. Functional enrichment analysis of differentially expressed genes (DEGs) was analyzed by the Database for Annotation, Visualization, and Integrated Discovery database (DAVID). Univariate and multivariate Cox regression analysis and least absolute shrinkage and selection operator (LASSO) regression analysis were performed for prognostic gene screening. The expression and prognostic value of these genes were further verified by KM-plotter database and the Human Protein Atlas (HPA) database. The correlation of the selected genes and the immune cell infiltration were analyzed by single sample gene set enrichment analysis (ssGSEA) algorithm and Tumor Immune Estimation Resource (TIMER). Results: Data analysis revealed that higher immune/stromal/Estimate scores were significantly associated with better survival benefits in HCC within 7 years, while the tumor purity showed a reverse trend. DEGs based on both immune and stromal scores primarily affected the cytokine-cytokine receptor interaction signaling pathway. Among the DEGs, three genes (CASKIN1, EMR3, and GBP5) were found most significantly associated with survival. Moreover, the expression levels of CASKIN1, EMR3, and GBP5 genes were significantly correlated with immune/stromal/Estimate scores or tumor purity and multiple immune cell infiltration. Among them, GBP5 genes were highly related to immune infiltration. Conclusion: This study identified three key genes which were related to the TME and had prognostic significance in HCC, which may be promising markers for predicting immunotherapy outcomes.

SUBMITTER: Xiang S

PROVIDER: S-EPMC8059369 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Purpose: The purpose of this study was to construct a novel risk scoring model with prognostic value that could elucidate tumor immune microenvironment of hepatocellular carcinoma (HCC). Samples and methods: Data were obtained through The Cancer Genome Atlas (TCGA) database. Univariate Cox analysis, least absolute shrinkage and selection operator (LASSO) analysis, and multivariate Cox analysis were carried out to screen for glycolysis-related long noncoding RNAs (lncRNAs) that could provide prognostic value. Finally, we established a risk score model to describe the characteristics of the model and verify its prediction accuracy. The receiver operating characteristic (ROC) curves of 1, 3, and 5 years of overall survival (OS) were depicted with risk score and some clinical features. ESTIMATE algorithm, single-sample gene set enrichment analysis (ssGSEA), and CIBERSORT analysis were employed to reveal the characteristics of tumor immune microenvironment in HCC. The nomogram was drawn by screening indicators with high prognostic accuracy. The correlation of risk signature with immune infiltration and immune checkpoint blockade (ICB) therapy was analyzed. After enrichment of related genes, active behaviors and pathways in high-risk groups were identified and lncRNAs related to poor prognosis were validated in vitro. Finally, the impact of MIR4435-2HG upon ICB treatment was uncovered. Results: After screening through multiple steps, four glycolysis-related lncRNAs were obtained. The risk score constructed with the four lncRNAs was found to significantly correlate with prognosis of samples. From the ROC curve of samples with 1, 3, and 5 years of OS, two indicators were identified with high prognostic accuracy and were used to draw a nomogram. Besides, the risk score significantly correlated with immune score, immune-related signature, infiltrating immune cells (i.e. B cells, etc.), and ICB key molecules (i.e. CTLA4,etc.). Gene enrichment analysis indicated that multiple biological behaviors and pathways were active in the high-risk group. In vitro validation results showed that MIR4435-2HG was highly expressed in the two cell lines, which had a significant impact on the OS of samples. Finally, we corroborated that MIR4435-2HG had intimate relationship with ICB therapy in hepatocellular carcinoma. Conclusion: We elucidated the crucial role of risk signature in immune cell infiltration and immunotherapy, which might contribute to clinical strategies and clinical outcome prediction of HCC.

Project description:AbstractGrowing evidence supports that the tumor microenvironment plays a key role in the development and progression of tumors. But immune microenvironment of hepatocellular carcinoma (HCC) has not yet been fully explored. In the present investigation, the clinical value and prognostic significance of immune-related genes in HCC were investigated.The immune and stromal scores of HCC were calculated through the application of Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression data Algorithm based on the Cancer Genome Atlas database. Differentially expressed genes were identified using the "edgeR" package of the R software. Functional annotation and pathway enrichment were performed using "ggplots2" and "clusterProfiler" packages in R software. Protein-protein interaction network was constructed using STRING, and the hub genes were identified through the Cytoscape. Survival analysis was performed using Kaplan-Meier methods. Tumor Immune Estimation Resource algorithm was used to view the immune landscape of the microenvironment in HCC.Firstly, the immune and stromal scores of HCC were calculated and we found that the immune and stromal scores of HCC were closely related to the patients' prognosis. Then the differentially expressed genes were identified respectively stratified by the median value of the immune and stromal scores, and the immune-related genes that related to the prognosis in HCC patients were further identified. Functional enrichment analysis and protein-protein interaction networks further showed that these genes mainly participated in immune-related biological process. In addition, dendritic cells were found to be the most abundant in the microenvironment of HCC through Tumor Immune Estimation Resource algorithm and were significantly associated with the patients' prognosis. To robust the results, the immune-related genes were validated in an independent dataset from the Gene Expression Omnibus database.We arrived at a more comprehensive understanding of the microenvironment of HCC and extracted 7 immune-related genes that were significantly associated with the recurrence survival of HCC.

Project description:BackgroundHepatocellular carcinoma (HCC) is a leading cause of cancer-related mortality worldwide, but effective early detection and prognostication methods are lacking.MethodsThe Cox regression model was built to stratify the HCC patients. The single-cell RNA sequencing data analysis and gene set enrichment analysis were employed to investigate the biological function of identified markers. PLCB1 gain- or loss-of-function experiments were performed, and obtained HCC samples were analyzed using quantitative real-time PCR and immunohistochemistry assay to validate the biological function of identified markers.ResultsIn this study, we developed a model using optimized markers for HCC recurrence prediction. Specifically, we screened out 8 genes through a series of data analyses, and built a multivariable Cox model based on their expression. The risk stratifications using the Eight-Gene Cox (EGC) model were closely associated with the recurrence-free survivals (RFS) in both training and three validation cohorts. We further demonstrated that this risk stratification could serve as an independent predictor in predicting HCC recurrence, and that the EGC model could outperform other models. Moreover, we also investigated the cell-type-specific expression patterns of the eight recurrence-related genes in tumor microenvironment using single-cell RNA sequencing data, and interpreted their functional roles from correlation and gene set enrichment analyses, in vitro and in vivo experiments. Particularly, PLCB1 and SLC22A7 were predominantly expressed in malignant cells, and they were predicted to promote angiogenesis and to help maintain normal metabolism in liver, respectively. In contrast, both FASLG and IL2RB were specifically expressed in T cells, and were highly correlated with T cell marker genes, suggesting that these two genes might assist in maintaining normal function of T cell-mediated immune response in tumor tissues.ConclusionIn conclusion, the EGC model and eight identified marker genes could not only facilitate the accurate prediction of HCC recurrence, but also improve our understanding of the mechanisms behind HCC recurrence.

Project description:Background: In recent years, immunotherapy has changed the therapeutic landscape of hepatocellular carcinoma (HCC). Since the efficacy of immunotherapy is closely related to the tumor microenvironment (TME), in this study, we constructed a prognostic model based on TME to predict the prognosis and immunotherapy effect of HCC patients. Methods: Transcriptome and follow-up data of 374 HCC patients were acquired from the TCGA Cancer Genome Atlas (TCGA) database. The immune/stromal/estimate scores (TME scores) and tumor purity were calculated using the ESTIMATE algorithm and the module most associated with TME scores were screened by the weighted gene co-expression network analysis (WGCNA). A TME score-related prognostic model was constructed and patients were divided into a high-risk group and a low-risk group. Kaplan-Meier survival curves and receiver operator characteristic curve (ROC) were used to evaluate the performance of the TME risk prognostic model and validated with the external database International Cancer Genome Consortium (ICGC) cohort. Combined with clinicopathologic factors, a prognostic nomogram was established. The nomogram's ability to predict prognosis was assessed by ROC, calibration curve, and the decision curve analysis (DCA). Gene Set Enrichment Analyses (GSEA) were conducted to explore the underlying biological functions and pathways of this risk signature. Moreover, the possible correlation of risk signature with TME immune cell infiltration, immune checkpoint inhibitor (ICI) treatment response, single-nucleotide polymorphisms (SNPs), and drug sensitivity were assessed. Finally, real-time PCR was used to verify the gene expression levels in normal liver cells and cancer cells. Results: KM survival analysis results indicated that high immune/stromal/estimate score groups were closely associated with a better prognosis, while the tumor purity showed a reverse trend (p < 0.01). WGCNA demonstrated that the yellow module was significantly correlated with the TME score. The 5-genes TME risk signature was built to predict the prognosis of patients with HCC including DAB2, IL18RAP, RAMP3, FCER1G, and LHFPL2. Patients with a low-risk score have higher levels of tumor-infiltrating immune cells and higher expression of immune checkpoints, which may be more sensitive to immunotherapy. Conclusion: It provided a theoretical basis for predicting the prognosis and personalized treatment of patients with HCC.

Project description:AimImmune cells that infiltrate the tumor microenvironment (TME) are associated with cancer prognosis. The aim of the current study was to identify TME related gene signatures related to the prognosis of sarcoma (SARC) by using the data from The Cancer Genome Atlas (TCGA).MethodsImmune and stromal scores were calculated by estimation of stromal and immune cells in malignant tumor tissues using expression data algorithms. The least absolute shrinkage and selection operator (lasso) based cox model was then used to select hub survival genes. A risk score model and nomogram were used to predict the overall survival of patients with SARC.ResultsWe selected 255 patients with SARC for our analysis. The Kaplan-Meier method found that higher immune (p = 0.0018) or stromal scores (p = 0.0022) were associated with better prognosis of SARC. The estimated levels of CD4+ (p = 0.0012) and CD8+ T cells (p = 0.017) via the tumor immune estimation resource were higher in patients with SARC with better overall survival. We identified 393 upregulated genes and 108 downregulated genes (p < 0.05, fold change >4) intersecting between the immune and stromal scores based on differentially expressed gene (DEG) analysis. The univariate Cox analysis of each intersecting DEG and subsequent lasso-based Cox model identified 11 hub survival genes (MYOC, NNAT, MEDAG, TNFSF14, MYH11, NRXN1, P2RY13, CXCR3, IGLV3-25, IGHV1-46, and IGLV2-8). Then, a hub survival gene-based risk score gene signature was constructed; higher risk scores predicted worse SARC prognosis (p < 0.0001). A nomogram including the risk scores, immune/stromal scores and clinical factors showed a good prediction value for SARC overall survival (C-index = 0.716). Finally, connectivity mapping analysis identified that the histone deacetylase inhibitors trichostatin A and vorinostat might have the potential to reverse the harmful TME for patients with SARC.ConclusionThe current study provided new indications for the association between the TME and SARC. Lists of TME related survival genes and potential therapeutic drugs were identified for SARC.

Dataset Information

Identification of Prognostic Genes in the Tumor Microenvironment of Hepatocellular Carcinoma.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets