A 4-gene prognostic signature predicting survival in hepatocellular carcinoma

A prediction model for prognosis of gastric adenocarcinoma based on six metabolism-related genes

Background

The study of tumor metabolism is of great value to elucidate the mechanism of tumorigenesis and predict the prognosis of patients. However, the prognostic role of metabolism-related genes (MRGs) in gastric adenocarcinoma (GAD) remains poorly understood.

Methods

We downloaded the gene chip dataset GSE79973 (n = 20) of GAD from the Gene Expression Omnibus (GEO) database to compare differentially expressed genes (DEGs) between normal and tumor tissues. We then extracted MRGs from these DEGs and systematically investigated the prognostic value of these differential MRGs for predicting patients' overall survival by univariable and multivariable Cox regression analysis. Six metabolic genes (ACOX3, APOE, DIO2, HSD17B4, NUAK1, and WHSC1L1) were identified as prognosis-associated hub genes, which were used to build a prognostic model in the training dataset GSE15459 (n = 200), and then validated in the dataset GSE62254 (n = 300).

Results

Patients were divided into high-risk and low-risk subgroups based on the model's risk score, and it was found that patients in the high-risk subgroup had shorter overall survival than those in the low-risk subgroup, both in the training and testing datasets. In addition, for the training and testing cohorts, the area under the ROC curve of the prognostic model for one-year survival prediction was 0.723 and 0.667, respectively, indicating that the model has good predictive performance. Furthermore, we established a nomogram based on tumor stage and risk score to effectively predict the overall survival (OS) of GAD patients. The expression of 6 MRGs at the protein level was confirmed by immunohistochemistry (IHC). Kaplan-Meier survival analysis further confirmed that their expression influenced OS in GAD patients.

Conclusion

Collectively, the 6 MRGs signature might be a reliable tool for assessing OS in GAD patients, with potential application value in clinical decision-making and individualized therapy.

The value of exosome-derived noncoding RNAs in colorectal cancer proliferation, metastasis, and clinical applications

Colorectal cancer (CRC) is the third most common cancer in the world today, and its incidence and mortality rates are increasing every year. The ease of proliferation and metastasis of CRC has long been an important reason for its high mortality rate. Exosomes serve as key mediators that mediate communication between tumor cells and various other cells. Non-coding RNAs (ncRNAs) have been shown to play a key role in apoptosis, immunosuppression and proliferation metastasis in cancer. ncRNAs are loaded on exosomes and initiate the onset of metastasis by promoting epithelial-mesenchymal transition (EMT) at the primary site of the tumor. Meanwhile, exosome-derived ncRNAs construct a pre-metastatic niche (PMN) for CRC metastasis by forming an inflammatory microenvironment in distant organs, immunosuppression, and promoting angiogenesis and remodeling of the extracellular matrix. Here, we summarize the specific mechanisms associated with exosome-derived ncRNAs promoting local invasion and metastasis in CRC. Finally, we focus on their value for clinical application in future CRC diagnosis and treatment.

Systematic analyses to explore immune gene sets-based signature in hepatocellular carcinoma, in which IGF2BP3 contributes to tumor progression

Tumor immune microenvironment (TIME) is of critical importance for the development and therapeutic response of hepatocellular carcinoma (HCC). However, limited studies have investigated immune-related indicators for clinical supervision and decision. The current study aimed to develop an improved prognostic signature based on TIME. HCC patients from TCGA and ICGC database were classified into three subtypes (Immunity High, Immunity Medium and Immunity Low) according to ssGSEA scores of 29 immune gene sets. Differentially expressed immune-related genes (DE IRGs) between Immune High and Low groups were screened with an adjusted P < 0.05. Weighted gene co-expression network analysis (WGCNA) was used to establish gene co-expression modules of differentially expressed genes (DEGs) between tumor and normal tissues. 45 survival-related immune genes (SRIGs) were identified at points of intersection between hub genes and DE IRGs. By performing Cox regression and LASSO analysis, 3 of the 45 SRIGs were screened to establish a prognostic model. Patients with high risk scores exhibited worse survival outcome and poorer response to chemotherapy. Potential mechanisms of chemotherapy resistance also have been discussed. More significantly, high -risk patients showed increased immune cell infiltration and checkpoints, which suggested a benefit of immunotherapy. In addition, knockdown of IGF2BP3 was determined to significantly inhibit cell proliferation and migration in HCC. Our immune-related model may be an effective tool for precise diagnosis and treatment of HCC. It may help to select patients suitable for chemotherapy, and immunotherapy.

Integrative Analysis Identifies Cell-Type-Specific Genes Within Tumor Microenvironment as Prognostic Indicators in Hepatocellular Carcinoma

Background

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related mortality worldwide, but effective early detection and prognostication methods are lacking.

Methods

The Cox regression model was built to stratify the HCC patients. The single-cell RNA sequencing data analysis and gene set enrichment analysis were employed to investigate the biological function of identified markers. PLCB1 gain- or loss-of-function experiments were performed, and obtained HCC samples were analyzed using quantitative real-time PCR and immunohistochemistry assay to validate the biological function of identified markers.

