Dataset Information

A six-miRNA signature as a novel biomarker for improving prediction of prognosis and patterns of immune infiltration in hepatocellular carcinoma.

ABSTRACT:

Background

Hepatocellular carcinoma (HCC) is one of the leading causes of tumor-related death. MicroRNAs (miRNAs) belong to a subfamily of functional non-coding RNAs (ncRNAs) and are essential regulators of tumorigenesis. They affect tumor-related therapeutic response, tumor metastasis, and clinical outcomes of several human malignant tumors. However, the prognostic value of miRNAs and their role in the tumor immune microenvironment (TIME) of HCC have not been clarified.

Materials and methods

Raw RNA-sequencing data (mRNA and miRNA) and clinicopathological characteristics of HCC samples were downloaded from the TCGA-GDC database. The Perl programming language, R software, Cytoscape software, and several online databases were used to clarify the clinical significance and biological functions of miRNAs and their target genes in HCC.

Results

A total of 424 mRNA-sequencing samples and 425 miRNA-sequencing samples were obtained from the TCGA database. There were 344 HCC cases with complete information in the TCGA dataset and they were randomly categorized into two subgroups. Six miRNAs were identified as independent prognostic biomarkers for HCC patients by univariate and multivariate Cox regression analysis. The constructed prognostic signature, which contains these six miRNAs, was significantly correlated with overall survival (OS). In addition, this prognostic signature is superior to single miRNA in predicting short-term prognosis of HCC patients. We also found that the prognostic signature was significantly associated with tumor-related immune cell infiltration, TIME, and immunotherapeutic response. Furthermore, a total of 4568 potential target genes of six miRNAs were identified. The miRNA-mRNA co-expression network, protein-protein interaction (PPI) network, and functional and pathway enrichment analysis demonstrated that these miRNA-related target genes have important biological effects during the initiation and progression of HCC.

Conclusions

This study demonstrates that the miRNA signature can accurately predict the prognosis of HCC patients and provide a basis for novel immunotherapy treatments.

SUBMITTER: Chen Z

PROVIDER: S-EPMC9274565 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundHepatocellular carcinoma (HCC) is a global health burden with poor prognosis. Anoikis, a novel programmed cell death, has a close interaction with metastasis and progression of cancer. In this study, we aimed to construct a novel bioinformatics model for evaluating the prognosis of HCC based on anoikis-related gene signatures as well as exploring the potential mechanisms.Materials and methodsWe downloaded the RNA expression profiles and clinical data of liver hepatocellular carcinoma from TCGA database, ICGC database and GEO database. DEG analysis was performed using TCGA and verified in the GEO database. The anoikis-related risk score was developed via univariate Cox regression, LASSO Cox regression and multivariate Cox regression, which was then used to categorize patients into high- and low-risk groups. Then GO and KEGG enrichment analyses were performed to investigate the function between the two groups. CIBERSORT was used for determining the fractions of 22 immune cell types, while the ssGSEA analyses was used to estimate the differential immune cell infiltrations and related pathways. The "pRRophetic" R package was applied to predict the sensitivity of administering chemotherapeutic and targeted drugs.ResultsA total of 49 anoikis-related DEGs in HCC were detected and 3 genes (EZH2, KIF18A and NQO1) were selected out to build a prognostic model. Furthermore, GO and KEGG functional enrichment analyses indicated that the difference in overall survival between risk groups was closely related to cell cycle pathway. Notably, further analyses found the frequency of tumor mutations, immune infiltration level and expression of immune checkpoints were significantly different between the two risk groups, and the results of the immunotherapy cohort showed that patients in the high-risk group have a better immune response. Additionally, the high-risk group was found to have higher sensitivity to 5-fluorouracil, doxorubicin and gemcitabine.ConclusionThe novel signature of 3 anoikis-related genes (EZH2, KIF18A and NQO1) can predict the prognosis of patients with HCC, and provide a revealing insight into personalized treatments in HCC.

