Skip to main content

A comprehensive computational analysis to explore the importance of SIGLECs in HCC biology

Abstract

Background

Hepatocellular carcinoma (HCC) is an aggressive, malignant cancer with a complex pathogenesis. However, effective therapeutic targets and prognostic biomarkers are limited. Sorafenib provides delaying cancer progression and survival improvement in advanced HCC. But despite 10 years of research on the clinical application of sorafenib, predictive markers for its therapeutic effect are lacking.

Methods

The clinical significance and molecular functions of SIGLEC family members were assessed by a comprehensive bioinformatic analysis. The datasets included in this study (ICGC-LIRI-JP, GSE22058 and GSE14520) are mainly based on patients with HBV infections or HBV-related liver cirrhosis. The TCGA, GEO, and HCCDB databases were used to explore the expression of SIGLEC family genes in HCC. The Kaplan–Meier Plotter database was used to evaluate relationships between the expression levels of SIGLEC family genes and prognosis. Associations between differentially expressed genes in the SIGLEC family and tumour-associated immune cells were evaluated using TIMER.

Results

The mRNA levels of most SIGLEC family genes were significantly lower in HCC than in normal tissues. Low protein and mRNA expression levels of SIGLECs were strongly correlated with tumour grade and clinical cancer stage in patients with HCC. Tumour-related SIGLEC family genes were associated with tumour immune infiltrating cells. High SIGLEC expression was significantly related to a better prognosis in patients with advanced HCC treated with sorafenib.

Conclusions

SIGLEC family genes have potential prognostic value in HCC and may contribute to the regulation of cancer progression and immune cell infiltration. More importantly, our results revealed that SIGLEC family gene expression may be used as a prognostic marker for HCC patients treated with sorafenib.

Peer Review reports

Background

Malignant tumours pose a serious threat to public health worldwide, including hepatocellular carcinoma (HCC), which is caused by various pathogenic factors, such as alcohol use and hepatitis [1]. HCC is currently the third-leading cause of cancer mortality [2]. Although some patients with HCC are cured by resection or liver transplant, overall survival (OS) is still low. The poor prognosis in HCC can be explained by the frequent diagnosis at a late stage [3]. The majority of HCC patients develop on the background of cirrhosis[4] and then generate a regional milieu that predisposes them to HCC[5]. Sorafenib, an inhibitor of multitarget tyrosine kinases, is used for the treatment of advanced-stage HCC [6]. However, the median life expectancy of patients with HCC treated with sorafenib is only 1 year. Despite 10 years of research on the clinical application of sorafenib, there are still no confirmed predictive factors for its therapeutic effect [7,8,9].

Sialic-acid-binding immunoglobulin-like lectins (SIGLECs) are a family of sialic acid-recognition proteins expressed primarily on leukocytes involved in the fine-tuning of leukocyte activities [10]. The roles of SIGLECs in neoplasms have already been described, and therapeutic agents directed at SIGLECs are under investigation [11]. SIGLEC are mainly expressed on immune cell surfaces. SIGLECs have been summarized in a review that they could play diverse regulatory functions in the tumor microenvironment (TME), while participating in HCC progression by various mechanisms, such as regulating cancer metastasis, and promoting cancer immune escape [12]. Thus, SIGLECs can become a target for regulating immunological processes [13]. The connection between SIGLEC expression and prognosis has been evaluated extensively. Previous studies have shown that SIGLECs family genes play important functions in HCC. SIGLEC2 is the best characterized member of the SIGLEC family [14], and its down-regulation has been found to predict worse overall survival of HBV-related HCC [15]. Wu et al. found that CD24-SIGLEC10 interaction may be involved in immune response, which could potentially become a therapeutic target for HCC patients [16]. Wang et al. reported that reduced SIGLEC7 expression could cause NK cell dysfunction in HCC patients [17]. Yamada reported an association between SIGLEC7 expression and prognosis in colorectal cancer (CRC) [18]. Wang et al. identified SIGLEC15 as a prospective target for tumour immunotherapy [19]. Ye et al. have shown that SIGLEC family genes may manipulate the tumour microenvironment via the chemokine axis and are associated with prognosis [20]. However, the function of SIGLEC family members in HCC remains uncertain, and comprehensive explorations of its biological effects and molecular pathways are essential for developing new prognostic indicators and therapeutic targets.

In this study, we explored this question by evaluating the expression of SIGLEC family members at the mRNA and protein levels using The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and Human Protein Atlas (HPA) databases. A multidimensional analysis was used to structure functional networks and genomic alterations associated with HCC, and its effect on tumour immunity was explored. Additionally, clinical characteristics associated with SIGLEC family genes and the prognostic value of these genes in HCC were analysed. Hence, this is a comprehensive computational analysis to explore the importance of SIGLECs in HCC biology and prognosis.

Methods

Differences in SIGLEC expression at the transcriptional level

Most HCC patients included in this study have a history of HBV infection or HBV-related liver cirrhosis. TCGA is a famous cancer gene project containing genomic information for 33 primary cancers and matching noncancer tissues [21]. GEPIA is a database that includes tumour and normal samples from TCGA and the Genotype-Tissue Expression (GTEx) dataset. GEPIA was used to explore the difference in SIGLEC expression between cancer and corresponding normal tissues in TCGA datasets [22]. UALCAN is a website for analysing tumour genomic data [23]. We used UALCAN to explore the correlations between the expression of SIGLEC family members and tumour staging and grading in HCC.