Results

In this study, we developed a model using optimized markers for HCC recurrence prediction. Specifically, we screened out 8 genes through a series of data analyses, and built a multivariable Cox model based on their expression. The risk stratifications using the Eight-Gene Cox (EGC) model were closely associated with the recurrence-free survivals (RFS) in both training and three validation cohorts. We further demonstrated that this risk stratification could serve as an independent predictor in predicting HCC recurrence, and that the EGC model could outperform other models. Moreover, we also investigated the cell-type-specific expression patterns of the eight recurrence-related genes in tumor microenvironment using single-cell RNA sequencing data, and interpreted their functional roles from correlation and gene set enrichment analyses, in vitro and in vivo experiments. Particularly, PLCB1 and SLC22A7 were predominantly expressed in malignant cells, and they were predicted to promote angiogenesis and to help maintain normal metabolism in liver, respectively. In contrast, both FASLG and IL2RB were specifically expressed in T cells, and were highly correlated with T cell marker genes, suggesting that these two genes might assist in maintaining normal function of T cell-mediated immune response in tumor tissues.

Conclusion

In conclusion, the EGC model and eight identified marker genes could not only facilitate the accurate prediction of HCC recurrence, but also improve our understanding of the mechanisms behind HCC recurrence.

Integrated analysis of single-cell RNA-seq dataset and bulk RNA-seq dataset constructs a prognostic model for predicting survival in human glioblastoma

BACKGROUND

Glioblastoma (GBM) is the most common primary malignant brain tumor in adults. For patients with GBM, the median overall survival (OS) is 14.6 months and the 5-year survival rate is 7.2%. It is imperative to develop a reliable model to predict the survival probability in new GBM patients. To date, most prognostic models for predicting survival in GBM were constructed based on bulk RNA-seq dataset, which failed to accurately reflect the difference between tumor cores and peripheral regions, and thus show low predictive capability. An effective prognostic model is desperately needed in clinical practice.

METHODS

We studied single-cell RNA-seq dataset and The Cancer Genome Atlas-glioblastoma multiforme (TCGA-GBM) dataset to identify differentially expressed genes (DEGs) that impact the OS of GBM patients. We then applied the least absolute shrinkage and selection operator (LASSO) Cox penalized regression analysis to determine the optimal genes to be included in our risk score prognostic model. Then, we used another dataset to test the accuracy of our risk score prognostic model.

RESULTS

We identified 2128 DEGs from the single-cell RNA-seq dataset and 6461 DEGs from the bulk RNA-seq dataset. In addition, 896 DEGs associated with the OS of GBM patients were obtained. Five of these genes (LITAF, MTHFD2, NRXN3, OSMR, and RUFY2) were selected to generate a risk score prognostic model. Using training and validation datasets, we found that patients in the low-risk group showed better OS than those in the high-risk group. We validated our risk score model with the training and validating datasets and demonstrated that it can effectively predict the OS of GBM patients.

CONCLUSION

We constructed a novel prognostic model to predict survival in GBM patients by integrating a scRNA-seq dataset and a bulk RNA-seq dataset. Our findings may advance the development of new therapeutic targets and improve clinical outcomes for GBM patients.

Non-coding RNAs as emerging regulators and biomarkers in colorectal cancer

CRC is the third most common cancer occurring worldwide and the second leading cause of cancer deaths. In the year 2020, 1,931,590 new cases of CRC and 935,173 deaths were reported. The last two decades have witnessed an intensive study of noncoding RNAs and their implications in various pathological conditions including cancer. Noncoding RNAs such as miRNAs, tsRNAs, piRNAs, lncRNAs, pseudogenes, and circRNAs have emerged as promising prognostic and diagnostic biomarkers in preclinical studies of cancer. Some of these noncoding RNAs have also been shown as promising therapeutic targets for cancer treatment. In this review, we have discussed the emerging roles of various types of noncoding RNAs in CRC and their future implications in colorectal cancer management and research.

A robust gene expression signature for NASH in liver expression data

Non-Alcoholic Fatty Liver Disease (NAFLD) is a progressive liver disease that affects up to 30% of worldwide population, of which up to 25% progress to Non-Alcoholic SteatoHepatitis (NASH), a severe form of the disease that involves inflammation and predisposes the patient to liver cirrhosis. Despite its epidemic proportions, there is no reliable diagnostics that generalizes to global patient population for distinguishing NASH from NAFLD. We performed a comprehensive multicohort analysis of publicly available transcriptome data of liver biopsies from Healthy Controls (HC), NAFLD and NASH patients. Altogether we analyzed 812 samples from 12 different datasets across 7 countries, encompassing real world patient heterogeneity. We used 7 datasets for discovery and 5 datasets were held-out for independent validation. Altogether we identified 130 genes significantly differentially expressed in NASH versus a mixed group of NAFLD and HC. We show that our signature is not driven by one particular group (NAFLD or HC) and reflects true biological signal. Using a forward search we were able to downselect to a parsimonious set of 19 mRNA signature with mean AUROC of 0.98 in discovery and 0.79 in independent validation. Methods for consistent diagnosis of NASH relative to NAFLD are urgently needed. We showed that gene expression data combined with advanced statistical methodology holds the potential to serve basis for development of such diagnostic tests for the unmet clinical need.

