scholarly journals A novel prognostic index model constructed with five autophagy-related genes may be a potential prognostic biomarker and therapeutic target in liver hepatocellular carcinoma

2020 ◽  
Author(s):  
Qingqing Zhu ◽  
Jia Wang ◽  
Menghan Liu ◽  
Qiujing Zhang ◽  
Miaomiao Shen ◽  
...  

Abstract Background Early diagnosis and effective treatment of liver hepatocellular carcinoma (LIHC) are keys to improving the prognosis of patients. Increasing evidences clarify that autophagy-related genes (ARGs) make great differences to the generation and progression of LIHC, and may serve as prognostic biomarkers for LIHC. Methods We randomly divided the LIHC patients in The Cancer Genome Atlas (TCGA) into the training and testing group. Next, use the training group to perform univariate Cox, LASSO and multivariate Cox analysis to construct our prognostic index (PI) model for LIHC; use the testing group and the whole TCGA set to make internal validations; use the International Cancer Genome Consortium to make external validations; and use the whole TCGA, GSE14520 and Oncomine to exam the expression patterns of the five ARGs. Then, we performed the ROC curve as well as univariate and multivariate analysis to evaluate the independent prognostic prediction power of the PI model, and made nomograms to estimate 1,3,5-year survival rate of LIHC patients. Besides, we conducted functional enrichment analyses of differentially expressed ARGs with GO, KEGG and GSEA, and made drug sensitivity analysis for the PI model via the GDSC database. Results A novel PI model which was composed of five key ARGs ( ATG9A , EIF2S1 , GRID1 , SAR1A and SQSTM1 ) succeeded to be constructed. All the internal and external validations testified that the PI model could well distinguish high-risk patients from low-risk ones, with AUC values > 0.60. Further comparison analysis showed that the PI model was no less than some common prognostic factors. People can estimate the 1,3,5-year survival rate of individual LIHC patient with the nomograms. Additionally, we obtained 62 differentially expressed ARGs and studied the potential mechanisms or pathways. Furthermore, we also found some potential targeted drugs associated to the five ARGs for LIHC patients. Conclusions The novel five-ARGs PI model has great potential to serve as a diagnostic or prognostic biomarker and therapeutic target in LIHC, which may guide future clinical applications to some extent and improve the outcome of LIHC patients.

2021 ◽  
Vol 12 ◽  
Author(s):  
Huaifeng Liu ◽  
Yu Gao ◽  
Shangshang Hu ◽  
Zhengran Fan ◽  
Xianggang Wang ◽  
...  

Liver Hepatocellular Carcinoma (LIHC), a malignant tumor with high incidence and mortality, is one of the most common cancers in the world. Multiple studies have found that the aberrant expression of rhythm genes is closely related to the occurrence of LIHC. This study aimed to use bioinformatics analysis to identify differentially expressed rhythm genes (DERGs) in LIHC. A total of 563 DERGs were found in LIHC, including 265 downregulated genes and 298 upregulated genes. KEGG pathway enrichment and GO analyses showed that DERGs were significantly enriched in rhythmic and metabolic processes. Survival analysis revealed that high expression levels of CNK1D, CSNK1E, and NPAS2 were significantly associated with the low survival rate in LIHC patients. Through cell experiment verification, the mRNA expression levels of CSNK1D, CSNK1E, and NPAS2 were found to be strongly upregulated, which was consistent with the bioinformatics analysis of LIHC patient samples. A total of 23 nodes and 135 edges were involved in the protein–protein interaction network of CSNK1D, CSNK1E, and NPAS2 genes. Clinical correlation analyses revealed that CSNK1D, CSNK1E, and NPAS2 expression levels were high-risk factors and independently connected with the overall survival rate in LIHC patients. In conclusion, the identification of these DERGs contributes to the exploration of the molecular mechanisms of LIHC occurrence and development and may be used as diagnostic and prognostic biomarkers and molecular targets for chronotherapy in LIHC patients in the future.


2020 ◽  
Vol 11 ◽  
Author(s):  
Xuanlong Du ◽  
Yewei Zhang

