- Research article
- Open Access
Distinct diagnostic and prognostic values of Glypicans gene expression in patients with hepatocellular carcinoma
BMC Cancer volume 21, Article number: 462 (2021)
In our current work, we aimed to investigate the expressions of glypican (GPC) family genes at the mRNA level and assess their prognostic significances in patients with hepatocellular carcinoma (HCC).
The pathological roles of GPC family genes were examined using bioinformatics analysis. The diagnostic values of GPC genes were explored with the Gene Expression Profiling Interactive Analysis. Moreover, the mRNA expression and prognostic values of GPC genes were assessed via the KM plotter database.
Our data showed that the expression of GPC-3 was dramatically increased in the liver tumor tissue. Moreover, the expressions of the other five GPC family members were not significantly different between the tumor and normal liver tissues (P > 0.05). Furthermore, the up-regulation of GPC-1 at the mRNA level was dramatically correlated to the reduced overall survival (OS) for all HCC patients (hazard ratio = 2.03, 95% confidence intervals =1.44–2.87, P = 4.1e-05) compared with its low-expression group. Besides, the prognosis of the Caucasians was related to most GPC family genes, while the prognosis of the Asian race was only related to the expression of GPC-2. Besides, for pathological factors, including stage, grade, AJCC, and vascular invasion, the higher the pathological grade and vascular invasiveness, the lower the expression levels of GPC family genes (P < 0.05). Finally, the expression levels of GPC-1, 2, and 3 in the hepatitis group were related to the poor prognosis of HCC in the risk factor (alcohol consumption and hepatitis) subgroup (P < 0.05).
Our findings indicated that GPC-3 was dysregulated in HCC compared with paracancerous tissues. The expression of GPC-1 could be used as a potent predictive index for the general prognosis of HCC. The pathology, patients, and risk factors might affect the prognostic value of GPC family genes in HCC.
Hepatic cancer is more frequently diagnosed in males compared with females, and a dramatically increasing number of patients with liver cancer have been reported in recent years. Most of the primary liver cancer is diagnosed as HCC . HCC ranks the second leading cause of tumor-associated mortality worldwide, accounting for more than 90% of all deaths from primary HCC. Although a great deal of effort has been made on HCC , the mechanism underlying the pathogenesis of HCC remains largely unexplored.
Glypicans (GPCs) are a group of heparan sulfate (HS) proteoglycans, which are identified to be associated with the exocytoplasmic surface of the plasma membrane through a glycosyl-phosphatidylinositol anchor. The GPC family consists of six members, namely GPC-1 ~ 6 [3,4,5]. These members show high expressions on the cell surface and in the extracellular matrix, which function mainly as mediators of growth factor signaling pathway. Recent investigations have shown that some of GPC genes are dysregulated in the pathogenesis of various tumors, which may probably participate in tumorigenesis, and these detectable proteins in the blood can be used as latent clinical indicators [6,7,8,9,10].
Like other tumors, the tumorigenesis of HCC can also be attributed to genetic and environmental origins. Moreover, dysregulated genes involved in the pathogenesis of tumors are the most prospective source of diagnostic and prognostic indices.
However, an in-depth investigation is required to comprehensively assess the diagnostic and prognostic values of GPC genes in HCC. In the present study, we aimed to explore the diagnostic and prognostic values of GPC family genes in HCC patients according to information derived from publicly accessible databases and bioinformatics assay.
Bioinformatics analysis of GPC family genes
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of GPCs were performed using the database for Annotation, Visualization and Integrated Discovery [11,12,13] (DAVID)v6.8 (accessed October 27, 2019; https://david.ncifcrf.gov/home.jsp;). The gene-gene interaction (GGI) network was established using the gene multiple association network integration algorithms [14,15,16] (GeneMANIA; http://www.genemania.org/; accessed October 27, 2019), and STRING [17, 18] (Search Tool for the Retrieval of Interacting Genes/Proteins) (STRINGv.10.0; https://string-db.org/;accessed October 27, 2019) was employed to build protein-protein interaction (PPI) network.
Correlation analysis and assessment of the diagnostic value
As a newly established interactive web server, GEPIA can evaluate the RNA sequencing data from the TCGA and GTEx projects , including 9736 tumors and 8587 normal samples, using a standard processing pipeline.
Construction of prognostic signature
The detailed information of GPC family genes was submitted to the online database Kaplan–Meier plotter [21, 22] (KM plotter; http://kmplot.com/analysis/index.php?p=service&cancer=liver_rnaseq), which can investigate the effects of 54,000 genes on overall survival (OS) in 21 types of cancers. The largest dataset includes liver, breast, ovarian, lung, and gastric cancers. The miRNA subsystems include 11,000 samples from 20 different types of cancer. The database includes gene chip and RNA-seq data sources for different databases, such as Gene Expression Omnibus (GEO), European Genome-phenome Archive (EGA), and The Cancer Genome Atlas (TCGA). Therefore, the KM plotter has been widely accepted for assessing the clinical impacts of individual genes on the survival of cancer patients, including HCC. In our current study, the gene expressions of GPC family genes and their prognostic values in HCC patients were explored via the KM plotter database.
The KM survival plots were drawn according to the medians of clinical parameters, including pathology factor (stage, grade, AJCC, vascular invasion), patient factor (sex, race, sorafenib treatment), and risk factor (alcohol consumption and hepatitis virus). Based on the above-mentioned clinical parameters, the patients were assigned to different subgroups.
SPSS v.22.0 software (IBM Corp, Armonk, NY, USA) was employed to conduct the statistical analysis. The data consisted of hazard ratio (HR), survival plot, 95% confidence intervals (CI), and log-rank. A p value of less than 0.05 was considered statistically significant.
Bioinformatics analysis of GPC family genes
Figure 1a shows the results of GO analysis. The functions of GPC family genes included anatomical structure morphogenesis, Wnt signaling pathway, planar cell polarity pathway, co-receptor activity involved in Wnt signaling pathway, planar cell polarity pathway, plasma membrane, integral component of plasma membrane, extracellular space, proteinaceous extracellular matrix, anchored component of membrane; Golgi lumen, lysosomal lumen, retinoid metabolic process, glycosaminoglycan biosynthetic process, glycosaminoglycan catabolic process, and HS proteoglycan binding. However, the KEGG pathway analysis using DAVID suggested that the functions of GPC family genes were involved in hsa05205: proteoglycans in cancer (Fig. 1b).
GeneMANIA was used to conduct a correlation analysis of GPC family members at the gene level, which revealed relationships in pathways, shared protein domains, co-localization, and co-expression between GPC1 to GPC6. A complex network was constructed with the GPC family genes and other related genes (Fig. 2a). GGI network suggested that the GPC family genes possessed a strong protein homology and co-expression pattern with each other. STRING analysis was performed to identify the interactions of GPC gene family members at the protein level. The PPI network indicated that GPC family genes contacted with each other directly or indirectly. GPC3 was shown to interact with GPC1, GPC2, GPC4, GPC5, and GPC6 in regards to gene co-occurrence, text-mining, co-expression, and protein homology. The GPC3 gene closely related to other GPC family genes was in the central position. There were relationships among GPC2, GPC4, GPC6, and GPC3 in texting, protein homology, and gene co-occurrence. Besides, there were relationships in text-mining and protein homology among GPC1, GPC5, and GPC3.
Correlation analysis and assessment of the diagnostic value
In the present study, we compared the expressions of GPC family genes between liver tumor tissues and paracancerous tissues using the GEPIA online tool. Our data showed that GPC-3 was significantly up-regulated in the liver tumor tissues (Fig. 3c). Moreover, the expressions of GPC-1, GPC-2, GPC-4, GPC-5, and GPC-6 were not significantly different between tumor and normal liver tissues (Fig. 3a, b, d, e, f).
Correlation analysis between the expressions of GPC family genes at the mRNA level and OS
Firstly, the correlation between the expressions of GPC-1 ~ 6 at the mRNA level and prognosis in HCC patients was investigated using on-line survival analysis database. The Affymetrix ID of GPC-1 in the KM plotter was 2817_at. We found that in overall HCC patients, the up-regulation of GPC-l at the mRNA level was dramatically correlated to a lower OS (HR = 2.03, 95%CI = 1.44–2.87, P = 4.1e-05, Fig. 4a). The up-regulation of GPC-l at the mRNA level was remarkably correlated to a lower OS in pathological stage (I ~ III), grade (I ~ III), AJCC_T (I ~ III) and none vascular invasion group (HR = 2.8, 95%CI = 1.48–5.27, P = 0.00092; HR = 2.59, 95% CI = 1.08–6.21, P = 0.0275; P = 0.0397, HR = 1.87, 95%CI = 1.02–3.43; HR = 3.1, 95%CI = 1.1–8.7, P = 0.0246; HR = 1.76, 95%CI = 1.04–2.96, P = 0.0318; HR = 2.63, 95%CI = 1.42–4.86, P = 0.0014; HR = 2.42, 95%CI = 1.32–4.42, P = 0.0032; HR = 2.5, 95%CI = 1.1–5.66, P = 0.0232; HR = 1.95, 95%CI = 1.04–3.67, P = 0.0342; HR = 2.13, 95%CI = 1.26–3.59, P = 0.0037, Fig. 9a-d). There were statistically significant differences in the males and Caucasians (HR = 2.46, 95%CI = 1.57–3.85, P = 4.5e-5; HR = 1.59, 95%CI = 1.01–2.51, P = 0.0434, Supplemental Table 1, 2). Our data also showed that the up-regulation of GPC-l at the mRNA level was remarkably correlated to a lower OS in non-alcohol consumption, non-hepatitis virus group (HR = 2.14, 95%CI = 1.35–3.39, P = 0.001; HR = 1.84, 95%CI = 1.18–2.89, P = 0.0066, Figs. 5, 7a) as well as alcohol consumption and hepatitis virus group (HR = 2.47, 95%CI = 1.25–4.88, P = 0.0069; HR = 2.38, 95%CI = 1.22–4.62, P = 0.0087, Fig. 6, 8a). These results were consistent with the general trend of HCC (P < 0.05, Fig. 4a).
The Affymetrix ID of GPC-2 in the KM plotter was 221914_at. The up-regulation of GPC-2 at the mRNA level was not correlated to OS for all HCC patients (HR = 1.36, 95%CI = 0.94–1.96, P = 0.1012, Fig. 4b). Figure 9a-d reveals that the pathological factors, including stage, grade, AJCC_T, and vascular invasion, affected the mRNA expression of GPC-2 and prognosis of HCC (P > 0.05). There were statistically significant differences in the Asian group (HR = 2.45, 95%CI = 1.33–4.5, P = 0.0029, Supplemental Table 2), but no significant difference in the sex group (P > 0.05, Supplemental Table 1). The up-regulation of GPC-2 at the mRNA level suggested a worse OS in the non-alcohol consumption group (HR = 2.08, 95% CI = 1.09–3.95, P = 0.022, Fig. 5b) and non-hepatitis virus group (HR = 1.63, 95% CI = 1.00–2.64, P = 0.047, Fig. 7b). Meanwhile, we analyzed such correlation in the alcohol consumption group and showed that its up-regulation implied a pernicious prognosis (HR = 0.49, 95% CI = 0.25–0.95, P = 0.031, Fig. 6b).
Figure 4 shows the prognostic value of GPC-3 (Affymetrix ID:2719_at) expression at the mRNA level. No significant correlation was found between the OS and GPC-3 expression for all HCC patients (P > 0.05, Fig. 4c). There was a significant correlation in stage (III), grade (II) and microvascular invasion pathological types (HR = 1.94, 95% CI = 1.08–3.51, P = 0.0252; HR = 1.74, 95% CI = 1.05–2.9, P = 0.0309; HR = 2.17, 95% CI = 1–4.72, P = 0.0442, Fig. 9a-b, d), but no significant difference between sex, race and sorafenib treatment subgroup (P > 0.05, Supplemental Table 1, 2, 3). We further analyzed the effects of risk factors on the prognosis of HCC. The up-regulation of GPC-3 at the mRNA level suggested a worse OS in the alcohol consumption group (HR = 2.78, 95% CI = 1.45–5.34, P = 0.0014, Fig. 6c) or non-hepatitis virus group (HR = 1.68, 95% CI = 1.05–2.70, P = 0.003, Fig. 7c). In contrast, the up-regulation of GPC-3 at the mRNA level was correlated to a better OS in the non-alcohol consumption group (HR = 0.59, 95% CI = 0.37–0.93, P = 0.023, Fig. 5c).
Subsequently, the effect of GPC-4 (Affymetrix ID:2239_at) was investigated. Figure 4 shows that the expression of GPC-4 exerted no effects on OS for all HCC patients (P > 0.05, Fig. 4d). In addition, the down-regulation of GPC-4 at the mRNA level suggested a better OS in the grade (II), AJCC_T (III) subgroup (HR = 1.69, 95% CI = 1.01–2.83, P = 0.0435; HR = 1.99, 95% CI = 1.04–3.83, P = 0.0343, Fig. 9b-c), whereas it indicated a worse OS in Caucasian HCC patients (HR = 1.56, 95% CI = 0.35–0.91, P = 0.0167, Supplemental Table 2). There were also no effects of risk factors, including alcohol consumption and hepatitis virus, on the prognosis of HCC, and no correlation between the expression of GPC-4 at the mRNA level and OS was found in different groups (P > 0.05, Figs. 5, 6, 7 and 8 D).
The Affymetrix ID of GPC-5 was 2262_at (Fig. 4e). The up-regulation of GPC-5 at the mRNA level had no relationship with the OS for HCC patients (P > 0.05, Fig. 4e), whereas such an up-regulation exhibited a favorable effect on OS in hepatitis virus HCC, sorafenib treatment group (HR = 0.42, 95% CI = 0.21–0.85, P = 0.0125, Fig. 8e; HR = 4.56, 95% CI = 0.94–22.27, P = 0.0439, Supplemental Table 3). The result also demonstrated that the expression of GPC-5 with female factors might affect the prognosis of HCC (HR = 1.59, 95% CI = 0.91–2.78, P = 0.0992, Supplemental Table 1).
We further assessed the prognostic value of GPC-6 expression at the mRNA level (Affymetrix ID:10082_at). The expression of GPC-6 had a significant impact on the OS in Caucasians and sorafenib-treated HCC subgroup (HR = 0.59, 95% CI = 0.36–0.97, P = 0.0364; HR = 0.15, 95% CI = 0.04–0.56, P = 0.0012, Supplemental Table 2, 3). The data showed that the expression of GPC-6 exerted no effect on the OS of general HCC patients as well as a pathological factor and risk factor subgroup (P > 0.05, Fig. 4f, Figs. 5, 6, 7 and 8 F, Fig. 9a-d).
In the present study, we aimed to assess the associations between GPC gene family members and HCC prognosis. We collected the clinical and pathological data of 364 HCC patients from GEO, EGA, and TCGA databases. The data were classified to discuss the prognostic value of GPC family genes by bioinformatics analysis. The study results showed that the expression level of GPC-1 was related to the overall prognosis of HCC. The lower the expression level of GPC-1, the longer the survival time of patients. Among them, GPC-1 was related to the prognosis of HCC in pathological stage, grade, AJCC, vascular invasion, sex, race, risk factors (alcohol consumption and hepatitis) subgroups, and the expression level of GPC-1 was negatively correlated with the OS of HCC. The lower the expression level of GPC-1, the better the prognosis. Besides, the mRNA expressions of GPC-2, 3, 4, 5, and 6 were not related to the overall prognosis of HCC. In the analysis of each subgroup, the mRNA expression of GPC-2 was related to OS with the impact factors of the race (Asian), alcohol consumption, and hepatitis. The mRNA expression of GPC-3 was related to OS in pathological stage (III), grade (II), vascular invasion (micro), sorafenib treatment, alcohol consumption, and hepatitis. The expression level of GPC-4 was related to poor prognosis OS of HCC in grade (II) and AJCC (III) and race (Caucasian and Asian). The expression level of GPC-5 was related to the poor prognosis of OS in sex (female), sorafenib treatment, and hepatitis. The expression level of GPC-6 was related to the poor prognosis of OS in the race (Caucasian) and sorafenib treatment. It could be seen that the significance of GPC family genes in the prognosis of HCC was not the same among different prognostic factors. Based on our findings, we speculated that GPC-1 was stable in predicting the clinical prognosis of each subgroup of HCC, and it could be a good prognostic biomarker. At present, there is no study on the relationship between GPC-1 and HCC prognosis. Our data provided valuable insights into future research and clinical work on HCC. Besides, the predictive effect of GPC-3 on the prognosis of HCC is not consistent among different centers [23,24,25,26]. Our results indicated that the expression of GPC-3 was not related to the overall prognosis of HCC, and this result was consistent with Chen IP’s report .
GPCs refer to a subset of cell-surface glycoproteins, among which HS glycosaminoglycan chains are covalently attached to a protein core. The GPC family genes are highly conserved across animal species, and they have critical functions in various biological processes . Moreover, GPC genes also play a fundamental role in some biological processes, and they have been considered to be the regulators of several cell signal transduction pathways [29, 30]. The GO and KEGG analyses also showed that GPC genes were mainly involved in Wnt signaling pathway, plasma membrane and planar cell polarity pathway, as well as glycosaminoglycan biosynthetic and catabolic processes. In our study, GGI and PPI networks showed that the GPC genes and other associated genes constructed an intricate network, in which they interacted with each other. Among them, GPC-3 played a core role in the network connection. The results suggested that GPC-3 also played a role in HCC.
GPC-1 is a cell surface HS proteoglycan that is over-expressed in a variety of solid tumors, while its expression is suppressed in most adult normal tissues. GPC-1 acts as a form of co-receptor for a range of signaling molecules, affecting signaling pathways, including Wnt, Hedgehog, TGF-β, and fibroblast growth factor . GPC-1 is an important clinical biomarker involved in the process of cancer onset, and it can be used as an important indicator of disease prognoses, such as breast cancer, lung cancer, colorectal cancer, glioma, pancreatic cancer, and esophageal cancer . Matsuda K et al. have found that GPC-1 may play a pivotal role in the ability of breast cancer cells to exhibit a mitogenic response to multiple heparin-binding growth factors and contribute to disease progression in this malignancy . Qian JY et al. have reported that GPC-1 in exosomes is identified as an early diagnostic biomarker for pancreatic adenocarcinoma, and the presence of GPC-1 in extracellular vesicles can serve as a predictor of RIAC outcome for patients with APC mutation . Although GPC-1 can be used as a diagnostic marker and a prognostic indicator in many cancers, the clinical value of GPC-1 in the diagnosis and prognosis of HCC has not been reported . Our study showed that the GPC-1 could be a good potential biomarker in HCC prognosis.
As an oncofetal glycoprotein, GPC-3 binds to the cell surface through a glycophosphatidylinositol anchor. To date, GPC-3 is the most well-documented GPC gene in its family in different cancers, including HCC, ovarian clear cell carcinoma, melanoma, squamous cell carcinoma of the lung, hepatoblastoma, nephroblastoma (Wilms tumor), yolk sac tumor, and some pediatric cancers . GPC-3 is expressed in over 70% of HCC cases . Capurro MI et al. have found that GPC-3 can trigger canonical Wnt signaling, leading to the accelerated growth of HCC . Besides its role as a biomarker, more and more attention has been paid to GPC-3 as a new therapeutic target molecule, and clinical trials targeting GPC-3 are in progress . Ortiz MV et al. have found that GPC-3-targeted immunotherapy is a key pediatric-specific consideration . Jeon Y et al. have concluded that GPC-3 expression is more frequently correlated to HCC with aggressive features in South Koreans . Zhang J et al. have revealed that the over-expression of GPC-3 can predict a poor OS in HCC patients. Although GPC-3 has been established as a vital prognostic index for HCC, the results remain controversial [23,24,25,26]. According to our research, the results showed that the expression of GPC-3 was increased in HCC tissues, but GPC-3 could not predict the overall prognosis of HCC. Chen IP et al. have studied 55 patients with early HCC who undergo initial hepatectomy between 1995 and 2010. They find that the GPC-3 expression is not significantly correlated with OS . This outcome is consistent with our result. The controversial conclusions might be attributed to the clinical cases. Although there are certain standards for data collection, some differences still exist in the data collection methods of each center. Each experimental sample contains different pathological, clinical and risk factors, leading to the biased results. Therefore, the data uniformity could not be guaranteed. Our data were derived from the public dataset and different centers. Information obtained from multiple centers can avoid error as much as possible.
A recent study has revealed that GPC-2 can positively regulate Wnt signaling in neuroblastoma, as silencing of GPC-2 inactivates Wnt/b-catenin signaling and reduces the expressions of target genes . As one of several mRNA transcripts, GPC-2 is highly expressed in several childhood cancers. Orentas RJ et al. have revealed that GPC-2 has a critical function in neurodevelopment, childhood cancer, and prostate cancer . Xu N et al. have reported that the GPC-2 plays a fundamental role in the progression and metastasis of prostate cancer and thus can be used as a candidate treatment target and a latent prognostic biomarker . Bosse KR has revealed that GPC-2 is a suitable tumor antigen in neuroblastoma . Although HCC is concerned, the diagnosis and prognosis associated with the GPC-2 expression at the mRNA level have not been reported. In our study, the result also revealed that the expression of GPC-2 was not related to the overall prognosis of HCC.
GPC-5 is a member of the GPC family, which is capable of binding to the external surface of the plasma membrane . GPC5 is a critical regulator of morphogens and growth factors during development . The GPC5 gene is located on chromosome 13q31.3, a region that is frequently mutated in many types of cancers . GPC-5 has been reported as a tumor suppressor gene in many cancers, such as pancreatic cancer, glioma, endometrial cancer, prostate cancer, and lung adenocarcinoma. Liu T et al. have shown that miR-709 enhances the invasion of HCC cells by mediating the GPC-5 expression . Hong X et al. have reported that GPC-5 plays a tumor-suppressive role in glioma and indicated that the inhibition of miR-301b suppresses the proliferation and invasion of glioma cells by increasing the expression of GPC-5 [50, 51]. Our study revealed that the expression level of GPC-5 was not related to the prognosis of OS in overall HCC patients.
Compared with GPC-1, GPC-2, GPC-3, and GPC-5, few studies have investigated the mechanism and prognostic significance of GPC-4 and GPC-6 expressions at the mRNA level in malignant tumors, and no related reports have been found in HCC. GPC-4 is a novel adipomyokine that enhances insulin signaling. Ussar et al. have initially shown that GPC-4, as a novel adipokine, is released from cells and adipose tissue of mouse explants, and this circulating GPC-4 interacts with the insulin receptor and enhances insulin receptor signaling and insulin sensitivity . Cao J et al. have demonstrated that GPC-4 participates in 5-fluorouracil (5-FU) resistance and pancreatic cancer stemness . Zhao D et al. have demonstrated that the T allele of GPC-4 may signify a risk factor for Epstein-Barr virus-associated and negative gastric carcinoma . As the newest member of the family, GPC-6 is most homologous to GPC4 and ubiquitously expressed . Karapetsas A et al. have reported that GPC-6 may function as a predictive indicator of CD8+ T-lymphocyte infiltration and a favorable prognostic index in early-stage ovarian cancer, showing significant clinical value in diagnosis, prognosis, and tumor immunobattling . Fan Cet al. have demonstrated that GPC-6 promotes proliferation, migration, and invasion of tumor cells in nasopharyngeal carcinoma . Our study also suggested that there was no difference in the expressions of GPC-4 and GPC-6 in normal liver and HCC, and these two genes had no predictive effect on the overall prognosis of HCC.
From our data graph, we could also see some interesting perspectives. Firstly, in terms of race, the prognosis of the Caucasian was related to most GPC family genes, while the prognosis of the Asian was only related to the expression of GPC-2. There were significant differences in the prognostic significance of GPC genes among different races. Secondly, for pathological factors, including stage, grade, AJCC, and vascular invasion, the higher the pathological grade and vascular invasiveness, the lower the expression levels of GPC family genes. The expression levels of GPC family genes could be used as an indicator of poor pathological prognosis, which was also consistent with common sense. Thirdly, in the risk factor (alcohol consumption and hepatitis) subgroup, the expression levels of GPC-1, 2, and 3 in the hepatitis group were related to the poor prognosis of HCC. According to our study, the reason for this opposite conclusion might be attributed to that the risk factor of drinking alcohol affected the expression levels of GPC-2 and GPC-3. The incidence, diagnosis, treatment, and medical status of HCC in each region were different, leading to the different prognosis of HCC. Consequently, the overall relationship between GPC genes and prognosis obtained from the public database would show some differences. We speculated that the pathology, patients, and risk factors might affect the GPC family’s prognostic value in HCC.
There were some limitations in this study that should be addressed. First, the study cohort was relatively small. Thus, larger studies are required to verify these findings. Besides, further studies from multiple centers with patients of various races are needed. Finally, our study was mainly based on the online website, which generates univariate analysis. Based on the present study, further explorations using multivariate analysis will aim to determine independent factors and potential molecular mechanisms. To address these issues, we are planning well-designed functional verification studies, including in vitro and in vivo models, shortly.
Collectively, this study aimed to assess the associations between HCC prognosis and the expression patterns of GPC family members. Our study found that GPC-3 was dysregulated in HCC compared with paracancerous tissues. Furthermore, the up-regulation of GPC-1 at the mRNA level was dramatically correlated to the reduced OS for overall HCC patients. Besides, the pathology, patients, and risk factors might affect the expressions of GPC genes and the clinical stage of HCC. Therefore, GPC-1 was a potentially prognostic biomarker for HCC. However, our findings still need to be further validated, and the prognostic values of other GPC genes still need to be prospectively confirmed in a larger number of patients.
Availability of data and materials
All analyzed data are included in this published article. The original data are available upon reasonable request to the corresponding author. The public access of links to all databases are open and available to all the visitors.
DAVID (https://david.ncifcrf.gov/home.jsp); GeneMANIA (http://www.genemania.org/); STRING (https://string-db.org/); GEPIA (http://gepia.cancer-pku.cn/index.html); Kaplan-Meier plotter (http://kmplot.com/analysis/index.php?p=service&cancer=liver_rnaseq);
Kyoto Encyclopedia of Genes and Genomes
Gene Expression Profiling Interactive Analysis
Gene Expression Omnibus
European Genome-phenome Archive
The Cancer Genome Atlas
Torbenson M. Hepatocellular carcinoma: making sense of morphological heterogeneity, growth patterns, and subtypes. Hum Pathol. 2020. https://doi.org/10.1016/j.humpath.2020.12.009.
Singal A, Lampertico P, Nahon P. Epidemiology and surveillance for hepatocellular carcinoma: new trends. J Hepatol. 2020;72(2):250–61. https://doi.org/10.1016/j.jhep.2019.08.025.
Li N, Spetz M, Ho M. The role of Glypicans in cancer progression and therapy. J Histochem Cytochem. 2020;68(12):841–62. https://doi.org/10.1369/0022155420933709.
Kaur S, Cummings B. Role of glypicans in regulation of the tumor microenvironment and cancer progression. Biochem Pharmacol. 2019;168:108–18. https://doi.org/10.1016/j.bcp.2019.06.020.
Li N, Gao W, Zhang Y, Ho M. Glypicans as cancer therapeutic targets. Trends Cancer. 2018;4(11):741–54. https://doi.org/10.1016/j.trecan.2018.09.004.
Xu D, Su C, Sun L, Gao Y, Li Y. Performance of serum Glypican 3 in diagnosis of hepatocellular carcinoma: a meta-analysis. Ann Hepatol. 2019;18(1):58–67. https://doi.org/10.5604/01.3001.0012.7863.
Zhang J, Zhang M, Ma H, Song X, He L, Ye X, et al. Overexpression of glypican-3 is a predictor of poor prognosis in hepatocellular carcinoma: an updated meta-analysis. Medicine. 2018;97(24):e11130. https://doi.org/10.1097/MD.0000000000011130.
De Robertis M, Arigoni M, Loiacono L, Riccardo F, Calogero R, Feodorova Y, et al. Novel insights into notum and glypicans regulation in colorectal cancer. Oncotarget. 2015;6(38):41237–57. https://doi.org/10.18632/oncotarget.5652.
Li Y, Li M, Shats I, Krahn J, Flake G, Umbach D, et al. Glypican 6 is a putative biomarker for metastatic progression of cutaneous melanoma. PLoS One. 2019;14(6):e0218067. https://doi.org/10.1371/journal.pone.0218067.
Tsuchiya N, Sawada Y, Endo I, Saito K, Uemura Y, Nakatsura T. Biomarkers for the early diagnosis of hepatocellular carcinoma. World J Gastroenterol. 2015;21(37):10573–83. https://doi.org/10.3748/wjg.v21.i37.10573.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13. https://doi.org/10.1093/nar/gkn923.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57. https://doi.org/10.1038/nprot.2008.211.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38(Web Server issue):W214–20.
Montojo J, Zuberi K, Rodriguez H, Kazi F, Wright G, Donaldson SL, et al. GeneMANIA Cytoscape plugin: fast gene function predictions on the desktop. Bioinformatics. 2010;26(22):2927–8. https://doi.org/10.1093/bioinformatics/btq562.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. https://doi.org/10.1101/gr.1239303.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–d613. https://doi.org/10.1093/nar/gky1131.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(null):D607–13.
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(W1):W98–w102. https://doi.org/10.1093/nar/gkx247.
Lonsdale J, Thomas J, Salvatore M, Phillips R, Lo E, Shad S, et al. The genotype-tissue expression (GTEx) project. Nat Genet. 2013;45(6):580–5.
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. https://doi.org/10.1098/rsos.181006.
Menyhárt O, Nagy Á, Győrffy B. Determining consistent prognostic biomarkers of overall survival and vascular invasion in hepatocellular carcinoma. R Soc Open Sci. 2018;5(12):181006. https://doi.org/10.1098/rsos.181006.
Moudi B, Heidari Z, Mahmoudzadeh-Sagheb H. Meta-analysis and systematic review of prognostic significance of Glypican-3 in patients with hepatitis B-related hepatocellular carcinoma. Virusdisease. 2019;30(2):193–200. https://doi.org/10.1007/s13337-019-00517-6.
Pan C, Wang X, Chen W, Tao C, Xu X, Jin L, et al. Reevaluation of glypican-3 as a prognostic marker in HCC using X-tile software. Med Oncol (Northwood, London, England). 2015;32(1):359.
Wang YL, Zhu ZJ, Teng DH, Yao Z, Gao W, Shen ZY. Glypican-3 expression and its relationship with recurrence of HCC after liver transplantation. World J Gastroenterol. 2012;18(19):2408–14. https://doi.org/10.3748/wjg.v18.i19.2408.
Yu MC, Lee YS, Lin SE, Wu HY, Chen TC, Lee WC, et al. Recurrence and poor prognosis following resection of small hepatitis B-related hepatocellular carcinoma lesions are associated with aberrant tumor expression profiles of glypican 3 and osteopontin. Ann Surg Oncol. 2012;19(Suppl 3):S455–63.
Yao M, Wang L, Fang M, Zheng W, Dong Z, Yao D. Advances in the study of oncofetal antigen glypican-3 expression in HBV-related hepatocellular carcinoma. Biosci Trends. 2016;10(5):337–43. https://doi.org/10.5582/bst.2016.01176.
Fico A, Maina F, Dono R. Fine-tuning of cell signaling by glypicans. Cell Mol Life Sci. 2011;68(6):923–9. https://doi.org/10.1007/s00018-007-7471-6.
Ho M, Kim H. Glypican-3: a new target for cancer immunotherapy. Eur J Cancer (Oxford, England : 1990). 2011;47(3):333–8.
Song HH, Shi W, Xiang YY, Filmus J. The loss of glypican-3 induces alterations in Wnt signaling. J Biol Chem. 2005;280(3):2116–25. https://doi.org/10.1074/jbc.M410090200.
Lin X. Functions of heparan sulfate proteoglycans in cell signaling during development. Development (Cambridge, England). 2004;131(24):6009–21.
Wang S, Qiu Y, Bai B. The expression, regulation, and biomarker potential of Glypican-1 in cancer. Front Oncol. 2019;9:614. https://doi.org/10.3389/fonc.2019.00614.
Matsuda K, Maruyama H, Guo F, Kleeff J, Itakura J, Matsumoto Y, et al. Glypican-1 is overexpressed in human breast cancer and modulates the mitogenic effects of multiple heparin-binding growth factors in breast cancer cells. Cancer Res. 2001;61(14):5562–9.
Qian JY, Tan YL, Zhang Y, Yang YF, Li XQ. Prognostic value of glypican-1 for patients with advanced pancreatic cancer following regional intra-arterial chemotherapy. Oncol Lett. 2018;16(1):1253–8. https://doi.org/10.3892/ol.2018.8701.
Xiao D, Dong Z, Zhen L, Xia G, Huang X, Wang T, et al. Combined Exosomal GPC1, CD82, and serum CA19-9 as multiplex targets: a specific, sensitive, and reproducible detection panel for the diagnosis of pancreatic cancer. Mol Cancer Res. 2020;18(2):300–10. https://doi.org/10.1158/1541-7786.MCR-19-0588.
Shimizu Y, Suzuki T, Yoshikawa T, Endo I, Nakatsura T. Next-generation cancer immunotherapy targeting Glypican-3. Front Oncol. 2019;9:248. https://doi.org/10.3389/fonc.2019.00248.
Baumhoer D, Tornillo L, Stadlmann S, Roncalli M, Diamantis EK, Terracciano LM. Glypican 3 expression in human nonneoplastic, preneoplastic, and neoplastic tissues: a tissue microarray analysis of 4,387 tissue samples. Am J Clin Pathol. 2008;129(6):899–906. https://doi.org/10.1309/HCQWPWD50XHD2DW6.
De Cat B, Muyldermans SY, Coomans C, Degeest G, Vanderschueren B, Creemers J, et al. Processing by proprotein convertases is required for glypican-3 modulation of cell survival, Wnt signaling, and gastrulation movements. J Cell Biol. 2003;163(3):625–35. https://doi.org/10.1083/jcb.200302152.
Nishida T, Kataoka H. Glypican 3-targeted therapy in hepatocellular carcinoma. Cancers. 2019;11(9):1139. https://doi.org/10.3390/cancers11091339.
Ortiz MV, Roberts SS, Glade Bender J, Shukla N, Wexler LH. Immunotherapeutic targeting of GPC3 in pediatric solid Embryonal tumors. Front Oncol. 2019;9:108. https://doi.org/10.3389/fonc.2019.00108.
Jeon Y, Kim H, Jang ES, Hong S, Kim JW, Yoon YS, et al. Expression profile and prognostic value of glypican-3 in post-operative south Korean hepatocellular carcinoma patients. APMIS Acta Pathologica, Microbiologica, et Immunologica Scandinavica. 2016;124(3):208–15. https://doi.org/10.1111/apm.12491.
Chen IP, Ariizumi S, Nakano M, Yamamoto M. Positive glypican-3 expression in early hepatocellular carcinoma predicts recurrence after hepatectomy. J Gastroenterol. 2014;49(1):117–25. https://doi.org/10.1007/s00535-013-0793-2.
Li N, Fu H, Hewitt SM, Dimitrov DS, Ho M. Therapeutically targeting glypican-2 via single-domain antibody-based chimeric antigen receptors and immunotoxins in neuroblastoma. Proc Natl Acad Sci U S A. 2017;114(32):E6623–e6631. https://doi.org/10.1073/pnas.1706055114.
Orentas RJ, Yang JJ, Wen X, Wei JS, Mackall CL, Khan J. Identification of cell surface proteins as potential immunotherapy targets in 12 pediatric cancers. Front Oncol. 2012;2:194.
Xu N, Wu YP, Yin HB, Xue XY, Gou X. Molecular network-based identification of competing endogenous RNAs and mRNA signatures that predict survival in prostate cancer. J Transl Med. 2018;16(1):274. https://doi.org/10.1186/s12967-018-1637-x.
Bosse KR, Raman P, Zhu Z, Lane M, Martinez D, Heitzeneder S, et al. Identification of GPC2 as an Oncoprotein and candidate immunotherapeutic target in high-risk neuroblastoma. Cancer Cell. 2017;32(3):295–309.e212.
Häcker U, Nybakken K, Perrimon N. Heparan sulphate proteoglycans: the sweet side of development. Nat Rev Mol Cell Biol. 2005;6(7):530–41. https://doi.org/10.1038/nrm1681.
De Cat B, David G. Developmental roles of the glypicans. Semin Cell Dev Biol. 2001;12(2):117–25. https://doi.org/10.1006/scdb.2000.0240.
Veugelers M, Vermeesch J, Reekmans G, Steinfeld R, Marynen P, David G. Characterization of glypican-5 and chromosomal localization of human GPC5, a new member of the glypican gene family. Genomics. 1997;40(1):24–30. https://doi.org/10.1006/geno.1996.4518.
Liu T, Zhang X, Sha K, Liu X, Zhang L. Wang B: miR-709 up-regulated in hepatocellular carcinoma, promotes proliferation and invasion by targeting GPC5. Cell Prolif. 2015;48(3):330–7. https://doi.org/10.1111/cpr.12181.
Hong X, Zhang Z, Pan L, Ma W, Zhai X, Gu C, et al. MicroRNA-301b promotes the proliferation and invasion of glioma cells through enhancing activation of Wnt/β-catenin signaling via targeting Glypican-5. Eur J Pharmacol. 2019;854:39–47. https://doi.org/10.1016/j.ejphar.2019.03.057.
Cao J, Ma J, Sun L, Li J, Qin T, Zhou C, et al. Targeting glypican-4 overcomes 5-FU resistance and attenuates stem cell-like properties via suppression of Wnt/β-catenin pathway in pancreatic cancer cells. J Cell Biochem. 2018;119(11):9498–512. https://doi.org/10.1002/jcb.27266.
Zhao D, Liu S, Sun L, Zhao Z, Liu S, Kuang X, et al. Glypican-4 gene polymorphism (rs1048369) and susceptibility to Epstein-Barr virus-associated and -negative gastric carcinoma. Virus Res. 2016;220:52–6. https://doi.org/10.1016/j.virusres.2016.04.005.
Veugelers M, De Cat B, Ceulemans H, Bruystens AM, Coomans C, Dürr J, et al. Glypican-6, a new member of the glypican family of cell surface heparan sulfate proteoglycans. J Biol Chem. 1999;274(38):26968–77. https://doi.org/10.1074/jbc.274.38.26968.
Karapetsas A, Giannakakis A, Dangaj D, Lanitis E, Kynigopoulos S, Lambropoulou M, et al. Overexpression of GPC6 and TMEM132D in early stage ovarian cancer correlates with CD8+ T-lymphocyte infiltration and increased patient survival. Biomed Res Int. 2015;2015:712438.
Fan C, Tu C, Qi P, Guo C, Xiang B, Zhou M, et al. GPC6 promotes cell proliferation, migration, and invasion in nasopharyngeal carcinoma. J Cancer. 2019;10(17):3926–32. https://doi.org/10.7150/jca.31345.
The authors wish to thank the curators of the Kaplan-Meier Plotter website for their valuable contributions to this work.
The study was supported by grants from the Sanming Project of Medicine in Shenzhen (No. SZSM201812055), Medical Scientific Research Foundation of Guangdong Province (A2020629), Shenzhen Key Medical Discipline Construction Fund (No. SZXK035), Shenzhen Fund for Guangdong Provincial High-level Clinical Key Specialties (No. SZXK035), National Nature Science Foundation of China (No.:81800489, 81560535, 81072321, 30760243, 30460143 and 30560133), 2009 Program for New Century Excellent Talents in University (NCET), Guangxi Nature Sciences Foundation (No.: GuiKeGong 1104003A-7), and Guangxi Health Ministry Medicine Grant (Key-Scientific Research-Grant Z201018). The present study is also partly supported by the Self-raised Scientific Research Fund of the Health and Family Planning Commission of Guangxi Zhuang Autonomous Region (Z2016318), the Basic Ability Improvement Project for Middle-aged and Young Teachers in Colleges and Universities in Guangxi (2018KY0110), Innovation Project of Guangxi Graduate Education (JGY2018037), and Research Institute of Innovative Think-tank in Guangxi Medical University (The gene-environment interaction in hepatocarcinogenesis in Guangxi HCCs and its translational applications in the HCC prevention). We would acknowledge the supported by the National Key Clinical Specialty Programs (General Surgery & Oncology) and the Key Laboratory of Early Prevention & Treatment for Regional High-Incidence-Tumor (Guangxi Medical University), Ministry of Education, China.
Ethics approval and consent to participate
Use of Not applicable.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The prognostic value of the mRNA expressions of GPC genes in different sex HCC patients.
The prognostic value of the mRNA expressions of GPC genes in different race HCC patients.
The prognostic value of the mRNA expressions of GPC genes in sorafenib treatment HCC patients.
About this article
Cite this article
Wang, JY., Wang, XK., Zhu, GZ. et al. Distinct diagnostic and prognostic values of Glypicans gene expression in patients with hepatocellular carcinoma. BMC Cancer 21, 462 (2021). https://doi.org/10.1186/s12885-021-08104-z
- Bioinformatics analysis