Immunogenomic Landscape Analysis of Prognostic Immune-Related Genes in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-related death. HBV infection is an important risk factor for the tumorigenesis of HCC, given that the inflammatory environment is closely related to morbidity and prognosis. Consequently, it is of urgent importance to explore the immunogenomic landscape to supplement the prognosis of HCC. The expression profiles of immune-related genes (IRGs) were integrated with 377 HCC patients to generate differentially expressed IRGs based on the Cancer Genome Atlas (TCGA) dataset. These IRGs were evaluated and assessed in terms of their diagnostic and prognostic values. A total of 32 differentially expressed immune-related genes resulted as significantly correlated with the overall survival of HCC patients. The Gene Ontology functional enrichment analysis revealed that these genes were actively involved in cytokine-cytokine receptor interaction. A prognostic signature based on IRGs (HSPA4, PSME3, PSMD14, FABP6, ISG20L2, TRAF3, NDRG1, NRAS, CSPG5, and IL17D) stratified patients into high-risk versus low-risk groups in terms of overall survival and remained as an independent prognostic factor in multivariate analyses after adjusting for clinical and pathologic factors. Several IRGs (HSPA4, PSME3, PSMD14, FABP6, ISG20L2, TRAF3, NDRG1, NRAS, CSPG5, and IL17D) of clinical significance were screened in the present study, revealing that the proposed clinical-immune signature is a promising risk score for predicting the prognosis of HCC.

Risk characteristics with seven epithelial-mesenchymal transition-related genes are used to predict the prognosis of patients with hepatocellular carcinoma

Background

Epithelial-mesenchymal transition (EMT)-related genes (ERGs) have been shown to play an important role in cancer invasion, tumor resistance, and tumor metastasis of hepatocellular carcinoma. This study sought to examine the prognostic value of ERGs and other pre-hepatoma genes.

Methods

Relevant data from The Cancer Genome Atlas (TCGA) were analyzed and synthesized. Specifically, 1,014 ERGs were downloaded and subject to a gene set enrichment analysis; 318 different EAG expressions were found, and the possible molecular mechanism of EAG was predicted by GO analysis and KEGG analysis. To determine the prediction of ERGS, a Cox regression model was used to establish a risk hypothesis. Based on risk patterns, patients were divided into high- or low-risk groups. Kaplan-Meier and receiver operating characteristic (ROC) curves confirmed the predictive value of the model.

Results

Seven prognostically relevant ERGs (i.e., ECT2, EZH2, MYCN, ROR2, SPP1, SQSTM1, and STC2) were identified. Using Cox's regression analysis method, appropriate cases were selected to establish a new risk prediction model. Under the risk model, the overall survival rate of the low-risk group samples was higher than that of the high-risk group samples (P<0.00001).

Conclusions

In short, we developed a risk model for liver cancer based on ERGs terminology. This model improve the postpartum treatment of patients with liver cancer.

Circulating non-coding RNAs as new biomarkers and novel therapeutic targets in colorectal cancer

Colorectal cancer (CRC) is one of the most common malignant tumors, and a large number of patients are diagnosed and die every year. Due to the lack of appropriate diagnosis, prediction and treatment, early diagnosis rate of CRC is low and the prognosis is poor. Studies have found that abnormally expressed non-coding RNAs (ncRNAs) (including microRNAs (miRNAs), circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs),etc.) play an important regulatory role in the occurrence and development of CRC. Some studies have shown that they are stable in the blood and can be detected repeatedly. They are expected to be non-invasive biomarkers for early diagnosis, prognosis evaluation, and prediction of drug sensitivity of CRC, as well as potential applications in the treatment of CRC.

Co-amplification of genes in chromosome 8q24: a robust prognostic marker in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is a leading cause of tumor-associated death worldwide, owing to its high 5-year postoperative recurrence rate and inter-individual heterogeneity. Thus, a prognostic model is urgently needed for patients with HCC. Several researches have reported that copy number amplification of the 8q24 chromosomal region is associated with low survival in many cancers. In the present work, we set out to construct a multi-gene model for prognostic prediction in HCC.

Methods

RNA sequencing and copy number variant data of tumor tissue samples of HCC from The Cancer Genome Atlas (n=328) were used to identify differentially expressed messenger RNAs of genes located on the chromosomal 8q24 region by the Wilcox test. Univariate Cox and Lasso-Cox regression analyses were carried out for the screening and construction of a prognostic multi-gene signature in The Cancer Genome Atlas cohort (n=119). The multi-gene signature was validated in a cohort from the International Cancer Genome Consortium (n=240). A nomogram for prognostic prediction was built, and the underpinning molecular mechanisms were studied by Gene Set Enrichment Analysis.

Results