BackgroundHepatocellular carcinoma (HCC) is a common malignant tumor with high mortality and poor prognoses around the world. Ferroptosis is a new form of cell death, and some studies have found that it is related to cancer immunotherapy. The aim of our research was to find immunity- and ferroptosis-related biomarkers to improve the treatment and prognosis of HCC by bioinformatics analysis.MethodsFirst, we obtained the original RNA sequencing (RNA-seq) expression data and corresponding clinical data of HCC from The Cancer Genome Atlas (TGCA) database and performed differential analysis. Second, we used immunity- and ferroptosis-related differentially expressed genes (DEGs) to perform a computational difference algorithm and Cox regression analysis. Third, we explored the potential molecular mechanisms and properties of immunity- and ferroptosis-related DEGs by computational biology and performed a new prognostic index based on immunity- and ferroptosis-related DEGs by multivariable Cox analysis. Finally, we used HCC data from International Cancer Genome Consortium (ICGC) data to perform validation.ResultsWe obtained 31 immunity (p < 0.001)- and 14 ferroptosis (p < 0.05)-related DEGs correlated with overall survival (OS) in the univariate Cox regression analysis. Then, we screened five immunity- and two ferroptosis-related DEGs (HSPA4, ISG20L2, NRAS, IL17D, NDRG1, ACSL4, and G6PD) to establish a predictive model by multivariate Cox regression analysis. Receiver operating characteristic (ROC) and Kaplan–Meier (K–M) analyses demonstrated a good performance of the seven-biomarker signature. Functional enrichment analysis including Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed that the seven-biomarker signature was mainly associated with HCC-related biological processes such as nuclear division and the cell cycle, and the immune status was different between the two risk groups.ConclusionOur results suggest that this specific seven-biomarker signature may be clinically useful in the prediction of HCC prognoses beyond conventional clinicopathological factors. Moreover, it also brings us new insights into the molecular mechanisms of HCC.


2020 ◽  
Author(s):  
Gaochen Lan ◽  
Xiaoling Yu ◽  
Yanna Zhao ◽  
Jinjian Lan ◽  
Wan Li ◽  
...  

Abstract Background: Breast cancer is the most common malignant disease among women. At present, more and more attention has been paid to long non-coding RNAs (lncRNAs) in the field of breast cancer research. We aimed to investigate the expression profiles of lncRNAs and construct a prognostic lncRNA for predicting the overall survival (OS) of breast cancer.Methods: The expression profiles of lncRNAs and clinical data with breast cancer were obtained from The Cancer Genome Atlas (TCGA). Differentially expressed lncRNAs were screened out by R package (limma). The survival probability was estimated by the Kaplan‑Meier Test. The Cox Regression Model was performed for univariate and multivariate analysis. The risk score (RS) was established on the basis of the lncRNAs’ expression level (exp) multiplied regression coefficient (β) from the multivariate cox regression analysis with the following formula: RS=exp a1 * β a1 + exp a2 * β a2 +……+ exp an * β an. Functional enrichment analysis was performed by Metascape.Results: A total of 3404 differentially expressed lncRNAs were identified. Among them, CYTOR, MIR4458HG and MAPT-AS1 were significantly associated with the survival of breast cancer. Finally, The RS could predict OS of breast cancer (RS=exp CYTOR * β CYTOR + exp MIR4458HG * β MIR4458HG + exp MAPT-AS1 * β MAPT-AS1). Moreover, it was confirmed that the three-lncRNA signature could be an independent prognostic biomarker for breast cancer (HR=3.040, P=0.000).Conclusions: This study established a three-lncRNA signature, which might be a novel prognostic biomarker for breast cancer.


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.


2021 ◽  
Author(s):  
Pejman Morovat ◽  
Saman Morovat ◽  
Arash M. Ashrafi ◽  
Shahram Teimourian

Abstract Hepatocellular carcinoma (HCC) is one of the most prevalent cancers worldwide, which has a high mortality rate and poor treatment outcomes with yet unknown molecular basis. It seems that gene expression plays a pivotal role in the pathogenesis of the disease. Circular RNAs (circRNAs) can interact with microRNAs (miRNAs) to regulate gene expression in various malignancies by acting as competitive endogenous RNAs (ceRNAs). However, the potential pathogenesis roles of the ceRNA network among circRNA/miRNA/mRNA in HCC are unclear. In this study, first, the HCC circRNA expression data were obtained from three Gene Expression Omnibus microarray datasets (GSE164803, GSE94508, GSE97332), and the differentially expressed circRNAs (DECs) were identified using R limma package. Also, the liver hepatocellular carcinoma (LIHC) miRNA and mRNA sequence data were retrieved from TCGA, and differentially expressed miRNAs (DEMIs) and mRNAs (DEGs) were determined using the R DESeq2 package. Second, CSCD website was used to uncover the binding sites of miRNAs on DECs. The DECs' potential target miRNAs were revealed by conducting an intersection between predicted miRNAs from CSCD and downregulated DEMIs. Third, some related genes were uncovered by intersecting targeted genes predicted by miRWalk and targetscan online tools with upregulated DEGs. The ceRNA network was then built using the Cytoscape software. The functional enrichment and the overall survival time of these potential targeted genes were analyzed, and a PPI network was constructed in the STRING database. Network visualization was performed by Cytoscape, and ten hub genes were detected using the CytoHubba plugin tool. Four DECs (hsa_circ_0000520, hsa_circ_0008616, hsa_circ_0070934, hsa_circ_0004315) were obtained and six miRNAs (hsa-miR-542-5p, hsa-miR-326, hsa-miR-511-5p, hsa-miR-195-5p, hsa-miR-214-3p, and hsa-miR-424-5p) which are regulated by the above DECs were identified. Then 543 overlapped genes regulated by six miRNAs mentioned above were predicted. Functional enrichment analysis showed that these genes are mostly associated with cancer regulation functions. Ten hub genes (TTK،AURKB, KIF20A، KIF23، CEP55، CDC6، DTL، NCAPG، CENPF، PLK4) have been screened from the PPI network of the 204 survival-related genes. KIF20A, NCAPG, TTK, PLK4, and CDC6 were selected for the highest significant p-values. In the end, a circRNA-miRNA-mRNA regulatory axis was established for five final selected hub genes. This study implies the potential pathogenesis of the obtained network and proposes that the two DECs (has_circ_0070934 and has_circ_0004315) may be important prognostic factor for HCC.


