Loss of MiR-192-5p Initiates a Hyperglycolysis and Stemness Positive Feedback in Hepatocellular Carcinoma.

Loss of MiR-192-5p Initiates a Hyperglycolysis and Stemness Positive Feedback in Hepatocellular Carcinoma. | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Loss of MiR-192-5p Initiates a Hyperglycolysis and Stemness Positive Feedback in Hepatocellular Carcinoma. Yuanzhuo Gu, Fubo Ji, Niya Liu, Yongzhi Zhao, Xiyang Wei, Shiyuan Hu, and 8 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-56821/v3 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 30 Nov, 2020 Read the published version in Journal of Experimental & Clinical Cancer Research → Version 3 posted 4 You are reading this latest preprint version Show more versions Abstract Background: Emerging studies revealed that cancer stem cells (CSCs) possessed peculiar metabolic properties, which however remained largely unknown in hepatocellular carcinoma (HCC). Genetic silencing of liver-abundant miR-192-5p was a key feature for multiple groups of CSC-positive HCCs. We thus aimed to investigate essential metabolic features of hepatic CSCs via using HCCs with miR-192-5p silencing as a model. Methods: Datasets from two independent HCC cohorts were used. Data integration analyses of miR-192-5p with metabolome and mRNA transcriptome data in HCC Cohort 1 were performed to investigate miR-192-5p related metabolic features, which was further validated in Cohort 2. Cellular and molecular assays were performed to examine whether and how miR-192-5p regulated the identified metabolic features . Co-culture systems consisting of HCC cells and LX2 (human hepatic stellate cell line) or THP1 (human monocyte cell line) were established to explore effects of the identified metabolic properties on stemness features of HCC cells via interacting with co-cultured non-tumor cells. Results: High levels of glycolysis-related metabolites and genes were present in HCCs with low miR-192-5p and CSC-positive HCCs in two independent HCC cohorts. miR-192-5p knockout cells displayed CSC features and miR-192-5p loss led to an enhanced glycolytic phenotype via upregulating three bona fide targets, GLUT1 and PFKFB3 (two glycolytic enzymes) and c-Myc (regulating glycolytic genes’ expression). Meanwhile, c-Myc suppressed miR-192-5p transcription, ensuring a low-miR-192-5p/high-c-Myc loop to maintain hyperglycolysis. Moreover, over-produced lactic acid from hyperglycolytic HCC cells stimulated the ERK phosphorylation of co-cultured LX2 and THP1 non-tumor cells partially via NDRG3 and MCT1, which in turn promoted cell malignancy and stemness of HCC cells. Consistently, HCC patients with low level of miR-192-5p in their tumor tissues and high level of NDRG3 or MCT1 in their non-tumor tissues had the shortest overall survival. Conclusions: In CSC-positive HCCs, miR-192-5p loss enhanced glycolysis and over produced lactate might further increase HCC malignant features via interacting with environmental non-tumor cells. Cancer Biology Hepatocellular carcinoma Cancer stem cell miR-192-5p Glycolysis c-Myc Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Background Abnormal cancer metabolism is one of the ten cancer hallmarks. Studies of cancer metabolism have revealed the important roles of metabolic reprogramming in cancer cells for proliferation, metastasis and drug resistance[ 1 , 2 ]. In tumors, a small population of cancer cells exhibit a high capacity of self-renewal and tumor initiation in NOD/SCID mice, which are referred to as cancer stem cells (CSCs). Currently, CSCs are considered to be responsible for tumor initiation, therapy resistance as well as recurrence, and rising attentions have been placed on CSC metabolism [ 3-5 ]. Recent studies indicate that CSCs have different metabolic properties when compared to the tumor bulk and that metabolic plasticity is dependent on the intrinsic demand of nutrients as well as the surrounding environment [ 6-9 ]. CSCs exhibit distinct metabolic phenotypes that include low mitochondrial respiration, high glycolytic activity, and high fatty acid oxidation based on tumor types and CSC biomarkers as well as isolation methods [ 6 , 10 , 11 ]. Exploiting metabolic vulnerability of CSCs may provide new effective cancer therapies to diminish tumor recurrence and metastasis. Primary liver cancer is the fourth lethal neoplasm worldwide, 90% of which are hepatocellular carcinoma (HCC) [ 12 ]. Hepatic CSCs are considered as one of the determining factors for HCC carcinogenesis and recurrence and multiple hepatic CSC biomarkers such as EpCAM, CD133, CD90, CD44, and CD24, have been used to enrich tumorigenic CSCs in both HCC cell lines and primary HCC tissues [ 13-19 ] . Although gene signatures and regulatory pathways of CSC + HCCs were thoroughly investigated, their metabolic properties remained largely unknown except for few recent studies. Three studies revealed a high glycolysis rate and a low oxygen consumption rate (OCR) in CD133 + hepatic CSCs[ 20 ], the importance of fatty acid oxidation in CD133 + CSCs[ 21 ], and increased lipid metabolites associated with stearoyl-CoA-desaturase in EpCAM + HCCs[ 22 ]. Hence, it is essential to systematically investigate the key metabolic features shared among different CSC biomarker positive HCCs. MiR-192-5p is the 2 nd most abundant miRNA in the liver [ 23 ]. Its genetic silencing frequently occurs in many groups of hepatic CSC marker-positive HCCs such as EpCAM + , CD90 + , CD133 + , CD44 + and CD24 + HCCs, as well as pluripotency marker-positive HCCs[ 23 ]. HCC cells with miR-192-5p loss display a highly invasive phenotype and essential CSC features partially mediated via the p53/miR-192-5p/PABPC4 pathway [ 23-25 ]. Thus, investigating key metabolic features in HCCs with miR-192-5p loss might allow us to explore the essential metabolic properties among multiple different groups of CSC + HCCs. Methods Omics dataset A total of five datasets from two HCC cohorts with 548 cases were used. In Cohort 1, there were a total of 176 Asian HCC cases, among which 22 cases had available metabolomics data and all cases had miRNA transcriptome (GSE6857) and mRNA transcriptome data (GSE14520) [ 22 , 23 , 26-29 ] (Fig. 1a). In Cohort 2, there were 372 HCC cases with different races. Their miRNA and mRNA sequencing data were downloaded from The Cancer Genome Atlas (TCGA) and available TP53 mutation status and MYC amplification information in 240 HCC cases were collected from www.cbioportal.org/index.do . Cell lines and miR-192-5p knockout HCC cells Human liver cancer cell lines including HLF, HLE, Huh7 and HepG2 cells; human embryonic kidney HEK293T cells; human hepatic stellate cell (HSC) line LX2 cells; human hepatocyte line HL7702 cells; and human leukemic monocyte cell line THP1 (with macrophage activation) were used in this study as described before[ 23 , 30 ]. HLF, HLE and Huh7 were originally from Japanese Collection of Research Biosources Cell Bank (JCRB). HepG2, 293T and THP1 were from American Type Culture Collection (ATCC), and LX2 and HL7702 were from Chinese Academy of Sciences (Shanghai, China). Cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% FBS and 100 U/mL penicillin–streptomycin, 1% L-glutamine. Cells were grown at 37°C in a 5% CO2 incubator and regularly tested to be free of mycoplasma contamination. MiR-192-5p knockout HCC cells were generated using the CRISPR/Cas9 system. Two single-guide RNAs (sgRNAs) were designed via an online web tool ( http://crispr.mit.edu ) to target the precursor mir-192. The synthesized sequences were cloned into pSpCas9(BB)-2A-Puro (PX459) vector, termed as PX459-sgMiR192. 48 hours after PX459-sgMiR192 transfection, HCC cells underwent puromycin selection (1.5µg/ml) for 5 days. Next, 1000 collected cells were seeded in a 10-cm dish and cultured for 10-12 days to select miR-192-5p knockout single cell clones. Positive clones were identified via PCR examination and sequencing validation. The sequences for guide RNAs and PCR primers are listed in Supplementary Table S1. Plasmids, siRNAs and cell treatments Lentiviral constructs pre-hsa-miR-control (pmiR-control) and pre-hsa-miR-192-5p (pmiR-192) with GFP as well as miRzip-control and miRZip-192-5p (miRZip-192) were obtained from SBI Biosciences and stored in our laboratory. pmiR-control with RFP (pmiR-ctrl/RFP) and pmiR-192 with RFP (pmiR-192/RFP) were constructed by replacing GFP with Ds-Red between the HindIII and NotI sites. Lentiviruses were packaged with plasmids psPAX2 and pMD2.G (Addgene) in 293T cells. For infection, 5 MOI of each lentivirus was used for all our studies. pMiR-Report-control (Luc-Ctrl), pmiR-192 reporter (192pos reporter), and miRZip-192 reporter (Zip192 reporter) were constructed before[ 23 ]. The miR-192-5p binding regions in the 3’UTR or coding region of GLUT1 , PFKPB3 , and MYC were inserted into the monoclonal sites (HindIII/SpeI) of the pMiR-Report plasmid to generate the corresponding luciferase reporters. The pT3-EF1α-c-Myc vector was originally modified from pT3-EF1α for c-Myc overexpression[ 31 ]. pGL-miR-192-5p (pGL-192) plasmids were constructed by inserting different lengths of the miR-192 promoter region to the monoclonal site (KpnI/XhoI) of pGL4.20-basic (Promega). PFKFB3 siRNAs, GLUT 1 siRNAs and negative control siRNA were purchased from RiboBio Co, Guangzhuo, China. TP53 siRNAs, MYC siRNAs, NDRG3 siRNAs, and MCT1 siRNAs were purchased from GenePharma Co., Shanghai, China. Lipofectamine 2000 (Invitrogen) reagent and Rfect siRNA Transfection Reagent (BIO-TRAN) were used for the transfections of plasmids and siRNAs, respectively. The detailed information of all primers and sequences is listed in Supplementary Table S1. Nutlin-3a (Selleck Chemicals) was dissolved in DMSO. When indicated, HCC cells were treated for 24 hours with 10 µM of Nutlin-3a. 2-DG (Sigma-Aldrich) was dissolved in ddH 2 O and HCC cells were treated for 72 hours with 5-10 µM of 2-DG. Extracellular acidification rate (ECAR) and oxygen consumption rate (OCR) assays In the XF8 Extracellular Flux Analyzer (Seahorse Bioscience), the ECAR and OCR were measured using the Seahorse XF Glycolysis Stress Test Kit (Agilent Technologies) and Seahorse XF Cell Mito Stress Test Kit (Agilent Technologies), respectively. Experiments were performed according to the manufacturer’s instructions. For the ECAR measurement, cells (1 × 10 4 cells/well) were plated in an XF8 cell-culture microplate for 12 hours and the culture medium was then replaced with XF assay medium supplemented plus 2 mM L-glutamine. After one-hour incubation at 37°C, the ECAR was measured by the sequential addition of glucose (10 mM), oligomycin (3.5 μM), and 2-deoxyglucose (2DG, 80 mM) in an XF8 flux analyzer. For the OCR measurement, cells (1 × 10 4 cells/well) were plated in an XF8 cell-culture microplate for 12 hours and then culture medium was replaced with XF assay medium supplemented with 10 mM D-glucose, 1 mM sodium pyruvate and 2 mM L-glutamine. The OCR was measured by sequential addition of oligomycin (2 μM), FCCP (0.5 μM), and antimycin A (1 μM). Data were analyzed by Seahorse XF Wave software. The results were normalized to the cell number. Co-culture system Two types of co-culture systems were used. In the first, HCC cells infected with corresponding RFP or GFP lentivirus were co-cultured with HL7702 or LX2 cells in the same dish. In this system, cells that displayed red or green under a fluorescence microscope or flow cytometry were regarded as HCC cells. The second co-culture system was based on a chamber system using Polyester (PET) Membrane Tissue Culture Plate Insert with 0.4 μm pores (JET biofil) according to the manufacturer’s protocol. Generally, LX2 or HL7702 were placed in the bottom layer of a 6-well plate while HCC cells were placed in the cell inserts of this 6-well plate. In addition, HCC cells were also co-cultured with THP1, a suspension cell line. RNA isolation, quantitative real-time PCR, and Western blot Total RNA was isolated using TRIzol (Invitrogen) following the manufacturer's instructions, and 1 μg of total RNA was reverse transcribed into cDNA using PrimeScript TM RT Reagent Kit with gDNA Eraser (TaKaRa). Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed with TB Green TM Premix EX Taq TM II (Tli RNaseH Plus) (TaKaRa). The expression of mature miRNAs was measured using TaqMan MiRNA Assays as described previously[ 28 , 32 , 33 ]. RNU6B was used as the reference gene for miRNAs, while 18S was the reference for measured mRNAs. Primers for TACSTD1, CD133, CD90, CD44, CD24, and UGT2B7 were used as before[ 23 ] . Primer sequences for GLUT1 , HK2, PFKP, PFKFB3, ALDOA, ENO2, PKM2, LDHA, MCT4, MCT1, NDRG3 and CYP1A2 are listed in Supplementary Table S1. The experiments were performed in triplicate. For the Western blot assay, cells were lysed and processed as previously described [ 28 , 33 ]. The membranes with transferred protein extracts were incubated with the indicated primary antibodies and secondary antibodies conjugated to horseradish peroxidase for enhanced chemiluminescence detection of the signals (Amersham, Arlington Height, IL). The detailed information of all antibodies is listed in Supplementary Table S2. Non-targeted metabolomics study The non-targeted metabolomics was performed by Metabo-Profile Biotechnology (Shanghai) Co., Ltd. Cell lysates and cell culture medium were collected from cells and used for metabolomics analysis. Metabolites were identified and quantified from gas chromatography/time-of-flight mass spectrometry (GC/TOF-MS) data, which was performed in the workflow of ADAP-GC 2.0. Each experiment was performed in triplicate. 2-NBDG uptake, glucose measurement and lactate measurement Glucose uptake of HCC cells was quantified by flow cytometry using the fluorescent D-glucose derivate 2-(N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino)-2-deoxy-D-glucose (2-NBDG). Briefly, the corresponding cells were cultured with 2-NBDG at 37°C for 1 h. The uptake of 2-NBDG was measured by flow cytometry. A glucose measurement kit (Shanghai Rongsheng Biotech Co., Ltd.) was used to measure the glucose concentration in cell culture medium. Extracellular lactate production was measured using a lactate assay kit (Nanjing Jiancheng Bioengineering Institute). All these assays were performed according to the manufacturer's protocols. Sphere formation assay and flow cytometry analysis Single-cell suspensions of 1000 cells were seeded in 6-well Ultra-Low Attachment Microplates (Corning, Corning, NY) for spheroid assays. The number of spheroids was measured 12 days after seeding. For flow cytometry analysis, cultured cells were trypsinized, washed, and resuspended in phosphate-buffered saline plus 0.5% bovine serum albumin. They were incubated with Allophycocyanin (APC)-conjugated antibodies on ice for 20 mins in the dark. Data were collected with a FACS Calibur flow cytometer (BD Biosciences) and analyzed using FlowJo software (Tree Star). The detailed information for all antibodies is listed in Supplementary Table S2. Luciferase reporter assay The pMiR-Report plasmids (with or without mature miR-192-5p binding sites) were transfected together with pRL-CMV vector containing Renilla luciferase. The Firefly and Renilla luciferase activities were measured 24 hours after transfection using Dual-Luciferase Reporter Assay (Promega, CA) with a PerkinElmer luminometer. Each experiment was performed in triplicate and repeated at least three times. To examine the promoter activity of miR-192-5p, HCC cells (24-well plate) were transfected with 300 ng pT3-EF1α vector or pT3-EF1α-cMyc vector, or 50 nM siRNAs (negative control or MYC siRNAs) on the first day and then transfected with 300ng pGL-192 promoter constructs and 10 ng pRL-CMV on the next day. The Firefly and Renilla luciferase activities were then measured 24 hours after transfection. Statistics Spearman’s rank correlation was performed to identify metabolites associated to miR-192-5p. Student’s t-test and Mann-Whitney rank test were used for statistical analysis of comparative data between groups. Two-way ANOVA was used to compare the glucose consumption, lactate production and cell viability of HCC cells at different time points. Hierarchical clustering analysis was performed by GENESIS software version 1.7.6 developed by Alexander Sturn (IBMT-TUG, Graz, Austria). Pearson correlation was used to identify genes correlated with miR-192-5p. Kaplan–Meier survival analysis was used to compare patient survival based on prediction results using GraphPad Prism V7.0 (San Diego, CA), and the p-value was generated by the Cox–Mantel log-rank test. Gene set enrichment analysis (GSEA) in the Molecular Signatures Database was performed using GSEA V3.0. All p-values were 2-sided. A p-value of ≤0.05 was regarded as statistically significant. Results Glycolysis-related metabolites and genes were highly expressed in HCC cases with low miR-192-5p expression We first investigated the metabolic features in HCCs with low miR-192-5p levels using HCC Cohort 1 with available metabolome and transcriptome data (Fig. 1a). We performed an integration analysis of miR-192-5p with the global metabolome in tumor tissues among 22 HCC patients, and found that 17 metabolites were significantly correlated with miR-192-5p (|r-value| >0.4, Fig. 1b, Supplementary Table S3). Among them, 7 metabolites presented |r-value| >0.5 and three of them were glycolysis-related metabolites, i.e., Glucose-6-Phosphate (G6P), Fructose-6-Phosphate (F6P), and Nicotinamide Adenine Dinucleotide Phosphate (NADPH). Meanwhile, 652 genes were significantly correlated with miR-192-5p with |r-value| >0.4, revealed by an integration analysis of miR-192-5p with mRNA transcriptome in tumor tissues among 176 HCC patients (Fig. 1c). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using these genes revealed 13 enriched metabolic features (p<0.001), three of which were associated with glycolysis and glycolysis-related pathways. These results suggest an altered glycolytic feature in HCC cases with low miR-192-5p expression. Available glycolysis-related metabolites and genes in our profiling data (Fig. 1d) were then compared between HCCs with high miR-192-5p levels (termed HCC 192High ) and HCCs with low miR-192-5p levels (termed HCC 192Low ), based on a miR-192-5p median cut-off in HCC tumors. Levels of G6P, F6P, and NADPH were significantly higher in tumors from HCC 192Low patients than HCC 192High patients, while no difference was found in their non-tumor tissues (Fig. 1e, Supplementary Fig. S1a). Consistently, many genes coding for key glycolytic enzymes such as GLUT1, HK2, PFKFB3, PFKP, and PKM2 were significantly upregulated in HCC 192Low tumors compared to HCC 192High tumors (Fig. 1e) but showed negligible alteration in their non-tumor tissues (Supplementary Fig. S1b). MCT1 was used as a negative control due to its main role in lactate import, but not in glycolysis[ 34 , 35 ]. These results demonstrate a hyperglycolytic metabolic feature in HCC cases with low miR-192-5p level. We further investigated the hyperglycolytic feature in CSC + HCC cases, i.e., cases with the top quartile expression of CSC markers as previously defined [ 23 ]. Hierarchical clustering analysis with glycolytic genes in Cohort 1 revealed two HCC subgroups with different expression levels of glycolytic genes. Consistently, in HCC subgroup with high expression levels of glycolysis-related genes, miR-192-5p level was low while various groups of CSC + HCC cases were enriched (Fig. 1f). Statistical analysis also showed that glycolytic genes were expressed at significantly higher levels in various groups of CSC + HCCs than in CSC - HCCs, but no difference was observed in the comparisons of their non-tumor tissues (Supplementary Fig. S1c-d). Comparable data were observed in Cohort 2 with 372 HCC patients (Supplementary Fig. S1e). Together, these data indicated that the hyperglycolytic feature was present in various groups of CSC + HCCs with low level of miR-192-5p. HCC cells with miR-192-5p loss were hyperglycolytic. We next investigated the role of miR-192-5p in regulating glycolysis. In Huh7, glycolytic genes with significant differential expression between HCC 192High and HCC 192Low patients were examined and eight out of nine genes showed a significant up-regulation after suppressing miR-192-5p (Supplementary Fig. S2a). Meanwhile, Huh7 cells with suppressed miR-192-5p exhibited a distinctly increased ECAR and a reduced OCR (Fig. 1g). The extracellular acid produced by cells is derived from lactate produced by glycolysis and CO 2 produced during respiration. OCR mainly represents mitochondrial respiration. Therefore, the increased ECAR in Huh7 cells with suppressed miR-192-5p was mainly due to lactate produced from glycolysis but not CO 2 from mitochondrial respiration. Consistently, over-expressed miR-192-5p in HLF and HLE cells lowered the extracellular acid production from glycolysis as shown by a reduced ECAR but an increased OCR (Fig. 1h). Most examined glycolytic genes were significantly reduced by miR-192-5p overexpression in both HLF and HLE cells (Supplementary Fig. S2b). These results demonstrate an important role of miR-192-5p in modulating a Warburg-like effect in HCC cells. To better elucidate the role of miR-192-5p in regulating glycolysis, we established miR-192-5p knockout (termed 192KO) clones from two human HCC cell lines HLF and HLE. A 69bp DNA fragment was deleted in the 192KO mixture clones as well as selected single 192KO clones (Supplementary Fig. S2c-d). One 192KO single clone from each HCC cell lines was used and miR-192-5p expression was undetectable in HLE-192KO and HLF-192KO cells (Fig. 2a). As a control, the expression of miR-194, a nearby miRNA of miR-192-5p, was not affected. As expected, HLF-192KO cells displayed significantly increased CSC features, such as increased populations of CD44 + , CD24 + and EpCAM + CSCs (Fig. 2b, Supplementary Fig. S2e); increased mRNA levels of multiple CSC biomarkers and reduced expression of a differentiation-related gene CYP1A2 (Supplementary Fig. S2f); and enlarged and more spheroid formation (Supplementary Fig. S2g). HLE-192KO cells displayed increased CSC features at a moderate level (Fig. 2b, Supplementary Fig. S2e-f). Consistently, these two 192KO lines also showed the hyperglycolytic features. As shown in Fig 2c, five glycolytic enzymes, i.e., GLUT1, HK2, PFKFB3, ALDOA, and PKM2, as well as c-Myc presented higher levels in 192KO cells than in wild-type cells. 192KO lines also exhibited increased ECARs but decreased OCRs (Fig. 2d), indicating that miR-192-5p loss largely increased the glycolysis-related extracellular acidification. Consistent data were also noticed in other 192KO clones of both HLF and HLE HCC cell lines (Supplementary Fig. S3a-b). Furthermore, overexpressed miR-192-5p in HLF-192KO cells significantly reduced CSC features and lactate accumulation in the culture medium (Fig. 2e). As a control, the intracellular lactate remained unchanged. Meanwhile, over culturing time, lactate gradually accumulated in the medium and was significantly higher in both HLF-192KO and HLE-192KO cells compared to their corresponding wild-type cells, which could be lowered by overexpressed miR-192-5p (Fig. 2f). A lower pH value was also observed in 192KO cells indicated by the orange/yellow medium vs. the pink medium of wild-type cells at 72 hours after seeding. Comparable data on lactate production were seen in other HLF and HLE 192KO clones (Supplementary Fig. S3c). Consistent data were also obtained through the detection of lactate using non-targeted metabolomics in HLF cells with different expression of miR-192-5p in both internal cells and culture medium (Supplementary Fig S3d). Together, miR-192-5p loss in HCC cells led to a hyperglycolytic phenotype. HCC cells with miR-192-5p loss had high glucose consumption We further examined glucose consumption among HCC cells with different levels of miR-192-5p as well as between HCC cells and their co-cultured non-tumor cells. As shown in Fig. 3a, both HLF-192KO and HLE-192KO cells exhibited significantly higher glucose consumption than HCC cells overexpressing miR-192-5p. In Huh7 cells, suppressing miR-192-5p increased their glucose usage (Supplementary Fig. S4a). Consistently, 192KO cells were more sensitive after exposure to 2-DG, a glucose analog, as shown by the significantly reduced cell viability compared to cells with miR-192-5p expression (Fig. 3b). Comparable data were obtained in HuH7 cells (Supplementary Fig. S4b). In co-culture systems of HCC cells with LX2, HL7702, and THP1, we further compared their glucose uptake via 2-NBDG uptake assay. HLF HCC cells infected with pmiR-ctrl/RFP or pmiR-192/RFP lentiviruses were used, and red fluorescent labeling efficiency was nearly 100% (Supplementary Fig. S4c). In this system, with or without co-culturing with other cells, HLF-192KO cells consistently showed higher 2-NBDG uptake than wild-type cells (Fig. 3c-d). In contrast, LX2 and HL7702 in co-culture with HLF-192KO cells exhibited lower 2-NBDG uptake compared to those in co-culture with HLF-WT cells. Moreover, forced-expression of miR-192-5p in HLF-192KO cells reduced the 2-NBDG uptake in HLF cells but increased 2-NBDG uptake in LX2 and HL7702 cells in the co-culture system. The alteration of 2-NBDG uptake was not observed in THP1 from our co-culture system (Fig. 3d, Supplementary Fig. S4d). Similar data were seen in HLE cells as well as in HLE cells co-cultured with LX2 and HL7702 (Supplementary Fig. S4e-f). These results demonstrate that HCC cells with loss of miR-192-5p actively utilize glucose from their environment to ensure a hyperglycolysis status. GLUT1, PFKFB3 and c-Myc were miR-192-5p bona fide targets and contributed to glycolytic and stemness features of HCC cells To investigate the target genes of miR-192-5p in regulating glycolysis flow, we assessed genes negatively correlated with miR-192-5p in 176 HCC cases (r0.2, p<0.01). Among these 554 genes, two main groups were observed. One group contained genes related to cell migration as we previously reported [ 23 ]. The other group included eight glycolysis-related genes (Fig. 4a) and three of them ( GLUT1, HK2 , and PKM2 ) were reported targets of c-Myc, an important regulator of glycolysis[ 36 , 37 ]. Next, using TargetScan and manual miRNA target prediction, we found that four of these glycolytic genes ( PFKFB3, GLUT1, MCT4 , and MYC ) contained miR-192-5p binding sites in their 3’UTR and/or coding regions (Fig. 4a). In HLF and HLE cells, miR-192-5p overexpression reduced the protein levels of PFKFB3, GLUT1, and c-Myc, but not that of MCT4 (Fig. 4b). Further, the predicted miR-192-5p binding regions of these three genes were cloned into a luciferase reporter and forced expression of miR-192-5p reduced the luciferase activities when the wild-type sequences for PFKFB3 and GLUT1 as well as the #2 binding site of MYC were present (Fig. 4c). These effects were significantly reduced when the corresponding miR-192-5p binding sites were mutated. Moreover, silencing PFKFB3, GLUT1, or c-Myc with 2 siRNAs for each gene reduced ECAR in both HLF-WT and HLF-192KO cells (Fig. 4d, Supplementary Fig. S5a). In HLF-192KO cells, silencing PFKFB3, GLUT1 or c-Myc notably reduced the ECAR rate to a level similar to that of HLF-WT cells with silencing of these genes. Comparable data were also noticed in HCC patients from Cohorts 1 and 2. PFKFB3 , GLUT1 , and MYC presented higher levels in HCC 192Low tumors compared to HCC 192High tumors (Supplementary Fig. S5b-c). These results indicate that PFKFB3, GLUT1, and c-Myc are miR-192-5p targets and are involved in the hyperglycolysis caused by miR-192-5p loss. PFKFB3, GLUT1, and c-Myc were reported to maintain stemness features in cancer at certain levels [ 38-41 ]. Consistently, si- PFKFB3, si- GLUT1 or double knockdown led to reduced levels of four CSC biomarkers, i.e., CD44 , CD24 , EpCAM and CD90 , as determined by RT-qPCR (Fig. 4e). Flow cytometry analysis also showed that si- PFKFB3 and GLUT1 reduced the populations of CD44 + and CD24 + CSCs (Fig. 4f). Meanwhile, si- MYC seemed to only reduce CD44 + CSCs moderately, but not the CD24 + CSC population (Supplementary Fig. S5d). Together, these data demonstrate that three glycolytic regulators, PFKFB3, GLUT1 and c-Myc were bona fide targets of miR-192-5p, and they contributed to both hyper-glycolysis and CSC features of HCCs caused by loss of miR-192-5p. c-Myc suppressed miR-192-5p transcription, ensuring a positive feedback of high c-Myc/low miR-192-5p in hyperglycolytic CSC + HCCs In our previous miRNA profiles of tumors and non-tumors from a hydrodynamic injection HCC FVB mouse model[ 31 ], miR-192-5p expression was significantly reduced in c-Myc-induced HCCs (Fig. 5a). Further, in a hydrodynamic injection HCC ICR mouse model, miR-192-5p level was also reduced >100 times in c-Myc induced HCCs but was not much in Ras-induced HCCs when compared to corresponding non-HCC liver tissues. In four different HCC cell lines, si- MYC led to an increased expression of miR-192-5p, while forced expression of c-Myc reduced the level of miR-192-5p (Fig. 5b). These data indicate that c-Myc might regulate miR-192-5p transcription. Consistently, among four different lengths of miR-192-5p promoter regions, the -266 nt to +186 nt region showed the strongest promoter activity (Fig. 5c) and the miR-192-5p promoter activity (-266 nt to +186 nt) was reduced by exogenous c-Myc, while enhanced by si- MYC (Fig. 5d). Wild-type p53 could bind to the miR-192-5p promoter region and induce its expression[ 23 , 42 ]. Consistently, in HepG2 cells with wild-type TP53, the expression of miR-192-5p was induced by p53 via exposure to Nutlin-3a (an MDM2 antagonist to stabilize p53) and reduced by silencing of p53 (Fig. 5e, Supplementary Fig. S6a-b). Over-expressed c-Myc significantly suppressed miR-192-5p expression in HCC cells with either activated p53 or silenced p53, indicating that c-Myc-mediated miR-192-5p down-regulation was independent on p53. Comparable data were noticed in HCC patients. In both HCC cohorts, hierarchical clustering analysis revealed two subgroups with distinct c-Myc activation status based on 76 c-Myc target genes from the online Human MYC Targets Profiler (Supplementary Fig. S6c-d). In Cohort 1, miR-192-5p expression in the c-Myc activation subgroup was significantly lower than that in c-Myc non-activation subgroup (Fig. 5f). In Cohort 2, miR-192-5p expression was always significantly lower in each c-Myc activation subgroup than in the corresponding non-activation subgroup, which was independent of TP53 mutation and mir-192 promoter methylation (Fig. 5f-g). In addition, different groups of CSC + HCCs consistently presented a low level of miR-192-5p, a high level of c-Myc activation and high frequency of MYC amplification (Supplementary Fig. S7). Together, c-Myc suppressed miR-192-5p transcription, which led to a positive feedback of high c-Myc/low miR-192-5p in CSC + HCC cells with glycolytic feature. Overproduced lactate from CSC + HCCs activated the ERK pathway in environmental non-tumor cells, and this effect further increased HCC cell stemness and malignancy features As the end product of glycolysis, the continuously produced lactate in hyperglycolytic miR-192-5p-loss HCC cells might affect their environment and contribute to HCC malignancy. The transport of lactate across the plasma membrane is mainly catalyzed by MCT1 and MCT4, with MCT1 typically involved in the import while MCT4 in export of lactate[ 34 , 35 ]. In HCC patients, the expression ratio of MCT1 vs. MCT4 showed no difference between tumor and non-tumor tissues of Cohort 1 but was significantly higher in non-tumor tissues than tumor tissues of Cohort 2 (Supplementary Fig. S8a), indicating the possibility of lactate uptake by environmental non-tumor cells. Lactate could also stabilize NDRG3, which in turn activated the ERK pathway to promote cell malignancy[ 43 ]. Consistently, lactate treatment stimulated ERK phosphorylation noticeably in non-tumor cells of HCC microenvironment, i.e., LX2, THP1 and HL7702 cells (Fig. 6a). In a chamber co-culture system, pERK level was increased in LX2 and THP1 cells when co-cultured with HLF-192KO cells compared to when co-cultured with HLF-WT (Fig. 6b). Moreover, pERK was further reduced in LX2 and THP1 when exposed to HLF-192KO cells with miR-192-5p overexpression (Fig. 6b). Similar results were observed from LX2 and THP1 cells co-cultured with HLE cells (Supplementary Fig. S8b). However, pERK was unaltered in HL7702 cells in this co-culture system. Thus, HCC cells with miR-192-5p loss could actively affect certain non-tumor cells via increased production of lactate. We have also found that lactate-induced pERK in environmental non-tumor cells partially relied on NDRG3 and MCT1 (Fig. 6c-e). LX2 and THP1 cells expressed relatively high levels of MCT1 and NDRG3 (Fig. 6c). Silencing NDRG3 or MCT1 in LX2 and THP1 cells reduced the level of lactate-induced pERK (Fig. 6d-e). Low expression levels of MCT1 and NDRG3 in HL7702 were consistent with its minor response to lactate (data not shown). We then explored the effects of an altered lactate/MCT1/NDRG3/pERK axis in LX2 or THP1 cells on the malignancy features in HCC cells (Fig. 7a). HLF-192KO/RFP cells were co-cultured with LX2 cells pre-transfected with si-Ctrl, or si- NDRG3 , or si- MCT1 (termed LX2 si-Ctrl , LX2 si- NDRG3 , and LX2 si- MCT1 , respectively). Wound-healing assay of red-fluorescence HCC cells showed that cell migration of HLF-192KO cells was slower upon co-culture with LX2 si- NDRG3 or LX2 si- MCT1 than with LX2 si-Ctrl (Fig. 7b). Consistent data were observed in HLF-192KO cells co-cultured with THP1 (Fig. 7b). Moreover, spheroid assays of HLF-192KO cells were performed in different conditioned medium settings. The number of spheroids of HLF-192KO cells was significantly lower under exposure to conditioned medium from co-culture of HLF-192KO with LX2 si- NDRG3 or LX2 si- MCT1 than from co-culture of HLF-192KO with LX2 si-Ctrl (Fig. 7c). CD44 + and CD24 + HLF-192KO populations were also significantly reduced when they were co-cultured with LX2 si- NDRG3 or LX2 si- MCT1 (Fig. 7d, Supplementary Fig. S8c). As a control, HLF-192KO with overexpressed miR-192-5p did not exhibit significant alteration of migration, spheroid formation and CSC populations when co-cultured with different LX2 cells (Fig. 7b-d). Consistent data were obtained using a second set of siRNAs for NDRG3 and MCT1 (Supplementary Fig. S9). Therefore, blocking the lactate/ERK pathway in HCC microenvironmental cells suppressed the malignancy and stemness features of HCC cells in co-culture experiments. In both HCC cohorts, patients were divided into four groups based on miR-192-5p expression in their tumor tissues (HCC 192Low and HCC 192High , medium cut-off) and levels of NDRG3 and MCT1 in non-tumor tissues (NT High_NDRG3 or MCT1 and NT Low_NDRG3 and MCT1 , medium cut-offs). There was no expression difference of NDRG3 or MCT1 in non-tumor tissues between HCC 192Low and HCC 192High patients (Fig. 7e, Supplementary Fig. S10a). In Cohort 1, patients with HCC 192High NT Low_NDRG3 and MCT1 had the best prognosis, as shown by a prolonged time to recurrence and overall survival. In HCC 192Low subgroup, patients with NT High_NDRG3 or MCT1 had worse prognosis compared to patients with NT Low_NDRG3 or MCT1 (Fig. 7f). Similar but less significant data were obtained in Cohort 2, which might be due to the limited number of patients (n=49) with available non-tumor mRNA data (Supplementary Fig. S10b). GSEA analysis in HCC 192Low patients revealed that several stem cell related gene-sets were enriched in patients with NT High_NDRG3 or MCT1 (Supplementary Fig. S10c-d). Together, HCC cells with miR-192-5p loss exhibited a highly malignant feature when they were surrounded by environmental non-tumors with high MCT1 or NDRG3 expression. Discussion HCC is one of the most malignant cancers worldwide. Although a trend of decreasing incidence has been observed for many cancers, the incidence of HCC is still rising in both developing and developed countries[ 44 ]. CSCs are thought to be responsible for tumorigenesis as well as tumor metastasis and targeting CSCs holds the hope of eliminating cancer[ 3 , 4 , 6 ]. In the last decade, researchers’ attention has been drawn to the field of hepatic CSCs, leading to the identification of multiple hepatic CSC biomarkers. Accordingly, accumulated studies have revealed gene signatures and regulatory signaling pathways of these identified hepatic CSCs [ 18 ]. In this study, we aimed to understand essential metabolic features of hepatic CSCs. MiR-192-5p, a liver-abundant miRNA, is functionally important in suppressing the stemness and malignancy features of HCC cells, and its genetic silencing frequently occurs in multiple groups of CSC + HCCs, i.e., EpCAM + , CD90 + , CD133 + , CD44 + , and CD24 + HCCs [ 23 ]. Here, we revealed that the miR-192-5p loss drove HCC cells to a hyperglycolytic metabolism status via targeting c-Myc and two glycolytic enzymes GLUT1 and PFKFB3, and a c-Myc/miR-192-5p positive feedback pathway. Five groups of CSC + HCCs presented a hyperglycolytic signature, high level of c-Myc and low miR-192-5p. Mechanistically, loss of miR-192-5p led to increased expression of its target genes, i.e., PFKFB3, GLUT1, and c-Myc, which facilitated a Warburg effect favored by cellular growth and stemness. Moreover, c-Myc also directly suppressed miR-192-5p expression, ensuring a positive feedback to drive hyperglycolytic features in CSC + HCCs. In addition, targeting PFKFB3, GLUT1, and c-Myc could also reduce hepatic CSC populations. Thus, hyperglycolysis was a key feature for CSC + HCCs and effective methods of targeting hyperglycolytic cells in combination with conventional anticancer methods might effectively eliminate CSCs, leading to a significant reduction of HCC recurrence and metastasis. However, it remains unclear whether liver cells gain their CSC features and hyperglycolysis features at the same time or if one is the consequence of the other. The efforts to answer these questions will further improve our understanding of HCC initiation and may open new avenues to prevent HCC tumor initiation or metastasis. The HCC microenvironment is complex, including but not limited to activated hepatic stellate cells, Kupffer cells, and hypoxia etc, and is believed to actively interact with HCC tumors as well as affect their malignancy status [ 7 , 45 , 46 ]. In our study, LX2 and THP1 (representing hepatic stellate cells and macrophages respectively) were used to co-culture with HCC cells in vitro. We found that hyperglycolytic HCC cells with miR-192-5p loss nourished LX2 and THP1 cells with overproduced lactate, which stimulated the pERK pathway in LX2 and THP1 cells via NDRG3/MCT1. Consistently, in the in vitro co-culture system, LX2 and THP1 cells with high levels of NDRG3 and/or MCT1 facilitated a more malignant phenotype of HCC cells with miR-192-5p loss. In vivo , HCC patients with low level of miR-192-5p in their tumor and high level of NDRG3/MCT1 in their non-tumor had the shorter overall survival when compared to other patient subgroups. Thus, hepatic CSC-related metabolism affected the CSC microenvironment via accumulated lactate, and lactate in turn interacted with environmental cellular components (LX2 and THP1) that contributed to HCC stemness features. This was consistent with reported studies demonstrating that the structural and cellular components of the CSC microenvironment played important roles in CSC plasticity, proliferation and invasion [ 8 , 9 ]. It remains undiscovered in our study how the pERK pathway in THP1 or LX2 cells promote HCC malignancy. It would be interesting to perform an in-depth study on cytokines or exosomes being produced by THP1 and LX2 cells through the lactate/ERK axis, and to explore their roles in promoting the malignancy of tumor cells. Both lactate and LX2/THP1 cells in the hepatic CSC environment contributed to the stemness and malignant features of HCC cells with miR-192-5p loss. Thus, we would need more insight on the interaction of tumor initiating cells and their environment to offer a better solution of cancer therapy regarding CSC elimination. Investigating the CSC-niche metabolism would be also important for the further development of metabolic targeting approaches, which might significantly impact cancer therapy. As a class of non-coding RNAs, miRNAs participate in many key events in hepatocarcinogenesis. As a liver-specific and abundant miRNA, miR-192-5p is significantly down-regulated in HCC tumor tissues due to a high frequency of TP53 mutation and its promoter hypermethylation [ 23 ]. In HCC cells, miR-192-5p loss not only led to increased populations of hepatic CSCs [ 23 ] and cell malignancy phenotype [ 23-25 ], but also a hyperglycolytic phenotype as revealed in this study (Fig 7g). These findings suggested a potentially important role of miR-192-5p in regulating hepatic carcinogenesis. It will be essential to establish conditional miR-192-5p knockout mice and further investigate whether and at which stages miR-192-5p regulates hepatic carcinogenesis, the metabolism of HCC tumor and tumor microenvironment, as well as roles of the miR-192-5p-related metabolic pathways in regulating hepatic carcinogenesis in the mouse model. Conclusions Taken together, loss of miR-192-5p promoted hyperglycolysis and stemness features in CSC + HCC cells through upregulating c-Myc and two glycolytic enzymes ( PFKFB3 and GLUT1), and through a miR-192-5p/c-Myc positive feedback circuit. Moreover, HCC cells with miR-192-5p loss could also actively interact with their environmental non-tumor cells to produce a more aggressive and stemness-related feature via the over-produced lactic acid by HCC cells and consequently activating a lactate/MCT1/NDRG3/pERK pathway in their environmental non-tumor cells (Fig. 7g). Abbreviations Cancer stem cell, CSC; Dulbecco’s modified Eagle’s medium, DMEM; Extracellular acidification rate, ECAR; Fructose-6-Phosphate, F6P; Gas chromatography/time-of-flight mass spectrometry, GC/TOF-MS; Gene Set Enrichment Analysis, GSEA; Glucose-6-Phosphate, G6P; Hepatocellular carcinoma, HCC; Hepatic stellate cell, HSC; Japanese Collection of Research Biosources Cell Bank, JCRB; Kyoto Encyclopedia of Genes and Genomes, KEGG; Oxygen consumption rate, OCR; Quantitative reverse transcription polymerase chain reaction, qRT-PCR; Single-guide RNAs, sgRNAs; Nicotinamide Adenine Dinucleotide Phosphate, NADPH; The Cancer Genome Atlas, TCGA; 2-deoxyglucose, 2-DG; 2-(N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino)-2-deoxy-D-glucose, 2-NBDG. Declarations Ethics approval and consent to participate Not applicable Consent for publication Not applicable Availability of data and materials miRNA transcriptome in HCC Cohort 1, GSE6857; mRNA transcriptome in HCC Cohort 1, GSE14520; metabolomics data in HCC Cohort 1 were referenced in Ref 12. miRNA sequencing and mRNA sequencing data were downloaded from TCGA. TP53 mutation status and MYC amplification information in HCC Cohort 2 were obtained from www.cbioportal.org/index.do. Competing interests The authors have no conflicts of interest with the material presented in this manuscript. Funding This work was supported by National Natural Science Foundation of China (No.81672905, 81874054, and 82073055) (J. Ji), National Key R&D Program of China (2018YFA0800504) (J. Ji), the Fundamental Research Funds for the Central Universities in China (J. Ji), and the Thousand Young Talents Plan of China (J. Ji). X.W.W. and A.B. were supported by the intramural program of the Center for Cancer Research, the U.S. National Cancer Institute. S.R. was supported by the German Research Foundation (DFG, Deutsche Forschungsgemeinschaft) within SFB/TR 209 Liver Cancer Project-ID 314905040. Authors’ contributions Conception and design: Y. Gu., and J. Ji.; Development of methodology: Y. Gu., F. Ji., N. Liu., J. Ji.; Acquisition of data: Y. Gu., N. Liu., Y. Zhao., S. Hu., J. Ji.; Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): Y. Gu., F. Ji., X. Wei., J. Ji.; Writing, review, and/or revision of the manuscript: Y. Gu., X. Wang., A. Budhu., S. Roessler., J. Ji.; Administrative, technical, or material support: X. Zheng., B. Zhao., X. Wang., A. Budhu., W. Jia., J. Ji.; Study supervision: J. Ji. Acknowledgments We thank Dr. Chaohui Yu for the HL7702 cell line, Dr. Long Zhang for the THP1 cell line, Dr. Xing Guo for PX459 plasmid and Dr. Hai Song for the anti-pERK and anti-ERK antibodies. The authors thank our core facility at Life Sciences Institute for their support on flow cytometry analysis and measurement using Sea Horse XF Analyzers. References Gorrini C, Mak TW: Glutathione Metabolism: An Achilles' Heel of ARID1A-Deficient Tumors. Cancer Cell. 2019;35:161-3. Ward PS, Thompson CB: Metabolic reprogramming: a cancer hallmark even warburg did not anticipate. Cancer Cell. 2012;21:297-308. Gupta PB, Chaffer CL, Weinberg RA: Cancer stem cells: mirage or reality? Nat Med. 2009;15:1010-2. Oikawa T: Cancer Stem cells and their cellular origins in primary liver and biliary tract cancers. Hepatology. 2016;64:645-51. Plaks V, Kong N, Werb Z: The cancer stem cell niche: how essential is the niche in regulating stemness of tumor cells? Cell Stem Cell. 2015;16:225-38. Deshmukh A, Deshpande K, Arfuso F, Newsholme P, Dharmarajan A: Cancer stem cell metabolism: a potential target for cancer therapy. Mol Cancer. 2016;15:69. Maman S, Witz IP: A history of exploring cancer in context. Nat Rev Cancer. 2018;18:359-76. Batlle E, Clevers H: Cancer stem cells revisited. Nat Med. 2017;23:1124-34. Vermeulen L, De Sousa EMF, van der Heijden M, Cameron K, de Jong JH, Borovski T et al : Wnt activity defines colon cancer stem cells and is regulated by the microenvironment. Nat Cell Biol. 2010;12:468-76. Ahmed N, Escalona R, Leung D, Chan E, Kannourakis G: Tumour microenvironment and metabolic plasticity in cancer and cancer stem cells: Perspectives on metabolic and immune regulatory signatures in chemoresistant ovarian cancer stem cells. Semin Cancer Biol. 2018;53:265-81. Peiris-Pages M, Martinez-Outschoorn UE, Pestell RG, Sotgia F, Lisanti MP: Cancer stem cell metabolism. Breast Cancer Res. 2016;18:55. Villanueva A: Hepatocellular Carcinoma. N Engl J Med. 2019;380:1450-62. Lee TK, Castilho A, Cheung VC, Tang KH, Ma S, Ng IO: CD24(+) liver tumor-initiating cells drive self-renewal and tumor initiation through STAT3-mediated NANOG regulation. Cell Stem Cell. 2011;9:50-63. Yamashita T, Ji J, Budhu A, Forgues M, Yang W, Wang HY et al : EpCAM-positive hepatocellular carcinoma cells are tumor-initiating cells with stem/progenitor cell features. Gastroenterology. 2009;136:1012-24. Yang ZF, Ho DW, Ng MN, Lau CK, Yu WC, Ngai P et al : Significance of CD90+ cancer stem cells in human liver cancer. Cancer Cell. 2008;13:153-66. Ma S, Tang KH, Chan YP, Lee TK, Kwan PS, Castilho A et al : miR-130b Promotes CD133(+) liver tumor-initiating cell growth and self-renewal via tumor protein 53-induced nuclear protein 1. Cell Stem Cell. 2010;7:694-707. Ma S, Chan KW, Hu L, Lee TK, Wo JY, Ng IO et al : Identification and characterization of tumorigenic liver cancer stem/progenitor cells. Gastroenterology. 2007;132:2542-56. Gu Y, Zheng X, Ji J: Liver cancer stem cells as a hierarchical society: yes or no? Acta Biochim Biophys Sin (Shanghai). 2020;52:723-35. Prasetyanti PR, Medema JP: Intra-tumor heterogeneity from a cancer stem cell perspective. Mol Cancer. 2017;16:41. Song K, Kwon H, Han C, Zhang J, Dash S, Lim K et al : Active glycolytic metabolism in CD133(+) hepatocellular cancer stem cells: regulation by MIR-122. Oncotarget. 2015;6:40822-35. Chen CL, Uthaya Kumar DB, Punj V, Xu J, Sher L, Tahara SM et al : NANOG Metabolically Reprograms Tumor-Initiating Stem-like Cells through Tumorigenic Changes in Oxidative Phosphorylation and Fatty Acid Metabolism. Cell Metab. 2016;23:206-19. Budhu A, Roessler S, Zhao X, Yu Z, Forgues M, Ji J et al : Integrated metabolite and gene expression profiles identify lipid biomarkers associated with progression of hepatocellular carcinoma and patient outcomes. Gastroenterology. 2013;144:1066-75 e1. Gu Y, Wei X, Sun Y, Gao H, Zheng X, Wong LL et al : miR-192-5p Silencing by Genetic Aberrations Is a Key Event in Hepatocellular Carcinomas with Cancer Stem Cell Features. Cancer Res. 2019;79:941-53. Lian J, Jing Y, Dong Q, Huan L, Chen D, Bao C et al : miR-192, a prognostic indicator, targets the SLC39A6/SNAIL pathway to reduce tumor metastasis in human hepatocellular carcinoma. Oncotarget. 2016;7:2672-83. Yang YM, Lee WH, Lee CG, An J, Kim ES, Kim SH et al : Galpha12 gep oncogene deregulation of p53-responsive microRNAs promotes epithelial-mesenchymal transition of hepatocellular carcinoma. Oncogene. 2015;34:2910-21. Budhu A, Jia HL, Forgues M, Liu CG, Goldstein D, Lam A et al : Identification of metastasis-related microRNAs in hepatocellular carcinoma. Hepatology. 2008;47:897-907. Ji J, Shi J, Budhu A, Yu Z, Forgues M, Roessler S et al : MicroRNA expression, survival, and response to interferon in liver cancer. N Engl J Med. 2009;361:1437-47. Ji J, Yamashita T, Budhu A, Forgues M, Jia HL, Li C et al : Identification of microRNA-181 by genome-wide screening as a critical player in EpCAM-positive hepatic cancer stem cells. Hepatology. 2009;50:472-80. Roessler S, Long EL, Budhu A, Chen Y, Zhao X, Ji J et al : Integrative genomic identification of genes on 8p associated with hepatocellular carcinoma progression and patient survival. Gastroenterology. 2012;142:957-66 e12. Sun Y, Ji F, Kumar MR, Zheng X, Xiao Y, Liu N et al : Transcriptome integration analysis in hepatocellular carcinoma reveals discordant intronic miRNA-host gene pairs in expression. Int J Biol Sci. 2017;13:1438-49. Tao J, Ji J, Li X, Ding N, Wu H, Liu Y et al : Distinct anti-oncogenic effect of various microRNAs in different mouse models of liver cancer. Oncotarget. 2015;6:6977-88. Ji J, Zhao L, Budhu A, Forgues M, Jia HL, Qin LX et al : Let-7g targets collagen type I alpha2 and inhibits cell migration in hepatocellular carcinoma. J Hepatol. 2010;52:690-7. Ji J, Zheng X, Forgues M, Yamashita T, Wauthier EL, Reid LM et al : Identification of microRNAs specific for epithelial cell adhesion molecule-positive tumor cells in hepatocellular carcinoma. Hepatology. 2015;62:829-40. Dimmer KS, Friedrich B, Lang F, Deitmer JW, Broer S: The low-affinity monocarboxylate transporter MCT4 is adapted to the export of lactate in highly glycolytic cells. Biochem J. 2000;350 Pt 1:219-27. Doherty JR, Cleveland JL: Targeting lactate metabolism for cancer therapeutics. J Clin Invest. 2013;123:3685-92. Sabnis HS, Somasagara RR, Bunting KD: Targeting MYC Dependence by Metabolic Inhibitors in Cancer. Genes (Basel). 2017;8. Stine ZE, Walton ZE, Altman BJ, Hsieh AL, Dang CV: MYC, Metabolism, and Cancer. Cancer Discov. 2015;5:1024-39. Li HM, Yang JG, Liu ZJ, Wang WM, Yu ZL, Ren JG et al : Blockage of glycolysis by targeting PFKFB3 suppresses tumor growth and metastasis in head and neck squamous cell carcinoma. J Exp Clin Cancer Res. 2017;36:7. Shibuya K, Okada M, Suzuki S, Seino M, Seino S, Takeda H et al : Targeting the facilitative glucose transporter GLUT1 inhibits the self-renewal and tumor-initiating capacity of cancer stem cells. Oncotarget. 2015;6:651-61. Zhang HL, Wang MD, Zhou X, Qin CJ, Fu GB, Tang L et al : Blocking preferential glucose uptake sensitizes liver tumor-initiating cells to glucose restriction and sorafenib treatment. Cancer Lett. 2017;388:1-11. Akita H, Marquardt JU, Durkin ME, Kitade M, Seo D, Conner EA et al : MYC activates stem-like cell potential in hepatocarcinoma by a p53-dependent mechanism. Cancer Res. 2014;74:5903-13. Pichiorri F, Suh SS, Rocci A, De Luca L, Taccioli C, Santhanam R et al : Downregulation of p53-inducible microRNAs 192, 194, and 215 impairs the p53/MDM2 autoregulatory loop in multiple myeloma development. Cancer Cell. 2010;18:367-81. Lee DC, Sohn HA, Park ZY, Oh S, Kang YK, Lee KM et al : A lactate-induced response to hypoxia. Cell. 2015;161:595-609. Bray F, Ferlay J, Soerjomataram I, Siegel RL, 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:394-424. Brodt P: Role of the Microenvironment in Liver Metastasis: From Pre- to Prometastatic Niches. Clin Cancer Res. 2016;22:5971-82. Craig AJ, von Felden J, Garcia-Lezana T, Sarcognato S, Villanueva A: Tumour evolution in hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol. 2020;17:139-52. Supplementary Files manuscriptSup1118.docx Cite Share Download PDF Status: Published Journal Publication published 30 Nov, 2020 Read the published version in Journal of Experimental & Clinical Cancer Research → Version 3 posted Editorial decision: Accept 18 Nov, 2020 Editor assigned by journal 17 Nov, 2020 Submission checks completed at journal 17 Nov, 2020 Editor invited by journal 17 Nov, 2020 You are reading this latest preprint version Show more versions Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-56821","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research","associatedPublications":[],"authors":[{"id":5028142,"identity":"87043ee8-ab46-4e24-bbcc-b8804ff2da8a","order_by":0,"name":"Yuanzhuo Gu","email":"","orcid":"","institution":"Zhejiang University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Yuanzhuo","middleName":"","lastName":"Gu","suffix":""},{"id":5028143,"identity":"ff8b7dfe-6010-4c20-846f-3c3e76884635","order_by":1,"name":"Fubo Ji","email":"","orcid":"","institution":"Zhejiang University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Fubo","middleName":"","lastName":"Ji","suffix":""},{"id":5028144,"identity":"0d9f99bf-2fe8-4f1c-ac41-3b77ecb7e43d","order_by":2,"name":"Niya Liu","email":"","orcid":"","institution":"Zhejiang University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Niya","middleName":"","lastName":"Liu","suffix":""},{"id":5028145,"identity":"49dc50e9-427d-4815-93ca-68ca416a33c8","order_by":3,"name":"Yongzhi Zhao","email":"","orcid":"","institution":"Zhejiang University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Yongzhi","middleName":"","lastName":"Zhao","suffix":""},{"id":5028146,"identity":"9bdf0511-35f2-4d0e-b5fb-20774379c755","order_by":4,"name":"Xiyang Wei","email":"","orcid":"","institution":"Zhejiang University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Xiyang","middleName":"","lastName":"Wei","suffix":""},{"id":5028147,"identity":"4cc059e2-51e5-4f88-b530-5e7b16fe405f","order_by":5,"name":"Shiyuan Hu","email":"","orcid":"","institution":"Zhejiang University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Shiyuan","middleName":"","lastName":"Hu","suffix":""},{"id":5028148,"identity":"7326092a-4239-4986-a94d-88b3c0c936cb","order_by":6,"name":"Wei Jia","email":"","orcid":"","institution":"Hong Kong Baptist University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Wei","middleName":"","lastName":"Jia","suffix":""},{"id":5028149,"identity":"a4f2601b-0a5c-4057-9f31-c6a3b2a05a18","order_by":7,"name":"Xin Wei Wang","email":"","orcid":"","institution":"National Cancer Institute","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Xin","middleName":"Wei","lastName":"Wang","suffix":""},{"id":5028150,"identity":"dd9a5691-0bc9-4c8c-b599-8b6f7f13c34b","order_by":8,"name":"Anuradha Budhu","email":"","orcid":"","institution":"National Cancer Institute","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Anuradha","middleName":"","lastName":"Budhu","suffix":""},{"id":5028151,"identity":"b5312d12-b372-4202-bdef-a85163910126","order_by":9,"name":"Juling Ji","email":"","orcid":"","institution":"Nantong University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Juling","middleName":"","lastName":"Ji","suffix":""},{"id":5028152,"identity":"66e2e578-0aeb-4aff-88e3-f4c340574600","order_by":10,"name":"Bin Zhao","email":"","orcid":"","institution":"Zhejiang University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Bin","middleName":"","lastName":"Zhao","suffix":""},{"id":5028153,"identity":"45879a6a-bbbc-45fe-8033-12a23eddee62","order_by":11,"name":"Stephanie Roessler","email":"","orcid":"","institution":"University Hospital Heidelberg","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Stephanie","middleName":"","lastName":"Roessler","suffix":""},{"id":5028154,"identity":"33906b28-af57-45b9-be01-d3777dcc69b2","order_by":12,"name":"Xin Zheng","email":"","orcid":"","institution":"EZKIT L.L.C.","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Xin","middleName":"","lastName":"Zheng","suffix":""},{"id":5028155,"identity":"b62f0b08-3d3c-4bb1-9df0-f96bca12de6c","order_by":13,"name":"Junfang Ji","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA6UlEQVRIiWNgGAWjYBACxgYGhg9AmpmNvYHBACpAUAvjDLAWngNEagEpmgGmJBLg1uIHzDNyDzb83FHLzif5/EExD4ON7IYDzM8e4LciL7Gx98xxZjbphARjHoY04w0H2MwN8GvJMX/A23YMpOUAUMvhxA0HeNgkCGgxbPwL0iJ5sAGo5T9xWpp522qY2SSYGYBaDhChpeeNYbNs2wFgIKcxGM4xSDaeeZjNDK8Ww3agw9621SXLtx9/ZvCmwk6273jzM/xaGsDU4WQgwWYAjkxmfOqBQB5C1dmB1D4goHgUjIJRMApGKAAAVahFygTCgGUAAAAASUVORK5CYII=","orcid":"https://orcid.org/0000-0002-3071-6494","institution":"Zhejiang University","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Junfang","middleName":"","lastName":"Ji","suffix":""}],"badges":[],"createdAt":"2020-08-10 12:25:46","currentVersionCode":3,"declarations":"","doi":"10.21203/rs.3.rs-56821/v3","doiUrl":"https://doi.org/10.21203/rs.3.rs-56821/v3","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s13046-020-01785-7","type":"published","date":"2020-11-30T15:02:01+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":3807128,"identity":"a56618d8-cee4-4165-929a-d561d2975d0b","added_by":"auto","created_at":"2020-11-24 21:26:42","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":212269,"visible":true,"origin":"","legend":"In HCCs, high expression of glycolysis-related metabolites and genes was associated with low miR-192-5p level.\na. Available omics dataset in Cohort 1 (n=176). b. Spearman correlation was performed between miR-192-5p and metabolome in 22 HCC cases of HCC Cohort 1. 17 metabolites were significantly correlated with miR-192-5p with |r-value| \u003e0.4 (p\u003c0.05). c. 652 genes significantly correlated with miR-192-5p (Pearson correlation, |r-value| \u003e0.4) in Cohort 1 were used for KEGG pathway enrichment analysis. d. Major metabolites and genes in glycolysis and glycolysis side pathways. e. Log2 ratios of glycolysis-related metabolites and genes in HCC192Low vs. HCC192High HCCs of Cohort 1. f. Hierarchical clustering analysis with glycolysis-related genes in Cohort 1. Red bar and green bar represent CSC biomarker positive and negative HCC cases, respectively. miR-192-5p expression levels were also shown. g. ECAR and OCR measurement in Huh7 cells infected with lentivirus miRZip-ctrl or miRZip-192. h. ECAR and OCR measurement in HLF and HLE cells infected with lentivirus pmiR-ctrl or pmiR-192.\n","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-56821/v3/aa0e98d7b7017abc1916e050.png"},{"id":3807129,"identity":"9ad20b83-aaf2-4a51-a403-c8f390db2e83","added_by":"auto","created_at":"2020-11-24 21:26:42","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":236027,"visible":true,"origin":"","legend":"MiR-192-5p knockout human HCC cell lines presented the increased CSC features and were highly glycolytic.\na. Schematic structure of miR-192-5p wild-type (WT) and knockout (KO) by CRISPR/Cas9 system. The knockout of miR-192-5p DNA fragment was assessed by Sanger sequencing. Levels of miR-192-5p and miR-194 was measured by RT-qPCR. b. Quantitative data of flow cytometry analysis using APC-conjugated antibodies against different CSC biomarkers in 192KO lines and WT lines from HLF and HLE cells. c. Western blot of glycolytic enzymes using 192KO lines and WT lines from HLF and HLE cells. d. ECAR and OCR were measured in 192KO lines and WT lines from HLF and HLE cells. e. RT-qPCR analysis of miR-192-5p, flow cytometry analysis of CD24+ and CD44+ populations, and lactate measurement using Lactate assay kit in HLF cells with different levels of miR-192-5p. f. Measurement of lactate in cultured medium from HLF and HLE cells at different time points. Representative images for color change of the corresponding cultured medium at 72hours after seeding are shown. \n","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-56821/v3/366d30f4b10a02cae5b1e951.png"},{"id":3807130,"identity":"8da78e9c-ad95-4b18-bfc5-8c4c469a1d03","added_by":"auto","created_at":"2020-11-24 21:26:43","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":115471,"visible":true,"origin":"","legend":"MiR-192-5p KO cells consumed more glucose from the environment.\na. Glucose concentration in culture medium from WT and 192KO lines of HLE and HLF cells infected with lentivirus pmiR-Ctrl or pmiR-192. b. Cell viability of WT and 192KO lines for HLE and HLF cells infected with lentivirus pmiR-Ctrl or pmiR-192 and exposed to different doses of 2-DG. c. 2-NBDG uptake by flow cytometry analysis of HLF-WT and -192KO cells infected with lentivirus pmiR-Ctrl (RFP) or pmiR-192 (RFP). d. 2-NBDG uptake for all cells in co-culture system. HLF cells were infected with lentivirus pmiR-Ctrl (RFP) or pmiR-192 (RFP) and co-cultured with LX2, HL7702 or THP1 cells. Two-way ANOVA analysis was performed for (a, b) and student t-test was performed for (c,d). *, p\u003c0.05. \n","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-56821/v3/127b7449331d947780383bc9.png"},{"id":3807131,"identity":"ca4be272-997c-4cbb-a908-a5628e39a59c","added_by":"auto","created_at":"2020-11-24 21:26:43","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":226511,"visible":true,"origin":"","legend":"c-Myc as well as two glycolytic enzymes were miR-192-5p targets. \na. Genes negatively correlated with miR-192-5p (r\u003c-0.3) and significantly altered in HCC192Low vs. HCC192High (Left panel) and predicted miR-192-5p binding sites in 3’-UTR or coding regions (CDS) of human PFKFB3, GLUT1, MYC and MCT4 (Right panel). b. Western blotting of PFKFB3, GLUT1, c-Myc and MCT4 in HLE and HLF cells infected with lentivirus pmiR-Ctrl or pmiR-192. c. Luciferase activities were measured using reporters with WT or MT miR-192-5p binding sites of three genes in the 3’-UTR or CDS of luciferase vector. 192pos reporter was used as the positive control. d. ECAR was measured in HLF-WT and HLF-192KO cells transfected with 50nM si-Ctrl, si-PFKFB3, si-GLUT1, and si-MYC respectively. e. RT-PCR analysis of five CSC surface markers in HLF-WT and HLF-192KO cells infected with si-Ctrl, si-PFKFB3, si-GLUT1, as well as the combination of si-PFKFB3 and si-GLUT1. UGT2B7 was measured as a negative control. f. Flow cytometry analysis of CD24+ and CD44+ populations in HLF-WT cells transfected with si-Ctrl, HLF-192KO cells transfected with si-Ctrl, and HLF-192KO cells transfected with the combined si-PFKFB3 and si-GLUT1. Student t-test was performed for (c.e.). *, p \u003c0.05.\n","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-56821/v3/ac36bbe59da38de4ca022a68.png"},{"id":3807132,"identity":"4b369c64-862c-4d8d-8f39-4fdd49c2d471","added_by":"auto","created_at":"2020-11-24 21:26:43","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":182801,"visible":true,"origin":"","legend":"C-Myc activation reduced miR-192-5p expression in HCC.\na. RT-qPCR analysis of miR-192-5p in murine oncogene-induced HCC samples and the corresponding liver controls. RNAs from 3-5 mice for each group were used. b. RT-qPCR analysis of MYC and miR-192-5p in HCC cell lines (Huh7, HLE, HLF and HepG2) transfected with Ctrl or MYC overexpressing plasmids, as well as with Ctrl or MYC siRNAs, respectively. c. Dual-luciferase assays for different miR-192-5p promoter fragments in HLE, HLF and HuH7 cells. d. Dual-luciferase assay using miR-192-5p promoter (-266nt, +186nt) in HLE, HLF and HuH7 cells transfected with Ctrl or MYC overexpressing plasmids, and with Ctrl or MYC siRNAs, respectively. e. Western blotting analysis for p53 and c-Myc, and RT-qPCR analysis for miR-192-5p in HepG2 treated with Nutlin-3a (10uM) and/or transfected with MYC. f. miR-192-5p expression levels were shown in c-Myc activation group and c-Myc non-activation group defined in Supplementary Fig. S6c-d in two cohorts. In Cohort 2, mir-192 expression levels were further shown in subgroups with different status of p53 mutation, mir-192 promoter methylation and c-Myc activation. Student t-test was performed. *, p\u003c0.05.\n","description":"","filename":"F5.png","url":"https://assets-eu.researchsquare.com/files/rs-56821/v3/e1f4ee5c82131d9d2c2ee803.png"},{"id":3807133,"identity":"f29a71a2-89bb-431d-8fda-c95b2ed1eced","added_by":"auto","created_at":"2020-11-24 21:26:43","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":351766,"visible":true,"origin":"","legend":"Overproduced lactate from HCC cells with miR-192-5p loss activated ERK pathway in environmental non-tumor cells.\na. Western blotting analysis for pERK and total ERK in LX2, THP1 and HL7702 cells treated with 0, 10, 20mM lactate. b. Western blotting analysis for pERK and total ERK in LX2, THP1 and HL7702 cells upon co-culture with HLF-WT or -192KO cells infected with lentivirus pmiR-Ctrl or pmiR-192 in a co-culture chamber system. c. RT-qPCR analysis of MCT1 and NDRG3 in different HCC and non-HCC cell lines. d. Western blot analysis of pERK and total ERK in LX2 cells transfected with si-NDRG3 or si-MCT1 and treated with different doses of Lactate. e. Western blot analysis of pERK and total ERK in THP1 cells transfected with si-NDRG3 or si-MCT1 and treated with different doses of lactate. Student t-test was performed. *, p\u003c0.05.\n","description":"","filename":"F6.png","url":"https://assets-eu.researchsquare.com/files/rs-56821/v3/0f8e873cb3c7d71362e51793.png"},{"id":3807134,"identity":"ddf9b546-a243-455a-b735-5e8682cee81c","added_by":"auto","created_at":"2020-11-24 21:26:43","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":239651,"visible":true,"origin":"","legend":"Silencing MCT1 and NDRG3 in LX2 or THP1 reduced malignancy and stemness features of co-cultured HCC cells.\na. The schematic diagram of examining malignancy features of HCC cells under a co-culture condition. b. For wound healing assay, scratches were generated in a confluent monolayer HLF cells infected with pmiR-ctrl/RFP or pmiR-192/RFP which were co-cultured with LX2 or THP1 pre-transfected with si-NDRG3 or si-MCT1. c. HLF-192KO infected with pmiR-ctrl/GFP or pmiR-192/GFP were used for spheroid formation assay. Conditioned medium was used for this assay and collected from corresponding HLF cells co-cultured with LX2 pre-transfected with si-NDRG3 or si-MCT1. d. Quantitative data of flow cytometry analysis of CD24+ and CD44+ populations in HLF-192KO cells or HLF-192KO infected with pmiR-192/GFP cells co-cultured with LX2 pre-transfected with si-NDRG3 or si-MCT1. e. Relative levels of MCT1 and NDRG3 in non-tumor samples from patients with different levels of miR-192-5p in their tumors. f. Kaplan–Meier curves of overall survival and time to recurrence in Cohort 1 according to miR-192-5p level in tumors as well as MCT1 or NDRG3 levels in non-tumors. g. The schematic model of the miR-192-5p regulatory pathway in glycolysis and hepatic CSC features. Student t-test was used for (c. d. and e). *, p\u003c0.05.\n","description":"","filename":"7.png","url":"https://assets-eu.researchsquare.com/files/rs-56821/v3/8cb78f76225d81c106f03faa.png"},{"id":13620898,"identity":"461c9dc1-cac5-4a27-86c7-4be6d6002626","added_by":"auto","created_at":"2021-09-17 07:07:39","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1620233,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-56821/v3/bac7a9f2-a738-4eda-a646-f3c032243bf8.pdf"},{"id":3807127,"identity":"e8bb7ec7-0f37-4a76-af3f-3fd58821ecf4","added_by":"auto","created_at":"2020-11-24 21:26:42","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":4842854,"visible":true,"origin":"","legend":"","description":"","filename":"manuscriptSup1118.docx","url":"https://assets-eu.researchsquare.com/files/rs-56821/v3/30e0f5fbf290a3970a56aa35.docx"}],"financialInterests":"","formattedTitle":"\u003cp\u003eLoss of MiR-192-5p Initiates a Hyperglycolysis and Stemness Positive Feedback in Hepatocellular Carcinoma.\u003c/p\u003e","fulltext":[{"header":"Background","content":"\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eAbnormal cancer metabolism is one of the ten cancer hallmarks. Studies of cancer metabolism have revealed the important roles of metabolic reprogramming in cancer cells for proliferation, metastasis and drug resistance[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_1\"\u003e1\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_2\"\u003e2\u003c/a\u003e]. In tumors, a small population of cancer cells exhibit a high capacity of self-renewal and tumor initiation in NOD/SCID mice, which are referred to as cancer stem cells (CSCs). Currently, CSCs are considered to be responsible for tumor initiation, therapy resistance as well as recurrence, and rising attentions have been placed on CSC metabolism [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_3\"\u003e3-5\u003c/a\u003e]. Recent studies indicate that CSCs have different metabolic properties when compared to the tumor bulk and that metabolic plasticity is dependent on the intrinsic demand of nutrients as well as the surrounding environment [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_6\"\u003e6-9\u003c/a\u003e]. CSCs exhibit distinct metabolic phenotypes that include low mitochondrial respiration, high glycolytic activity, and high fatty acid oxidation based on tumor types and CSC biomarkers as well as isolation methods [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_6\"\u003e6\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_10\"\u003e10\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_11\"\u003e11\u003c/a\u003e]. Exploiting metabolic vulnerability of CSCs may provide new effective cancer therapies to diminish tumor recurrence and metastasis.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003ePrimary liver cancer is the fourth lethal neoplasm worldwide, 90% of which are hepatocellular carcinoma (HCC) [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_12\"\u003e12\u003c/a\u003e]. Hepatic CSCs are considered as one of the determining factors for HCC carcinogenesis and recurrence and multiple hepatic CSC biomarkers such as EpCAM, CD133, CD90, CD44, and CD24, have been used to enrich tumorigenic CSCs in both HCC cell lines and primary HCC tissues [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_13\"\u003e13-19\u003c/a\u003e] . Although gene signatures and regulatory pathways of CSC\u003csup\u003e+\u003c/sup\u003eHCCs were thoroughly investigated, their metabolic properties remained largely unknown except for few recent studies. Three studies revealed a high glycolysis rate and a low oxygen consumption rate (OCR) in CD133\u003csup\u003e+\u003c/sup\u003e hepatic CSCs[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_20\"\u003e20\u003c/a\u003e], the importance of fatty acid oxidation in CD133\u003csup\u003e+\u003c/sup\u003e CSCs[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_21\"\u003e21\u003c/a\u003e], and increased lipid metabolites associated with stearoyl-CoA-desaturase in EpCAM\u003csup\u003e+\u003c/sup\u003e HCCs[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_22\"\u003e22\u003c/a\u003e]. Hence, it is essential to systematically investigate the key metabolic features shared among different CSC biomarker positive HCCs.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eMiR-192-5p is the 2\u003csup\u003end\u003c/sup\u003e most abundant miRNA in the liver [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e]. Its genetic silencing frequently occurs in many groups of hepatic CSC marker-positive HCCs such as EpCAM\u003csup\u003e+\u003c/sup\u003e, CD90\u003csup\u003e+\u003c/sup\u003e, CD133\u003csup\u003e+\u003c/sup\u003e, CD44\u003csup\u003e+\u003c/sup\u003e and CD24\u003csup\u003e+\u003c/sup\u003eHCCs, as well as pluripotency marker-positive HCCs[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e]. HCC cells with miR-192-5p loss display a highly invasive phenotype and essential CSC features partially mediated via the p53/miR-192-5p/PABPC4 pathway [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23-25\u003c/a\u003e]. Thus, investigating key metabolic features in HCCs with miR-192-5p loss might allow us to explore the essential metabolic properties among multiple different groups of CSC\u003csup\u003e+\u003c/sup\u003eHCCs.\u003c/span\u003e\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eOmics dataset\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eA total of five datasets from two HCC cohorts with 548 cases were used. In Cohort 1, there were a total of 176 Asian HCC cases, among which 22 cases had available metabolomics data and all cases had miRNA transcriptome (GSE6857) and mRNA transcriptome data (GSE14520) [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_22\"\u003e22\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_26\"\u003e26-29\u003c/a\u003e] (Fig. 1a). In Cohort 2, there were 372 HCC cases with different races. Their miRNA and mRNA sequencing data were downloaded from The Cancer Genome Atlas (TCGA) and available \u003cem\u003eTP53 \u003c/em\u003emutation status and \u003cem\u003eMYC\u003c/em\u003e amplification information in 240 HCC cases were collected from \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"http://www.cbioportal.org/index.do\"\u003ewww.cbioportal.org/index.do\u003c/a\u003e.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eCell lines and miR-192-5p knockout HCC cells \u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eHuman liver cancer cell lines including HLF, HLE, Huh7 and HepG2 cells; human embryonic kidney HEK293T cells; human hepatic stellate cell (HSC) line LX2 cells; human hepatocyte line HL7702 cells; and human leukemic monocyte cell line THP1 (with macrophage activation) were used in this study as described before[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_30\"\u003e30\u003c/a\u003e]. HLF, HLE and Huh7 were originally from Japanese Collection of Research Biosources Cell Bank (JCRB). HepG2, 293T and THP1 were from American Type Culture Collection (ATCC), and LX2 and HL7702 were from Chinese Academy of Sciences (Shanghai, China). Cells were cultured in Dulbecco\u0026rsquo;s modified Eagle\u0026rsquo;s medium (DMEM) supplemented with 10% FBS and 100\u0026nbsp;U/mL penicillin\u0026ndash;streptomycin, 1% L-glutamine. Cells were grown at 37\u0026deg;C in a 5% CO2 incubator and regularly tested to be free of mycoplasma contamination.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eMiR-192-5p knockout HCC cells were generated using the CRISPR/Cas9 system. Two single-guide RNAs (sgRNAs) were designed via an online web tool (\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"http://crispr.mit.edu\"\u003ehttp://crispr.mit.edu\u003c/a\u003e) to target the precursor mir-192. The synthesized sequences were cloned into pSpCas9(BB)-2A-Puro (PX459) vector, termed as PX459-sgMiR192. 48 hours after PX459-sgMiR192 transfection, HCC cells underwent puromycin selection (1.5\u0026micro;g/ml) for 5 days. Next, 1000 collected cells were seeded in a 10-cm dish and cultured for 10-12 days to select miR-192-5p knockout single cell clones. Positive clones were identified via PCR examination and sequencing validation. The sequences for guide RNAs and PCR primers are listed in Supplementary Table S1.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003ePlasmids, siRNAs and cell treatments\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eLentiviral constructs pre-hsa-miR-control (pmiR-control) and pre-hsa-miR-192-5p (pmiR-192) with GFP as well as miRzip-control and miRZip-192-5p (miRZip-192) were obtained from SBI Biosciences and stored in our laboratory. pmiR-control with RFP (pmiR-ctrl/RFP) and pmiR-192 with RFP (pmiR-192/RFP) were constructed by replacing GFP with Ds-Red between the HindIII and NotI sites. Lentiviruses were packaged with plasmids psPAX2 and pMD2.G (Addgene) in 293T cells. For infection, 5 MOI of each lentivirus was used for all our studies.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003epMiR-Report-control (Luc-Ctrl), pmiR-192 reporter (192pos reporter), and miRZip-192 reporter (Zip192 reporter) were constructed before[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e]. The miR-192-5p binding regions in the 3\u0026rsquo;UTR or coding region of \u003cem\u003eGLUT1\u003c/em\u003e, \u003cem\u003ePFKPB3\u003c/em\u003e, and \u003cem\u003eMYC \u003c/em\u003ewere inserted into the monoclonal sites (HindIII/SpeI) of the pMiR-Report plasmid to generate the corresponding luciferase reporters. The pT3-EF1\u0026alpha;-c-Myc vector was originally modified from pT3-EF1\u0026alpha; for c-Myc overexpression[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_31\"\u003e31\u003c/a\u003e]. pGL-miR-192-5p (pGL-192) plasmids were constructed by inserting different lengths of the miR-192 promoter region to the monoclonal site (KpnI/XhoI) of pGL4.20-basic (Promega).\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cem\u003ePFKFB3\u003c/em\u003e siRNAs, \u003cem\u003eGLUT\u003c/em\u003e1 siRNAs and negative control siRNA were purchased from RiboBio Co, Guangzhuo, China.\u003cem\u003e TP53 \u003c/em\u003esiRNAs, \u003cem\u003eMYC\u003c/em\u003e siRNAs, \u003cem\u003eNDRG3\u003c/em\u003e siRNAs, and\u003cem\u003e MCT1\u003c/em\u003e siRNAs were purchased from GenePharma Co., Shanghai, China. Lipofectamine 2000 (Invitrogen) reagent and Rfect siRNA Transfection Reagent (BIO-TRAN) were used for the transfections of plasmids and siRNAs, respectively. The detailed information of all primers and sequences is listed in Supplementary Table S1. \u0026nbsp;\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eNutlin-3a (Selleck Chemicals) was dissolved in DMSO. When indicated, HCC cells were treated for 24 hours with 10 \u0026micro;M of Nutlin-3a. 2-DG (Sigma-Aldrich) was dissolved in ddH\u003csub\u003e2\u003c/sub\u003eO and HCC cells were treated for 72 hours with 5-10 \u0026micro;M of 2-DG.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eExtracellular acidification rate (ECAR) \u003c/strong\u003e\u003cstrong\u003eand oxygen consumption rate (OCR) \u003c/strong\u003e\u003cstrong\u003eassays\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eIn the XF8 Extracellular Flux Analyzer (Seahorse Bioscience), the ECAR and OCR were measured using the Seahorse XF Glycolysis Stress Test Kit (Agilent Technologies) and Seahorse XF Cell Mito Stress Test Kit (Agilent Technologies), respectively. Experiments were performed according to the manufacturer\u0026rsquo;s instructions. For the ECAR measurement, cells (1 \u0026times; 10\u003csup\u003e4\u003c/sup\u003e cells/well) were plated in an XF8 cell-culture microplate for 12 hours and the culture medium was then replaced with XF assay medium supplemented plus 2\u0026nbsp;mM L-glutamine. After one-hour incubation at 37\u0026deg;C, the ECAR was measured by the sequential addition of glucose (10\u0026nbsp;mM), oligomycin (3.5\u0026nbsp;\u0026mu;M), and 2-deoxyglucose (2DG, 80\u0026nbsp;mM) in an XF8 flux analyzer. For the OCR measurement, cells (1 \u0026times; 10\u003csup\u003e4\u003c/sup\u003e cells/well) were plated in an XF8 cell-culture microplate for 12 hours and then culture medium was replaced with XF assay medium supplemented with 10\u0026nbsp;mM D-glucose, 1\u0026nbsp;mM sodium pyruvate and 2\u0026nbsp;mM L-glutamine. The OCR was measured by sequential addition of oligomycin (2\u0026nbsp;\u0026mu;M), FCCP (0.5\u0026nbsp;\u0026mu;M), and antimycin A (1\u0026nbsp;\u0026mu;M). Data were analyzed by Seahorse XF Wave software. The results were normalized to the cell number.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eCo-culture system \u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eTwo types of co-culture systems were used. In the first, HCC cells infected with corresponding RFP or GFP lentivirus were co-cultured with HL7702 or LX2 cells in the same dish. In this system, cells that displayed red or green under a fluorescence microscope or flow cytometry were regarded as HCC cells. The second co-culture system was based on a chamber system using Polyester (PET) Membrane Tissue Culture Plate Insert with 0.4 \u0026mu;m pores (JET biofil) according to the manufacturer\u0026rsquo;s protocol. Generally, LX2 or HL7702 were placed in the bottom layer of a 6-well plate while HCC cells were placed in the cell inserts of this 6-well plate. In addition, HCC cells were also co-cultured with THP1, a suspension cell line.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eRNA isolation, quantitative real-time PCR, and Western blot\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eTotal RNA was isolated using TRIzol (Invitrogen) following the manufacturer's instructions, and 1 \u0026mu;g of total RNA was reverse transcribed into cDNA using PrimeScript\u003csup\u003eTM\u003c/sup\u003e RT Reagent Kit with gDNA Eraser (TaKaRa). Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed with TB Green\u003csup\u003eTM\u003c/sup\u003e Premix EX Taq\u003csup\u003eTM\u003c/sup\u003e II (Tli RNaseH Plus) (TaKaRa). The expression of mature miRNAs was measured using TaqMan MiRNA Assays as described previously[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_28\"\u003e28\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_32\"\u003e32\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_33\"\u003e33\u003c/a\u003e]. RNU6B was used as the reference gene for miRNAs, while 18S was the reference for measured mRNAs. Primers for\u003cem\u003e TACSTD1, CD133, CD90, CD44, CD24,\u003c/em\u003e and \u003cem\u003eUGT2B7 \u003c/em\u003ewere used as before[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e]\u003cem\u003e. \u003c/em\u003ePrimer sequences for \u003cem\u003eGLUT1\u003c/em\u003e, \u003cem\u003eHK2, PFKP, PFKFB3, ALDOA, ENO2, PKM2, LDHA, MCT4, MCT1, NDRG3 \u003c/em\u003eand\u003cem\u003e CYP1A2\u003c/em\u003e are listed in Supplementary Table S1. The experiments were performed in triplicate.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eFor the Western blot assay, cells were lysed and processed as previously described [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_28\"\u003e28\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_33\"\u003e33\u003c/a\u003e]. The membranes with transferred protein extracts were incubated with the indicated primary antibodies and secondary antibodies conjugated to horseradish peroxidase for enhanced chemiluminescence detection of the signals (Amersham, Arlington Height, IL). The detailed information of all antibodies is listed in Supplementary Table S2.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eNon-targeted metabolomics study\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eThe non-targeted metabolomics was performed by Metabo-Profile Biotechnology (Shanghai) Co., Ltd. Cell lysates and cell culture medium were collected from cells and used for metabolomics analysis. Metabolites were identified and quantified from gas chromatography/time-of-flight mass spectrometry (GC/TOF-MS) data, which was performed in the workflow of ADAP-GC 2.0. Each experiment was performed in triplicate.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003e2-NBDG uptake, glucose measurement and lactate measurement \u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eGlucose uptake of HCC cells was quantified by flow cytometry using the fluorescent D-glucose derivate 2-(N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino)-2-deoxy-D-glucose (2-NBDG). Briefly, the corresponding cells were cultured with 2-NBDG at 37\u0026deg;C for 1 h. The uptake of 2-NBDG was measured by flow cytometry. A glucose measurement kit (Shanghai Rongsheng Biotech Co., Ltd.) was used to measure the glucose concentration in cell culture medium. Extracellular lactate production was measured using a lactate assay kit (Nanjing Jiancheng Bioengineering Institute). All these assays were performed according to the manufacturer's protocols.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eSphere formation assay and flow cytometry analysis\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eSingle-cell suspensions of 1000 cells were seeded in 6-well Ultra-Low Attachment Microplates (Corning, Corning, NY) for spheroid assays. The number of spheroids was measured 12 days after seeding.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eFor flow cytometry analysis, cultured cells were trypsinized, washed, and resuspended in phosphate-buffered saline plus 0.5% bovine serum albumin. They were incubated with Allophycocyanin (APC)-conjugated antibodies on ice for 20 mins in the dark. Data were collected with a FACS Calibur flow cytometer (BD Biosciences) and analyzed using FlowJo software (Tree Star). The detailed information for all antibodies is listed in Supplementary Table S2.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eLuciferase reporter assay\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eThe pMiR-Report plasmids (with or without mature miR-192-5p binding sites) were transfected together with pRL-CMV vector containing Renilla luciferase. The Firefly and Renilla luciferase activities were measured 24 hours after transfection using Dual-Luciferase Reporter Assay (Promega, CA) with a PerkinElmer luminometer. Each experiment was performed in triplicate and repeated at least three times.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eTo examine the promoter activity of miR-192-5p, HCC cells (24-well plate) were transfected with 300 ng pT3-EF1\u0026alpha; vector or pT3-EF1\u0026alpha;-cMyc vector, or 50 nM siRNAs (negative control or \u003cem\u003eMYC\u003c/em\u003e siRNAs) on the first day and then transfected with 300ng pGL-192 promoter constructs and 10 ng pRL-CMV on the next day. The Firefly and Renilla luciferase activities were then measured 24 hours after transfection.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eStatistics\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eSpearman\u0026rsquo;s rank correlation was performed to identify metabolites associated to miR-192-5p. Student\u0026rsquo;s t-test and Mann-Whitney rank test were used for statistical analysis of comparative data between groups. Two-way ANOVA was used to compare the glucose consumption, lactate production and cell viability of HCC cells at different time points. Hierarchical clustering analysis was performed by GENESIS software version 1.7.6 developed by Alexander Sturn (IBMT-TUG, Graz, Austria). Pearson correlation was used to identify genes correlated with miR-192-5p. Kaplan\u0026ndash;Meier survival analysis was used to compare patient survival based on prediction results using GraphPad Prism V7.0 (San Diego, CA), and the p-value was generated by the Cox\u0026ndash;Mantel log-rank test. Gene set enrichment analysis (GSEA) in the Molecular Signatures Database was performed using GSEA V3.0. All p-values were 2-sided. A p-value of \u0026le;0.05 was regarded as statistically significant.\u003c/span\u003e\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eGlycolysis-related metabolites and genes were highly expressed in HCC cases with low miR-192-5p expression \u0026nbsp;\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eWe first investigated the metabolic features in HCCs with low miR-192-5p levels using HCC Cohort 1 with available metabolome and transcriptome data (Fig. 1a). We performed an integration analysis of miR-192-5p with the global metabolome in tumor tissues among 22 HCC patients, and found that 17 metabolites were significantly correlated with miR-192-5p (|r-value| \u0026gt;0.4, Fig. 1b, Supplementary Table S3). Among them, 7 metabolites presented |r-value| \u0026gt;0.5 and three of them were glycolysis-related metabolites, i.e., Glucose-6-Phosphate (G6P), Fructose-6-Phosphate (F6P), and Nicotinamide Adenine Dinucleotide Phosphate (NADPH). Meanwhile, 652 genes were significantly correlated with miR-192-5p with |r-value| \u0026gt;0.4, revealed by an integration analysis of miR-192-5p with mRNA transcriptome in tumor tissues among 176 HCC patients (Fig. 1c). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using these genes revealed 13 enriched metabolic features (p\u0026lt;0.001), three of which were associated with glycolysis and glycolysis-related pathways. These results suggest an altered glycolytic feature in HCC cases with low miR-192-5p expression.\u0026nbsp;\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eAvailable glycolysis-related metabolites and genes in our profiling data (Fig. 1d) were then compared between HCCs with high miR-192-5p levels (termed HCC\u003csup\u003e192High\u003c/sup\u003e) and HCCs with low miR-192-5p levels (termed HCC\u003csup\u003e192Low\u003c/sup\u003e), based on a miR-192-5p median cut-off in HCC tumors. Levels of G6P, F6P, and NADPH were significantly higher in tumors from HCC\u003csup\u003e192Low\u003c/sup\u003e patients than HCC\u003csup\u003e192High \u003c/sup\u003epatients, while no difference was found in their non-tumor tissues (Fig. 1e, Supplementary Fig. S1a). Consistently, many genes coding for key glycolytic enzymes such as \u003cem\u003eGLUT1, HK2, PFKFB3, PFKP, and PKM2\u003c/em\u003e were significantly upregulated in HCC\u003csup\u003e192Low\u003c/sup\u003e tumors compared to HCC\u003csup\u003e192High\u003c/sup\u003e tumors (Fig. 1e) but showed negligible alteration in their non-tumor tissues (Supplementary Fig. S1b). \u003cem\u003eMCT1\u003c/em\u003e was used as a negative control due to its main role in lactate import, but not in glycolysis[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_34\"\u003e34\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_35\"\u003e35\u003c/a\u003e]. These results demonstrate a hyperglycolytic metabolic feature in HCC cases with low miR-192-5p level.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eWe further investigated the hyperglycolytic feature in CSC\u003csup\u003e+\u003c/sup\u003e HCC cases, i.e., cases with the top quartile expression of CSC markers as previously defined [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e]. Hierarchical clustering analysis with glycolytic genes in Cohort 1 revealed two HCC subgroups with different expression levels of glycolytic genes. Consistently, in HCC subgroup with high expression levels of glycolysis-related genes, miR-192-5p level was low while various groups of CSC\u003csup\u003e+\u003c/sup\u003e HCC cases were enriched (Fig. 1f). Statistical analysis also showed that glycolytic genes were expressed at significantly higher levels in various groups of CSC\u003csup\u003e+\u003c/sup\u003e HCCs than in CSC\u003csup\u003e- \u003c/sup\u003eHCCs, but no difference was observed in the comparisons of their non-tumor tissues (Supplementary Fig. S1c-d). Comparable data were observed in Cohort 2 with 372 HCC patients (Supplementary Fig. S1e). Together, these data indicated that the hyperglycolytic feature was present in various groups of CSC\u003csup\u003e+\u003c/sup\u003eHCCs with low level of miR-192-5p.\u0026nbsp; \u0026nbsp;\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eHCC cells with miR-192-5p loss were hyperglycolytic.\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eWe next investigated the role of miR-192-5p in regulating glycolysis. In Huh7, glycolytic genes with significant differential expression between HCC\u003csup\u003e192High\u003c/sup\u003e and HCC\u003csup\u003e192Low\u003c/sup\u003e patients were examined and eight out of nine genes showed a significant up-regulation after suppressing miR-192-5p (Supplementary Fig. S2a). Meanwhile, Huh7 cells with suppressed miR-192-5p exhibited a distinctly increased ECAR and a reduced OCR (Fig. 1g). The extracellular acid produced by cells is derived from lactate produced by glycolysis and CO\u003csub\u003e2\u003c/sub\u003e produced during respiration. OCR mainly represents mitochondrial respiration. Therefore, the increased ECAR in Huh7 cells with suppressed miR-192-5p was mainly due to lactate produced from glycolysis but not CO\u003csub\u003e2\u003c/sub\u003e from mitochondrial respiration. Consistently, over-expressed miR-192-5p in HLF and HLE cells lowered the extracellular acid production from glycolysis as shown by a reduced ECAR but an increased OCR (Fig. 1h). Most examined glycolytic genes were significantly reduced by miR-192-5p overexpression in both HLF and HLE cells (Supplementary Fig. S2b). These results demonstrate an important role of miR-192-5p in modulating a Warburg-like effect in HCC cells.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eTo better elucidate the role of miR-192-5p in regulating glycolysis, we established miR-192-5p knockout (termed 192KO) clones from two human HCC cell lines HLF and HLE. A 69bp DNA fragment was deleted in the 192KO mixture clones as well as selected single 192KO clones (Supplementary Fig. S2c-d). One 192KO single clone from each HCC cell lines was used and miR-192-5p expression was undetectable in HLE-192KO and HLF-192KO cells (Fig. 2a). As a control, the expression of miR-194, a nearby miRNA of miR-192-5p, was not affected. As expected, HLF-192KO cells displayed significantly increased CSC features, such as increased populations of CD44\u003csup\u003e+\u003c/sup\u003e, CD24\u003csup\u003e+\u003c/sup\u003e and EpCAM\u003csup\u003e+\u003c/sup\u003e CSCs (Fig. 2b, Supplementary Fig. S2e); increased mRNA levels of multiple CSC biomarkers and reduced expression of a differentiation-related gene \u003cem\u003eCYP1A2\u003c/em\u003e (Supplementary Fig. S2f); and enlarged and more spheroid formation (Supplementary Fig. S2g). HLE-192KO cells displayed increased CSC features at a moderate level (Fig. 2b, Supplementary Fig. S2e-f). Consistently, these two 192KO lines also showed the hyperglycolytic features. As shown in Fig 2c, five glycolytic enzymes, i.e., GLUT1, HK2, PFKFB3, ALDOA, and PKM2, as well as c-Myc presented higher levels in 192KO cells than in wild-type cells. 192KO lines also exhibited increased ECARs but decreased OCRs (Fig. 2d), indicating that miR-192-5p loss largely increased the glycolysis-related extracellular acidification. Consistent data were also noticed in other 192KO clones of both HLF and HLE HCC cell lines (Supplementary Fig. S3a-b).\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eFurthermore, overexpressed miR-192-5p in HLF-192KO cells significantly reduced CSC features and lactate accumulation in the culture medium (Fig. 2e). As a control, the intracellular lactate remained unchanged. Meanwhile, over culturing time, lactate gradually accumulated in the medium and was significantly higher in both HLF-192KO and HLE-192KO cells compared to their corresponding wild-type cells, which could be lowered by overexpressed miR-192-5p (Fig. 2f). A lower pH value was also observed in 192KO cells indicated by the orange/yellow medium vs. the pink medium of wild-type cells at 72 hours after seeding. Comparable data on lactate production were seen in other HLF and HLE 192KO clones (Supplementary Fig. S3c). Consistent data were also obtained through the detection of lactate using non-targeted metabolomics in HLF cells with different expression of miR-192-5p in both internal cells and culture medium (Supplementary Fig S3d). Together, miR-192-5p loss in HCC cells led to a hyperglycolytic phenotype.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eHCC cells with miR-192-5p loss had high glucose consumption\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eWe further examined glucose consumption among HCC cells with different levels of miR-192-5p as well as between HCC cells and their co-cultured non-tumor cells. As shown in Fig. 3a, both HLF-192KO and HLE-192KO cells exhibited significantly higher glucose consumption than HCC cells overexpressing miR-192-5p. In Huh7 cells, suppressing miR-192-5p increased their glucose usage (Supplementary Fig. S4a). Consistently, 192KO cells were more sensitive after exposure to 2-DG, a glucose analog, as shown by the significantly reduced cell viability compared to cells with miR-192-5p expression (Fig. 3b). Comparable data were obtained in HuH7 cells (Supplementary Fig. S4b).\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eIn co-culture systems of HCC cells with LX2, HL7702, and THP1, we further compared their glucose uptake via 2-NBDG uptake assay. HLF HCC cells infected with pmiR-ctrl/RFP or pmiR-192/RFP lentiviruses were used, and red fluorescent labeling efficiency was nearly 100% (Supplementary Fig. S4c). In this system, with or without co-culturing with other cells, HLF-192KO cells consistently showed higher 2-NBDG uptake than wild-type cells (Fig. 3c-d). In contrast, LX2 and HL7702 in co-culture with HLF-192KO cells exhibited lower 2-NBDG uptake compared to those in co-culture with HLF-WT cells. Moreover, forced-expression of miR-192-5p in HLF-192KO cells reduced the 2-NBDG uptake in HLF cells but increased 2-NBDG uptake in LX2 and HL7702 cells in the co-culture system. The alteration of 2-NBDG uptake was not observed in THP1 from our co-culture system (Fig. 3d, Supplementary Fig. S4d). Similar data were seen in HLE cells as well as in HLE cells co-cultured with LX2 and HL7702 (Supplementary Fig. S4e-f). These results demonstrate that HCC cells with loss of miR-192-5p actively utilize glucose from their environment to ensure a hyperglycolysis status.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eGLUT1, PFKFB3 and c-Myc were miR-192-5p \u003c/strong\u003e\u003cstrong\u003e\u003cem\u003ebona fide\u003c/em\u003e targets and contributed to glycolytic and stemness features of HCC cells\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eTo investigate the target genes of miR-192-5p in regulating glycolysis flow, we assessed genes negatively correlated with miR-192-5p in 176 HCC cases (r\u0026lt;-0.3) and significantly up-regulated in HCC\u003csup\u003e192Low\u003c/sup\u003e tumors versus HCC\u003csup\u003e192High\u003c/sup\u003e tumors (log\u003csub\u003e2\u003c/sub\u003efold \u0026gt;0.2, p\u0026lt;0.01). Among these 554 genes, two main groups were observed. One group contained genes related to cell migration as we previously reported [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e]. The other group included eight glycolysis-related genes (Fig. 4a) and three of them (\u003cem\u003eGLUT1, HK2\u003c/em\u003e, and \u003cem\u003ePKM2\u003c/em\u003e) were reported targets of c-Myc, an important regulator of glycolysis[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_36\"\u003e36\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_37\"\u003e37\u003c/a\u003e].\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eNext, using TargetScan and manual miRNA target prediction, we found that four of these glycolytic genes (\u003cem\u003ePFKFB3, GLUT1, MCT4\u003c/em\u003e, and\u003cem\u003e MYC\u003c/em\u003e) contained miR-192-5p binding sites in their 3\u0026rsquo;UTR and/or coding regions (Fig. 4a). In HLF and HLE cells, miR-192-5p overexpression reduced the protein levels of PFKFB3, GLUT1, and c-Myc, but not that of MCT4 (Fig. 4b). Further, the predicted miR-192-5p binding regions of these three genes were cloned into a luciferase reporter and forced expression of miR-192-5p reduced the luciferase activities when the wild-type sequences for \u003cem\u003ePFKFB3\u003c/em\u003e and\u003cem\u003e GLUT1\u003c/em\u003e as well as the #2 binding site of \u003cem\u003eMYC\u003c/em\u003e were present (Fig. 4c). These effects were significantly reduced when the corresponding miR-192-5p binding sites were mutated. Moreover, silencing PFKFB3, GLUT1, or c-Myc with 2 siRNAs for each gene reduced ECAR in both HLF-WT and HLF-192KO cells (Fig. 4d, Supplementary Fig. S5a). In HLF-192KO cells, silencing PFKFB3, GLUT1 or c-Myc notably reduced the ECAR rate to a level similar to that of HLF-WT cells with silencing of these genes. Comparable data were also noticed in HCC patients from Cohorts 1 and 2. \u003cem\u003ePFKFB3\u003c/em\u003e, \u003cem\u003eGLUT1\u003c/em\u003e, and \u003cem\u003eMYC\u003c/em\u003e presented higher levels in HCC\u003csup\u003e192Low\u003c/sup\u003e tumors compared to HCC\u003csup\u003e192High\u003c/sup\u003e tumors (Supplementary Fig. S5b-c). These results indicate that PFKFB3, GLUT1, and c-Myc are miR-192-5p targets and are involved in the hyperglycolysis caused by miR-192-5p loss.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003ePFKFB3, GLUT1, and c-Myc were reported to maintain stemness features in cancer at certain levels [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_38\"\u003e38-41\u003c/a\u003e]. Consistently, si-\u003cem\u003ePFKFB3,\u003c/em\u003e si-\u003cem\u003eGLUT1\u003c/em\u003e or double knockdown led to reduced levels of four CSC biomarkers, i.e., \u003cem\u003eCD44\u003c/em\u003e, \u003cem\u003eCD24\u003c/em\u003e, \u003cem\u003eEpCAM\u003c/em\u003e and \u003cem\u003eCD90\u003c/em\u003e, as determined by RT-qPCR (Fig. 4e). Flow cytometry analysis also showed that si-\u003cem\u003ePFKFB3\u003c/em\u003e and \u003cem\u003eGLUT1\u003c/em\u003e reduced the populations of CD44\u003csup\u003e+\u003c/sup\u003e and CD24\u003csup\u003e+\u003c/sup\u003e CSCs (Fig. 4f). Meanwhile, si-\u003cem\u003eMYC\u003c/em\u003e seemed to only reduce CD44\u003csup\u003e+\u003c/sup\u003e CSCs moderately, but not the CD24\u003csup\u003e+\u003c/sup\u003e CSC population (Supplementary Fig. S5d). Together, these data demonstrate that three glycolytic regulators, PFKFB3, GLUT1 and c-Myc were \u003cem\u003ebona fide\u003c/em\u003e targets of miR-192-5p, and they contributed to both hyper-glycolysis and CSC features of HCCs caused by loss of miR-192-5p.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003ec-Myc suppressed miR-192-5p transcription, ensuring a positive feedback of high c-Myc/low miR-192-5p in hyperglycolytic CSC\u003csup\u003e+\u003c/sup\u003eHCCs\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eIn our previous miRNA profiles of tumors and non-tumors from a hydrodynamic injection HCC FVB mouse model[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_31\"\u003e31\u003c/a\u003e], miR-192-5p expression was significantly reduced in c-Myc-induced HCCs (Fig. 5a). Further, in a hydrodynamic injection HCC ICR mouse model, miR-192-5p level was also reduced \u0026gt;100 times in c-Myc induced HCCs but was not much in Ras-induced HCCs when compared to corresponding non-HCC liver tissues. In four different HCC cell lines, si-\u003cem\u003eMYC\u003c/em\u003e led to an increased expression of miR-192-5p, while forced expression of c-Myc reduced the level of miR-192-5p (Fig. 5b). These data indicate that c-Myc might regulate miR-192-5p transcription.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eConsistently, among four different lengths of miR-192-5p promoter regions, the -266 nt to +186 nt region showed the strongest promoter activity (Fig. 5c) and the miR-192-5p promoter activity (-266 nt to +186 nt) was reduced by exogenous c-Myc, while enhanced by si-\u003cem\u003eMYC\u003c/em\u003e (Fig. 5d). Wild-type p53 could bind to the miR-192-5p promoter region and induce its expression[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_42\"\u003e42\u003c/a\u003e]. Consistently, in HepG2 cells with wild-type \u003cem\u003eTP53,\u003c/em\u003e the expression of miR-192-5p was induced by p53 via exposure to Nutlin-3a (an MDM2 antagonist to stabilize p53) and reduced by silencing of p53 (Fig. 5e, Supplementary Fig. S6a-b). Over-expressed c-Myc significantly suppressed miR-192-5p expression in HCC cells with either activated p53 or silenced p53, indicating that c-Myc-mediated miR-192-5p down-regulation was independent on p53.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eComparable data were noticed in HCC patients. In both HCC cohorts, hierarchical clustering analysis revealed two subgroups with distinct c-Myc activation status based on 76 c-Myc target genes from the online Human MYC Targets Profiler (Supplementary Fig. S6c-d). In Cohort 1, miR-192-5p expression in the c-Myc activation subgroup was significantly lower than that in c-Myc non-activation subgroup (Fig. 5f). In Cohort 2, miR-192-5p expression was always significantly lower in each c-Myc activation subgroup than in the corresponding non-activation subgroup, which was independent of \u003cem\u003eTP53\u003c/em\u003e mutation and mir-192 promoter methylation (Fig. 5f-g). In addition, different groups of CSC\u003csup\u003e+\u003c/sup\u003eHCCs consistently presented a low level of miR-192-5p, a high level of c-Myc activation and high frequency of \u003cem\u003eMYC\u003c/em\u003e amplification (Supplementary Fig. S7). Together, c-Myc suppressed miR-192-5p transcription, which led to a positive feedback of high c-Myc/low miR-192-5p in CSC\u003csup\u003e+\u003c/sup\u003eHCC cells with glycolytic feature.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003e\u003cstrong\u003eOverproduced lactate from CSC\u003csup\u003e+\u003c/sup\u003eHCCs activated the ERK pathway in environmental non-tumor cells, and this effect further \u003c/strong\u003e\u003cstrong\u003eincreased HCC cell stemness and malignancy features\u003c/strong\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eAs the end product of glycolysis, the continuously produced lactate in hyperglycolytic miR-192-5p-loss HCC cells might affect their environment and contribute to HCC malignancy. The transport of lactate across the plasma membrane is mainly catalyzed by MCT1 and MCT4, with MCT1 typically involved in the import while MCT4 in export of lactate[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_34\"\u003e34\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_35\"\u003e35\u003c/a\u003e]. In HCC patients, the expression ratio of \u003cem\u003eMCT1 vs. MCT4\u003c/em\u003e showed no difference between tumor and non-tumor tissues of Cohort 1 but was significantly higher in non-tumor tissues than tumor tissues of Cohort 2 (Supplementary Fig. S8a), indicating the possibility of lactate uptake by environmental non-tumor cells. Lactate could also stabilize NDRG3, which in turn activated the ERK pathway to promote cell malignancy[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_43\"\u003e43\u003c/a\u003e]. Consistently, lactate treatment stimulated ERK phosphorylation noticeably in non-tumor cells of HCC microenvironment, i.e., LX2, THP1 and HL7702 cells (Fig. 6a).\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eIn a chamber co-culture system, pERK level was increased in LX2 and THP1 cells when co-cultured with HLF-192KO cells compared to when co-cultured with HLF-WT (Fig. 6b). Moreover, pERK was further reduced in LX2 and THP1 when exposed to HLF-192KO cells with miR-192-5p overexpression (Fig. 6b). Similar results were observed from LX2 and THP1 cells co-cultured with HLE cells (Supplementary Fig. S8b). However, pERK was unaltered in HL7702 cells in this co-culture system. Thus, HCC cells with miR-192-5p loss could actively affect certain non-tumor cells via increased production of lactate.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eWe have also found that lactate-induced pERK in environmental non-tumor cells partially relied on NDRG3 and MCT1 (Fig. 6c-e). LX2 and THP1 cells expressed relatively high levels of \u003cem\u003eMCT1\u003c/em\u003e and \u003cem\u003eNDRG3 \u003c/em\u003e(Fig. 6c). Silencing NDRG3 or MCT1 in LX2 and THP1 cells reduced the level of lactate-induced pERK (Fig. 6d-e). Low expression levels of MCT1 and NDRG3 in HL7702 were consistent with its minor response to lactate (data not shown).\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eWe then explored the effects of an altered lactate/MCT1/NDRG3/pERK axis in LX2 or THP1 cells on the malignancy features in HCC cells (Fig. 7a). HLF-192KO/RFP cells were co-cultured with LX2 cells pre-transfected with si-Ctrl, or si-\u003cem\u003eNDRG3\u003c/em\u003e, or si-\u003cem\u003eMCT1\u003c/em\u003e (termed LX2\u003csup\u003esi-Ctrl\u003c/sup\u003e, LX2\u003csup\u003esi-\u003cem\u003eNDRG3\u003c/em\u003e\u003c/sup\u003e, and LX2\u003csup\u003esi-\u003cem\u003eMCT1\u003c/em\u003e\u003c/sup\u003e\u003cem\u003e,\u003c/em\u003e respectively). Wound-healing assay of red-fluorescence HCC cells showed that cell migration of HLF-192KO cells was slower upon co-culture with LX2\u003csup\u003esi-\u003cem\u003eNDRG3\u003c/em\u003e\u003c/sup\u003e or LX2\u003csup\u003esi-\u003cem\u003eMCT1\u003c/em\u003e\u003c/sup\u003e than with LX2\u003csup\u003esi-Ctrl\u003c/sup\u003e (Fig. 7b). Consistent data were observed in HLF-192KO cells co-cultured with THP1 (Fig. 7b). Moreover, spheroid assays of HLF-192KO cells were performed in different conditioned medium settings. The number of spheroids of HLF-192KO cells was significantly lower under exposure to conditioned medium from co-culture of HLF-192KO with LX2\u003csup\u003esi-\u003cem\u003eNDRG3\u003c/em\u003e\u003c/sup\u003e or LX2\u003csup\u003esi-\u003cem\u003eMCT1 \u003c/em\u003e\u003c/sup\u003ethan from co-culture of HLF-192KO with LX2\u003csup\u003esi-Ctrl\u003c/sup\u003e (Fig. 7c). CD44\u003csup\u003e+\u003c/sup\u003e and CD24\u003csup\u003e+\u003c/sup\u003e HLF-192KO populations were also significantly reduced when they were co-cultured with LX2\u003csup\u003esi-\u003cem\u003eNDRG3\u003c/em\u003e\u003c/sup\u003e or LX2\u003csup\u003esi-\u003cem\u003eMCT1 \u003c/em\u003e\u003c/sup\u003e(Fig. 7d, Supplementary Fig. S8c). As a control, HLF-192KO with overexpressed miR-192-5p did not exhibit significant alteration of migration, spheroid formation and CSC populations when co-cultured with different LX2 cells (Fig. 7b-d). Consistent data were obtained using a second set of siRNAs for \u003cem\u003eNDRG3\u003c/em\u003e and \u003cem\u003eMCT1\u003c/em\u003e (Supplementary Fig. S9). Therefore, blocking the lactate/ERK pathway in HCC microenvironmental cells suppressed the malignancy and stemness features of HCC cells in co-culture experiments.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eIn both HCC cohorts, patients were divided into four groups based on miR-192-5p expression in their tumor tissues (HCC\u003csup\u003e192Low\u003c/sup\u003e and HCC\u003csup\u003e192High\u003c/sup\u003e, medium cut-off) and levels of NDRG3 and MCT1 in non-tumor tissues (NT\u003csup\u003eHigh_NDRG3 or MCT1\u003c/sup\u003e and NT\u003csup\u003e Low_NDRG3 and MCT1\u003c/sup\u003e, medium cut-offs). There was no expression difference of NDRG3 or MCT1 in non-tumor tissues between HCC\u003csup\u003e192Low\u003c/sup\u003e and HCC\u003csup\u003e192High\u003c/sup\u003e patients (Fig. 7e, Supplementary Fig. S10a). In Cohort 1, patients with HCC\u003csup\u003e192High\u003c/sup\u003e NT\u003csup\u003eLow_NDRG3 and MCT1\u003c/sup\u003e had the best prognosis, as shown by a prolonged time to recurrence and overall survival. In HCC\u003csup\u003e192Low\u003c/sup\u003e subgroup, patients with NT\u003csup\u003eHigh_NDRG3 or MCT1\u003c/sup\u003e had worse prognosis compared to patients with NT\u003csup\u003eLow_NDRG3 or MCT1\u003c/sup\u003e (Fig. 7f). Similar but less significant data were obtained in Cohort 2, which might be due to the limited number of patients (n=49) with available non-tumor mRNA data (Supplementary Fig. S10b). GSEA analysis in HCC\u003csup\u003e192Low\u003c/sup\u003e patients revealed that several stem cell related gene-sets were enriched in patients with NT\u003csup\u003eHigh_NDRG3 or MCT1\u003c/sup\u003e (Supplementary Fig. S10c-d). Together, HCC cells with miR-192-5p loss exhibited a highly malignant feature when they were surrounded by environmental non-tumors with high MCT1 or NDRG3 expression.\u003c/span\u003e\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eHCC is one of the most malignant cancers worldwide. Although a trend of decreasing incidence has been observed for many cancers, the incidence of HCC is still rising in both developing and developed countries[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_44\"\u003e44\u003c/a\u003e]. CSCs are thought to be responsible for tumorigenesis as well as tumor metastasis and targeting CSCs holds the hope of eliminating cancer[\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_3\"\u003e3\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_4\"\u003e4\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_6\"\u003e6\u003c/a\u003e]. In the last decade, researchers\u0026rsquo; attention has been drawn to the field of hepatic CSCs, leading to the identification of multiple hepatic CSC biomarkers. Accordingly, accumulated studies have revealed gene signatures and regulatory signaling pathways of these identified hepatic CSCs [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_18\"\u003e18\u003c/a\u003e]. In this study, we aimed to understand essential metabolic features of hepatic CSCs.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eMiR-192-5p, a liver-abundant miRNA, is functionally important in suppressing the stemness and malignancy features of HCC cells, and its genetic silencing frequently occurs in multiple groups of CSC\u003csup\u003e+\u003c/sup\u003e HCCs, i.e., EpCAM\u003csup\u003e+\u003c/sup\u003e, CD90\u003csup\u003e+\u003c/sup\u003e, CD133\u003csup\u003e+\u003c/sup\u003e, CD44\u003csup\u003e+\u003c/sup\u003e, and CD24\u003csup\u003e+ \u003c/sup\u003eHCCs [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e]. Here, we revealed that the miR-192-5p loss drove HCC cells to a hyperglycolytic metabolism status via targeting c-Myc and two glycolytic enzymes GLUT1 and PFKFB3, and a c-Myc/miR-192-5p positive feedback pathway. Five groups of CSC\u003csup\u003e+\u003c/sup\u003eHCCs presented a hyperglycolytic signature, high level of c-Myc and low miR-192-5p. Mechanistically, loss of miR-192-5p led to increased expression of its target genes, i.e., PFKFB3, GLUT1, and c-Myc, which facilitated a Warburg effect favored by cellular growth and stemness. Moreover, c-Myc also directly suppressed miR-192-5p expression, ensuring a positive feedback to drive hyperglycolytic features in CSC\u003csup\u003e+\u003c/sup\u003eHCCs. In addition, targeting PFKFB3, GLUT1, and c-Myc could also reduce hepatic CSC populations. Thus, hyperglycolysis was a key feature for CSC\u003csup\u003e+\u003c/sup\u003e HCCs and effective methods of targeting hyperglycolytic cells in combination with conventional anticancer methods might effectively eliminate CSCs, leading to a significant reduction of HCC recurrence and metastasis. However, it remains unclear whether liver cells gain their CSC features and hyperglycolysis features at the same time or if one is the consequence of the other. The efforts to answer these questions will further improve our understanding of HCC initiation and may open new avenues to prevent HCC tumor initiation or metastasis.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eThe HCC microenvironment is complex, including but not limited to activated hepatic stellate cells, Kupffer cells, and hypoxia etc, and is believed to actively interact with HCC tumors as well as affect their malignancy status [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_7\"\u003e7\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_45\"\u003e45\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_46\"\u003e46\u003c/a\u003e]. In our study, LX2 and THP1 (representing hepatic stellate cells and macrophages respectively) were used to co-culture with HCC cells in vitro. We found that hyperglycolytic HCC cells with miR-192-5p loss nourished LX2 and THP1 cells with overproduced lactate, which stimulated the pERK pathway in LX2 and THP1 cells via NDRG3/MCT1. Consistently, in the \u003cem\u003ein vitro\u003c/em\u003e co-culture system, LX2 and THP1 cells with high levels of NDRG3 and/or MCT1 facilitated a more malignant phenotype of HCC cells with miR-192-5p loss. \u003cem\u003eIn vivo\u003c/em\u003e, HCC patients with low level of miR-192-5p in their tumor and high level of \u003cem\u003eNDRG3/MCT1\u003c/em\u003e in their non-tumor had the shorter overall survival when compared to other patient subgroups. Thus, hepatic CSC-related metabolism affected the CSC microenvironment via accumulated lactate, and lactate in turn interacted with environmental cellular components (LX2 and THP1) that contributed to HCC stemness features. This was consistent with reported studies demonstrating that the structural and cellular components of the CSC microenvironment played important roles in CSC plasticity, proliferation and invasion [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_8\"\u003e8\u003c/a\u003e, \u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_9\"\u003e9\u003c/a\u003e]. It remains undiscovered in our study how the pERK pathway in THP1 or LX2 cells promote HCC malignancy. It would be interesting to perform an in-depth study on cytokines or exosomes being produced by THP1 and LX2 cells through the lactate/ERK axis, and to explore their roles in promoting the malignancy of tumor cells. Both lactate and LX2/THP1 cells in the hepatic CSC environment contributed to the stemness and malignant features of HCC cells with miR-192-5p loss. Thus, we would need more insight on the interaction of tumor initiating cells and their environment to offer a better solution of cancer therapy regarding CSC elimination. Investigating the CSC-niche metabolism would be also important for the further development of metabolic targeting approaches, which might significantly impact cancer therapy.\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cspan style=\"color: #000000;\"\u003eAs a class of non-coding RNAs, miRNAs participate in many key events in hepatocarcinogenesis. As a liver-specific and abundant miRNA, miR-192-5p is significantly down-regulated in HCC tumor tissues due to a high frequency of \u003cem\u003eTP53\u003c/em\u003e mutation and its promoter hypermethylation [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e]. In HCC cells, miR-192-5p loss not only led to increased populations of hepatic CSCs [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23\u003c/a\u003e] and cell malignancy phenotype [\u003ca style=\"color: #000000; text-decoration: underline;\" href=\"#_ENREF_23\"\u003e23-25\u003c/a\u003e], but also a hyperglycolytic phenotype as revealed in this study (Fig 7g). These findings suggested a potentially important role of miR-192-5p in regulating hepatic carcinogenesis. It will be essential to establish conditional miR-192-5p knockout mice and further investigate whether and at which stages miR-192-5p regulates hepatic carcinogenesis, the metabolism of HCC tumor and tumor microenvironment, as well as roles of the miR-192-5p-related metabolic pathways in regulating hepatic carcinogenesis in the mouse model. \u0026nbsp;\u003c/span\u003e\u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003eTaken together, loss of miR-192-5p promoted hyperglycolysis and stemness features in CSC\u003csup\u003e+\u003c/sup\u003eHCC cells through upregulating c-Myc and two glycolytic enzymes \u003cu\u003e(\u003c/u\u003e\u003cu\u003ePFKFB3 and GLUT1), and through\u003c/u\u003e a miR-192-5p/c-Myc positive feedback circuit. Moreover, HCC cells with miR-192-5p loss could also actively interact with their environmental non-tumor cells \u003cu\u003eto produce\u003c/u\u003e a more aggressive and stemness-related feature via the over-produced lactic acid by HCC cells and consequently activating a lactate/MCT1/NDRG3/pERK pathway in their environmental non-tumor cells (Fig. 7g).\u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cp\u003eCancer stem cell, CSC; Dulbecco\u0026rsquo;s modified Eagle\u0026rsquo;s medium, DMEM; Extracellular acidification rate, ECAR; Fructose-6-Phosphate, F6P; Gas chromatography/time-of-flight mass spectrometry, GC/TOF-MS; Gene Set Enrichment Analysis, GSEA; Glucose-6-Phosphate, G6P; Hepatocellular carcinoma, HCC; Hepatic stellate cell, HSC; Japanese Collection of Research Biosources Cell Bank, JCRB; Kyoto Encyclopedia of Genes and Genomes, KEGG; Oxygen consumption rate, OCR; Quantitative reverse transcription polymerase chain reaction, qRT-PCR; Single-guide RNAs, sgRNAs; Nicotinamide Adenine Dinucleotide Phosphate, NADPH; The Cancer Genome Atlas, TCGA; 2-deoxyglucose, 2-DG; 2-(N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino)-2-deoxy-D-glucose, 2-NBDG.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003emiRNA transcriptome in HCC Cohort 1, GSE6857; mRNA transcriptome in HCC Cohort 1, GSE14520; metabolomics data in HCC Cohort 1 were referenced in Ref 12. miRNA sequencing and mRNA sequencing data were downloaded from TCGA. TP53 mutation status and MYC amplification information in HCC Cohort 2 were\u003cu\u003e obtained\u003c/u\u003e from www.cbioportal.org/index.do.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors have no conflicts of interest with the material presented in this manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was supported by National Natural Science Foundation of China (No.81672905, 81874054, and 82073055) (J. Ji), National Key R\u0026amp;D Program of China (2018YFA0800504) (J. Ji), the Fundamental Research Funds for the Central Universities in China (J. Ji), and the Thousand Young Talents Plan of China (J. Ji). X.W.W. and A.B. were supported by the intramural program of the Center for Cancer Research, the U.S. National Cancer Institute. S.R. was supported by the German Research Foundation (DFG, Deutsche Forschungsgemeinschaft) within SFB/TR 209 Liver Cancer Project-ID 314905040.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026rsquo; contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eConception and design: Y. Gu., and J. Ji.; Development of methodology: Y. Gu., F. Ji., N. Liu., J. Ji.; Acquisition of data: Y. Gu., N. Liu., Y. Zhao., S. Hu., J. Ji.; Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): Y. Gu., F. Ji., X. Wei., J. Ji.; Writing, review, and/or revision of the manuscript: Y. Gu., X. Wang., A. Budhu., S. Roessler., J. Ji.; Administrative, technical, or material support: X. Zheng., B. Zhao., X. Wang., A. Budhu., W. Jia., J. Ji.; Study supervision: J. Ji.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgments\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe thank Dr. Chaohui Yu for the HL7702 cell line, Dr. Long Zhang for the THP1 cell line, Dr. Xing Guo for PX459 plasmid and Dr. Hai Song for the anti-pERK and anti-ERK antibodies. The authors thank our core facility at Life Sciences Institute for their support on flow cytometry analysis and measurement using Sea Horse XF Analyzers.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eGorrini C, Mak TW: Glutathione Metabolism: An Achilles' Heel of ARID1A-Deficient Tumors. Cancer Cell. 2019;35:161-3.\u003c/li\u003e\n\u003cli\u003eWard PS, Thompson CB: Metabolic reprogramming: a cancer hallmark even warburg did not anticipate. Cancer Cell. 2012;21:297-308.\u003c/li\u003e\n\u003cli\u003eGupta PB, Chaffer CL, Weinberg RA: Cancer stem cells: mirage or reality? Nat Med. 2009;15:1010-2.\u003c/li\u003e\n\u003cli\u003eOikawa T: Cancer Stem cells and their cellular origins in primary liver and biliary tract cancers. Hepatology. 2016;64:645-51.\u003c/li\u003e\n\u003cli\u003ePlaks V, Kong N, Werb Z: The cancer stem cell niche: how essential is the niche in regulating stemness of tumor cells? Cell Stem Cell. 2015;16:225-38.\u003c/li\u003e\n\u003cli\u003eDeshmukh A, Deshpande K, Arfuso F, Newsholme P, Dharmarajan A: Cancer stem cell metabolism: a potential target for cancer therapy. Mol Cancer. 2016;15:69.\u003c/li\u003e\n\u003cli\u003eMaman S, Witz IP: A history of exploring cancer in context. Nat Rev Cancer. 2018;18:359-76.\u003c/li\u003e\n\u003cli\u003eBatlle E, Clevers H: Cancer stem cells revisited. Nat Med. 2017;23:1124-34.\u003c/li\u003e\n\u003cli\u003eVermeulen L, De Sousa EMF, van der Heijden M, Cameron K, de Jong JH, Borovski T\u003cem\u003e et al\u003c/em\u003e: Wnt activity defines colon cancer stem cells and is regulated by the microenvironment. Nat Cell Biol. 2010;12:468-76.\u003c/li\u003e\n\u003cli\u003eAhmed N, Escalona R, Leung D, Chan E, Kannourakis G: Tumour microenvironment and metabolic plasticity in cancer and cancer stem cells: Perspectives on metabolic and immune regulatory signatures in chemoresistant ovarian cancer stem cells. Semin Cancer Biol. 2018;53:265-81.\u003c/li\u003e\n\u003cli\u003ePeiris-Pages M, Martinez-Outschoorn UE, Pestell RG, Sotgia F, Lisanti MP: Cancer stem cell metabolism. Breast Cancer Res. 2016;18:55.\u003c/li\u003e\n\u003cli\u003eVillanueva A: Hepatocellular Carcinoma. N Engl J Med. 2019;380:1450-62.\u003c/li\u003e\n\u003cli\u003eLee TK, Castilho A, Cheung VC, Tang KH, Ma S, Ng IO: CD24(+) liver tumor-initiating cells drive self-renewal and tumor initiation through STAT3-mediated NANOG regulation. Cell Stem Cell. 2011;9:50-63.\u003c/li\u003e\n\u003cli\u003eYamashita T, Ji J, Budhu A, Forgues M, Yang W, Wang HY\u003cem\u003e et al\u003c/em\u003e: EpCAM-positive hepatocellular carcinoma cells are tumor-initiating cells with stem/progenitor cell features. Gastroenterology. 2009;136:1012-24.\u003c/li\u003e\n\u003cli\u003eYang ZF, Ho DW, Ng MN, Lau CK, Yu WC, Ngai P\u003cem\u003e et al\u003c/em\u003e: Significance of CD90+ cancer stem cells in human liver cancer. Cancer Cell. 2008;13:153-66.\u003c/li\u003e\n\u003cli\u003eMa S, Tang KH, Chan YP, Lee TK, Kwan PS, Castilho A\u003cem\u003e et al\u003c/em\u003e: miR-130b Promotes CD133(+) liver tumor-initiating cell growth and self-renewal via tumor protein 53-induced nuclear protein 1. Cell Stem Cell. 2010;7:694-707.\u003c/li\u003e\n\u003cli\u003eMa S, Chan KW, Hu L, Lee TK, Wo JY, Ng IO\u003cem\u003e et al\u003c/em\u003e: Identification and characterization of tumorigenic liver cancer stem/progenitor cells. Gastroenterology. 2007;132:2542-56.\u003c/li\u003e\n\u003cli\u003eGu Y, Zheng X, Ji J: Liver cancer stem cells as a hierarchical society: yes or no? Acta Biochim Biophys Sin (Shanghai). 2020;52:723-35.\u003c/li\u003e\n\u003cli\u003ePrasetyanti PR, Medema JP: Intra-tumor heterogeneity from a cancer stem cell perspective. Mol Cancer. 2017;16:41.\u003c/li\u003e\n\u003cli\u003eSong K, Kwon H, Han C, Zhang J, Dash S, Lim K\u003cem\u003e et al\u003c/em\u003e: Active glycolytic metabolism in CD133(+) hepatocellular cancer stem cells: regulation by MIR-122. Oncotarget. 2015;6:40822-35.\u003c/li\u003e\n\u003cli\u003eChen CL, Uthaya Kumar DB, Punj V, Xu J, Sher L, Tahara SM\u003cem\u003e et al\u003c/em\u003e: NANOG Metabolically Reprograms Tumor-Initiating Stem-like Cells through Tumorigenic Changes in Oxidative Phosphorylation and Fatty Acid Metabolism. Cell Metab. 2016;23:206-19.\u003c/li\u003e\n\u003cli\u003eBudhu A, Roessler S, Zhao X, Yu Z, Forgues M, Ji J\u003cem\u003e et al\u003c/em\u003e: Integrated metabolite and gene expression profiles identify lipid biomarkers associated with progression of hepatocellular carcinoma and patient outcomes. Gastroenterology. 2013;144:1066-75 e1.\u003c/li\u003e\n\u003cli\u003eGu Y, Wei X, Sun Y, Gao H, Zheng X, Wong LL\u003cem\u003e et al\u003c/em\u003e: miR-192-5p Silencing by Genetic Aberrations Is a Key Event in Hepatocellular Carcinomas with Cancer Stem Cell Features. Cancer Res. 2019;79:941-53.\u003c/li\u003e\n\u003cli\u003eLian J, Jing Y, Dong Q, Huan L, Chen D, Bao C\u003cem\u003e et al\u003c/em\u003e: miR-192, a prognostic indicator, targets the SLC39A6/SNAIL pathway to reduce tumor metastasis in human hepatocellular carcinoma. Oncotarget. 2016;7:2672-83.\u003c/li\u003e\n\u003cli\u003eYang YM, Lee WH, Lee CG, An J, Kim ES, Kim SH\u003cem\u003e et al\u003c/em\u003e: Galpha12 gep oncogene deregulation of p53-responsive microRNAs promotes epithelial-mesenchymal transition of hepatocellular carcinoma. Oncogene. 2015;34:2910-21.\u003c/li\u003e\n\u003cli\u003eBudhu A, Jia HL, Forgues M, Liu CG, Goldstein D, Lam A\u003cem\u003e et al\u003c/em\u003e: Identification of metastasis-related microRNAs in hepatocellular carcinoma. Hepatology. 2008;47:897-907.\u003c/li\u003e\n\u003cli\u003eJi J, Shi J, Budhu A, Yu Z, Forgues M, Roessler S\u003cem\u003e et al\u003c/em\u003e: MicroRNA expression, survival, and response to interferon in liver cancer. N Engl J Med. 2009;361:1437-47.\u003c/li\u003e\n\u003cli\u003eJi J, Yamashita T, Budhu A, Forgues M, Jia HL, Li C\u003cem\u003e et al\u003c/em\u003e: Identification of microRNA-181 by genome-wide screening as a critical player in EpCAM-positive hepatic cancer stem cells. Hepatology. 2009;50:472-80.\u003c/li\u003e\n\u003cli\u003eRoessler S, Long EL, Budhu A, Chen Y, Zhao X, Ji J\u003cem\u003e et al\u003c/em\u003e: Integrative genomic identification of genes on 8p associated with hepatocellular carcinoma progression and patient survival. Gastroenterology. 2012;142:957-66 e12.\u003c/li\u003e\n\u003cli\u003eSun Y, Ji F, Kumar MR, Zheng X, Xiao Y, Liu N\u003cem\u003e et al\u003c/em\u003e: Transcriptome integration analysis in hepatocellular carcinoma reveals discordant intronic miRNA-host gene pairs in expression. Int J Biol Sci. 2017;13:1438-49.\u003c/li\u003e\n\u003cli\u003eTao J, Ji J, Li X, Ding N, Wu H, Liu Y\u003cem\u003e et al\u003c/em\u003e: Distinct anti-oncogenic effect of various microRNAs in different mouse models of liver cancer. Oncotarget. 2015;6:6977-88.\u003c/li\u003e\n\u003cli\u003eJi J, Zhao L, Budhu A, Forgues M, Jia HL, Qin LX\u003cem\u003e et al\u003c/em\u003e: Let-7g targets collagen type I alpha2 and inhibits cell migration in hepatocellular carcinoma. J Hepatol. 2010;52:690-7.\u003c/li\u003e\n\u003cli\u003eJi J, Zheng X, Forgues M, Yamashita T, Wauthier EL, Reid LM\u003cem\u003e et al\u003c/em\u003e: Identification of microRNAs specific for epithelial cell adhesion molecule-positive tumor cells in hepatocellular carcinoma. Hepatology. 2015;62:829-40.\u003c/li\u003e\n\u003cli\u003eDimmer KS, Friedrich B, Lang F, Deitmer JW, Broer S: The low-affinity monocarboxylate transporter MCT4 is adapted to the export of lactate in highly glycolytic cells. Biochem J. 2000;350 Pt 1:219-27.\u003c/li\u003e\n\u003cli\u003eDoherty JR, Cleveland JL: Targeting lactate metabolism for cancer therapeutics. J Clin Invest. 2013;123:3685-92.\u003c/li\u003e\n\u003cli\u003eSabnis HS, Somasagara RR, Bunting KD: Targeting MYC Dependence by Metabolic Inhibitors in Cancer. Genes (Basel). 2017;8.\u003c/li\u003e\n\u003cli\u003eStine ZE, Walton ZE, Altman BJ, Hsieh AL, Dang CV: MYC, Metabolism, and Cancer. Cancer Discov. 2015;5:1024-39.\u003c/li\u003e\n\u003cli\u003eLi HM, Yang JG, Liu ZJ, Wang WM, Yu ZL, Ren JG\u003cem\u003e et al\u003c/em\u003e: Blockage of glycolysis by targeting PFKFB3 suppresses tumor growth and metastasis in head and neck squamous cell carcinoma. J Exp Clin Cancer Res. 2017;36:7.\u003c/li\u003e\n\u003cli\u003eShibuya K, Okada M, Suzuki S, Seino M, Seino S, Takeda H\u003cem\u003e et al\u003c/em\u003e: Targeting the facilitative glucose transporter GLUT1 inhibits the self-renewal and tumor-initiating capacity of cancer stem cells. Oncotarget. 2015;6:651-61.\u003c/li\u003e\n\u003cli\u003eZhang HL, Wang MD, Zhou X, Qin CJ, Fu GB, Tang L\u003cem\u003e et al\u003c/em\u003e: Blocking preferential glucose uptake sensitizes liver tumor-initiating cells to glucose restriction and sorafenib treatment. Cancer Lett. 2017;388:1-11.\u003c/li\u003e\n\u003cli\u003eAkita H, Marquardt JU, Durkin ME, Kitade M, Seo D, Conner EA\u003cem\u003e et al\u003c/em\u003e: MYC activates stem-like cell potential in hepatocarcinoma by a p53-dependent mechanism. Cancer Res. 2014;74:5903-13.\u003c/li\u003e\n\u003cli\u003ePichiorri F, Suh SS, Rocci A, De Luca L, Taccioli C, Santhanam R\u003cem\u003e et al\u003c/em\u003e: Downregulation of p53-inducible microRNAs 192, 194, and 215 impairs the p53/MDM2 autoregulatory loop in multiple myeloma development. Cancer Cell. 2010;18:367-81.\u003c/li\u003e\n\u003cli\u003eLee DC, Sohn HA, Park ZY, Oh S, Kang YK, Lee KM\u003cem\u003e et al\u003c/em\u003e: A lactate-induced response to hypoxia. Cell. 2015;161:595-609.\u003c/li\u003e\n\u003cli\u003eBray F, Ferlay J, Soerjomataram I, Siegel RL, 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:394-424.\u003c/li\u003e\n\u003cli\u003eBrodt P: Role of the Microenvironment in Liver Metastasis: From Pre- to Prometastatic Niches. Clin Cancer Res. 2016;22:5971-82.\u003c/li\u003e\n\u003cli\u003eCraig AJ, von Felden J, Garcia-Lezana T, Sarcognato S, Villanueva A: Tumour evolution in hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol. 2020;17:139-52.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"journal-of-experimental-and-clinical-cancer-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"jecc","sideBox":"Learn more about [Journal of Experimental \u0026 Clinical Cancer Research](http://jeccr.biomedcentral.com)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/jecc/default.aspx","title":"Journal of Experimental \u0026 Clinical Cancer Research","twitterHandle":"@OncoBioMed","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Hepatocellular carcinoma, Cancer stem cell, miR-192-5p, Glycolysis, c-Myc","lastPublishedDoi":"10.21203/rs.3.rs-56821/v3","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-56821/v3","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground:\u003c/strong\u003e Emerging studies revealed that cancer stem cells (CSCs) possessed peculiar metabolic properties, which however remained largely unknown in hepatocellular carcinoma (HCC). Genetic silencing of liver-abundant miR-192-5p was a key feature for multiple groups of CSC-positive HCCs. We thus aimed to investigate essential metabolic features of hepatic CSCs via using HCCs with miR-192-5p silencing as a model.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eMethods:\u003c/strong\u003e Datasets from two independent HCC cohorts were used. Data integration analyses of miR-192-5p with metabolome and mRNA transcriptome \u003cu\u003edata\u003c/u\u003e in HCC Cohort 1 were performed to investigate miR-192-5p related metabolic features, which was further validated in Cohort 2. Cellular and molecular assays were performed to examine whether and how miR-192-5p regulated the identified metabolic \u003cu\u003efeatures\u003c/u\u003e. Co-culture systems consisting of HCC cells and LX2 (human hepatic stellate cell line) or THP1 (human monocyte cell line) were established to explore effects of the identified metabolic \u003cu\u003eproperties\u003c/u\u003e on stemness features of HCC cells via interacting with co-cultured non-tumor cells.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eResults: \u003c/strong\u003eHigh levels of glycolysis-related metabolites and genes \u003cu\u003ewere present\u003c/u\u003e in HCCs with low miR-192-5p and CSC-positive HCCs in two independent HCC cohorts. miR-192-5p knockout cells displayed CSC features and miR-192-5p loss led to an enhanced glycolytic phenotype via upregulating three \u003cem\u003ebona fide \u003c/em\u003etargets, GLUT1 and PFKFB3 (two glycolytic enzymes) and c-Myc (regulating glycolytic genes’ expression). Meanwhile, c-Myc suppressed miR-192-5p transcription, ensuring a low-miR-192-5p/high-c-Myc loop to maintain hyperglycolysis. Moreover, over-produced lactic acid from hyperglycolytic HCC cells stimulated the ERK phosphorylation of co-cultured LX2 and THP1 non-tumor cells partially via NDRG3 and MCT1, which in turn promoted cell malignancy and stemness of HCC cells. Consistently, HCC patients with low level of miR-192-5p in their tumor tissues and high level of NDRG3 or MCT1 in their non-tumor tissues had the shortest overall survival.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eConclusions:\u003c/strong\u003e In CSC-positive HCCs, miR-192-5p loss enhanced glycolysis and over produced lactate might further increase HCC malignant features via interacting with environmental non-tumor cells.\u003c/p\u003e","manuscriptTitle":"Loss of MiR-192-5p Initiates a Hyperglycolysis and Stemness Positive Feedback in Hepatocellular Carcinoma.","msid":"","msnumber":"","nonDraftVersions":[{"code":3,"date":"2020-11-24 21:26:40","doi":"10.21203/rs.3.rs-56821/v3","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Accept","date":"2020-11-19T00:00:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2020-11-18T00:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2020-11-17T23:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2020-11-17T23:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"journal-of-experimental-and-clinical-cancer-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"jecc","sideBox":"Learn more about [Journal of Experimental \u0026 Clinical Cancer Research](http://jeccr.biomedcentral.com)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/jecc/default.aspx","title":"Journal of Experimental \u0026 Clinical Cancer Research","twitterHandle":"@OncoBioMed","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}},{"code":2,"date":"2020-11-05 15:47:46","doi":"10.21203/rs.3.rs-56821/v2","editorialEvents":[{"type":"communityComments","content":0},{"type":"editorInvitedReview","content":"","date":"2020-11-05T00:00:00+00:00","index":3,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"decision","content":"Minor revision","date":"2020-11-05T00:00:00+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2020-11-02T00:00:00+00:00","index":1,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"editorInvitedReview","content":"","date":"2020-11-02T00:00:00+00:00","index":2,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"reviewerAgreed","content":"","date":"2020-11-01T01:00:00+00:00","index":3,"fulltext":""},{"type":"reviewerAgreed","content":"","date":"2020-11-01T00:00:00+00:00","index":2,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2020-10-31T00:00:00+00:00","index":"","fulltext":""},{"type":"reviewerAgreed","content":"","date":"2020-10-31T00:00:00+00:00","index":1,"fulltext":""},{"type":"editorAssigned","content":"","date":"2020-10-30T00:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2020-10-29T23:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2020-10-29T23:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"journal-of-experimental-and-clinical-cancer-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"jecc","sideBox":"Learn more about [Journal of Experimental \u0026 Clinical Cancer Research](http://jeccr.biomedcentral.com)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/jecc/default.aspx","title":"Journal of Experimental \u0026 Clinical Cancer Research","twitterHandle":"@OncoBioMed","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}},{"code":1,"date":"2020-08-17 16:52:43","doi":"10.21203/rs.3.rs-56821/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Major revision","date":"2020-08-22T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2020-08-20T12:00:00+00:00","index":3,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"editorInvitedReview","content":"","date":"2020-08-13T12:00:00+00:00","index":1,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"editorInvitedReview","content":"","date":"2020-08-13T12:00:00+00:00","index":2,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"reviewerAgreed","content":"","date":"2020-08-11T12:00:00+00:00","index":3,"fulltext":""},{"type":"editorAssigned","content":"","date":"2020-08-10T12:00:00+00:00","index":"","fulltext":""},{"type":"reviewersInvited","content":"","date":"2020-08-10T12:00:00+00:00","index":"","fulltext":""},{"type":"reviewerAgreed","content":"","date":"2020-08-10T12:00:00+00:00","index":1,"fulltext":""},{"type":"reviewerAgreed","content":"","date":"2020-08-10T12:00:00+00:00","index":2,"fulltext":""},{"type":"submitted","content":"","date":"2020-08-09T12:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2020-08-09T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2020-08-09T12:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"journal-of-experimental-and-clinical-cancer-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"jecc","sideBox":"Learn more about [Journal of Experimental \u0026 Clinical Cancer Research](http://jeccr.biomedcentral.com)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/jecc/default.aspx","title":"Journal of Experimental \u0026 Clinical Cancer Research","twitterHandle":"@OncoBioMed","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"89c66e14-4f82-4e13-a650-8aab5c3b7d84","owner":[],"postedDate":"November 24th, 2020","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":1200897,"name":"Cancer Biology"}],"tags":[],"updatedAt":"2020-12-06T15:03:00+00:00","versionOfRecord":{"articleIdentity":"rs-56821","link":"https://doi.org/10.1186/s13046-020-01785-7","journal":{"identity":"journal-of-experimental-and-clinical-cancer-research","isVorOnly":false,"title":"Journal of Experimental \u0026 Clinical Cancer Research"},"publishedOn":"2020-11-30 15:02:01","publishedOnDateReadable":"November 30th, 2020"},"versionCreatedAt":"2020-11-24 21:26:40","video":"","vorDoi":"10.1186/s13046-020-01785-7","vorDoiUrl":"https://doi.org/10.1186/s13046-020-01785-7","workflowStages":[]},"version":"v3","identity":"rs-56821","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-56821","identity":"rs-56821","version":["v3"]},"buildId":"-HB7Z8yhvgn0wM9Nzuekk","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}