Immune- and metabolism-related gene signature analysis uncovers the prognostic and immune microenvironments of hepatocellular carcinoma.

Metabolic reprogramming is an emerging hallmark that influences the tumour microenvironment (TME) by regulating the behavior of cancer cells and immune cells. The relationship between metabolism and immunity remains elusive. The purpose of this study was to explore the predictive value of immune- and metabolism-related genes in hepatocellular carcinoma (HCC) and their intricate interplay with TME. We established the immune- and metabolism-related signature (IMRPS) based on the LIHC cohort from The Cancer Genome Atlas (TCGA) dataset. Kaplan-Meier analysis, receiver operating characteristic (ROC) curve analysis and Cox regression analysis confirmed the prognostic value of IMRPS. We investigated differences in immune cell infiltration, clinical features, and therapeutic response between risk groups. The quantitative real-time PCR (qPCR) was used to confirm the expression of signature genes. Immunohistochemical staining was performed to evaluate immune infiltration features in HCC tissue samples. We conducted cell experiments including gene knockout, cell counting kit-8 (CCK-8), and flow cytometry to explore the role of the IMRPS key gene UCK2 in HCC. RNA-seq was used to further investigate the potential underlying mechanism involved. The IMRPS, composed of four genes, SMS, UCK2, PFKFB4 and MAPT, exhibited significant correlations with survival, immune cell infiltration, clinical features, immune checkpoints and therapeutic response. The IMRPS was shown to be an excellent predictor of HCC prognosis. It could stratify patients appropriately and characterize the TME accurately. The high-risk HCC group exhibited an immunosuppressive microenvironment with abundant M2-like macrophage infiltration, which was confirmed by the immunohistochemistry results. The results of qPCR revealed that the expression of signature genes in 20 HCC tissues was significantly greater than that in adjacent normal tissues. After the key gene UCK2 was knocked out, the proliferation of the Huh7 cell line was significantly inhibited, and monocyte-derived macrophages polarized towards an M1-like phenotype in the coculture system. RNA-seq and GSEA suggested that the phenotypes were closely related to the negative regulation of growth and regulation of macrophage chemotaxis. This study established a new IMRS for the accurate prediction of patient prognosis and the TME, which is also helpful for identifying new targets for the treatment of HCC.

Gu Y ，Ma E ，Jiang S ，Shan Z ，Xia G ，Ma R ，Fu J ，Wang Z ... -  《-》

Predictive value of a stemness-based classifier for prognosis and immunotherapy response of hepatocellular carcinoma based on bioinformatics and machine-learning strategies.

Significant advancements have been made in hepatocellular carcinoma (HCC) therapeutics, such as immunotherapy for treating patients with HCC. However, there is a lack of reliable biomarkers for predicting the response of patients to therapy, which continues to be challenging. Cancer stem cells (CSCs) are involved in the oncogenesis, drug resistance, and invasion, as well as metastasis of HCC cells. Therefore, in this study, we aimed to create an mRNA expression-based stemness index (mRNAsi) model to predict the response of patients with HCC to immunotherapy. We retrieved gene expression and clinical data of patients with HCC from the GSE14520 dataset and the Cancer Genome Atlas (TCGA) database. Next, we used the "one-class logistic regression (OCLR)" algorithm to obtain the mRNAsi of patients with HCC. We performed "unsupervised consensus clustering" to classify patients with HCC based on the mRNAsi scores and stemness subtypes. The relationships between the mRNAsi model, clinicopathological features, and genetic profiles of patients were compared using various bioinformatic methods. We screened for differentially expressed genes to establish a stemness-based classifier for predicting the patient's prognosis. Next, we determined the effect of risk scores on the tumor immune microenvironment (TIME) and the response of patients to immune checkpoint blockade (ICB). Finally, we used qRT-PCR to investigate gene expression in patients with HCC. We screened CSC-related genes using various bioinformatics tools in patients from the TCGA-LIHC cohort. We constructed a stemness classifier based on a nine-gene (PPARGC1A, FTCD, CFHR3, MAGEA6, CXCL8, CABYR, EPO, HMMR, and UCK2) signature for predicting the patient's prognosis and response to ICBs. Further, the model was validated in an independent GSE14520 dataset and performed well. Our model could predict the status of TIME, immunogenomic expressions, congenic pathway, and response to chemotherapy drugs. Furthermore, a significant increase in the proportion of infiltrating macrophages, Treg cells, and immune checkpoints was observed in patients in the high-risk group. In addition, tumor cells in patients with high mRNAsi scores could escape immune surveillance. Finally, we observed that the constructed model had a good expression in the clinical samples. The HCC tumor size and UCK2 genes expression were significantly alleviated and decreased, respectively, by treatments of anti-PD1 antibody. We also found knockdown UCK2 changed expressions of immune genes in HCC cell lines. The novel stemness-related model could predict the prognosis of patients and aid in creating personalized immuno- and targeted therapy for patients in HCC.