We successfully established a 7-gene prognostic signature model to predict the prognosis of patients with HCC. Using the model, we divided individuals into high-risk and low-risk sets, which showed a significant difference in overall survival in the training dataset (HR =0.17, 95% CI: 0.1-0.28; P<0.001) and in the testing dataset (HR = 0.42, 95% CI: 0.23-0.74; P=0.002). Multivariate Cox regression analysis showed the signature to be an independent prognostic factor of HCC survival. A nomogram including the prognostic signature was constructed and showed a better predictive performance in short-term (1 and 3 years) than in long-term (5 years) survival. Furthermore, Gene Set Enrichment Analysis identified several pathways of significance, which may aid in explaining the underlying molecular mechanism.

Conclusions

Our 7-gene signature is a reliable prognostic marker for HCC, which may provide meaningful information for therapeutic customization and treatment-related decision making.

Development of a prognostic model for hepatocellular carcinoma using genes involved in aerobic respiration

Hepatocellular carcinoma (HCC) is the second leading cause of cancer-related death worldwide. Currently, recent risk stratification has only focused on liver function and tumor characteristics. Thus, the purpose of this study was to develop a prognostic model based on genes involved in aerobic respiration. Matched tumor and normal tissues from TCGA and ICGC cohorts were analyzed to identify 15 overlapping differential expressed genes. Cox univariate analysis of the 15 genes in the TCGA cohort revealed they were all associated with disease-specific survival (DSS) in HCC patients. Using LASSO estimation and the optimal value for penalization coefficient lambda 12 genes were selected for the prognostic model, and then HCC patients in the TCGA cohort were dichotomized into low-risk and high-risk groups. Univariate and multivariate Cox analysis demonstrated patients in low-risk group had better survival. Validation of the risk score model with the ICGC cohort produces results consistent with those of the TCGA cohort. In conclusion, this study developed and validated a prognostic model of HCC through a comprehensive analysis of genes involved in aerobic respiration. This model may help develop personalized treatments for patients with HCC.

Identification of Key Genes Related to the Prognosis of Esophageal Squamous Cell Carcinoma Based on Chip Re-Annotation

Esophageal cancer (EC) is one of the deadliest cancers worldwide. However, reliable biomarkers for early diagnosis, or those for the prognosis of therapy, remain unfulfilled goals for its subtype esophageal squamous cell carcinoma (ESCC). The purpose of this study was to identify reliable biomarkers for the diagnosis and prognosis of ESCC by gene chip re-annotation technique and downstream bioinformatics analysis. In our research, the GSE53624 dataset was downloaded from the GEO database. Then, we reannotated the gene expression probe and obtained the gene expression matrix. Differential expressed genes (DEGs) were found by R packages and they were subjected to Gene Ontology enrichment analysis and protein–protein interaction (PPI) network construction. As a result, a total of 28,885 mRNA probes were reannotated, among which 210 down-regulated and 80 up-regulated DEGs were screened out. By combining these genes set in clinical prognosis information and Western blot analysis, we found four genes with diagnostic and prognostic significance, including MMP13, SPP1, MMP10, and COL1A1. Furthermore, markers of infiltrating immune cells exhibited different DEG-related immune infiltration patterns.

Identification of a Novel Epithelial-to-Mesenchymal-related Gene Signature in Predicting Survival of Patients with Hepatocellular Carcinoma

BACKGROUND

Epithelial-mesenchymal transformation (EMT) promotes cancer metastasis including hepatocellular carcinoma. Therefore, EMT-related gene signature was explored.

OBJECTIVE

The present study was designed to develop an EMT-related gene signature for predicting the prognosis of patients with hepatocellular carcinoma.

METHODS

We conducted an integrated gene expression analysis based on tumor data of the patients with hepatocellular carcinoma from The Cancer Genome Atlas (TCGA), HCCDB18 and GSE14520 dataset. An EMT-related gene signature was constructed by least absolute shrinkage and selection operator (LASSO) and COX regression analysis of univariate and multivariate survival.

RESULTS

A 3-EMT gene signature was developed and validated based on gene expression profiles of hepatocellular carcinoma from three microarray platforms. Patients with a high risk score had a significantly worse overall survival (OS) than those with low risk scores. The EMT-related gene signature showed a high performance in accurately predicting prognosis and in examining the clinical characteristics and immune score analysis. Univariate and multivariate Cox regression analyses confirmed that the EMT-related gene signature was an independent prognostic factor for predicting survival in hepatocellular carcinoma patients. Compared with the existing models, our EMT-related gene signature reached higher area under curve (AUC).

CONCLUSION

Our findings provide novel insight into understanding EMT and help identify hepatocellular carcinoma patients with poor prognosis.

A signature of 13 autophagy‑related gene pairs predicts prognosis in hepatocellular carcinoma

