Skip to main content

B4GALT5 high expression associated with poor prognosis of hepatocellular carcinoma

Abstract

Background

B4GALT5 is postulated to be an important protein in sugar metabolism that catalyzes the synthesis of lactosylceramide (LacCer). However, its role in hepatocellular carcinoma (HCC) remains unknown.

Method

We characterized the expression of B4GALT5 in HCC tissue compared to normal tissue, and explored its function of B4GALT5 in HCC by enrichment analysis based on its co-expressed gene set. Next, we checked whether B4GALT5 expression is correlated to immune infiltration level and clinical prognosis in hepatocellular carcinoma. Finally, we verified the expression of B4GALT5 using clinical samples evaluated by RT-PCR, and conducted in vitro experiments with B4GALT5-knockdown HCC cells to investigate the function of B4GALT5 in the HCC cell proliferation, migration and invasion.

Results

We found B4GALT5 mRNA and protein expression levels were significantly high in HCC tissue compared to normal tissue. The enrichment analysis of the gene sets that co-expressed with B4GALT5 showed specificity in HCC-related pathways and functions. Also, the expression pattern of B4GALT5 was significantly related to the immune infiltration level, especially CD4+ T cell and macrophage cells. B4GALT5 higher mRNA expression was associated with poor overall survival (OS) in HCC patients. Furthermore, In vitro experiments showed that depletion of B4GALT5 significantly inhibited HCC cell proliferation, migration and invasion. This study revealed the function and its mediated pathways of B4GALT5 in HCC, indicating that B4GALT5 may serve as a prognostic biomarker of HCC.

Peer Review reports

Background

Hepatocellular carcinoma (HCC) is the most common form of liver cancer and highly malignant with poor prognosis [1]. HCC has become a threat to global health, according to the global cancer report in 2018 [2]. More than 800 thousand new cases of HCC are reported worldwide each year, more than 87% HCC patients died and approximately half of them are Chinese [3]. It is a pressing demand to find novel prognostic biomarkers and therapeutic targets of HCC [4].Yet we know little about how it plays its role in HCC.

Aberrant glycosylation of receptors on cell surface often causes oncogenic transformation involved in the development and progression of tumor, including tumor cell proliferation, invasion, metastasis and angiogenesis [5]. Altered glycosylation of proteins is frequently attributed to abnormal expression of glycosyltransferases [6, 7]. As one of the seven beta-1,4-galactosyltransferase genes, B4GALT5 (beta-1,4-galactosyltransferase 5) catalyzes the synthesis of lactosylceramide (LacCer) via the transfer of galactose from UDP-galactose to glucosylceramide (GlcCer) [8, 9]. B4GALT5 is present on the cell surface and located in the Golgi complex, similar to other glycosyltransferases that function as adhesion molecules involved in matrix interactions, cell spreading and migration, and signal transduction cascades [10, 11]. Many studies have shown that glycosyltransferase is closely related to tumor occurrence and development [12, 13]. The relationship between B4GALT5 and several tumors, such as gynecological tumor and embryonic tumor, has been confirmed. However, the role of B4GALT5 in HCC has not been investigated to date.

In this study, we explored the expression and function of B4GALT5 in HCC, and examined the results using clinical samples and in vitro experiments. We found that B4GALT5 is significantly upregulated in HCC, and its overexpression is related to the poor prognosis of HCC patients. Also, the expression level of B4GALT5 is significantly related to the immune infiltration level in HCC, especially CD4+ T cell and macrophage cells. Our in vitro experiments showed that depletion of B4GALT5 significantly inhibited HCC cell proliferation, migration and invasion. In conclusion, B4GALT5 may be a potential biomarker for prognosis of patients with hepatocellular carcinoma.

Results

B4GALT5 is highly expressed in HCC and associated to poor prognosis

We identified the differentially expressed genes (DEGs) between HCC tissues and normal tissues based on the GSE14520 dataset (T=225, N=220) and drew the volcano plot, as shown in Fig. 1 (a). The differential expression gene set included 2,273 upregulated genes and 1,458 downregulated genes (|FC|>1.5). Among the DEGs identified above, B4GALT5 was one of the genes with drastically increased expression in HCC tissues compared to normal ones. The survival analysis between high and low expression groups showed that high B4GALT5 expression indicate poor prognosis (Fig. 1 b-c).

Fig. 1
figure 1

Differential expression and survival analysis of B4GALT5 in HCC. (a) Volcano plot of DEGs in HCC tissues and normal tissues in GSE14520 dataset. (b) B4GALT5 expression levels in HCC and normal tissues in GSE14520 dataset (p<0.05). (c) Kaplan-CMeier curves of overall survival of GSE14520 cohort grouped by B4GALT5 median expression level (p<0.05). (d) B4GALT5 expression levels in HCC and normal tissues in TCGA dataset (p<0.05). (e-f) Kaplan-CMeier curves of overall survival (OS) and disease-free survival (DFS) of TCGA cohort grouped by B4GALT5 median expression level (p<0.05)