Chen E ，Zou Z ，Wang R ，Liu J ，Peng Z ，Gan Z ，Lin Z ，Liu J ... -  《Frontiers in Immunology》

A Cuproptosis-Related LncRNA Risk Model for Predicting Prognosis and Immunotherapeutic Efficacy in Patients with Hepatocellular Carcinoma.

Cuproptosis is a novel programmed cell death pathway that is initiated by direct binding of copper to lipoylated tricarboxylic acid (TCA) cycle proteins. Recent studies have demonstrated that cuproptosis-related genes regulate tumorigenesis. However, the potential role and clinical significance of cuproptosis-related long noncoding RNAs (lncRNAs) in hepatocellular carcinoma (HCC) have not been established. We performed a bioinformatics analyses of RNA-sequencing data of HCC patients extracted from The Cancer Genome Atlas (TCGA) dataset to identify and validate a cuproptosis-related lncRNA prognostic signature. Furthermore, we analyzed the clinical significance of the prognostic signature of cuproptosis-related lncRNA in predicting the immunotherapeutic efficacy and the status of the tumor immune microenvironment. The RNA-sequencing data, genomic mutations, and clinical information were downloaded for 374 HCC samples and 50 normal liver samples from TCGA-Liver Hepatocellular Carcinoma (TCGA-LIHC) dataset. Co-expression analysis of Gene-lncRNA pairs with 49 known cuproptosis-related prognostic genes was used to define cuproptosis-related prognostic lncRNAs. We performed the LASSO algorithm and univariate and multivariate Cox regression analysis, respectively, to gradually identify the prognostic risk models of cuproptosis-related lncRNA based on the TCGA-LIHC dataset. Subsequently, the predictive performance of the model was evaluated using receiver operation characteristic (ROC) curves, Kaplan-Meier survival curves, and prognostic nomogram. The analysis of gene-lncRNA co-expression with 49 known cuproptosis-related genes identified 1359 cuproptosis-related lncRNAs in the TCGA-LIHC data set. A prognostic model was constructed with nine cuproptosis-related prognostic lncRNAs (AC007998.3, AC003086.1, AC009974.2, IQCH-AS1, LINC0256 1, AC105345.1, ZFPM2-AS1, AL353708.1 and WAC-AS1) using LASSO regression and Cox regression analyses. Risk scores were calculated for all HCC patient samples based on the four cuproptosis-related lncRNA prognostic models. All HCC patients were divided into high-risk and low-risk subgroups according to a 1:1 ratio column. The Kaplan-Meier survival curve analysis showed that the overall survival rate (OS) of the high-risk group patients was significantly lower than that of the low-risk group. The principal component analysis (PCA) confirmed that the prognostic lncRNA model accurately distinguished between high- and low-risk HCC patients. Furthermore, regression analysis as well as ROC curves confirmed the prognostic value of the risk score. A nomogram with risk scores and other clinicopathological characteristics was constructed. The nomogram accurately predicted the probability of 1-, 3-, and 5-year OS in HCC patients. Tumor mutation burden (TMB) scores were higher for high-risk patients than for low-risk patients. HCC patients in the low-risk group showed lower TIDE scores and greater sensitivity to antitumor drugs than those in the high-risk group. Tumor immune responses and tumor immune cell infiltration were significantly different between the high-risk and low-risk groups of patients with HCC. Our study identified a 9-cuproptosis-related lncRNA signature that accurately predicted prognosis, immunotherapeutic efficacy, and the status of the tumor immune microenvironment in HCC patients. Therefore, this cuproptosis-related lncRNA risk model is a potential prognostic biometric feature in HCC and shows high clinical value in identifying HCC patients who are potentially responsive to immunotherapy.

