Establishment and Validation of a Novel Risk Score for Hepatocellular Carcinoma Based on Bile Acid and Bile Salt Metabolism-Related Genes

Liver cancer is a public disease burden with an increasing incidence rate globally. Bile acid and bile salt’s metabolic pathways participate in liver tumorigenesis and regulate the tumor microenvironment. However, there still remains a lack of systematic analysis of the genes related to bile acid and bile salt metabolic pathways in hepatocellular carcinoma (HCC). The mRNA expression data and clinical follow-up information of patients with HCC were obtained from public databases, including The Cancer Genome Atlas, Hepatocellular Carcinoma Database, Gene Expression Omnibus, and IMvigor210. The bile acid and bile salt metabolism-related genes were extracted from Molecular Signatures Database. Univariate Cox and logistic least absolute shrinkage and selection operator regression analyses were conducted to establish the risk model. Single sample gene set enrichment analysis, Estimation of STromal and Immune cells in MAlignant Tumour tissues using Expression data, and Tumor Immune Dysfunction and Exclusion were adopted to analyze immune status. The efficiency of the risk model was tested using a decision tree and a nomogram. We determined two molecular subtypes based on bile acid and bile salt metabolism-related genes, with the prognosis of the S1 subtype being markedly superior to the S2 subtype. Next, we established a risk model based on the differentially expressed genes between the two molecular subtypes. The high-risk and low-risk groups showed significant differences in the biological pathways, immune score, immunotherapy response, and drug susceptibility. Our results demonstrated the good predictive performance of the risk model in immunotherapy datasets and established that it could be an essential factor affecting the prognosis of HCC. In conclusion, we identified two molecular subtypes based on bile acid and bile salt metabolism-related genes. The risk model established in our study could effectively predict the prognosis of patients with HCC and their immunotherapeutic response, which may contribute to targeted immunotherapy in HCC.

[1]  Yufeng Liu,et al.  Identification of a tumour immune barrier in the HCC microenvironment that determines the efficacy of immunotherapy. , 2023, Journal of hepatology.

[2]  A. Jemal,et al.  Cancer statistics, 2023 , 2023, CA: a cancer journal for clinicians.

[3]  Chen Xue,et al.  Prediction of hepatocellular carcinoma prognosis and immunotherapeutic effects based on tryptophan metabolism-related genes , 2022, Cancer Cell International.

[4]  N. Mukhopadhyay,et al.  Melanoma differentiation associated gene‐9/syndecan binding protein promotes hepatocellular carcinoma , 2022, Hepatology.

[5]  J. Ferlay,et al.  Global burden of primary liver cancer in 2020 and predictions to 2040 , 2022, Journal of hepatology.

[6]  Genhao Zhang,et al.  Endoplasmic Reticulum Stress-Related Signature for Predicting Prognosis and Immune Features in Hepatocellular Carcinoma , 2022, Journal of immunology research.

[7]  Jian Yang,et al.  Immune infiltration and clinical significance analyses of the coagulation-related genes in hepatocellular carcinoma , 2022, Briefings Bioinform..

[8]  D. Rozman,et al.  The role of bile acids in carcinogenesis , 2022, Cellular and Molecular Life Sciences.

[9]  K. Schoonjans,et al.  Metabolic Messengers: bile acids , 2022, Nature Metabolism.

[10]  P. Galle,et al.  Emerging immunotherapy for HCC: A guide for hepatologists , 2022, Hepatology.

[11]  Zhiyong Zhang,et al.  Loss of SIRT5 promotes bile acid-induced immunosuppressive microenvironment and hepatocarcinogenesis. , 2022, Journal of hepatology.

[12]  Z. Zhao,et al.  Potential biomarkers in the fibrosis progression of nonalcoholic steatohepatitis (NASH) , 2022, Journal of Endocrinological Investigation.

[13]  M. Trauner,et al.  Role of bile acids and their receptors in gastrointestinal and hepatic pathophysiology , 2022, Nature Reviews Gastroenterology & Hepatology.

[14]  Wanqing Chen,et al.  Cancer statistics in China and United States, 2022: profiles, trends, and determinants , 2022, Chinese medical journal.

