Downregulated Acetyl-CoA Acyltransferase 2 Promoted the Progression of Hepatocellular Carcinoma and Participated in the Formation of Immunosuppressive Microenvironment

Abstract

Background

The aim of this study is to explore the role of acetyl-CoA acyltransferase 2 (ACAA2) in the progression of hepatocellular carcinoma (HCC).

Methods

Bulk RNA data and single-cell RNA data were acquired from The Cancer Genome Atlas and Gene Expression Omnibus. Both in vitro and in vivo studies were used to determine the effect of ACAA2 on the progression of HCC, and RNA sequencing analysis was performed to explore the mechanism.

Results

We found downregulation of ACAA2 was involved in the malignant progression of HCC. The patient with low ACAA2 level had an immunosuppressive microenvironment in the HCC and predicted to have a poor prognosis. Decreased ACAA2 facilitated HCC proliferation and metastasis by activating the nuclear factor-κB (NFκB) signaling pathway. And increased CXCL1 induced by NFκB signaling pathway might be responsible for low level of ACAA2 related immunosuppressive microenvironment. Furthermore, the expression of ACAA2 was also detected in immune cells. The expression of ACAA2 in CD4⁺TCF7⁺T, CD4⁺FOXP3⁺T, CD8⁺GZMK⁺T, and CD8⁺KLRD1⁺T cells was inversely correlated with the composition of CD8⁺PDCD1⁺T cells in HCC. This effect might be due to the CCL5-CCRs and HLA-E-KLRCs ligand-receptor networks.

Conclusion

In a conclusion, downregulated ACAA2 promoted the progression of hepatocellular carcinoma and might be participated in the formation of immunosuppressive microenvironment. ACAA2 could be served as a favorable indicator for the prognosis of HCC and an ideal biomarker for immunotherapy.

Keywords:

• hepatocellular carcinoma
• acetyl-CoA acyltransferase 2
• immunotherapy
• nuclear factor-κB signaling pathway
• biomarker

Introduction

Hepatocellular carcinoma (HCC) remains a global challenge, accounting for 90% of liver malignancies.¹ Owing to late diagnosis and limited treatment, the overall survival (OS) rate of HCC is unsatisfactory, and HCC ranks third in cancer-related deaths worldwide.² Until 2017, only tyrosine kinase inhibitors (TKIs) were used for treating advanced-stage HCC.³ In the past few years, immune checkpoint blockades (ICBs) have exhibited remarkable effects during cancer therapy and have become a first-line option for patients with advanced HCC who fail TKIs, based on the last version of the National Comprehensive Cancer Network (NCCN Guidelines Version 1.2022 Hepatocellular Carcinoma). ICBs are monoclonal antibodies that block the interaction between checkpoint proteins and their ligands to improve the immune response in HCC.⁴ The objective remission rates of ICBs in HCC are 15–20% (including 1–5% complete responses) in a cohort of patients with HCC treated with a single agent.⁵ Responders to ICBs can gain durable benefits and prolong survival.⁵ In addition, approximately 20% of patients receiving pembrolizumab remained free from disease progression for more than a year compared to less than 7% of controls.⁶

Metabolism reprogramming has thought to play key role in the progression of HCC.^7,8 Recent studies have also revealed that oncogene-induced metabolic reprogramming in tumor cells can affect the tumor microenvironment to limit immune responses and present barriers to immunotherapy.⁹ The poor vascular condition and demand for uncontrolled proliferation forces tumor cells to compete with immune cells for the necessary nutrients.¹⁰ Little is known about the link between the metabolic features of tumor cells and immunotherapy.

In this study, we aimed to explore underlying metabolic biomarkers and their biological functions in HCC. Acetyl-CoA acyltransferase 2 (ACAA2), which decreases during HCC initiation, was identified as a gene of interest.¹¹ ACAA2 is involved in mitochondrial fatty acid elongation and degradation by catalyzing the last step of the respective β-oxidation pathway.¹² Decreased ACAA2 is essential for the development of hepatic steatosis by diminishing hepatic fatty acid β-oxidation.¹³ Besides, upregulated ACAA2 could also diminish the hepatotoxicity induced by acetaminophen, indicating the liver protect effect of ACAA2.¹⁴ However, there is fewer studies to demonstrate the role of ACAA2 in the HCC.

In the present study, we aimed to investigate the role of ACAA2 in HCC progression and subsequently, we evaluated its potential role in the tumor microenvironment.

Materials and Methods

Datasets and Cell Lines