Growing evidences suggest that autophagy plays a momentous part in the tumorigenesis and development of hepatocellular carcinoma (HCC). However, there are not many researches to predict the prognosis of HCC using autophagy-related genes. Therefore, based on the clinical information and RNA-Seq data of The Cancer Genome Atlas data portal (TCGA), 13 autophagy‑related gene pairs were screened to build the autophagy‑related signature to predict the prognosis by least absolute shrinkage and selection operator (LASSO) regression analysis. Besides, the International Cancer Genome Consortium (ICGC) cohort was further applied to verify the autophagy‑related prognostic signature. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene set enrichment analysis (GSEA) were also used to predict the relevant function of the autophagy-related gene pairs signature. As shown in the results, the autophagy-related gene pairs were mainly involved in process utilizing autophagic mechanism, autophagy, macroautophagy and cellular response to oxidative stress. The immune cell levels in the high-risk and low-risk group were explored, which showed that the three immune cells were obviously increased in the high-risk group, while the five immune cells were obviously increased in the low-risk group. In conclusion, an autophagy‑related prognostic signature was established to predict the prognosis of HCC patients with great accuracy and we found that autophagy‑related prognostic signature was related to infiltrating immune cells.

Identifying Prognostic Significance of RCL1 and Four-Gene Signature as Novel Potential Biomarkers in HCC Patients

BACKGROUND

Hepatocellular carcinoma (HCC) is a highly malignant disease, and it is characterized by rapid progression and low five-year survival rate. At present, there are no effective methods for monitoring the treatment and prognosis of HCC.

METHODS

The transcriptome and gene expression profiles of HCC were obtained from the Cancer Genome Atlas (TCGA) program, International Cancer Genome Consortium (ICGC), and Gene Expression Omnibus (GEO) databases. The random forest method was applied to construct a four-gene prognostic model based on RNA terminal phosphate cyclase like 1 (RCL1) expression. The Kaplan-Meier method was performed to evaluate the prognostic value of RCL1, long noncoding RNAs (AC079061, AL354872, and LINC01093), and four-gene signature (SPP1, MYBL2, TRNP1, and FTCD). We examined the relationship between RCL1 expression and immune cells infiltration, tumor mutation burden (TMB), and microsatellite instability (MSI).

RESULTS

The results of multiple databases indicated that the aberrant expression of RCL1 was associated with clinical outcome, immune cells infiltration, TMB, and MSI in HCC patients. Meanwhile, we found that long noncoding RNAs (AC079061, AL354872, and LINC01093) and RCL1 were significantly coexpressed in HCC patients. We also confirmed that the four-gene signature was an independent prognostic factor for HCC patients. Ferroptosis potential index, immune checkpoint molecules, and clinical feature were found to have obvious correlations with risk score. The area under the receiver operating characteristic curve values for the model were 0.7-0.8 in the training set and the validation set, suggesting high robustness of the four-gene signature. We then built a nomogram for facilitating the use in clinical practice.

CONCLUSION

Our study demonstrated that RCL1 and a novel four-gene signature can be used as prognostic biomarkers for predicting clinical outcome in HCC patients; and this model may assist in individualized treatment monitoring of HCC patients in clinical practice.

Development of a macrophages-related 4-gene signature and nomogram for the overall survival prediction of hepatocellular carcinoma based on WGCNA and LASSO algorithm

BACKGROUND

Immune system instability and poor prognosis are the two major clinical performance of hepatocellular carcinoma (HCC). Abnormal expression of MiR-424-5p has been reported to accelerate the progression of liver cancer, but it mediated immune cell infiltration imbalance is still unknown. We aim to mine the immune-related genes (IRGs) targeted by miR-424-5p and construct a multi-gene signature to improve the prognostic prediction of HCC.

METHODS

The HCC-related data of the cancer genome atlas (TCGA) database and the GSE14520 dataset of the gene expression omnibus (GEO) database were downloaded as the discovery dataset and the validation dataset, respectively. Weighted gene co-expression network analysis (WGCNA), the deconvolution algorithm of CIBERSORT and LASSO algorithm participated in the identification of IRGs and the development of prognostic signature and nomogram.

RESULTS

Our study found that the abundance of macrophages M0, M1 and M2 are all drastically changed during the cancerous process. A total of 920 macrophages infiltration-related LRGs were identified and a novel 4-gene signature (CDCA8, CBX2, UCK2 and SOCS2) with superior prognostic independence was established. The prognostic signature based risk score has superior capability to identify high-risk patients and predict overall survival (p < 0.001; AUC = 0.798 for 1 year; AUC = 0.748 for 3 years; AUC = 0.721 for 5 years). And it (C-index = 0.726) has a better prognostic potential than the TNM stage (C-index = 0.619), which is widely adopted in clinical practice. Additionally, the nomogram formed by combining the risk score and TNM stage further improved the accuracy of survival prediction (C-index = 0.733).

CONCLUSION

In summary, the immune landscape with abnormal infiltration of macrophages may be one of the prelude to the cancerous process. The novel macrophages-related 4-gene signature is expected to become a potential prognostic marker in liver cancer.

KPNA2-Associated Immune Analyses Highlight the Dysregulation and Prognostic Effects of GRB2, NRAS, and Their RNA-Binding Proteins in Hepatocellular Carcinoma

