Identification of Key Genes in the HBV-Related HCC Immune Microenvironment Using Integrated Bioinformatics Analysis

Ding, Wei; Zhang, Zheng; Ye, Nianyuan; Zhuang, Ling; Yuan, Zhiping; Xue, Wenbo; Tan, Yulin; Xu, Xuezhong

doi:https://doi.org/10.1155/2022/2797033

Journal of Oncology

On this page

Abstract Introduction Methods and Materials Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Novel Roles for Genomic Markers in Potential Targeted Therapy for Digestive Cancers

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 2797033 | https://doi.org/10.1155/2022/2797033

Identification of Key Genes in the HBV-Related HCC Immune Microenvironment Using Integrated Bioinformatics Analysis

Wei Ding,^1,2,3Zheng Zhang,^1,2Nianyuan Ye,^1,2Ling Zhuang,^1,2Zhiping Yuan,⁴Wenbo Xue,^1,2Yulin Tan ,^1,2and Xuezhong Xu^1,2

Academic Editor: Yanqing Liu

Received08 Sept 2022

Accepted28 Sept 2022

Published15 Oct 2022

Abstract

Purpose. Hepatocellular carcinoma (HCC) has poor prognosis and high mortality among gastrointestinal tumors because of its insidious onset and strong invasiveness. However, there was little understanding of their pathogenesis. The purpose of this study was to use bioinformatics analysis to identify genes associated with the immune microenvironment in HBV-related HCC and to develop new therapeutic targets to prevent and treat cancer. Methods. RNA-seq data of HBV-related HCC cases were downloaded from TCGA-LIHC database. ESTIMATE and Deseq2 algorithms were used to screen out differentially expressed genes (DEGs). WGCNA was used to construct gene coexpression networks. In key modules, functional enrichment analysis was performed. Protein-protein interaction (PPI) was used to screen hub genes, and survival analysis was conducted to assess their prognostic significance. Following, we search for key genes differentially expressed between cancerous and paracancerous tissues in GSE136247 and GSE121248 datasets. Reveal the potential links between key genes in immune infiltration by using TIMER. Finally, in TCGA-LIHC database, integration of key genes with clinical data were used to further validate their correlation with prognosis. Results. In the cohort of HBV-related HCC patients, immune/stromal/ESTIMATE scores were not significantly associated with patient prognosis. After bioinformatics analysis, screening out five key genes was significantly related to the prognosis of HBV-related HCC. Downregulation of SLAMF1 and TRAF3IP3 suggested poor prognosis and was related to a variety of immune cell infiltration. Furthermore, compared with adjacent nontumor tissues, TRAF3IP3 and SLAMF1 were highly expressed in tumor tissues and were linked to tumor recurrences. Conclusion. In conclusion, SLAMF1 and TRAF3IP3 were identified with higher expression in tumor tissues and associated with tumor recurrence. It will be a new research direction of tumor progress and treatment.

1. Introduction

About 90% of the pathological types of liver cancer were hepatocellular carcinoma (HCC) in clinical. Both morbidity and mortality are far higher than other types of tumors [1]. Gender differences (predominant in males) and geographic differences (mainly East Asia) influence the incidence of HCC [2]. The main causative factors associated with HCC are virus (chronic hepatitis B and C), metabolism (diabetes and nonalcoholic fatty liver disease), toxicity (alcohol and aflatoxins), and immune system-related diseases [3]. Affected by these factors, the morbidity of HCC is rising continuously year by year. According to statistics data, more than half of HCC patients in the world are infected with hepatitis B virus (HBV). It is the main risk factor for human [4]. Mortality associated with HCC is also increasing. Recent studies have shown that there were 85% of patients with cirrhosis infected with HBV, and less than 20% of people survived more than five years [5]. Reassuringly, the incidence of HCC was significantly reduced in the middle-aged population aged 30-59, largely due to the global hepatitis B virus vaccination program [6]. Liver transplantation and surgical resection are treatment options for HCC in early-stage HCC and when the tumor size is <5 cm [7]. However, because the disease is mostly asymptomatic in its early stages, most symptomatic patients are diagnosed at an advanced stage. Currently, patients with advanced HCC are mostly treated with radiofrequency ablation (RFA), transhepatic arterial chemotherapy and embolization (TACE), tyrosine kinase inhibitor (TKI), and immunotherapy, but with the emergence of drug resistance and disease recurrence, these modalities do not significantly prolong lifespan [8]. As research progresses in depth, new and diverse avenues for the treatment of HCC are being discovered.

