Parts of this article's content have been modified or rectified in:
Erratum: Construction of a tumor immune infiltration macrophage signature for predicting prognosis and immunotherapy response in liver cancer
Anmin Huang1,2,3†
Bei Lv4†
Yunjie Zhang1,2,3
Junhui Yang1,2,3Jie Li1,2,3Chengjun Li1,2,3
Zhijie Yu5*
Jinglin Xia1,6,7*- 1Key Laboratory of Diagnosis and Treatment of Severe Hepato-Pancreatic Diseases of Zhejiang Province, The First Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
- 2Translational Medicine Laboratory, The First Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
- 3The First Clinical Medical College, Wenzhou Medical University, Wenzhou, China
- 4Department of Radiation Oncology, Zhongshan Hospital, Fudan University, Shanghai, China
- 5Wenzhou Key Laboratory of Hematology, The First Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
- 6Department of Intervention, The First Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
- 7Liver Cancer Institute, Zhongshan Hospital, Fudan University, Shanghai, China
Liver cancer is an extraordinarily heterogeneous malignant disease. The tumor microenvironment (TME) and tumor-associated macrophages (TAMs) are the major drivers of liver cancer initiation and progression. It is critical to have a better understanding of the complicated interactions between liver cancer and the immune system for the development of cancer immunotherapy. Based on the gene expression profiles of tumor immune infiltration cells (TIICs), upregulated genes in TAMs and downregulated genes in other types of immune cells were identified as macrophage-specific genes (MSG). In this study, we combined MSG, immune subtypes, and clinical information on liver cancer to develop a tumor immune infiltration macrophage signature (TIMSig). A four-gene signature (S100A9, SLC22A15, TRIM54, and PPARGC1A) was identified as the TAM-related prognostic genes for liver cancer, independent of multiple clinicopathological parameters. Survival analyses showed that patients with low TIMSig had a superior survival rate than those with high TIMSig. Additionally, clinical immunotherapy response and TIMSig was observed as highly relevant. In addition, TIMSig could predict the response to chemotherapy. Collectively, the TIMSig could be a potential tool for risk-stratification, clinical decision making, treatment planning, and oncology immunotherapeutic drug development.
Introduction
The development is regulated by the dynamic tumor microenvironment (TME) which is comprised of a complex network of multipotent stromal cells, fibroblasts, blood vessels, endothelial cell precursors, immune cells, and various secreted factors such as cytokines (Anderson and Simon, 2020). Tumor immune infiltration cells (TIICs) have been demonstrated to be important in tumor proliferation, invasion, metastasis, and drug resistance through complex immune interaction with malignant cells (Binnewies et al., 2018; Lei et al., 2020). Recent studies have shown that TIICs are directly linked to the prognosis of cancer patients (Fridman et al., 2017; Giraldo et al., 2019). The temporal and spatial heterogeneity of TME in hepatocellular carcinoma (HCC) suggests that the disorder and imbalance of TME caused by various factors, such as DNA methylation and chromatin architecture, may be one of the most critical mechanisms of tumorigenesis and progression (Kurebayashi et al., 2018; Fu et al., 2019; Johnstone et al., 2020). At present, most studies on the occurrence and development of primary liver carcinoma focus on HCC, while little attention has been paid to cholangiocarcinoma (CCA) and combined hepatocellular-cholangiocarcinoma (cHCC-CCA). Herein, we integrate the information from all pathological types of liver cancer to establish a novel risk assessment model to evaluate the risk of liver cancer.
Macrophages that reside within the TME are known as tumor-associated macrophages (TAMs). As the predominant infiltrated TIICs, TAMs promote tumor development at multiple levels, such as accelerating genomic mutation, cultivating tumor stem cells, paving the way for metastasis, and taming the immune system (Wynn et al., 2013; Mantovani et al., 2017). There are two types of macrophages classically activated pro-inflammatory (M1) and alternatively activated anti-inflammatory (M2) phenotypes. TAMs can be polarized into M2 phenotype macrophages or M1 phenotype macrophages by different stimuli (Vitale et al., 2019; Yoon et al., 2019). M1 macrophages release various pro-inflammatory cytokines, active oxygen, and nitrogen oxide that drive tumor-killing activities (Qian and Pollard, 2010). In contrast, factors such as IL-10, TGF-β, and VEGF secreted by M2 macrophages promote tumor progression, angiogenesis, metastasis, and suppression of anti-tumor immunity (Lujambio et al., 2013). Previous studies on TAMs have suggested that high M2 macrophage infiltration is associated with a poor prognosis in cancers (Martínez et al., 2017; Väyrynen et al., 2021). Therefore, TAMs are correlated to the development of liver tumors and the clinical prognosis of patients, and a TAM-based signature could have potential application in predicting clinical outcomes and immunotherapeutic responses.
Immune checkpoint inhibitors (ICIs) have made an indelible mark in the field of cancer immunotherapy (Vaddepally et al., 2020). ICI immunotherapy covers a series of monoclonal antibodies designed to block the binding of immune checkpoints expressed on the surface of immune cells to their ligands, which eliminates the immunosuppression caused by immune checkpoints and revitalizes the function of T cells (Dyck and Mills, 2017). ICIs such as anti-cytotoxic Tlymphocyte associated protein 4 (CTLA4, CTLA-4), anti-programmed cell death protein 1 (PD1, PD-1), and anti-programmed cell death one ligand 1 (PDL1, PD-L1) have been approved for the clinical application of advanced malignancies, including melanoma, non-small-cell lung cancer, urothelial carcinoma, gastric cancer, liver cancer, and Hodgkin’s lymphoma (Wei et al., 2018; Perez-Ruiz et al., 2019). Although ICIs have demonstrated their efficiency and durability in the treatment of solid tumors, a great number of patients have limited benefits in terms of response and survival (Syn et al., 2017; Schoenfeld and Hellmann, 2020; Bagchi et al., 2021). At present, some research perspectives are that this phenomenon may be related to T cell depletion and mechanical factors of TME, which place restrictions on the efficacy of ICI treatment in cancer patients (Sacks et al., 2018; Wong et al., 2021). Therefore, new molecules and prognostic models should be explored to predict or improve the clinical response and application of ICI therapy.
This study developed a computational algorithm framework for identifying prognostic signatures with TAMs (Figure 1). We selected different immune datasets from the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/database) database, which contains bulk RNA-sequencing data of cell lines, primary animal cells, and primary human cells, for screening the tumor immune infiltration macrophage genes (TIMGs). We performed a systematic and thorough biomarker analysis and validation to locate and construct a risk predictive signature for the prognosis of liver cancer by taking into account multiple components of TME. Herein, we report a tumor immune infiltration macrophage signature (TIMSig), which can predict the prognosis and immunotherapy of liver cancer individuals and reflect the cellular functions of macrophages.
FIGURE 1
FIGURE 1. Workflow of TIMSig generation and validation.
Methods
Data collection
The transcriptome sequencing information and corresponding clinical data (n = 404) of liver cancer samples were obtained from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/). The gene expression data were given in log2 (x+1) transformed and multiple imputations with R package “mice” prior to analyses. Transcriptional profiles of macrophages from three datasets (GSE158792, GSE56755, and GSE75829) and 21 other immune cell profiles from 15 datasets (GSE23371, GSE27291, GSE157737, GSE27838, GSE37750, GSE8059, GSE51540, GSE149425, GSE104852, GSE52156, GSE155148, GSE83441, GSE42058, GSE106932, and GSE28726) were obtained from the GEO database, these profiles were used to perform differential gene analyses using GEO2R online tools with R package “limma”. The intersection of upregulated genes in macrophages (cut-off was set at logFC >1.5, p < 0.05) and downregulated genes (cut-off was set at logFC < -1.5, p < 0.05) in other immune cells was considered as macrophage-specific genes (MSG).
Identification and prognostic assessment of immune subtypes in liver cancer
CIBERSORT algorithm was performed to characterize 22 immune cell subtypes of immune landscape using the deconvolution strategy. The 404 liver cancer patients from TCGA were clustered according to immune landscape by using the R package “ConsensusClusterPlus”, and a consensus matrix was constructed to define the immune subtypes. The prognostic potential of each immune subtype was estimated by Kaplan–Meier (K-M) curves. Then, the association between the MSG and different immune subtypes of cancer patients was investigated by using the R package “limma”, with |logFC| > 1.5 and p < 0.05 considered as significantly different.
Construction and validation of the TIMSig
Cancer cases from the TCGA database were randomized into the training cohort or test cohort for the construction and validation of risk scores at the ratio of 1:1 (202 in the training set, 202 in the test set). In the training cohort, we used univariate Cox proportional hazards regression analysis and LASSO regression analysis to investigate the prognostic performance of candidate TIMGs. Then, TIMGs were determined by multivariate Cox proportional hazards regression analysis, and the risk score based on TIMGs was constructed. The samples were divided into the high-risk group and low-risk group based on the median risk score. The difference in overall survival (OS) between the high-risk group and low-risk group was performed using the log-rank test. Various statistical methods, including K–M curves, univariate Cox proportional hazards regression analysis, multivariate Cox proportional hazards regression analysis, and time-dependent receiver-operating characteristic (ROC) curves were performed to identify the association between TIMSig and OS.
Validation of the TIMSig in the GSE14520 cohort and ICI clinical samples
Another sample information of the GSE14520 dataset included information from 242 tumor samples, which was downloaded from the GEO database. Due to different platforms, only part of four TIMGs in the TILSig was covered by the other dataset. Therefore, the risk score is only based on part of the TIMGs. Moreover, we collected 24 PD-1/PD-L1 clinical treated cases with HCC from GSE140901 with available progression-free survival (PFS) records based on the nCounter PanCancer Immune Profiling Panel platform, which came from the National Taiwan University Hospital for estimation of the potential of TIMGs as immunotherapy biomarkers. The association between TIMSig and PD-1/PD-L1 immunotherapy response was investigated by paired t-test. ROC curve was performed to estimate the accuracy of TIMSig response to predict ICI treatment response. K–M curve was performed to evaluate the risk-stratification capability of TIMSig and PD-1/PD-L1 immunotherapy response.
Tumor immune microenvironment analyses
ESTIMATE algorithm was performed to calculate the scores of immune infiltration and tumor purity based on gene expression profile. The ESTIMATE analysis method is integrated into the “estimate” R package. Tumor Immune Estimation Resource 2.0 (TIMER, http://timer.comp-genomics.org/) is an excellent tool to assess the level of immune infiltration, applied to analyze the association between TIMGs and abundance of TIICs for estimation of the capability of TIMGs in predicting TIICs. In addition, the expression differences of TIMGs between cancer and normal samples were visualized.
Functional analyses
We carried out the R package “ggpubr” to investigate the function of macrophages in secreting cytokines by analyzing the association of the TIMSig and the expression level of corresponding cytokine molecules. Then, The Gene Set Enrichment Analysis (GSEA, ver. 4.1.0) was applied to analyze the organism function and cell development with TIMSig through the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The Human Protein Atlas (HPA, https://www.proteinatlas.org/) provides immunohistochemistry (IHC) information for human cancers, which applied the HPA database to evaluate the protein expression levels for TIMGs.
Statistical analyses
All statistical analyses were performed by R version 4.1.1, in addition to corresponding R packages and online tools. A p-value < 0.05 was considered statistically significant. Univariate Cox proportional hazards regression analysis and multivariate Cox proportional hazards regression analysis methods were integrated in the R package “survival” and “survminer”. LASSO Cox proportional regression analysis was performed with the R package “glmnet”. ROC curves were performed, and the area under curves (AUC) was applied to estimate the prognostic performance of TIMSig using the R package “survivalROC”. The relationship between the TIMSig and immune checkpoints was evaluated to verify the potential prognostic capacity of TIMSig for ICI immunotherapy. Drug sensitivity analysis was applied with the R package “pRRophetic”.
Results
Identification of MSG
To characterize the gene expression pattern in various immune cell subtypes, we retrieved 18 expression profiles of immune cells from the GEO database. Differential expression analyses were performed on each dataset with “limma” algorithm, including macrophages (GSE158792, GSE56755, and GSE75829) and other immune cells (GSE23371, GSE27291, GSE157737, GSE27838, GSE37750, GSE8059, GSE51540, GSE149425, GSE104852, GSE52156, GSE155148, GSE83441, GSE42058, GSE106932, and GSE28726). The platform information of these datasets was shown in Supplementary Table S1. Subsequently, 994 upregulated macrophage genes and 6,057 downregulated genes were considered as statistically significantly different in TIICs. Then, 422 dysregulated genes were identified that were both highly upregulated in macrophages and downregulated in other immune cells, demonstrating their expression specificity for TAMs rather than other types of immune cells. Finally, 404 MSG were screened and combined with clinical features for the construction of a prognostic risk model (Figure 2A).
FIGURE 2
FIGURE 2. Immune subtype analysis of TCGA cohort. (A) Venn diagram plot among TCGA, macrophage, and other immune cell genes. (B) Correlation heatmap for TCGA cancer samples based on the CIBERSORT immune infiltration algorithm. (C–E) Consensus score matrix for cancer samples when k = 5. A higher consensus score between two samples indicates that they are more likely to be assigned to the same cluster in different iterations. (F) A K-M survival curve of the five immune subtypes.
Identification and prognostic assessment of immune subtypes in liver cancer
Immunotyping can be used to mirror the immune status in TME and help identify suitable genes for cancer immunotherapy. The immune landscape based on the CIBERSORT algorithm was applied to investigate the immune subtypes. We compared median survival rates among these subtypes in the TCGA cohorts by setting p < 0.05 as the threshold for screening. The correlation heatmap of the CIBERSORT immune infiltration profile was shown in Figure 2B. We chose k = 5 where immune infiltration appeared to be stably clustered (Figures 2C–E), and obtained five immune subtypes that had a significant difference for OS (Figure 2F, p = 0.010). These results support a distinct immune infiltration feature of five immune subtypes, reflecting clinical outcomes due to immune infiltration disorder and its predictive value in classifying liver cancer patients. The association of MSG and immune subtypes was performed with R package “limma”. The cut-off was set at |logFC| >1.5, p < 0.05, adjust p < 0.05.
Identification and verification of the prognostic capability of the TIMSig
The demographic and clinical characteristics of patients in the two cohorts were similar (Supplementary Table S2). In the training cohort (n = 202), we used univariate Cox proportional hazards regression analysis, LASSO regression analysis (Supplementary Figures S1A,B), and multivariate Cox proportional hazards regression analysis to develop TIMSig. The following formula was used for calculation: risk score = (0.170 × S100A9) + (0.090 × SLC22A15) + (0.143 × TRIM54) + (-0.088 × PPARGIA). Cancer patients were classified into different groups according to respective median risk scores. With the increase in the risk scores, the number of liver cancer patients in the high-risk group increased (Figures 3C,D, Figures 4C,D, and Figures 5C,D). In the training cohort (n = 202), the K–M survival curve demonstrated that patients in the low-risk group had a significantly longer OS compared to those in the high-risk group (Figure 3A, p < 0.001, HR: 3.164, 95% HR CI: 1.908–5.245). The AUC of ROC curves for 1, 2, and 3 years was higher than 0.7 in the training cohort (Figure 3B). Similar trends were found in the test cohort (n = 202), the high-risk group exhibited a shorter survival time in liver cancer patients compared to those in the low-risk group (Figure 4A, p = 0.019, HR: 1.755, 95% HR CI: 1.090–2.825). According to the results of ROC curves, the AUC for 1, 2, and 3 years was 0.672, 0.616, and 0.595 in the test cohort, respectively (Figure 4B). To further verify the prognosis prediction capability of TIMSig, the prognostic performance of the TIMSig was further tested using another completely independent GEO cohort with the GSE14520 dataset (n = 242). As shown in the Figure 5A, the K–M curve for OS showed a statistically significant difference between the two groups (p = 0.037, HR: 1.600, 95% HR CI: 1.025–2.497). Patients tend to have a better prognosis in the low-risk group compared to the high-risk group. In Figure 5B, the AUC of the ROC curve was 0.641 for 3 years, which was significantly higher than other clinical features.
FIGURE 3
FIGURE 3. Survival analysis of the training cohort. (A) K–M curve for the training cohort. Patients in the low-risk group represented a better OS. (B) The 1, 2, and 3 years ROC curves of the risk score in the training cohort. (C) Distribution of risk score of liver cancer patients in the training cohort. (D) Distribution of survival status of liver cancer patients in the training cohort. (E) Heatmap of the TIMGs in the training cohort.
FIGURE 4
FIGURE 4. Survival analysis of the test cohort. (A) K–M curve for the test cohort. Patients in the low-risk group represented a better OS. (B) The 1, 2, and 3 years ROC curves of the risk score in the test cohort. (C) Distribution of risk scores of liver cancer patients in the test cohort. (D) Distribution of survival status of liver cancer patients in the test cohort. (E) Heatmap of the TIMGs in the test cohort.
FIGURE 5
FIGURE 5. Validation of GSE14520 cohort. (A) K–M curve for the GSE14520 cohort. Patients in the low-risk group represented a better OS. (B) Comparison of the 3 years' ROC curve with other clinical characteristics. (C) Distribution of risk score of liver cancer patients in the GSE14520 cohort. (D) Distribution of survival status of liver cancer patients in the GSE14520 cohort.
The TIMSig is an independent prognostic factor in liver cancer patients
To investigate the prognostic power of our risk model, the TIMSig was performed as an independent prognostic factor in univariable and multivariable Cox analyses for each cohort (Table 1). For the training cohort, univariate Cox regression analysis showed that risk score (p < 0.001, HR: 2.717, 95% HR CI: 1.963–3.760) and stage (p = 0.003, HR: 1.443, 95% HR CI: 1.133–1.837) were significantly related to the prognosis of liver cancer. Multivariate Cox regression analysis confirmed that the risk score (p < 0.001, HR: 2.562, 95% HR CI: 1.854–3.540) still was an independent prognostic factor after adjusting for other clinicopathologic factors. Then, the role of TIMSig was validated by using the test cohort, the results of univariate Cox regression analysis (p = 0.006, HR: 1.529, 95% HR CI: 1.132–2.067) and multivariate Cox regression analysis (p = 0.014, HR: 1.493, 95% HR CI: 1.084–2.057) for OS represented that risk score were significantly relevant to the clinical outcomes of liver cancer. For the GSE14520 cohort, the TIMSig still were an independent prognostic predictor for OS in the univariate Cox regression analysis (p = 0.022, HR: 2.179, 95% HR CI: 1.118–4.249) and multivariate Cox regression analysis (p = 0.019, HR: 2.282, 95% HR CI: 1.146–4.545). Taken together, the TIMSig can be applied to predict the prognosis of patients with liver cancer and has a great accuracy, which is consistent across different cohorts. Moreover, patients with a low expression value of each TIMGs had longer survival times than patients in the high-risk group (Figures 6A–D).
TABLE 1
TABLE 1. Univariate Cox analysis and multivariate Cox analysis of OS in each cohort.
FIGURE 6
FIGURE 6. (A–D) K–M curves for S100A9, SLC22A15, TRIM54, and PPARGC1A in the TCGA cohort.
Potential of the TIMSig as an indicator of immunotherapy response in liver cancer
Given the success of anti-PD1, anti-PDL1, and anti-CTLA4 treatments in cancer patients, we compared the difference in the expression value of immune checkpoints (including PD-1, PD-L1, and CTLA-4), and the different subgroups categorized by the median of TIMSig. In the training cohort, patients with high TIMSig tend to show higher expression value of immune checkpoints compared with low TIMSig (Figure 7A, p < 0.001, p = 0.004, p < 0.001 for comparison of PD-1, PD-L1, and CTLA-4 with each group). Also, the patients with the high TIMSig still show higher expression value of immune checkpoints compared with low TIMSig in the test cohort, (Figure 7B, p < 0.001, p = 0.017, p < 0.001 for comparison of PD-1, PD-L1, CTLA-4 with each group).
FIGURE 7
FIGURE 7. Immunotherapy response with TIMSig. (A,B) Comparison of the expression value of immune checkpoints (PD1, PDL1, and CTLA4) with different groups of TIMSig in each cohort. (C–E) The Pearson correlation between immune checkpoints and the risk scores in the TCGA cohort. (F–H) K-M survival curves of OS among four patient groups stratified by the TIMSig and PD1, PDL1, and CTLA4 in the TCGA cohort. (I) Comparison of the risk score with different groups of PD1/PDL1 therapy responses, paired t-test was used as the significance test. (J) The ROC curve to estimate the sensitivity of TIMSig to PD1/PDL1 therapy responses. (K) K–M survival curves of PFS among four patient groups stratified by the TIMSig and PD1/PDL1 therapy reponses.
To verify whether the TIMSig has an impact on OS in liver cancer patients with the homologous trend of immune checkpoints and different risk groups. Based on the TIMSig and the expression value of immune checkpoints, the survival distribution of total TCGA patients was partitioned into four groups and compared by log-rank test. We found that TIMSig was positively related to PD-1 (Figure 7C, R = 0.34, p < 0.001), PD-L1 (Figure 7D, R = 0.23, p < 0.001), CTLA-4 (Figure 7E, R = 0.43, p < 0.001). As shown in Figure 7F, patients with a low level of PD-1 and high TIMSig would have significantly worse OS than the other three groups (p < 0.001), whereas patients with high TIMSig and a low level of PD-L1 were inclined to be the worst clinical outcomes relative to other three groups (Figure 7G, p < 0.001). The same statistical difference with OS was repeated using TIMSig and CTLA-4 (Figure 7H, p < 0.001). It is indicated that patients with liver cancer stratified by TIMSig and PD-L1 or CTLA-4 exhibited OS analogous to PD-1.
The 24 clinical cases of anti-PD1/PDL1 clinical treatment for HCC patients were obtained from GSE140901. A total of 13 PD-1/PD-L1 responses (54.2%) occurred, while 11 cases (45.8%) have nonresponses on PD-1/PD-L1 therapy. Based on the results of different analyses by paired t-test (Figure 7I, p = 0.012), TIMSig is significantly associated with PD-1/PD-L1 therapy response. The AUC value of the model established in the ICI clinical cohort is 0.622 and is a great predicted value for the estimation of clinical immunotherapeutic efficacy (Figure 7J). In addition, TIMSig combined PD-1/PD-L1 response can significantly stratify PFS (Figure 7K, p < 0.0001). These observed associations between the TIMSig and immune checkpoints confirmed our hypothesis that the TIMSig may be a great predictive biomarker for cancer immunotherapy response.
The pertinence of TIMGs and immune infiltration level in liver cancer
ESTIMATE algorithm gave scores for 404 TCGA cancer samples, containing the immune score, ESTIMATE score, and tumor purity. We compared the distribution of these scores with TIMSig to immune characteristics. Pearson correlation tests were applied while each specific score and TIMSig was inspected and analyzed as shown in Figures 8A–C. TIMSig was positively associated with immune score (R = 0.29, p < 0.001) and ESTIMATE score (R = 0.21, p < 0.001), while TIMSig and tumor purity was negatively correlated (R = -0.21, p < 0.001). To evaluate whether TIMGs can accurately predict the distribution of TIICs, the association of TIMGs expression level with TIICs’ abundance was evaluated using the TIMER2.0 database. It suggested that high expression levels of S100A9, SLC22A15, TRIM54, and PPARGC1A were significantly correlated with increased infiltration degree of macrophages, DCs, B cells, and NK cells (Figures 8D–G). Therefore, TIMGs may directly or indirectly present the immune cells in triggering an immune response.
FIGURE 8
FIGURE 8. Tumor immune microenvironment analyses. (A) Correlation between the immune score and risk score in TCGA cancer samples. (B) Correlation between the ESTIMATE score and risk score in TCGA cancer samples. (C) Correlation between the tumor purity and risk score in TCGA cancer samples. (D–G) Identification of TIMGs associated with TIICs.
Functional analyses for TIMGs
The main role of TAMs is to accelerate tumor growth by secreting many cytokines. To analyze the function of macrophages in liver cancer tissues, we analyzed the relationship between TIMSig and cytokine gene expression. As shown in Figure 9A, eight cytokines (IL-1β, IL-6, IL-10, TGF-β, CCL2, CXCL8, CSF-1, and VEGF) were identified to have higher expression values in the high-risk group than in the low-risk group, which was relevant to the progression, invasion, and metastasis of cancers. The GO analysis of the biological process (BP), molecular function (MF), and cell component (CC) showed that most of the enriched terms were related to cellular material transport and material metabolism (Figure 9B). The KEGG pathway enrichment analysis illustrated that the TIMSig were mainly enriched in NOD-like receptor signaling pathway, cytokine–cytokine receptor interaction, and Fc gamma R-mediated phagocytosis, most of which are related to immunity and metabolism (Figure 9C). Concurrently, the protein level of S100A9, SLC22A15, and TRIM54 was found to be much higher in the cancer tissues or the cells around the blood sinus compared to the normal tissues (Supplementary Figure S2).
FIGURE 9
FIGURE 9. Functional analysis for TIMGs. (A) Comparison of the expression value of cytokines (IL-1β, IL-6, IL-10, TGF-β, CCL2, CXCL8, CSF-1, and VEGF) with different groups of TIMSig in the TCGA cohort. (B) GO function analysis of TIMSig in the TCGA cohort. (C) KEGG pathway analysis of TIMSig in the TCGA cohort.
Screening of sensitive chemotherapy drugs
Based on the pRRophetic algorithm, we explored the relationship between TIMSig and drug chemoresistance by calculating the half-maximal inhibitory concentration (IC50) of six common chemotherapeutic drugs (cytarabine, rapamycin, cisplatin, sunitinib, erlotinib, and methotrexate) for liver cancer. In Figures 10A–F, we observed that cancer patients in the high-risk group were more resistant to erlotinib. On the contrary, the patients with a high score of TIMSig were more sensitive to cytarabine, rapamycin, cisplatin, sunitinib, and methotrexate.
FIGURE 10
FIGURE 10. Chemotherapeutic drug sensitivity analysis based on TIMSig. (A–F) Comparison of the IC50 levels (cytarabine, rapamycin, cisplatin, sunitinib, erlotinib, and methotrexate) with different groups of TIMSig in the TCGA cohort.
Discussion
Liver cancer is the leading cause of cancer-related death in the world, due to its characteristics of late diagnosis, poor prognosis, and high heterogeneity (Heinrich et al., 2021; Sung et al., 2021). Hence, identifying reliable and effective biomarkers for liver cancer prognosis is of great importance. The tumor's immune status largely influences the effectiveness of cancer immunotherapy. Accumulating evidence has suggested that distinct molecular subtypes of tumors are positively associated with OS (Zhang et al., 2020; Hu et al., 2021; Liu et al., 2021). In this study, we defined five immune infiltration subtypes and screened TAM-related molecular biomarkers, which could help to predict the clinical response of ICI treatments.
TAMs can promote malignant cell proliferation by interacting with cancer cells by secreting ingredient exosomes and cytokines. Accumulating evidence has indicated that TAMs in TME is associated with the prognosis and immunological characteristics of a variety of cancers (Lam et al., 2019; Li et al., 2019; Xiao N et al., 2021). TAMs can mediate PD-1 drug resistance in HCC through the PD-L1 pathway and regulation of T cells (Pu and Ji, 2022). Molecular biomarkers of TAMs could be used for risk-stratification and response prediction of cancer treatments, though the fundamental cognition of cellular structures and molecular landscapes in liver cancer remains difficult to define. Traditional strategies for clinical models perform investigation mainly at the single component level and have inherent limitations in providing precise prognostic information on complex component cells residing in a highly multicomponent TME. In previous investigations, the signatures for prognostic prediction based on immune-related genes or immune microenvironment have been described in many kinds of cancers, and several mRNAs or lncRNA-associated signatures have also been developed to predict the clinical outcomes of HCC patients (Xiao B et al., 2021; Kang et al., 2021; Pi et al., 2021). However, the potential ability of TAMs for immunotherapy prediction has been ignored. In this study, we used the transcriptome profile of TAM to construct a risk model that could reflect the function of macrophages.
In our results, S100A9, SLC22A15, TRIM54, and PPARGC1A were identified as TIMGs, which were identified as prognostic biomarkers for liver cancer sufferers. The expression of these TIMGs in the training cohort and test cohort was shown in Figure 3E and Figure 4E. In addition, based on the pan-cancer analyses, TIMGs are differentially expressed not only in patients with liver cancer, but also in many other cancers (Supplementary Figure S3). In most cancer types, TIMGs exhibited an upregulated tendency in cancer samples compared to normal samples. To date, studies on S100A9 and PPARGC1A have been suggested as a new potential biomarker for liver cancer. S100A9, a secreted protein related to the inflammatory immune microenvironment and the functional phenotype of macrophages, is significantly increased in TAMs of HCC (Ganta et al., 2019). As a downstream regulator of VEGFR1 macrophage polarization, S100A9 is a promoter of M2 polarization. It has been found that after the knockout of S100A9, the levels and activities of CX3CR1 and Nur77 in macrophages decreased significantly, which resulted in the decrease of efferocytosis in macrophages and the accumulation of necrotic cells in tissues (Marinković et al., 2020; Willers et al., 2020). The existence of S100A9 positive macrophages in tumor tissues, a key gene in the growth and metastasis of HCC, was related to the shorter survival time and the poor treatment of PD-1 antibody in metastatic cancer patients (Duan et al., 2018; Kwak et al., 2020; Wei et al., 2021). PPARGC1A is a mitochondrial regulator, which can regulate mitochondrial biogenesis in macrophages and play a regulatory role in the growth and metastasis of liver cancer. The up-regulation of PPARGC1A will enhance the oxidation of fatty acids in mitochondria to reduce the accumulation of free fatty acids, which leads to abnormal mitochondrial function (Soliman et al., 2020). PPARGC1A mediates YAP to reprogram cell metabolism, shifting substrates from gluconeogenesis to growth anabolism (Hu et al., 2017). At the same time, it effectively inhibits aerobic glycolysis by adjusting Wnt/β-catenin/pdk1 axis, thus inhibiting the migration and invasion of HCC. Meanwhile, PPARGC1A interacts with TNFAIP3 and HSPA12A, produces nuclear translocation, induces AOAH expression, participates in mitochondrial regulation and the expression of inflammatory genes such as NLRP3, and plays a vital role in the homeostasis of telomeres and macrophages’ mitochondria (Kang et al., 2018; Liu et al., 2020).
TAM promotes intravasation, extravasation, and metastasis of tumor cells by secreting various pro-tumoral molecule proteins. Macrophages change their functional state by responding to signaling molecules in the TME, and secrete various cytokines (such as IL-1β, IL-6, CXCL-8, IL-10, CCL2, and CSF-1) to interact with various types of cells in the TME, thus playing a momentous role in the regulation of tumor invasiveness, metastasis and drug resistance (Galdiero et al., 2018; Lin et al., 2019; Korbecki et al., 2020). A common feature of the nine cytokines in Figure 9 is the enhancement of tumor progression via efferocytosis of TAMs which is essential for understanding the immune status of TME and improving the prognosis of liver cancer patients. It is proved that TIMGs are significantly associated with the biological functions of the immune system and will be a promising tool reflecting the function of TAMs in liver cancer. Thus, we consider that TIMSig plays a role in the crosstalk between the macrophages and the TME of liver cancer. For further studies, the connection among TIMGs, cytokines, and corresponding pathways could be an interesting direction.
Many previous studies have identified meaningful immune infiltration signatures for cancer patients (Chen et al., 2021; Lin et al., 2021; Xu et al., 2021). However, the immune subtype or the interaction between TAMs and immune cells in cancer metabolism was rarely considered. Compared with other similar prognostic models based on bioinformatics, our signature passed the test of clinical PD-1/PD-L1 treated cohort and has good relevance with PFS, which is the biggest advantage of our model. In addition, our study not only screened the immune subtypes of liver cancer based on immune infiltration, but also investigate the potential TIMSig-based regulatory mechanisms where the immune system participated and allowed robust risk-stratification, thus enhancing a broader notion of the TME-based prognostic model. Finally, based on the drug resistance algorithm, we explored six chemotherapeutic drugs that were related to TIMSig, which could provide guidelines for clinical cancer treatments. Except for macrophages, our algorithm is also suitable for identifying the specific genes of other immune cells in the prognosis of tumor patients.
There are several limitations to our study. The components of tumor tissues were complicated, it is not certain whether the collected cases and algorithm could accurately reflect the function of TAMs on the survival of liver cancer sufferers. In addition, our prognostic model is only based on bulk RNA-sequencing data and retrospective clinical cohort, and validation of this model by single-cell sequencing or flow cytometry might give a higher potential for its application. Finally, the mechanism of TRIM54 and SLC22A15 affecting liver cancer development through TAMs is still unclear, and further research was warranted.
Conclusion
Overall, S100A9, SLC22A15, TRIM54, and PPARGC1A were screened as TIMGs that can be used for prognostic prediction and be the potential targets of the ICI treatments for patients with liver cancer. The prediction model integrated specific genes of immune cells in the TME will help clinicians not only make rational immunotherapeutic decisions but also understand the driving nodes in the machinery of liver cancer.
Data availability statement
Publicly available datasets were analyzed in this study. These data can be found at: GEO (https://www.ncbi.nlm.nih.gov/geo/database), TCGA (https://portal.gdc.cancer.gov/), CIBERSORT (https://cibersort.stanford.edu/), TIMER2.0 (http://timer.comp-genomics.org/), GSEA (http://www.gsea-msigdb.org/gsea/index.jsp), HPA (https://www.proteinatlas.org/).
Author contributions
JX and ZY were responsible for the study design. AH, BL, JY, and CL developed the methodology. AH and YZ was mainly responsible for data collection and data analysis. AH, BL, and JL was responsible for manuscript draft and review. All authors read and approved the final manuscript.
Funding
This research was supported by the major scientific and technological innovation project of Wenzhou Science and Technology Bureau (grant no. ZY2021009).
Acknowledgments
All authors gratefully acknowledge TCGA, GEO, CIBERSORT, TIMER, GSEA, and HPA databases for providing such excellent platforms.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.983840/full#supplementary-material
References
Anderson, N. M., and Simon, M. C. (2020). The tumor microenvironment. Curr. Biol. 30 (16), R921. Epub 2020/08/19PubMed PMID: 32810447; PubMed Central PMCID: PMCPMC8194051. doi:10.1016/j.cub.2020.06.081
Bagchi, S., Yuan, R., and Engleman, E. G. (2021). Immune checkpoint inhibitors for the treatment of cancer: Clinical impact and mechanisms of response and resistance. Annu. Rev. Pathol. 16, 223–249. Epub 2020/11/17PubMed PMID: 33197221. doi:10.1146/annurev-pathol-042020-042741
Binnewies, M., Roberts, E. W., Kersten, K., Chan, V., Fearon, D. F., Merad, M., et al. (2018). Understanding the tumor immune microenvironment (TIME) for effective therapy. Nat. Med. 24 (5), 541–550. Epub 2018/04/25PubMed PMID: 29686425; PubMed Central PMCID: PMCPMC5998822. doi:10.1038/s41591-018-0014-x
Chen, X., Yan, L., Jiang, F., Lu, Y., Zeng, N., Yang, S., et al. (2021). Identification of a ferroptosis-related signature associated with prognosis and immune infiltration in adrenocortical carcinoma. Int. J. Endocrinol. 2021, 4654302. Epub 2021/08/03PubMed PMID: 34335745; PubMed Central PMCID: PMCPMC8318759. doi:10.1155/2021/4654302
Duan, L., Wu, R., Zhang, X., Wang, D., You, Y., Zhang, Y., et al. (2018). HBx-induced S100A9 in NF-κB dependent manner promotes growth and metastasis of hepatocellular carcinoma cells. Cell Death Dis. 9 (6), 629. PubMed PMID: 29795379. doi:10.1038/s41419-018-0512-2
Dyck, L., and Mills, K. H. G. (2017). Immune checkpoints and their inhibition in cancer and infectious diseases. Eur. J. Immunol. 47 (5), 765–779. Epub 2017/04/11PubMed PMID: 28393361. doi:10.1002/eji.201646875
Fridman, W. H., Zitvogel, L., Sautès-Fridman, C., and Kroemer, G. (2017). The immune contexture in cancer prognosis and treatment. Nat. Rev. Clin. Oncol. 14 (12), 717–734. Epub 2017/07/26PubMed PMID: 28741618. doi:10.1038/nrclinonc.2017.101
Fu, Y., Liu, S., Zeng, S., and Shen, H. (2019). From bench to bed: the tumor immune microenvironment and current immunotherapeutic strategies for hepatocellular carcinoma. J. Exp. Clin. Cancer Res. 38 (1), 396. Epub 2019/09/11PubMed PMID: 31500650; PubMed Central PMCID: PMCPMC6734524. doi:10.1186/s13046-019-1396-4
Galdiero, M. R., Marone, G., and Mantovani, A. (2018). Cancer inflammation and cytokines. Cold Spring Harb. Perspect. Biol. 10 (8), a028662. Epub 2017/08/06PubMed PMID: 28778871; PubMed Central PMCID: PMCPMC6071493. doi:10.1101/cshperspect.a028662
Ganta, V. C., Choi, M., Farber, C. R., and Annex, B. H. (2019). Antiangiogenic VEGF165b regulates macrophage polarization via S100A8/S100A9 in peripheral artery disease. Circulation 139 (2), 226–242. PubMed PMID: 30586702. doi:10.1161/CIRCULATIONAHA.118.034165
Giraldo, N. A., Sanchez-Salas, R., Peske, J. D., Vano, Y., Becht, E., Petitprez, F., et al. (2019). The clinical role of the TME in solid cancer. Br. J. Cancer 120 (1), 45–53. Epub 2018/11/11PubMed PMID: 30413828; PubMed Central PMCID: PMCPMC6325164. doi:10.1038/s41416-018-0327-z
Heinrich, S., Craig, A. J., Ma, L., Heinrich, B., Greten, T. F., and Wang, X. W. (2021). Understanding tumour cell heterogeneity and its implication for immunotherapy in liver cancer using single-cell analysis. J. Hepatol. 74 (3), 700–715. Epub 2020/12/04PubMed PMID: 33271159. doi:10.1016/j.jhep.2020.11.036
Hu, Y., Shin, D-J., Pan, H., Lin, Z., Dreyfuss, J. M., Camargo, F. D., et al. (2017). YAP suppresses gluconeogenic gene expression through PGC1α. Hepatol. Baltim. Md) 66 (6), 2029–2041. PubMed PMID: 28714135. doi:10.1002/hep.29373
Hu, Y., Zheng, M., Wang, S., Gao, L., Gou, R., Liu, O., et al. (2021). Identification of a five-gene signature of the RGS gene family with prognostic value in ovarian cancer. Genomics 113 (4), 2134–2144. Epub 2021/04/13PubMed PMID: 33845140. doi:10.1016/j.ygeno.2021.04.012
Johnstone, S. E., Reyes, A., Qi, Y., Adriaens, C., Hegazi, E., Pelka, K., et al. (2020). Large-scale topological changes restrain malignant progression in colorectal cancer. Cell 182 (6), 1474–1489. Epub 2020/08/26PubMed PMID: 32841603; PubMed Central PMCID: PMCPMC7575124. doi:10.1016/j.cell.2020.07.030
Kang, Y., Zhang, H., Zhao, Y., Wang, Y., Wang, W., He, Y., et al. (2018). Telomere dysfunction disturbs macrophage mitochondrial metabolism and the NLRP3 inflammasome through the PGC-1α/TNFAIP3 Axis. Cell Rep. 22 (13), 3493–3506. PubMed PMID: 29590618. doi:10.1016/j.celrep.2018.02.071
Kang, Z., Li, W., Yu, Y. H., Che, M., Yang, M. L., Len, J. J., et al. (2021). Identification of immune-related genes associated with bladder cancer based on immunological characteristics and their correlation with the prognosis. Front. Genet. 12, 763590. Epub 2021/12/14PubMed PMID: 34899848; PubMed Central PMCID: PMCPMC8664377. doi:10.3389/fgene.2021.763590
Korbecki, J., Grochans, S., Gutowska, I., Barczak, K., and Baranowska-Bosiacka, I. (2020). CC chemokines in a tumor: A review of pro-cancer and anti-cancer properties of receptors CCR5, CCR6, CCR7, CCR8, CCR9, and CCR10 ligands. Int. J. Mol. Sci. 21 (20), E7619. Epub 2020/10/21PubMed PMID: 33076281; PubMed Central PMCID: PMCPMC7590012. doi:10.3390/ijms21207619
Kurebayashi, Y., Ojima, H., Tsujikawa, H., Kubota, N., Maehara, J., Abe, Y., et al. (2018). Landscape of immune microenvironment in hepatocellular carcinoma and its additional impact on histological and molecular classification. Hepatology 68 (3), 1025–1041. Epub 2018/04/01PubMed PMID: 29603348. doi:10.1002/hep.29904
Kwak, T., Wang, F., Deng, H., Condamine, T., Kumar, V., Perego, M., et al. (2020). Distinct populations of immune-suppressive macrophages differentiate from monocytic myeloid-derived suppressor cells in cancer. Cell Rep. 33 (13), 108571. PubMed PMID: 33378668. doi:10.1016/j.celrep.2020.108571
Lam, J. H., Ng, H. H. M., Lim, C. J., Sim, X. N., Malavasi, F., Li, H., et al. (2019). Expression of CD38 on macrophages predicts improved prognosis in hepatocellular carcinoma. Front. Immunol. 10, 2093. Epub 2019/09/26PubMed PMID: 31552039; PubMed Central PMCID: PMCPMC6738266. doi:10.3389/fimmu.2019.02093
Lei, X., Lei, Y., Li, J. K., Du, W. X., Li, R. G., Yang, J., et al. (2020). Immune cells within the tumor microenvironment: Biological functions and roles in cancer immunotherapy. Cancer Lett. 470, 126–133. Epub 2019/11/16 PubMed PMID: 31730903. doi:10.1016/j.canlet.2019.11.009
Li, X., Liu, R., Su, X., Pan, Y., Han, X., Shao, C., et al. (2019). Harnessing tumor-associated macrophages as aids for cancer immunotherapy. Mol. Cancer 18 (1), 177. Epub 2019/12/07PubMed PMID: 31805946; PubMed Central PMCID: PMCPMC6894344. doi:10.1186/s12943-019-1102-3
Lin, Y., Xu, J., and Lan, H. (2019). Tumor-associated macrophages in tumor metastasis: biological roles and clinical therapeutic applications. J. Hematol. Oncol. 12 (1), 76. Epub 2019/07/14PubMed PMID: 31300030; PubMed Central PMCID: PMCPMC6626377. doi:10.1186/s13045-019-0760-3
Lin, R., Fogarty, C. E., Ma, B., Li, H., Ni, G., Liu, X., et al. (2021). Identification of ferroptosis genes in immune infiltration and prognosis in thyroid papillary carcinoma using network analysis. BMC Genomics 22 (1), 576. Epub 2021/07/29PubMed PMID: 34315405; PubMed Central PMCID: PMCPMC8314640. doi:10.1186/s12864-021-07895-6
Liu, J., Du, S., Kong, Q., Zhang, X., Jiang, S., Cao, X., et al. (2020). HSPA12A attenuates lipopolysaccharide-induced liver injury through inhibiting caspase-11-mediated hepatocyte pyroptosis via PGC-1α-dependent acyloxyacyl hydrolase expression. Cell Death Differ. 27 (9), 2651–2667. PubMed PMID: 32332915. doi:10.1038/s41418-020-0536-x
Liu, C., Papukashvili, D., Dong, Y., Wang, X., Hu, X., Yang, N., et al. (2021). Identification of tumor antigens and design of mRNA vaccine for colorectal cancer based on the immune subtype. Front. Cell Dev. Biol. 9, 783527. Epub 2022/02/08PubMed PMID: 35127707; PubMed Central PMCID: PMCPMC8811447. doi:10.3389/fcell.2021.783527
Lujambio, A., Akkari, L., Simon, J., Grace, D., Tschaharganeh, D. F., Bolden, J. E., et al. (2013). Non-cell-autonomous tumor suppression by p53. Cell 153 (2), 449–460. Epub 2013/04/09PubMed PMID: 23562644; PubMed Central PMCID: PMCPMC3702034. doi:10.1016/j.cell.2013.03.020
Mantovani, A., Marchesi, F., Malesci, A., Laghi, L., and Allavena, P. (2017). Tumour-associated macrophages as treatment targets in oncology. Nat. Rev. Clin. Oncol. 14 (7), 399–416. Epub 2017/01/25PubMed PMID: 28117416; PubMed Central PMCID: PMCPMC5480600. doi:10.1038/nrclinonc.2016.217
Marinković, G., Koenis, D. S., de Camp, L., Jablonowski, R., Graber, N., de Waard, V., et al. (2020). S100A9 links inflammation and repair in myocardial infarction. Circ. Res. 127 (5), 664–676. PubMed PMID: 32434457. doi:10.1161/CIRCRESAHA.120.315865
Martínez, V. G., Rubio, C., Martínez-Fernández, M., Segovia, C., López-Calderón, F., Garín, M. I., et al. (2017). BMP4 induces M2 macrophage polarization and favors tumor progression in bladder cancer. Clin. Cancer Res. 23 (23), 7388–7399. Epub 2017/09/21PubMed PMID: 28928159. doi:10.1158/1078-0432.Ccr-17-1004
Perez-Ruiz, E., Minute, L., Otano, I., Alvarez, M., Ochoa, M. C., Belsue, V., et al. (2019). Prophylactic TNF blockade uncouples efficacy and toxicity in dual CTLA-4 and PD-1 immunotherapy. Nature 569 (7756), 428–432. Epub 2019/05/03PubMed PMID: 31043740. doi:10.1038/s41586-019-1162-y
Pi, Y. N., Guo, J. N., Lou, G., and Cui, B. B. (2021). Comprehensive analysis of prognostic immune-related genes and drug sensitivity in cervical cancer. Cancer Cell Int. 21 (1), 639. Epub 2021/12/03PubMed PMID: 34852825; PubMed Central PMCID: PMCPMC8638517. doi:10.1186/s12935-021-02333-9
Pu, Y., and Ji, Q. (2022). Tumor-associated macrophages regulate PD-1/PD-L1 immunosuppression. Front. Immunol. 13, 874589. Epub 2022/05/21PubMed PMID: 35592338; PubMed Central PMCID: PMCPMC9110638. doi:10.3389/fimmu.2022.874589
Qian, B. Z., and Pollard, J. W. (2010). Macrophage diversity enhances tumor progression and metastasis. Cell 141 (1), 39–51. Epub 2010/04/08PubMed PMID: 20371344; PubMed Central PMCID: PMCPMC4994190. doi:10.1016/j.cell.2010.03.014
Sacks, D., Baxter, B., Campbell, B. C. V., Carpenter, J. S., Cognard, C., Dippel, D., et al. (2018). Multisociety consensus quality improvement revised consensus statement for endovascular therapy of acute ischemic stroke. Int. J. Stroke 13 (6), 612–632. Epub 2018/05/23PubMed PMID: 29786478. doi:10.1177/1747493018778713
Schoenfeld, A. J., and Hellmann, M. D. (2020). Acquired resistance to immune checkpoint inhibitors. Cancer Cell 37 (4), 443–455. Epub 2020/04/15PubMed PMID: 32289269; PubMed Central PMCID: PMCPMC7182070. doi:10.1016/j.ccell.2020.03.017
Soliman, E., Elhassanny, A. E. M., Malur, A., McPeek, M., Bell, A., Leffler, N., et al. (2020). Impaired mitochondrial function of alveolar macrophages in carbon nanotube-induced chronic pulmonary granulomatous disease. Toxicology 445, 152598. PubMed PMID: 32976959. doi:10.1016/j.tox.2020.152598
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Ca. Cancer J. Clin. 71 (3), 209–249. Epub 2021/02/05PubMed PMID: 33538338. doi:10.3322/caac.21660
Syn, N. L., Teng, M. W. L., Mok, T. S. K., and Soo, R. A. (2017). De-novo and acquired resistance to immune checkpoint targeting. Lancet. Oncol. 18 (12), e731–e741. Epub 2017/12/07PubMed PMID: 29208439. doi:10.1016/s1470-2045(17)30607-1
Vaddepally, R. K., Kharel, P., Pandey, R., Garje, R., and Chandra, A. B. (2020). Review of indications of FDA-approved immune checkpoint inhibitors per NCCN guidelines with the level of evidence. Cancers (Basel) 12 (3), E738. Epub 2020/04/05PubMed PMID: 32245016; PubMed Central PMCID: PMCPMC7140028. doi:10.3390/cancers12030738
Väyrynen, J. P., Haruki, K., Lau, M. C., Väyrynen, S. A., Zhong, R., Dias Costa, A., et al. (2021). The prognostic role of macrophage polarization in the colorectal cancer microenvironment. Cancer Immunol. Res. 9 (1), 8–19. Epub 2020/10/08PubMed PMID: 33023967; PubMed Central PMCID: PMCPMC7785652. doi:10.1158/2326-6066.Cir-20-0527
Vitale, I., Manic, G., Coussens, L. M., Kroemer, G., and Galluzzi, L. (2019). Macrophages and metabolism in the tumor microenvironment. Cell Metab. 30 (1), 36–50. Epub 2019/07/04PubMed PMID: 31269428. doi:10.1016/j.cmet.2019.06.001
Wei, S. C., Duffy, C. R., and Allison, J. P. (2018). Fundamental mechanisms of immune checkpoint blockade therapy. Cancer Discov. 8 (9), 1069–1086. Epub 2018/08/18PubMed PMID: 30115704. doi:10.1158/2159-8290.Cd-18-0367
Wei, R., Zhu, W-W., Yu, G-Y., Wang, X., Gao, C., Zhou, X., et al. (2021). S100 calcium-binding protein A9 from tumor-associated macrophage enhances cancer stem cell-like properties of hepatocellular carcinoma. Int. J. Cancer 148 (5), 1233–1244. PubMed PMID: 33205453. doi:10.1002/ijc.33371
Willers, M., Ulas, T., Vãllger, L., Vogl, T., Heinemann, A. S., Pirr, S., et al. (2020). S100A8 and S100A9 are important for postnatal development of gut microbiota and immune system in mice and infants. Gastroenterology 159 (6), 2130–2145. PubMed PMID: 32805279. doi:10.1053/j.gastro.2020.08.019
Wong, W. K., Yin, B., Lam, C. Y. K., Huang, Y., Yan, J., Tan, Z., et al. (2021). The interplay between epigenetic regulation and CD8(+) T cell differentiation/exhaustion for T cell immunotherapy. Front. Cell Dev. Biol. 9, 783227. Epub 2022/01/29PubMed PMID: 35087832; PubMed Central PMCID: PMCPMC8787221. doi:10.3389/fcell.2021.783227
Wynn, T. A., Chawla, A., and Pollard, J. W. (2013). Macrophage biology in development, homeostasis and disease. Nature 496 (7446), 445–455. Epub 2013/04/27PubMed PMID: 23619691; PubMed Central PMCID: PMCPMC3725458. doi:10.1038/nature12034
Xiao, N., Li, K., Zhu, X., Xu, B., Liu, X., Lei, M., et al. (2021). CD74(+) macrophages are associated with favorable prognosis and immune contexture in hepatocellular carcinoma. Cancer Immunol. Immunother. 71, 57. Epub 2021/05/20PubMed PMID: 34009409. doi:10.1007/s00262-021-02962-z
Xiao, B., Liu, L., Li, A., Wang, P., Xiang, C., Li, H., et al. (2021). Identification and validation of immune-related lncRNA prognostic signatures for melanoma. Immun. Inflamm. Dis. 9 (3), 1044–1054. Epub 2021/06/03PubMed PMID: 34077998; PubMed Central PMCID: PMCPMC8342236. doi:10.1002/iid3.468
Xu, Q., Chen, S., Hu, Y., and Huang, W. (2021). Landscape of immune microenvironment under immune cell infiltration pattern in breast cancer. Front. Immunol. 12, 711433. Epub 2021/09/14PubMed PMID: 34512634; PubMed Central PMCID: PMCPMC8429934. doi:10.3389/fimmu.2021.711433
Yoon, J., Um, H. N., Jang, J., Bae, Y. A., Park, W. J., Kim, H. J., et al. (2019). Eosinophil activation by toll-like receptor 4 ligands regulates macrophage polarization. Front. Cell Dev. Biol. 7, 329. Epub 2020/01/11PubMed PMID: 31921842; PubMed Central PMCID: PMCPMC6933835. doi:10.3389/fcell.2019.00329
Keywords: macrophage, tumor microenvironment, immune infiltration, prognosis, immunotherapy response
Citation: Huang A, Lv B, Zhang Y, Yang J, Li J, Li C, Yu Z and Xia J (2022) Construction of a tumor immune infiltration macrophage signature for predicting prognosis and immunotherapy response in liver cancer. Front. Mol. Biosci. 9:983840. doi: 10.3389/fmolb.2022.983840
Received: 01 July 2022; Accepted: 03 August 2022;
Published: 02 September 2022.
Edited by:
Zhe-Sheng Chen, St. John’s University, United StatesReviewed by:
Siu Hong Dexter Wong, Hong Kong Polytechnic University, Hong Kong SAR, ChinaLishun Wang, Fudan University, China
Copyright © 2022 Huang, Lv, Zhang, Yang, Li, Li, Yu and Xia. 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: Zhijie Yu, emhpamllX3l1QHd6aG9zcGl0YWwuY24=; Jinglin Xia, eGlhamluZ2xpbkB3emhvc3BpdGFsLmNu
†These authors have contributed equally to this work