scholarly journals Identification of Immune-Related Genes for Establishment of Prognostic Index in Hepatocellular Carcinoma

Author(s):  
Yinfang Li ◽  
Ling Zou ◽  
Xuejun Liu ◽  
Judong Luo ◽  
Hui Liu

Background: Immune checkpoint inhibitor (ICI) therapy has been proved to be a promising therapy to many types of solid tumors. However, effective biomarker for estimating the response to ICI therapy and prognosis of hepatocellular carcinoma (HCC) patients remains underexplored. The aim of this study is to build a novel immune-related prognostic index based on transcriptomic profiles.Methods: Weighted gene co-expression network analysis (WGCNA) was conducted to identify immune-related hub genes that are differentially expressed in HCC cohorts. Next, univariate Cox regression analysis and least absolute shrinkage and selection operator (LASSO) analysis were used to detect hub genes associated to overall survival (OS). To validate the immune-related prognostic index, univariate and multivariate Cox regression analysis were performed. CIBERSORT and ESTIMATE were used to explore the tumor microenvironment and immune infiltration level.Results: The differential expression analysis detected a total of 148 immune-related genes, among which 25 genes were identified to be markedly related to overall survival in HCC patients. LASSO analysis yielded 10 genes used to construct the immune-related gene prognostic index (IRGPI), by which a risk score is computed to estimate low vs. high risk indicating the response to ICI therapy and prognosis. Further analysis confirmed that this immune-related prognostic index is an effective indicator to immune infiltration level, response to ICI treatment and OS. The IRGPI low-risk patients had better overall survival (OS) than IRGPI high-risk patients on two independent cohorts. Moreover, we found that IRGPI high-risk group was correlated with high TP53 mutation rate, immune-suppressing tumor microenvironment, and these patients acquired less benefit from ICI therapy. In contrast, IRGPI-low risk group was associated with low TP53 and PIK3CA mutation rate, high infiltration of naive B cells and T cells, and these patients gained relatively more benefit from ICI therapy.

2020 ◽  
Author(s):  
Li Liu ◽  
She Tian ◽  
Zhu Li ◽  
Yongjun Gong ◽  
Hao Zhang

Abstract Background : Hepatocellular carcinoma (HCC) is one of the most common clinical malignant tumors, resulting in high mortality and poor prognosis. Studies have found that LncRNA plays an important role in the onset, metastasis and recurrence of hepatocellular carcinoma. The immune system plays a vital role in the development, progression, metastasis and recurrence of cancer. Therefore, immune-related lncRNA can be used as a novel biomarker to predict the prognosis of hepatocellular carcinoma. Methods : The transcriptome data and clinical data of HCC patients were obtained by using The Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA‑LIHC), and immune-related genes were extracted from the Molecular Signatures Database (IMMUNE RESPONSE M19817 and IMMUNE SYSTEM PROCESS M13664). By constructing the co-expression network and Cox regression analysis, 13 immune-lncRNAs was identified to predict the prognosis of HCC patients. Patients were divided into high risk group and low risk group by using the risk score formula, and the difference in overall survival (OS) between the two groups was reflected by Kaplan-Meier survival curve. The time - dependent receiver operating characteristics (ROC) analysis and principal component analysis (PCA) were used to evaluate 13 immune -lncRNAs signature. Results : Through TCGA - LIHC extracted from 343 cases of patients with hepatocellular carcinoma RNA - Seq data and clinical data, 331 immune-related genes were extracted from the Molecular Signatures Database , co-expression networks and Cox regression analysis were constructed, 13 immune-lncRNAs signature was identified as biomarkers to predict the prognosis of patients. At the same time using the risk score median divided the patients into high risk and low risk groups, and through the Kaplan-Meier survival curve analysis found that high-risk group of patients' overall survival (OS) less low risk group of patients. The AUC value of the ROC curve is 0.828, and principal component analysis (PCA) results showed that patients could be clearly divided into two parts by immune-lncRNAs, which provided evidence for the use of 13 immune-lncRNAs signature as prognostic markers. Conclusion : Our study identified 13 immune-lncRNAs signature that can effectively predict the prognosis of HCC patients, which may be a new prognostic indicator for predicting clinical outcomes.


2020 ◽  
Author(s):  
YuPing Bai ◽  
Wenbo Qi ◽  
Le Liu ◽  
Jing Zhang ◽  
Lan Pang ◽  
...  

