- 1Nanjing Drum Tower Hospital, Medicine School of Southeast University, Nanjing, China
- 2School of Life Science and Technology, Southeast University, Nanjing, China
- 3Department of General Surgery, Affiliated Drum Tower Hospital, Medical School, Nanjing University, Nanjing, China
Background: Hepatocellular carcinoma (HCC) is one of the malignant tumors with high morbidity and mortality worldwide. Immunotherapy has emerged as an increasingly important cancer treatment modality. However, the potential relationship between immune genes and HCC still needs to be explored. The purpose of this study is to construct a new prognostic risk signature to predict the prognosis of HCC patients based on the expression of immune-related genes (IRGs) and explore its potential mechanism.
Methods: We analyzed the gene expression data of 332 HCC patient samples and 46 adjacent normal tissues samples (Solid Tissue Normal including cirrhotic tissue) in The Cancer Genome Atlas (TCGA) database and clinical characteristics. We analyzed the gene expression data, identified differentially expressed IRGs in HCC tissues, filtered IRGs with prognostic value to construct an IRG signature, and classified patients into high and low gene expression groups based on the expression of IRGs in their tumor tissues. We also investigated the potential molecular mechanisms of IRGs through a bioinformatics approach using Protein-Protein Interaction (PPI) network, Kyoto Encyclopedia of Genes and Genomes (KEGG) database analysis and Gene Ontology (GO) database analysis. Differentially expressed IRGs associated with significant clinical outcomes (SIRGs) were identified by univariate Cox regression analysis. An immune-related risk score model (IRRSM) was established based on Lasso Cox regression analysis and multivariate Cox regression analysis. Based on the IRRSM, the immune score of the patients was calculated, and the patients were divided into high-risk and low-risk patients according to the median score, and the differences in survival between the two groups were compared. Then, the correlation analysis between the IRRSM and clinical characteristics was performed, and the IRRSM was validated using the International Cancer Genome Consortium (ICGC) database.
Results: The IRRSM was eventually constructed and confirmed to be an independent prognostic model for HCC patients. The IRRSM was shown to be positively correlated with the infiltration of four types of immune cells.
Conclusion: Our results showed that some SIRGs have potential value for predicting the prognosis and clinical outcomes of HCC patients. IRGs affect the prognosis of HCC patients by regulating the tumor immune microenvironment (TIME). This study provides a new insight for immune research and treatment strategies in HCC patients.
Introduction
Hepatocellular carcinoma (HCC), the most common primary liver cancer (accounting for over 80% of cases), is one of the malignant tumors with high morbidity and mortality worldwide. The most recent epidemiological data reveal that liver cancer ranked sixth in incidence, with 906,000 new cases, and third in cancer deaths, with 830,000 deaths, worldwide in 2020. Currently, surgical resection remains the only radical treatment for HCC. However, the surgical resection success rate is low, and the postoperative recurrence rate is high (1, 2), which severely limits the beneficial effect of surgical treatment.
HCC is the result of alterations in the interaction of a variety of cellular and molecular factors, including genetic, epigenetic, transcriptomic, and metabolic, that ultimately affect interactions in the tumor microenvironment. In particular, the tumor immune microenvironment (TIME) plays a very important role in the occurrence and treatment of HCC. The formation and development of tumors appears to be the result of the interaction between hepatocytes or HCC cells and immune cells. The immune system can recognize cancer cells and activate the immune response to eliminate tumor cells (3). Immunotherapy has become a highly promising therapeutic modality for the treatment of various cancers, including HCC (4). Concurrently, immune checkpoint inhibitors have become the most promising strategies in cancer immunotherapy (5, 6).
As a treatment modality for cancer, immunotherapy offers a promising treatment method for patients with advanced tumor. In particular, HCC is a highly heterogeneous tumor, which is closely related to immunity. Studies have revealed that some components of the TIME play key roles in the occurrence, development, invasion, metastasis and drug resistance of tumors, such as infiltrating immune cells, secreted cytokines/chemokines and other components. Research on the relationship of immunity and the development of HCC is very important (7). In the past few decades, immune-related genes (IRGs) have been shown to play a critical role in the occurrence, development and prognosis of various tumors (8, 9). IRGs can better predict the prognosis of patients and also provide a reference for early treatment of tumors. Although numerous studies have examined the correlation between IRGs and the prognosis of HCC patients, and there are differences in predictors in different studies. Therefore, additional prospective studies and basic experiments are still needed to verify the correlation between IRGs and the prognosis of HCC patients.
Accordingly, we designed this study to further investigate the clinical value of TIME and IRGs in the assessment of the prognosis of HCC patients. We obtained and identified HCC-related differentially expressed IRGs in tumor tissue samples, as well as the corresponding clinicopathological characteristics of the patients, and further evaluated the relationship between IRGs and overall survival (OS). We undertook to develop an IRGs-based immune-related risk score model (IRRSM) to predict the prognosis of HCC patients. This study revealed the likely mechanism of the role of IRGs in the progression of HCC, and established a suitable and accurate model providing a fresh perspective for clinical decision-making.
Methods
Data Downloading and Preprocessing
Total RNA‐sequencing data of 368 HCC tissue samples and 50 adjacent normal tissues samples were obtained from The Cancer Genome Atlas (TCGA) database (https://cancergenome.nih.gov/). After excluding patients with a shorter survival time (the patients with OS ≤90 days were excluded as these patients probably died of some unpredictable factors, such as severe infections, hemorrhage and liver failure), 332 hepatocellular carcinoma tissues and 46 adjacent tissues were eventually included. RNA expression sequencing data and patient clinical information of 237 HCC samples were also obtained from the International Cancer Genome Consortium (ICGC) database (https://dcc.icgc.org/) for validation.
The Analysis of Differential Genes and Differential IRGs
The limma package in the R software (https://bioconductor.Org/packages/release/bioc/html/limma.html) was used to screen differentially expressed genes (DEGs) between HCC and adjacent non-tumor tissues. The threshold used to screen for DEGs was set as follows: “False Discovery Rate (FDR) < 0.05 and log2| FC |> 1”. We obtained 192 IRGs from the ImmPort database (https://immport.niaid.nih.gov), and through the intersection with the identified DEGs, a total of 102 IRGs were ultimately obtained. In order to investigate the interaction between these genes, we mapped the Protein-Protein Interaction (PPI) network of these immune genes using the STRING online database. The CytoHubba in Cytoscape software version 3.8.2 is used to display the PPI network results. Additionally, we perform gene enrichment analysis using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases to investigate the potential molecular mechanism of the roles of these IRGs. The GO and KEGG pathways were identified using the cluster profiler, org.Hs.eg.db and enrichplot packages in the R software.
Construction of the IRG Signature
We performed a single-factor Cox regression analysis of the IRGs to determine the relationship between the expression of IRGs and the OS of HCC patients. In order to avoid overfitting and delete highly related genes, Lasso Cox regression was performed using the survival and glmnet package in the R software. Genes detected using the Lasso algorithm were evaluated by stepwise multivariate Cox regression analysis. Risk scores were obtained based on genes expression multiplied by a linear combination of the regression coefficient obtained by the multivariate Cox regression analysis. Patients were categorized into high-risk and low-risk groups according to the median risk score. The Kaplan–Meier analysis was performed to compare OS between high-risk and low-risk groups using the survival package in R.
We verified the model internally and externally, and subsequently tested the feasibility of the model through the clinical decision-making curve. The CIBERSORT (http://cibersortx.stanford.edu/) database was used to assess the levels of tumor infiltrating immune cells, which contains immune-related scores for various genes. Therefore, we scored according to the number of IRGs expressed in each sample, and eventually obtained the immune score of each sample. We determined the correlation between IRRSM and immune cell infiltration.
In order to confirm the relationship between IRGs associated with significant clinical outcomes (SIRGs) and the clinicopathological characteristics of HCC patients, we analyzed the correlation between the IRRSM and clinicopathological characteristics. The Tumor, Nodes and Metastasis (TNM) staging system is a well-recognized and the most commonly used system for classifying the extent of tumor cancers. Cancer stage defined according to the TNM classification system takes into account overall tumor burden and metastatic extent. Likewise, microvascular invasion and histopathological grading are closely associated with patient prognosis and treatment response. In order to explore the potential molecular mechanism SIRGs, we construct a regulatory network between SIRGs and transcription factors (TFs) based on NetworkAnalyst 3.0 (www.networkanalyst.ca).
Statistical Analysis
In order to predict the median survival and 1-, 3-, 5-year survival probability of HCC patients, we used the rms package in the R software to construct a nomogram. To verify the independent predictive ability of the nomogram model, the receiver operating characteristic (ROC) curve was drawn using the survival ROC package in the R software, and the models of different single clinical indicators and nomogram were used to evaluate the prognostic predictive ability of patients. We also showed the degree of calibration of the model through the calibration curve. Univariate Cox regression analysis, Lasso regression analysis and multivariate regression analysis were used to confirm SIRGs. The Kaplan-Meier curve was used to estimate the OS of the high-risk and low-risk groups and verify the independent prognostic factors of HCC patients. A radargram was drawn by the fmsb package in the R software, which displays the relationship of the IRRSM and level of immune cell infiltration. All statistical analyses were performed using SPSS26.0 software (IBM Corporation, Armonk, USA). P < 0.05 is considered statistically significant. All experiments were performed 3-4 times to reduce the error.
Results
Differential Gene Expression Analysis
We acquired the RNA-seq gene expression data and clinicopathological data of a total of 418 HCC patients with tumor tissues or adjacent normal tissues from TCGA database. In addition, after excluding patients with incomplete clinicopathological data and survival time of less than 90 days (36 HCC tissue samples and 4 adjacent normal tissues), patients-related data were extracted for the remaining 378 samples. Ultimately, the expression data of 332 HCC tissue samples and 46 adjacent normal tissues and their respective clinicopathological characteristics were included in this study. After conversion of 34,428 Ensembl human gene ID into gene names, 2,024 DEGs between HCC tissues and adjacent normal tissues, with threshold set as: (log2| FC |>1 and FDR <0.05), were identified using the limma package, of which 961 were downregulated and 1,063 were upregulated (Figure 1A). The top 20 upregulated and downregulated mRNAs of DEGs were visualized by heatmaps as shown in Figure 1B. We found 102 IRGs, of which 67 were downregulated and 35 were upregulated genes (Figure 1C). The expression level of IRGs is illustrated by a heatmap (Figure 1D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766072/fonc-11-766072-HTML/image_m/fonc-11-766072-g001.jpg)
Figure 1 Differentially expressed HCC genes and immune-related genes. Volcano plot (A) and heatmap (B) illustrated the differentially expressed genes between HCC tissues and para-carcinoma tissues. Immune-related genes (IRGs) with different expression levels were showed in the volcano plot (C) and heatmap (D). FDR < 0.05 and log2FC |> 1.
The Characteristics of IRGs
To investigate the mutual interaction among IRGs, a PPI network was constructed, as shown in Figure 2A. In this PPI network, FPR2, PMCH, RLN3, PNOC, OPRK1, PF4, CCR8, CXCL5, PPBP, CXCR5 were regarded as the hub genes among the 102 IRGs (Figure 2B). We performed GO term, KEGG pathway, and functional enrichment analyses to explore the potential biological functions of the common immune-related DEGs. As shown in the Figure 2C, “response to molecule of bacterial origin”, “secretory granule lumen” and “receptor ligand activity” were the most enriched terms in the biological processes (BP), cellular components (CC) and molecular functions (MF) categories, respectively. “Cytokine−cytokine receptor interaction” was identified to be the most enriched among the KEGG pathways of the IRGs (Figure 2D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766072/fonc-11-766072-HTML/image_m/fonc-11-766072-g002.jpg)
Figure 2 The functional enrichment analysis of differentially expressed IRGs. PPI network (A) of IRGs and the hub IRGs (B). The top pathways of IRGs were shown in biological process (BP), cellular component (CC), molecular function (MF) (C), and top KEGG pathways (D).
Clinical Outcomes of the High-Risk and the Low-Risk Groups
In order to establish the IRRSM, we identified 8 statistically significant SIRGs (including six high-risk IRGs and two low-risk IRGs) among the 21 SIRGs of the univariate Cox regression analysis (Table 1) using Lasso regression analysis and multivariate Cox regression analysis (Table 2). According to the IRRSM, the HCC patients were divided into the high-risk group and the low-risk group (Figure 3A). The mortality of patients with the higher risk scores was significantly higher than that of patients with the lower risk scores (Figure 3B). There were distinct differences in the expression levels of the 8 genes between the high-risk and low-risk groups (Figure 3C). In addition, in the IRRSM, the survival probability of the low-risk group was significantly higher than that of the high-risk group (Figure 3D). To further verify the predictive capacity of the IRG signature, 236 HCC patients from the ICGC database were used for external validation. To this end, we divided the patients into high-risk and low-risk groups. Consistent with the result obtained with TCGA dataset, the K‐M analysis showed that patients in the high‐risk group had a worse prognosis than those in the low‐risk group (Figure 3E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766072/fonc-11-766072-HTML/image_m/fonc-11-766072-g003.jpg)
Figure 3 IRRSM was created by SIRGs (A–C). The distribution of high-risk group and low-risk group (A). Survival status of the high-risk and the low-risk group (B). The heatmap of the expression levels of twelve SIRGs (C); Survival curve for the high risk and low risk groups of HCC patients in the TCGA database (D). Survival curve for the high risk and low risk groups of HCC patients for validation in ICGC database (E).
Development and Validation of Prognostic Signature
The Child Pugh score, histologic grade, TMN stage (tumor, metastasis, node), Risk score and microvascular invasion degree were included in TCGA cohort. Histologic grade, TMN stage and Risk score were included in ICGC cohort. A total of 332 patients in TCGA cohort and 236 patients in the ICGC cohort with complete clinical information were included. The hazard ratio (HR) was calculated and expressed in a forest plot (Figure 4A). According to the multivariate independent prognostic analysis results of TCGA cohort, a prognostic nomogram including significant factors (P-value < 0.05) was developed and constructed with the rms package in the R software (Figure 4B). Receiver operating characteristic (ROC) curves of the nomogram model, risk score model, histologic grade model, and TMN stage model were plotted using survival ROC package in the R software (Figures 4C, D). The area under the ROC curve (AUC) was also provided. Additionally, we compared the calibration and the discrimination inside and outside the model (Figures 4E, F). The calibration curves were also plotted and revealed the good prognostic prediction efficacy of the model (Figures 4G–L). The assessment of the Clinical Impact of Risk Prediction Models with Decision Curves is shown in Figures 4M, N.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766072/fonc-11-766072-HTML/image_m/fonc-11-766072-g004.jpg)
Figure 4 Independent prognostic prediction analysis (A). The clinical features, including Child Pugh class, histologic grade, TMN stage and microvascular invasion, were combined with Risk Score to construct the nomogram model (B). The survival predicted ROC curves of four distinct models, namely, Nomogram model, TNM stage model, Risk score model and histologic grade model were compared in TCGA and ICGC database (C, D). The calibration curve diagrams of the training set (E) and the validation set (F) have good agreement between the predicted probability and the actual probability, both S: p>0.05. The calibration curves of 1,2,3-year survival (G–L). The dotted line represented the ideal prediction model, and the red solid line represented the observed model in TCGA database (G–I). The calibration curves of 1,2,3-year survival. The dotted line represented the ideal prediction model, and the red solid line represented the observed model in ICGC database (J–lL). Decision Curve Analysis model: A clinical decision-making model based on the final multiple risk factors in TCGA and ICGC database (M, N).
Relationship of the IRG Signature and Clinical Variables
We also examined the correlations between clinical factors and risk scores calculated from the IRG signature. As shown in Figure 5, the results indicated that higher TNM stage, microvascular invasion and histologic grade tended to correlate with higher risk scores. The correlations between the eight IRGs and clinical tumor histologic grade were also assessed. As shown in Supplementary Figure 1, higher expression of NR0B1 and IL20RA were associated with higher tumor histologic grade (Supplementary Figures 1E, F), and lower expression of ESRRG were associated with higher tumor histologic grade (Supplementary Figure 1G).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766072/fonc-11-766072-HTML/image_m/fonc-11-766072-g005.jpg)
Figure 5 Relationships between risk score and clinical characteristics. The figures demonstrated that higher of TNM stage (A), microvascular invasion (B) and histologic grade (C) tended to correlated with higher risk scores.
The Analysis of the Immune Status of the High-Risk and Low-Risk Groups
The functional enrichment analysis of the IRG signature was performed by gene set enrichment analysis (GSEA) to compare the high‐risk and the low‐risk groups in HCC. The KEGG pathway enrichment analysis revealed that “cytokine-cytokine receptor interaction”, “MAPK signaling pathway”, “neuroactive ligand receptor interaction”, “olfactory transduction”, “pathway in cancer”, “regulation of actin cytoskeleton” were significantly related to the high‐risk group (Figures 6A–F). In order to further examine whether the IRRSM accurately reflects the TIME, we analyzed the relationships between the IRRSM and immune cells infiltration (Figure 6G). We found that activated dendritic cells (Figure 6H), resting dendritic cells (Figure 6I), non-activated (M0) macrophages (Figure 6J) and CD8 T cells (Figure 6K) showed positive correlation with risk score, whereas naïve B cells (Figure 6L), classically activated (M1) macrophages (Figure 6M), alternatively activated (M2) macrophages (Figure 6N), activated mast cells (Figure 6O), resting mast cells (Figure 6P) and gamma delta T cells (Figure 6Q) showed the opposite results. These results motivated us to further pursue the research on the underlying functions and mechanisms in future studies.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766072/fonc-11-766072-HTML/image_m/fonc-11-766072-g006.jpg)
Figure 6 The most significantly enriched KEGG pathways between high- and low-risk patients. KEGG: cytokine-cytokine receptor interaction (A); MAPK signaling pathway (B); neuroactive ligand receptor interaction (C); olfactory transduction (D); pathway in cancer (E); regulation of actin cytoskeleton (F). Relationships between the IRRSM and immune cells. The relationships between IRRSM and immune cells infiltration: radar plot (G); dendritic cells activated (H); dendritic cells resting (I); macrophages M0 (J); T cells CD8 (K); B cells naive (L); macrophages M1 (M); macrophages M2 (N); mast cells activated (O); mast cells resting (P); T cells gamma delta (Q).
Mechanisms of the Immune‐Related Signature
We explored the potential regulatory mechanisms of the eight immune‐related genes that made up the signature, which could reflect the regulatory mechanisms of the signature. Based eight immune‐related genes, we constructed a relevant regulatory network between TFs and immune‐related genes by NetworkAnalyst 3.0 (www.networkanalyst.ca) (Figure 7). In terms of the figure, we can see that GIP, GNRH2, FGF9, CHGA and IL20RA were included in the regulatory network, while there was no regulatory relationship between the RAET1E, NR0B1, ESRRG and TFs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766072/fonc-11-766072-HTML/image_m/fonc-11-766072-g007.jpg)
Figure 7 Mechanisms of the immune‐related risk signature. Transcription factors‐based regulatory network with the eight immune‐related genes in the signature.
Discussion
With the rapid development of next-generation sequencing technologies, our understanding of transcriptional alterations in HCC has greatly improved. The molecular landscape of synchronous HCC reveals genetic heterogeneity and predicts the prognosis of patients, thereby providing an innovative and effective way to identify promising markers and targets for HCC diagnosis and treatment (10–14). According to the survival prediction of patients with prognostic biomarkers, individualized treatment can be better developed in clinical practice (15, 16). Compared with a single gene, a multigene-based model has been shown to be more robust and precise in diagnosis and prognosis in many cancers (17, 18). Thus, it provides new options for prognostic prediction and treatment of patients.
Since generally people have an in depth understanding of the importance of immune activity in tumorigenesis, progression and prognosis, the emerging cancer immunotherapy approaches have always had an attention-catching position in the field of cancer treatment (19). Over the past few decades, numerous targeted cancer and immunotherapeutic drugs have been and are currently being developed, with many already used in routine clinical care, including pembrolizumab, ipilimumab and nivolumab, which have been approved for clinical application and have achieved good results (20, 21). Therefore, the discovery of more promising and sensitive immune-related biomarkers, such as IRGs is necessary to improve the efficacy of immunotherapy. In this study, we developed a novel prediction method based on IRGs in HCC patients, which may also independently predict the OS and median survival of HCC patients. Most importantly, we also investigated the molecular mechanisms related with the IRG signature using the bioinformatics tools and analysis of the correlations between our IRG signature and immune cell infiltrations to evaluate immune cells infiltration. These findings indicate that some clinical signature is of great significance for predicting the prognosis of HCC patients and exploring potential immunotherapy targets in HCC.
In this study, a novel IRG signature was identified to classify HCC patients according to the risk score. We integrated multiple genes into a single signature through multivariate Cox regression analysis modeling. The novel IRG signature was verified to be a significant predictor both in TCGA and ICGC cohorts. The eight genes included in the identified IRG signature were CHGA, RAET1E, FGF9, GIP, NR0B1, IL20RA, ESRRG and GNRH2. The heatmaps were plotted to illustrate the risk of the identified genes based on gene expression profile. As shown in the heatmaps, the expression of the twelve altered genes increased together with the risk score in both cohorts.
The signature was constructed by eight IRGs that were able to predict prognosis of HCC patients. In this study, all eight TRGs were shown to be independent prognostic markers for HCC patients.
Expression of CHGA has been reported to be associated with prognosis in colorectal cancer (22–24), and is generally recognized as the main biomarker for neuroendocrine neoplasms (25, 26). Recently, CHGA has been proposed as an early diagnosis biomarker for gastric cancer (27), prostate cancer (28), and pancreatic neuroendocrine tumors (29). However, there is no study about CHGA expression and early diagnosis in HCC. ESRRG has been identified as a tumor suppressor gene in several cancers (30), which is consistent with our findings that ESRRG is lowly expressed in HCC tissue and can predict the prognosis of HCC patients, indicating that ESRRG could regulate the occurrence and development of HCC. FGF9 was first identified from the secretions of human glioma MCF-G1 cells, and found to promote NIH-3T3 cell line malignant transformation, suggesting that FGF9 may be a tumor-promoting factor (31). Indeed, FGF9 has been found to be highly expressed in HCC, and significantly associated with the proliferation of cancer cells (32–34). Several alterations in FGF-signaling have been found to affect liver carcinogenesis (35). Aberrant expression of the endocrine FGF19 and its high affinity FGFR4 receptor contributes to HCC progression (36). GIPR is found in many cancers, including the majority of pancreatic, ileal, and bronchial neuroendocrine tumors (NETs). The high GIPR content in a broad spectrum of cancers, makes GIPR a potential target for cancer diagnostic and therapeutic purposes (37, 38). Expression of the gonadotropin-releasing hormone homologue peptides GNRH2 and their receptor GNRHR has been demonstrated in various malignancies. Studies have shown that GNRH2 is closely related to the occurrence of multiple cancers, and to affect the prognosis of patients. Recent studies have provided evidence that IL20RA signaling regulates the development of certain cancers. Another study found that IL20RA expression was elevated in breast cancer and colorectal cancer (39). Epigenetic modification of NR0B1 leads to its ectopic activation in Ewing’s sarcoma and lung cancer, enabling it to promote cancer cell proliferation (40). Additionally, it can also inhibit the proliferation of liver cancer cells by suppressing the transcriptional activity of β-catenin (41). It has been reported that RAET1E might impair NKG2D-mediated NK cell cytotoxicity to liver cancer cells (42).
This study established a signature for calculating the risk score based on the expression levels of eight genes, which predicted the median survival and OS of HCC patients. Moreover, this study also provides an innovative analysis approach to identify significant biomarkers. Also, through comprehensive data analysis, reliable results with excellent repeatability can be obtained, which can be verified with external data and clinical outcomes.
Conclusion
In this paper, we analyzed the roles of SIRGs in the prediction and evaluation of the clinical prognosis of HCC patients and verified the predictive value of some SIRGs. We establish a reliable and accurate IRRSM to predict the prognosis of HCC patients, which offer a new perspective for treatments by early intervention.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
DY, JP, and QX designed the study. DY and YL analyzed the data and wrote the manuscript. JP and QX provided technical expertise and support. The manuscript was prepared and reviewed by DY, YL, GL, and BS. 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.
Acknowledgments
We sincerely thank the researchers and study participants for their contributions toward this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.766072/full#supplementary-material.
Supplementary Figure 1 | Relationships between clinical tumor histologic grade with expression of the eight immune‐related genes. CHGA (A), RAET1E (B), FGF9(C), GIP (D), NR0B1 (E), IL20RA (F), ESRRG (G), GNRH2 (H).
References
1. Allaire M, Goumard C, Lim C, Le Cleach A, Wagner M, Scatton O. New Frontiers in Liver Resection for Hepatocellular Carcinoma. JHEP Rep (2020) 2:100134. doi: 10.1016/j.jhepr.2020.100134
2. Gunasekaran G, Bekki Y, Lourdusamy V, Schwartz M. Surgical Treatments of Hepatobiliary Cancers. Hepatology (2021) 73(Suppl 1):128–36. doi: 10.1002/hep.31325
3. Zhang Q, Lou Y, Bai X-L, Liang T-B. Immunometabolism: A Novel Perspective of Liver Cancer Microenvironment and Its Influence on Tumor Progression. World J Gastroenterol (2018) 24:3500–12. doi: 10.3748/wjg.v24.i31.3500
4. Donisi C, Puzzoni M, Ziranu P, Lai E, Mariani S, Saba G, et al. Immune Checkpoint Inhibitors in the Treatment of HCC. Front Oncol (2020) 10:601240. doi: 10.3389/fonc.2020.601240
5. Singh S, Hassan D, Aldawsari HM, Molugulu N, Shukla R, Kesharwani P. Immune Checkpoint Inhibitors: A Promising Anticancer Therapy. Drug Discovery Today (2020) 25:223–9. doi: 10.1016/j.drudis.2019.11.003
6. Le Saux O, Lounici Y, Wajda P, Barrin S, Caux C, Dubois B, et al. Neoadjuvant Immune Checkpoint Inhibitors in Cancer, Current State of the Art. Crit Rev Oncol Hematol (2021) 157:103172. doi: 10.1016/j.critrevonc.2020.103172
7. Ruf B, Heinrich B, Greten TF. Immunobiology and Immunotherapy of HCC: Spotlight on Innate and Innate-Like Immune Cells. Cell Mol Immunol (2021) 18:112–27. doi: 10.1038/s41423-020-00572-w
8. Gnjatic S, Bronte V, Brunet LR, Butler MO, Disis ML, Galon J, et al. Identifying Baseline Immune-Related Biomarkers to Predict Clinical Outcome of Immunotherapy. J Immunother Cancer (2017) 5:44. doi: 10.1186/s40425-017-0243-4
9. Shi X, Li R, Dong X, Chen AM, Liu X, Lu D, et al. IRGS: An Immune-Related Gene Classifier for Lung Adenocarcinoma Prognosis. J Transl Med (2020) 18:55. doi: 10.1186/s12967-020-02233-y
10. Xu Q, Wang Y, Huang W. Identification of Immune-Related lncRNA Signature for Predicting Immune Checkpoint Blockade and Prognosis in Hepatocellular Carcinoma. Int Immunopharmacol (2021) 92:107333. doi: 10.1016/j.intimp.2020.107333
11. Zhang R, Ye J, Huang H, Du X. Mining Featured Biomarkers Associated With Vascular Invasion in HCC by Bioinformatics Analysis With TCGA RNA Sequencing Data. BioMed Pharmacother (2019) 118:109274. doi: 10.1016/j.biopha.2019.109274
12. Zhang J, Chen G, Zhang J, Zhang P, Ye Y A. Construction of a Prognostic Model Based on Nine Immune-Related Genes and Identification of Small Molecule Drugs for Hepatocellular Carcinoma (HCC). Am J Transl Res (2020) 12:5108–30.
13. Li Z, Lin Y, Cheng B, Zhang Q, Cai Y. Identification and Analysis of Potential Key Genes Associated With Hepatocellular Carcinoma Based on Integrated Bioinformatics Methods. Front Genet (2021) 12:571231. doi: 10.3389/fgene.2021.571231
14. Dong S, Huang F, Zhang H, Chen Q. Overexpression of BUB1B, CCNA2, CDC20, and CDK1 in Tumor Tissues Predicts Poor Survival in Pancreatic Ductal Adenocarcinoma. Biosci Rep (2019) 39:BSR20182306. doi: 10.1042/BSR20182306
15. Sun Y, Wu L, Zhong Y, Zhou K, Hou Y, Wang Z, et al. Single-Cell Landscape of the Ecosystem in Early-Relapse Hepatocellular Carcinoma. Cell (2021) 184:404–421. doi: 10.1016/j.cell.2020.11.041
16. Hu X, Bao M, Huang J, Zhou L, Zheng S. Identification and Validation of Novel Biomarkers for Diagnosis and Prognosis of Hepatocellular Carcinoma. Front Oncol (2020) 10:541479. doi: 10.3389/fonc.2020.541479
17. Chen L, Li Z, Zeng T, Zhang Y-H, Li H, Huang T, et al. Predicting Gene Phenotype by Multi-Label Multi-Class Model Based on Essential Functional Features. Mol Genet Genomics (2021) 296:905–18. doi: 10.1007/s00438-021-01789-8
18. Gu X, Guan J, Xu J, Zheng Q, Chen C, Yang Q, et al. Model Based on Five Tumour Immune Microenvironment-Related Genes for Predicting Hepatocellular Carcinoma Immunotherapy Outcomes. J Transl Med (2021) 19:26. doi: 10.1186/s12967-020-02691-4
19. Kather JN, Halama N, Jaeger D. Genomics and Emerging Biomarkers for Immunotherapy of Colorectal Cancer. Semin Cancer Biol (2018) 52:189–97. doi: 10.1016/j.semcancer.2018.02.010
20. Tintelnot J, Stein A. Immunotherapy in Colorectal Cancer: Available Clinical Evidence, Challenges and Novel Approaches. World J Gastroenterol (2019) 25:3920–8. doi: 10.3748/wjg.v25.i29.3920
21. Ganesh K, Stadler ZK, Cercek A, Mendelsohn RB, Shia J, Segal NH, et al. Immunotherapy in Colorectal Cancer: Rationale, Challenges and Potential. Nat Rev Gastroenterol Hepatol (2019) 16:361–75. doi: 10.1038/s41575-019-0126-x
22. Kleist B, Poetsch M. Neuroendocrine Differentiation: The Mysterious Fellow of Colorectal Cancer. World J Gastroenterol (2015) 21:11740–7. doi: 10.3748/wjg.v21.i41.11740
23. Zeng Y-J, Lai W, Liu L, Wu H, Luo X-X, Wang J, et al. Prognostic Significance of Neuroendocrine Differentiation in Colorectal Adenocarcinoma After Radical Operation: A Meta-Analysis. J Gastrointest Surg (2014) 18:968–76. doi: 10.1007/s11605-014-2480-x
24. Dost Gunay FS, Kırmızı BA, Ensari A, İcli F, Akbulut H. Tumor-Associated Macrophages and Neuroendocrine Differentiation Decrease the Efficacy of Bevacizumab Plus Chemotherapy in Patients With Advanced Colorectal Cancer. Clin Colorectal Cancer (2019) 18:e244–50. doi: 10.1016/j.clcc.2018.12.004
25. Wang L, Li H, Li Z, Li M, Tang Q, Wu C, et al. Smoothened Loss is a Characteristic of Neuroendocrine Prostate Cancer. Prostate (2021) 81:508–20. doi: 10.1002/pros.24122
26. Galvin SG, Larraufie P, Kay RG, Pitt H, Bernard E, McGavigan AK, et al. Peptidomics of Enteroendocrine Cells and Characterisation of Potential Effects of a Novel Preprogastrin Derived-Peptide on Glucose Tolerance in Lean Mice. Peptides (2021) 140:170532. doi: 10.1016/j.peptides.2021.170532
27. Wen L, Han Z, Du Y. Identification of Gene Biomarkers and Immune Cell Infiltration Characteristics in Rectal Cancer. J Gastrointest Oncol (2021) 12:964–80. doi: 10.21037/jgo-21-255
28. Mather RL, Parolia A, Carson SE, Venalainen E, Roig-Carles D, Jaber M, et al. The Evolutionarily Conserved Long Non-Coding RNA LINC00261 Drives Neuroendocrine Prostate Cancer Proliferation and Metastasis via Distinct Nuclear and Cytoplasmic Mechanisms. Mol Oncol (2021) 15:1921–41. doi: 10.1002/1878-0261.12954
29. Dugnani E, Sordi V, Pellegrini S, Chimienti R, Marzinotto I, Pasquale V, et al. Gene Expression Analysis of Embryonic Pancreas Development Master Regulators and Terminal Cell Fate Markers in Resected Pancreatic Cancer: A Correlation With Clinical Outcome. Pancreatology (2018) 18:945–53. doi: 10.1016/j.pan.2018.09.006
30. Kim J-H, Choi Y-K, Byun J-K, Kim M-K, Kang YN, Kim SH, et al. Estrogen-Related Receptor γ is Upregulated in Liver Cancer and Its Inhibition Suppresses Liver Cancer Cell Proliferation via Induction of P21 and P27. Exp Mol Med (2016) 48:e213. doi: 10.1038/emm.2015.115
31. Naruo K, Seko C, Kuroshima K, Matsutani E, Sasada R, Kondo T, et al. Novel Secretory Heparin-Binding Factors From Human Glioma Cells (Glia-Activating Factors) Involved in Glial Cell Growth. Purification and Biological Properties. J Biol Chem (1993) 268:2857–64. doi: 10.1016/S0021-9258(18)53852-7
32. Cui G, Shao M, Gu X, Guo H, Zhang S, Lu J, et al. The Value of as a Novel Biomarker in the Diagnosis of Prostate Cancer. Artif Cells Nanomed Biotechnol (2019) 47:2241–5. doi: 10.1080/21691401.2019.1620250
33. He X, Chen S-Y, Yang Z, Zhang J, Wang W, Liu M-Y, et al. miR-4317 Suppresses Non-Small Cell Lung Cancer (NSCLC) by Targeting Fibroblast Growth Factor 9 (FGF9) and Cyclin D2 (Ccnd2). J Exp Clin Cancer Res (2018) 37:230. doi: 10.1186/s13046-018-0882-4
34. Wang R, Sun Y, Yu W, Yan Y, Qiao M, Jiang R, et al. Downregulation of miRNA-214 in Cancer-Associated Fibroblasts Contributes to Migration and Invasion of Gastric Cancer Cells Through Targeting FGF9 and Inducing EMT. J Exp Clin Cancer Res (2019) 38:20. doi: 10.1186/s13046-018-0995-9
35. Paur J, Valler M, Sienel R, Taxauer K, Holzmann K, Marian B, et al. Interaction of FGF9 With FGFR3-IIIb/IIIc, a Putative Driver of Growth and Aggressive Behaviour of Hepatocellular Carcinoma. Liver Int (2020) 40:2279–90. doi: 10.1111/liv.14505
36. Raja A, Park I, Haq F, Ahn S-M. FGF19- Signaling in Hepatocellular Carcinoma. Cells (2019) 8:536. doi: 10.3390/cells8060536
37. Waser B, Rehmann R, Sanchez C, Fourmy D, Reubi JC. Glucose-Dependent Insulinotropic Polypeptide Receptors in Most Gastroenteropancreatic and Bronchial Neuroendocrine Tumors. J Clin Endocrinol Metab (2012) 97:482–8. doi: 10.1210/jc.2011-2454
38. Sherman SK, Carr JC, Wang D, O'Dorisio MS, O'Dorisio TM, Howe JR. Gastric Inhibitory Polypeptide Receptor (GIPR) is a Promising Target for Imaging and Therapy in Neuroendocrine Tumors. Surgery (2013) 154:1206–14. doi: 10.1016/j.surg.2013.04.052
39. Gao W, Wen H, Liang L, Dong X, Du R, Zhou W, et al. IL20RA Signaling Enhances Stemness and Promotes the Formation of an Immunosuppressive Microenvironment in Breast Cancer. Theranostics (2021) 11:2564–80. doi: 10.7150/thno.45280
40. Lu Y, Liu Y, Liao S, Tu W, Shen Y, Yan Y, et al. Epigenetic Modifications Promote the Expression of the Orphan Nuclear Receptor NR0B1 in Human Lung Adenocarcinoma Cells. Oncotarget (2016) 7:43162–76. doi: 10.18632/oncotarget.9012
41. Jiang H-L, Xu D, Yu H, Ma X, Lin G-F, Ma D-Y, et al. DAX-1 Inhibits Hepatocellular Carcinoma Proliferation by Inhibiting β-Catenin Transcriptional Activity. Cell Physiol Biochem (2014) 34:734–42. doi: 10.1159/000363038
Keywords: hepatocellular carcinoma, immune-related genes, immune-related risk score, prognosis, nomogram
Citation: Ye D, Liu Y, Li G, Sun B, Peng J and Xu Q (2021) A New Risk Score Based on Eight Hepatocellular Carcinoma- Immune Gene Expression Can Predict the Prognosis of the Patients. Front. Oncol. 11:766072. doi: 10.3389/fonc.2021.766072
Received: 28 August 2021; Accepted: 01 November 2021;
Published: 19 November 2021.
Edited by:
Maen Abdelrahim, Houston Methodist Research Institute, United StatesReviewed by:
Antonio Giovanni Solimando, University of Bari Aldo Moro, ItalyJuan Pablo Cerapio Arroyo, INSERM U1037 Centre de Recherche en Cancérologie de Toulouse (INSERM), France
Copyright © 2021 Ye, Liu, Li, Sun, Peng and Xu. 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: Qingxiang Xu, 2841958311@qq.com; Jin Peng, pengjin119@163.com