Skip to main content

CD97 serves as a novel biomarker of immune cell infiltration in hepatocellular carcinoma

Abstract

Background

CD97 is the most widely expressed G protein-coupled receptor in the epidermal growth factor seven-span transmembrane family. It plays a vital role in cell adhesion, migration, and cell connection regulation. We explored the role of CD97 in hepatocellular carcinoma (HCC).

Methods

We evaluated CD97 mRNA expression in HCC using TNMplot and the Gene Expression Omnibus database. The clinical prognostic significance of CD97 in HCC patients was evaluated by gene expression profiling interactive analysis, the Kaplan–Meier plotter, and the UALCAN database. The Tumor Immune Estimation Resource (TIMER) and CIBERSORT databases were used to analyze the relationships among CD97, genes positively related with CD97, and tumor-infiltrating immune cells.

Results

CD97 was highly expressed in HCC tissues and was associated with an adverse prognosis. CD97 and genes positively related with CD97 were positively correlated with the abundance of tumor-infiltrating immune cells and strongly correlated with tumor-infiltrating macrophages (all r ≥ 0.513, P < 0.001). CD97 was positively correlated with M2 macrophage and tumor-associated macrophage markers (both r ≥ 0.464, P < 0.001). CD97 was found to be an immune-related gene in HCC and positively correlated with the TOX, PD-L1, PD-L2, CTLA4, and PD-1 immune checkpoint genes. CD97 copy number alterations affect the level of immune cell infiltration and mRNA expression.

Conclusions

CD97 can be used as a potential molecular marker of prognosis in HCC, which is associated with immune cell infiltration.

Background

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors worldwide. The incidence and mortality rate of HCC rank sixth and third, respectively [1], posing a severe threat to human health. HCC is increasing in incidence in China and is the primary causative factor in chronic hepatitis B virus infection [2]. At present, hepatitis B accounts for the largest HCC burden globally [3]. Recently, with in-depth tumor research and rapid developments in science and technology, the early detection rate and treatment efficacy of HCC have improved [4]. With extensive studies, an increasing number of molecular mechanisms that promote HCC are being discovered, including the inactivation of tumor suppressor genes (such as p53), abnormal activation of oncogenes (such as K-ras), multiple signaling pathways (PI3K, MAPK, JAK/STAT, etc.), and abnormal regulation of epigenetic factors (such as microRNAs) are all involved in HCC development and progression [5]. However, HCC patients still have an extremely shorter survival. Thus, it is very important to explore the novel diagnosis tools, therapy strategy, and prognosis biomarkers for HCC.

Adhesion G protein-coupled receptors (GPCRs) are seven-span transmembrane receptors widely expressed in various cells that play vital roles in cell adhesion and signal transduction [6]. The epidermal growth factor seven-span transmembrane family comprises adherent GPCRs with multiple functions, constituting the largest family of cell surface proteins [7]. CD97/ADGRE5 belongs to class B GPCRs, consisting of an extracellular α-subunit and a transmembrane β-subunit [8]. CD97 is the most widely expressed member of the epidermal growth factor seven-span transmembrane family and is found on the surface of lymphocytes, monocytes, macrophages, dendritic cells, granulocytes, and smooth muscle [9]. It plays a role in cell adhesion, migration, and regulation of intercellular junctions [9]. Previous studies have investigated the relationship between CD97 expression and tumors. Ward et al. found that CD97 and lysophosphatidic acid receptor 1 heterodimerize and function synergistically, mediate Gα12/13 to increase RHO-GTP levels in thyroid and prostate cancers, and promote tumor invasion [10, 11]. Aust et al. found high expression of CD97 in gastric, pancreatic, and esophageal carcinomas [12]. Staining of CD97 was stronger in aggressive tumor cells in > 50% of gastric cancer (GC) patients, and a similar expression pattern was also observed in colorectal cancer. The high expression of CD97 was associated with a worse clinical stage and lymphatic invasion [13]. In a study on rectal adenocarcinoma, the expression of CD97 was significantly higher in patients with recurrence or metastasis than in patients without recurrence or metastasis [14]. Chidambaram et al. found that CD97 increased the migration ability of tumor cells in glioblastoma, but with no effect on proliferation, and the survival of patients with CD97 overexpression was significantly shorter than that in patients with lower expression [15]. Liu et al. found that CD97 was related to the invasion depth and TNM stage of human GC and that CD97 promoted GC cell proliferation and invasion via the MAPK signaling pathway, mediated by exosomes [16, 17]. Yin et al. found that the synergy between CD97 and GPCR kinase 6 (GRK6) promoted epithelial–mesenchymal transition by regulating the expression of downstream MMP2/MMP9, leading to metastasis of HCC [18]. The function of CD97 in tumors remains unknown. Although many studies have shown that CD97 plays a key role in tumor invasion, metastasis, and angiogenesis, the specific regulatory mechanisms remain unclear. With this background, it is especially appealing to explore the role of CD97 in the pathogenesis and development of tumors.

Here, we found that CD97 was correlated with immune cell infiltration and can be used as a prognostic marker in HCC. First, we revealed that CD97 and genes co-expressed with CD97 were abnormally expressed and related to prognosis in HCC. Second, we analyzed the signaling pathways that CD97 is involved in. Furthermore, we showed the relationships among CD97, tumor-infiltrating immune cells (TIICs), and immune-related genes. This study aims to provide a new perspective on the underlying mechanisms of hepatocellular carcinogenesis and to help identify new potential targets for diagnosing and treating HCC.

Materials and methods

CD97 gene expression analysis