Abstract Background: Hepatocellular carcinoma is ranked fifth among the most common cancer worldwide. Hypoxia can induce tumor growth, but the relationship with HCC prognosis remains unclear. Our study aims to construct a hypoxia-related multigene model to predict the prognosis of HCC. Methods: RNA-seq expression data and related clinical information were download from TCGA database and ICGC database, respectively. Univariate/multivariate Cox regression analysis was used to construct prognostic models. KM curve analysis, and ROC curve were used to evaluate the prognostic models, which were further verified in the clinical traits and ICGC database. GSEA analyzed pathway enrichment in high-risk groups. Nomogram was constructed to predict the personalized treatment of patients. Finally, real-time fluorescence quantitative PCR(RT-qPCR) was used to detect the expressions of KDELR3 and SCARB1 in normal hepatocytes and 4 hepatocellular carcinoma cells. Results: Through a series of analyses, 7 prognostic markers related to HCC survival were constructed. HCC patients were divided into the high and low risk group, and the results of KM curve showed that there was a significant difference between the two groups. Stratified analysis,found that there were significant differences in risk values of different ages, genders, stages and grades, which could be used as independent predictors. In addition, we assessed the risk value in the clinical traits analysis and found that it could accelerate the progression of cancer, while the results of GSEA enrichment analysis showed that the high-risk group patients were mainly distributed in the cell cycle and other pathways. Then, Nomogram was constructed to predict the overall survival of patients. Finally, RT-qPCR showed that KDELR3 and SCARB1 were highly expressed in HepG2 and L02, respectively. Conclusion: This study provides a potential diagnostic indicator for HCC patients, and help clinicians to deepen the comprehension in HCC pathogenesis so as to make personalized medical decisions.


2021 ◽  
Author(s):  
Dandong Luo ◽  
Qiang Tao ◽  
HuiJuan Jiang ◽  
JunLing Zhu ◽  
Ning Zhang ◽  
...  

Abstract Background Hepatocellular carcinoma (HCC) is characterized by widespread epidemiology and extraordinary heterogeneity, with challenging prognosis prediction. Ferroptosis is a regulatory cell death driven by iron-dependent lipid peroxidation. The main aim of this study was to determine the predictive value of ferroptosis-related genes (FRGs) in HCC. Methods Herein, the data of HCC patients were downloaded from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) public databases. In the ICGC cohort, a multigenic signature was constructed using the LASSO Cox regression model. Next, patients in the TCGA cohort were used to verify the reliability of the model. Results Results showed that 30.07% of the differentially expressed genes (DEGs) in the ICGC cohort were associated with ferroptosis. Among them, 35 genes were identified as intersected genes associated with overall survival in both cohorts. Moreover, an 8-gene signature for prediction of HCC patients was constructed and the patients were divided it into low-risk and high-risk groups. The results indicated that the overall survival (OS) of patients in the high-risk group was lower than OS of patients in the low-risk group (P < 0.001 in both cohorts). Multivariate Cox regression analysis indicated that the risk score was an independent predictor of OS (HR > 1, P < 0.001). Receiver operating curves (ROC) demonstrated the predictive power of the signature. Furthermore, functional enrichment analysis revealed the existence of significantly correlated immune-related pathways, and their immune states were different between groups. Conclusions In summary, the genetic signature described in this study was associated with ferroptosis and it can be used to predict the prognosis of HCC. Therefore, targeted treatment of ferroptosis may be an alternative treatment option for HCC.


2021 ◽  
Vol 12 ◽  
Author(s):  
Yuping Bai ◽  
Wenbo Qi ◽  
Le Liu ◽  
Jing Zhang ◽  
Lan Pang ◽  
...  