Data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) were acquired using the SangerBox platform (http://vip.sangerbox.com/). Mutant-allele tumor heterogeneity (MATH) and stemness scores were obtained from the SangerBox platform. All the cell lines were purchased from Shanghai Cell Bank of the Chinese Academy of Sciences. All cell lines were cultured in Dulbecco’s modified Eagle medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and 1% antibiotics (100 U/mL penicillin and 100 µg/mL streptomycin). All the cell lines were cultured at 37 °C in a 5% CO₂ incubator.

Determination of Immune Cell Infiltration Levels and Immune Related Functional Scores

Immune cell infiltration level was obtained from the tumor and immune system interaction database (TISIDB), an integrated repository portal for tumor-immune system interactions.¹⁵ Tumor immune dysfunction and exclusion (TIDE) was generated as previously described.^16,17 Single-cell RNA analysis was performed using the GSE125449 and GSE140228 datasets. The cell types were identified using tumor immune single-cell hub (TISCH) (http://tisch.comp-genomics.org/home/) and the scTIME portal (http://sctime.sklehabc.com/unicellular/home).

Lentivirus, Stable Cell Line Construction, and Quantitative Real-Time PCR (qPCR)

The lentiviral vector system (LV) and empty vector was purchased from the GeneChem Corporation (Shanghai, China). Stable cell lines expressing the target gene or negative control were selected by adding 0.5 μg/mL puromycin into the medium. Total RNA was extracted and quantified using an RNA purification kit (Thermo Fisher, US) according to the manufacturer’s protocol. Quantitative real-time PCR was performed using the SYBR green power master mix (Promega, Madison, WI, USA). Primers used for the qRT-PCR experiments are listed:

F: 5′ - CTGCTCCGAGGTGTGTTTGTA-3′,

R: 5′ - GGCAGCAAATTCAGACAAGTCA-3′ for ACAA2;

F: 5′ - CTCGAGGCCCCTGGGGCAGAAGCCTC -3′,

R: 5′- GATATCGGGGCTCAGCAGGCGGGTCT -3′ for CXCL1;

F: 5′-AGAAGGCTGGGGCTCATTTG-3′,

R: 5′-AGGGGCCATCCACAGTCTTC-3′ for GAPDH.

Cell Proliferation and Cell Cycle Assay

For cell proliferation, cells were seeded in 96-well plates at 2×10³ cells per well, and cell viability was measured using a cell counting kit-8 (CCK-8) assay (Dojindo Molecular Technologies, Inc., Japan) at different time points. For cell cycle tests, 4×10⁴ cells were stained according to the protocol of the Cycle TESTTM PLUS DNA reagent kit (BD Biosciences San Jose, CA) and analyzed using flow cytometry (Beckman Coulter FC500).

Migration and Invasion

The migration and invasion of HCC cells were determined using Transwell assays. Chambers with or without precoated Matrigel were placed in a 24-well plate. DMEM supplemented with 10% FBS was added to the lower chamber. Cells (5×10⁴) suspended in 100 µL of FBS-free DMEM were seeded in the upper chamber and cultured for 24 h.

Western Blot

Cells were lysed in radioimmunoprecipitation assay buffer containing protease and phosphatase inhibitors. The samples were separated using SDS–PAGE and transferred to polyvinylidene fluoride membranes. After blocking with 5% skimmed milk at room temperature for 1 h, the membranes were incubated with primary antibodies overnight. The membranes were then incubated with the secondary antibodies for 1 h at room temperature, and the bands were detected using an enhanced chemiluminescence kit. The primary antibodies used in this study were ACAA2 (ab128929, ABCAM), pNFKB (3033, Cell Signaling Technology), and NFKB (8242, Cell Signaling Technology).

In vivo Tumor Growth

All experimental protocols involving animals were approved by the Animal Ethics Committee of the First Affiliated Hospital of the Harbin Medical University (2019049). Four–six-week-old male BALB/c nude mice were purchased from Shanghai Animal Centre (Shanghai, China). Flank subcutaneous xenografts were established by subcutaneous injection of 1×10⁶ cells suspended in 100 μL of phosphate-buffered saline. After 6 weeks, tumor growth was monitored using bioluminescence signals. For the liver orthotopic xenograft implantation model, subcutaneous tumors were dissected into 1 mm³ sections and then implanted into the liver parenchyma. Tumor growth was monitored using bioluminescence signals, and the mice were sacrificed after 6 weeks.

Statistical Analysis

The Student’s t-test was used for comparisons between the two groups. Pearson’s correlation analysis was used to determine the linear relationship between the two groups. Multivariate analysis was performed using the Cox regression model. Kaplan Meier curves were used to compare survival, and the Log rank test was used to compare survival between the different groups. The threshold was defined as P <0.05.

Results

Downregulated ACAA2 Was Associated with HCC Progression

ACAA2 was frequently downregulated in different tumor tissues, including HCC (Figure 1A). Furthermore, ACAA2 gradually decreased during HCC initiation (Figure 1B). Patients with lower levels of ACAA2 were predicted to have a larger tumor size, more vascular invasion, higher tumor grade, and higher inflammation levels (Figure 1C–F). Gene set enrichment analysis (GSEA) results revealed that the cell cycle-related gene set was enriched in patients with low levels of ACAA2 (Figure 1G). The MATH value was used to determine HCC heterogeneity. The expression of ACAA2 was negatively associated with MATH, implying that the heterogeneity of HCC gradually increased as ACAA2 decreased (Figure 1H). The level of stemness of HCC was determined by RNA-based stemness scores. The results showed that the gene expression of ACAA2 was negatively associated with the stemness level of HCC (Figure 1I). The above results showed that downregulation of ACAA2 was involved in the malignant progression of HCC.

Figure 1 Downregulated ACAA2 was associated with HCC progression. (A) Gene expression of ACAA2 in different tumor types; blue: normal tissues; red: tumor tissues; (B) Gene expression of ACAA2 in different disease stages, data derived from GSE6764; (C–F) Compared gene expression of ACAA2 between different groups, data derived from TCGA-LIHC dataset; (G) GSEA data revealed the correlation between gene expression of ACAA2 and cell cycle; (H) The correlation between the expression of ACAA2 and MATH score; (I) The correlation between the expression of ACAA2 and stemness. **p<0.01; ***p<0.001; ****p<0.0001.

ACAA2 Served as a Biomarker for Immunotherapy

We then investigated the association between ACAA2 and the tumor microenvironment (TME). We found that ACAA2 was positively correlated with the infiltration of natural killer (NK) cells in HCC tissues but was inversely associated with regulatory T cells (Tregs) and myeloid-derived suppressor cells (MDSCs)(Figure 2A). In addition, ACAA2 was negatively correlated with immune checkpoints such as PDCD1, CTLA4, TIGHT, and LAG3 (Figure 2B). These results demonstrate that a low level of ACAA2 is accompanied by an immunosuppressive microenvironment. A TIDE score was generated to assess the response of the tumor to immunotherapy. Linear regression analysis revealed a reverse correlation between the expression of ACAA2 and TIDE score (Figure 2C). Moreover, patients with lower ACAA2 levels were predicted to have a poor response to ICBs, based on TIDE scores (Figure 2D). Patients with lower levels of ACAA2 were predicted to have shorter overall survival (OS) time and disease-free survival (DFS) time (Figure 2E–F). Immune checkpoints that were significantly correlated with ACAA2 were filtered. The risk scores of the ACAA2 related immune checkpoints were generated using a stepwise multivariate Cox regression analysis. The patient with a high-risk score was predicted to have a poor prognosis (Figure 2G). The combined effect of ACAA2 and risk score was generated using Kaplan-Meier curves. In the low-risk score group, the expression of ACAA2 did not influence the survival time, whereas in the high-risk score group, patients with low ACAA2 expression had a shorter OS (Figure 2H). The above data suggest that patients with low levels of ACAA2 had a low response to ICBs, and ACAA2 could be a biomarker for immunotherapy.

Figure 2 ACAA2 served as a biomarker for immunotherapy. (A) The correlation between ACAA2 and immune cell types and MHC; (B) Correlation between ACAA2 and immune check points; (C) Correlation between ACAA2 and TIDE score; (D) Compared the expression of ACAA2 between different groups; (E) Kaplan-Meier curves were used to determine the overall survival between different groups; (F) Kaplan-Meier curves were used to determine the disease-free survival between different groups; (G) Kaplan-Meier curves were used to determine the prognostic value of risk score; (H) Kaplan-Meier curves were used to determine the combine effect of ACAA2 and risk score on the prognosis of HCC; *p<0.05; **p<0.01; ****p<0.0001.

Abbreviation: RS, risk score.

Downregulated ACAA2 Promoted Proliferation and Metastasis of HCC

To determine the effect of ACAA2 on HCC, its mRNA expression levels were evaluated in normal (THLE-3) and different HCC cell lines. The HCC cell lines exhibited a relatively lower expression of ACAA2 compared with that noted for THLE-3 (Figure 3A). Subsequently, a stably transfected cell line was established by overexpressing ACAA2 in Huh7 cells, which had a relatively lower expression of ACAA2. In addition, silencing of ACAA2 was performed in the HCCLM3 cell line, which exhibited relatively higher expression of ACAA2. The efficiency of LV transfection was assessed using real-time PCR and Western blotting, and Lv-shRNA-3 (59%) was selected as the best candidate for silencing ACAA2 expression (Figure 3B–D). Elevated ACAA2 suppressed proliferation and induced cell cycle arrest in Huh7 cells (Figure 3E and F). Knockdown of ACAA2 accelerated proliferation and promoted the cell cycle process in the HCCLM3 cell line (Figure 3G and H). Overexpression of ACAA2 inhibited the migration and invasion ability of Huh7 cells, and decreased ACAA2 facilitated the migration and invasion of the HCCLM3 cell line (Figure 3I and J). The in vivo results were consistent with those of in vitro studies. ACAA2 knockdown resulted in larger tumor volumes in the subcutaneous tumor model (Figure 3K, Supplementary Figure 1). Decreased ACAA2 levels also promoted tumor growth in the xenograft model (Figure 3L). These results demonstrated that downregulation of ACAA2 promotes the proliferation and metastasis of HCC.

Figure 3 Downregulated ACAA2 promoted proliferation and metastasis of HCC. (A) Gene expression of ACAA2 in different cell lines; (B and C) Gene expression of ACAA2 in the indicated cell lines; (D) The protein expression of ACAA2 was examined by Western blot in the indicated cell lines; (E and G) CCK8 assays was used to determine the proliferation of the indicated cell lines and the value of OD450 was related to control; (F and H) Representative images for cell cycle assays in indicated cell lines in the left panel, and statistical analysis in the right panel. (I and J) The migration and invasion abilities of the indicated cell lines were detected by Transwell assays in the left panel, and the statistical analysis was placed in the right panel. Bar=50µm (K) Representative images of subcutaneous tumors derived from indicated cell lines in left panel; the statistical analysis of bioluminescence signal in the right panel; (L) Representative images of liver xenograft derived from indicated cell lines in left panel; the statistical analysis of bioluminescence signal in the right panel. Five mice were used for each group *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001.

Downregulated ACAA2 Facilitated HCC Progression via the NF-κB Signaling Pathway

To further explore the mechanism by which of downregulated ACAA2 promoting HCC, an mRNA sequencing analysis was performed. Differentially expressed genes (|fold change| > 2, p < 0.05) were filtered and enriched in the signaling pathway (Figure 4A). KEGG enrichment analysis revealed that dysregulated signaling pathways induced by ACAA2 knockdown were majorly related to human diseases and organismal systems (Figure 4B). Hallmark enrichment analysis showed that downregulation of ACAA2 increased the adaptation of HCC to hypoxia and promoted epithelial mesenchymal transition (Figure 4C). In addition, diminished ACAA2 could activate the NF-κB signaling pathway (Figure 4C). GO enrichment analysis suggested that reduced ACAA2 promoted the biosynthesis of macromolecules and RNA, intracellular delivery, and transcription (Figure 4D–F). Activation of the NF-κB signaling pathway was examined using Western blotting. The results showed that elevated ACAA2 decreased the expression of pNFκB, whereas decreased ACAA2 increased the level of pNFκB (Figure 4G). The expression variation of ACAA2 had no effect on total NFκB (Figure 4G). BAY 11–7082 was used to inhibit NF-κB signaling. The activation of the NF-κB signaling pathway by downregulating ACAA2 was significantly suppressed by 5 μM BAY 11–7082 (Figure 4H). Blocking the NF-κB signaling pathway attenuated the proliferation induced by ACAA2 knockdown (Figure 4I). These results revealed that downregulation of ACAA2 facilitated HCC progression via the NF-κB signaling pathway.

Figure 4 Downregulated ACAA2 facilitated HCC progression via the NF-κB signaling pathway. (A) Volcano plot exhibited the gene expression altered after ACAA2 knockdown; (B) Different expression genes were enriched in KEGG pathways; (C) Different expression genes were enriched in Hallmark pathways; (D–F) Different expression genes were enriched in GO pathway; (G and H) The protein expression of pNKFB and NFKB were examined by Western blot in the indicated cell lines. 7.5µM BAY11-7082 was used to inhibit the NFKB signaling pathway; (I) CCK8 assays was used to determine the proliferation of the indicated cell lines and the value of OD450 was related to control; **p<0.01; ***p<0.001.

Abbreviations: BP, biological process; CC, cellular component; MF, molecular function.

CXCL1 Contributed to the Low Level of ACAA2 Related Immune Suppression Microenvironment

Downregulation of ACAA2 increased the expression of CXCL1 via the NF-κB signaling pathway (Figure 5A). CXCL1, a secreted growth factor, plays a key role in tumor microenvironment modification. The GSEA results demonstrated that ACAA2 and CXCL1 were inversely involved in several metabolic signaling pathways in HCC (Figure 5B). Single-cell RNA (scRNA) analysis was performed using GSE125449. The cell types were identified using TISCH (Figure 5C). Both ACAA2 and CXCL1 were highly expressed in malignant cells (Figure 5D). The ACAA2–CXCL1 axis demonstrated that metabolic reprogramming of tumor cells could modify the tumor microenvironment in a paracrine manner. High level of CXCL1 in HCC tissues were positively correlated with the infiltration of MDSCs, neutrophils, and Tregs as well as the expression of CTLA4, but not PDCD1 (Figure 5E). The patient with a high CXCL1 level showed a poor response to immunotherapy (Figure 5F). Linear regression analysis revealed a reverse correlation between ACAA2 and CXCL1 expression in HCC (Figure 5G). These results suggest that CXCL1 is responsible for the low level of ACAA2 related immune suppression microenvironment.

Figure 5 CXCL1 contributed to the low level of ACAA2 related immune suppression microenvironment. (A) Gene expression of CXCL1 in the indicated cell lines in the upper panel; protein expression of CXCL1 in the indicated cell lines in the lower panel; (B) The top 10 common signaling pathways that both ACAA2 and CXCL1 were involved in; (C) The cell types were identified; single-cell RNA (scRNA) data was derived from GSE125449; (D) The expression of ACAA2 and CXCL1 in different cell types; (E) The correlation between ACAA2 and immune cell types and immune check points; (F) Compared the expression of ACAA2 in different groups; (G) Linear regression analysis was used to determine the relationship between ACAA2 and CXCL1 in HCC; data was derived from TCGA-LIHC. **p<0.01; ***p<0.001; ****p<0.0001.

Exploring the Role of ACAA2 in the Immune Cells

The scRNA data were derived from GSE140228, and the cell types were identified using the scTIME portal. We found that ACAA2 was also gently expressed in the immune cells (Figure 6A). Basing on the expression of ACAA2, the immune cell types were divided into ACAA2 positive group and ACAA2 negative group (Figure 6B). CD8⁺PDCD1⁺T cells belonged to the ACAA2 negative group (Figure 6B). A ligand-receptor (LR) network among the immune cell types was generated, and 28 cell types interacted with CD8+PDCD1+T cells including 12 ACAA2 positive cell types (Figure 6C and D). Among these cell types, the expression of ACAA2 in CD4⁺TCF7⁺T, CD4⁺FOXP3⁺T, CD8⁺GZMK⁺T, and CD8⁺KLRD1⁺T cells was strongly negatively correlated with the composition of CD8⁺PDCD1⁺T cells (Figure 6E). The top 10 LR pairs are shown (Figure 6F). All four cell types existed CCL5–CCRs and HLA-E-KLRCs LR pairs with CD8+PDCD1+ T cells. These results suggested that ACAA2 in immune cells decreased the composition of CD8⁺PDCD1⁺T cells by affecting the LR network. The lack of ACAA2 receded the antitumor immune response and increased resistance to immunotherapy.

Figure 6 Exploring the role of ACAA2 in the immune cells. (A) The cell types were identified; single-cell RNA (scRNA) data was derived from GSE140228; (B) The expression of ACAA2 in different cell types; (C) Ligand-receptor networks among the cell types; (D) Ligand-receptor networks between CD8-C8-PDCD1 cell type and other cell types; (E) The correlation between the expression of ACAA2 in the four cell types and composition of CD8-CD-PDCD1 cell type in HCC; (F) The top 10 ligand-receptor networks for cell types.

Discussion

Immunotherapy is the first-line treatment for advanced HCC.¹⁸ Immunotherapy is associated with improved OS compared with chemotherapy in advanced HCC. The combined efficacy and safety of immunotherapy and antiangiogenic therapy for unresectable HCC have been confirmed in multiple studies.¹⁹ IFNα and anti-PD1 cotreatment as a novel combination strategy for HCC has also been suggested.²⁰ A recent study demonstrated that perioperative camrelizumab plus apatinib displays promising efficacy and manageable toxicity in patients with resectable HCC.²¹ However, the benefits of immunotherapy are limited in certain patients. Therefore, identifying an effective biomarker is essential for immunotherapy.

Metabolic reprogramming as a hallmark of cancer can limit immune responses and present barriers to cancer therapy, and growing evidence suggests promising indications for targeting metabolism to improve the antitumor microenvironment.²² Tumor-associated macrophages promote tumor progression and resistance to immunotherapy by affecting the metabolic profile of the TME.²³ Reshaping lipid metabolism could inhibit effector T cell senescence induced by tumor cells and Tregs, and enhance tumor immunotherapy.²⁴ Dysregulated metabolism in HCC tissues is manifested mainly by decreased gene expression.²⁵ Lack of key regulators of metabolism can satisfy the demand for uncontrolled proliferation of tumor cells, and excess metabolites also participate in the creation of an immunosuppressive microenvironment.^25,26 Therefore, metabolic regulators may serve as ideal biomarkers for immunotherapy.

In the present study, we found that ACAA2 is downregulated during HCC development. Decreased ACAA2 levels are accompanied by an increased immunosuppressive tumor microenvironment. Knockdown of ACAA2 in tumor cells promotes progression of HCC malignancy. The lack of ACAA2 elevated the expression of CXCL1 via the NF-κB signaling pathway. CXCL1 promotes tumor progression and induces immune suppression by activating different signaling pathways.^27–29 GSEA revealed that ACAA2 and CXCL1 were inversely involved in several metabolic pathways. Our results demonstrate that, in addition to acquiring nutrients and increasing metabolites, metabolic reprogramming of cancer cells could modulate the TME in a paracrine manner. However, the key mechanism for ACAA2 regulating NF-κB signaling pathway is still unclear. ACAA2 is involved in lipid metabolism. Lipid metabolism not only provides the energetic needs of the cells but also provides the raw material for signaling molecular for many oncogenic signaling pathways.³⁰ Downregulated ACAA2 could impede the lipid metabolism and lead to the accumulation of intermediate metabolites, which might mediate decreased ACAA2 inducing the activation of NF-κB signaling pathway.

ACAA2 is also expressed in immune cells. The expression of ACAA2 in immune cells was negatively correlated with the composition of CD8⁺ T cells with PDCD1 positive in the tumor microenvironment. This effect may be mediated by CCL5–CCRs and the HLA-E–KLRCs LR network. A mouse experiment reported that CCL5 expression could restore immune surveillance in HCC.³¹ These results suggest that loss of ACAA2 facilitates the immune suppression microenvironment in a dual manner. Based on the TIDE scores, we found that the expression of ACAA2 was higher in the patient response to immunotherapy group than in the no response group.

Conclusion

In a conclusion, downregulated ACAA2 promotes the progression of hepatocellular carcinoma and participates in the formation of immunosuppressive microenvironment. We thought that ACAA2 could be an ideal biomarker for prognosis and immunotherapy of HCC.

Institutional Review Board Statement

All experiments were approved by the Ethics Committee of the first affiliated hospital of Harbin Medical University. The welfare of the laboratory animals was directed by the following guidelines: Regulations for the management of laboratory animals by State Science and Technology Commission of China, Guidance on the good treatment of laboratory animals by Ministry of Science and Technology of China, and Regulations of Heilongjiang Province on the management of laboratory animals by the Government of Heilongjiang Province.

Author Contributions

All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

Disclosure

The authors declare that they have no competing interests to disclose.

Acknowledgments

We thank for the technical supports of Shanghai Biotechnology Corporation.

Data Sharing Statement

The datasets presented in this study can be found in online datasets. The names of the datasets can be found in the article material. The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.

Additional information

Funding

This study was supported by the National Natural Science Foundation of China (81803001), Heilongjiang Province Youth Innovating Talents Training Plan Fund for Regular Undergraduate Colleges (UNPYSCT-2017066), Heilongjiang Province Postdoctoral Science Foundation (LBH-Q21034, 21042220024), Scientific Foundation of the First Affiliated Hospital of Harbin Medical University (HYD2020YQ0011, HYD2020JQ0007), Haiyan Research Fund of Harbin Medical University Cancer Hospital (JJMS2022-01).