Dataset Information

Identification of Prognostic Glycolysis-Related lncRNA Signature in Tumor Immune Microenvironment of Hepatocellular Carcinoma.

ABSTRACT: 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.

SUBMITTER: Bai Y

PROVIDER: S-EPMC8100457 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background: An accumulating body of evidence suggests that long non-coding RNAs (lncRNAs) can serve as potential cancer prognostic factors. However, the utility of lncRNA combinations in estimating overall survival (OS) for hepatocellular carcinoma (HCC) remains to be elucidated. This study aimed to construct a powerful lncRNA signature related to the OS for HCC to enhance prognostic accuracy. Methods: The expression patterns of lncRNAs and related clinical data of 371 HCC patients were obtained based on The Cancer Genome Atlas (TCGA). Differentially expressed lncRNAs (DElncRNAs) were acquired by comparing tumors with adjacent normal samples. lncRNAs displaying significant association with OS were screened through univariate Cox regression analysis and the least absolute shrinkage and selection operator (LASSO) algorithm. All cases were classified into the validation or training group at the ratio of 3:7 to validate the constructed lncRNA signature. Data from the Gene Expression Omnibus (GEO) were used for external validation. We conducted real-time polymerase chain reaction (PCR) and assays for Transwell invasion, migration, CCK-8, and colony formation to determine the biological roles of lncRNA. Gene set enrichment analysis (GSEA) of the lncRNA model risk score was also conducted. Results: We identified 1292 DElncRNAs, among which 172 were significant in univariate Cox regression analysis. In the training group (n = 263), LASSO regression analysis confirmed 11 DElncRNAs including AC010547.1, AC010280.2, AC015712.7, GACAT3 (gastric cancer associated transcript 3), AC079466.1, AC089983.1, AC051618.1, AL121721.1, LINC01747, LINC01517, and AC008750.3. The prognostic risk score was calculated, and the constructed risk model showed significant correlation with HCC OS (log-rank P-value of 8.489e-9, hazard ratio of 3.648, 95% confidence interval: 2.238-5.945). The area under the curve (AUC) for this lncRNA model was up to 0.846. This risk model was confirmed in the validation group (n = 108), the entire cohort, and the external GEO dataset (n = 203). GACAT3 was highly expressed in HCC tissues and cell lines. Based on online databases, GACAT3 expression independently affects both OS and disease-free survival in HCC patients. Silencing GACAT3 in vitro significantly suppressed HCC cell proliferation, invasion, and migration. Moreover, pathways related to the lncRNA model risk score were confirmed by GSEA. Conclusion: The lncRNA signature established in this study can be used to predict HCC prognosis, which could provide novel clinical evidence to guide targeted HCC treatment.

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:Objective: The present study aimed to establish a prognostic signature based on the autophagy-related long non-coding RNAs (lncRNAs) analysis in patients with hepatocellular carcinoma (HCC). Methods: Patients with HCC from The Cancer Genome Atlas (TCGA) were taken as the training cohort, and patients from the International Cancer Genome Consortium (ICGC) were treated as the validation cohort. Autophagy-related lncRNAs were obtained via a co-expression network analysis. According to univariate and multivariate analyses, a multigene prognostic signature was constructed in the training cohort. The predictive power of the signature was confirmed in both cohorts. The detailed functions were investigated using functional analysis. The single-sample gene set enrichment analysis (ssGSEA) score was used to evaluate the tumor microenvironment. The expression levels of immunotherapy and targeted therapy targets between the two risk groups were compared. Finally, a nomogram was constructed by integrating clinicopathological parameters with independently predictive value and the risk score. Results: Four autophagy-related lncRNAs were identified to establish a prognostic signature, which separated patients into high- and low-risk groups. Survival analysis showed that patients in the high-risk group had a shorter survival time in both cohorts. A time-independent receiver-operating characteristic (ROC) curve and principal component analysis (PCA) confirmed that the prognostic signature had a robust predictive power and reliability in both cohorts. Functional analysis indicated that the expressed genes in the high-risk group are mainly enriched in autophagy- and cancer-related pathways. ssGSEA revealed that the different risk groups were associated with the tumor microenvironment. Moreover, the different risk groups had positive correlations with the expressions of specific mutant genes. Multivariate analysis showed that the risk score also exhibited excellent predictive power irrespective of clinicopathological characteristics in both cohorts. A nomogram was established. The nomogram showed good discrimination, with Harrell's concordance index (C-index) of 0.739 and good calibration. Conclusion: The four autophagy-related lncRNAs could be used as biological biomarkers and therapeutic targets. The prognostic signature and nomogram might aid clinicians in individual treatment optimization and clinical decision-making for patients with HCC.

Project description: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.