BackgroundHepatocellular carcinoma (HCC) is ranked fifth among the most common cancer worldwide. Hypoxia can induce tumor growth, but the relationship with HCC prognosis remains unclear. Our study aims to construct a hypoxia-related multigene model to predict the prognosis of HCC.MethodsRNA-seq expression data and related clinical information were download from TCGA database and ICGC database, respectively. Univariate/multivariate Cox regression analysis was used to construct prognostic models. KM curve analysis, and ROC curve were used to evaluate the prognostic models, which were further verified in the clinical traits and ICGC database. GSEA analyzed pathway enrichment in high-risk groups. Nomogram was constructed to predict the personalized treatment of patients. Finally, real-time fluorescence quantitative PCR (RT-qPCR) was used to detect the expressions of KDELR3 and SCARB1 in normal hepatocytes and 4 HCC cells. The expressions of SCARB1 in hepatocellular carcinoma tissue in 46 patients were detected by immunohistochemistry, and the correlation between its expressions and disease free survival of patient was calculated.ResultsThrough a series of analyses, seven prognostic markers related to HCC survival were constructed. HCC patients were divided into the high and low risk group, and the results of KM curve showed that there was a significant difference between the two groups. Stratified analysis, found that there were significant differences in risk values of different ages, genders, stages and grades, which could be used as independent predictors. In addition, we assessed the risk value in the clinical traits analysis and found that it could accelerate the progression of cancer, while the results of GSEA enrichment analysis showed that the high-risk group patients were mainly distributed in the cell cycle and other pathways. Then, Nomogram was constructed to predict the overall survival of patients. Finally, RT-qPCR showed that KDELR3 and SCARB1 were highly expressed in HepG2 and L02, respectively. Results of IHC staining showed that SCARB1 was highly expressed in cancer tissues compared to adjacent normal liver tissues and its expression was related to hepatocellular carcinoma differentiation status. The Kaplan-Meier survival showed a poor percent survival in the SCARB1 high group compared to that in the SCARB1 low group.ConclusionThis study provides a potential diagnostic indicator for HCC patients, and help clinicians to deepen the comprehension in HCC pathogenesis so as to make personalized medical decisions.


2020 ◽  
Author(s):  
Ye Liu ◽  
Zhixiang Qin ◽  
Hai Yang ◽  
Yang Gu ◽  
Kun Li

Abstract Background Hepatocellular carcinoma (HCC) represents one of the deadliest malignancies worldwide. Despite significant advances in diagnosis and treatment, the mortality rate from HCC persists at a substantial level. This research strives to establish a prognostic model based on the RNA binding proteins (RBPs) that can predict HCC patients’ OS. Methods There was an RNA-seq data set derived from the Cancer Genome Atlas (TCGA) databank which was included in our research as well as a Microarray data set (GSE14520). The differentially expressed RBPs between HCC and normal tissues were investigated in TCGA dataset. Subsequently, the TCGA data set was randomly split into a training and a testing cohort. The prognostic model of the training cohort was developed by applying univariate Cox regression and lasso Cox regression analyses and multivariate Cox regression analysis. In order to evaluate the prognostic value of the model, a comprehensive survival assessment was conducted. Results A total of 133 differentially expressed RBPs were identified. Five RBPs (RPL10L, EZH2, PPARGC1A, ZNF239, IFIT1) were used to construct the model. The model accurately predicted the prognosis of liver cancer patients in both the TCGA cohort and the GSE14520 validation cohort. HCC patients could be assigned into a high-risk group and a low-risk group by this model, and the overall survival of these two groups was significantly different. Furthermore, the risk scores obtained by our model were highly correlated with immune cell infiltration. . Conclusions Five RBPs-related prognostic models were constructed and validated to predict OS reliably in HCC individuals. It helps to identify patients at high risk of mortality with the risk prediction score, which optimizes personalized therapeutic decision-making.


PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e11273
Author(s):  
Lei Yang ◽  
Weilong Yin ◽  
Xuechen Liu ◽  
Fangcun Li ◽  
Li Ma ◽  
...  