We also used the LIHC dataset obtained from TCGA (T=366, N=160) to verify the association between B4GALT5 expression level and prognosis. Consistently, this also showed that the expression of B4GALT5 was higher in tumor samples than in normal samples (Fig. 1 d), suggesting its close relation with HCC. In addition, the patients with higher expression of B4GALT5 were found to have shorter overall survival (OS) and progression-free survival (PFS) than those in the normal group (Fig. 1 e-f). These results suggested that B4GALT5 high expression is unfavorable indicator of liver cancer patients.

B4GALT5 conjugated with clinicopathological parameters improve prognosis accuracy

The mean expression level and statistic tests of B4GALT5 in different subgroup divided by other clinicopathological variables of HCC were listed in Table 1 on GSE14520 and TCGA, respectively. It can be found that the B4GALT5 only shows statistical difference in clinical stage subgroup (I-II vs. III-IV, p-value < 0.01) and gender subgroup (male vs. female, p-value < 0.05) on TCGA HCC cohort. On GSE14520 cohort, B4GALT5 only has different expression level in AFP subgroups (high vs. low, p-value < 0.05). The results showed that B4GALT5 is another independent biomarker of prognosis in HCC.

Table 1 Statistical correlation between B4GALT5 and clinicopathological parameters

Next, we checked whether B4GALT5 expression level conjugated with another clinicopathological variable can improve the prognosis accuracy. Survivals of GSE14520 cohort were shown in Fig. 2 (a-d), the overall survival status can be better speculated using B4GALT5 expression level for the age ≥ 60 patients. B4GALT5 expression level in female patients yield to more accurate prognosis. Especially,we observed statistically significant difference of survival status between the high and low B4GALT5 expression in white subgroup and early stage (I-II) subgroup. For TCGA cohort, B4GALT5 conjugated with age ≥ 60, male, low ALT and cirrhosis acquired significantly as shown in Fig. 2 (e-h). We supposed that B4GALT5 may be a promising biomarker of prognosis especially in certain subgroup of HCC patients.

Fig. 2
figure 2

Survival analysis of B4GALT5 combined with clinicopathological parameters on TCGA cohort. (a-d) Kaplan-CMeier curves according to B4GALT5 signature with age ≥ 60, female, white, stage I-II patients in GSE14520 cohort(p<0.05). (e-h) Kaplan-CMeier curves of B4GALT5 signature with age ≥ 60, male, low ALT and cirrhosis patients in TCGA cohort(p<0.05)

B4GALT5 co-expressed genes are enriched in HCC-related pathways and functions

We identified the genes that were closely correlated with B4GALT5 based on the TCGA LIHC dataset. The genes with extremely low expression (TPM<0.5) in %50 samples were filtered out. Heatmap of the top 20 most positively and negatively correlated genes was shown in (Fig. 3 a). The set of co-expressed genes identified above were further used to explore the disease associations, gene patterns, GO functional enrichment and transcriptional factors by Metascape [14].

Fig. 3
figure 3

Enrichment analysis based on the set of B4GALT5 co-expressed genes in HCC. (a) The heatmap of the top 20 positively and negatively correlated genes of B4GALT5. (b-c) Bar plots of enriched terms in the DisGeNET10, PaGenBase, GO and TF-target knowledge repository of the set of B4GALT5 co-expressed genes

It can be found that B4GALT5 co-expressed genes are mainly enriched in the glucose-related diseases and lipid metabolism closely linked to HCC. Also, the set of co-expressed genes showed enrichment in liver tissue and HEPG2 cell line (Fig. 3 b). GO biological process were mainly concentrated in localization, metabolic process, growth and cellular process. Transcription factor targets included NR1H4, ZMYM2 and ZNF436, which might be metabolism-related and cancer-related factors (Fig. 3 c). These results further demonstrated that B4GALT5 is involved in the oncogenic pathways of HCC.

B4GALT5 expression is correlated with immune infiltration level in HCC

We conducted immune infiltration analysis on TCGA LIHC dataset and explored the correlation between B4GALT5 expression and immune infiltration levels (n=371). As shown in Fig. 4 (a-b), the heatmap and vioplot showed that T-reg cells have high infiltration level in HCC tissues than in normal tissue (p-value < 0.05). The Correlation heatmap showed the dendritic resting cell, NK activated cells B memory cells showed significantly positive correlation (Fig. 4c).

Fig. 4
figure 4

Immune cell counts in HCC tissues and association of B4GALT5 with immune infiltration level. (a) Heatmap of infiltrated immune cells in TCGA HCC samples. (b) Vioplot of expression of immune cells in HCC samples and normal samples. (c) Correlation heatmap of immune cells in HCC samples. (d) Correlation between B4GALT5 expression and immune infiltration levels evaluated by CIBERSORT and TIMER

