scholarly journals A prognostic model for hepatocellular carcinoma based on apoptosis-related genes

Author(s):  
Renjie Liu ◽  
Guifu Wang ◽  
Chi Zhang ◽  
Dousheng Bai

Abstract Background: Dysregulation of the balance between proliferation and apoptosis is the basis for human hepatocarcinogenesis. In many malignant tumors, such as hepatocellular carcinoma (HCC), there is a correlation between apoptotic dysregulation and poor prognosis. However, the prognostic values of apoptosis-related genes (ARGs) in HCC have not been elucidated. Methods: To screen for differentially expressed ARGs, the expression levels of 161 ARGs from The Cancer Genome Atlas (TCGA) database(https://cancergenome.nih.gov/) were analyzed. Gene Ontology (GO) enrichment and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to evaluate the underlying molecular mechanisms of differentially expressed ARGs in HCC. The prognostic values of ARGs were established using Cox regression, and subsequently, a prognostic risk model for scoring patients was developed. Kaplan-Meier (K-M) and receiver operating characteristic (ROC) curves were plotted to determine the prognostic value of the model. Results: Compared to normal tissues, 43 highly up-regulated and 8 down-regulated ARGs in HCC tissues were screened. GO analysis results revealed that these 51 genes are indeed related to the apoptosis function. KEGG analysis revealed that these 51 genes were correlated with MAPK, P53, TNF, and PI3K-AKT signaling pathways, while Cox regression revealed that 5 ARGs (PPP2R5B, SQSTM1, TOP2A, BMF, and LGALS3) were associated with prognosis and were, therefore, obtained to develop the prognostic model. Based on the median risk scores, patients were categorized into high-risk and low-risk groups. Patients in the low-risk groups exhibited significantly elevated two-year or five-year survival probabilities (p < 0.0001). The risk model had a better clinical potency than the other clinical characteristics, with the area under the ROC curve (AUC = 0.741). The prognosis of HCC patients was established from a plotted nomogram. Conclusion: Based on the differential expression of ARGs, we established a novel risk model for predicting HCC prognosis. This model can also be used to inform the individualized treatment of HCC patients.

2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Renjie Liu ◽  
Guifu Wang ◽  
Chi Zhang ◽  
Dousheng Bai

Abstract Background Dysregulation of the balance between proliferation and apoptosis is the basis for human hepatocarcinogenesis. In many malignant tumors, such as hepatocellular carcinoma (HCC), there is a correlation between apoptotic dysregulation and poor prognosis. However, the prognostic values of apoptosis-related genes (ARGs) in HCC have not been elucidated. Methods To screen for differentially expressed ARGs, the expression levels of 161 ARGs from The Cancer Genome Atlas (TCGA) database (https://cancergenome.nih.gov/) were analyzed. Gene Ontology (GO) enrichment and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to evaluate the underlying molecular mechanisms of differentially expressed ARGs in HCC. The prognostic values of ARGs were established using Cox regression, and subsequently, a prognostic risk model for scoring patients was developed. Kaplan–Meier (K-M) and receiver operating characteristic (ROC) curves were plotted to determine the prognostic value of the model. Results Compared with normal tissues, 43 highly upregulated and 8 downregulated ARGs in HCC tissues were screened. GO analysis results revealed that these 51 genes are indeed related to the apoptosis function. KEGG analysis revealed that these 51 genes were correlated with MAPK, P53, TNF, and PI3K-AKT signaling pathways, while Cox regression revealed that 5 ARGs (PPP2R5B, SQSTM1, TOP2A, BMF, and LGALS3) were associated with prognosis and were, therefore, obtained to develop the prognostic model. Based on the median risk scores, patients were categorized into high-risk and low-risk groups. Patients in the low-risk groups exhibited significantly elevated 2-year or 5-year survival probabilities (p < 0.0001). The risk model had a better clinical potency than the other clinical characteristics, with the area under the ROC curve (AUC = 0.741). The prognosis of HCC patients was established from a plotted nomogram. Conclusion Based on the differential expression of ARGs, we established a novel risk model for predicting HCC prognosis. This model can also be used to inform the individualized treatment of HCC patients.


2021 ◽  
Author(s):  
Renjie Liu ◽  
Guifu Wang ◽  
Chi Zhang ◽  
Dousheng Bai

Abstract Background Dysregulation of the balance between proliferation and apoptosis is the basis in human hepatocarcinogenesis. There is a possible association of apoptosis dysregulation with poor prognosis in many malignant tumors, such as hepatocellular carcinoma (HCC). However, the prognostic effect of Apoptosis-related genes (ARGs) on HCC is still unclear. Methods A total of 161 ARGs expression levels were analyzed based on The Cancer Genome Atlas (TCGA) database(https://cancergenome.nih.gov/) to screen for differentially expressed ARGs. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to determine the underlying molecular mechanisms of screened ARGs in HCC. ARGs prognostic values were identified using Cox regression to subsequently establish a prognostic risk model and scoring in patients. Kaplan-Meier (K-M) and receiver operating characteristic (ROC) curves were plotted to determine the prognostic value in the model. Results Compared to normal specimens, 43 highly upregulated and 8 downregulated ARGs respectively and their normal counterparts in HCC specimens were screened. KEGG analysis demonstrated pathways correlated with these 51 genes which included MAPK, P53, TNF, PI3K-Akt signaling pathways. With Cox regression, 5 prognostic correlated with ARGs (PPP2R5B, SQSTM1, TOP2A, BMF, and LGALS3) were obtained to develop the prognosis model. According to the median of risk scores, patients were categorized into high-risk and low-risk groups. Patients in low-risk groups had a significantly higher two-year or five-year survival probability (p < 0.0001). The risk model had better potency than other clinical characteristics, with the area under the ROC curve (AUC = 0.741). Prognosis of HCC patients was established from a plotted nomogram. Conclusion This present study established a novel prognostic risk model for predicting HCC according to the expression of ARGs. The present advancement can potentially contribute to prediction prognosis and individualized treatment of HCC patients.


2021 ◽  
Vol 2021 ◽  
pp. 1-15
Author(s):  
Hao Guo ◽  
Jing Zhou ◽  
Yanjun Zhang ◽  
Zhi Wang ◽  
Likun Liu ◽  
...  

Background. Hypoxia closely relates to malignant progression and appears to be prognostic for outcome in hepatocellular carcinoma (HCC). Our research is aimed at mining the hypoxic-related genes (HRGs) and constructing a prognostic predictor (PP) model on clinical prognosis in HCC patients. Methods. RNA-sequencing data about HRGs and clinical data of patients with HCC were obtained from The Cancer Genome Atlas (TCGA) database portal. Differentially expressed HRGs between HCC and para-carcinoma tissue samples were obtained by applying the Wilcox analysis in R statistical software. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes were used for gene functional enrichment analyses. Then, the patients who were asked to follow up for at least one month were enrolled in the following study. Cox proportional risk regression model was applied to obtain key HRGs which related to overall survival (OS) in HCC. PP was constructed and defined, and the accuracy of PP was validated by constructing the signature in a training set and validation set. Connectivity map (CMap) was used to find potential drugs, and gene set cancer analysis (GSCA) was also performed to explore the underlying molecular mechanisms. Results. Thirty-seven differentially expressed HRGs were obtained. It contained 28 upregulated and 9 downregulated genes. After the univariate Cox regression model analysis, we obtained 27 prognosis-related HRGs. Of these, 25 genes were risk factors for cancer, and 2 genes were protective factors. The PP was composed by 12 key genes (HDLBP, SAP30, PFKP, DPYSL4, SLC2A1, HMOX1, PGK1, ERO1A, LDHA, ENO2, SLC6A6, and TPI1). GSCA results showed the overall activity of these 12 key genes in 10 cancer-related pathways. Besides, CMap identified deferoxamine, crotamiton, talampicillin, and lycorine might have effects with HCC. Conclusions. This study firstly reported 12 prognostic HRGs and constructed the model of the PP. This comprehensive research of multiple databases helps us gain insight into the biological properties of HCC and provides deferoxamine, crotamiton, talampicillin, and lycorine as potential drugs to fight against HCC.


2020 ◽  
Author(s):  
Xiang Zhou ◽  
Keying Zhang ◽  
Fa Yang ◽  
Chao Xu ◽  
Jianhua Jiao ◽  
...  

Abstract Background: Hepatocellular carcinoma (HCC) is a disease with higher morbidity, mortality, and poor prognosis in the whole world. Understanding the crosslink between HCC and the immune system is essential for people to uncover a few potential and valuable therapeutic strategies. This study aimed to reveal the correlation between HCC and immune-related genes and establish a clinical evaluation model. Methods: We had analyzed the clinical information consisted of 373 HCC and 49 normal samples from the cancer genome atlas (TCGA). The differentially expressed genes (DEGs) were selected by the Wilcoxon test and the immune-related differentially expressed genes (IRDEGs) in DEGs were identified by matching DEGs with immune-related genes downloaded from the ImmPort database. Furthermore, the univariate Cox regression analysis and multivariate Cox regression analysis were performed to construct a prognostic risk model. Then, twenty-two types of tumor immune-infiltrating cells (TIICs) were downloaded from Tumor Immune Estimation Resource (TIMER) and were used to construct the correlational graphs between the TIICs and risk score by the CIBERSORT. Subsequently, the transcription factors (TFs) were gained in the Cistrome website and the differentially expressed TFs (DETFs) were achieved. Finally, the KEGG pathway analysis and GO analysis were performed to further understand the molecular mechanisms between DETFs and PDIRGs.Results: In our study, 5839 DEGs, 326 IRDEGs, and 31 prognosis-related IRDEGs (PIRDEGs) were identified. And 8 optimal PIRDEGs were employed to construct a prognostic risk model by multivariate Cox regression analysis. The correlation between risk genes and clinical characterizations and TIICs has verified that the prognostic model was effective in predicting the prognosis of HCC patients. Finally, several important immune-related pathways and molecular functions of the eight PIRDEGs were significantly enriched and there was a distinct association between the risk IRDEGs and TFs. Conclusion: The prognostic risk model showed a more valuable predicting role for HCC patients, and produced many novel therapeutic targets and strategies for HCC.


BMC Cancer ◽  
2020 ◽  
Vol 20 (1) ◽  
Author(s):  
Min Wang ◽  
Shan Huang ◽  
Zefeng Chen ◽  
Zhiwei Han ◽  
Kezhi Li ◽  
...  

Abstract Background Hepatocellular carcinoma (HCC) is among the deadliest forms of cancer. While RNA-binding proteins (RBPs) have been shown to be key regulators of oncogenesis and tumor progression, their dysregulation in the context of HCC remains to be fully characterized. Methods Data from the Cancer Genome Atlas - liver HCC (TCGA-LIHC) database were downloaded and analyzed in order to identify RBPs that were differentially expressed in HCC tumors relative to healthy normal tissues. Functional enrichment analyses of these RBPs were then conducted using the GO and KEGG databases to understand their mechanistic roles. Central hub RBPs associated with HCC patient prognosis were then detected through Cox regression analyses, and were incorporated into a prognostic model. The prognostic value of this model was then assessed through the use of Kaplan-Meier curves, time-related ROC analyses, univariate and multivariate Cox regression analyses, and nomograms. Lastly, the relationship between individual hub RBPs and HCC patient overall survival (OS) was evaluated using Kaplan-Meier curves. Finally, find protein-coding genes (PCGs) related to hub RBPs were used to construct a hub RBP-PCG co-expression network. Results In total, we identified 81 RBPs that were differentially expressed in HCC tumors relative to healthy tissues (54 upregulated, 27 downregulated). Seven prognostically-relevant hub RBPs (SMG5, BOP1, LIN28B, RNF17, ANG, LARP1B, and NR0B1) were then used to generate a prognostic model, after which HCC patients were separated into high- and low-risk groups based upon resultant risk score values. In both the training and test datasets, we found that high-risk HCC patients exhibited decreased OS relative to low-risk patients, with time-dependent area under the ROC curve values of 0.801 and 0.676, respectively. This model thus exhibited good prognostic performance. We additionally generated a prognostic nomogram based upon these seven hub RBPs and found that four other genes were significantly correlated with OS. Conclusion We herein identified a seven RBP signature that can reliably be used to predict HCC patient OS, underscoring the prognostic relevance of these genes.


Author(s):  
Qianqian Wu ◽  
Sutian Jiang ◽  
Tong Cheng ◽  
Manyu Xu ◽  
Bing Lu

Hepatocellular carcinoma (HCC) is the second most lethal malignant tumor because of its significant heterogeneity and complicated molecular pathogenesis. Novel prognostic biomarkers are urgently needed because no effective and reliable prognostic biomarkers currently exist for HCC patients. Increasing evidence has revealed that pyroptosis plays a role in the occurrence and progression of malignant tumors. However, the relationship between pyroptosis-related genes (PRGs) and HCC patient prognosis remains unclear. In this study, 57 PRGs were obtained from previous studies and GeneCards. The gene expression profiles and clinical data of HCC patients were acquired from public data portals. Least absolute shrinkage and selection operator (LASSO) Cox regression analysis was performed to establish a risk model using TCGA data. Additionally, the risk model was further validated in an independent ICGC dataset. Our results showed that 39 PRGs were significantly differentially expressed between tumor and normal liver tissues in the TCGA cohort. Functional analysis confirmed that these PRGs were enriched in pyroptosis-related pathways. According to univariate Cox regression analysis, 14 differentially expressed PRGs were correlated with the prognosis of HCC patients in the TCGA cohort. A risk model integrating two PRGs was constructed to classify the patients into different risk groups. Poor overall survival was observed in the high-risk group of both TCGA (p &lt; 0.001) and ICGC (p &lt; 0.001) patients. Receiver operating characteristic curves demonstrated the accuracy of the model. Furthermore, the risk score was confirmed as an independent prognostic indicator via multivariate Cox regression analysis (TCGA cohort: HR = 3.346, p &lt; 0.001; ICGC cohort: HR = 3.699, p &lt; 0.001). Moreover, the single-sample gene set enrichment analysis revealed different immune statuses between high- and low-risk groups. In conclusion, our new pyroptosis-related risk model has potential application in predicting the prognosis of HCC patients.


2021 ◽  
Author(s):  
Ju Kun Wang ◽  
Ke Han ◽  
Chao Zhang ◽  
Xin Chen ◽  
Yu Li ◽  
...  

Purpose: ADME genes are genes involved in drug absorption, distribution, metabolism, and excretion (ADME). Previous studies report that expression levels of ADME-related genes correlate with prognosis of hepatocellular carcinoma (HCC) patients. However, the role of ADME gene expression on HCC prognosis has not been fully explored. This study sought to construct a prediction model using ADME-related genes for prognosis of HCC. Methods: Transcriptome and clinical data were retrieved from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC), which were used as training and validation cohorts, respectively. A prediction model was constructed using univariate Cox regression and LASSO analysis. Patients were divided into high- and low-risk groups based on the median risk score. The predictive ability of the risk signature was estimated through bioinformatics analyses. Results: Six ADME-related genes (CYP2C9, ABCB6, ABCC5, ADH4, DHRS13, and SLCO2A1) were used to construct the prediction model with a good predictive ability. Univariate and multivariate Cox regression analyses showed the risk signature was an independent predictor of overall survival. A single-sample gene set enrichment analysis (ssGSEA) strategy showed a significant relationship between risk signature and immune status. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed differentially expressed genes in the high- and low-risk groups were enriched in biological process associated with metabolic and cell cycle pathways. Conclusion: A prediction model was constructed using six ADME-related genes for prediction of HCC prognosis. This signature can be used to improve HCC diagnosis, treatment, and prognosis in clinical use.


2020 ◽  
Author(s):  
Min wang ◽  
Shan Huang ◽  
Zefeng Chen ◽  
Zhiwei Han ◽  
Kezhi Li ◽  
...  

Abstract Background: Hepatocellular carcinoma (HCC) is among the deadliest forms of cancer. While RNA-binding proteins (RBPs) have been shown to be key regulators of oncogenesis and tumor progression, their dysregulation in the context of HCC remains to be fully characterized. Methods: Data from the Cancer Genome Atlas - liver HCC (TCGA-LIHC) database were downloaded and analyzed in order to identify RBPs that were differentially expressed in HCC tumors relative to healthy normal tissues. Functional enrichment analyses of these RBPs were then conducted using the GO and KEGG databases to understand their mechanistic roles. Central hub RBPs associated with HCC patient prognosis were then detected through Cox regression analyses, and were incorporated into a prognostic model. The prognostic value of this model was then assessed through the use of Kaplan-Meier curves, time-related ROC analyses, univariate and multivariate Cox regression analyses, and nomograms. Lastly, the relationship between individual hub RBPs and HCC patient overall survival (OS) was evaluated using Kaplan-Meier curves. Results: In total, we identified 81 RBPs that were differentially expressed in HCC tumors relative to healthy tissues (54 upregulated, 27 downregulated). Seven prognostically-relevant hub RBPs (SMG5, BOP1, LIN28B, RNF17, ANG, LARP1B, and NR0B1) were then used to generate a prognostic model, after which HCC patients were separated into high- and low-risk groups based upon resultant risk score values. In both the training and test datasets, we found that high-risk HCC patients exhibited decreased OS relative to low-risk patients, with time-dependent area under the ROC curve values of 0.801 and 0.676, respectively. This model thus exhibited good prognostic performance. We additionally generated a prognostic nomogram based upon these seven hub RBPs and found that four other genes were significantly correlated with OS. Conclusion: We herein identified a seven RBP signature that can reliably be used to predict HCC patient OS, underscoring the prognostic relevance of these genes.


Author(s):  
Yue Li ◽  
Ruoyi Shen ◽  
Anqi Wang ◽  
Jian Zhao ◽  
Jieqi Zhou ◽  
...  

BackgroundLung adenocarcinoma (LUAD) originates mainly from the mucous epithelium and glandular epithelium of the bronchi. It is the most common pathologic subtype of non-small cell lung cancer (NSCLC). At present, there is still a lack of clear criteria to predict the efficacy of immunotherapy. The 5-year survival rate for LUAD patients remains low.MethodsAll data were downloaded from The Cancer Genome Atlas (TCGA) database. We used Gene Set Enrichment Analysis (GSEA) database to obtain immune-related mRNAs. Immune-related lncRNAs were acquired by using the correlation test of the immune-related genes with R version 3.6.3 (Pearson correlation coefficient cor = 0.5, P &lt; 0.05). The TCGA-LUAD dataset was divided into the testing set and the training set randomly. Based on the training set to perform univariate and multivariate Cox regression analyses, we screened prognostic immune-related lncRNAs and given a risk score to each sample. Samples were divided into the high-risk group and the low-risk group according to the median risk score. By the combination of Kaplan–Meier (KM) survival curve, the receiver operating characteristic (ROC) (AUC) curve, the independent risk factor analysis, and the clinical data of the samples, we assessed the accuracy of the risk model. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed on the differentially expressed mRNAs between the high-risk group and the low-risk group. The differentially expressed genes related to immune response between two risk groups were analyzed to evaluate the role of the model in predicting the efficacy and effects of immunotherapy. In order to explain the internal mechanism of the risk model in predicting the efficacy of immunotherapy, we analyzed the differentially expressed genes related to epithelial-mesenchymal transition (EMT) between two risk groups. We extracted RNA from normal bronchial epithelial cell and LUAD cells and verified the expression level of lncRNAs in the risk model by a quantitative real-time polymerase chain reaction (qRT-PCR) test. We compared our risk model with other published prognostic signatures with data from an independent cohort. We transfected LUAD cell with siRNA-LINC0253. Western blot analysis was performed to observed change of EMT-related marker in protein level.ResultsThrough univariate Cox regression analysis, 24 immune-related lncRNAs were found to be strongly associated with the survival of the TCGA-LUAD dataset. Utilizing multivariate Cox regression analysis, 10 lncRNAs were selected to establish the risk model. The K-M survival curves and the ROC (AUC) curves proved that the risk model has a fine predictive effect. The GO enrichment analysis indicated that the effect of the differentially expressed genes between high-risk and low-risk groups is mainly involved in immune response and intercellular interaction. The KEGG enrichment analysis indicated that the differentially expressed genes between high-risk and low-risk groups are mainly involved in endocytosis and the MAPK signaling pathway. The expression of genes related to the efficacy of immunotherapy was significantly different between the two groups. A qRT-PCR test verified the expression level of lncRNAs in LUAD cells in the risk model. The AUC of ROC of 5 years in the independent validation dataset showed that this model had superior accuracy. Western blot analysis verified the change of EMT-related marker in protein level.ConclusionThe immune lncRNA risk model established by us could better predict the prognosis of patients with LUAD.


2020 ◽  
Author(s):  
Bangyou Zuo ◽  
Haitao Zhao ◽  
Jin Bian ◽  
Junyun Long ◽  
Xu Yang ◽  
...  

Abstract Background The function of exosome includes cell-to-cell communication, neovascularization, and metastasis of cancer cell and drug resistance, which plays an important part in the occurrence and progression of hepatocellular carcinoma (HCC). Because the mechanism in this area is less studied, our goal is to identify exosome-related genes in HCC, establish a reliable prognostic model for liver cancer patients, and explore its underlying mechanisms. Methods The exoRbase database and The Cancer Genome Atlas (TCGA) database were used to analyze differentially expressed genes (DEGs). Cox regression and LASSO analysis were applied to determine DEGs closely related to overall survival (OS). Then the exosome-related prognostic model was constructed in TCGA and validated in the database of International Cancer Genome Consortium (ICGC). Nomogram graph was performed to predict the survival. CIBERSORT was used to estimate the score of different type of immune cells. DEGs related to immunotherapy are used to predict the effect of immunotherapy. Results 48 exosome-related DEGs were obtained and five genes (XPO1, IFI30, FBXO16, CALM1, MORC3) among them were selected to construct predictive model. Then we divided the HCC patients into low-risk and high-risk groups by the best cut-off value according to the X-tile software. The high-risk related to exosome were significantly associated with a poor prognosis. Moreover, the features related to exosome could positively regulate immune response. At the same time, the proportion of T cell regulatory factors (Tregs) and macrophages M2 is higher in the high-risk group, and high-risk group exhibited higher expression of immune checkpoint molecular including PD-L1, PD-L2, TIGIT, and IDO1. Conclusions Overall, our research showed that markers related to exosomes were potential biomarkers for the prognosis of HCC, providing an immunological perspective for the development of precision treatment.


Sign in / Sign up

Export Citation Format

Share Document