Background Hepatocellular carcinoma (HCC) is considered to be a malignant tumor with a high incidence and a high mortality. Accurate prognostic models are urgently needed. The present study was aimed at screening the critical genes for prognosis of HCC. Methods The GSE25097, GSE14520, GSE36376 and GSE76427 datasets were obtained from Gene Expression Omnibus (GEO). We used GEO2R to screen differentially expressed genes (DEGs). A protein-protein interaction network of the DEGs was constructed by Cytoscape in order to find hub genes by module analysis. The Metascape was performed to discover biological functions and pathway enrichment of DEGs. MCODE components were calculated to construct a module complex of DEGs. Then, gene set enrichment analysis (GSEA) was used for gene enrichment analysis. ONCOMINE was employed to assess the mRNA expression levels of key genes in HCC, and the survival analysis was conducted using the array from The Cancer Genome Atlas (TCGA) of HCC. Then, the LASSO Cox regression model was performed to establish and identify the prognostic gene signature. We validated the prognostic value of the gene signature in the TCGA cohort. Results We screened out 10 hub genes which were all up-regulated in HCC tissue. They mainly enrich in mitotic cell cycle process. The GSEA results showed that these data sets had good enrichment score and significance in the cell cycle pathway. Each candidate gene may be an indicator of prognostic factors in the development of HCC. However, hub genes expression was weekly associated with overall survival in HCC patients. LASSO Cox regression analysis validated a five-gene signature (including CDC20, CCNB2, NCAPG, ASPM and NUSAP1). These results suggest that five-gene signature model may provide clues for clinical prognostic biomarker of HCC.


2021 ◽  
Vol 20 ◽  
pp. 153303382110414
Author(s):  
Xiaoyong Li ◽  
Jiaqong Lin ◽  
Yuguo pan ◽  
Peng Cui ◽  
Jintang Xia

Background: Liver progenitor cells (LPCs) play significant roles in the development and progression of hepatocellular carcinoma (HCC). However, no studies on the value of LPC-related genes for evaluating HCC prognosis exist. We developed a gene signature of LPC-related genes for prognostication in HCC. Methods: To identify LPC-related genes, we analyzed mRNA expression arrays from a dataset (GSE57812 & GSE 37071) containing LPCs, mature hepatocytes, and embryonic stem cell samples. HCC RNA-Seq data from The Cancer Genome Atlas (TCGA) were used to explore the differentially expressed genes (DEGs) related to prognosis through DEG analysis and univariate Cox regression analysis. Lasso and multivariate Cox regression analyses were performed to construct the LPC-related gene prognostic model in the TCGA training dataset. This model was validated in the TCGA testing set and an external dataset (International Cancer Genome Consortium [ICGC] dataset). Finally, we investigated the relationship between this prognostic model with tumor-node-metastasis stage, tumor grade, and vascular invasion of HCC. Results: Overall, 1770 genes were identified as LPC-related genes, of which 92 genes were identified as DEGs in HCC tissues compared with normal tissues. Furthermore, we randomly assigned patients from the TCGA dataset to the training and testing cohorts. Twenty-six DEGs correlated with overall survival (OS) in the univariate Cox regression analysis. Lasso and multivariate Cox regression analyses were performed in the TCGA training set, and a 3-gene signature was constructed to stratify patients into 2 risk groups: high-risk and low-risk. Patients in the high-risk group had significantly lower OS than those in the low-risk group. Receiver operating characteristic curve analysis confirmed the signature's predictive capacity. Moreover, the risk score was confirmed to be an independent predictor for patients with HCC. Conclusion: We demonstrated that the LPC-related gene signature can be used for prognostication in HCC. Thus, targeting LPCs may serve as a therapeutic alternative for HCC.


2021 ◽  
Vol 11 ◽  
Author(s):  
Fen Liu ◽  
Zongcheng Yang ◽  
Lixin Zheng ◽  
Wei Shao ◽  
Xiujie Cui ◽  
...  

BackgroundGastric cancer is a common gastrointestinal malignancy. Since it is often diagnosed in the advanced stage, its mortality rate is high. Traditional therapies (such as continuous chemotherapy) are not satisfactory for advanced gastric cancer, but immunotherapy has shown great therapeutic potential. Gastric cancer has high molecular and phenotypic heterogeneity. New strategies for accurate prognostic evaluation and patient selection for immunotherapy are urgently needed.MethodsWeighted gene coexpression network analysis (WGCNA) was used to identify hub genes related to gastric cancer progression. Based on the hub genes, the samples were divided into two subtypes by consensus clustering analysis. After obtaining the differentially expressed genes between the subtypes, a gastric cancer risk model was constructed through univariate Cox regression, least absolute shrinkage and selection operator (LASSO) regression and multivariate Cox regression analysis. The differences in prognosis, clinical features, tumor microenvironment (TME) components and immune characteristics were compared between subtypes and risk groups, and the connectivity map (CMap) database was applied to identify potential treatments for high-risk patients.ResultsWGCNA and screening revealed nine hub genes closely related to gastric cancer progression. Unsupervised clustering according to hub gene expression grouped gastric cancer patients into two subtypes related to disease progression, and these patients showed significant differences in prognoses, TME immune and stromal scores, and suppressive immune checkpoint expression. Based on the different expression patterns between the subtypes, we constructed a gastric cancer risk model and divided patients into a high-risk group and a low-risk group based on the risk score. High-risk patients had a poorer prognosis, higher TME immune/stromal scores, higher inhibitory immune checkpoint expression, and more immune characteristics suitable for immunotherapy. Multivariate Cox regression analysis including the age, stage and risk score indicated that the risk score can be used as an independent prognostic factor for gastric cancer. On the basis of the risk score, we constructed a nomogram that relatively accurately predicts gastric cancer patient prognoses and screened potential drugs for high-risk patients.ConclusionsOur results suggest that the 7-gene signature related to tumor progression could predict the clinical prognosis and tumor immune characteristics of gastric cancer.