2021 ◽  
Vol 11 ◽  
Author(s):  
Peng Liu ◽  
Jinhong Wei ◽  
Feiyu Mao ◽  
Zechang Xin ◽  
Heng Duan ◽  
...  

Hepatocellular carcinoma (HCC) is one of the most common types of cancer worldwide and its incidence continues to increase year by year. Endoplasmic reticulum stress (ERS) caused by protein misfolding within the secretory pathway in cells and has an extensive and deep impact on cancer cell progression and survival. Growing evidence suggests that the genes related to ERS are closely associated with the occurrence and progression of HCC. This study aimed to identify an ERS-related signature for the prospective evaluation of prognosis in HCC patients. RNA sequencing data and clinical data of patients from HCC patients were obtained from The Cancer Genome Atlas (TCGA) and The International Cancer Genome Consortium (ICGC). Using data from TCGA as a training cohort (n=424) and data from ICGC as an independent external testing cohort (n=243), ERS-related genes were extracted to identify three common pathways IRE1, PEKR, and ATF6 using the GSEA database. Through univariate and multivariate Cox regression analysis, 5 gene signals in the training cohort were found to be related to ERS and closely correlated with the prognosis in patients of HCC. A novel 5-gene signature (including HDGF, EIF2S1, SRPRB, PPP2R5B and DDX11) was created and had power as a prognostic biomarker. The prognosis of patients with high-risk HCC was worse than that of patients with low-risk HCC. Multivariate Cox regression analysis confirmed that the signature was an independent prognostic biomarker for HCC. The results were further validated in an independent external testing cohort (ICGC). Also, GSEA indicated a series of significantly enriched oncological signatures and different metabolic processes that may enable a better understanding of the potential molecular mechanism mediating the progression of HCC. The 5-gene biomarker has a high potential for clinical applications in the risk stratification and overall survival prediction of HCC patients. In addition, the abnormal expression of these genes may be affected by copy number variation, methylation variation, and post-transcriptional regulation. Together, this study indicated that the genes may have potential as prognostic biomarkers in HCC and may provide new evidence supporting targeted therapies in HCC.


2021 ◽  
Vol 12 ◽  
Author(s):  
Xingte Chen ◽  
Lei Wang ◽  
Liang Hong ◽  
Zhixiong Su ◽  
Xiaohong Zhong ◽  
...  

Background: Aging is a well-studied concept, but no studies have comprehensively analyzed the association between aging-related genes (AGs) and hepatocellular carcinoma (HCC) prognosis.Methods: Gene candidates were selected from differentially expressed genes and prognostic genes in The Cancer Genome Atlas (TCGA) database. A gene risk score for overall survival prediction was established using the least absolute shrinkage and selection operator (LASSO) regression analysis, and this was validated using data from the International Cancer Genome Consortium (ICGC) database. Functional analysis was conducted using gene ontology enrichment, Kyoto Encyclopedia of Genes and Genomes analysis, gene set enrichment analysis, and immune microenvironment and tumor stemness analyses.Results: Initially, 72 AGs from the TCGA database were screened as differentially expressed between normal and tumor tissues and as genes associated with HCC prognosis. Then, seven AGs (POLA1, CDK1, SOCS2, HDAC1, MAPT, RAE1, and EEF1E1) were identified using the LASSO regression analysis. The seven AGs were used to develop a risk score in the training set, and the risk was validated to have a significant prognostic value in the ICGC set (p < 0.05). Patients with high risk scores had lower tumor differentiation, higher stage, and worse prognosis (all p < 0.05). Multivariate Cox regression analyses also confirmed that the risk score was an independent prognostic factor for HCC in both the TCGA and ICGC sets (all p < 0.05). Further analysis showed that a high risk score was correlated with the downregulation of metabolism and tumor immunity.Conclusion: The risk score predicts HCC prognosis and could thus be used as a biomarker not only for predicting HCC prognosis but also for deciding on treatment.


