- Open Access
Integrated analysis reveals critical glycolytic regulators in hepatocellular carcinoma
Cell Communication and Signaling volume 18, Article number: 97 (2020)
Cancer cells primarily utilize aerobic glycolysis for energy production, a phenomenon known as the Warburg effect. Increased aerobic glycolysis supports cancer cell survival and rapid proliferation and predicts a poor prognosis in cancer patients.
Molecular profiles from The Cancer Genome Atlas (TCGA) cohort were used to analyze the prognostic value of glycolysis gene signature in human cancers. Gain- and loss-of-function studies were performed to key drivers implicated in hepatocellular carcinoma (HCC) glycolysis. The molecular mechanisms underlying Osteopontin (OPN)-mediated glycolysis were investigated by real-time qPCR, western blotting, immunohistochemistry, luciferase reporter assay, and xenograft and diethyl-nitrosamine (DEN)-induced HCC mouse models.
Increased glycolysis predicts adverse clinical outcome in many types of human cancers, especially HCC. Then, we identified a handful of differentially expressed genes related to HCC glycolysis. Gain- and loss-of-function studies showed that OPN promotes, while SPP2, LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD inhibit HCC cell glycolysis as revealed by glucose utilization, lactate production, and extracellular acidification ratio. These glycolysis-related genes exhibited significant tumor-promoting or tumor suppressive effect on HCC cells and these effects were glycolysis-dependent. Mechanistically, OPN enhanced HCC glycolysis by activating the αvβ3-NF-κB signaling. Genetic or pharmacological blockade of OPN-αvβ3 axis suppressed HCC glycolysis in xenograft tumor model and hepatocarcinogenesis induced by DEN.
Our findings reveal crucial determinants for controlling the Warburg metabolism in HCC cells and provide a new insight into the oncogenic roles of OPN in HCC.
Hepatocellular carcinoma (HCC), the main type of liver cancer, is an increasingly prevalent clinical problem. HCC is the fifth most common cancer worldwide and the third most frequent cause of cancer-related death worldwide . Many risk factors contributes for HCC initiation, such as infection with either hepatitis B virus or hepatitis C virus, nonalcoholic steatohepatitis, alcoholic cirrhosis, and exposure to environmental toxins [2, 3]. Despite many progresses in the diagnosis and treatment of HCC, the prognosis remains dismal due to high rates of recurrence and intrahepatic metastasis after surgical resection . Therefore, uncovering the mechanisms underlying HCC progression to improve clinical outcomes and to develop better therapeutic strategies is of paramount importance .
Accumulating evidences have shown that reprogrammed energy metabolism is a hallmark of tumor cells and a critical contributor to tumor development [6, 7]. Different from most normal cells, cancer cells metabolize glucose to lactate even in the presence of sufficient oxygen, a phenomenon termed aerobic glycolysis, also known as “Warburg effect” . In the past decade, there is a considerable resurgence of interest in the oncogenic role of aerobic glycolysis in cancers. The aerobic glycolysis is characterized by a much higher rate of glucose uptake, consumption and lactate release in cancer cells. Increased glycolysis facilitates cancer cells to rapid utilization of glucose to produce abundant ATP. Through the pentose phosphate pathway, glucose metabolism provides cancer cells with abundant cellular buildings for biosynthetic pathways, including nucleotides, lipids, and nonessential amino acids. Therefore, the glycolysis benefits both bioenergetics and biosynthesis necessary for cancer cell survival and rapid proliferation [8, 9]. Moreover, glycolysis-derived lactate can lead to an acidic tumor microenvironment, which is profoundly implicated in tumor progression by modulating tumor metastasis and immune response . Thus, targeting tumor glycolysis holds distinct promise for development of therapeutic strategies for cancer patients . However, the molecular mechanism for the Warburg effect in HCC is far from explored.
Increased glycolysis is commonly viewed as a consequence of oncogenic event that drives the aggressive phenotypes of cancer cells. During liver malignant transformation, a series of genetic alterations will be selected, including activation of oncogenes and inactivation of tumor suppressors. For example, highly expressed PARP14 can regulate PKM2 activity to promote the Warburg effect and cell survival in HCC ; phosphorylation of ATG4B can impair mitochondrial activity and enhance the Warburg effect ; inactivation of HELLS leads to metabolic reprogramming and reverses the Warburg effect . Notably, loss of FBP1 facilitates aggressive features of HCC cells through the Warburg effect . Several glycolytic enzymes are also upregulated in cirrhotic livers and significantly associated with an increased risk for developing HCC . Previously, Jiang et al. have documented the glycolysis gene expression profiling in HCC and identified six mRNAs (DPYSL4, HOMER1, ABCB6, CENPA, CDK1, and STMN1) significantly associated with overall survival for HCC . Although several modulators involved in HCC glycolysis have been revealed, a comprehensive characterization of glycolysis-related genes is timely needed.
In this study, by leveraging large-scale molecular profiles from The Cancer Genome Atlas (TCGA) cohort, we identify many differentially expressed glycolysis-related genes in HCC. Functional experiments show that Osteopontin (OPN) is a positive regulator, while SPP2, LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD are negative regulator for HCC glycolysis. OPN, encoded by SPP1 (secreted phosphoprotein1) gene, plays a crucial role in HCC initiation and progression and can be acts as a biomarker for HCC [18,19,20,21,22]. Previously, intense investigations have well characterized the oncogenic roles of OPN in HCC. Here, we provide new evidence that OPN is also profoundly implicated in HCC glycolysis by activating the αvβ3-NF-κB signaling.
Materials and methods
The RNA-sequencing data of HCC and corresponding adjacent non-tumor liver tissues were downloaded from The Cancer Genome Atlas (TCGA, https://gdc.cancer.gov/) database. The glycolysis score was calculated based on the mRNA expression value of SLC2A1, HK2, GPI, PFKL, ALDOA, GAPDH, PGK1, PGAM1, ENO1, PKM2, and LDHA. Differentially expressed genes were identified by estimating an exact test P-value.
Cell culture and reagents
The human HCC cell lines (Huh7, HCC-LM3, SMCC-7721, MHCC-97H, Hep3B, and MHCC-97H), NIH3T3 cells, and MEFs were obtained from American Type Culture Collection (ATCC, UK). All cells were cultured with Dulbecco’s modified Eagle’s medium (DMEM, Gibco, Shanghai) supplemented with 10% fetal bovine serum (FBS), 2 mM L-glutamine, 1% penicillin-streptomycin (Sigma, USA), and maintained in a humidified incubator with 5% CO2. All cells were tested negative for cross-contamination of other human cells and mycoplasma contamination. The glycolysis inhibitor 2-Deoxy-D-glucose (2-DG) was obtained from Sigma Aldrich (D6134, St Louis, MO, USA). Galactose was purchased from Sigma-Aldrich (G5388, St Louis, MO, USA). The recombinant OPN protein was purchased from R&D systems (1433-OP-050, Shanghai, China).
The genetic modulation in this study was performed by lentiviral system. The shRNA for OPN and scramble control was obtained from Genepharma Biotechnology (Shanghai, China). The full length sequences of genes for overexpression strategy were synthesized from Sangon Biotech (Shanghai, China). For transfection, cells were grown on 6-well culture plates with 70–80% confluence. Transfection was performed with Lipofectamine 2000 (#11668019, Invitrogen, Carlsbad, CA, USA) according to the procedures of the manufacturer. Lentiviral particles were packaged in 293 T cells using second-generation packaging plasmids psPAX2 (Addgene plasmid 12,260) and pMD2.G (Addgene plasmid 12,259). The cells supernatant was collected 72 h after transfection. The virus was concentrated using ultracentrifuge at 20,000 rpm for 3 h and added to target cells with 5 μg/ml of polybrene (Millipore, USA).
Detection of lactate and glucose level
HCC cells seeded into 6-well plates at a density of 5 × 105 cells were cultured in fresh phenol red-free medium and incubated for 24 h before the culture medium were collected. The lactate and glucose levels in the culture medium were determined by the glucose assay kit (Abcam, Cambridge, MA, USA) and lactate assay kit according to the manufacturer’s instruction (Biovision), and normalized with cell number.
Measurement of extracellular acidification rate
The Seahorse XF96 Flux Analyser (Seahorse Bioscience) was used to monitor extracellular acidification rate (ECAR) according to the manufacturer’s instructions. In brief, HCC cells seeded into 96-well plates at a density of 1.5 × 104 per well were cultured in complete medium. Before experiment, the culture medium was replaced by test buffer (Seahorse Bioscience). For ECAR test, cells were incubated with buffered medium followed by a sequential injection of 10 mM glucose, 1 mM oligomycin and 50 mM 2-DG. The final glycolytic ability was normalized to cell number.
Colony formation assay
For anchorage-dependent growth assay, cells were seeded onto a 6-well plate (1000 cells/well). The culture medium was replaced every 2 days. After 10–12 days, the colonies were fixed with 4% paraformaldehyde and stained with 0.5% crystal violet. The numbers of colonies greater than 100 μm in diameter were counted.
Enzyme linked immunosorbent assay (ELISA) analysis
HCC cell supernatants were collected, cleared by centrifugation and used immediately. The amount of supernatants used was normalized to cell number. The OPN ELISA kits were purchased from R&D systems (DOST00, Shanghai, China) and used according to the manufacturer’s instructions.
Quantitative real-time PCR
Total RNA was extracted from HCC cell lines using Trizol Reagent (Invitrogen, CA, USA). Total RNA (1 μg) was reverse transcribed using PrimeScript® reverse transcriptase Master Mix (TaKaRa, Dalian, China) according to the manufacturer’s instructions. Real-time qPCR was carried out on an Applied Biosystems 7500 apparatus using SYBR-Green Master mix (TaKaRa, Dalian, China). All samples including the template controls were tested in triplicate. The relative expression of target transcripts was normalized to human ACTB gene transcripts. The delta Ct method was used to calculate the relative expression. The primer sequences used in this study were provided in Supplementary Table 1.
Western blotting analysis
Total protein was extracted from HCC cell lines using a total protein extraction buffer (Beyotime, Shanghai, China) with phosphatase and proteinase inhibitors (1 mM EDTA, 1 mM sodium orthovanadate, 10 mM sodium pyrophosphate, 100 mM NaF, 10 mg/ml leupeptin, 10 units/ml aprotinin, 1 mM phenylmethylsulfonyl fluoride) and protein concentration was measured using a BCA Protein Assay Kit (Pierce Biotechnology, USA). Protein samples were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred onto polyvinylidene fluoride membranes. After blocked with 5% non-fat milk/TBST, the membrane was incubated with the primary antibody against OPN (Abcam, ab8448), Akt (Cell Signaling Technology, #4685), p-Akt (Cell Signaling Technology, #4060), Erk1/2 (Cell Signaling Technology, #4695), p-Erk1/2 (Cell Signaling Technology, #9101), p65 (Cell Signaling Technology, #8242), p-p65 (Cell Signaling Technology, #3033), and β-actin (Abcam, ab8227). The signals were acquired after incubation with the species-specific secondary antibodies using an Odyssey Infrared Imaging System (LI-COR, Lincoln, NE).
Paraffin-embedded liver tissues were mounted on silicon-coated slides, dewaxed, and rehydrated, and antigen retrieval was then performed. The slides were incubated with anti-OPN antibody (1:100; Abcam; ab8448)). Quantitation of immunoreactivity was calculated based on the intensity and number of positive cells.
Luciferase reporter assay
HCC cells were transfected with NF-κB-driven luciferase plasmid together with the pRL-TK, which contains the herpes simplex virus thymidine kinase (HSK-TK) promoter to provide low to moderate levels of Renilla luciferase expression by using Lipofectamine 2000 (#11668019, Invitrogen, Carlsbad, CA, USA). After 24 h, the cells were collected and subjected to luciferase activity detection using the Dual Luciferase Reporter Assay System (Promega, USA).
Six-week-old male BALB/c nude mice (Experimental Animal Centre, SIBS) were used in our study. Wild type or modified HCC-LM3 cells were injected subcutaneously into the right flank of these mice to establish xenograft model. Mouse experiments were conducted in accordance with the NIH Guidelines for the Care and Use of Laboratory Animals. The study procedures were approved by the Fifth Affiliated Hospital of Wenzhou Medical University. For pharmacological inhibition experiment, mice were randomly divided into two groups when the tumor volume reached to 200 mm3. Mice were treated with Cilengitide (5 mg/kg) three times a week, while control group was treated with saline with 0.01% DMSO. Tumor volume (mm3) was estimated by the formula: volume = length × width2/2. The tumor volumes data are presented as means ± SD. For DEN-induced HCC model, male C57BL/6-Spp1tm1Blh (OPN−/−) (OPN KO) mice were purchased from Jackson Laboratory (Bar Harbor, ME, USA). The mice at 2 weeks old were injected with 25 mg/kg of DEN (Sigma, Shanghai, China) or vehicle intraperitoneally once to induce hepatic carcinogenesis. Mice in the control groups were not subjected to DEN injection. All mice were sacrificed at 35 weeks after the injection of DEN.
All experiments were conducted independently for three times. Statistical analyses were performed using the SPSS 16.0 (SPSS, Chicago, IL, USA) or Prism 5.0 (GraphPad, La Jolla, CA, USA). The differences between groups were analyzed using Student’s t-test or one-way ANOVA with Dunnett’s multiple comparisons. The Spearman’s rank correlation coefficient and chi-square tests were used to analysis the correlation between indicated parameters. Kaplan-Meier survival analysis was performed to compare patient survival data. P < 0.05 was considered to be statistically significant. *P < 0.05, **P < 0.01, and ***P < 0.001.
Prognostic analyses of glycolysis gene signature in human cancers
Emerging evidences suggest that increased glycolysis contributes to tumor malignant phenotypes and patients’ poor prognosis [9, 23]. More than ten genes encoding the glucose transporter and glycolytic enzymes are directly responsible for the aerobic glycolysis  (Fig. 1a). To identify the correlation between glycolysis and patients’ clinical outcome in human cancers, we calculated a glycolysis score of each tumor sample from The Cancer Genome Atlas (TCGA) cohort based on the expression level of glycolytic components (glucose transporter and glycolytic enzymes). Using the median glycolysis as a cutoff, we analyzed the prognostic value of glycolysis gene signature across 33 types of human cancers (Fig. 1b). As a result, tumor glycolysis was closely correlated with poor prognosis in many types of human cancers, especially for cervical squamous cell carcinoma and endocervical adenocarcinoma (CSES), kidney chromophobe (KICH), liver hepatocellular carcinoma (LIHC), and uveal melanoma (UVM), which had a hazard ratio lager than 2 (Fig. 1c). Specifically, increased glycolysis predicted a better clinical outcome in KIRC patients (Fig. 1b).
Identification of differentially expressed genes for HCC glycolysis
Aerobic glycolysis in cancer cells are largely mediated by activated oncogenes and/or loss of tumor suppressors. In this study, we aimed to uncover the critical modulators responsible for HCC glycolysis. To achieve this, we compared the gene expression pattern among TCGA HCC samples with high and low glycolysis score and identified a handful of differentially expressed genes (DEGs). The top 10 up-regulated and down-regulated DEGs were shown in Fig. 2a. In the top 10 up-regulated genes, S100P and SPP1 were highly expressed in HCC tissues compared with normal liver tissues, whereas PLA2G2A and APOA4 expression were down-regulated in HCC (Fig. 2b); highly expressed SPP1, TMEM92, and EGLN3 predicted a poor prognosis in HCC patients (Fig. 2c and Supplementary Fig. 1). In the top 10 down-regulated genes, SPP2, LECT2, SLC10A1, CYP2A6, CYP3A4, HSD17B13, CYP2A7, and IYD expression were significantly down-regulated, while CYP7A1 expression was up-regulated in tumor tissues in comparison to normal liver tissues (Fig. 2b); notably, SPP2, LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD predicted a better prognosis in HCC patients (Fig. 2c and Supplementary Fig. 1). Combined these data above, we hypothesized that SPP1, SPP2, LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD might be critical modulators for HCC glycolysis. Indeed, there genes had a close correlation with the glycolysis gene signature (Supplementary Fig. 2). Consistent with the results in HCC-LM3, NIH3T3 cells and MEFs expressing OPN exhibited increased glycolysis (Fig. 3e-g), excluding the possibility that the observed differences in glycolysis could be attributed to a cell type-specific phenomenon.
Certification of glycolysis-related genes in HCC cells
Previous studies have well documented diverse aggressive phenotypes induced by OPN (SPP1) in HCC, including tumor growth, metastasis, epithelial-mesenchymal transition, stemness, drug resistance, immunosuppression, and etc. [25,26,27,28]. However, no report is available about its regulatory role in HCC glycolysis. To test whether up-regulated OPN contributes to HCC glycolysis, loss-of-function study was performed in HCC-LM3 cells, which preserve higher OPN at both mRNA and secreted level (Supplementary Fig. 3A and 3B). Two specific shRNA against OPN led to significant reduction in OPN protein level (Fig. 3a). In HCC-LM3 cells, genetic silencing of OPN significantly inhibited glucose utilization and lactate production (Fig. 3b). In addition, treatment with recombinant OPN protein increased glucose utilization and lactate production in Huh7 and Hep3B cells in a dose-dependent manner (Supplementary Fig. 3C and D). To further confirm this observation, Seahorse Extracellular Flux analysis was performed to determine the cellular bioenergetic activity. As a result, the extracellular acidification ratio (ECAR) was markedly reduced by OPN knockdown (Fig. 3c). Moreover, blocking of OPN with its neutralizing antibody also attenuated the glycolytic ability of HCC-LM3 cells as revealed by reduced glucose uptake and lactate production (Fig. 3d). Consistent with the results in HCC-LM3, NIH3T3 cells and MEFs expressing OPN exhibited increased glycolysis (Fig. 3e-g), excluding the possibility that the observed differences in glycolysis could be attributed to a cell type-specific phenomenon. Similarly, we evaluated the potential regulatory role of SPP2, LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD in HCC glycolysis by gain-of-function strategy. The overexpression efficiency of indicated genes was verified by western blotting (Fig. 4a) and real-time qPCR (Fig. 4b), respectively. As displayed in Fig. 4c-e, overexpression of SPP2, LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD differentially suppressed HCC glycolytic capacity as demonstrated by reduced glucose uptake, lactate release, and ECAR.
Effects of glycolysis-related genes on HCC tumor growth
Warburg effect can provide abundant cellular blocks for cancer cell rapid proliferation . Therefore, we investigated whether stimulating or inhibiting the Warburg effect is an important mechanism regulating HCC tumorigenesis by genetic modulation of OPN, SPP2, LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD. The colony formation assay showed that OPN knockdown led to drastic reduction in HCC-LM3 cell anchorage-dependent growth (Fig. 5a), whereas recombinant OPN protein clearly facilitated Hep3B cell anchorage-dependent proliferation (Fig. 5b). Expectedly, overexpression of SPP2, LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD diversely suppressed Huh7 cell anchorage-dependent growth (Fig. 5c). To further confirm this observation, we added 5 mM 2-Deoxy-D-glucose (2-DG) in the culture medium to inhibit glycolysis. The results showed that 2-DG greatly inhibited anchorage-dependent growth of HCC-LM3 and Huh7 cells (Fig. 5d). Notably, 2-DG largely compromised the effects induced by OPN recombinant protein or LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD knockdown (Fig. 5d). Moreover, we replaced the glucose in the culture medium with galactose, which occurs at a much lower rate than glucose entry into glycolysis. Likewise, galactose mirrored the effect induced by 2-DG (Fig. 5e). Collectively, these data above strongly suggest that OPN, SPP2, LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD in regulating the aerobic glycolysis affect HCC tumorigenesis.
OPN promotes HCC glycolysis by modulating αvβ3-NF-κB signaling
Next, we pursued the mechanism by which OPN promotes HCC glycolysis. Integrin receptor αvβ3 is well documented as a critical receptor for OPN-mediated oncogenic activities in HCC . Indeed, pharmacological inhibition of αvβ3 with Cilengitide significantly inhibited HCC-LM3 cell glucose utilization, lactate production, and ECAR (Fig. 6a). Upon binding to integrins and CD44 receptor family, OPN activates various signaling cascade, including PI3K/Akt, and MAPK, and nuclear factor kappa B (NF-κB) in HCC. Western blotting showed that OPN knockdown drastically inhibited the phosphorylation levels of p65, while the phosphorylation levels of Akt and Erk1/2 were slightly reduced. Consistently, OPN treatment significantly increased phosphorylation levels of p65 in Hep3B cells, indicating the prominent regulatory role of OPN in NF-κB activity (Fig. 6b). Interestingly, small molecule inhibitors against NF-κB (Bay11–7082) but not PI3K/Akt (LY294002) and MAPK (U0126) largely abrogated the glycolysis induced by OPN (Fig. 6c). By luciferase reporter assay, we found that OPN knockdown inhibited, whereas recombinant OPN protein increased NF-κB-responsive reporter activity in HCC cells (Fig. 6d). Notably, ectopic expression of constitutively active I-kappa B kinase complex β (CA-IKKβ), at least, to some extent, rescued the reduced glycolysis induced by OPN knockdown (Fig. 6e). Taken together, these data provide mechanism support for the regulatory role of OPN in HCC glycolysis by activating NF-κB signaling.
Inhibition of OPN-αvβ3 axis suppresses HCC tumor growth and glycolyis
Furthermore, we evaluated the therapeutic potential of targeting OPN-αvβ3 axis in HCC. Consistent with previous report, tumor growth of xenografts formed from sh-OPN HCC-LM3 cells was much slower compared with that formed from sh-Ctrl HCC-LM3 cells (Fig. 7a). Moreover, pharmacological blockade of αvβ3 with Cilengitide remarkably inhibited tumor growth of HCC-LM3 xenografts (Fig. 7b). By measuring the lactate level in xenograft tumor tissues, we found that both OPN knockdown and αvβ3 inhibition reduced lactate concentration in vivo (Fig. 7c). Moreover, real-time qPCR analysis showed that OPN knockdown significantly down-regulated the expression of glucose transporter (SLC2A1) and glycolytic enzymes (HK2, PFKL, PKM2, and LDHA) (Fig. 7d). Similar result was also noticed in the pharmacological inhibition studies (Fig. 7d). Moreover, we generated DEN-induced HCC model in OPN-KO and wide type (WT) mice. As shown in Fig. 7e, the nodules from the OPN KO mice at 35 weeks were composed of well-differentiated tumor cells or vacuolated cells that formed trabeculae or nests, whereas the nodules from the WT mice showed a sessile and solid growth pattern. Real-time qPCR resulted revealed that expression of the glycolytic enzymes was markedly reduced in liver tissues of OPN-KO mice compared to that in WT mice (Fig. 7f). Collectively, these data indicate that targeting OPN-αvβ3 axis is sufficient to block HCC tumor growth and glycolysis.
Expression pattern of OPN in clinical samples
Consistent with previous reports, OPN is highly expressed in human HCC samples and correlated a poor overall survival rate and disease-free survival rate in HCC patients (Supplementary Fig. 4). In this study, we further investigated the correlation between OPN and the Warburg metabolism in clinical HCC samples. In a cohort of 60 HCC patients, the glycolytic enzymes were highly expressed in OPN-high samples compared to OPN-low HCC samples (Fig. 8a). Moreover, in a cohort of 20 HCC patients who received preoperative 18F-FDG PET/CT, we found that the SUVmax was much higher in the group with high expression of OPN (Fig. 8b). This result further supports the critical regulatory of the OPN in HCC glycolysis.
Many molecular targeted therapies and diagnostic methods have been developed for the treatment of HCC; however, the clinical outcome of HCC is still unsatisfactory. Recent data show that increased glycolysis is closely associated with poor prognosis in cancer patients, and targeting cancer glycolysis metabolism represents an attractive strategy for cancer therapy. Therefore, it is very important to uncover critical modulators that regulate HCC glycolysis. In this study, we found that tumor samples with high glycolysis gene signature had a significant reduced overall survival time in many cancer types, especially in HCC. By leveraged TCGA HCC RNA-seq data, we identified seven differentially expressed genes (OPN, SPP2, LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD) associated with HCC glycolysis. In addition, these genes had a close correlation with the glycolysis score and patients’ prognosis in HCC. Notably, none of these genes have been reported to play a role in the regulation in cancer metabolism. Our experimental data demonstrated glycolytic alterations induced by these genes are linked to tumor growth. The functional validation of glycolysis-related genes provides candidate therapeutic targets and associated biomarkers for this deadly disease.
OPN is a multifunctional secreted phosphoprotein with intense studies in liver disease, including acute liver injury, autoimmune and viral hepatitis, alcoholic liver diseases, non-alcoholic fatty fiver diseases, liver fibrosis and HCC . Importantly, OPN is involved in numerous oncogenic activities in HCC, including the cell survival, proliferation, stemness, angiogenesis, invasion, and metastasis [21, 28, 30, 31]. High OPN level is significantly associated with poor overall survival and OPN expression is positively associated with stage and tumor size . OPN is also markedly increased in the plasma of HCC patients and emerged as a diagnostic biomarker that improves AFP performance in HCC surveillance among patients with HBV or HCV-related cirrhosis [33, 34]. In line with previous report, we found that OPN is overexpressed in HCC tissues and predicts a poor prognosis. By loss-of-function study, we certified that OPN is an important glycolysis regulator in HCC. Consistently, blockade of OPN with neutralizing antibody significantly suppressed HCC cell glycolytic activity. Interesting, tumor glycolysis-induced acidic microenvironment can activate hepatic stellate cells, which produce OPN under acidic condition and in turn promote the migration of HCC cells . Previously, Shi et al. showed that OPN-a isoform increases the cellular glucose levels, and OPN-c utilizes this glucose to generate energy in breast cancer cells . These data together suggest the regulatory role of OPN in tumor glycolysis. In HCC, OPN acts through multiple adhesion receptor binding motifs including integrins and CD44 receptors to promote tumor progression and metastasis [29, 37]. In this study, we revealed that αvβ3-NF-κB signaling is responsible for OPN-induced glcolysis in HCC. Given other integrins and CD44 receptor are also profoundly implicated in reprogrammed energy metabolism of cancer cells, we cannot fully rule out their contributions in HCC [38, 39]. OPN has been well studied as a potential therapeutic target in HCC treatment . Consistent observation with in HCC cells with lower OPN level, glycolytic capacity was also increased by overexpression of OPN in normal MEF and NIH3T3 cells, suggesting the regulatory role of OPN on glycolysis is not cell type-specific phenomenon. Given that OPN is highly expressed in the HCC tumor tissues and can be merely detected in normal liver, there should be a therapeutic window for targeting OPN for HCC treatment. Therefore, our study further broadens the mechanism regarding the therapeutic potential by targeting cancer glucose metabolism.
Apart from OPN, the role of LECT2, SLC10A1, CYP3A4, and HSD17B13 are also characterized in HCC. LECT2 (Leukocyte cell-derived chemotoxin 2) is originally identified as a chemotactic factor for neutrophils and stimulates the growth of chondrocytes and osteoblasts and acts as a tumor suppressor in HCC by multiple mechanisms, such as inactivation of MET, controlling inflammatory monocytes, and inhibition of VEGF165/VEGFR2-dependent signaling [41,42,43]. Down-regulated SLC10A1 is correlated with poor post-surgery survival rate and larger tumor tissue mass in HCC patients and ectopic expression of NTCP significantly suppresses HCC cell proliferation . Down-regulation of CYP3A4 gene is an independent predictor for HCC early recurrence . HSD17B13 expression in peritumoral tissues is associated with worse recurrence free survival and overall survival of HCC patients and the HSD17B13 rs72613567 loss of function variant is protective of HCC development in patients with alcoholic liver disease [46, 47]. The expression pattern and functional role of IYD in cancers are not known. In this study, we for the first time showed that LECT2, SLC10A1, CYP3A4, and HSD17B13 can suppress HCC tumor growth by inhibiting Warburg effect. However, the detailed mechanism by which these genes regulate glycolysis warrants further investigations.
Our comprehensive analysis provides essential insights into HCC glycolysis and uncovers several critical glycolysis-related genes. OPN promotes HCC glycolysis through activation of αvβ3-NF-κB signaling. LECT2, SLC10A1, CYP3A4, HSD17B13, and IYD can suppress HCC cell proliferation by down-regulating glycolytic flux. Therefore, our results open new avenues for designing potential therapeutic strategies for human HCC.
Availability of data and materials
Source data and reagents are available from the corresponding author upon reasonable request.
The Cancer Genome Atlas (TCGA)
American Type Culture Collection
Dulbecco’s modified Eagle’s medium
Fetal bovine serum
Extracellular acidification rate
Liver hepatocellular carcinoma
Differentially expressed genes
Leukocyte cell-derived chemotoxin 2
Forner A, Reig M, Bruix J. Hepatocellular carcinoma. Lancet. 2018;391:1301–14.
Xu XF, Xing H, Han J, Li ZL, Lau WY, Zhou YH, Gu WM, et al. Risk factors, patterns, and outcomes of late recurrence after liver resection for hepatocellular carcinoma: a multicenter study from China. JAMA Surg. 2018;154:209–17.
Fujiwara N, Friedman SL, Goossens N, Hoshida Y. Risk factors and prevention of hepatocellular carcinoma in the era of precision medicine. J Hepatol. 2018;68:526–49.
Kulik L, El-Serag HB. Epidemiology and Management of Hepatocellular Carcinoma. Gastroenterology. 2019;156:477–491 e471.
Braillon A. Hepatocellular carcinoma screening: seeking robust evidence. Gastroenterology. 2019;156:288–9.
Cairns RA, Harris IS, Mak TW. Regulation of cancer cell metabolism. Nat Rev Cancer. 2011;11:85–95.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144:646–74.
Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009;324:1029–33.
Gatenby RA, Gillies RJ. Why do cancers have high aerobic glycolysis? Nat Rev Cancer. 2004;4:891–9.
Corbet C, Feron O. Tumour acidosis: from the passenger to the driver's seat. Nat Rev Cancer. 2017;17:577–93.
Martinez-Outschoorn UE, Peiris-Pages M, Pestell RG, Sotgia F, Lisanti MP. Cancer metabolism: a therapeutic perspective. Nat Rev Clin Oncol. 2017;14:113.
Iansante V, Choy PM, Fung SW, Liu Y, Chai JG, Dyson J, Del Rio A, et al. PARP14 promotes the Warburg effect in hepatocellular carcinoma by inhibiting JNK1-dependent PKM2 phosphorylation and activation. Nat Commun. 2015;6:7882.
Ni Z, He J, Wu Y, Hu C, Dai X, Yan X, Li B, et al. AKT-mediated phosphorylation of ATG4B impairs mitochondrial activity and enhances the Warburg effect in hepatocellular carcinoma cells. Autophagy. 2018;14:685–701.
Law CT, Wei L, Tsang FH, Chan CY, Xu IM, Lai RK, Ho DW, et al. HELLS regulates chromatin remodeling and epigenetic silencing of multiple tumor suppressor genes in human hepatocellular carcinoma. Hepatology. 2019;69:2013–30.
Yang J, Wang C, Zhao F, Luo X, Qin M, Arunachalam E, Ge Z, et al. Loss of FBP1 facilitates aggressive features of hepatocellular carcinoma cells through the Warburg effect. Carcinogenesis. 2017;38:134–43.
Lee NCW, Carella MA, Papa S, Bubici C. High expression of glycolytic genes in cirrhosis correlates with the risk of developing liver cancer. Front Cell Dev Biol. 2018;6:138.
Jiang L, Zhao L, Bi J, Guan Q, Qi A, Wei Q, He M, et al. Glycolysis gene expression profilings screen for prognostic risk signature of hepatocellular carcinoma. Aging (Albany NY). 2019;11:10861–82.
Zhu Y, Yang J, Xu D, Gao XM, Zhang Z, Hsu JL, Li CW, et al. Disruption of tumour-associated macrophage trafficking by the osteopontin-induced colony-stimulating factor-1 signalling sensitises hepatocellular carcinoma to anti-PD-L1 blockade. Gut. 2019;68:1653–66.
Shang S, Plymoth A, Ge S, Feng Z, Rosen HR, Sangrajrang S, Hainaut P, et al. Identification of osteopontin as a novel marker for early hepatocellular carcinoma. Hepatology. 2012;55:483–90.
Xue YH, Zhang XF, Dong QZ, Sun J, Dai C, Zhou HJ, Ren N, et al. Thrombin is a therapeutic target for metastatic osteopontin-positive hepatocellular carcinoma. Hepatology. 2010;52:2012–22.
Sun BS, Dong QZ, Ye QH, Sun HJ, Jia HL, Zhu XQ, Liu DY, et al. Lentiviral-mediated miRNA against osteopontin suppresses tumor growth and metastasis of human hepatocellular carcinoma. Hepatology. 2008;48:1834–42.
Zhao J, Dong L, Lu B, Wu G, Xu D, Chen J, Li K, et al. Down-regulation of osteopontin suppresses growth and metastasis of hepatocellular carcinoma via induction of apoptosis. Gastroenterology. 2008;135:956–68.
Kitamura K, Hatano E, Higashi T, Narita M, Seo S, Nakamoto Y, Yamanaka K, et al. Proliferative activity in hepatocellular carcinoma is closely correlated with glucose metabolism but not angiogenesis. J Hepatol. 2011;55:846–57.
Li L, Liang Y, Kang L, Liu Y, Gao S, Chen S, Li Y, et al. Transcriptional regulation of the Warburg effect in cancer by SIX1. Cancer Cell. 2018;33:368–85 e367.
Wen Y, Jeong S, Xia Q, Kong X. Role of osteopontin in liver diseases. Int J Biol Sci. 2016;12:1121–8.
Dong Q, Zhu X, Dai C, Zhang X, Gao X, Wei J, Sheng Y, et al. Osteopontin promotes epithelial-mesenchymal transition of hepatocellular carcinoma through regulating vimentin. Oncotarget. 2016;7:12997–3012.
Cao L, Fan X, Jing W, Liang Y, Chen R, Liu Y, Zhu M, et al. Osteopontin promotes a cancer stem cell-like phenotype in hepatocellular carcinoma cells via an integrin-NF-kappaB-HIF-1alpha pathway. Oncotarget. 2015;6:6627–40.
Liu G, Fan X, Tang M, Chen R, Wang H, Jia R, Zhou X, et al. Osteopontin induces autophagy to promote chemo-resistance in human hepatocellular carcinoma cells. Cancer Lett. 2016;383:171–82.
Rangaswami H, Bulbule A, Kundu GC. Osteopontin: role in cell signaling and cancer progression. Trends Cell Biol. 2006;16:79–87.
Takafuji V, Forgues M, Unsworth E, Goldsmith P, Wang XW. An osteopontin fragment is essential for tumor cell invasion in hepatocellular carcinoma. Oncogene. 2007;26:6361–71.
Gao X, Sheng Y, Yang J, Wang C, Zhang R, Zhu Y, Zhang Z, et al. Osteopontin alters DNA methylation through up-regulating DNMT1 and sensitizes CD133+/CD44+ cancer stem cells to 5 azacytidine in hepatocellular carcinoma. J Exp Clin Cancer Res. 2018;37:179.
Sun T, Li P, Sun D, Bu Q, Li G. Prognostic value of osteopontin in patients with hepatocellular carcinoma: a systematic review and meta-analysis. Medicine (Baltimore). 2018;97:e12954.
Salem M, Atti SA, Raziky ME, Darweesh SK, Sharkawy ME. Clinical significance of plasma Osteopontin level as a biomarker of hepatocellular carcinoma. Gastroenterol Res. 2013;6:191–9.
Kim J, Ki SS, Lee SD, Han CJ, Kim YC, Park SH, Cho SY, et al. Elevated plasma osteopontin levels in patients with hepatocellular carcinoma. Am J Gastroenterol. 2006;101:2051–9.
Song J, Ge Z, Yang X, Luo Q, Wang C, You H, Ge T, et al. Hepatic stellate cells activated by acidic tumor microenvironment promote the metastasis of hepatocellular carcinoma via osteopontin. Cancer Lett. 2015;356:713–20.
Shi Z, Wang B, Chihanga T, Kennedy MA, Weber GF. Energy metabolism during anchorage-independence. Induction by osteopontin-c. PLoS One. 2014;9:e105675.
Shirasaki T, Honda M, Yamashita T, Nio K, Shimakami T, Shimizu R, Nakasyo S, et al. The osteopontin-CD44 axis in hepatic cancer stem cells regulates IFN signaling and HCV replication. Sci Rep. 2018;8:13143.
Ata R, Antonescu CN. Integrins and cell metabolism: an intimate relationship impacting cancer. Int J Mol Sci. 2017;18:189.
Nam K, Oh S, Shin I. Ablation of CD44 induces glycolysis-to-oxidative phosphorylation transition via modulation of the c-Src-Akt-LKB1-AMPKalpha pathway. Biochem J. 2016;473:3013–30.
Cheng F, Ma C, Wang X, Zhai C, Wang G, Xu X, Mu J, et al. Effect of traditional Chinese medicine formula Sinisan on chronic restraint stress-induced nonalcoholic fatty liver disease: a rat study. BMC Complement Altern Med. 2017;17:203.
Chen CK, Yu WH, Cheng TY, Chen MW, Su CY, Yang YC, Kuo TC, et al. Inhibition of VEGF165/VEGFR2-dependent signaling by LECT2 suppresses hepatocellular carcinoma angiogenesis. Sci Rep. 2016;6:31398.
L'Hermitte A, Pham S, Cadoux M, Couchy G, Caruso S, Anson M, Crain-Denoyelle AM, et al. Lect2 controls inflammatory monocytes to constrain the growth and progression of hepatocellular carcinoma. Hepatology. 2019;69:160–78.
Chen CK, Yang CY, Hua KT, Ho MC, Johansson G, Jeng YM, Chen CN, et al. Leukocyte cell-derived chemotaxin 2 antagonizes MET receptor activation to suppress hepatocellular carcinoma vascular invasion by protein tyrosine phosphatase 1B recruitment. Hepatology. 2014;59:974–85.
Kang J, Wang J, Cheng J, Cao Z, Chen R, Li H, Liu S, et al. Down-regulation of NTCP expression by cyclin D1 in hepatitis B virus-related hepatocellular carcinoma has clinical significance. Oncotarget. 2017;8:56041–50.
Ashida R, Okamura Y, Ohshima K, Kakuda Y, Uesaka K, Sugiura T, Ito T, et al. CYP3A4 gene is a novel biomarker for predicting a poor prognosis in hepatocellular carcinoma. Cancer Genomics Proteomics. 2017;14:445–53.
Yang J, Trepo E, Nahon P, Cao Q, Moreno C, Letouze E, Imbeaud S, et al. A 17-Beta-Hydroxysteroid dehydrogenase 13 variant protects from hepatocellular carcinoma development in alcoholic liver disease. Hepatology. 2019;70:231–40.
Chen J, Zhuo JY, Yang F, Liu ZK, Zhou L, Xie HY, Xu X, et al. 17-beta-hydroxysteroid dehydrogenase 13 inhibits the progression and recurrence of hepatocellular carcinoma. Hepatobiliary Pancreat Dis Int. 2018;17:220–6.
This study was supported by National Natural Science Foundation of China (No.81803778), The Key Research and development Project of Zhejiang Province (No. 2018C0302), The Medical and Health Care Key Project of Zhejiang Province (No.WKJ-ZJ-1629), Natural Science Foundation of Zhejiang Province (No. LQ17H180001), The Public Welfare Project of Zhejiang Province (No.2016C37101, 2017C33216 and LGF18H160035), The Science and Technology Development Project of Lishui City (No. 2017ZDXK07), and The Analysis Test Project of Zhejiang Province (No.2018C37039).
Ethics approval and consent to participate
The study has been examined and certified by the Fifth Affiliated Hospital of Wenzhou Medical University, and informed consent was obtained from all participants included in the study, in agreement with institutional guidelines.
Consent for publication
The Central Hospital of Zhejiang Lishui.
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 glycolysis-related genes in HCC. Supplementary Fig. 2 The correlation between glycolysis-related genes and the glycolysis signature. Supplementary Fig. 3 Expression pattern and roles of OPN in HCC cells. Supplementary Fig. 4 Expression pattern and prognostic value of OPN in HCC tissues. Supplementary Table 1: Primers used in this study. Supplementary Table 2: Group information.
About this article
Cite this article
Lu, C., Fang, S., Weng, Q. et al. Integrated analysis reveals critical glycolytic regulators in hepatocellular carcinoma. Cell Commun Signal 18, 97 (2020). https://doi.org/10.1186/s12964-020-00539-4
- Liver cancer
- Tumor metabolism
- Energy metabolism