Hepatocellular carcinoma (HCC) is one of the most common types of cancers and a global health challenge with a low early diagnosis rate and high mortality. The coagulation cascade plays an important role in the tumor immune microenvironment (TME) of HCC. In this study, based on the coagulation pathways collected from the KEGG database, two coagulation-related subtypes were distinguished in HCC patients. We demonstrated the distinct differences in immune characteristics and prognostic stratification between two coagulation-related subtypes. A coagulation-related risk score prognostic model was developed in the Cancer Genome Atlas (TCGA) cohort for risk stratification and prognosis prediction. The predictive values of the coagulation-related risk score in prognosis and immunotherapy were also verified in the TCGA and International Cancer Genome Consortium cohorts. A nomogram was also established to facilitate the clinical use of this risk score and verified its effectiveness using different approaches. Based on these results, we can conclude that there is an obvious correlation between the coagulation and the TME in HCC, and the risk score could serve as a robust prognostic biomarker, provide therapeutic benefits for chemotherapy and immunotherapy and may be helpful for clinical decision making in HCC patients.

BackgroundTranscatheter arterial chemoembolization LIHC, Liver hepatocellular carcinoma; (TACE) is a valid therapeutic method for hepatocellular carcinoma (HCC). However, many patients respond poorly to TACE, thus leading to an adverse outcome. Therefore, finding new biomarkers for forecasting TACE refractoriness occurrence and prognosis becomes one of the current research priorities in the field of HCC treatment.Materials and MethodsBased on microarray datasets and a high-throughput sequencing dataset, the TACE refractoriness–related genes (TRGs) were identified by differential expression analysis. LASSO and Cox regression were applied to construct TACE refractoriness diagnostic score (TRD score) and prognostic score (TRP score) and validated their accuracy in external datasets. Functional correlation of TRP score was analyzed by gene set variation analysis and Gene Ontology. CIBERSORT and IMMUNCELL AI algorithms were performed to understand the correlation between the two scores and immune activity. We further carried out the efficacy analysis of immunotherapy and targeted drugs in the different TRP score groups. Furthermore, a nomogram was built by integrating various independent prognostic factors and validated its effectiveness in different datasets.ResultsWe identified 487 TRGs combined with GSE104580 and TCGA datasets. Then four novel TRGs (TTK, EPO, SLC7A11, and PON1) were screened out to construct TRD score and TRP score models, and both two scores had good predictive ability in external datasets. Tumors with high TRP score show an immunosuppressive phenotype with more infiltrations of regulatory T cells and macrophages. Immunotherapy and chemotherapy response evaluation revealed patients with a high TRP score demonstrated well reactions to immune checkpoint inhibitors (ICIs) and sorafenib. TRP score, TNM stage, and cancer type were brought into the combined nomogram with optimum prediction.ConclusionsOur research provided dependable and simplified methods for patients with HCC to assess tumors’ susceptibility to TACE refractoriness and prognosis and guide patients’ clinical therapy choices.

Background: Liver hepatocellular carcinoma (LIHC) remains a global health challenge with a low early diagnosis rate and high mortality. Therefore, finding new biomarkers for diagnosis and prognosis is still one of the current research priorities.Methods: Based on the variation of gene expression patterns in different stages, the LIHC-development genes (LDGs) were identified by differential expression analysis. Then, prognosis-related LDGs were screened out to construct the LIHC-unfavorable gene set (LUGs) and LIHC-favorable gene set (LFGs). Gene set variation analysis (GSVA) was conducted to build prognostic scoring models based on the LUGs and LFGs. ROC curve analysis and univariate and multivariate Cox regression analysis were carried out to verify the diagnostic and prognostic utility of the two GSVA scores in two independent datasets. Additionally, the key LCGs were identified by the intersection analysis of the PPI network and univariate Cox regression and further evaluated their performance in expression level and prognosis prediction. Single-sample GSEA (ssGSEA) was performed to understand the correlation between the two GSVA enrichment scores and immune activity.Result: With the development of LIHC, 83 LDGs were gradually upregulated and 247 LDGs were gradually downregulated. Combining with LIHC survival analysis, 31 LUGs and 32 LFGs were identified and used to establish the LIHC-unfavorable GSVA score (LUG score) and LIHC-favorable GSVA score (LFG score). ROC curve analysis and univariate/multivariate Cox regression analysis suggested the LUG score and LFG score could be great indicators for the early diagnosis and prognosis prediction. Four genes (ESR1, EHHADH, CYP3A4, and ACADL) were considered as the key LCGs and closely related to good prognosis. The frequency of TP53 mutation and copy number variation (CNV) were high in some LCGs. Low-LFG score patients have active metabolic activity and a more robust immune response. The high-LFG score patients characterized immune activation with the higher infiltration abundance of type I T helper cells, DC, eosinophils, and neutrophils, while the high-LUG score patients characterized immunosuppression with the higher infiltration abundance of type II T helper cells, TRegs, and iDC. The high- and low-LFG score groups differed significantly in immunotherapy response scores, immune checkpoints expression, and IC50 values of common drugs.Conclusion: Overall, the LIHC-progression characteristic genes can be great diagnostic and prognostic signatures and the two GSVA score systems may become promising indices for guiding the tumor treatment of LIHC patients.