Next, we investigated whether B4GALT5 expression is correlated with immune infiltration levels in HCC. Figure 4 (d) showed that B4GALT5 expression level had statistical correlation with the infiltration level of typical immune cells, particularly the NK resting cell, myeloid dendritic resting cell, mast activated cell. B4GALT5 expression level had significantly positive correlation with the infiltration level of macrophage cell, CD4+T cell and neutrophil. Table 2 results showed that B4GALT5 expression level was associated with immune cell infiltration. Especially, some types of immune cell infiltration in HCC had significantly higher correlation with B4GALT5 expression contrast to normal tissue, such as CD4 T, Monocyte and Th1 immune cells. These observations indicated that B4GALT5 played specific role in immune infiltration in HCC.

Table 2 Correlation between B4GALT5 expression and immune cell infiltration level on TCGA samples

Immunohistochemistry and qRT-PCR validate higher B4GALT5 expression in HCC tissues

The protein level of B4GALT5 in HCC tissue and matched adjacent tissue samples were analyzed using immunohistochemistry experiments. Figure 5 (a) showed different results of immunohistochemical staining. According to the total scores analysis, HCC tissues exhibited higher B4GALT5 protein expression level than paracarcinoma tissues (Fig. 5 b). Similarly, B4GALT5 mRNA expression levels were analyzed in the HCC tissues and paracarcinoma tissues by PCR, the results showed that the mRNA expression of B4GALT5 HCC tissues was higher than that in paracarcinoma tissues (Fig. 5 c).

Fig. 5
figure 5

B4GALT5 protein abundance of B4GALT5 in HCC tissues. (a) B4GALT5 protein abundance in normal, negative, neutral and positive expression in HCC tissue by IHC. (b) B4GALT5 protein level measured by IHC in HCC tissues and adjacent normal tissues(p<0.05). (c) B4GALT5 protein level measured by PCR in HCC tissues and adjacent normal tissues (p<0.05)

Knockdown of B4GALT5 significantly reduces proliferation, migration and invasion of HCC cells

To explore the function of B4GALT5 in HCC cells, we transfected siRNA into HCC cells and employed RT-qPCR and western blot to verify that the mRNA level of B4GALT5 was significantly decreased following si-B4GALT5 #3 transfection compared with si-NC transfection (Fig. 6 a). We then used these transfected cells for cell functional experiments. Cells transfected with #3 si-B4GALT5 exhibit significantly decreased proliferation compared with those in the control group (si-NC) and blank group (Fig. 6 b). Transwell and scratch assays showed that the number of invaded cells and the migration capacity of transfected cells are less than those cells in the control group, indicating that the invasive and migratory capacities of the cells with knockdown B4GALT5 were reduced (Fig. 6 c-e).

Fig. 6
figure 6

siRNA-induced Knockdown of B4GALT5 affect the capacity of proliferation, migration and invasion of HCC cells by in vitro experiments. (a) B4GALT5 transfection efficiency of liver cells evaluated by PCR and western blot. (b) Proliferation curve of CCK8 assay between B4GALT5 knockdown and control group (p<0.05). (c) Snapshot of transwell assay and scratch assay among B4GALT5 knockdown, control and blank group. (d-e) Invasion and migration efficiency of B4GALT5 knockdown cells by transwell assay (p<0.05)

Discussion

The malignant degree of HCC is high while surgi- cal treatment is only applicable for early-stage patients. Although various treatments, such as transcatheter arterial chemoembolization/(TACE), radiotherapy and locoregional therapy [15], have been developed, HCC still has a high mortality rate. The advanced stage and frequent recurrence after surgical resection are the main reasons of high mortality rate of HCC. Another reason lie in that the majority of the patients were limited to surgery or chemotherapy [16]. A few factors that are responsible for the initiation of HCC has been investigated, the underlying mechanisms for the conversion of healthy hepatic cells to neoplastic cells are still undefined. Therefore, early diagnosis and prognostic biomarkers are the focus of liver cancer research persistently [17].

Studies have shown that abnormal changes of glycosyltransferase are closely related to the occurrence and development of tumors [18]. In tumor cells the function of glycosyltransferase is aberrant because of genetic alterations and/or epigenetic modifications, which results in abnormal extracellular matrix and cell-cell interactions [19]. As aberrant glycans is involved in promoting tumor progression and metastasis [20], these enzymes are considered as therapeutic targets. In fact, genetic alterations in some of these glycotransferases have been gradually confirmed to be closely associated with liver cancer. For example, FUT8 activates PI3K-Akt-NF- κB signaling pathway and promotes the proliferation of HCC cells, it also stimulates the expression of drug-resistant proteins. Inhibition of FUT8 can weaken hepatocyte epidermal growth factor (EGF) and hepatocyte growth factor (HGF), and thus reduce the incidence of liver cancer [21, 22]. miR-23a may affect the formation of N-glycochain branches on the cell surface through glucose transferase MGAT3, thereby increasing the metastasis potential of liver cancer, which provides a new idea for the mechanism of action of MGAT3 in tumor metastasis [23, 24]. B4GalT5 participates in the synthesis of both N-linked oligosaccharides and various glycolipids, it involves in extraembryonic development and tumorigenesis in many kinds of tumors, such as astrocytoma, glioma, breast cancer and cervical cancer [9, 2527]. In fact, B4GalT5 has been reported to be one of the target genes connected to adverse clinical outcomes in patients with HCC [28].

This study sets about to clarify the role of B4GALT5 expression and its prognostic significance in HCC, by combining large-scale cohort analysis and in vitro experiments. On the one hand, glycosylation is closely related to glucose and lipid metabolism and insulin resistance, which promotes the occurrence and development of liver cancer [29]. Both our results and previous studies indicate this point. For instance, the metabolic status in mice after B4GALT5 gene knockout is improved [8, 30]. Moreover, glycosylation is related to tumor immune microenvironment which leads to poor prognosis [31]. Our results show that the degree of immune infiltration level in HCC is significantly correlated with B4GALT5 expression, especially in macrophages and neutrophils. Our finding is actually consistent to the conclusions of previous reports. For instance, the precursor of B4GALT5 β-GlcCer and its product LacCer are shown to be associated with tumor immune microenvironment. β-GlcCer was identified as an endogenous ligand for Mincle and possess immunostimulatory activity [32], while LacCer and kinase Lyn that form microdomains rich in spingolipids on the plasma membrane of neutrophils play an important role in innate immune function of neutrophils [33, 34]. In fact, many studies have suggested that B4GALT5 may be an inflammatory factor and a regulator of M1 infiltration, which increases inflammation and insulin resistance [35, 36]. As the regulation and functional mechanism of B4GALT5 still remains unclear, our future work will focus on the inflammatory reaction and immune regulation to reveal the role of B4GALT5 in HCC.

Conclusion

In this study, we demonstrated for the first time that B4GALT5 is over-expressed in HCC and positively associates with poor survival of HCC patients. We also found that B4GALT5 combined with another clinicopathological parameter can improve prognosis accuracy. The gene set enrichment analysis showed that B4GALT5 co-expressed genes were enriched in HCC-related pathways and functions. We also conducted immune infiltration level analysis and confirmed the association of B4GALT5 with tumor immune microenvironment in HCC. Next, we used clinical samples and cell experiments to verify the function of B4GALT5 in HCC cells. Our findings was actually consistent with the conclusions of previous studies that high B4GALT5 expression is a biomarker of poor prognosis of HCC patients. As a result, we suggested that inhibition of B4GALT5 may be a promising therapeutic target to treat HCC.

Methods

Expression and survival analysis of B4GALT5 in HCC

The gene expression dataset of human hepatocellular carcinoma was obtained from GEO GSE14520 [37]. This dataset included 239 HCC tumor samples and 249 normal samples. The R package limma [38] was used to run differential expression analysis. The volcano plot were produced by limma. The survival analysis for differentially expressed gene set were analyzed using Graphpad Prism (version 8.2.1 Windows version, GraphPad Software, San Diego). Another dataset obtained from TCGA LIHC project [39] contained 371 tumor samples and 50 normal samples was used to verify gene expression and prognosis.

Survival analysis based on B4GALT5 coupled with clinicopathological parameters

The clinical features of HCC patients were downloaded from GEO and TCGA. Each cohort was grouped according to different features, and then the mean expression levels of B4GALT5 in each subgroups were computed. The p-value was obtained by group t-test. A three-line table was used to show the relationships between clinical feature and B4GALT5 expression. Next, for each subgroup of patients divided by a clinical feature, survival analysis was run on the subgroups stratified by mean B4GALT5 expression level. The curves of survival analysis with statistical significance were presented.

Enrichment analysis of B4GALT5 co-expressed genes

The co-expressed genes of B4GALT5 were selected based on Spearman correlation coefficient, using the criteria of p-value less than 0.01. The set of co-expressed genes were filtered and clustered based on their membership similarities. Top 20 positively and 20 negatively co-expressed genes were used for further enrichment analysis. Metascape [14] was used for functional analysis on the set of co-expressed genes. The results included gene set enrichment analysis in disease association in DisGeNET database [40], transcription factor targets, and tissue/cell type association in PaGenBase database [41].

Immune infiltration level analysis

The immune-related gene expressions in HCC were screened based on gene expression matrix. The heatmap, vioplot and correlation plot of immune cells were drawn by R package. To further identify the relationship between B4GALT5 and immune infiltration level, the Spearman correlation of B4GALT5 expression with typical immune cells, including B cells, CD4+ T cells, CD8+ T cells, neutrophils and macrophages, was calculated.

The R package immunedeconv [42] was used to evaluate the scores of immune cells with different B4GALT5 expression. CIBERSORT [43] and TIMER [44] algorithms were both used to estimate immune cell count in HCC tissues. TIMER considered tissue specificity in the estimation of immune cell counts, while CIBERSORT can assess T cell characteristics in more detail (* p<0.05, ** p<0.01, *** p<0.001). The correlation between B4GALT5 expression and immune cells were presented in tabular form. p-value less than 0.5 is considered as statistical significance between the infiltration level of immune cell and B4GALT5.

Immunohistochemistry experiments

A total of 24 hepatocellular carcinoma tissue samples and matched adjacent normal tissue samples were acquired from the Affiliated Changzhou No. 2 People’s Hospital of Nanjing Medical University (Jiangsu, China). All tissues were obtained by surgery. The histology of the 24 hepatocellular carcinoma tissue specimens was confirmed by the senior pathologist of the Department of Pathology at the hospital. First, tissue sections were treated for antigen retrieval using EDTA buffer, and then the tissues were incubated with a primary anti-B4GALT5 antibody (1:200, Biorbyt) at 4C overnight. Next, the sections were incubated with a secondary antibody and then visualized using diaminobenzidine (DAB) and hematoxylin counterstaining after washing with PBS. The immunohistochemical results were evaluated by pathologists and scored as follows: 0, negative; 1, +; 2, ++; and 3, +++. The positive staining rate was defined according to the proportion of positively stained cancer cells: 0, negative; 1, 1-20%; 2, 21-40%; 3, 41-60%; 4, 61-80%; and 5, 81-100%. The staining intensity score plus the staining positive rate score is counted as the total score.

Quantitative real-time PCR

Total RNA was isolated from cells using an RNA extraction kit (centrifugal column assay) according to the manufacturer’s instructions. Quantitative RT-PCR (q RT-PCR) was performed on an ABI ViiA7 Series PCR instrument (Applied Biosystem, ThermoFisher, USA) using FastStart Universal SYBR Green Master Mix (ROX) (Roche, Sigma-Aldrich). The primer sequences were as follows: B4GALT5, forward 5’-TCCTCGCTGCTGTACTTCG-3’ and reverse 5’-AATGCCTTGGGCTTGCATCA-3’; 18S (serving as the internal reference), forward 5’-TCCTCGCTGCTGTACTTCG-3’ and reverse 5’-TTACAGGGCCTCGAAAGAGTCC-3’. Reactions were carried out at 95C for 30s, followed by 40 cycles at 95C for 5s and 60C for 30s. The fluorescence data were collected in the 60C extension phase. Each sample was measured in 3 technical replicates.

Western blot

Western blot was used to verify the transfection efficiency.The total protein of cells was extracted using radioimmunoprecipitation assay (RIPA) lysis buffer. The intact protein was separated by polyacrylamide gel electrophoresis (PAGE) and transferred onto a nitrocellulose membrane using a wet transfer method. The membrane was blocked with 5% BSA for 1h and then probed at 4C overnight with primary antibodies specific for B4GALT5 (1:1000, Bioworld, 45 KD) or Alpha tubulin (1:2000, 55 KD); Alpha tubulin served as the internal reference. After washing with TBST, the membrane was incubated with horseradish peroxidase (HRP)-labeled goat anti-rabbit IgG (1:10,000), washed with TBST, developed with an enhanced chemiluminescence (ECL) solution, and imaged for analysis,and the greater the amount of target protein, the less obvious the knockdown efficiency.

Cell culture and transfection

SK-Hep-1 human hepatocellular carcinoma cells were cultured in DMEM (Gibco, USA) supplemented with 10% FBS (Gibco, USA) in 5% CO2 at 37C. B4GALT5 silencing was performed using three custom-made siRNAs targeting the B4GALT5 mRNA region (siB4GALT5 #1 sense: GCUGCUGUACUUCGUCUAUTT, siB4GALT5 #1 antisense: AUAGACGAAGUACAGCAGCTT; siB4GALT5 #2 sense: GGAAGCCUUCUGAUUGCAUTT, siB4GALT5 #2 antisense: AUGCAAUCAGAAGGCUUCCTT; and siB4GALT5 #3 sense: CCAGUUUCUUGGAAGGUAUTT, siB4GALT5 #3 antisense: AUACCUUCCAAGAAACUGGTT) and Nc control siRNA (GenePharma, China). Cells in the exponential growth phase were plated in a six-well plate(5 × 105 cells/well) and transfected with siRNA after 24h by Lipofectamine 3000 Transfection Reagent following the manufacturer’s instructions. For validation, we extracted total RNA, on which the described RT-PCR was performed, and Western blotting of the total proteins from transfected cells was performed.

Cell proliferation, invasion and migration assays

Transfected cells and control cells in the logarithmic growth phase were harvested. They were seeded in 96-well plates (1000 cells/well) with five replicates per sample and incubated at 37C, and were quantified every 24h using Cell Counting Kit-8 (Dojindo, Kumamoto, Japan). The absorbance was measured at 450 nm for proliferation analysis. A transwell assay was used to evaluate cell invasion. A total of 5 × 104 cells were seeded in the upper chamber with serum-free culture medium (200 μl) after the upper chamber was coated with Matrigel (BD Biosciences, San Jose, CA), and the lower chamber was filled with 20% FBS medium. After culturing for 24h, the cells were fixed with 4% paraformaldehyde and stained with crystal violet for 15 min.The migrated cells were counted in random fields. A scratch assay was used to analyze the migratory ability of cells. After 24h of transfection, cells were seeded in a 6-well plate at 1 × 106 cells/well. When the cells were approximately 95% confluent, they were vertically and linearly scratched with a 20 μl micropipette. After washing to remove the nonadherent cells, the cells were cultured with serum-free culture medium and imaged at 0 h, 24 h and 48 h. The cell migration distance was measured using Image-Pro Plus Analysis software.

Statistical analysis

Survival curves were plotted by the Kaplan-Meier (KM) method,univariate and multivariate analyses were performed using the Cox proportional hazards model. We used the Chi-square/Fisher’s exact test for categorical variables and the two-sample t-test/Wilcoxon rank-sum test for continuous variables to observe the statistical significance. Unless otherwise specified, the 2-tailed p-value less than 0.05 was considered statistically significant.

Availability of data and material

The datasets generated and analysed during the current study are available in the GitHub repository at https://github.com/hliu2016/B4GALT5.

References

  1. Kucukoglu O, Sowa J-P, Mazzolini GD, Syn W-K, Canbay A. Hepatokines and adipokines in nash-related hepatocellular carcinoma. J Hepatol. 2021; 74(2):442–57.

    CAS  PubMed  Article  Google Scholar 

  2. Bray F, Ferlay J, Soerjomataram I, Siegel R, Torre LA, Jemal A. Global cancer statistics 2018: Globocan estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018; 68(6):394–424.

    PubMed  Article  Google Scholar 

  3. Tanaka M, Katayama F, Kato H, Tanaka H, Wang J, Qiao YL, Inoue M. Hepatitis b and c virus infection and hepatocellular carcinoma in china: a review of epidemiology and control measures. J Epidemiol. 2011; 21(6):401–16.

    PubMed  PubMed Central  Article  Google Scholar 

  4. Rimassa L, Personeni N, Czauderna C, Foerster F, Galle P. Systemic treatment of hcc in special populations. J Hepatol. 2021; 74(4):931–43.

    PubMed  Article  Google Scholar 

  5. Fuster MM, Esko JD. The sweet and sour of cancer: glycans as novel therapeutic targets. Nat Rev Cancer. 2005; 5(7):526–42.

    CAS  PubMed  Article  Google Scholar 

  6. Yates B, Braschi B, Gray KA, Seal RL, Tweedie S, Bruford EA. Genenames. org: the hgnc and vgnc resources in 2017. Nucleic Acids Res. 2017; 45(D1):D619–25.

    CAS  PubMed  Article  Google Scholar 

  7. De Matteis S, Ragusa A, Marisi G, De Domenico S, Casadei Gardini A, Bonafè M, Giudetti AM. Aberrant metabolism in hepatocellular carcinoma provides diagnostic and therapeutic opportunities. Oxidative Med Cell Longev. 2018; 11(4):7512159.

    Google Scholar 

  8. Li S-F, Zhu C-S, Wang Y-M, Xie X-X, Xiao L-L, Zhang Z-C, Tang Q-Q, Li X. Downregulation of β1, 4-galactosyltransferase 5 improves insulin resistance by promoting adipocyte commitment and reducing inflammation. Cell Death Dis. 2018; 9(2):1–14.

    Google Scholar 

  9. Tang W, Li M, Qi X, Li J. β1, 4-galactosyltransferase v modulates breast cancer stem cells through wnt/ β-catenin signaling pathway. Cancer Res Treat: Off J Korean Cancer Assoc. 2020; 52(4):1084.

    CAS  Google Scholar 

  10. Rodeheffer C, Shur BD. Targeted mutations in β1, 4-galactosyltransferase i reveal its multiple cellular functions. Biochim Biophys Acta (BBA) - Gen Subj. 2002; 1573(3):258–70.

    CAS  Article  Google Scholar 

  11. Youakim A, Dubois DH, Shur BD. Localization of the long form of beta-1, 4-galactosyltransferase to the plasma membrane and golgi complex of 3t3 and f9 cells by immunofluorescence confocal microscopy. Proc Natl Acad Sci. 1994; 91(23):10913–17.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. Silsirivanit A. Glycosylation markers in cancer. Adv Clin Chem. 2019; 89:189–213.

    CAS  PubMed  Article  Google Scholar 

  13. Kang X, Wang N, Pei C, Sun L, Sun R, Chen J, Liu Y. Glycan-related gene expression signatures in human metastatic hepatocellular carcinoma cells. Exp Ther Med. 2012; 3(3):415–22.

    CAS  PubMed  Article  Google Scholar 

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

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  15. Dutta R, Mahato RI. Recent advances in hepatocellular carcinoma therapy. Pharmacol Ther. 2017; 173:106–17.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. Ruiz E, Rojas TR, Berrospi F, Chávez I, Luque C, Cano L, Doimi F, Pineau P, Deharo E, Bertani S. Hepatocellular carcinoma surgery outcomes in the developing world: A 20-year retrospective cohort study at the national cancer institute of peru. Heliyon. 2016; 2(1):00052.

    Article  Google Scholar 

  17. De Stefano F, Chacon E, Turcios L, Marti F, Gedaly R. Novel biomarkers in hepatocellular carcinoma. Dig Liver Dis. 2018; 50(11):1115–23.

    CAS  PubMed  Article  Google Scholar 

  18. Sun Y, Ma W, Yang Y, He M, Li A, Bai L, Yu B, Yu Z. Cancer nanotechnology: Enhancing tumor cell response to chemotherapy for hepatocellular carcinoma therapy. Asian J Pharm Sci. 2019; 14(6):581–94.

    PubMed  PubMed Central  Article  Google Scholar 

  19. 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.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. Cooper J, Giancotti FG. Integrin signaling in cancer: mechanotransduction, stemness, epithelial plasticity, and therapeutic resistance. Cancer Cell. 2019; 35(3):347–67.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. Wang Y, Fukuda T, Isaji T, Lu J, Im S, Hang Q, Gu W, Hou S, Ohtsubo K, Gu J. Loss of αl, 6-fucosyltransferase inhibits chemical-induced hepatocellular carcinoma and tumorigenesis by down-regulating several cell signaling pathways. FASEB J. 2015; 29(8):3217–27.

    CAS  PubMed  Article  Google Scholar 

  22. Li S, Liu X-Y, Pan Q, Wu J, Liu Z-H, Wang Y, Liu M, Zhang X-L. Hepatitis c virus-induced fut8 causes 5-fu drug resistance in human hepatoma huh7. 5.1 cells. Viruses. 2019; 11(4):378.

    CAS  PubMed Central  Article  Google Scholar 

  23. Huang H, Liu Y, Yu P, Qu J, Guo Y, Li W, Wang S, Zhang J. Mir-23a transcriptional activated by runx2 increases metastatic potential of mouse hepatoma cell via directly targeting mgat3. Sci Rep. 2018; 8(1):1–11.

    Google Scholar 

  24. Lu J, Isaji T, Im S, Fukuda T, Kameyama A, Gu J. Expression of n-acetylglucosaminyltransferase iii suppresses α2, 3-sialylation, and its distinctive functions in cell migration are attributed to α2, 6-sialylation levels. J Biol Chem. 2016; 291(11):5708–20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  25. Yoshihara T, Satake H, Nishie T, Okino N, Hatta T, Otani H, Naruse C, Suzuki H, Sugihara K, Kamimura E, et al. Lactosylceramide synthases encoded by b4galt5 and 6 genes are pivotal for neuronal generation and myelin formation in mice. PLoS Genet. 2018; 14(8):1007545.

    Article  CAS  Google Scholar 

  26. Danolic D, Heffer M, Wagner J, Skrlec I, Alvir I, Mamic I, Susnjar L, Banovic M, Danolić D, Puljiz M. Role of ganglioside biosynthesis genetic polymorphism in cervical cancer development. J Obstet Gynaecol. 2020; 40(8):1127–32.

    CAS  PubMed  Article  Google Scholar 

  27. Narayan G, Murty VV. Integrative genomic approaches in cervical cancer: implications for molecular pathogenesis. Futur Oncol. 2010; 6(10):1643–52.

    CAS  Article  Google Scholar 

  28. Wang D, Zhu Z-Z, Jiang H, Zhu J, Cong W-M, Wen B-J, He S-Q, Liu S-F. Multiple genes identified as targets for 20q13. 12–13.33 gain contributing to unfavorable clinical outcomes in patients with hepatocellular carcinoma. Hep Intl. 2015; 9(3):438–46.

    Article  Google Scholar 

  29. Parker BL, Thaysen-Andersen M, Fazakerley DJ, Holliday M, Packer NH, James DE. Terminal galactosylation and sialylation switching on membrane glycoproteins upon tnf-alpha-induced insulin resistance in adipocytes. Mol Cell Proteomics. 2016; 15(1):141–53.

    CAS  PubMed  Article  Google Scholar 

  30. Kumagai T, Tanaka M, Yokoyama M, Sato T, Shinkai T, Furukawa K. Early lethality of β-1, 4-galactosyltransferase v-mutant mice by growth retardation. Biochem Biophys Res Commun. 2009; 379(2):456–59.

    CAS  PubMed  Article  Google Scholar 

  31. Sperandio M, Gleissner CA, Ley K. Glycosylation in immune cell trafficking. Immunol Rev. 2009; 230(1):97–113.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  32. Nagata M, Izumi Y, Ishikawa E, Kiyotake R, Doi R, Iwai S, Omahdi Z, Yamaji T, Miyamoto T, Bamba T, et al. Intracellular metabolite β-glucosylceramide is an endogenous mincle ligand possessing immunostimulatory activity. Proc Natl Acad Sci. 2017; 114(16):3285–94.

    Article  CAS  Google Scholar 

  33. Zhang T, de Waard AA, Wuhrer M, Spaapen RM. The role of glycosphingolipids in immune cell functions. Front Immunol. 2019; 10:90.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  34. Yoshizaki F, Nakayama H, Iwahara C, Takamori K, Ogawa H, Iwabuchi K. Role of glycosphingolipid-enriched microdomains in innate immunity: microdomain-dependent phagocytic cell functions. Biochim Biophys Acta (BBA) - Gen Subj. 2008; 1780(3):383–92.

    CAS  Article  Google Scholar 

  35. Iwabuchi K. Gangliosides in the immune system: Role of glycosphingolipids and glycosphingolipid-enriched lipid rafts in immunological functions. Gangliosides. 2018; 1804:83–95.

    CAS  Article  Google Scholar 

  36. Mayo L, Trauger SA, Blain M, Nadeau M, Patel B, Alvarez JI, Mascanfroni ID, Yeste A, Kivisäkk P, Kallas K, et al. B4galt6 regulates astrocyte activation during cns inflammation. Nat Med. 2014; 20(10):1147.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

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

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. Ritchie ME, Belinda P, Wu D, Hu Y, Law CW, Wei S, Smyth GK. limma powers differential expression analyses for rna-sequencing and microarray studies. Nucleic Acids Res. 2015; 43(7):47.

    Article  CAS  Google Scholar 

  39. Li P, Xiu WB, Di KL, Xu C, Wang GM, Xia QY, Xiong Q. Large-scale rna-seq transcriptome analysis of 4043 cancers and 548 normal tissue controls across 12 tcga cancer types. Sci Rep. 2015; 5:13413.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  40. Janet P, Núria Q, àlex B, Jordi DP, Anna BM, Martin B, Ferran S, Furlong LI. Disgenet: a discovery platform for the dynamical exploration of human diseases and their genes. Database. 2015; 4(15):bav028.

    Google Scholar 

  41. Pan JB, Hu SC, Shi D, Cai MC, Li YB, Zou Q, Ji ZL, Arnar P. Pagenbase: A pattern gene database for the global and dynamic understanding of gene function. PLoS ONE. 2013; 8(12):80747.

    Article  CAS  Google Scholar 

  42. Sturm G, Finotello F, List M. Immunedeconv: An r package for unified access to computational methods for estimating immune cell fractions from bulk rna-sequencing data. 2020.

  43. Chen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling tumor infiltrating immune cells with cibersort. Methods Mol Biol. 2018; 1711:243.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. Timer: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017; 77(21):108.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

Not applicable.

Funding

The data analysis and in vitro experiments were supported by National Natural Science Foundation of China (82073339; 62072058), Natural Science Foundation of Jiangsu Province (BK20191157), Funding of CSCO-Linghang (Y-2019AZMS-0002).

Author information

Affiliations

Authors

Contributions

HY did the biological experiments and drafted the manuscript. LZ collected the data and performed the statistical analysis. LH and QW helped to improve the idea and the manuscript. SZQ conceived the main idea and the framework of the manuscript. WY reviewed drafts of the paper. LJD supervised the study and provided funding. All authors read and commented on the manuscript. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yong Wu or Judong Luo.

Ethics declarations

Ethics approval and consent to participate

Approval from the ethical board for this study was not required because of the public nature of all the data. All methods were performed in accordance with the relevant guidelines and regulations.

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.

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, visithttp://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

Verify currency and authenticity via CrossMark

Cite this article

Han, Y., Li, Z., Wu, Q. et al. B4GALT5 high expression associated with poor prognosis of hepatocellular carcinoma. BMC Cancer 22, 392 (2022). https://doi.org/10.1186/s12885-022-09442-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-022-09442-2

Keywords

  • Hepatocellular carcinoma
  • Prognostic biomarker
  • Immune infiltration level
  • Expression level
  • Survival analysis