[15]  W. Fang,et al.  Single-cell immune signature for detecting early-stage HCC and early assessing anti-PD-1 immunotherapy efficacy , 2022, Journal for ImmunoTherapy of Cancer.

[16]  Lanjuan Li,et al.  Crosstalk between circRNAs and the PI3K/AKT signaling pathway in cancer progression , 2021, Signal Transduction and Targeted Therapy.

[17]  V. Mazzaferro,et al.  Immunotherapies for hepatocellular carcinoma , 2021, Nature Reviews Clinical Oncology.

[18]  J. Bruix,et al.  BCLC strategy for prognosis prediction and treatment recommendation Barcelona Clinic Liver Cancer (BCLC) staging system. The 2022 update. , 2021, Journal of hepatology.

[19]  R. Quinn,et al.  Review: microbial transformations of human bile acids , 2021, Microbiome.

[20]  G. Giannelli,et al.  Epithelial–Mesenchymal Transition (EMT) Induced by TGF-β in Hepatocellular Carcinoma Cells Reprograms Lipid Metabolism , 2021, International journal of molecular sciences.

[21]  P. Zhao,et al.  Construction of TME and Identification of crosstalk between malignant cells and macrophages by SPP1 in hepatocellular carcinoma , 2021, Cancer Immunology, Immunotherapy.

[22]  Ze Li,et al.  PIVKA-II serves as a potential biomarker that complements AFP for the diagnosis of hepatocellular carcinoma , 2021, BMC cancer.

[23]  B. Sangro,et al.  Advances in immunotherapy for hepatocellular carcinoma , 2021, Nature Reviews Gastroenterology & Hepatology.

[24]  J. Zucman‐Rossi,et al.  Hepatocellular carcinoma , 1998, Nature Reviews Disease Primers.

[25]  Guangxing Zhang,et al.  Upregulated lncRNA DLX6-AS1 underpins hepatocellular carcinoma progression via the miR-513c/Cul4A/ANXA10 axis , 2020, Cancer Gene Therapy.

[26]  Wei Jia,et al.  Targeting the alternative bile acid synthetic pathway for metabolic diseases , 2020, Protein & cell.

[27]  K. Arihiro,et al.  Senescent hepatic stellate cells caused by deoxycholic acid modulates malignant behavior of hepatocellular carcinoma , 2020, Journal of Cancer Research and Clinical Oncology.

[28]  Y. Saeys,et al.  Osteopontin Expression Identifies a Subset of Recruited Macrophages Distinct from Kupffer Cells in the Fatty Liver , 2020, Immunity.

[29]  Zemin Zhang,et al.  The history and advances in cancer immunotherapy: understanding the characteristics of tumor-infiltrating immune cells and their therapeutic implications , 2020, Cellular & Molecular Immunology.

[30]  Jeffrey S. Damrauer,et al.  Comprehensive Analysis of Genetic Ancestry and Its Molecular Correlates in Cancer , 2020, Cancer Cell.

[31]  A. Lok,et al.  Utility of Liquid Biopsy Analysis in Detection of Hepatocellular Carcinoma, Determination of Prognosis, and Disease Monitoring: a Systematic Review. , 2020, Clinical gastroenterology and hepatology : the official clinical practice journal of the American Gastroenterological Association.

[32]  Zihai Li,et al.  Next-generation immuno-oncology agents: current momentum shifts in cancer immunotherapy , 2020, Journal of Hematology & Oncology.

[33]  Qiao Wu,et al.  FGF15 Activates Hippo Signaling to Suppress Bile Acid Metabolism and Liver Tumorigenesis. , 2019, Developmental cell.

[34]  Rachel S. Riley,et al.  Delivery technologies for cancer immunotherapy , 2019, Nature Reviews Drug Discovery.

[35]  J. Zucman‐Rossi,et al.  Lect2 Controls Inflammatory Monocytes to Constrain the Growth and Progression of Hepatocellular Carcinoma , 2018, Hepatology.

[36]  Xiaosheng Wang,et al.  Classification of triple-negative breast cancers based on Immunogenomic profiling , 2018, Journal of experimental & clinical cancer research : CR.

[37]  X. Liu,et al.  Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response , 2018, Nature Medicine.