Project description:Hepatocellular carcinoma (HCC) is one of the malignant tumors with high mortality and a worse prognosis globally. Necroptosis is a programmed death mediated by receptor-interacting Protein 1 (RIP1), receptor-interacting Protein 1 (RIP3), and Mixed Lineage Kinase Domain-Like (MLKL). Our study aimed to create a new Necroptosis-related lncRNAs (NRlncRNAs) risk model that can predict survival and tumor immunity in HCC patients. The RNA expression and clinical data originated from the TCGA database. Pearson correlation analysis was applied to identify the NRlncRNAs. The LASSO-Cox regression analysis was employed to build the risk model. Next, the ROC curve and the area under the Kaplan-Meier curve were utilized to evaluate the accuracy of the risk model. In addition, based on the two groups of risk model, we performed the following analysis: clinical correlation, differential expression, PCA, TMB, GSEA analysis, immune cells infiltration, and clinical drug prediction analysis. Plus, qRT-PCR was applied to test the expression of genes in the risk model. Finally, a prognosis model covering six necroptosis-related lncRNAs was constructed to predict the survival of HCC patients. The ROC curve results showed that the risk model possesses better accuracy. The 1, 3, and 5-years AUC values were 0.746, 0.712, and 0.670, respectively. Of course, we also observed that significant differences exist in the following analysis, such as functional signaling pathways, immunological state, mutation profiles, and medication sensitivity between high-risk and low-risk groups of HCC patients. The result of qRT-PCR confirmed that three NRlncRNAs were more highly expressed in HCC cell lines than in the normal cell line. In conclusion, based on the bioinformatics analysis, we constructed an NRlncRNAs associated risk model, which predicts the prognosis of HCC patients. Although our study has some limitations, it may greatly contribute to the treatment of HCC and medical progression.

Project description:BackgroundThe tumor microenvironment (TME) regulates the proliferation and metastasis of solid tumors and the effectiveness of immunotherapy against them. We investigated the prognostic role of TME-related genes based on transcriptomic data of bladder urothelial carcinoma (BLCA) and formulated a prediction model of TME-related signatures.MethodsMolecular subtypes were identified using the non-negative matrix factorization (NMF) algorithm based on TME-related genes from the TCGA database. TME-related genes with prognostic significance were screened with univariate Cox regression analysis and lasso regression. Nomogram was developed based on risk genes. Receiver operating characteristic (ROC) curve and decision curve analysis (DCA) were used for inner and outer validation of the model. Risk scores (RS) of patients were calculated and divided into high-risk group (HRG) and low-risk group (LRG) to compare the differences in clinical characteristics and PD-L1 treatment responsiveness between HRG and LRG.ResultsWe identified two molecular subtypes (C1 and C2) according to the NMF algorithm. There were significant differences in overall survival (OS) (p<0.05), progression-free survival (PFS) (p<0.05), and immune cell infiltration between the two subtypes. A total of eight TME-associated genes (CABP4, ZNF432, BLOC1S3, CXCL11, ANO9, OAS1, FBN2, CEMIP) with independent prognostic significance were screened to build prognostic risk models. Age (p<0.001), grade (p<0.001), and RS (p<0.001) were independent predictors of survival in BLCA patients. The developed RS nomogram was able to predict the prognosis of BLCA patients at 1, 3, and 5 years more potentially than the models of other investigators according to ROC and DCA. RS showed significantly higher values (p = 0.047) in patients with stable disease (SD)/progressive disease (PD) compared to patients with complete response (CR)/partial response (PR).ConclusionsWe successfully clustered and constructed predictive models for TME-associated genes and helped guide immunotherapy strategies.

Project description:BackgroundVarious studies highlighted that immune cell-mediated inflammatory processes play crucial roles in the progression and treatment of hepatocellular carcinoma (HCC). However, the immune microenvironment of HCC is still poorly characterized. Exploring the role of immune-related genes (IRGs) and describing the immune landscape in HCC would provide insights into tumor-immune co-evolution along HCC progression.MethodsWe integrated the datasets with complete prognostic information from the Cancer Genome Atlas (TCGA) database and GEO DataSets (GSE14520, GSE76427, and GSE54236) to construct a novel immune landscape based on the Cibersort algorithm and reveal the prognostic signature in HCC patients.ResultsTo describe the tumor microenvironment (TME) in HCC, immune infiltration patterns were defined using the CIBERSORT method, and a prognostic signature contains 5 types of immune cells, including 3 high-risk immune cells (T.cells. CD4. memory. resting, Macrophages.M0, Macrophages.M2) and 2 low-risk immune cells (Plasma. cells, T.cells.CD8), were finally constructed. A novel prognostic index, based on prognostic immune risk score (pIRG), was developed using the univariate Cox regression analyses and LASSO Cox regression algorithm. Furthermore, the ROC curve and KM curve showed that the TME signatures had a stable value in predicting the prognosis of HCC patients in the internal training cohort, internal validation, and external validation cohort. Differential genes analysis and qPCR experiment showed that the expression levels of AKR1B10, LAPTM4B, MMP9, and SPP1 were significantly increased in high-risk patients, while the expression of CD5L was lower. Further analysis found that AKR1B10 and MMP9 were associated with higher M0 macrophage infiltration, while CD5L was associated with higher plasma cell infiltration.ConclusionsTaken together, we performed a comprehensive evaluation of the immune landscape of HCC and constructed a novel and robust prognostic prediction model. AKR1B10, LAPTM4B, MMP9, SPP1, and CD5L were involved in important processes in the HCC tumor microenvironment and were expected to become HCC prediction markers and potential targets of treatment.

Dataset Information

Identification of Prognostic Glycolysis-Related lncRNA Signature in Tumor Immune Microenvironment of Hepatocellular Carcinoma.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets