scholarly journals A Novel Five-Gene Signature for Prognosis Prediction in Hepatocellular Carcinoma

2021 ◽  
Vol 11 ◽  
Author(s):  
Lisa Su ◽  
Genhao Zhang ◽  
Xiangdong Kong

Hepatocellular carcinoma (HCC) has been a global health issue and attracted wide attention due to its high incidence and poor outcomes. In this study, our purpose was to explore an effective prognostic marker for HCC. Five cohort profile datasets from GEO (GSE25097, GSE36376, GSE62232, GSE76427 and GSE101685) were integrated with TCGA-LIHC and GTEx dataset to identify differentially expressed genes (DEGs) between normal and cancer tissues in HCC patients, then 5 upregulated differentially expressed genes and 32 downregulated DEGs were identified as common DEGs in total. Next, we systematically explored the relationship between the expression of 37 common DEGs in tumor tissues and overall survival (OS) rate of HCC patients in TCGA and constructed a novel prognostic model composed of five genes (AURKA, PZP, RACGAP1, ACOT12 and LCAT). Furthermore, the predicted performance of the five-gene signature was verified in ICGC and another independent clinical samples cohort, and the results demonstrated that the signature performed well in predicting the OS rate of patients with HCC. What is more, the signature was an independent hazard factor for HCC patients when considering other clinical factors in the three cohorts. Finally, we found the signature was significantly associated with HCC immune microenvironment. In conclusion, the prognostic five-gene signature identified in our present study could efficiently classify patients with HCC into subgroups with low and high risk of longer overall survival time and help clinicians make decisions for individualized treatment.

2020 ◽  
Author(s):  
Cuiyun Wu ◽  
Yaosheng Luo ◽  
Yinghui Chen ◽  
Hongling Qu ◽  
Lin Zheng ◽  
...  

Abstract Background: Accurate prediction of overall survival is important for prognosis and the assignment of appropriate personalized clinical treatment in hepatocellular carcinoma (HCC) patients. The aim of the present study was to establish an optimal gene model for the independent prediction of prognosis associated with common clinical patterns.Methods: Gene expression profiles and the corresponding clinical information of the LIHC cohort were obtained from The Cancer Genome Atlas. Differentially expressed genes were found using the R package “limma”. Subsequently, a prognostic gene signature was developed using the LASSO Cox regression model. Kaplan–Meier, log-rank, and receiver operating characteristic (ROC) analyses were performed to verify the predictive accuracy of the prognostic model. Finally, a nomogram and calibration plot were created using the “rms” package.Results: Differentially expressed genes were screened with threshold criteria (FDR < 0.01 and |log FC|>3) and 563 differentially expressed genes were obtained, including 448 downregulated and 115 upregulated genes. Using the LASSO Cox regression model, a prognostic gene signature was developed based on nine genes,IQGAP3, BIRC5, PTTG1, STC2, CDKN3, PBK, EXO1, NEIL3, and HOXD9, the expression levels of which were quantitated using RT-qPCR. According to the risk scores, patients were separated into high-risk and low-risk groups. Patients with lower risk scores generally had a better survival rate than those with higher risk scores. The mortality rate in the high-risk group was 42.02%, while that in the low-risk group was 12.50%. Results of the log-rank test showed significant differences in mortality between the two groups (HR: 4.86; 95% CI: 2.72–8.69; P = 1.01E-08). Subsequently, we assessed the prognostic accuracy of the gene signature using an ROC curve and the results show good sensitivity and specificity, with an average area under the curve (AUC) of 0.81 at 5 years (P < 0.01). Following multivariate adjustment for conventional clinical patterns, the prognostic gene signature remained a powerful and independent factor (HR: 4.70; 95% CI: 2.61–8.38; P = 2.06E-07), confirming its robust predictive ability of overall survival in HCC patients. Finally, a nomogram was established based on the gene signature and four clinicopathological features, which demonstrated an advantageous discriminating ability with the potential to facilitate clinical decision-making in HCC.Conclusion: Our prognostic gene signature can be used as a combined biomarker for the independent prediction of overall survival in HCC patients. Moreover, we created a nomogram that can be used to infer prognosis and aid individualized decisions regarding treatment and surveillance.


2020 ◽  
Vol 27 (1) ◽  
pp. 107327482097711
Author(s):  
Jiasheng Lei ◽  
Dengyong Zhang ◽  
Chao Yao ◽  
Sheng Ding ◽  
Zheng Lu