Karyopherin α2 (KPNA2) was reported to be overexpressed and have unfavorable prognostic effects in many malignancies including hepatocellular carcinoma (HCC). Although its contributions to inflammatory response were reported in many studies, its specific associations with immune infiltrations and immune pathways during cancer progression were unclear. Here, we aimed to identify new markers for HCC diagnosis and prognosis through KPNA2-associated immune analyses. RNA-seq expression data of HCC datasets were downloaded from The Cancer Genome Atlas and International Cancer Genome Consortium. The gene expressions were counts per million normalized. The infiltrations of 24 kinds of immune cells in the samples were evaluated with ImmuCellAI (Immune Cell Abundance Identifier). The Spearman correlations of the immune infiltrations with KPNA2 expression were investigated, and the specific positive correlation of B-cell infiltration with KPNA2 expression in HCC tumors was identified. Fifteen genes in KEGG (Kyoto Encyclopedia of Genes and Genomes) B-cell receptor signaling pathway presented significant correlations with KPNA2 expression in HCC. Among them, GRB2 and NRAS were indicated to be independent unfavorable prognostic factors for HCC overall survival. Clinical Proteomic Tumor Analysis Consortium HCC dataset was investigated to validate the results at protein level. The upregulation and unfavorable prognostic effects of KPNA2 and GRB2 were confirmed, whereas, unlike its mRNA form, NRAS protein was presented to be downregulated and have favorable prognostic effects. Through receiver operating characteristic curve analysis, the diagnostic potential of the three proteins was shown. The RNA-binding proteins (RBPs) of KPNA2, NRAS, and GRB2, downloaded via The Encyclopedia of RNA Interactomes, were investigated for their clinical significance in HCC at protein level. An eight-RBP signature with independent prognostic value and dysregulations in HCC was identified. All the RBPs were significantly correlated with MKI67 expression and at least one of KPNA2, GRB2, and NRAS at protein level in HCC, indicating their roles in HCC progression and the regulation of the three proteins. We concluded that KPNA2, GRB2, NRAS, and their RBPs might have coordinating roles in HCC immunoregulation and progression. They might be new markers for HCC diagnosis and prognosis predication and new targets for HCC immunotherapy.

A novel prognostic models for identifying the risk of hepatocellular carcinoma based on epithelial-mesenchymal transition-associated genes

Several epithelial-mesenchymal transition (EMT)-associated genes (EAGs) have been confirmed to correlate with the prognosis of hepatocellular carcinoma (HCC) patients. Herein, we explored the value of EAGs in the prognosis of HCC relying on data from The Cancer Genome Atlas (TCGA) database. A total of 200 EMT-associated genes were downloaded from the Gene set enrichment analysis (GSEA) website. Moreover, 96 differentially expressed EAGs were identified. Using Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, we forecasted the potential molecular mechanisms of EAGs. To identify prognostic EAGs, Cox regression was used in developing a prognostic risk model. Then, the Kaplan-Meier and receiver operating characteristic (ROC) curves were plotted to validate the prognostic significance of the model. A total of 5 prognostic correlated EAGs (P3H1, SPP1, MMP1, LGALS1, and ITGB5) were screened via Cox regression, which provided the basis for developing a novel prognostic risk model. Based on the risk model, patients were subdivided into high-risk and low-risk groups. The overall survival of the low-risk group was better compared to the high-risk group (P < 0.00001). The ROC curve of the risk model showed a higher AUC (Area under Curve) (AUC = 0.723) compared to other clinical features (AUC ≤ 0.511). A nomogram based on this model was constructed to predict the 1-year, 2-year, and 3-year overall survival rates (OS) of patients. Conclusively, we developed a novel HCC prognostic risk model based on the expression of EAGs, which help advance the prognostic management of HCC patients. Abbreviations: HCC: hepatocellular carcinoma; TCGA: The Cancer Genome Atlas; EMT: epithelial-mesenchymal transition; EAGs: EMT-associated genes; GSEA: gene set enrichment analysis; GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; PPI: protein-protein interaction; TF: transcription factor; ROC: receiver operating characteristic; K-M: Kaplan-Meier; AUC: the area under the ROC curve; FDR: false discovery rate; TNM: Tumor size/lymph nodes/distance metastasis.

Development and validation of a RNA binding protein gene pair-associated prognostic signature for prediction of overall survival in hepatocellular carcinoma

Background

Increasing evidence has demonstrated the correlation between hepatocellular carcinoma (HCC) prognosis and RNA binding proteins (RBPs) dysregulation. Thus, we aimed to develop and validate a reliable prognostic signature that can estimate the prognosis for HCC.

Methods

Gene expression profiling and clinical information of 374 HCC patients were derived from the TCGA data portal. The survival-related RBP pairs were determined using univariate cox-regression analysis and the signature was built based on LASSO analysis. All patients were divided patients into high-and low-risk groups according to the optimal cut off of the signature score determined by time-dependent receiver operating characteristic (ROC) curve analysis. The predictive value of the signature was further validated in an independent cohort.