GEO is a public genome data storage centre for MIAME-compliant data. Biological and medical researchers can query and download gene expression profiles of interest [24]. GSE14520 [25] and GSE22058 [26] datasets were obtained from the GEO database. In addition, the International Cancer Genome Consortium (ICGC) was used to explore differences in SIGLEC gene expression at the transcriptional level. The specific dataset was ICGC-LIRI-JP, which contains 231 liver cancer samples, RNA-seq data, and clinical information [27].

Differential SIGLEC expression at the protein level

To investigate SIGLEC expression in HCC at the protein level, HPA was used to explore the immunohistochemical data for SIGLEC protein expression in HCC and normal tissues [28]. The SIGLEC protein expression levels in normal hepatic tissues were explored in the tissue module, which contains abundant tissue protein expression information for healthy humans.

Molecular network structure and functional clustering of SIGLECs in HCC

GeneMANIA is a platform for gene function speculation, analyses of gene lists, and gene prioritization [29]. In this study, SIGLEC family members were submitted to GeneMANIA to explore correlations with other genes, including coexpressed genes, associated protein domains, interactions, colocalizing genes, and genetic interactions.

WebGestalt [30] (http://www.webgestalt.org) is a tool for functional enrichment analyses. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway [31,32,33] and Gene Ontology (GO) functional analyses of SIGELCs and correlated genes were performed using WebGestalt. GO functional enrichment was evaluated in the biological process (BP), cellular component (CC), and molecular function (MF) categories. The enriched pathways were evaluated by KEGG pathway analysis. Metascape is a platform used for inclusive gene list annotation and source analysis for investigational biologists. We further analysed gene enrichment of the SIGLEC family using DisGeNET [34] in Metascape [35] (https://metascape.org/gp/index.html).

Analysis of tumour immune invasion

Tumour immune invasion was explored by using Tumour Immune Estimation Resource (TIMER) [36], which can be used for integrated analyses of cancer-infiltrating immune cells. SIGLECs-related immune infiltration was explored in HCC using TIMER. SIGLEC mRNA expression in HCC and relationships between SIGLEC expression and immune cell abundances were explored, including CD4+ T cells, B cells, and CD8+ T cells. The purity-corrected partial Spearman’s rho values were obtained. Using TIMER, the relationships between SIGLEC family member expression levels and PD1, PD-L1, and CTLA4 levels were evaluated.

Survival and mutation analysis of SIGLEC in HCC

cBioPortal is a cancer database for data integration, data mining, and visualization [37]. Gene alterations in SIGLECs were visualized using OncoPrint. cBioPortal was used to analyse SIGLEC alterations in TCGA-HCC cases. The analytical parameters included mRNA expression, copy number variants, mutations, and survival.

In this study, the prognostic value of SIGLEC family gene expression in HCC was explored using Kaplan–Meier plotter. This database is used to analyse correlations of the expression levels of various genes and survival in multiple kinds of cancer. The data sources include GEO, TCGA, and EGA. A Kaplan–Meier survival analysis was performed to contrast two patient cohorts. The HR with 95% confidence interval and log-rank P value were assessed. Clinical information in the database is regularly updated and maintained [38].

Steps of data extraction and analysis in different databases

HCCDB [39] (http://lifeome.net/database/hccdb/home.html). Time of login: April 2022. Version: HCCDB 2018. Data source of HCCDB: GEO data series, ICGC and TCGA-LIHC. The box plot of HCCDB18 (ICGC-LIRI-JP) included the SIGLECs differential expression was shown in Fig. 3A. And Fig. 3B, C shows the results of HCCDB 1 (GSE22058) and HCCDB 16 (GSE14520). A heatmap of SIGLEC family expression in tissue and prognosis was generated (Fig. 3D). The data of the datasets contained in the HCCDB database could be downloaded, including the number of adjacent samples and the number of HCC samples. This is where the content for Table 1 comes from.

Table 1 Information on HCCDB datasets

The Human Protein Atlas [28] (https://www.proteinatlas.org). Time of login: August 2021. Version: HPA 21.1. Data source of The Human Protein Atlas: A Swedish-based program. The immunohistochemical images of liver cancer tissue and normal liver tissue were obtained and the protein expression of SIGELC family members in HCC (Fig. 5) can be generated in HPA.

GeneMANIA [29] (http://genemania.org). Time of login: August 2021. Version: GeneMANIA 2.6. Data source of GeneMANIA: Genomics and proteomics data based on instantly updated publicly available databases. The network of SIGLECs can be generated in GeneMANIA (Fig. 6A).

PINA [40, 41] (https://omics.bjcancer.org/pina/home.action). Time of login: August 2021. Version: PINA v3.0. Data source of PINA: TCGA-LIHC. The network of SIGLECs in HCC can be generated in PINA (Fig. 6B).

WebGestalt [30] (http://www.webgestalt.org). Time of login: August 2021. Version: WebGestalt 2019. Data source of WebGestalt: WebGestalt supports 342 gene identifiers and 155 175 functional categories, as well as user-uploaded functional databases. The data were derived from large-scale omics studies. The functional enrichment analysis of SIGLEC genes (Fig. 6C–F). The multiple test correction was perfomed in Gene set enrichment analysis by WebGestalt.

Metascape [35] (https://metascape.org/gp/index.html). Time of login: August 2021. Version: Metascape 2019. Data source of Metascape: A broad set of instantly updated biological databases. Summary of enrichment analysis in DisGeNET (Fig. 6G).

TIMER [36] (https://cistrome.shinyapps.io/timer). Time of login: August 2021. Version: TIMER 1.0. Data source of TIMER: TCGA-LIHC. The correlation results between PD-1, PD-L1, CTLA-4 and SIGLEC genes (Fig. 7). The correlation results between T-cell infiltration, B-cell infiltration and SIGLEC genes (Table 4).

TISIDB [42] (http://cis.hku.hk/TISIDB/index.php). Time of login: August 2021. Version: TISIDB 2019. Data source of TISIDB: TCGA-LIHC. The correlations between SIGLEC mRNA expression levels and tumour grade in HCC (Fig. 8).

GEPIA [22] (http://gepia.cancer-pku.cn). Time of login: August 2021. Version: GEPIA 1.0. Data source of GEPIA: TCGA-LIHC. The correlations between SIGLEC mRNA expression levels and cancer stage in HCC (Figs. 4, 9).

cBioPortal [43] (https://www.cbioportal.org). Time of login: August 2021. Version: cBioPortal v5.1.0. Data source of cBioPortal: Multidimensional cancer genomics datasets. The overview of 279 SIGLEC gene mutations (Fig. 10A). The SIGLEC gene alteration frequencies (Fig. 10B). The survival in SIGLEC mutation group (Fig. 10C, D).

Kaplan–Meier Plotter [38] (http://kmplot.com). Time of login: August 2021. Version: Kaplan–Meier Plotter 2021. Data source of Kaplan–Meier Plotter: TCGA-LIHC, GSE20017 and GSE9843. The survival results correlated with SIGLEC expression (Fig. 11). The effect of SIGLEC expression on survival and prognosis in patients with advanced HCC treated with sorafenib (Fig. 12).

Statistical analysis

SIGLECs expression levels were compared by t-tests. Kaplan–Meier analyses were used to evaluate differences in prognosis. OS was defined as the time interval from diagnosis to death. Recurrence-free survival (RFS) was defined as the time interval from diagnosis to recurrence or last follow‐up. Progression-free survival (PFS) was defined as the time from random assignment in a clinical trial to disease progression or death from any cause. Log-rank P < 0.05 revealed a significant difference in prognosis. One-way ANOVA, Wilcoxon signed-rank tests, and logistic regression were used to analyse the correlations between clinicopathologic features and SIGLEC expression levels. The relationships between SIGLEC gene expression levels and immune scores were calculated by using Spearman methods. The data visualization were produced using the graphics package ggplot2 [44, 45] (Fig. 1).

Fig. 1
figure 1

Flow chart of analytical pipeline

Results

SIGLECs family genes are expressed at low levels in HCC

An overview of the research process is shown in Fig. 2.

Fig. 2
figure 2

Flow chart of the analysis process in this study

To explore SIGLEC family member expression in HCC, the ICGC and GEO databases were used for comparisons of expression levels between cancer tissues and correlated normal tissues. First, in the ICGC-LIRI-JP datasets, the levels of SIGLEC family members were significantly lower in HCC tissues than in adjacent nontumor tissues (Fig. 3A). Furthermore, based on the GSE22058 dataset, the levels of SIGLEC family members were lower in HCC tissues than in paired adjacent nontumor tissues (Fig. 3B). Similar results were obtained for the GSE14520 dataset (Fig. 3C). These results indicate that SIGLEC family mRNA expression levels were significantly downregulated in HCC tissues. A heatmap showed that the most highly differentially expressed genes included SIGLEC1 and SIGLEC7 (Fig. 3D). SIGLECs family members showed decreased expression in cancer tissues in different HCCDB datasets, and the differences were significant. Information about the datasets is summarized in Tables 1, 2 and 3.

Fig. 3
figure 3

Transcriptional expression of SIGLEC family members in ICGC and GEO databases. Expression levels of SIGLEC family genes in cancer and normal tissues based on data from the A ICGC database, B GEO database (GSE22058), and C GEO database (GSE14520). D Heatmap representation of differential expression. Diff indicates the number of differentially expressed genes, blue indicates downregulation, and red indicates upregulation (HCCDB)

Table 2 The clinical information of TCGA-HCC patients
Table 3 Details of the primary databases used in this study

The mRNA expression levels of SIGLEC family genes were explored using TCGA data. The expression levels of most SIGLEC genes were significantly lower in HCC tissues than in normal samples, and genes with statistically significant differences included SIGLEC1, 4, 5, 7, 8, 9, 11, 14, and 16. These results were consistent with those of the GEO analysis (Fig. 4). Since the majority of cases develop on the background of cirrhosis, SIGLECs expression in cirrhotic controls, normal healthy controls, adjacent normal controls and HCC (GSE25097) were analysed (Additional file 1: Fig. S6).

Fig. 4
figure 4

Expression levels of SIGLEC family genes in normal liver tissues and HCC tissues (TCGA database). The mRNA expression levels of most SIGLEC genes were significantly reduced in HCC (AN). ***P < 0.001

Protein expression of SIGELC family members in HCC

SIGLECs protein expression in HCC was investigated using data from the HPA. The results were similar to those of the mRNA-level analysis. SIGLECs protein expression levels were low in HCC tissues (Fig. 5). SIGLEC3, 4, 5, 8, 9, and 14 were expressed at significantly lower levels in HCC tissues than in normal tissues. SIGLEC1, 2, 6, 10, and 11 protein expression levels were low in both HCC tissues and normal liver tissues (Fig. 5). Overall, these results indicated that SIGLEC expression is significantly reduced in HCC.

Fig. 5
figure 5

Immunohistochemical results for SIGLECs in normal liver tissues and HCC tissues (HPA database)

Functional enrichment analysis of SIGLEC genes in HCC

A network of SIGLECs and correlated genes was created using GeneMANIA (Fig. 6A) and PINA (Fig. 6B). With the functional network diagram, specific negatively and positively correlated genes that interact with SIGLEC genes were obtained.

Fig. 6
figure 6

Functional enrichment of SIGLEC genes in HCC. A A gene network of SIGLEC and correlated genes was generated using GeneMANIA. B Interaction network analysis of SIGLEC family members using the PINA platform. C Cellular component. D Biological process. E Molecular function. F KEGG pathway analysis [31,32,33]. G Gene list enrichment

WebGestalt was used to explore the GO and KEGG pathway enrichment. Various biological processes, such as cell adhesion, biological adhesion, and neutrophil activation, involved in the immune response were significantly related to SIGLECs in HCC (Fig. 6C). Cellular components, including tertiary granules, tertiary granule membranes, and secretory granule membranes, were identified (Fig. 6D). SIGLECs were also enriched for terms in the molecular function category (Fig. 6E), such as sialic acid binding, CD4 receptor binding, and glycosphingolipid binding.

A KEGG analysis showed that the SIGELC family was involved in various pathways [31,32,33], including the haematopoietic cell lineage, cell adhesion molecules (CAMs), and the B-cell receptor signalling pathway (Fig. 6F). A gene list enrichment analysis of members of the SIGLEC family was performed using Metascape. The SIGLEC family was obviously enriched in tumour immunity, Guillain–Barre syndrome, immune suppression, and other functions (Fig. 6G).

Relationships between SIGELC family mRNA expression levels and the degree of immune cell infiltration in HCC

We used the TIMER database to analyse the correlations between SIGLEC family members and immune infiltration in HCC. The mRNA expression levels of most SIGLECs were significantly correlated with B-cell frequencies. The mRNA expression levels of SIGLEC1, 2, 3, 5, 6, 7, 8, 9, 10, 11, 14, 15, and 16 were obviously related to CD4+ T cells. In addition, the mRNA expression levels of SIGLEC1, 2, 3, 5, 6, 7, 8, 9, 10, 11, 14, and 16 were closely related to the level of CD8+ T-cell infiltration in HCC (Table 4). Correlations between gene expression levels and levels of PD-1 (Fig. 7A), PD-L1 (Fig. 7B), and CTLA-4 (Fig. 7C) were evaluated.

Table 4 Association between the expression levels of SIGLEC family members and the levels of T-cell infiltration and B-cell infiltration in HCC (TCGA database)
Fig. 7
figure 7

Correlations between the SIGELC family and the immune-related genes PD-1, PD-L1, and CTLA4 (TCGA database). A Correlations between gene expression levels of SIGLECs and PD-1. B Correlations between the gene expression levels of SIGLECs and PD-L1. (C) Correlations between gene expression levels of SIGLECs and CTLA-4

Relationships between SIGELC family gene expression levels and clinicopathological characteristics of patients with HCC

TCGA data were downloaded from GEPIA, UALCAN, and TISIDB. These data were used to analyse the correlations between SIGLEC mRNA expression levels and clinicopathological characteristics (including tumour grade and individual stage) in patients with HCC. SIGLECs family mRNA expression levels were correlated with cancer stage. Lower mRNA expression levels of SIGLECs corresponded with a worse tumour stage. Similarly, a worse tumour grade was related to lower mRNA expression levels of SIGLECs (Figs. 8, 9 and Additional file 1: Figs. S1–S2).

Fig. 8
figure 8

Correlations between SIGLEC mRNA expression levels and tumour grade in HCC. (TCGA database) A–L. *P < 0.05; **P < 0.01; ***P < 0.001

Fig. 9
figure 9

Correlations between SIGLEC mRNA expression levels and cancer stage in HCC. (TCGA database) A–N. *P < 0.05; **P < 0.01; ***P < 0.001

Overall, the mRNA expression levels of most SIGLEC family genes were correlated with the clinicopathological characteristics of patients with HCC.

Mutations in SIGLEC family genes in HCC and correlations with prognosis

As shown in Fig. 10A, the type and frequency of SIGLEC family gene mutations in eight HCC datasets (containing 1507 samples) were explored using cBioPortal. The alteration frequencies in TCGA-Firehose Legacy, TCGA-PanCancer Atlas, INSERM-Nat Genet 2015, and AMC-Hepatology were 11.9%, 11.7%, 8%, 5.2%, and 2.1%, respectively (Fig. 10).

Fig. 10
figure 10

Type and frequency of SIGLEC gene mutations in HCC (cBioPortal database). A Overview of SIGLEC gene mutations. B SIGLEC gene alteration frequencies. C Significantly shortened disease-specific survival in the SIGLEC mutation group. D Significantly shortened overall survival in the SIGLEC mutation group

The SIGLEC family gene mutations were remarkably correlated with prognosis in HCC; the SIGLEC mutation group showed a significantly shorter OS and DSS than the wild-type group.

Overall, there was an important connection between SIGLEC family gene mutations and poor prognosis in HCC. Patients with HCC with normal expression of SIGLEC family genes are expected to have a better prognosis.

Value of SIGLEC mRNA expression levels in predicting the prognostic of advanced HCC patients treated with sorafenib

Kaplan–Meier Plotter was used to analyse the correlations between the mRNA expression levels of SIGLEC family genes and prognosis. Lower SIGLEC mRNA expression levels were significantly correlated with shorter RFS and PFS in advanced HCC patients (Fig. 11).

Fig. 11
figure 11

Forest plots of HRs for RFS and PFS with SIGLEC expression (Kaplan–Meier Plotter)

In view of the lack of markers for the response to sorafenib, we analysed the expression levels of SIGLECs in patients with advanced HCC treated with sorafenib. Intriguingly, we found that the expression levels of SIGLECs were significantly correlated with the prognosis of patients with advanced HCC treated with sorafenib. High levels of SIGLEC expression were significantly associated with a better prognosis (Fig. 12 and Additional file 1: Figs. S3–S5).

Fig. 12
figure 12

Effect of SIGLEC expression on survival and prognosis in patients with HCC treated with sorafenib (Kaplan–Meier Plotter)

The above results indicate that SIGLEC expression is correlated with prognosis in patients with advanced HCC treated with sorafenib. In the future, the SIGLEC family may become a predictive index for evaluating the efficacy of sorafenib, which is very important for the prognostic evaluation of patients with advanced HCC treated with sorafenib.

Discussion

HCC has a poor prognosis and a high mortality worldwide, with an increasing number of cases every year [46, 47]. It is possible that the early stages of HCC could be cured. However, the diagnosis of early-stage HCC is difficult owing to the lack of obvious signs and diagnostic biomarkers [48]. Sorafenib prolongs the survival period of end-stage HCC by only several months [49]. Immunotherapy for HCC is still in its infancy compared to its application to other tumours [50]. Accordingly, HCC treatment is far from perfect. Therefore, it is imperative to explore therapeutic targets.

The functions of SIGLECs in the regulation of immune cell infiltration and cancer have been described by Macauley et al. [51]. In fact, SIGLECs are a family of sialic acid immunoglobulin receptors involved in immune discrimination. These receptors can regulate immune cells, such as T cells and B cells, in the adaptive and innate immune systems via their glycan ligands [51]. More recently, a report proposed that SIGLEC family genes in tumour cells are related to carcinogenesis [20]. Immune cell infiltration plays crucial functions in HCC progression [52]. Another study showed that CD24-SIGLEC10 interaction may be involve in immune responses[16] and reduced SIGLEC7 expression could cause NK cell dysfunction in HCC patients[17]. We observed strong correlations between the expression levels of SIGLECs and tumour-infiltrating lymphocytes. Taken together, these results indicate that SIGLECs are potential targets for cancer immunotherapy. The mechanisms by which SIGLECs interact with T cells and B cells in immune surveillance and immunoediting have been summarized previously [53]. Additionally, SIGLEC15 has been identified as an emerging target for immunotherapy [54]. Barkal et al. reported that CD24 is highly expressed and is related to the inhibitory receptor SIGLEC10 expressed by tumour‐associated macrophages to promote immune evasion in breast cancer and ovarian cancer [55]. Additionally, SIGLEC9 is involved in the innate immune response to cancer [56]. There are many studies of the relationships between SIGLEC family members and prognosis in patients with malignant tumours. Kensuke et al. observed a relationship between SIGLEC7 and prognosis in CRC; they found that the expression of SIGLEC7 in macrophages may become a novel prognostic biomarker for the efficacy of immunotherapy against metastatic CRC [18]. An in-depth understanding of the roles of SIGLECs in immune cell functions is expected to accelerate the development of cancer inhibitors [51].

Our results showed that the mRNA expression levels of the majority of SIGLECs genes were remarkably reduced in HCC tissues based on TCGA and GEO datasets. We obtained similar results at the protein level based on data from the HPA. Furthermore, we identified additional proteins, e.g., TNR, PTPN11, PLP1, CD19, PSG1, PSG2, and PSG7, related to SIGLEC by network analysis using GeneMANIA. The predicted biological functions of SIGLECs mainly included biological adhesion, neutrophil activation, and primary cell adhesion. These results were consistent with the molecular pathway correlated with HCC development. In an analysis of the associations of SIGLEC expression with clinicopathological features of patients with HCC, lower gene expression levels were correlated with a worse tumour grade.

High SIGELC family gene expression was remarkably correlated with a better prognosis in patients with HCC undergoing treatment with sorafenib. This improved prognosis may be explained by the relationship between high SIGLEC expression and the activation of the immune system. Sorafenib has several targets, including VEGFR, PDGFR, and RAF [57]. Additionally, SIGLEC5 may be related to VEGFR [58]. However, the interaction between these factors in improving prognosis needs to be further explored. In addition, the aetiology of liver disease, presence of cirrhosis, alpha-fetoprotein levels, tumour burden, and other factors affect the response to sorafenib. For example, high AFP and NLR have been identified as prognostic factors for poor OS in patients treated with sorafenib [59]. Therefore, future research aimed at the development of effective indicators should incorporate these variables.

Combination strategies involving sorafenib and immunotherapy are promising methods and are gaining attention in HCC research. It is now clear that SIGLEC family members function in the regulation of natural immunity. Our results reveal the molecular mechanism and prognostic significance of SIGLECs in HCC, providing a basis for the development of therapeutic and diagnostic strategies. Surprisingly, we found that SIGLECs are effective prognostic markers for the response to sorafenib in patients with advanced HCC; therefore, our findings have important clinical significance. Therefore, we explored the dual effects of SIGLEC family members to overcome the current limitations of HCC therapy. However, the predictive value of SIGLECs for the response to new first-line therapies (atezolizumab + bevacizumab) requires further research and exploration.

Conclusion

In general, our results contribute substantially to research on the roles of SIGLECs as prognostic markers and anti-HCC therapeutic targets. Our study also provides an important theoretical basis for the clinical value of SIGLECs as a prognostic marker for HCC patients treated with sorafenib.

Availability of data and materials

The datasets analysed in this research are available from the corresponding author. For the usage of R language, please contact the corresponding author. The accessible links of all the publicly available databases used in the study are as follows. The TCGA database is available at https://portal.gdc.cancer.gov/. The GEPIA database is available at http://gepia.cancer-pku.cn/. The Kaplan–Meier plotter is available at http://kmplot.com/analysis/. cBioPortal is available at https://www.cbioportal.org/. TIMER is available at https://cistrome.shinyapps.io/timer/. Linkedomics is available at http://linkedomics.org/. The HCCDB is available at http://lifeome.net/database/hccdb/.

Abbreviations

AFP:

Alpha fetoprotein

BP:

Biological process

CAMs:

Cell adhesion molecules

CC:

Cellular component

CRC:

Colorectal cancer

GEPIA:

Gene Expression Profiling Interactive Analysis

GEO:

Gene Expression Omnibus

GO:

Gene Ontology

GTEx:

Genotype-Tissue Expression

HBV:

Hepatitis B virus

HPA:

Human Protein Atlas

HCC:

Hepatocellular carcinoma

HCCDB:

Integrative Molecular Database of Hepatocellular Carcinoma

ICGC:

International Cancer Genome Consortium

KEGG:

Kyoto Encyclopedia of Genes and Genomes

MF:

Molecular function

OS:

Overall survival

ORA:

Over-Representation analysis

PFS:

Progression-free survival

RFS:

Recurrence free survival

SIGLECs:

Sialic-acid-binding Ig-like lectins

TIMER:

Tumour Immune Estimation Resource

TCGA:

The Cancer Genome Atlas

References

  1. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71(1):7–33.

    Article  PubMed  Google Scholar 

  2. Zhu C, Wang A, Zhang H, Wan X, Yang X, Chen S, Zhao H. Research progress and prospects of markers for liver cancer stem cells. World J Gastroenterol. 2015;21(42):12190–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Sun N, Lee Y, Zhang R, Kao R, Teng P, Yang Y, Yang P, Wang J, Smalley M, Chen P, et al. Purification of HCC-specific extracellular vesicles on nanosubstrates for early HCC detection by digital scoring. Nat Commun. 2020;11(1):4489.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Zhao J, Adams A, Roberts B, O’Neil M, Vittal A, Schmitt T, Kumer S, Cox J, Li Z, Weinman SA, et al. Protein arginine methyl transferase 1- and Jumonji C domain-containing protein 6-dependent arginine methylation regulate hepatocyte nuclear factor 4 alpha expression and hepatocyte proliferation in mice. Hepatology. 2018;67(3):1109–26.

    Article  CAS  PubMed  Google Scholar 

  5. Guo L, Sharma SD, Debes JD, Beisang D, Rattenbacher B, Vlasova-St Louis I, Wiesner DL, Cameron CE, Bohjanen PR. The hepatitis C viral nonstructural protein 5A stabilizes growth-regulatory human transcripts. Nucleic Acids Res. 2018;46(5):2537–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Poh W, Wong W, Ong H, Aung M, Lim S, Chua B, Ho H. Klotho-beta overexpression as a novel target for suppressing proliferation and fibroblast growth factor receptor-4 signaling in hepatocellular carcinoma. Mol Cancer. 2012;11:14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Bayard Q, Meunier L, Peneau C, Renault V, Shinde J, Nault JC, Mami I, Couchy G, Amaddeo G, Tubacher E, et al. Cyclin A2/E1 activation defines a hepatocellular carcinoma subclass with a rearrangement signature of replication stress. Nat Commun. 2018;9(1):5235.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Marisi G, Cucchetti A, Ulivi P, Canale M, Cabibbo G, Solaini L, Foschi FG, De Matteis S, Ercolani G, Valgiusti M, et al. Ten years of sorafenib in hepatocellular carcinoma: are there any predictive and/or prognostic markers? World J Gastroenterol. 2018;24(36):4152–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Brunetti O, Gnoni A, Licchetta A, Longo V, Calabrese A, Argentiero A, Delcuratolo S, Solimando A, Casadei-Gardini A, Silvestris N. Predictive and prognostic factors in HCC patients treated with Sorafenib. Medicina. 2019;55(10):707.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Huang P, Low P, Wang I, Hsu S, Angata T. Soluble Siglec-14 glycan-recognition protein is generated by alternative splicing and suppresses myeloid inflammatory responses. J Biol Chem. 2018;293(51):19645–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Angata T. Associations of genetic polymorphisms of Siglecs with human diseases. Glycobiology. 2014;24(9):785–93.

    Article  CAS  PubMed  Google Scholar 

  12. Jiang KY, Qi LL, Kang FB, Wang L. The intriguing roles of Siglec family members in the tumor microenvironment. Biomarker Res. 2022;10(1):22.

    Article  Google Scholar 

  13. Bornhöfft K, Viergutz T, Kühnle A, Galuska S. Nanoparticles equipped with α2,8-linked sialic acid chains inhibit the release of neutrophil extracellular traps. Nanomaterials. 2019;9(4):610.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Abdu-Allah HH, Watanabe K, Completo GC, Sadagopan M, Hayashizaki K, Takaku C, Tamanaka T, Takematsu H, Kozutsumi Y, Paulson JC, et al. CD22-antagonists with nanomolar potency: the synergistic effect of hydrophobic groups at C-2 and C-9 of sialic acid scaffold. Bioorg Med Chem. 2011;19(6):1966–71.

    Article  CAS  PubMed  Google Scholar 

  15. Ren X, Ji Y, Jiang X, Qi X. Down-regulation of siglec-2 (CD22) predicts worse overall survival from HBV-related early-stage hepatocellular carcinoma: a preliminary analysis from Gene Expression Omnibus. 2018. Biosci Rep. https://doi.org/10.1042/BSR20181423.

  16. Zhang P, Lu X, Tao K, Shi L, Li W, Wang G, Wu K. Siglec-10 is associated with survival and natural killer cell dysfunction in hepatocellular carcinoma. J Surg Res. 2015;194(1):107–13.

    Article  CAS  PubMed  Google Scholar 

  17. Tao L, Wang S, Yang L, Jiang L, Li J, Wang X. Reduced Siglec-7 expression on NK cells predicts NK cell dysfunction in primary hepatocellular carcinoma. Clin Exp Immunol. 2020;201(2):161–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Yamada K, Hazama S, Suzuki N, Xu M, Nakagami Y, Fujiwara N, Tsunedomi R, Yoshida S, Tomochika S, Matsukuma S, et al. Siglec-7 is a predictive biomarker for the efficacy of cancer vaccination against metastatic colorectal cancer. Oncol Lett. 2021;21(1):10.

    CAS  PubMed  Google Scholar 

  19. Wang J, Sun J, Liu L, Flies D, Nie X, Toki M, Zhang J, Song C, Zarr M, Zhou X, et al. Siglec-15 as an immune suppressor and potential target for normalization cancer immunotherapy. Nat Med. 2019;25(4):656–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Chen Z, Yu M, Guo L, Zhang B, Liu S, Zhang W, Zhou B, Yan J, Ma Q, Yang Z, et al. Tumor derived SIGLEC family genes may play roles in tumor genesis, progression, and immune microenvironment regulation. Front Oncol. 2020;10:586820.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Tomczak K, Czerwińska P, Wiznerowicz M. The Cancer Genome Atlas (TCGA): an immeasurable source of knowledge. Contemp Oncol. 2015;19:A68-77.

    Google Scholar 

  22. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45:W98–102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Chandrashekar D, Bashel B, Balasubramanya S, Creighton C, Ponce-Rodriguez I, Chakravarthi B, Varambally S. 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 

  24. Barrett T, Wilhite S, Ledoux P, Evangelista C, Kim I, Tomashevsky M, Marshall K, Phillippy K, Sherman P, Holko M, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2013;41:D991-995.

    Article  CAS  PubMed  Google Scholar 

  25. Roessler S, Jia H, Budhu A, Forgues M, Ye Q, Lee J, Thorgeirsson S, Sun Z, Tang Z, Qin L, et al. A unique metastasis gene signature enables prediction of tumor relapse in early-stage hepatocellular carcinoma patients. Can Res. 2010;70(24):10202–12.

    Article  CAS  Google Scholar 

  26. Burchard J, Zhang C, Liu A, Poon R, Lee N, Wong K, Sham P, Lam B, Ferguson M, Tokiwa G, et al. microRNA-122 as a regulator of mitochondrial metabolic gene network in hepatocellular carcinoma. Mol Syst Biol. 2010;6:402.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Zhang J, Baran J, Cros A, Guberman J, Haider S, Hsu J, Liang Y, Rivkin E, Wang J, Whitty B, et al. International Cancer Genome Consortium Data Portal—a one-stop shop for cancer genomics data. Database J Boil Databases Curation. 2011;2011:bar026.

    Google Scholar 

  28. Thul P, Åkesson L, Wiking M, Mahdessian D, Geladaki A, Ait Blal H, Alm T, Asplund A, Björk L, Breckels L, et al. A subcellular map of the human proteome. Science. 2017. https://doi.org/10.1126/science.aal3321.

    Article  PubMed  Google Scholar 

  29. Warde-Farley D, Donaldson S, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes C, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38:W214-220.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B. WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res. 2019;47(W1):W199-w205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545-d551.

    Article  CAS  PubMed  Google Scholar 

  34. Piñero J, Saüch J, Sanz F, Furlong LI. The DisGeNET cytoscape app: exploring and visualizing disease genomics data. Comput Struct Biotechnol J. 2021;19:2960–7.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, Benner C, Chanda SK. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu J, Li B, Liu X. 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 

  37. Cerami E, Gao J, Dogrusoz U, Gross B, Sumer S, Aksoy B, Jacobsen A, Byrne C, Heuer M, Larsson E, 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 

  38. Nagy Á, Munkácsy G, Győrffy B. Pancancer survival analysis of cancer hallmark genes. Sci Rep. 2021;11(1):6047.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Lian Q, Wang S, Zhang G, Wang D, Luo G, Tang J, Chen L, Gu J. HCCDB: a database of hepatocellular carcinoma expression atlas. Genomics Proteomics Bioinform. 2018;16(4):269–75.

    Article  Google Scholar 

  40. Du Y, Cai M, Xing X, Ji J, Yang E, Wu J. PINA 3.0: mining cancer interactome. Nucleic Acids Res. 2021;49(D1):D1351–7.

    Article  CAS  PubMed  Google Scholar 

  41. Wu J, Vallenius T, Ovaska K, Westermarck J, Mäkelä TP, Hautaniemi S. Integrated network analysis platform for protein-protein interactions. Nat Methods. 2009;6(1):75–7.

    Article  CAS  PubMed  Google Scholar 

  42. Ru B, Wong CN, Tong Y, Zhong JY, Zhong SSW, Wu WC, Chu KC, Wong CY, Lau CY, Chen I, et al. TISIDB: an integrated repository portal for tumor-immune system interactions. Bioinformatics. 2019;35(20):4200–2.

    Article  CAS  PubMed  Google Scholar 

  43. Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Luo W, Brouwer C. Pathview: an R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics. 2013;29(14):1830–1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Asrani S, Devarbhavi H, Eaton J, Kamath P. Burden of liver diseases in the world. J Hepatol. 2019;70(1):151–71.

    Article  PubMed  Google Scholar 

  47. Sia D, Villanueva A, Friedman S, Llovet J. Liver cancer cell of origin, molecular class, and effects on patient prognosis. Gastroenterology. 2017;152(4):745–61.

    Article  CAS  PubMed  Google Scholar 

  48. Jiang Y, Tie C, Wang Y, Bian D, Liu M, Wang T, Ren Y, Liu S, Bai L, Chen Y, et al. Upregulation of Serum Sphingosine (d18:1)-1-P potentially contributes to distinguish HCC including AFP-negative HCC from cirrhosis. Front Oncol. 2020;10:1759.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Xia J, Ozaki I, Matsuhashi S, Kuwashiro T, Takahashi H, Anzai K, Mizuta T. Mechanisms of PKC-mediated enhancement of HIF-1α activity and its inhibition by vitamin K2 in hepatocellular carcinoma cells. Int J Mol Sci. 2019;20(5):1022.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Johnston M, Khakoo S. Immunotherapy for hepatocellular carcinoma: current and future. World J Gastroenterol. 2019;25(24):2977–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Macauley M, Crocker P, Paulson J. Siglec-mediated regulation of immune cell function in disease. Nat Rev Immunol. 2014;14(10):653–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Yang C, Huang X, Liu Z, Qin W, Wang C. Metabolism-associated molecular classification of hepatocellular carcinoma. Mol Oncol. 2020;14(4):896–913.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Muenst S, Läubli H, Soysal SD, Zippelius A, Tzankov A, Hoeller S. The immune system and cancer evasion strategies: therapeutic concepts. J Intern Med. 2016;279(6):541–62.

    Article  CAS  PubMed  Google Scholar 

  54. Sun J, Lu Q, Sanmamed MF, Wang J. Siglec-15 as an emerging target for next-generation cancer immunotherapy. Clin Cancer Res. 2021;27(3):680–8.

    Article  CAS  PubMed  Google Scholar 

  55. Barkal A, Brewer R, Markovic M, Kowarsky M, Barkal S, Zaro B, Krishnan V, Hatakeyama J, Dorigo O, Barkal L, et al. CD24 signalling through macrophage Siglec-10 is a target for cancer immunotherapy. Nature. 2019;572(7769):392–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Hsu Y, Hsu F, Chiang M, Tsai D, Li F, Angata T, Crocker P, Chau L. Siglec-E retards atherosclerosis by inhibiting CD36-mediated foam cell formation. J Biomed Sci. 2021;28(1):5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Liu H, Nazmun N, Hassan S, Liu X, Yang J. BRAF mutation and its inhibitors in sarcoma treatment. Cancer Med. 2020;9(14):4881–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Antille C, Frei M, Sorg O, Tran C, Kaya G, Masouyé I, Fontao L, Prins C, Saurat JH. Amicrobial pustulosis of the folds associated with auto-immune disorders. A case report with an analysis of cytokine expression profile in skin lesions of cutaneous neutrophilic lupus. Dermatology (Basel, Switzerland). 2008;216(4):324–9.

    Article  CAS  PubMed  Google Scholar 

  59. Bruix J, Cheng AL, Meinhardt G, Nakajima K, De Sanctis Y, Llovet J. Prognostic factors and predictors of sorafenib benefit in patients with hepatocellular carcinoma: analysis of two phase III studies. J Hepatol. 2017;67(5):999–1008.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We recognize the GEO and TCGA platforms and contributors for providing imperative data. We sincerely thank Kanehisa Laboratories for offering the formal permission of KEGG pathway maps. We also thank the reviewers for providing valuable suggestions, which have substantially improved our research.

Funding

This work was endorsed by the National Natural Science Foundation of China (Nos. 81872504; 81372487).

Author information

Authors and Affiliations

Authors

Contributions

SY designed the investigation. WC, TC, HZ, ZB, and XZ compiled the data. SY, YJ, XY, and LP wrote the manuscript. HC modified the manuscript. All authors contributed to this research and authorized this version of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Huaidong Cheng.

Ethics declarations

Ethics approval and consent to participate

This research was conducted in accordance with the ethical guidelines of the Declaration of Helsinki. All procedures were performed in accordance with relevant guidelines. All of the above data were sourced from an open source database, ethics approval and consent to participate have been acquired by the corresponding database construction team, and the use and acquisition of these data were performed according to relevant data access policies and publication guidelines.

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.

Further analysis of SIGLEC family expression and its relationship with patient survival.

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

Yao, S., Chen, W., Chen, T. et al. A comprehensive computational analysis to explore the importance of SIGLECs in HCC biology. BMC Gastroenterol 23, 42 (2023). https://doi.org/10.1186/s12876-023-02672-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12876-023-02672-z

Keywords