Background: Hepatocellular carcinoma (HCC) remains the third leader cancer-associated cause of death globally, but the etiological basis for this complex disease remains poorly clarified. The present study was thus conceptualized to define a prognostic immune-related gene (IRG) signature capable of predicting immunotherapy responsiveness and overall survival (OS) in patients with HCC. Methods: Five differentially expressed IRG associated with HCC were established the immune-related risk model through univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses. Patients were separated at random into training and testing cohorts, after which the association between the identified IRG signature and OS was evaluated using the “survival” R package. In addition, maftools was leveraged to assess mutational data, with tumor mutation burden (TMB) scores being calculated as follows: (total mutations/total bases) × 106. Immune-related risk term abundance was quantified via “ssGSEA” algorithm using the “gsva” R package. Results: HCC patients were successfully stratified into low-risk and high-risk groups based upon a signature composed of 5 differentially expressed IRGs, with overall survival being significantly different between these 2 groups in training cohort, testing cohort and overall patient cohort ( P = 1.745e-06, P = 1.888e-02, P = 4.281e-07). No association was observed between TMB and this IRG risk score in the overall patient cohort ( P = 0.461). Notably, 19 out of 29 immune-related risk terms differed substantially in the overall patient dataset. These risk terms mainly included checkpoints, human leukocyte antigens, natural killer cells, dendritic cells, and major histocompatibility complex class I. Conclusion: In summary, an immune-related prognostic gene signature was successfully developed and used to predict survival outcomes and immune system status in patients with HCC. This signature has the potential to help guide immunotherapeutic treatment planning for patients affected by this deadly cancer.


Author(s):  
Dafeng Xu ◽  
Yu Wang ◽  
Jincai Wu ◽  
Yuliang Zhang ◽  
Zhehao Liu ◽  
...  

Background: The prognosis of patients with hepatocellular carcinoma (HCC) is negatively affected by the lack of effective prognostic indicators. The change of tumor immune microenvironment promotes the development of HCC. This study explored new markers and predicted the prognosis of HCC patients by systematically analyzing immune characteristic genes.Methods: Immune-related genes were obtained, and the differentially expressed immune genes (DEIGs) between tumor and para-cancer samples were identified and analyzed using gene expression profiles from TCGA, HCCDB, and GEO databases. An immune prognosis model was also constructed to evaluate the predictive performance in different cohorts. The high and low groups were divided based on the risk score of the model, and different algorithms were used to evaluate the tumor immune infiltration cell (TIIC). The expression and prognosis of core genes in pan-cancer cohorts were analyzed, and gene enrichment analysis was performed using clusterProfiler. Finally, the expression of the hub genes of the model was validated by clinical samples.Results: Based on the analysis of 730 immune-related genes, we identified 64 common DEIGs. These genes were enriched in the tumor immunologic related signaling pathways. The first 15 genes were selected using RankAggreg analysis, and all the genes showed a consistent expression trend across multi-cohorts. Based on lasso cox regression analysis, a 5-gene signature risk model (ATG10, IL18RAP, PRKCD, SLC11A1, and SPP1) was constructed. The signature has strong robustness and can stabilize different cohorts (TCGA-LIHC, HCCDB18, and GSE14520). Compared with other existing models, our model has better performance. CIBERSORT was used to assess the landscape maps of 22 types of immune cells in TCGA, GSE14520, and HCCDB18 cohorts, and found a consistent trend in the distribution of TIIC. In the high-risk score group, scores of Macrophages M1, Mast cell resting, and T cells CD8 were significantly lower than those of the low-risk score group. Different immune expression characteristics, lead to the different prognosis. Western blot demonstrated that ATG10, PRKCD, and SPP1 were highly expressed in cancer tissues, while IL18RAP and SLC11A1 expression in cancer tissues was lower. In addition, IL18RAP has a highly positive correlation with B cell, macrophage, Neutrophil, Dendritic cell, CD8 cell, and CD4 cell. The SPP1, PRKCD, and SLC11A1 genes have the strongest correlation with macrophages. The expression of ATG10, IL18RAP, PRKCD, SLC11A1, and SPP1 genes varies among different immune subtypes and between different T stages.Conclusion: The 5-immu-gene signature constructed in this study could be utilized as a new prognostic marker for patients with HCC.


2021 ◽  
pp. 153537022110535
Author(s):  
Nan Li ◽  
Kai Yu ◽  
Zhong Lin ◽  
Dingyuan Zeng