2021 ◽  
Vol 12 ◽  
Author(s):  
Susu Zheng ◽  
Xiaoying Xie ◽  
Xinkun Guo ◽  
Yanfang Wu ◽  
Guobin Chen ◽  
...  

Pyroptosis is a novel kind of cellular necrosis and shown to be involved in cancer progression. However, the diverse expression, prognosis and associations with immune status of pyroptosis-related genes in Hepatocellular carcinoma (HCC) have yet to be analyzed. Herein, the expression profiles and corresponding clinical characteristics of HCC samples were collected from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Then a pyroptosis-related gene signature was built by applying the least absolute shrinkage and selection operator (LASSO) Cox regression model from the TCGA cohort, while the GEO datasets were applied for verification. Twenty-four pyroptosis-related genes were found to be differentially expressed between HCC and normal samples. A five pyroptosis-related gene signature (GSDME, CASP8, SCAF11, NOD2, CASP6) was constructed according to LASSO Cox regression model. Patients in the low-risk group had better survival rates than those in the high-risk group. The risk score was proved to be an independent prognostic factor for overall survival (OS). The risk score correlated with immune infiltrations and immunotherapy responses. GSEA indicated that endocytosis, ubiquitin mediated proteolysis and regulation of autophagy were enriched in the high-risk group, while drug metabolism cytochrome P450 and tryptophan metabolism were enriched in the low-risk group. In conclusion, our pyroptosis-related gene signature can be used for survival prediction and may also predict the response of immunotherapy.


2021 ◽  
Vol 7 ◽  
Author(s):  
Xiaoyu Deng ◽  
Qinghua Bi ◽  
Shihan Chen ◽  
Xianhua Chen ◽  
Shuhui Li ◽  
...  

Although great progresses have been made in the diagnosis and treatment of hepatocellular carcinoma (HCC), its prognostic marker remains controversial. In this current study, weighted correlation network analysis and Cox regression analysis showed significant prognostic value of five autophagy-related long non-coding RNAs (AR-lncRNAs) (including TMCC1-AS1, PLBD1-AS1, MKLN1-AS, LINC01063, and CYTOR) for HCC patients from data in The Cancer Genome Atlas. By using them, we constructed a five-AR-lncRNA prognostic signature, which accurately distinguished the high- and low-risk groups of HCC patients. All of the five AR lncRNAs were highly expressed in the high-risk group of HCC patients. This five-AR-lncRNA prognostic signature showed good area under the curve (AUC) value (AUC = 0.751) for the overall survival (OS) prediction in either all HCC patients or HCC patients stratified according to several clinical traits. A prognostic nomogram with this five-AR-lncRNA signature predicted the 3- and 5-year OS outcomes of HCC patients intuitively and accurately (concordance index = 0.745). By parallel comparison, this five-AR-lncRNA signature has better prognosis accuracy than the other three recently published signatures. Furthermore, we discovered the prediction ability of the signature on therapeutic outcomes of HCC patients, including chemotherapy and immunotherapeutic responses. Gene set enrichment analysis and gene mutation analysis revealed that dysregulated cell cycle pathway, purine metabolism, and TP53 mutation may play an important role in determining the OS outcomes of HCC patients in the high-risk group. Collectively, our study suggests a new five-AR-lncRNA prognostic signature for HCC patients.


Sign in / Sign up

Export Citation Format

Share Document