Results

A 37-RBP pairs signature consisting of 61 unique genes was constructed which was significantly associated with the survival. The RBP-related signature accurately predicted the prognosis of HCC patients, and patients in high-risk groups showed poor survival in two cohorts. The novel signature was an independent prognostic factor of HCC in two cohorts (all P < 0.001). Furthermore, the C-index of the prognostic model was 0.799, which was higher than that of many established risk models. Pathway and process enrichment analysis showed that the 61 unique genes were mainly enriched in translation, ncRNA metabolic process, RNA splicing, RNA modification, and translational termination.

Conclusion

The novel proposed RBP-related signature based on relative expression orderings could serve as a promising independent prognostic biomarker for patients with HCC, and could improve the individualized survival prediction in HCC.

A Novel Ten-Gene Signature Predicting Prognosis in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) has a dismal long-term outcome. We aimed to construct a multi-gene model for prognosis prediction to inform HCC management. The cancer-specific differentially expressed genes (DEGs) were identified using RNA-seq data of paired tumor and normal tissue. A prognostic signature was built by LASSO regression analysis. Gene set enrichment analysis (GSEA) was performed to further understand the underlying molecular mechanisms. A 10-gene signature was constructed to stratify the TCGA and ICGC cohorts into high- and low-risk groups where prognosis was significantly worse in the high-risk group across cohorts (P < 0.001 for all). The 10-gene signature outperformed all previously reported models for both C-index and the AUCs for 1-, 3-, 5-year survival prediction (C-index, 0.84 vs 0.67 to 0.73; AUCs for 1-, 3- and 5-year OS, 0.84 vs 0.68 to 0.79, 0.81 to 0.68 to 0.80, and 0.85 vs 0.67 to 0.78, respectively). Multivariate Cox regression analysis revealed risk group and tumor stage to be independent predictors of survival in HCC. A nomogram incorporating tumor stage and signature-based risk group showed better performance for 1- and 3-year survival than for 5-year survival. GSEA revealed enrichment of pathways related to cell cycle regulation among high-risk samples and metabolic processes in the low-risk group. Our 10-gene model is robust for prognosis prediction and may help inform clinical management of HCC.

Bioinformatics analysis reveals meaningful markers and outcome predictors in HBV-associated hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common type of malignant neoplasm of the liver with high morbidity and mortality. Extensive research into the pathology of HCC has been performed; however, the molecular mechanisms underlying the development of hepatitis B virus-associated HCC have remained elusive. Thus, the present study aimed to identify critical genes and pathways associated with the development and progression of HCC. The expression profiles of the GSE121248 dataset were downloaded from the Gene Expression Omnibus database and the differentially expressed genes (DEGs) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) analyses were performed by using the Database for Annotation, Visualization and Integrated Discovery. Subsequently, protein-protein interaction (PPI) networks were constructed for detecting hub genes. In the present study, 1,153 DEGs (777 upregulated and 376 downregulated genes) were identified and the PPI network yielded 15 hub genes. GO analysis revealed that the DEGs were primarily enriched in ‘protein binding’, ‘cytoplasm’ and ‘extracellular exosome’. KEGG analysis indicated that DEGs were accumulated in ‘metabolic pathways’, ‘chemical carcinogenesis’ and ‘fatty acid degradation’. After constructing the PPI network, cyclin-dependent kinase 1, cyclin B1, cyclin A2, mitotic arrest deficient 2 like 1, cyclin B2, DNA topoisomerase IIα, budding uninhibited by benzimidazoles (BUB)1, TTK protein kinase, non-SMC condensin I complex subunit G, NDC80 kinetochore complex component, aurora kinase A, kinesin family member 11, cell division cycle 20, BUB1B and abnormal spindle microtubule assembly were identified as hub genes based on the high degree of connectivity by using Cytoscape software. In addition, overall survival (OS) and disease-free survival (DFS) analyses were performed using the Gene Expression Profiling Interactive Analysis online database, which revealed that the increased expression of all hub genes were associated with poorer OS and DFS outcomes. Receiver operating characteristic curves were constructed using GraphPad prism 7.0 software. The results confirmed that 15 hub genes were able to distinguish HCC form normal tissues. Furthermore, the expression levels of three key genes were analyzed in tumor and normal samples of the Human Protein Atlas database. The present results may provide further insight into the underlying mechanisms of HCC and potential therapeutic targets for the treatment of this disease.

Identification of Potentially Therapeutic Target Genes of Hepatocellular Carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is a major threat to public health. However, few effective therapeutic strategies exist. We aimed to identify potentially therapeutic target genes of HCC by analyzing three gene expression profiles.

METHODS

The gene expression profiles were analyzed with GEO2R, an interactive web tool for gene differential expression analysis, to identify common differentially expressed genes (DEGs). Functional enrichment analyses were then conducted followed by a protein-protein interaction (PPI) network construction with the common DEGs. The PPI network was employed to identify hub genes, and the expression level of the hub genes was validated via data mining the Oncomine database. Survival analysis was carried out to assess the prognosis of hub genes in HCC patients.