Cervical cancer mortality is the second highest in gynecological cancers. This study developed a new model based on copy number variation data and mRNA data for overall survival prediction of cervical cancer. Differentially expressed genes from The Cancer Genome Atlas dataset detected by univariate Cox regression analysis were further simplified to six by least absolute shrinkage and selection operator (Lasso) and stepwise Akaike information criterion (stepAIC). The study developed a six-gene signature, which was further verified in independent dataset. Association between immune infiltration and risk score was investigated by immune score. The relation between the signature and functional pathways was examined by gene set enrichment analysis. Ninety-nine differentially expressed genes were detected, and C11orf80, FOXP3, GSN, HCCS, PGAM5, and RIBC2 were identified as key genes to construct a six-gene signature. The prognostic signature showed a significant correlation with overall survival (hazard ratio, HR = 3.45, 95% confidence interval (CI) = 2.08–5.72, p <  0.00001). Immune score showed a negative correlation with the risk score calculated by the signature ( p <  0.05). Four immune-related pathways were closely associated with risk score ( p <  0.0001). The six-gene prognostic signature was an effective tool to predict overall survival of cervical cancer. In conclusion, the newly identified six genes may be considered as new drug targets for cervical cancer treatment.


2021 ◽  
Vol 12 ◽  
Author(s):  
Jing Jia ◽  
Yuhan Dai ◽  
Qing Zhang ◽  
Peiyu Tang ◽  
Qiang Fu ◽  
...  

BackgroundGrowing evidence has revealed the crucial roles of stromal cells in the microenvironment of various malignant tumors. However, efficient prognostic signatures based on stromal characteristics in colon cancer have not been well-established yet. The present study aimed to construct a stromal score-based multigene prognostic prediction model for colon cancer.MethodsStromal scores were calculated based on the expression profiles of a colon cancer cohort from TCGA database applying the ESTIMATE algorithm. Linear models were used to identify differentially expressed genes between low-score and high-score groups by limma R package. Univariate, LASSO, and multivariate Cox regression models were used successively to select the prognostic gene signature. Two independent datasets from GEO were used as external validation cohorts.ResultsLow stromal score was demonstrated to be a favorable factor to the overall survival of colon cancer patients in TCGA cohort (p = 0.0046). Three hundred and seven stromal score-related differentially expressed genes were identified. Through univariate, LASSO, and multivariate Cox regression analyses, a gene signature consisting of LEP, NOG, and SYT3 was recognized to build a prognostic prediction model. Based on the predictive values estimated by the established integrated model, patients were divided into two groups with significantly different overall survival outcomes (p &lt; 0.0001). Time-dependent Receiver operating characteristic curve analyses suggested the satisfactory predictive efficacy for the 5-year overall survival of the model (AUC value = 0.733). A nomogram with great predictive performance combining the multigene prediction model and clinicopathological factors was developed. The established model was validated in an external cohort (AUC value = 0.728). In another independent cohort, the model was verified to be of significant prognostic value for different subgroups, which was demonstrated to be especially accurate for young patients (AUC value = 0.763).ConclusionThe well-established model based on stromal score-related gene signature might serve as a promising tool for the prognostic prediction of colon cancer.


2021 ◽  
Author(s):  
Lianxiang Luo ◽  
Xinyue Yao ◽  
Fangfang Huang ◽  
Hui Luo

Abstract Ferroptosis is a novel type of cell death depending on iron, which has been confirmed strongly related with the development of tumor. Hepatocellular carcinoma (HCC) is a malignancy with high incidence. Despite some reports demonstrated the relation between ferroptosis-related genes and HCC, more details have not been excavated. In present study, we analyzed ferroptosis-related genes with their clinical information from TCGA-LIHC project to find out prognostic genes in HCC. And four genes (GPX2, MT3, PRDX1 and SRXN1) were established as a prognostic model after differentially expressed analysis, Cox regression analyses and LASSO approach. High-risk group separating by cutoff value with poor prognosis was proved, and risk score regarded as an independent prognostic factor. Subsequently, enrichment analyses were processed to the differentially expressed genes, we found that genes generally enriched in ferroptosis-related functions and pathways, also in some immune cells and functions with different immune status. Eventually, we constructed ferroptosis potential index (FPI) to reveal the functional roles of ferroptosis in tumor tissues. The immunohistochemistry performed prognostic genes expression in normal and tumor tissues. In conclusion, these results demonstrated the four-gene signature can be a biomarker for predicting HCC prognosis and its members can be drug target genes for HCC treatment.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Kai Yu ◽  
Huan Yang ◽  
Qiao-li Lv ◽  
Li-chong Wang ◽  
Zi-long Tan ◽  
...  