[38]  Lei Chen,et al.  HCCDB: A Database of Hepatocellular Carcinoma Expression Atlas , 2018, Genom. Proteom. Bioinform..

[39]  Casey M. Theriot,et al.  Gut microbiome–mediated bile acid metabolism regulates liver cancer via NKT cells , 2018, Science.

[40]  R. Bourgon,et al.  TGF-β attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells , 2018, Nature.

[41]  D. Wasserman,et al.  The liver , 2017, Current Biology.

[42]  Pornpimol Charoentong,et al.  Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade , 2016, bioRxiv.

[43]  N. Chavez-Tapia,et al.  The Role of Bile Acids in Glucose Metabolism and Their Relation with Diabetes. , 2017, Annals of hepatology.

[44]  W. Qiu,et al.  Phenotypic and Signaling Consequences of a Novel Aberrantly Spliced Transcript FGF Receptor-3 in Hepatocellular Carcinoma. , 2016, Cancer research.

[45]  Hanns-Ulrich Marschall,et al.  Intestinal Crosstalk between Bile Acids and Microbiota and Its Impact on Host Metabolism. , 2016, Cell metabolism.

[46]  G. Chen,et al.  Preoperative serum α‐fetoprotein and prognosis after hepatectomy for hepatocellular carcinoma , 2016, The British journal of surgery.

[47]  Andrea Bild,et al.  Alternative preprocessing of RNA-Sequencing data in The Cancer Genome Atlas leads to improved analysis results , 2015, Bioinform..

[48]  G. Getz,et al.  Inferring tumour purity and stromal and immune cell admixture from expression data , 2013, Nature Communications.

[49]  J. Boyer,et al.  Bile formation and secretion. , 2013, Comprehensive Physiology.

[50]  M. Trauner,et al.  Bile acid transporters and regulatory nuclear receptors in the liver and beyond , 2013, Journal of hepatology.

[51]  Sean R. Davis,et al.  NCBI GEO: archive for functional genomics data sets—update , 2012, Nucleic Acids Res..

[52]  Justin Guinney,et al.  GSVA: gene set variation analysis for microarray and RNA-Seq data , 2013, BMC Bioinformatics.

[53]  Helga Thorvaldsdóttir,et al.  Molecular signatures database (MSigDB) 3.0 , 2011, Bioinform..

[54]  S. Rittling Osteopontin in macrophage function , 2011, Expert Reviews in Molecular Medicine.

[55]  A. Bhattacharyya,et al.  Carcinogenicity of deoxycholate, a secondary bile acid , 2011, Archives of Toxicology.

[56]  H. Imamura,et al.  AFP, AFP-L3, DCP, and GP73 as markers for monitoring treatment response and recurrence and as surrogate markers of clinicopathological variables of HCC , 2010, Journal of Gastroenterology.

[57]  Ben S. Wittner,et al.  Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1 , 2009, Nature.

[58]  H. Bernstein,et al.  Bile acids as endogenous etiologic agents in gastrointestinal cancer. , 2009, World journal of gastroenterology.

[59]  Kenji Ikeda,et al.  Gene expression in fixed tissues and outcome in hepatocellular carcinoma. , 2008, The New England journal of medicine.

[60]  T. Rikiyama,et al.  Bile acids repress E‐cadherin through the induction of Snail and increase cancer invasiveness in human hepatobiliary carcinoma , 2008, Cancer science.

[61]  H. Shin,et al.  SPP1 polymorphisms associated with HBV clearance and HCC occurrence. , 2007, International journal of epidemiology.

[62]  A. Hofmann Biliary secretion and excretion in health and disease: current concepts. , 2007, Annals of hepatology.

[63]  Mary Ann Comunale,et al.  GP73, a resident Golgi glycoprotein, is a novel serum marker for hepatocellular carcinoma. , 2005, Journal of hepatology.

[64]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[65]  K. Zatloukal,et al.  Hepatobiliary transporter expression in human hepatocellular carcinoma , 2005, Liver international : official journal of the International Association for the Study of the Liver.

[66]  H. Bernstein,et al.  Bile acids as carcinogens in human gastrointestinal cancers. , 2005, Mutation research.

[67]  R. Tibshirani Regression Shrinkage and Selection via the Lasso , 1996 .