The current study suggests that HCC is caused by HBV-induced DNA damage that triggers hepatocyte regeneration and chronic inflammation in the liver [9, 10]. The nucleocapsid of HBV-infected hepatocytes allows the virus to replicate stealthily without being recognized by type I IFN [11]. It is now widely believed that the immune pathogenic mechanism of HCC is mainly that HBV, as a noncytopathic virus, promotes the disorder of the liver immune system and causes liver damage through abnormal immune attack. It is increasingly believed that immune pathogenesis significantly influenced the development of HBV-related HCC [12]. Although HBV was thought to contribute to HCC, there was still no clear understanding of the mechanism.

In this study, RNA-seq data and clinical feature information of HBV-related HCC patients were accessed by TCGA-LIHC. The prognosis-related DEGs and modules were screened by Sangerbox and WGCNA. In addition, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched for these DEGs and module genes and construct the PPI network to search the prognosis-related hub genes. Through the gene expression information of GSE136247 and GSE121248 in the GEO database, the possible key genes (SLAMF1 and TRAF3IP3) were finally locked. Through the TIMER database, the association between SLAMF1 and TRAF3IP3 and immune cell infiltration was analyzed. Finally, by analyzing the clinical characteristics of HBV-related HCC patients, it was confirmed that SLAMF1 and TRAF3IP3 were negatively correlated with the recurrence of patients.

2. Methods and Materials

2.1. Data Acquisition