Abstract Background Glioblastoma is the most common primary malignant brain tumor. Because of the limited understanding of its pathogenesis, the prognosis of glioblastoma remains poor. This study was conducted to explore potential competing endogenous RNA (ceRNA) network chains and biomarkers in glioblastoma by performing integrated bioinformatics analysis. Methods Transcriptome expression data from The Cancer Genome Atlas database and Gene Expression Omnibus were analyzed to identify differentially expressed genes between glioblastoma and normal tissues. Biological pathways potentially associated with the differentially expressed genes were explored by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis, and a protein-protein interaction network was established using the STRING database and Cytoscape. Survival analysis using Gene Expression Profiling Interactive Analysis was based on the Kaplan–Meier curve method. A ceRNA network chain was established using the intersection method to align data from four databases (miRTarBase, miRcode, TargetScan, and lncBace2.0), and expression differences and correlations were verified by quantitative reverse-transcription polymerase chain reaction analysis and by determining the Pearson correlation coefficient. Additionally, an MTS assay and the wound-healing and transwell assays were performed to evaluate the effects of complement C1s (C1S) on the viability and migration and invasion abilities of glioblastoma cells, respectively. Results We detected 2842 differentially expressed (DE) mRNAs, 2577 DE long non-coding RNAs (lncRNAs), and 309 DE microRNAs (miRNAs) that were dysregulated in glioblastoma. The final ceRNA network consisted of six specific lncRNAs, four miRNAs, and four mRNAs. Among them, four DE mRNAs and one DE lncRNA were correlated with overall survival (p < 0.05). C1S was significantly correlated with overall survival (p= 0.015). In functional assays, knockdown of C1S inhibited the proliferation and invasion of glioblastoma cell lines. Conclusions We established four ceRNA networks that may influence the occurrence and development of glioblastoma. Among them, the MIR155HG/has-miR-129-5p/C1S axis is a potential marker and therapeutic target for glioblastoma. Knockdown of C1S inhibited the proliferation, migration, and invasion of glioblastoma cells. These findings clarify the role of the ceRNA regulatory network in glioblastoma and provide a foundation for further research.


2021 ◽  
Author(s):  
Zhigang Wang ◽  
Leyu Pan ◽  
Deliang Guo ◽  
Xiaofeng Luo ◽  
Jie Tang ◽  
...  

Author(s):  
Heng Cao ◽  
Peng Guo ◽  
Xiaohui Wu ◽  
Jiankun Li ◽  
Chenlong Ge ◽  
...  

Background: Hepatocellular carcinoma (HCC) is one of the most common malignant tumors of digestive tract in the world. Therefore, it is important to carry out studies on the molecular mechanisms of early diagnosis and treatment of HCC to reduce mortality. Methods: Bioinformatic analysis was performed to explore the significant role of GCSF on the occurrence and development of neoplasm. Differently expressed genes (DEGs) were screened, and the significant hub genes related with GCSF were identified by the multiple algorithms of Cytoscape. Functional annotation for DEGs, pathological stage and overall survival analysis were implemented. In addition, the verification for the role of GCSF on HCC was made via the clinical samples. A total of 70 participates diagnosed as HCC were recruited from November 2014 to November 2019. The immunohistochemistry assay, qRT-PCR, receiver operating characteristic (ROC) curves, and overall survival analysis were carried out. Results: GCSF was related with the tumor size, and the expression of GCSF was up-regulated in hepatocellular carcinoma tissues. The enrichment results of GO and KEGG analysis were mainly enriched in “Inflammatory response”, “Protein binding”, “Metabolic pathways”, and “Proteasome”. The tumor diameter (P < 0.001), and survival time (P < 0.001) were significantly associated with expression of GCSF via the verification of clinical data. The univariate and multivariate Cox proportional regression analysis manifested that high expression of GCSF in patients with HCC was related to poor OS. Conclusion: The expression level of GCSF is significantly associated with the prognostic survival of HCC, and it is expected to become a new prognostic marker of HCC, providing a novel idea for future basic research as well as targeted therapy.


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.


Sign in / Sign up

Export Citation Format

Share Document