Wang S ，Bai H ，Fei S ，Miao B ... -  《-》

Histamine-related genes participate in the establishment of an immunosuppressive microenvironment and impact the immunotherapy response in hepatocellular carcinoma.

Chronic inflammation is pivotal in the pathogenesis of hepatocellular carcinoma (HCC). Histamine is a biologically active substance that amplifies the inflammatory and immune response and serves as a neurotransmitter. However, knowledge of histamine's role in HCC and its effects on immunotherapy remains lacking. We focused on histamine-related genes to investigate their potential role in HCC. The RNA-seq data and clinical information regarding HCC were obtained from The Cancer Genome Atlas (TCGA). After identifying the differentially expressed genes, we constructed a signature using the univariate Cox proportional hazard regression and least absolute shrinkage and selection operator (LASSO) analyses. The signature's predictive performance was evaluated using a receiver operating characteristic curve (ROC) analysis. Furthermore, drug sensitivity, immunotherapy effects, and enrichment analyses were conducted. Histamine-related gene expression in HCC was confirmed using quantitative real-time polymerase chain reaction (qRT-PCR). A histamine-related gene prognostic signature (HRGPS) was developed in TCGA. Time-dependent ROC and Kaplan-Meier survival analyses demonstrated the signature's strong predictive power. Importantly, patients in high-risk groups exhibited a higher frequency of TP53 mutations, elevated immune checkpoint-related gene expression, and increased infiltration of immunosuppressive cells-indicating a potentially favorable response to immunotherapy. In addition, drug sensitivity analysis revealed that the signature could effectively predict chemotherapy efficacy and sensitivity. qRT-PCR results validated histamine-related gene overexpression in HCC. Our findings demonstrate that inhibiting histamine-related genes and signaling pathways can impact the therapeutic effect of anti-PD-1/PD-L1. The precise predictive ability of our signature in determining the response to different therapeutic options highlights its potential clinical significance.

Zhang X ，Zheng P ，Meng B ，Zhuang H ，Lu B ，Yao J ，Han F ，Luo S ... -  《-》

Construction of a risk score prognosis model based on hepatocellular carcinoma microenvironment.

Hepatocellular carcinoma (HCC) is a common cancer with a poor prognosis. Previous studies revealed that the tumor microenvironment (TME) plays an important role in HCC progression, recurrence, and metastasis, leading to poor prognosis. However, the effects of genes involved in TME on the prognosis of HCC patients remain unclear. Here, we investigated the HCC microenvironment to identify prognostic genes for HCC. To identify a robust gene signature associated with the HCC microenvironment to improve prognosis prediction of HCC. We computed the immune/stromal scores of HCC patients obtained from The Cancer Genome Atlas based on the ESTIMATE algorithm. Additionally, a risk score model was established based on Differentially Expressed Genes (DEGs) between high- and low-immune/stromal score patients. The risk score model consisting of eight genes was constructed and validated in the HCC patients. The patients were divided into high- or low-risk groups. The genes (Disabled homolog 2, Musculin, C-X-C motif chemokine ligand 8, Galectin 3, B-cell-activating transcription factor, Killer cell lectin like receptor B1, Endoglin and adenomatosis polyposis coli tumor suppressor) involved in our risk score model were considered to be potential immunotherapy targets, and they may provide better performance in combination. Functional enrichment analysis showed that the immune response and T cell receptor signaling pathway represented the major function and pathway, respectively, related to the immune-related genes in the DEGs between high- and low-risk groups. The receiver operating characteristic (ROC) curve analysis confirmed the good potency of the risk score prognostic model. Moreover, we validated the risk score model using the International Cancer Genome Consortium and the Gene Expression Omnibus database. A nomogram was established to predict the overall survival of HCC patients. The risk score model and the nomogram will benefit HCC patients through personalized immunotherapy.

Zhang FP ，Huang YP ，Luo WX ，Deng WY ，Liu CQ ，Xu LB ，Liu C ... -  《-》