- 1Department of Infectious Diseases and Liver Diseases, Ningbo Medical Center Lihuili Hospital, Affiliated Lihuili Hospital of Ningbo University, Ningbo, China
- 2Department of Infectious Diseases, Shaoxing People's Hospital, Shaoxing, China
Background: Immune-related genes (IRGs) are closely connected to the occurrence and development of tumors. Their influence on the prognosis of patients with HCC, however, remains unclear.
Methods: From the TCGA database, we integrated 365 liver cancer tissues and 50 normal tissues to identify differential immune genes related to prognosis. Multivariate COX analysis was used to establish a new prognostic index on account of IRGs, whereby risk score = (Expression level of HSPA4*0.022) + (Expression level of PSMD14*0.042) + (Expression level of RBP2*0.019) + (Expression level of MAPT*0.197) + (Expression level of TRAF3*0.146) + (Expression level of NDRG1*(0.006) + (Expression level of NRAS*0.027) + (Expression level of IL17D*0.075).
Results: The risk score was clearly correlated with an unfavorable survival rate and with clinical characteristics. By integrating the immune-related risk score model with clinical features, a nomogram was constructed to predict the survival rate of HCC patients (1-, 3- and 5-year AUC of 0.721, 0.747 and 0.781, respectively).
Conclusion: We have established a valuable prognostic risk score for HCC patients that may be a better predictor of survival than the present method. With the risk score's strong predictive value for immune cells and functions, it may provide clinical guidance for the diagnosis and prognosis of different immunophenotypes, and provide multiple therapeutic targets for the treatment of HCC patients based on subtype-specific immune molecules.
Introduction
Hepatocellular carcinoma (HCC) is the second prime cause of cancer-related deaths internationally (1). Approximately 782,500 new cases of HCC are diagnosed worldwide each year (2); poor prognosis remains the most common feature in HCC patients due to their late diagnosis (3). HCC has several known risk factors, including hepatitis viral infection, alcohol insobriety, autoimmune hepatitis, adiposis and grievous metabolic diseases (4). The occurrence of HCC is closely related to an inflammatory microenvironment caused by virus infection, obesity and other liver damage (5). Tumor-associated immune cells are the main component of the inflammatory microenvironment of liver cancer. These cells respond to environmental signals and participate in the development of tumors (6).
Current clinical treatment options for HCC are varied and include surgical resection, liver transplantation, targeted therapy and immunotherapy. However, the prognosis of patients with liver cancer remains far from satisfactory (7). Consequently, the confirmation of prognostic biomarkers and novel drug targets is essential in providing HCC patients with better prognostic information and more efficient personalized treatments (8). Immune-related genes (IRGs) play a key role in regulating the systemic immune response (9). Improving our knowledge of IRGs is critical to understanding the mechanism of cancer immunotherapy. Some key IRGs could also serve as biomarkers to predict the outcome of cancer patients (10).
The aims of this research are to evaluate the prognostic significance of IRGs in HCC, to construct an IRG-based prognostic score, and to investigate the possible mechanism by which these immune genes affect the development of liver cancer. Ideally, the predictive scoring model could help clinicians estimate the risk of HCC in clinical practice and provide important information for personalized immunotherapy of HCC.
Methods
TCGA data and IRG identification
Transcriptome RNA sequencing data and clinical information for HCC was downloaded from the TCGA database (https://cancergenome.nih.gov/). There were 374 liver cancer tissue samples and 50 normal liver tissue samples. Overall, 365 patients were included, while HCC patients with zero overall survival were excluded. IRGs were identified from the ImmPort database (https://www.immport.org/), with 2,483 genes included in this study.
Differentially expressed IRG
The R LIMMA software package (FDRFilter = 0.05, logFCFilter = 1 and Wilcox Test) was used to screen for differentially expressed genes between liver cancer and normal liver tissues. Data for these differentially expressed IRGs was then extracted. The R heatmap software package was used to display data as heat maps and volcano maps.
Construction and validation of IRG-based prognostic model
To estimate the relationship between IRG expression and the overall survival rate of patients, we applied the univariate Cox proportional hazard regression analysis. A predictive model for IRGs was constructed depending upon their expression and the contribution coefficient (β) from the multivariate Cox proportional hazard regression model. The risk score was calculated for each patient and used to classify them into two groups (low- and high-risk) based on the optimal cut-off value for gene expression. The ICGC database (243 HCC samples) was applied to validate the precision of the risk score model.
The relationship of Risk score with clinical information and immune cells
The clinical information available included gender, age, survival time, AJCC (American Joint Committee on Cancer) tumor grade and TNM (Tumor, Node and Metastasis) staging system. Tumor-infiltrating immune cells involved B cells, CD4 T cells, CD8 T cells, macrophages, neutrophils and dendritic cells. The TIMER online database analyzed and visualized the number of tumor-infiltrating immune cells (https://cistrome.shinyapps.io/timer/). We also used the Tumor Immunity Estimation Resource to analyze the correlation between risk score and 6 immune infiltrating cell types. The purity-corrected partial Spearman's correlation coefficient was used to evaluate the relationship between them. The flow chart for this research is shown in Figure 1.
Construction and verification of the nomogram
Employing clinical information from the TCGA-LIHC cohort, multivariate Cox regression analysis was used to identify significant clinical variables. The nomogram was constructed using the R “rms” package combining the clinicopathological characteristics and risk scores of the patients. The ROC curve was used to evaluate the nomogram's 1-, 3- and 5-year prognosis, and the calibration curve was constructed by the bootstrap method (1,000 cycles) displaying the deviation between the predicted values and the actual probability of occurrence.
Statistical analysis
Statistical analyses were performed using R software (version 3.6). Pearson correlation coefficient was used to evaluate correlations among genes in coexpression analysis. The patients were divided into high or low group according to the median expression of the gene in the TCGA and ICGC database, respectively. Kaplan–Meier curves were constructed and the log-rank test was performed using R package survival software. Univariate and multivariate Cox regression analyses were used to determine independent prognostic factors. ROC analysis with area under the curve (AUC) was used to evaluate the accuracy of the Risk score. The nomogram was developed by incorporating the clinicopathological features and risk scores of patients using the R “rms” package. A P-value of <0.05 indicated statistical significance.
Results
Identification of differentially expressed IRG
Data from the TCGA database was extracted from 365 cases of liver cancer and 50 cases of non-diseased liver. 7,667 differentially expressed genes (DEG) were verified between the two groups, including 7,273 up-regulated and 394 down-regulated genes (Figures 2A,B). We then cross-matched these DEGs against the list of IRGs to identify differentially expressed immune genes. Among the total DEG, we identified 335 differentially expressed IRGs, of which 267 were up-regulated and 62 were down-regulated (Figures 2C,D).
Construction of risk score model
The most important indicator of prognosis is patient survival rate. Accordingly, we investigated for associations between IRG expression and patient survival using univariate Cox proportional hazard regression analysis. The results showed that 31s IRG were remarkably related with the survival of HCC patients (Figure 3A). We also conducted enrichment analysis of the KEGG pathways for these 31 IRGs (Figure 3B). Cox regression analysis of these genes was applied to construct a risk score model. The calculation formula for the risk score contained eight IRGs, the coefficients and HR values of which are shown in Table 1. The calculation formula for the Risk score is as follows: Risk score = (Expression level of HSPA4*0.022) + (Expression level of PSMD14*0.042) + (Expression level of RBP2*0.019) + (Expression level of MAPT*0.197) + (Expression level of TRAF3*0.146) + (Expression level of NDRG1*(0.006) + (Expression level of NRAS*0.027) + (Expression level of IL17D*0.075).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819491/fsurg-09-819491-HTML/image_m/fsurg-09-819491-g003.jpg)
Figure 3. Analysis of survival-related differentially expressed IRGs. (A) Hazard ratio forest chart shows the predicted value of IRGs. (B) Analysis of KEGG pathway of survival-related IRGs.
Evaluation and validation of risk score performance
HCC patients were divided into a high-risk group (risk score ≥7.65) or low-risk group (risk score <7.65) according to the median in TCGA-LIHC. The difference in IRG expression profiles between the two groups is demonstrated in Figure 4A, and risk score ranking is demonstrated in Figure 4B. The mortality rate observed for the high-risk group was remarkably higher than for the low-risk group (Figure 4C), while the AUROC of the risk score for 1- and 3-year mortality were 0.765 and 0.665, respectively (Figure 4D). The risk score was validated using the ICGC database, which includes 243 LIHC samples, HCC patients were divided into a high-risk group (risk score ≥ 1.24) or low-risk group (risk score <1.24) according to the median. The difference in IRG expression profiles between the two groups is demonstrated in Figure 4E and risk score ranking is demonstrated in Figure 4F. The mortality rate observed for the high-risk group was remarkably higher than for the low-risk group (Figure 4G); the AUROC of the risk score for 1- and 3-year mortality were 0.642 and 0.674, respectively (Figure 4H).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819491/fsurg-09-819491-HTML/image_m/fsurg-09-819491-g004.jpg)
Figure 4. Prognostic evaluation and validation based on IRG risk score in TCGA-LIHC and ICGC-LIHC database. (A) Heatmap of gene expression profiles in high-risk and low-risk populations. (B) The grade and distribution of prognostic indicators for high-risk and low-risk populations. (C) Comparison of survival time between high-risk and low-risk groups. (D) Risk score for predicting the 1- and 3-year overall survival (OS) rates of HCC patients. (E) Heatmap of gene expression profiles in high-risk and low-risk populations. (F) The grade and distribution of prognostic indicators for high-risk and low-risk populations. (G) Comparison of survival time between high-risk and low-risk groups. (H) Risk score for predicting the 1- and 3-year OS rates of HCC patients. (A–D from the TCGA-LIHC database and E–H from the ICGC-LIHC database).
Construction and verification of the nomogram
With clinical information from the TCGA-LIHC cohort, multivariate Cox regression analysis was used to identify significant clinical variables (Figure 5A). A nomogram was constructed by integrating risk scores and clinicopathological features to predict 1-, 3-, and 5-year overall survival in HCC patients (Figure 5B). The AUCs were 0.721, 0.747 and 0.781 for predicting in 1-, 3- and 5-years OS, respectively (Figure 5C). The calibration curve shows that the predicted and actual results of our nomogram are significantly consistent (Figure 5D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819491/fsurg-09-819491-HTML/image_m/fsurg-09-819491-g005.jpg)
Figure 5. Construction of the nomogram and analysis of ssGSEA scores. (A) Results of the multivariate Cox regression analyses regarding OS in the TCGA cohort. (B) The nomogram for predicting 1-, 3- and 5-years overall survival in TCGA-LIHC cohort. (C) ROC analyses of the nomogram predicting of 1-, 3- and 5-years overall survival. (D) Calibration curve for validation of the nomogram. (E) Comparison of the ssGSEA scores between different risk groups in the TCGA cohort.
Relationship between IRG-based Risk score and clinical features
To further explore the correlation of risk scores with immune status, we quantified enriched fractional subsets and associated functions of different immune cells by ssGSEA, including 16 immune cells and 13 immune-related function analyses. Interestingly, there were significant differences in both immune cell subsets and cellular function between the high- and low-risk groups. The scores for type-II IFN response, mast-cells and natural killer cells were lower in the high-risk group, while those for the macrophages, Tregs and Th2-Cells showed an opposite trend (Figure 5E). Univariate and multivariate Cox analyses were used to evaluate the prognostic value of the risk score (Table 2). We further analyzed the correlation between the IRG-based risk score and clinical information. The risk score showed significant positive correlations with liver cancer patient grade, AJCC stage and TNM stage. The risk score-derived immune gene HSPA4 showed significant correlation with grade, IL17D with gender and grade, and TRAF3 with stage and TNM (Figure 6). We further analyzed the correlation between risk score and immune cell infiltration. The number of CD8 T cells, dendritic cells, macrophages and neutrophils showed significant positive correlation with risk score, whereas CD4 T cells and B cells did not show any significant correlation (Figure 7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819491/fsurg-09-819491-HTML/image_m/fsurg-09-819491-g006.jpg)
Figure 6. The correlation between risk score and clinical information. The risk score showed significant positive correlation with liver cancer patient grade (A), stage (B) and TNM stage (C). The risk score-derived immune gene TRAF3 showed significant correlation with stage and TNM (D,E), HSPA4 with grade (F), IL17D with gender and grade (G,H).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819491/fsurg-09-819491-HTML/image_m/fsurg-09-819491-g007.jpg)
Figure 7. Relationship between the risk score and infiltration abundance of six immune cells. We further analyzed the correlation between risk score and immune cell infiltration. The number of macrophages (A), neutrophils (B), dendritic cells (C) and CD8 T cells (D) showed significant positive correlation with risk score, whereas CD4 T cells (E) and B cells (F) did not show any significant correlation.
Transcription factors and regulatory networks
Transcription factors (TFs) directly control gene expression levels. A total of 318 TFs were present in the Cistrome Cancer database (http://cistrome.org/). Comparison of the two risk score groups (high and low) revealed 117 differentially expressed TFs. Cytoscape (V3.6.1) was used to assemble regulatory networks of the survival-related IRGs for these TFs (Figure 8).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819491/fsurg-09-819491-HTML/image_m/fsurg-09-819491-g008.jpg)
Figure 8. TF-mediated regulatory network. Heatmap (A) and volcano (B) plot of differentially expressed TFs. (C) Construct a regulatory network based on differentially expressed TFs and survival-related IRGs (green triangles indicate TFs and pink polygons indicate IRGs).
Further validation of the IRG-based Risk score
The expression of IRGs that make up the risk score showed a significant positive correlation with patient prognosis and related clinical information. We further evaluated whether the expression of these IRGs was different between normal individuals and liver cancer patients. The protein expression of these genes in HCC using human protein atlas database. However, expression of these immune genes was remarkably higher in the hepatic tissue of liver cancer patients compared to non-liver cancer individuals (Figure 9).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/819491/fsurg-09-819491-HTML/image_m/fsurg-09-819491-g009.jpg)
Figure 9. IRGs based on risk score were validated in liver tissue using human protein atlas database. CA: the hepatic tissue of liver cancer patients, Nor: the hepatic tissue of non-liver cancer patients.
Discussion
There is increasing evidence that expression of immune genes is closely associated with the occurrence and progression of tumors (11, 12). Consequently, immunotherapy is now commonly used against various cancer types (13, 14), including HCC (15). Providing patients with a suitable treatment plan has become the core element of treatment (16). Although immunotherapy has been used to treat HCC, nevertheless, not all patients benefit from it. Thus, it is essential to develop an IRG-based score that can accurately predict the response to cancer immunotherapy. This will allow the most appropriate patients to be selected for treatment. Currently, several IRGs have been used to assess whether immunotherapy may be suitable and to predict the outcome of patients (17).
In our study, eight IRGs (HSPA4, PSMD14, RBP2, MAPT, TRAF3, NDRG1, NRAS and IL17D) with prognostic value in HCC patients were identified by Cox regression analysis and LASSO regression modeling. These IRGs, which were used to construct a risk score to evaluate the prognosis of HCC patients, showed significant positive correlations with HCC patient grade, AJCC stage and TNM. In other words, a high risk score correlated with worse patient prognosis. PSMD14 expression was significantly up-regulated in HCC tissues and was significantly associated with vascular invasion, tumor number, tumor recurrence and poor overall survival in HCC patients (18). High levels of NDRG1 expression correlated with poor prognosis in HCC patients (19), while NRAS overexpression was also associated with poor survival in HCC patients and proliferation in vivo (20).
HSPA4 expression was positively correlated with the grade of liver cancer in our study. Previous studies have shown that HSPA4 was overexpressed in liver cancer and was involved in tumor migration, invasion, transformation and recurrence (21, 22). TRAF3 has been reported to promote the proliferation of lung adenocarcinoma (23), but has not been well studied in liver cancer. Recent work has shown significant improvements in obesity, insulin resistance, liver steatosis and inflammatory responses in liver-specific TRAF3 knockout mice (24). Conversely, transgenic mice that overexpress TRAF3 in liver cells showed increased obesity. We found here that TRAF3 expression correlated significantly with the stage and TNM of HCC patients, suggesting it may be an important factor in the development of HCC.
Inflammatory processes contribute to immune tolerance as well as to tumor progression and metastasis. The role of IL17 in carcinogenesis, tumor metastasis and resistance to chemotherapy in several solid cancer types (25), together with high expression of IL17 in colorectal cancers, indicates an association with poor prognosis (26). Our results showed that IL17 expression correlated with tumor grade in HCC patients and that its expression was also significantly higher in female patients compared to males. IL17-mediated inflammation shows gender differences in a variety of diseases, including allergic asthma and urinary tract infection (27); to our knowledge, however, this has not been reported in liver cancer patients. HSPA4 was positively associated with immune cell infiltration and immune checkpoints (PD-1 and CTLA-4) in HCC (28). We validated the value of HSPA4, TRAF3 and IL17 in the diagnosis and prediction of prognosis of HCC. These proteins may be involved not only in the occurrence and development of HCC, but in its immune regulation as well. Therefore, they can be a potential diagnostic and prognostic biomarker as well as a therapeutic target for HCC. Although direct associations between HCC and the RBP2 and MAPT genes have yet to be reported, we believe these potential associations warrant further experimental study.
Furthermore, the high-risk group had higher proportions of macrophages, regulatory T cells (Tregs) and Th2 cells, whereas the activity of type II IFN responses, Mast cells and NK cells were lower in the high-risk group, which was associated with impaired antitumor immunity. Tumor-infiltrating macrophages are one of the most abundant stromal cell types in the HCC tumor microenvironment, not only suppressing anti-tumor immunity by inducing extracellular matrix remodeling, angiogenesis, tumor metastasis, and therapy resistance, but also secreting various inflammatory mediators that promote tumor cell progression (29). Tregs impair cancer immune surveillance by creating an immunosuppressive environment that promotes tumor cell survival (30). Disturbed immune response caused by up-regulation of Th2 cytokines and down-regulation of Th1 cytokines is a key predictor of HCC metastasis (31). The activity of type II IFN response, together with reduction of Mast and NK cells, was significantly associated with tumor immune escape and poor immunotherapy efficacy (32, 33). Therefore, weakened antitumor immunity in high-risk patients may be the key factor for the poor prognosis.
The tumor microenvironment coexists and interacts with a variety of immune cells to maintain the growth of HCC (34). The IRGs screened in the present study play an important role in the immune response to HCC. We evaluated the relationship between risk score and different types of immune cells and found an association with macrophages and neutrophils, in particular. The risk score developed herein was found to be valuable for predicting the outcome of HCC patients. Eight clinically significant IRGs in HCC were screened using bioinformatics, helping to expand our knowledge of the molecular mechanisms relating to HCC immunology. One of the limitations of our study was the small number of patients included in the validation group. We suggest that more HCC patients should be investigated in prospective studies to further understand the underlying mechanisms. In conclusion, we have established a valuable prognostic risk score for HCC patients, which may serve as a more accurate predictor of survival than the present method. The risk score's improved predictions for immune cells and functions can provide clinical guidance for the diagnosis and prognosis of different immunophenotypes, and provide multiple therapeutic targets for the treatment of HCC patients based on subtype-specific immune molecules.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author/s.
Ethics statement
Our study is based on open source data, so there are no ethical issues and other conflicts of interest. The authors state that no human or animal experiments were adopted in this study.
Author contributions
QC contributed to the conception and design. JD, LD contributed to the provision of the study materials or patients. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Dong Z, Chen Y, Yang C, Zhang M, Chen A, Yang J, et al. STAT Gene family mRNA expression and prognostic value in hepatocellular carcinoma. Onco Targets Ther. (2019) 12:7175–91. doi: 10.2147/ott.S202122
2. Forner A, Llovet JM, Bruix J. Hepatocellular carcinoma. Lancet. (2012) 379(9822):1245–55. doi: 10.1016/s0140-6736(11)61347-0
3. Kobayashi T, Aikata H, Kobayashi T, Ohdan H, Arihiro K, Chayama K. Patients with early recurrence of hepatocellular carcinoma have poor prognosis. Hepatobiliary Pancreat Dis Int. (2017) 16(3):279–88. doi: 10.1016/s1499-3872(16)60181-9
4. Ren Z, Li A, Jiang J, Zhou L, Yu Z, Lu H, et al. Gut microbiome analysis as a tool towards targeted non-invasive biomarkers for early hepatocellular carcinoma. Gut. (2019) 68(6):1014–23. doi: 10.1136/gutjnl-2017-315084
5. Wu Q, Li J, Li Z, Sun S, Zhu S, Wang L, et al. Exosomes from the tumour-adipocyte interplay stimulate beige/brown differentiation and reprogram metabolism in stromal adipocytes to promote tumour progression. J Exp Clin Cancer Res. (2019) 38(1):223. doi: 10.1186/s13046-019-1210-3
6. Wang L, Wang FS. Clinical immunology and immunotherapy for hepatocellular carcinoma: current progress and challenges. Hepatol Int. (2019) 13(5):521–33. doi: 10.1007/s12072-019-09967-y
7. Rimassa L, Pressiani T, Merle P. Systemic treatment options in hepatocellular carcinoma. Liver Cancer. (2019) 8(6):427–46. doi: 10.1159/000499765
8. Couri T, Pillai A. Goals and targets for personalized therapy for HCC. Hepatol Int. (2019) 13(2):125–37. doi: 10.1007/s12072-018-9919-1
9. Chen W, Ou M, Tang D, Dai Y, Du W. Identification and validation of immune-related gene prognostic signature for hepatocellular carcinoma. J Immunol Res. (2020) 2020:5494858. doi: 10.1155/2020/5494858
10. Masucci GV, Cesano A, Hawtin R, Janetzki S, Zhang J, Kirsch I, et al. Validation of biomarkers to predict response to immunotherapy in cancer: volume I – pre-analytical and analytical validation. J Immunother Cancer. (2016) 4(1):76. doi: 10.1186/s40425-016-0178-1
11. Berraondo P, Minute L, Ajona D, Corrales L, Melero I, Pio R. Innate immune mediators in cancer: between defense and resistance. Immunol Rev. (2016) 274(1):290–306. doi: 10.1111/imr.12464
12. Elola MT, Ferragut F, Méndez-Huergo SP, Croci DO, Bracalente C, Rabinovich GA. Galectins: multitask signaling molecules linking fibroblast, endothelial and immune cell programs in the tumor microenvironment. Cell Immunol. (2018) 333:34–45. doi: 10.1016/j.cellimm.2018.03.008
13. Banerjee K, Kumar S, Ross KA, Gautam S, Poelaert B, Nasser MW, et al. Emerging trends in the immunotherapy of pancreatic cancer. Cancer Lett. (2018) 417:35–46. doi: 10.1016/j.canlet.2017.12.012
14. Sanmamed MF, Chen L. A paradigm shift in cancer immunotherapy: from enhancement to normalization. Cell. (2018) 175(2):313–26. doi: 10.1016/j.cell.2018.09.035
15. Iñarrairaegui M, Melero I, Sangro B. Immunotherapy of hepatocellular carcinoma: facts and hopes. Clin Cancer Res. (2018) 24(7):1518–24. doi: 10.1158/1078-0432.Ccr-17-0289
16. Finkelmeier F, Waidmann O, Trojan J. Nivolumab for the treatment of hepatocellular carcinoma. Expert Rev Anticancer Ther. (2018) 18(12):1169–75. doi: 10.1080/14737140.2018.1535315
17. Ge P, Wang W, Li L, Zhang G, Gao Z, Tang Z, et al. Profiles of immune cell infiltration and immune-related genes in the tumor microenvironment of colorectal cancer. Biomed Pharmacother. (2019) 118:109228. doi: 10.1016/j.biopha.2019.109228
18. Lv J, Zhang S, Wu H, Lu J, Lu Y, Wang F, et al. Deubiquitinase PSMD14 enhances hepatocellular carcinoma growth and metastasis by stabilizing GRB2. Cancer Lett. (2020) 469:22–34. doi: 10.1016/j.canlet.2019.10.025
19. Luo Q, Wang CQ, Yang LY, Gao XM, Sun HT, Zhang Y, et al. FOXQ1/NDRG1 Axis exacerbates hepatocellular carcinoma initiation via enhancing crosstalk between fibroblasts and tumor cells. Cancer Lett. (2018) 417:21–34. doi: 10.1016/j.canlet.2017.12.021
20. Dietrich P, Gaza A, Wormser L, Fritz V, Hellerbrand C, Bosserhoff AK. Neuroblastoma RAS viral oncogene homolog (NRAS) is a novel prognostic marker and contributes to sorafenib resistance in hepatocellular carcinoma. Neoplasia. (2019) 21(3):257–68. doi: 10.1016/j.neo.2018.11.011
21. Yang MH, Chang SY, Chiou SH, Liu CJ, Chi CW, Chen PM, et al. Overexpression of NBS1 induces epithelial-mesenchymal transition and co-expression of NBS1 and Snail predicts metastasis of head and neck cancer. Oncogene. (2007) 26(10):1459–67. doi: 10.1038/sj.onc.1209929
22. Wu CY, Lin CT, Wu MZ, Wu KJ. Induction of HSPA4 and HSPA14 by NBS1 overexpression contributes to NBS1-induced in vitro metastatic and transformation activity. J Biomed Sci. (2011) 18(1):1. doi: 10.1186/1423-0127-18-1
23. Du X, Wang S, Liu X, He T, Lin X, Wu S, et al. MiR-1307-5p targeting TRAF3 upregulates the MAPK/NF-κB pathway and promotes lung adenocarcinoma proliferation. Cancer Cell Int. (2020) 20:502. doi: 10.1186/s12935-020-01595-z
24. Wang PX, Zhang XJ, Luo P, Jiang X, Zhang P, Guo J, et al. Hepatocyte TRAF3 promotes liver steatosis and systemic insulin resistance through targeting TAK1-dependent signalling. Nat Commun. (2016) 7:10592. doi: 10.1038/ncomms10592
25. Fabre J, Giustiniani J, Garbar C, Antonicelli F, Merrouche Y, Bensussan A, et al. Targeting the tumor microenvironment: the protumor effects of IL-17 related to cancer type. Int J Mol Sci. (2016) 17(9):1433. doi: 10.3390/ijms17091433
26. Hurtado CG, Wan F, Housseau F, Sears CL. Roles for interleukin 17 and adaptive immunity in pathogenesis of colorectal cancer. Gastroenterology. (2018) 155(6):1706–15. doi: 10.1053/j.gastro.2018.08.056
27. Zychlinsky Scharff A, Rousseau M, Lacerda Mariano L, Canton T, Consiglio CR, Albert ML, et al. Sex differences in IL-17 contribute to chronicity in male versus female urinary tract infection. JCI Insight. (2019) 5(13):e122998. doi: 10.1172/jci.insight.122998
28. Ma HY, Yamamoto G, Xu J, Liu X, Karin D, Kim JY, et al. IL-17 signaling in steatotic hepatocytes and macrophages promotes hepatocellular carcinoma in alcohol-related liver disease. J Hepatol. (2020) 72(5):946–59. doi: 10.1016/j.jhep.2019.12.016
29. Hou PP, Luo LJ, Chen HZ, Chen QT, Bian XL, Wu SF, et al. Ectosomal PKM2 promotes HCC by inducing macrophage differentiation and remodeling the tumor microenvironment. Mol Cell. (2020) 78(6):1192–1206.e10. doi: 10.1016/j.molcel.2020.05.004
30. Crispe IN, Giannandrea M, Klein I, John B, Sampson B, Wuensch S. Cellular and molecular mechanisms of liver tolerance. Immunol Rev. (2006) 213:101–18. doi: 10.1111/j.1600-065X.2006.00435.x
31. Couzin-Frankel J. Breakthrough of the year 2013. Cancer immunotherapy. Science. (2013) 342(6165):1432–3. doi: 10.1126/science.342.6165.1432
32. Alisa A, Ives A, Pathan AA, Navarrete CV, Williams R, Bertoletti A, et al. Analysis of CD4+ T-cell responses to a novel alpha-fetoprotein-derived epitope in hepatocellular carcinoma patients. Clin Cancer Res. (2005) 11(18):6686–94. doi: 10.1158/1078-0432.Ccr-05-0382
33. Lee JH, Lee JH, Lim YS, Yeon JE, Song TJ, Yu SJ, et al. Adjuvant immunotherapy with autologous cytokine-induced killer cells for hepatocellular carcinoma. Gastroenterology. (2015) 148(7):1383–91.e6. doi: 10.1053/j.gastro.2015.02.055
Keywords: Immune-related genes, hepatocellular carcinoma, prognostic risk score, immunotherapy, predictor
Citation: Cai Q, Duan J and Ding L (2022) Prognostic model of immune-related genes for patients with hepatocellular carcinoma. Front. Surg. 9:819491. doi: 10.3389/fsurg.2022.819491
Received: 21 November 2021; Accepted: 30 June 2022;
Published: 21 July 2022.
Edited by:
Alex Giakoustidis, Aristotle University of Thessaloniki, GreeceReviewed by:
Long Liu, First Affiliated Hospital of Zhengzhou University, ChinaYuting He, First Affiliated Hospital of Zhengzhou University, China
© 2022 Cai, Duan and Ding. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qun Cai quncai@zju.edu.cn
Specialty Section: This article was submitted to Surgical Oncology, a section of the journal Frontiers in Surgery