We used the pan-cancer analysis in the TNMplot database to evaluate the mRNA levels of CD97 in various tumors. The TNMplot database [19] is a collection of 56,938 high-quality tissue samples, including RNA-seq data from the Genotype-Tissue Expression (GTEx), The Cancer Genome Atlas (TCGA), and the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) databases, as well as GeneChip data from the Gene Expression Omnibus (GEO) database. Next, we selected two datasets (GSE14323 [20] and GSE64041 [21]) containing HCC and normal hepatic tissue gene expression data from the GEO database [22]. Differences in CD97 mRNA expression between HCC and normal hepatic tissues were analyzed, and plots were created using GraphPad Prism (version 8.0.0) software. We also used the Human Protein Atlas (HPA) database to analyze differences in CD97 protein expression between HCC and normal hepatic tissues.

Immunohistochemistry (IHC) staining and scoring

For IHC analysis of CD97 protein expression, the study enrolled 20 patients with HCC from January 2021 to December 2021 at The Affiliated Calmette Hospital of Kunming Medical University (Yunnan, China). For IHC staining, thick serial sections of paraffin-embedded tumor fragment from surgically resected HCC were used to prepare the slides. Antigen retrieval was performed with citrate antigen retrieval solution, and blocked with the endogenous peroxidase enzyme blocking buffer and normal goat serum, CD97 rabbit monoclonal antibody (#GTX108192; GeneTex) at a dilution of 1:200 was incubated with the slides. The slides were observed and photographed using an inverted microscope. The IHC staining score of CD97 in each sample was as follows: one slide was evaluated for each sample, and five regions were randomly selected for each slide for scoring and the average is obtained; the staining extent score was on a scale of 0–4, corresponding to the percentage of stained tumor cells (0%, 1–5%, 6–25%, 26–75%, and 76–100%, respectively), while staining intensity was scored as negative (score = 0), weak (score = 1), moderate (score = 2), and strong (score = 3); the final IHC staining score of CD97 was calculated by multiplying the staining extent score with the intensity score. IHC staining was assessed by two independent pathologists, and a plot was created using GraphPad Prism (version 8.0.0) software.

Correlation analysis between CD97 and clinicopathological features

We downloaded the clinicopathological data of primary HCC patients from the UCSC Xena browser [23]. The HCC cohort in TCGA contains 115 clinical parameters. We selected age, sex, tumor grade, invasion depth, lymph node metastasis, distant metastasis, TNM stage, and survival status as the clinicopathological parameters. After excluding the HCC patients with missing clinicopathological data, we classified the remaining 365 patients into two groups according to the median expression level of CD97: high and low CD97 expression groups. The chi-square test was used to analyze the correlations between CD97 expression and clinicopathological parameters (SPSS version 25.0 software).

Prognostic analysis of CD97 expression in HCC

The Kaplan–Meier plotter (KM plotter) online database [24] was used to evaluate the prognostic value of specific genes in HCC patients (n = 364). UALCAN is a web analysis tool [25] of TCGA that is also used to analyze the prognostic impact of individual genes. We used the KM plotter and UALCAN databases to analyze the relationship between CD97 expression and prognosis in HCC. Furthermore, we used a KM plotter to analyze the clinicopathological features of patients with HCC, and these were mapped using the “forestplot” package (version 2.0.1).

Significant prognostic marker analysis

We downloaded prognostic data of patients with primary HCC from the UCSC Xena browser. We used SPSS (version 25.0) software to determine the independent prognostic factors by univariate and multivariate Cox regression analyses. The results are presented as hazard ratios (HRs), 95% confidence intervals (CIs), and P-values.

Gene set enrichment analysis

We downloaded processed data of patients with HCC from the UCSC Xena browser and divided the patients into high and low CD97 expression groups according to the median expression of CD97. Gene set enrichment analysis (GSEA) (version 4.1.0) software was used to identify the enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in both groups. In the GSEA analysis, the number of genes parameters was 1000, and the nominal P-value, false discovery rate, and normalized enrichment score were used to evaluate the enriched pathways.

Profiling of genes positively related with CD97

To investigate the molecular mechanism of CD97, we used the UALCAN web analysis tool to screen genes positively related with CD97 and select CD97-associated proteins. We selected the top six genes most positively correlated with CD97 for further in-depth research. The relationship between these genes was verified using the Tumor Immune Estimation Resource (TIMER) database [26]. Moreover, we used TNMplot and UALCAN web analysis tools to explore expression differences among genes positively related with CD97. Finally, we used the KM plotter and UALCAN to analyze the prognostic value of the genes positively related with CD97.

Correlation analysis between CD97 or CD97 co-expressed gene expression and immune infiltrating cells

We used the TIMER database to analyze TIICs, including CD8+ T cells, CD4+ T cells, B cells, macrophages, neutrophils, and dendritic cells. We also explored the relationships of the expression of CD97 and its co-expressed genes (IQGAP1, GMIP, MOBKL2A, ARPC2, CD68, and FAM102B) with TIICs. Furthermore, we used the CIBERSORT website [27] to identify the LM22 gene in the HCC samples (including 22 immune cell types) to explore the proportion of immune cells in the HCC high and low CD97 expression groups further. Next, we used the “vioplot” package (version 0.3.7) to create the figures. The TIMER database was used to analyze the relationships between CD97 and gene markers of M1 macrophages, M2 macrophages, and tumor-associated macrophages (TAMs) in HCC.

Correlation analysis between CD97 and tumor immunity-related parameters

We evaluated changes in the tumor microenvironment (TME) in HCC based on the CIBERSORT analysis results and used the “ggplot2” package (version 3.3.3), “heatmap” package (version 1.0.12), and “corrplot” package (version 0.92) in the RStudio (version 3.5.2) to generate the corresponding images. Next, we used the TIMER database to analyze the relationship between the expression of CD97 and that of five immune checkpoints (TOX, PD-L1, PD-L2, CTLA4, and PD-1). We also analyzed the relationship between CD97 copy number changes and the level of immune cell infiltration in HCC. Finally, we used the cBioportal database [28] to analyze the relationship between CD97 copy number alternations and mRNA level changes.

Data analysis

We used RStudio (version 3.5.2) and SPSS (version 25.0) for the statistical analyses, and the values are presented as means ± standard deviation. We used independent and paired sample t-tests to analyze differences in CD97 mRNA expression between HCC and normal hepatic tissues based on the TNMplot and GEO databases. The Pearson chi-squared test was used to analyze the association between CD97 expression and characteristic clinical variables. We then used the KM plotter and UALCAN database to generate survival curves. Univariate and multivariate analyses were performed using Cox proportional hazards regression. Statistically significant differences were defined by P < 0.05.

Results

Elevated expression of CD97 in HCC

We evaluated CD97 mRNA levels in different tumors by analyzing RNA-seq and GeneChip data in the TNMplot database (Fig. 1a), and these results indicate that CD97 is abnormally expressed in various tumor tissues. We showed that CD97 mRNA expression was significantly higher in HCC tissue than in normal hepatic tissue (Fig. 1b). Next, we evaluated CD97 mRNA expression in HCC and normal hepatic tissues using the GEO database and found higher expression in HCC tissue than in normal hepatic tissue (Fig. 1c, d). Moreover, IHC data from the HPA database showed a significantly higher CD97 protein level in HCC than normal hepatic tissues (P = 0.0497; Fig. 2a). Similarly, in our study, the IHC staining score of CD97 in HCC tissues was significantly higher than adjacent normal hepatic tissues (P < 0.001; Fig. 2b). These results suggest higher CD97 expression in HCC than normal hepatic tissues.

Fig. 1
figure 1

CD97 mRNA expression in HCC. a CD97 mRNA levels in different tumors analyzed in TNMplot. b CD97 mRNA levels in HCC tissues and normal hepatic tissues in TNMplot. CD97 mRNA levels in HCC tissues and normal hepatic tissues in the GSE14323 (c) and GSE64041 (d) datasets. *P < 0.05, **P < 0.01, and ***P < 0.001

Fig. 2
figure 2

CD97 protein levels in HCC. Based on the HPA database (a) and our study (b), IHC staining analysis of CD97 protein levels in patients with HCC. Corresponding images and statistics of CD97 levels in adjacent normal hepatic tissues and HCC tissues are shown. *P < 0.05 and ***P < 0.001

Relationship between CD97 expression and clinicopathological features in HCC patients

We assessed the relationship between CD97 mRNA expression and clinicopathological features of HCC patients by using HCC data from TCGA. Based on the median CD97 mRNA level, we assigned 183 HCC patients to the CD97 high expression group and 182 patients to the CD97 low expression group. As shown in Table 1, there were no significant differences in clinicopathological characteristics, including survival status (P = 0.052), invasion depth (P = 0.099), and TNM stage (P = 0.099), between the two groups.

Table 1 Relationship between CD97 expression level and clinicopathological variables in HCC patients

The prognostic value of CD97 expression in HCC

First, we analyzed the association between CD97 expression and the prognosis of HCC patients. The prognosis of HCC patients with high CD97 expression was significantly worse in the UALCAN (P < 0.001; Fig. S1a) and KM plotter (HR = 1.8, P < 0.001; Fig. S1b) databases compared with the patients with low expression. Kaplan–Meier survival analysis showed that male sex (HR = 2.37, P < 0.001), Asian ethnicity (HR = 3, P < 0.001), alcohol consumption (HR = 1.96, P = 0.036), no alcohol consumption (HR = 1.89, P = 0.0068), and no hepatitis viral infection (HR = 1.86, P = 0.0068) were associated with a poor prognosis in the high CD97 expression group (Fig. S1c). The prognosis of HCC patients with high CD97 expression was significantly worse at TNM stages III and IV, pathological grades II and III, and American Joint Committee on Cancer stage III (Fig. S1c). These results indicate that the prognosis of HCC patients with high CD97 expression is correlated with clinicopathological features, especially for those with advanced HCC. Furthermore, Cox univariate survival analysis showed that grade (P < 0.001), invasion depth (P < 0.001), distant metastasis (P = 0.026), TNM stage (P < 0.001), and CD97 expression (P = 0.048) were significantly associated with the prognosis of HCC patients (Table 2). Cox multivariate survival analysis revealed that TNM stage and CD97 expression were independent risk factors for adverse prognosis in HCC (all P < 0.05; Table 3). In summary, our results suggest that CD97 expression predicts an adverse outcome and acts as an independent prognostic marker in HCC.

Table 2 Univariate analysis of the prognostic factors in HCC patients using a Cox regression model
Table 3 Multivariate analysis of the prognostic factors in HCC patients using a Cox regression model

CD97-related signaling pathways in HCC identified by GSEA

We explored the molecular mechanism of CD97 in the development of HCC using GSEA and found 131 KEGG signaling pathways associated with the high CD97 expression group, of which 19 were significantly enriched in the high CD97 expression group (nominal P < 0.05, false discovery rate < 0.1, and normalized enrichment score > 1.6; Table S1). These signaling pathways were positively correlated with immunity and inflammation, including “cytokine receptor interaction,” “T cell receptor signaling pathway,” “B cell receptor signaling pathway,” “graft versus host disease,” “primary immunodeficiency,” “NOD-like receptor signaling pathway,” “intestinal immune network for IgA production,” “allograft,” “chemokine signaling pathway,” “toll-like receptor signaling pathway,” and “FC gamma R mediated phagocytosis,” and to cell adhesion, including “cell adhesion molecules CAMs” (Fig. 3).

Fig. 3
figure 3

GSEA pathways enriched in HCC samples with high CD97 expression. The GSEA results showed that the terms “cytokine receptor interaction” (a), “T cell receptor signaling pathway” (b), “B cell receptor signaling pathway” (c), “graft versus host disease” (d), “primary immunodeficiency” (e), “NOD-like receptor signaling pathway” (f), “intestinal immune network for IgA production” (g), “allograft” (h), “chemokine signaling pathway” (i), “toll-like receptor signaling pathway” (j), “FC gamma R mediated phagocytosis” (k), and “cell adhesion molecules CAMs” (l) were enriched in HCC samples with high CD97 expression. NES, normalized enrichment score

Co-expression analysis of CD97 in HCC

We used the UALCAN analysis tool to mine for genes co-expressed with CD97 to further evaluate the effect of CD97 on HCC development. We selected the top six genes co-expressed with CD97 (r ≥ 0.66) (Fig. 4). Next, we used the TIMER database to verify the co-expression relationships between CD97 and these genes. The results showed that CD97 expression was positively related with IQGAP1 (r = 0.679, P < 0.001), GMIP (r = 0.673, P < 0.001), MOBKL2A (r = 0.642, P < 0.001), ARPC2 (r = 0.607, P < 0.001), CD68 (r = 0.568, P < 0.001), and FAM102B (r = 0.558, P < 0.001) (Fig. 5g–l). UALCAN analysis showed that only CD68 expression was decreased in HCC, while the other five genes had high expression (Fig. S2). The TNMplot database (Fig. S3) mirrored these findings. We also explored the correlations of the six genes with the prognosis of HCC patients using the UALCAN database and found that HCC patients with high GMIP, MOBKL2A, ARPC2, CD68, or FAM102B expression had a significantly worse prognosis (Fig. S4). The KM plotter database showed the same results: patients with high GMIP (HR = 1.46, P = 0.031), MOBKL2A (HR = 1.49, P = 0.022), ARPC2 (HR = 2.03, P = 0.002), CD68 (HR = 1.44, P = 0.046), or FAM102B (HR = 1.57, P = 0.011) expression had a shorter survival (Fig. S5). These results suggested that CD97 and genes positively related with CD97 are associated with the development and poor prognosis of HCC.

Fig. 4
figure 4

The top six genes co-expressed with CD97 in HCC. af The genes co-expressed with CD97 in HCC (r ≥ 0.66) were assessed using the UALCAN database. gl CD97 was positively correlated with IQGAP1 (r = 0.679, P < 0.001), GMIP (r = 0.673, P < 0.001), MOBKL2A (r = 0.642, P < 0.001), ARPC2 (r = 0.607, P < 0.001), CD68 (r = 0.568, P < 0.001), and FAM102B (r = 0.558, P < 0.001) in HCC in the TIMER database

Fig. 5
figure 5

Correlations of CD97 and its co-expressed genes with immune cell filtration levels in HCC. Correlations of CD97 (a), IQGAP1 (b), GMIP (c), CDH11 (d), MOBKL2A (e), ARPC2 (f), and FAM102B (g) expression with TIICs in HCC

Significant correlations of CD97 and its co-expressed genes with TIICs in HCC

We used the TIMER database to explore the relationship between CD97 expression and TIICs in the HCC microenvironment. The expression of CD97 was positively correlated with different TIICs, including B cells (r = 0.406, P < 0.001), CD8+ T cells (r = 0.482, P < 0.001), CD4+ T cells (r = 0.405, P < 0.001), neutrophils (r = 0.543, P < 0.001), and dendritic cells (r = 0.531, P < 0.001), and was significantly related to macrophages (r = 0.605, P < 0.001) (Fig. 5a). Further analysis of the relationships between CD97-co-expressed genes and TIICs in HCC showed that IQGAP1, GMIP, MOBKL2A, ARPC2, CD68, and FAM102B were significantly positively correlated with TIICs, especially macrophages (all r ≥ 0.281, P < 0.001; Fig. 5b–g). These results suggest that CD97 and its co-expressed genes are involved in the immune response in the TME by affecting TIICs, especially macrophages.

CIBERSORT analysis showed that CD97 expression was significantly associated with TIICs, including CD8+ T cells (P < 0.001), naive CD4 + T cells (P = 0.01), regulatory T cells (P = 0.007), resting NK cells (P = 0.005), monocytes (P = 0.007), M0 macrophages (P = 0.001), activated dendritic cells (P = 0.044), and resting mast cells (P < 0.001) (Fig. 6). In short, CD97 and its co-expressed genes showed close correlations with TIICs in HCC.

Fig. 6
figure 6

Violin plots visualizing the distributions of TIICs according to CD97 expression using the CIBERSORT database. Blue and red represent the CD97 low and high expression groups, respectively. There were significant differences in the numbers of CD8+ T cells (P < 0.001), naive CD4 + T cells (P = 0.01), regulatory T cells (P = 0.007), resting NK cells (P = 0.005), monocytes (P = 0.007), M0 macrophages (P = 0.001), activated dendritic cells (P = 0.044), and resting mast cells (P < 0.001)

CD97 is associated with M2 macrophage polarization

There was a significant positive correlation between CD97 and the infiltration level of macrophages (Fig. 5a). We explored the associations between CD97 expression and expression of markers of M1 macrophages, M2 macrophages, and TAMs in HCC using the TIMER database to determine the relationships between CD97 and the different macrophage subtypes. We found that the expression of the M1 macrophage gene marker IL1A (r = 0.278, P < 0.001) was significantly but weakly correlated with the expression of CD97 (Fig. 7a), while that of the M2 macrophage gene markers CD163 (r = 0.464, P < 0.001), MS4A4A (r = 0.536, P < 0.001), and VSIG4 (r = 0.525, P < 0.001) and the TAM markers CD86 (r = 0.661, P < 0.001), CCL2 (r = 0.479, P < 0.001), and IL10 (r = 0.508, P < 0.001) showed stronger positive correlations with CD97 expression (Fig. 7b, c). These findings suggest that high CD97 expression promotes M2 macrophage polarization and differentiation into TAMs, contributing to hepatocellular carcinogenesis.

Fig. 7
figure 7

Correlation of CD97 expression with macrophage polarization in HCC. a Correlations between CD97 expression and gene markers of M1 macrophages (NOS2, IL1A, and ARG2). b Correlations between CD97 expression and gene markers of M2 macrophages (CD163, MS4A4A, and VSIG4). c Correlations between CD97 expression and gene markers of TAM (CD86, CCL2, and IL10)

CD97 acts as an immune-related gene in HCC

We used the CIBERSORT algorithm to evaluate immune cell infiltration in the TME of HCC. We obtained data from 373 primary HCC and 50 normal hepatic samples from TCGA and found significant differences in the proportions of TIICs between HCC and normal hepatic tissues (Fig. 8a, b). There was an apparent discrepancy in the proportion of TIICs between the CD97 high and low expression groups in HCC (Fig. 8c). A heatmap showed that the proportions of different TIICs in HCC were strongly correlated with each other. There were significant negative correlations between M0 and M2 macrophages (r = −0.44), between CD8+ T cells and resting memory CD4+ T cells (r = −0.42), and between resting and activated NK cells (r = −0.42) and a significant positive correlation between CD8+ T cells and activated memory CD4+ T cells (r = 0.38) (Fig. 8d). We used the TIMER database to evaluate the relationships between the proportions of diverse types of TIICs and prognosis in HCC but found no significant associations (Fig. 9a). We also analyzed the relationship between CD97 and the immune checkpoint genes TOX, PD-L1 (or CD274), PD-L2 (or PDCD1LG2), CTLA4, and PD-1 (or PDCD1), which are critical targets of immunotherapy. The results showed that CD97 expression was positively correlated with the expression of the immune checkpoint genes and negatively correlated with tumor purity (Fig. 9b). In addition, we explored whether CD97 copy number alterations had a significant impact on the infiltration level of immune cells and CD97 mRNA expression (Fig. 10). Our results indicate that the copy number of CD97 may affect immune cell infiltration and thus the prognosis of HCC. Therefore, CD97 has a potential application in HCC immunotherapy.

Fig. 8
figure 8

Landscape of immune cell infiltration in HCC by CIBERSORT. a, b The composition of 22 types of immune cells in HCC and normal hepatic tissues. c Heatmap showing the differences in the infiltration levels of 22 types of immune cells between the CD97 high and low expression groups in HCC. d Correlation heatmap visualizing the correlations among 22 types of immune cells in HCC

Fig. 9
figure 9

Correlation analysis between CD97 expression and tumor immunity-related parameters. The associations between immune cells and prognosis in the TIMER database (a). The associations between CD97 and immune checkpoint expression (b)

Fig. 10
figure 10

The levels of CD97 and immune cell infiltration were associated with CD97 copy number alternations. ac The relationship between CD97 copy number and its mRNA level. d the relationships between CD97 copy number and immune cell infiltration levels, *P < 0.05

Discussion

HCC is a common malignant tumor of the digestive system [1]. Although there have been substantial advancements in the diagnosis and treatment of HCC over recent decades, patient prognosis remains poor [4]. Abnormal gene expression may be involved in tumor progression and is positively correlated with patient prognosis [5]. This study showed that CD97 is abnormally highly expressed in HCC tissues. Furthermore, CD97 overexpression is significantly associated with the prognosis and overall survival (OS) of HCC patients, suggesting that CD97 is an important factor that participates in HCC development. However, our knowledge of the role of CD97 molecules in HCC is limited. We aimed to comprehensively analyze the biological functions of CD97 in HCC using public databases and to reveal its related regulatory pathways.

Transcriptomic analysis of multiple databases revealed significantly higher CD97 expression in HCC tissues than in normal hepatic tissues. Furthermore, univariate and multivariate Cox survival analyses of HCC data in TCGA showed that CD97 was an independent adverse prognostic risk factor in HCC. We speculate that CD97 is highly expressed in HCC and can be used as a molecular marker for diagnosis and prognosis, but more clinical cohorts are needed for verification.

Based on the GSEA analysis, we identified 19 upregulated signaling pathways predominantly related to immune and inflammatory responses. CD97 is widely expressed in various immune cells, which plays an essential role in T cell activation [29]. The expression of CD97 is rapidly upregulated during the activation of lymphocytes and myeloid cells, which may promote adhesion and migration of these cells to inflammatory sites [8]. Kop et al. used a mouse model of collagen-induced arthritis to demonstrate that CD97 inhibitors reduce arthritis-related joint injury and inflammation by neutralizing antibodies [30]. These discoveries demonstrate that CD97 is positively correlated with immune cell infiltration. However, we need to further explore the relationship between CD97 expression and TIICs in HCC. Moreover, the significantly upregulated “cell adhesion molecules CAMs” signaling pathway participates in cell adhesion. In the current study, cell adhesion based on E-cadherin is important for maintaining epithelial integrity, which regulates the motility that spreads to the surrounding and distant tissues, affecting HCC progression [31]. CD97 has been proven to play a vital role in intercellular connections, which involves in the adhesion and migration of tumor cells and mediates the invasion of human tumors. Together, these results indicate that CD97 participates in tumorigenesis and immune-related KEGG signaling pathways. We showed that high CD97 expression in HCC predicted adverse outcomes. Therefore, CD97 may play a role in the development of HCC by affecting these signaling pathways.

Co-expression and TIMER analyses showed that CD97 expression was positively correlated with IQGAP1, GMIP, MOBKL2A, ARPC2, CD68, and FAM102B expression. These genes are involved in inflammation, the immune response, or tumorigenesis [32,33,34,35,36]. Genes positively related with CD97 were highly expressed in HCC tissues, except for CD68, and were associated with a poor prognosis in HCC patients, except for IQGAP1. TIICs are positively correlated with the progression and prognosis of HCC [37]. The mRNA expression levels of CD97 and its co-expressed genes were significantly and positively correlated with B cells, CD8+ T cells, CD4+ T cells, neutrophils, dendritic cells, and macrophages. At the same time, we evaluated the proportion of TIICs in HCC based on the differential expression of CD97 using CIBERSORT analysis and found significant differences in the proportions of eight types of immune cells: CD8+ T cells, naive CD4+ T cells, regulatory T cells, resting NK cells, monocytes, M0 macrophages, activated dendritic cells, and resting mast cells. The results further suggest that CD97 is strongly associated with TIICs in HCC and that CD97 and genes positively related with CD97 may be involved in the immune response to HCC development, leading to a poor prognosis.

We identified that CD97 expression in HCC was strongly associated with marker genes of M2 macrophages and TAMs, but there were only weak associations with the marker genes of M1 macrophages. These results indicate that high CD97 expression in HCC can promote the polarization of macrophages to M2 macrophages that eventually differentiate into TAMs, contributing to HCC development. The development and progression of HCC are complex and are associated with cumulative gene mutations in tumor cells but are regulated by the TME [38]. Scientific evidence shows that the TME plays a key role in tumor stimulation and progression. TAMs are vital components of the TME in HCC [38] and are composed of M1 and M2 macrophages. M1 macrophages, known as anti-tumor macrophages, activate host defenses and kill tumor cells. M2 macrophages, known as tumor-promoting macrophages, can promote tumor progression [39]. Previous studies have shown that tumor-promoting factors affect the polarization of macrophages to M2 macrophages in the TME. TAMs are similar to M2 macrophages and have similar characteristics [40]. TAMs play important roles in tumor growth by enhancing tumor cell resistance to chemotherapy and radiotherapy and promoting tumor angiogenesis (including invasion, infiltration, and metastasis) and immunosuppression [40]. Therefore, we speculate that high CD97 expression promotes the differentiation of macrophages into TAMs, accelerating their progression and leading to a poor prognosis in HCC patients.

There was a significant positive correlation between the expression of CD97 and immune checkpoint genes such as TOX, PD-L1, PD-L2, CTLA4, and PD-1. Immune checkpoint inhibitors have made significant breakthroughs in cancer treatment as new tumor therapies [41]. In HCC, the influence of immune checkpoints on the prognosis of patients is extensive and complex; increased expression of PD-1, PD-L1, and tumor-infiltrating lymphocytes in HCC induces immunosuppression [42]. Furthermore, high expression of CTLA-4 in HCC patients inhibits activation of T cell antigen presentation and proliferation and induces apoptosis, leading to immunosuppression and HCC development and progression [43]. Therefore, we believe that CD97 may affect the progression of HCC and the prognosis of patients by altering the expression of immune checkpoint genes such as TOX, PD-L1, PD-L2, CTLA4, and PD-1.

This study has certain limitations. First, our research data were obtained from the TCGA, GTEx, GEO, and TARGET databases, and our results must be verified using larger sample sizes. Moreover, the results obtained by the algorithm based on RNA sequences may not be accurate. Additional in vivo and in vitro experiments are required to explore the potential physiological mechanism of CD97 in HCC, including its role in tumor–immune interactions.

Conclusion

In summary, our results suggested that CD97 has prognostic value in HCC and affects tumor immunity. Patients with high CD97 expression had poorer clinicopathological features and worse prognoses compared with patients with low expression. High CD97 expression is associated with TIICs and can promote the conversion of macrophages to TAMs. Therefore, CD97 can be used as a prognostic marker in HCC.

Availability of data and materials

The raw data of this study are derived from the TCGA database (https://portal.gdc.cancer.gov/), GEO database (https://www.ncbi.nlm.nih.gov/geo/), and UCSC Xena database (http://xena.ucsc.edu), which are publicly available databases.

Abbreviations

GPCRs:

G protein-coupled receptors

HCC:

Hepatocellular carcinoma

GEO:

Gene Expression Omnibus

TCGA:

The Cancer Genome Atlas

KM:

Kaplan–Meier

GSEA:

Gene Set Enrichment Analysis

TIMER:

Tumor Immune Estimation Resource

GC:

Gastric cancer

GRK6:

G protein-coupled receptor kinase 6

EMT:

Epithelial–mesenchymal transition

KEGG:

Kyoto Encyclopedia of Genes and Genomes

GTEx:

Genotype-Tissue Expression

TARGET:

Therapeutically Applicable Research to Generate Effective Treatments

HPA:

The Human Protein Atlas

HR:

Hazard ratio

NES:

Normalized enrichment score

TME:

Tumor microenvironment

OS:

Overall survival

TIICs:

Tumor-infiltrating immune cells

TAMs:

Tumor-associated macrophages

References

  1. Forner A, Reig M, Bruix J. Hepatocellular carcinoma. Lancet. 2018;391(10127):1301–14.

    Article  PubMed  Google Scholar 

  2. Chen Y, Tian Z. HBV-induced immune imbalance in the development of HCC. Front Immunol. 2019;10:2048.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Levrero M, Zucman-Rossi J. Mechanisms of HBV-induced hepatocellular carcinoma. J Hepatol. 2016;64(1 Suppl):S84–S101.

    Article  CAS  PubMed  Google Scholar 

  4. Couri T, Pillai A. Goals and targets for personalized therapy for HCC. Hepatol Int. 2019;13(2):125–37.

    Article  PubMed  Google Scholar 

  5. Jiang Y, Han QJ, Zhang J. Hepatocellular carcinoma: mechanisms of progression and immunotherapy. World J Gastroenterol. 2019;25(25):3151–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Nieto Gutierrez A, McDonald PH. GPCRs: emerging anti-cancer drug targets. Cell Signal. 2018;41:65–74.

    Article  CAS  PubMed  Google Scholar 

  7. Yona S, Lin HH, Stacey M. Immunity and adhesion-GPCRs. Adv Exp Med Biol. 2010;706:121–7.

    Article  CAS  PubMed  Google Scholar 

  8. Hamann J, Eichler W, Hamann D, Kerstens HM, Poddighe PJ, Hoovers JM, et al. Expression cloning and chromosomal mapping of the leukocyte activation antigen CD97, a new seven-span transmembrane molecule of the secretion receptor superfamily with an unusual extracellular domain. J Immunol. 1995;155(4):1942–50.

    CAS  PubMed  Google Scholar 

  9. Eichler W, Hamann J, Aust G. Expression characteristics of the human CD97 antigen. Tissue Antigens. 1997;50(5):429–38.

    Article  CAS  PubMed  Google Scholar 

  10. Ward Y, Lake R, Martin PL, Killian K, Salerno P, Wang T, et al. CD97 amplifies LPA receptor signaling and promotes thyroid cancer progression in a mouse model. Oncogene. 2013;32(22):2726–38.

    Article  CAS  PubMed  Google Scholar 

  11. Ward Y, Lake R, Yin JJ, Heger CD, Raffeld M, Goldsmith PK, et al. LPA receptor heterodimerizes with CD97 to amplify LPA-initiated RHO-dependent signaling and invasion in prostate cancer cells. Cancer Res. 2011;71(23):7301–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Aust G, Steinert M, Schutz A, Boltze C, Wahlbuhl M, Hamann J, et al. CD97, but not its closely related EGF-TM7 family member EMR2, is expressed on gastric, pancreatic, and esophageal carcinomas. Am J Clin Pathol. 2002;118(5):699–707.

    Article  CAS  PubMed  Google Scholar 

  13. Steinert M, Wobus M, Boltze C, Schutz A, Wahlbuhl M, Hamann J, et al. Expression and regulation of CD97 in colorectal carcinoma cell lines and tumor tissues. Am J Pathol. 2002;161(5):1657–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Han SL, Xu C, Wu XL, Li JL, Liu Z, Zeng QQ. The impact of expressions of CD97 and its ligand CD55 at the invasion front on prognosis of rectal adenocarcinoma. Int J Color Dis. 2010;25(6):695–702.

    Article  Google Scholar 

  15. Chidambaram A, Fillmore HL, Van Meter TE, Dumur CI, Broaddus WC. Novel report of expression and function of CD97 in malignant gliomas: correlation with Wilms tumor 1 expression and glioma cell invasiveness. J Neurosurg. 2012;116(4):843–53.

    Article  CAS  PubMed  Google Scholar 

  16. Liu Y, Chen L, Peng SY, Chen ZX, Hoang-Vu C. Role of CD97(stalk) and CD55 as molecular markers for prognosis and therapy of gastric carcinoma patients. J Zhejiang Univ Sci B. 2005;6(9):913–8.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Li C, Liu DR, Li GG, Wang HH, Li XW, Zhang W, et al. CD97 promotes gastric cancer cell proliferation and invasion through exosome-mediated MAPK signaling pathway. World J Gastroenterol. 2015;21(20):6215–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Yin Y, Xu X, Tang J, Zhang W, Zhangyuan G, Ji J, et al. CD97 promotes tumor aggressiveness through the traditional G protein-coupled receptor-mediated signaling in hepatocellular carcinoma. Hepatology. 2018;68(5):1865–78.

    Article  CAS  PubMed  Google Scholar 

  19. Bartha Á, Győrffy B. TNMplot.com: a web tool for the comparison of gene expression in normal, tumor and metastatic tissues. Int J Mol Sci. 2021;22(5):2622.

  20. Mas VR, Maluf DG, Archer KJ, Yanek K, Kong X, Kulik L, et al. Genes involved in viral carcinogenesis and tumor initiation in hepatitis C virus-induced hepatocellular carcinoma. Mol Med. 2009;15(3-4):85–94.

    Article  CAS  PubMed  Google Scholar 

  21. Makowska Z, Boldanova T, Adametz D, Quagliata L, Vogt JE, Dill MT, et al. Gene expression analysis of biopsy samples reveals critical limitations of transcriptome-based molecular classifications of hepatocellular carcinoma. J Pathol Clin Res. 2016;2(2):80–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Cline MS, Craft B, Swatloski T, Goldman M, Ma S, Haussler D, et al. Exploring TCGA Pan-Cancer data at the UCSC Cancer Genomics Browser. Sci Rep. 2013;3:2652.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Menyhart O, Nagy A, Gyorffy B. Determining consistent prognostic biomarkers of overall survival and vascular invasion in hepatocellular carcinoma. R Soc Open Sci. 2018;5(12):181006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4.

    Article  PubMed  Google Scholar 

  29. Jaspars LH, Vos W, Aust G, Van Lier RA, Hamann J. Tissue distribution of the human CD97 EGF-TM7 receptor. Tissue Antigens. 2001;57(4):325–31.

    Article  CAS  PubMed  Google Scholar 

  30. Kop EN, Adriaansen J, Smeets TJ, Vervoordeldonk MJ, van Lier RA, Hamann J, et al. CD97 neutralisation increases resistance to collagen-induced arthritis in mice. Arthritis Res Ther. 2006;8(5):R155.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Giannelli G, Koudelkova P, Dituri F, Mikulits W. Role of epithelial to mesenchymal transition in hepatocellular carcinoma. J Hepatol. 2016;65(4):798–808.

    Article  CAS  PubMed  Google Scholar 

  32. Wei T, Lambert PF. Role of IQGAP1 in carcinogenesis. Cancers (Basel). 2021;13(16):3940.

  33. Hartmann EM, Campo E, Wright G, Lenz G, Salaverria I, Jares P, et al. Pathway discovery in mantle cell lymphoma by integrated analysis of high-resolution gene expression and copy number profiling. Blood. 2010;116(6):953–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Cheng Z, Wei W, Wu Z, Wang J, Ding X, Sheng Y, et al. ARPC2 promotes breast cancer proliferation and metastasis. Oncol Rep. 2019;41(6):3189–200.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Yamaguchi T, Fushida S, Yamamoto Y, Tsukada T, Kinoshita J, Oyama K, et al. Tumor-associated macrophages of the M2 phenotype contribute to progression in gastric cancer with peritoneal dissemination. Gastric Cancer. 2016;19(4):1052–65.

    Article  CAS  PubMed  Google Scholar 

  36. Baker KF, Skelton AJ, Lendrem DW, Scadeng A, Thompson B, Pratt AG, et al. Predicting drug-free remission in rheumatoid arthritis: a prospective interventional cohort study. J Autoimmun. 2019;105:102298.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Yu S, Wang Y, Hou J, Li W, Wang X, Xiang L, et al. Tumor-infiltrating immune cells in hepatocellular carcinoma: Tregs is correlated with poor overall survival. PLoS One. 2020;15(4):e0231003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Baig MS, Roy A, Rajpoot S, Liu D, Savai R, Banerjee S, et al. Tumor-derived exosomes in the regulation of macrophage polarization. Inflamm Res. 2020;69(5):435–51.

    Article  CAS  PubMed  Google Scholar 

  39. Deng L, He K, Pan Y, Wang H, Luo Y, Xia Q. The role of tumor-associated macrophages in primary hepatocellular carcinoma and its related targeting therapy. Int J Med Sci. 2021;18(10):2109–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Pathria P, Louis TL, Varner JA. Targeting tumor-associated macrophages in cancer. Trends Immunol. 2019;40(4):310–27.

    Article  CAS  PubMed  Google Scholar 

  41. Darvin P, Toor SM, Sasidharan Nair V, Elkord E. Immune checkpoint inhibitors: recent progress and potential biomarkers. Exp Mol Med. 2018;50(12):1–11.

    Article  PubMed  Google Scholar 

  42. Xu F, Jin T, Zhu Y, Dai C. Immune checkpoint therapy in liver cancer. J Exp Clin Cancer Res. 2018;37(1):110.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Zongyi Y, Xiaowu L. Immunotherapy for hepatocellular carcinoma. Cancer Lett. 2020;470:8–17.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors gratefully acknowledge the TNMplot, TIMER, TCGA, GEO, Kaplan–Meier plotter, UALCAN, and HPA database, which made the data available.

Funding

This research was supported by the National Natural Science Foundation of China under the grant [Grant No. NSF82060127], the Yunnan Provincial Organ Transplantation Clinical Medical Center under the grant [Grant No. 2020SYZ-Z-044], and the Postgraduate Innovation Fund of Kunming Medical University under the grant [Grant No. 2022B16].

Author information

Authors and Affiliations

Authors

Contributions

QS and JR conceived and designed the study. QS, LL, and XL performed the data curation and analysis. QS, JR, WL, and XZ analyzed and interpreted the results. QS, JR, YD, LH, and DW drafted and reviewed the manuscript. All authors reviewed and approved the final manuscript.

Corresponding author

Correspondence to Jianghua Ran.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the Medical Ethics Committee of the Affiliated Calmette Hospital of Kunming Medical University. Signed written informed consent was obtained from the patients and/or guardians.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Prognostic value of CD97 expression in HCC. a Correlation between CD97 expression and OS in HCC patients using the UALCAN database. b Correlation between CD97 expression and OS in HCC patients using the KM plotter. c OS analysis of HCC patients using the KM plotter database. OS, overall survival.

Additional file 2: Figure S2.

The mRNA expression of CD97 co-expressed genes in HCC using the UALCAN database. a-f the mRNA expression levels of IQGAP1, GMIP, MOBKL2A, ARPC2, CD68, and FAM102B in HCC in the UALCAN database.

Additional file 3: Figure S3.

The mRNA expression of CD97 co-expressed genes in HCC using the TNMplot database. a-f the mRNA expression levels of IQGAP1, GMIP, MOBKL2A, ARPC2, CD68, and FAM102B in HCC in the TNMplot database.

Additional file 4: Figure S4.

The prognostic value of CD97 co-expressed genes in HCC using the UALCAN database. a-f Correlations between OS and the mRNA levels of IQGAP1, GMIP, MOBKL2A, ARPC2, CD68, and FAM102B in HCC in the UALCAN database.

Additional file 5: Figure S5.

The prognostic value of CD97 co-expressed genes in HCC using the KM plotter database. a-f Correlations between OS and the mRNA levels of IQGAP1, GMIP, MOBKL2A, ARPC2, CD68, and FAM102B in HCC in the KM plotter database.

Additional file 6: Table S1.

GSEA pathways upregulated due to high expression of CD97.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Su, Q., Li, L., Li, X. et al. CD97 serves as a novel biomarker of immune cell infiltration in hepatocellular carcinoma. World J Surg Onc 20, 382 (2022). https://doi.org/10.1186/s12957-022-02829-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12957-022-02829-2

Keywords