RESULTS

A total of 51 common up-regulated DEGs and 201 down-regulated DEGs were obtained after gene differential expression analysis of the profiles. Functional enrichment analyses indicated that these common DEGs are linked to a series of cancer events. We finally identified 10 hub genes, six of which (OIP5, ASPM, NUSAP1, UBE2C, CCNA2, and KIF20A) are reported as novel HCC hub genes. Data mining the Oncomine database validated that the hub genes have a significant high level of expression in HCC samples compared normal samples (t-test, p < 0.05). Survival analysis indicated that overexpression of the hub genes is associated with a significant reduction (p < 0.05) in survival time in HCC patients.

CONCLUSIONS

We identified six novel HCC hub genes that might be therapeutic targets for the development of drugs for some HCC patients.

Genome-wide screening diagnostic biomarkers and the construction of prognostic model of hepatocellular carcinoma

Although methods in diagnosis and therapy of hepatocellular carcinoma (HCC) have made significant progress in decades, the overall survival (OS) of HCC remains dissatisfactory, so it is particularly important to find better diagnostic and prognostic biomarkers. In this study, we found a more reliable potential diagnostic biomarkers and constructed a more accurate prognostic evaluation model based on integrated transcriptome sequencing analysis of multiple independent data sets. First, we performed quality evaluation and differential analysis on seven Gene Expression Omnibus (GEO) data sets, and then comprehensively analyzed the differentially expressed genes with a robust rank aggregation algorithm. Next, Least absolute shrinkage and selection operator (LASSO) regression was used to establish an 8-gene prognostic risk score (RS) model. Finally, the prognostic model was further validated in the GEO data set. Also, RS has independence on other clinicopathological characteristics but has similarities in prognostic assessment compared with the T stage. Moreover, the combination of T stage and prognostic RS model based on the 8-gene had a better prognostic evaluation effect. In brief, our research suggest that the prognostic risk model of 8 genes has important clinical significance in HCC patients, and can further enrich the prognostic guidance value of the traditional T stage.

DNA methylation-driven genes for constructing diagnostic, prognostic, and recurrence models for hepatocellular carcinoma

In this study, we performed a comprehensively analysis of gene expression and DNA methylation data to establish diagnostic, prognostic, and recurrence models for hepatocellular carcinoma (HCC).

Methods: We collected gene expression and DNA methylation datasets for over 1,200 clinical samples. Integrated analyses of RNA-sequencing and DNA methylation data were performed to identify DNA methylation-driven genes. These genes were utilized in univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses to build a prognostic model. Recurrence and diagnostic models for HCC were also constructed using the same genes.

Results: A total of 123 DNA methylation-driven genes were identified. Two of these genes (SPP1 and LCAT) were chosen to construct the prognostic model. The high-risk group showed a markedly unfavorable prognosis compared to the low-risk group in both training (HR = 2.81; P < 0.001) and validation (HR = 3.06; P < 0.001) datasets. Multivariate Cox regression analysis indicated the prognostic model to be an independent predictor of prognosis (P < 0.05). Also, the recurrence model successfully distinguished the HCC recurrence rate between the high-risk and low-risk groups in both training (HR = 2.22; P < 0.001) and validation (HR = 2; P < 0.01) datasets. The two diagnostic models provided high accuracy for distinguishing HCC from normal samples and dysplastic nodules in the training and validation datasets, respectively.

Conclusions: We identified and validated prognostic, recurrence, and diagnostic models that were constructed using two DNA methylation-driven genes in HCC. The results obtained by integrating multidimensional genomic data offer novel research directions for HCC biomarkers and new possibilities for individualized treatment of patients with HCC.

Identification of a four-gene metabolic signature predicting overall survival for hepatocellular carcinoma

While hundreds of consistently altered metabolic genes had been identified in hepatocellular carcinoma (HCC), the prognostic role of them remains to be further elucidated. Messenger RNA expression profiles and clinicopathological data were downloaded from The Cancer Genome Atlas-Liver Hepatocellular Carcinoma and GSE14520 data set from the Gene Expression Omnibus database. Univariate Cox regression analysis and lasso Cox regression model established a novel four-gene metabolic signature (including acetyl-CoA acetyltransferase 1, glutamic-oxaloacetic transaminase 2, phosphatidylserine synthase 2, and uridine-cytidine kinase 2) for HCC prognosis prediction. Patients in the high-risk group shown significantly poorer survival than patients in the low-risk group. The signature was significantly correlated with other negative prognostic factors such as higher α-fetoprotein. The signature was found to be an independent prognostic factor for HCC survival. Nomogram including the signature shown some clinical net benefit for overall survival prediction. Furthermore, gene set enrichment analyses revealed several significantly enriched pathways, which might help explain the underlying mechanisms. Our study identified a novel robust four-gene metabolic signature for HCC prognosis prediction. The signature might reflect the dysregulated metabolic microenvironment and provided potential biomarkers for metabolic therapy and treatment response prediction in HCC.