We gathered gene expression RNA-seq and accompanying medical data of HBV-related HCC patients from TCGA-LIHC database (https://portal.gdc.cancer.gov/) [13]. HBV-infected and noninfected patients were differentiated based on the patient’s past infection history. Download gene expression data from GSE136247 and GSE121248 datasets from GEO database (https://www.ncbi.nlm.nih.gov/geo/) [14]. The GSE136247 dataset contained 39 HCC tissues (25 with HBV infection) and 30 noncancerous normal tissues (19 with HBV infection) [15]. The GSE121248 dataset contained cancer and normal tissues from 37 HCC cases, and these patients had a history of HBV infection [16].

2.2. ESTIMATE Scores

The amount of tumor cells, immune cells, and stromal cells was determined using ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data) based on the transcriptional profile of cancer samples. According to the stromal signature (stromal signature gene) and immune signature (immune cell signature gene), the stromal score and immune score were calculated by ssGSEA, respectively. Finally, the two scores were combined to generate an ESTIMATE score, which was used to analyze tumor purity. From the official website (https://bioinformatics.mdanderson.org/estimate/), matrix, immune, and ESTIMATE scores were downloaded for each sample in TCGA-LIHC cohort, and non-HBV-infected samples were excluded. In addition, this score was compared with tumor patient survival in a correlation analysis.

2.3. Acquisition of DEGs Based on Immune and Matrix Scores

All HBV-related HCC patients were divided into two groups (high vs. low) with positive and negative values. Data analysis was performed on Sangerbox [17] using the “Deseq2” package. The filter range for DEGs were determined to be , .

2.4. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Function Analysis

Analyze the biological functions of DEGs by using the GO enrichment analysis (including BP, CC, and MF) and KEGG pathway enrichment analysis from DAVID online website tools (database annotation, visualization, and comprehensive discovery, https://david.ncifcrf.gov/tools.jsp) [18].

2.5. Weighted Correlation Network Analysis (WGCNA)

WGCNA is an analytical method for analyzing gene expression patterns of multiple samples, which can cluster the similar expression gene and investigate the association between specific traits and phenotypes in modules. It will help us find relevant biomarker genes and therapeutic targets. The “WGCNA” package was used to build the DEGs coexpression network on Sangerbox to identify the modules related to prognosis.

2.6. Construction of PPI Network and Filtration of Hub Gene

The protein information and PPI network information of key modules were analyzed using String database (https://cn.string-db.org/) [19]. MCODE was a plugin for constructing functional modules of gene (protein) network clustering in Cytoscape 3.8.0. According to the analysis results, the hub gene can be determined [20].

2.7. Survival Analysis

Plot the Kaplan-Meier survival curves of these hub genes, and screen out the hub genes significantly correlated with overall survival () by log-rank test.

2.8. Verification of the Expression of Hub Genes

Compare the previous hub genes with the expression data in GSE136247 and GSE121248 datasets to obtain the final key genes. We used “TIMER” in Sangerbox to analyze the correlation of key genes with 6 tumor-infiltrating immune cells in HBV-related HCC tissues.

2.9. Clinical Features in Patients with HBC-Related HCC

For comparison, the basic information and medical data of HBV-related HCC patients were separated into high and low groups based on the expression of key genes.

2.10. Statistical Analysis

Analysis in the present study were conducted using the R package on Sangerbox and GraphPad prism 8.0.2. We used log-rank tests and chi-square tests for data analysis. A statistically significant difference was considered to be less than 0.05. The whole process of bioinformatics analysis was shown in Figure 1.

3. Results

3.1. Scores of the Immune System and Stroma Correlated with Overall Survival

Based on TCGA database, the statistical data of 104 HBV-related HCC patients were gained. Patients ranged from 23 to 83 years of age. 85 (81.7%) were male, and 19 (19.3%) were female. For each sample, ESTIMATE scores were calculated based on matrix, immune, and ESTIMATE scores. Stromal scores ranged from -1731.43 to 261.96, immune scores ranged from -964.97 to 2311.6, and ESTIMATE scores ranged from -2488.91 to 2306.2. In order to probe the possibility of the connection between immune/stromal/ESTIMATE scores and patient survival, we categorize HBV-related HCC patients into low and high groups on the basis of 0 scores. There were no positive results between the two groups (Figures S1A, S1B, and S1C).

3.2. Identification of DEGs in HBV-Related HCC

For expounding the connection between gene expression profiles and immune status, we used “DESeq2” package to identify. Genes were significantly differential expression among the three groups of scores. and were as screening criteria. As shown in Figure 2(a), 571 downregulated genes and 1,845 upregulated genes were detected in the immune score group; in the stromal score group, 1,457 downregulated genes and 1,014 upregulated genes were detected; in the ESTIMATE score group, 1,052 were detected downregulated genes and 1,584 upregulated genes. According to the heat map, there were significant differences between the three groups in the differential genes (Figure 2(b)). Through further data screening, the differential express gene in all three groups were obtained, including 111 upregulated genes and 322 downregulated genes as shown in Figures 2(c) and 2(d).

(a)

(b)

(c)

(d)

3.3. Functional Enrichment Analysis

DAVID website was used for GO and KEGG analyses. As a result of the enrichment analysis, cellular components (CC), molecular functions (MF), and biological processes (BP) were enriched by GO enrichment analysis (Figure S2). For BP, DEGs were mainly enriched in external encapsulating structure organization, biological adhesion, and collagen fibril organization. For CC, DEGs were mainly enriched in collagen containing extracellular matrix, external encapsulating structure, and T cell receptor complex. For MF, DEGs were mainly enriched in extracellular matrix structural constituents, glycosaminoglycan binding, and heparin binding. For KEGG, DEGs were mainly enriched in the regulation of hematopoietic cell lineage, cytokine-cytokine receptor interaction, and viral protein interaction with cytokine-cytokine receptor (Figure S2).

3.4. Weighted Correlation Network Analysis

The role network of DEGs was constructed by WGCNA analysis. With the network’s soft threshold set at 16, coexpression networks resembled scale-free networks most closely (Figures 3(a)–3(c)). According to different functions, DEGs can be divided into 7 modules. Turquoise was the module with the highest significant difference over survival (OS time) (Figure 3(d)). The module contained a total of 50 genes.

(a)

(b)

(c)

(d)

3.5. Functional Enrichment Analysis of Genes in Turquoise

These genes were mostly associated with T cell activation, lymphocyte activation, and leukocyte differentiation in BP, based on GO enrichment analysis (Figure 4(a)). For CC, these genes were mainly enriched in immunological synapse, external side of plasma membrane, and side of membrane (Figure 4(b)). For MF, these genes were mainly enriched in cytokine receptor activity, immune receptor activity, and C-C chemokine binding (Figure 4(c)). For KEGG, these genes were mainly involved in the regulation of T cell receptor signaling pathway, cytokine-cytokine receptor interaction, and primary immunodeficiency (Figure 4(d)).

(a)

(b)

(c)

(d)

3.6. Filtration of PPI Network and Identification of Prognostic-Related Genes

Through the String database, a PPI network was constructed using 50 genes (Figure 5(a)). Then, use Cytoscape 3.8.0 to further optimize the obtained PPI network, and use the MCODE plugin to draw important subnetworks (Figures 5(b) and 5(c)). There were 16 central genes (CD53, TAGAP, IKZF1, CARD11, WDFY4, PTPRC, PTPN22, CYTIP, TRAF3IP3, CCR7, ITK, IL7R, CD40LG, SLAMF1, CD5, and SPN) in the protein interaction network.

(a)

(b)

(c)

3.7. Survival Analysis in Blue Module

Identifying genes associated with overall survival in patients with HBV-related HCC was the purpose of this research. We constructed Kaplan-Meier survival curves of these genes using the prognostic information in TCGA-LIHC. Among them, the expression levels of CCR7, CD5, SLAMF1, SPN, and TRAF3IP3 were significantly associated with the prognosis of patients (Figure 6 and Figure S3).

3.8. Validation of the Analysis in the GEO Database

In addition, we determined the use of GSE136247 and GSE121248 to explore the expression of these genes in cancerous and paracancerous tissues. As shown in Figure 7(a), in GSE136247, CCR7 was expressed significantly upregulated in cancer tissues (whether or not infected with HBV) relative to adjacent tissues, whereas TRAF3IP3 was completely the opposite. In HBV-related HCC patients, SLAMF1 and SPN in cancer tissues were significantly decreased, but this phenomenon was not observed in patients without HBV infection. In addition, the expression of CD5 was significantly decreased in cancer tissues relative to adjacent tissues of HBV-uninfected patients, while in HBV-related HCC patients, the two groups did not differ significantly. As shown in Figure 7(b), in GSE121248, the expressions of SLAMF1 and TRAF3IP3 were significantly decreased in HBV-related HCC tissues relative to paracancerous tissues, while the expressions of CCR7, CD5, and SPN were not significantly different. Based on the above results, we defined SLAMF1 and TRAF3IP3 as key genes for follow-up studies.

(a)

(b)

3.9. Connection between Key Genes and Immune Infiltration

In the present study, we explored possible associations between key gene expression and the infiltration of immune cells using TIMER. A positive correlation was found between SLAMF1 and TRAF3IP3, but not between neutrophils and macrophages, with the infiltration of B cells, CD4+ T cells, CD8+ T cells, and dendritic cells (Figure 8). In light of this, key genes may play an important role in regulating immune cells.

(a)

(b)

3.10. Relationship between Key Genes and Clinical Features

Based on TCGA-LIHC database, we examined the relationship between SLAMF1 and TRAF3IP3 and HBV-related HCC clinical characteristics. The results showed that SLAMF1 and TRAF3IP3 were inversely associated with tumor recurrence, regardless of gender, age, tumor stage (T), lymph node stage (N), and metastasis stage (M) (Table 1).

4. Discussion

HCC has a poor prognosis, killing more than 800,000 people worldwide each year [21]. In European and American countries, the prevalence of nonalcoholic fatty liver disease is rising rapidly every year, leading to the subsequent development of HCC and HCC-related death; while in developing countries in Asia, hepatitis and cirrhosis caused by viral infection are the main causes of HCC [22]. As most HCC cases are secondary to hepatitis (hepatitis B, hepatitis C, or alcoholic and nonalcoholic liver disease) or cirrhosis, HCC is now gradually considered to be the inflammatory cancer induced by chronic liver injury [23–26]. Patients with advanced HCC lack access to surgery and rely mainly on immunization or chemotherapy, example for sorafenib, a kind of the tyrosine kinase inhibitor. In recent years, several treatment options (lenvatinib, regorafenib, cabozantinib, and ramucirumab) have emerged for the various aspects’ treatment of advanced HCC [27].

Cellular components of the HCC immune microenvironment (tumor cells, immune cells, stromal cells, endothelial cells, and cancer-associated fibroblasts) are critical for the response to immunotherapy [28]. Through the portal vein, antigen-rich blood from the gut is constantly exposed to the liver, which acts as a central immune organ. In order to reduce inflammatory stimulation and tissue damage from the blood and the liver, establish an immune-tolerant microenvironment which has a strong resistance to hit and self-cleaning ability. The homeostasis of this immune microenvironment is also disrupted when hepatitis or cirrhosis or even HCC develops [29]. The TME in HCC is the hallmark of tumor, which has an important influence on tumor growth, invasion, and drug resistance [30, 31].

First, we screened out HCC patients with HBV infection from TCGA-LIHC. In the immune microenvironment of HBV-related HCC, DEGs were identified based on immunity, stroma, and ESTIMATE scores. Activation of CD4+ T cells, CD8+ T cells, NK cells, NKT cells, monocytes/macrophages, and HSCs occurs in chronic hepatitis caused by HBV. Hepatitis is further aggravated, and HCC is further encouraged by the simultaneous production of TNF-α, IFN-γ, IL-12, IL-4, and IL-13 [32, 33]. Additionally, several immunosuppressive cells, including Treg, Breg, MDSC, and Kupffer cells, inhibit immune cell activity by producing cytokines such as TGF-β and IL-10 and inducing key factors in CD8+ T and NK cell depletion, leading to immune escape of HBV and HCC tumor cell [34–37].

Additionally, we analyzed DEG enrichment. These DEGs have various biological properties and participate in various signaling pathways, such as external encapsulating structure organization, biological adhesion, T cell receptor complex, cytokine-cytokine receptor interaction, and viral protein interaction with cytokine-cytokine receptor. All of these confirmed that their involvement was in the regulation of the immune microenvironment in HBV-related HCC [38, 39]. Coexpression networks were constructed, with turquoise modules identified as key modules by WGCNA. Their main functions are to activate T cells, activate lymphocytes, and differentiate leukocytes; they are mainly located at the immune synapse, the outer and membrane sides of the plasma membrane; they mainly regulate cytokine receptor activity, immune receptor activity, and C-C chemokine binding. They are also involved in the regulation of T cell receptor signaling pathway cytokine-cytokine receptor interactions and primary immunodeficiency. A close correlation can be found between the immune regulation of HBV-related HCC and the genes of this module.

Through the PPI construction and prognostic-related genes analysis of this module, we identified five HBV-related hub genes for the prognosis of HCC patients, namely, CCR7, CD5, SLAMF1, SPN, and TRAF3IP3. On the basis of GSE136247 and GSE121248, the expression of each gene in cancer tissue and normal tissue was verified, and two key genes were finally obtained, namely, SLAMF1 and TRAF3IP3.

The SLAMF1/CD150 receptor is a member of the cell surface receptor signaling lymphocyte activation molecule (SLAM) family and is considered a marker of activated T cells, B cells, monocytes, and DCs [40, 41]. SLAMF1 is actively involved in the regulation of different types of immune responses as well as keeping the tissue microenvironment [42]. Recent studies have demonstrated that the expression level of SLAMF1 is significantly increased in liver tissue of NASH compared with non-NASH controls and that the level of SLAMF1 was dramatically related to the seriousness of the NASH phenotype. This study was the first to identify the role of SLAMF1 in the mediating of hepatocyte death in NASH and as a measure of NASH in humans [43]. In another study, the concentration of SLAMF1 has a profound effect on the formation of cirrhosis in the plasma. But no significant difference was found between HCC and cirrhosis [44]. TRAF3IP3 (TRAF3-interacting protein 3) was identified as a TRAF3-interacting protein in original [45]. Recent studies have shown that TRAF3IP3 is involved in B and T cell development and for maintaining the functional stability of regulatory T cells [46, 47]. TRAF3IP3 has been shown to function as an oncogene in melanoma and glioma [48, 49].

Our study showed that SLAMF1 and TRAF3IP3 were lowly expressed in HBV-related HCC and positively related with the infiltration of B cells, CD4+ T cells, CD8+ T cells, and dendritic cells, but not neutrophils and macrophages. Taken together, SLAMF1 and TRAF3IP3 may contribute to the pathogenesis of HBV-related HCC. Through their effect on the immune-suppressive microenvironment, furthermore, we found that SLAMF1 and TRAF3IP3 were also associated with the recurrence of HBV-related HCC.

5. Conclusion

We used bioinformatics to comprehensively analyze the expression of immune microenvironment-related genes in HBV-associated HCC patients in TCGA. Further study of the screened DEGs yielded two genes related with prognosis. We explained that SLAMF1 and TRAF3IP3 were low-expressed in HBV-associated HCC tissues and were correlated with tumor recurrence. Our findings had clear implications for SLAMF1 and TRAF3IP3 as biomarkers for predicting the prognosis of HBV-related HCC patients and provide new research directions and diagnosis and treatment options for HBV-related HCC. However, follow-up clinical studies are required to confirm these opinions.

Data Availability

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

DW and XWB conceived and designed the project; DW, ZL, and YZP acquired the data; DW, TYL, and XXZ analyzed and interpreted the data; DW, ZZ, and YNY wrote the paper. All authors contributed to the study and approved the submitted version. Wei Ding, Zheng Zhang, Nianyuan Ye, and Zhiping Yuan contributed equally to this work.

Acknowledgments

This work was supported by the Changzhou Sci&Tech Program (CJ20210013 and CJ20220008), Young Talent Development Plan of Changzhou Health Commission (CZQM2020118 and CZQM2021028), the Development Foundation of Affiliated Hospital of Xuzhou Medical University (XYFY2020016), Medical Research Project of Jiangsu Health Commission (no. Z2019027), and Changzhou High-Level Medical Talents Training Project (2022CZBJ105).

Supplementary Materials

Figure S1: the relationship between immune status and overall survival in HBV-related HCC. (A) Kaplan-Meier curve shows the overall survival of the high and low immune score groups. (B) Kaplan-Meier curve shows the overall survival of the high and low stromal score groups. (C) Kaplan-Meier curve shows the overall survival of the high and low ESTIMATE score groups. Figure S2: chord diagram demonstrates GO and KEGG analysis of DEGs. Biological processes (BP), cellular components (CC), molecular functions (MF), and KEGG pathways. Figure S3: survival analysis. The relations between the expression levels of CD53, TAGAP, IKZF1, CARD11, WDFY4, PTPRC, PTPN22, CYTIP, ITK, IL7R, and CD40LG and OS in HBV-related HCC. (Supplementary Materials)

References

J. Ferlay, I. Soerjomataram, R. Dikshit et al., “Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012,” International Journal of Cancer, vol. 136, no. 5, pp. E359–E386, 2015.
View at: Publisher Site | Google Scholar
E. Chakraborty and D. Sarkar, “Emerging therapies for hepatocellular carcinoma (HCC),” Cancers (Basel), vol. 14, no. 11, p. 2798, 2022.
View at: Publisher Site | Google Scholar
L. Rinaldi, M. Guarino, A. Perrella et al., “Role of liver stiffness measurement in predicting HCC occurrence in direct-acting antivirals setting: a real-life experience,” Digestive Diseases and Sciences, vol. 64, no. 10, pp. 3013–3019, 2019.
View at: Publisher Site | Google Scholar
A. P. Venook, C. Papandreou, J. Furuse, and L. Ladrón de Guevara, “The incidence and epidemiology of hepatocellular carcinoma: a global and regional perspective,” The Oncologist, vol. 15, Suppl 4, pp. 5–13, 2010.
View at: Publisher Site | Google Scholar
N. E. Rich, A. C. Yopp, A. G. Singal, and C. C. Murphy, “Hepatocellular carcinoma incidence is decreasing among younger adults in the United States,” Clinical Gastroenterology and Hepatology, vol. 18, no. 1, pp. 242–248.e5, 2020.
View at: Publisher Site | Google Scholar
M. Yarchoan, P. Agarwal, A. Villanueva et al., “Recent developments and therapeutic strategies against hepatocellular carcinoma,” Cancer Research, vol. 79, no. 17, pp. 4326–4330, 2019.
View at: Publisher Site | Google Scholar
S. Medavaram and Y. Zhang, “Emerging therapies in advanced hepatocellular carcinoma,” Experimental Hematology & Oncology, vol. 7, no. 1, p. 17, 2018.
View at: Publisher Site | Google Scholar
J. Musa, J. Li, and T. G. Grünewald, “Hepatitis B virus large surface protein is priming for hepatocellular carcinoma development via induction of cytokinesis failure,” The Journal of Pathology, vol. 247, no. 1, pp. 6–8, 2019.
View at: Publisher Site | Google Scholar
M. Levrero and J. Zucman-Rossi, “Mechanisms of HBV-induced hepatocellular carcinoma,” Journal of Hepatology, vol. 64, no. 1, pp. S84–s101, 2016.
View at: Publisher Site | Google Scholar
B. Rehermann and R. Thimme, “Insights from antiviral therapy into immune responses to hepatitis B and C virus infection,” Gastroenterology, vol. 156, no. 2, pp. 369–383, 2019.
View at: Publisher Site | Google Scholar
Y. Chen and Z. Tian, “HBV-induced immune imbalance in the development of HCC,” Frontiers in Immunology, vol. 10, p. 2048, 2019.
View at: Publisher Site | Google Scholar
Z. Wang, M. A. Jensen, and J. C. Zenklusen, “A practical guide to The Cancer Genome Atlas (TCGA),” Methods in Molecular Biology, vol. 1418, pp. 111–141, 2016.
View at: Publisher Site | Google Scholar
E. Clough and T. Barrett, “The Gene Expression Omnibus database,” Methods in Molecular Biology, vol. 1418, pp. 93–110, 2016.
View at: Publisher Site | Google Scholar
J. P. Cerapio, A. Marchio, L. Cano et al., “Global DNA hypermethylation pattern and unique gene expression signature in liver cancer from patients with Indigenous American ancestry,” Oncotarget, vol. 12, no. 5, pp. 475–492, 2021.
View at: Publisher Site | Google Scholar
S. M. Wang, L. L. Ooi, and K. M. Hui, “Identification and validation of a novel gene signature associated with the recurrence of human hepatocellular carcinoma,” Clinical Cancer Research, vol. 13, no. 21, pp. 6275–6283, 2007.
View at: Publisher Site | Google Scholar
W. Shen, Z. Song, X. Zhong et al., “Sangerbox: a comprehensive, interaction-friendly clinical bioinformatics analysis platform,” iMeta, vol. 1, no. 3, article e36, 2022.
View at: Publisher Site | Google Scholar
B. Sherman, M. Hao, J. Qiu et al., “DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update),” Nucleic Acids Research, vol. 50, no. W1, pp. W216–W221, 2022.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. Franceschini, S. Wyder et al., “STRING v10: protein-protein interaction networks, integrated over the tree of life,” Nucleic Acids Research, vol. 43, no. D1, pp. D447–D452, 2015.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: a Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
R. L. Siegel, K. D. Miller, and A. Jemal, “Cancer statistics, 2020,” CA: a Cancer Journal for Clinicians, vol. 70, no. 1, pp. 7–30, 2020.
View at: Publisher Site | Google Scholar
C. Trépo, H. L. Chan, and A. Lok, “Hepatitis B virus infection,” Lancet, vol. 384, no. 9959, pp. 2053–2063, 2014.
View at: Publisher Site | Google Scholar
B. Gao and R. Bataller, “Alcoholic liver disease: pathogenesis and new therapeutic targets,” Gastroenterology, vol. 141, no. 5, pp. 1572–1585, 2011.
View at: Publisher Site | Google Scholar
D. Y. Zhang and S. L. Friedman, “Fibrosis-dependent mechanisms of hepatocarcinogenesis,” Hepatology, vol. 56, no. 2, pp. 769–775, 2012.
View at: Publisher Site | Google Scholar
Q. M. Anstee, H. L. Reeves, E. Kotsiliti, O. Govaere, and M. Heikenwalder, “From NASH to HCC: current concepts and future challenges,” Nature Reviews. Gastroenterology & Hepatology, vol. 16, no. 7, pp. 411–428, 2019.
View at: Publisher Site | Google Scholar
B. Ruf, B. Heinrich, and T. F. Greten, “Immunobiology and immunotherapy of HCC: spotlight on innate and innate-like immune cells,” Cellular & Molecular Immunology, vol. 18, no. 1, pp. 112–127, 2021.
View at: Publisher Site | Google Scholar
F. Heymann and F. Tacke, “Immunology in the liver--from homeostasis to disease,” Nature Reviews. Gastroenterology & Hepatology, vol. 13, no. 2, pp. 88–110, 2016.
View at: Publisher Site | Google Scholar
J. P. Böttcher, P. A. Knolle, and D. Stabenow, “Mechanisms balancing tolerance and immunity in the liver,” Digestive Diseases, vol. 29, no. 4, pp. 384–390, 2011.
View at: Publisher Site | Google Scholar
D. F. Quail and J. A. Joyce, “Microenvironmental regulation of tumor progression and metastasis,” Nature Medicine, vol. 19, no. 11, pp. 1423–1437, 2013.
View at: Publisher Site | Google Scholar
M. Binnewies, E. W. Roberts, K. Kersten et al., “Understanding the tumor immune microenvironment (TIME) for effective therapy,” Nature Medicine, vol. 24, no. 5, pp. 541–550, 2018.
View at: Publisher Site | Google Scholar
J. Wang, W. Zhao, L. Cheng et al., “CD137-mediated pathogenesis from chronic hepatitis to hepatocellular carcinoma in hepatitis B virus-transgenic mice,” Journal of Immunology, vol. 185, no. 12, pp. 7654–7662, 2010.
View at: Publisher Site | Google Scholar
X. Hou, X. Hao, M. Zheng et al., “CD205-TLR9-IL-12 axis contributes to CpG-induced oversensitive liver injury in HBsAg transgenic mice by promoting the interaction of NKT cells with Kupffer cells,” Cellular & Molecular Immunology, vol. 14, no. 8, pp. 675–684, 2017.
View at: Publisher Site | Google Scholar
Q. Wang, W. Pan, Y. Liu et al., “Hepatitis B virus-specific CD8+ T cells maintain functional exhaustion after antigen reexposure in an acute activation immune environment,” Frontiers in Immunology, vol. 9, p. 219, 2018.
View at: Publisher Site | Google Scholar
N. Trehanpati and A. K. Vyas, “Immune regulation by T regulatory cells in hepatitis B virus-related inflammation and cancer,” Scandinavian Journal of Immunology, vol. 85, no. 3, pp. 175–181, 2017.
View at: Publisher Site | Google Scholar
Q. F. Zhang, W. W. Yin, Y. Xia et al., “Liver-infiltrating CD11b⁻CD27⁻ NK subsets account for NK-cell dysfunction in patients with hepatocellular carcinoma and are associated with tumor progression,” Cellular & Molecular Immunology, vol. 14, no. 10, pp. 819–829, 2017.
View at: Publisher Site | Google Scholar
M. Li, R. Sun, L. Xu et al., “Kupffer cells support hepatitis B virus-mediated CD8+ T cell exhaustion via hepatitis B core antigen-TLR2 interactions in mice,” Journal of Immunology, vol. 195, no. 7, pp. 3100–3109, 2015.
View at: Publisher Site | Google Scholar
P. Qin, M. Zhang, X. Liu, and Z. Dong, “Immunogenomic landscape analysis of prognostic immune-related genes in hepatocellular carcinoma,” Journal of Healthcare Engineering, vol. 2021, Article ID 3761858, 13 pages, 2021.
View at: Publisher Site | Google Scholar
S. Xiang, J. Li, J. Shen et al., “Identification of prognostic genes in the tumor microenvironment of hepatocellular carcinoma,” Frontiers in Immunology, vol. 12, article 653836, 2021.
View at: Publisher Site | Google Scholar
V. G. Pinchouk, S. P. Sidorenko, D. F. Gluzman, E. P. Vetrova, A. G. Berdova, and L. N. Schlapatskaya, “Monoclonal antibodies IPO-3 and IPO-10 against human B cell differentiation antigens,” Anticancer Research, vol. 8, no. 6, pp. 1377–1380, 1988.
View at: Google Scholar
S. P. Sidorenko, E. P. Vetrova, O. V. Yurchenko, A. G. Berdova, L. N. Shlapatskaya, and D. F. Gluzman, “Monoclonal antibodies of IPO series against B cell differentiation antigens in leukemia and lymphoma immunophenotyping,” Neoplasma, vol. 39, no. 1, pp. 3–9, 1992.
View at: Google Scholar
I. Gordiienko, L. Shlapatska, L. Kovalevska, and S. P. Sidorenko, “SLAMF1/CD150 in hematologic malignancies: silent marker or active player?” Clinical Immunology, vol. 204, pp. 14–22, 2019.
View at: Publisher Site | Google Scholar
O. Gomez-Torres, S. Amatya, L. Kamberov et al., “SLAMF1 is expressed and secreted by hepatocytes and the liver in nonalcoholic fatty liver disease,” American Journal of Physiology. Gastrointestinal and Liver Physiology, vol. 323, no. 3, pp. G177–g187, 2022.
View at: Publisher Site | Google Scholar
I. Argirion, R. M. Pfeiffer, T. K. Lam et al., “Association between immunologic markers and cirrhosis in individuals with chronic hepatitis B,” Scientific Reports, vol. 11, no. 1, p. 21194, 2021.
View at: Publisher Site | Google Scholar
H. Dadgostar, S. E. Doyle, A. Shahangian, D. E. Garcia, and G. Cheng, “T3JAM, a novel protein that specifically interacts with TRAF3 and promotes the activation of JNK1,” FEBS Letters, vol. 553, no. 3, pp. 403–407, 2003.
View at: Publisher Site | Google Scholar
X. Yu, X. L. Teng, F. Wang et al., “Metabolic control of regulatory T cell stability and function by TRAF3IP3 at the lysosome,” The Journal of Experimental Medicine, vol. 215, no. 9, pp. 2463–2476, 2018.
View at: Publisher Site | Google Scholar
S. Peng, K. Wang, Y. Gu et al., “TRAF3IP3, a novel autophagy up-regulated gene, is involved in marginal zone B lymphocyte development and survival,” Clinical and Experimental Immunology, vol. 182, no. 1, pp. 57–68, 2015.
View at: Publisher Site | Google Scholar
G. Yang, S. Tang, J. Zhang, and L. Qin, “High TRAF3IP3 level predicts poor prognosis of patients with gliomas,” World Neurosurgery, vol. 148, pp. e436–e449, 2021.
View at: Publisher Site | Google Scholar
P. Nasarre, I. V. Bonilla, J. S. Metcalf, E. G. Hilliard, and N. Klauber-DeMore, “TRAF3-interacting protein 3, a new oncotarget, promotes tumor growth in melanoma,” Melanoma Research, vol. 28, no. 3, pp. 185–194, 2018.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Wei Ding et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

406

Downloads

522

Citations