Project description:Background:The prognosis of hepatocellular carcinoma (HCC) patients remains poor. Identifying prognostic markers to stratify HCC patients might help to improve their outcomes. Methods:Six gene expression profiles (GSE121248, GSE84402, GSE65372, GSE51401, GSE45267 and GSE14520) were obtained for differentially expressed genes (DEGs) analysis between HCC tissues and non-tumor tissues. To identify the prognostic genes and establish risk score model, univariable Cox regression survival analysis and Lasso-penalized Cox regression analysis were performed based on the integrated DEGs by robust rank aggregation method. Then Kaplan-Meier and time-dependent receiver operating characteristic (ROC) curves were generated to validate the prognostic performance of risk score in training datasets and validation datasets. Multivariable Cox regression analysis was used to identify independent prognostic factors in liver cancer. A prognostic nomogram was constructed based on The Cancer Genome Atlas (TCGA) dataset. Finally, the correlation between DNA methylation and prognosis-related genes was analyzed. Results:A twelve-gene signature including SPP1, KIF20A, HMMR, TPX2, LAPTM4B, TTK, MAGEA6, ANX10, LECT2, CYP2C9, RDH16 and LCAT was identified, and risk score was calculated by corresponding coefficients. The risk score model showed a strong diagnosis performance to distinguish HCC from normal samples. The HCC patients were stratified into high-risk and low-risk group based on the cutoff value of risk score. The Kaplan-Meier survival curves revealed significantly favorable overall survival in groups with lower risk score (P?<?0.0001). Time-dependent ROC analysis showed well prognostic performance of the twelve-gene signature, which was comparable or superior to AJCC stage at predicting 1-, 3-, and 5-year overall survival. In addition, the twelve-gene signature was independent with other clinical factors and performed better in predicting overall survival after combining with age and AJCC stage by nomogram. Moreover, most of the prognostic twelve genes were negatively correlated with DNA methylation in HCC tissues, which SPP1 and LCAT were identified as the DNA methylation-driven genes. Conclusions:We identified a twelve-gene signature as a robust marker with great potential for clinical application in risk stratification and overall survival prediction in HCC patients.

Project description:Hepatocellular carcinoma (HCC), the most frequently occurring type of cancer, is strongly associated with metabolic disorders. In this study, we aimed to characterize the metabolic features of HCC and normal tissue adjacent to the tumor (NAT). By using samples from The Cancer Genome Atlas (TCGA) liver cancer cohort and comparing 85 well-defined metabolic pathways obtained from the Kyoto Encyclopedia of Genes and Genomes (KEGG), 70 and 7 pathways were found to be significantly downregulated and upregulated, respectively, in HCC, revealing that tumor tissue lacks the ability to maintain normal metabolic levels. Through unsupervised hierarchical clustering of metabolic pathways, we found that metabolic heterogeneity correlated with prognosis in HCC samples. Thus, using the least absolute shrinkage and selection operator (LASSO) and filtering independent prognostic genes by the Cox proportional hazards model, a six-gene-based metabolic score model was constructed to enable HCC classification. This model showed that high expression of LDHA and CHAC2 was associated with an unfavorable prognosis but that high ADPGK, GOT2, MTHFS, and FTCD expression was associated with a favorable prognosis. Patients with higher metabolic scores had poor prognoses (p value = 2.19e-11, hazard ratio = 3.767, 95% CI = 2.555-5.555). By associating the score level with clinical features and genomic alterations, it was found that NAT had the lowest metabolic score and HCC with tumor stage III/IV the highest. qRT‒PCR results for HCC patients also revealed that tumor samples had higher score levels than NAT. Regarding genetic alterations, patients with higher metabolic scores had more TP53 gene mutations than those with lower metabolic scores (p value = 8.383e-05). Validation of this metabolic score model was performed using another two independent HCC cohorts from the Gene Expression Omnibus (GEO) repository and other TCGA datasets and achieved good performance, suggesting that this model may be used as a reliable tool for predicting the prognosis of HCC patients.