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.

BackgroundIt is well known that patients with systemic lupus erythematosus (SLE) had a high risk of venous thromboembolism (VTE). This study aimed to identify the crosstalk genes between SLE and VTE and explored their clinical value and molecular mechanism initially.MethodsWe downloaded microarray datasets of SLE and VTE from the Gene Expression Omnibus (GEO) dataset. Differential expression analysis was applied to identify the crosstalk genes (CGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on the shared genes. The shared diagnostic biomarkers of the two diseases were further screened from CGs using least absolute shrinkage and selection operator (Lasso) regression. Two risk scores for SLE and VTE were constructed separately to predict the likelihood of illness according to the diagnostic biomarkers using a logical regression algorithm. The immune infiltration levels of SEL and VTE were estimated via the CIBERSORT algorithm and the relationship of CGs with immune cell infiltration was investigated. Finally, we explored potential phenotype subgroups in SLE and VTE based on the expression level of CGs through the consensus clustering method and studied immune cell infiltration in different subtypes.ResultA total of 171 CGs were obtained by the intersection of differentially expressed genes (DEGs) between SLE and VTE cohorts. The functional enrichment shown these CGs were mainly related to immune pathways. After screening by lasso regression, we found that three hub CGs (RSAD2, HSP90AB1, and FPR2) were the optimal shared diagnostic biomarkers for SLE and VTE. Based on the expression level of RSAD2 and HSP90AB1, two risk prediction models for SLE and VTE were built by multifactor logistic regression and systemically validated in internal and external validation datasets. The immune infiltration results revealed that CGs were highly correlated with multiple infiltrated immunocytes. Consensus clustering was used to respectively regroup SLE and VTE patients into C1 and C2 clusters based on the CGs expression profile. The levels of immune cell infiltration and immune activation were higher in C1 than in C2 subtypes.ConclusionIn our study, we further screen out diagnostic biomarkers from crosstalk genes SLE and VTE and built two risk scores. Our findings reveal a close relationship between CGs and the immune microenvironment of diseases. This provides clues for further exploring the common mechanism and interaction between the two diseases.