2020 ◽  
Vol 40 (7) ◽  
Author(s):  
Weidong Shi ◽  
Lanyun Feng ◽  
Shu Dong ◽  
Zhouyu Ning ◽  
Yongqiang Hua ◽  
...  

Abstract The present study aimed to screen the immune-related genes (IRGs) in patients with liver hepatocellular carcinoma (LIHC) and construct a synthetic index for indicating the prognostic outcomes. The bioinformatic analysis was performed on the data of 374 cancer tissues and 50 normal tissues, which were downloaded from TCGA database. We observed that 17 differentially expressed IRGs were significantly associated with survival in LIHC patients. These LIHC-specific IRGs were validated with function analysis and molecular characteristics. Cox analysis was applied for constructing a RiskScore for predicting the survival. The RiskScore involved six IRGs and corresponding coefficients, which was calculated with the following formula: RiskScore = [Expression level of FABP5 *(0.064)] + [Expression level of TRAF3 * (0.198)] + [Expression level of CSPG5 * (0.416)] + [Expression level of IL17D * (0.197)] + [Expression level of STC2 * (0.036)] + [Expression level of BRD8 * (0.140)]. The RiskScore was positively associated with the poor survival, which was verified with the dataset from ICGC database. Further analysis revealed that the RiskScore was independent of any other clinical feature, while it was linked with the infiltration levels of six types of immune cells. Our study reported the survival-associated IRGs in LIHC and then constructed IRGs-based RiskScore as prognostic indicator for screening patients with high risk of short survival. Both the screened IRGs and IRGs-based RiskScore were clinically significant, which may be informative for promoting the individualized immunotherapy against LIHC.


2021 ◽  
Vol 12 ◽  
Author(s):  
Rong Deng ◽  
Xiaohan Cui ◽  
Yuxiang Dong ◽  
Yanqiu Tang ◽  
Xuewen Tao ◽  
...  

BackgroundCircular RNAs (circRNAs) are now under hot discussion as novel promising biomarkers for patients with hepatocellular carcinoma (HCC). The purpose of our study is to identify several competing endogenous RNA (ceRNA) networks related to the prognosis and progression of HCC and to further investigate the mechanism of their influence on tumor progression.MethodsFirst, we obtained gene expression data related to liver cancer from The Cancer Genome Atlas (TCGA) database (http://www.portal.gdc.cancer.gov/), including microRNA (miRNA) sequence, RNA sequence, and clinical information. A co-expression network was constructed through the Weighted Correlation Network Analysis (WGCNA) software package in R software. The differentially expressed messenger RNAs (DEmRNAs) in the key module were analyzed with the Database for Annotation Visualization and Integrated Discovery (DAVID) (https://david.ncifcrf.gov/summary.jsp) to perform functional enrichment analysis including Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO). The data of miRNA expression and clinical information downloaded from TCGA were utilized for survival analysis to detach the prognostic value of the DEmiRNAs of the key module.ResultsThe 201 differentially expressed miRNAs (DEmiRNAs) and 3,783 DEmRNAs were preliminarily identified through differential expression analysis. The co-expression networks of DEmiRNAs and DEmRNAs were constructed with WGCNA. Further analysis confirmed four miRNAs in the most significant module (blue module) were associated with the overall survival (OS) of patients with liver cancer, including hsa-miR-92b-3p, hsa-miR-122-3p, hsa-miR-139-5p, and hsa-miR-7850-5p. DAVID was used for functional enrichment analysis of 286 co-expressed mRNAs. The GO analysis results showed that the top enriched GO terms were oxidation–reduction process, extracellular exosome, and iron ion binding. In KEGG pathway analysis, the top three enriched terms included metabolic pathways, fatty acid degradation, and valine, leucine, and isoleucine degradation. In addition, we intersected the miRNA–mRNA interaction prediction results with the differentially expressed and prognostic mRNAs. We found that hsa-miR-92b-3p can be related to CPEB3 and ACADL. By overlapping the data of predicted circRNAs by circBank and differentially expressed circRNAs of GSE94508, we screened has_circ_0077210 as the upstream regulatory molecule of hsa-miR-92b-3p. Hsa_circ_0077210/hsa-miR-92b-3p/cytoplasmic polyadenylation element binding protein-3 (CPEB3) and acyl-Coenzyme A dehydrogenase, long chain (ACADL) were validated in HCC tissue.ConclusionOur research provides a mechanistic elucidation of the unknown ceRNA regulatory network in HCC. Hsa_circ_0077210 might serve a momentous therapeutic role to restrain the occurrence and development of HCC.


Sign in / Sign up

Export Citation Format

Share Document