Activation of miR-500a-3p/CDK6 axis suppresses aerobic glycolysis and colorectal cancer progression
Journal of Translational Medicine volume 20, Article number: 106 (2022)
Colorectal cancer (CRC) is one of the lethal cancers with a high mortality rate worldwide and understanding the mechanisms behind its progression is critical for improving patients’ prognosis and developing therapeutics. MiR-500a-3p has been demonstrated to be involved in the progression of several human cancers but its role in CRC remains unclear. The aim of this study is to uncover the expression pattern and mechanisms of action of miR-500a-3p during the CRC progression.
The expression of miR-500a-3p and Cyclin-dependent kinases 6 (CDK6) in 134 CRC tissues were tested by quantitative PCR (qPCR) and immunohistochemistry staining (IHC), respectively. The effect of miR-500a-3p on cell proliferation was explored in vitro and in vivo. The glycolysis of CRC cells was determined by Mass Spectrometry and Seahorse XF 96 Extracellular Flux Analyzer. A dual-luciferase reporter assay was performed to validate the relationship between miR-500a-3p and CDK6.
miR-500a-3p was abnormally downregulated in CRC tissues and cell lines and was negatively associated with a worse prognosis. miR-500a-3p mimics impeded CRC cell proliferation in vitro and in vivo. miR-500a-3p inhibited glucose consumption, lactate and ATP production, and down-regulated the expression of hexokinase2 (HK2). In silico prediction combined with western blot and luciferase assay confirmed that CDK6 is a direct target of miR-500a-3p. Overexpression of CDK6 phenotypically rescued the inhibitory effect of miR-500a-3p on the proliferation and glycolysis of CRC cells.
Our study revealed a potential tumor-suppressive role of miR-500a-3p in CRC, specifically targeting CDK6 and inhibiting cancer cell proliferation and aerobic glycolysis, which may provide new insights into novel prognostic biomarkers and therapeutic targets for CRC.
Colorectal cancer (CRC) is the third most commonly diagnosed cancer in the world with the second-highest mortality rate . The stage status of tumors often determine the prognosis of CRC patients. The 5-year survival rate of CRC patients ranges from greater than 90% in patients with stage I disease to slightly greater than 10% in patients with stage IV disease . Therefore, understanding the mechanisms behind the CRC progression is the key to improving patients’ prognosis.
MicroRNAs (miRNAs) are a class of single-stranded non-coding RNAs with a length of 21–25 nucleotides. miRNAs regulate the expression of approximately 60% of the protein-coding genes . Aerobic glycolysis, one of the central contributors to cancer progression , help cancer cell meet the energy demand, accumulate glycolytic intermediates for cancer biomass synthesis and create an acidic microenvironment . Existing studies have proved that microRNAs play an important role in the progression of CRC by regulating key glycolytic genes directly or indirectly [6, 7]. For example, miR-34a-5p inhibited HK1-mediated glycolysis while the miR-122-PKM2 axis contributed to attenuated glycolysis in CRC [8, 9]. miR-500a is located within the p11 locus on the X chromosome and plays various roles in cancer. miR-500a is upregulated in prostate cancer and promotes cell proliferation and invasion . miR-500a promotes cell stemness via different pathways in gastric and liver cancer [11, 12]. However, miR-500a suppresses cell proliferation, migration and invasion in non-small cell lung cancer . According to the data of a cohort of 551 CRC patients derived from the cBioportal database, miR-500a has the most significant correlation with prognosis (with the lowest P value). miR-500a has two arms, miR-500a-5p and miR-500a-3p. miR-500a-5p is downregulated and acts as a tumor suppresser in CRC. miR-500a-5p may attenuate epithelial‑mesenchymal transition through targeting the transforming growth factor‑β signaling pathway  and suppress cell proliferation by targeting the p300/YY1/HDAC2 axis . However, the role of miR-500a-3p has not been elucidated in CRC.
Cyclin-dependent kinases 6 (CDK6) is known as a classic cell cycle kinase that facilitates the proliferation of cells through the early G1 phase of the cell cycle. CDK6 participates in the process of cancer progression through the kinase-dependent or non-kinase-dependent function . Recently, the role of CDK6 in metabolic regulation has been reported. Wang et al. found that CDK6 regulates the catalytic activity of two key enzymes in the glycolytic pathway, 6-phosphofructokinase and pyruvate kinase M2 . Nevertheless, the role of CDK6 in regulating CRC metabolism remains to be revealed.
In the current study, we investigated the effects of miR-500a-3p on the cell proliferation and glycolysis in CRC using patient samples and cell lines and uncovered its potential link to CDK6.
miR-500a-3p is down regulated in tumor tissue and its low expression predicts poor prognosis
According to the data of a cohort of 551 patients from the cBioportal database, we found that the expression of has-mir-500a was significantly associated with the prognosis of CRC patients (Additional file 1: Table S1) Patients who had low miR-500a expression had significantly poor overall survival (P = 0.008, Additional file 2: Fig. S1A) and progression-free survival (P < 0.001, Additional file 2: Fig. S1B). We firstly examined the miR-500a-3p expression by qPCR in CRC cell lines, including HCT116, HT29, RKO, SW620, SW48 and SW480, and normal human colon epithelial cell line NCM460. The expression level of miR-500a-3p was significantly down-regulated in CRC cells than that in NCM460 cells (Fig. 1A). In addition, it was significantly downregulated in CRC tissues compared to adjacent non-tumor tissues (Fig. 1B). Next, we examined miR-500a-3p expression in the cohort that contained 134 continuous CRC specimens. We divided the patients into low or high-expression groups by a cut-off of the median expression value. We found that the expression level of miR-500a-3p was negatively correlated with tumor size (P = 0.022) and differentiation (P = 0.043) (Table 1). Survival analysis showed that low miR-500a-3p expression was significantly correlated with poor overall survival (P = 0.0026, Fig. 1C) and disease-free survival (P = 0.0276, Fig. 1D). To further evaluate the prognostic value of miR-500a-3p, both univariate and multivariate analyses were performed. The results indicated that high miR-500a-3p expression is an independent protective factor for CRC [HR, 0.109; 95% confidence interval (CI), 0.014–0.875; P = 0.037]. (Additional file 1: Table S2). We also analyzed the correlation of miR-500a-3p expression with the biomarkers of cell proliferation (Ki-67) and glucose uptake (Standard Uptake Value, SUV) in CRC tissues. The expression of miR-500a-3p was significantly lower in the high Ki67 group than that in the low Ki67 group (Fig. 1E). Besides, in patients receiving PETCT preoperatively (n = 52), there was a negative correlation between the expression of miR-500a-3p and SUV with the P value of 0.038 and r value of -0.262 (Fig. 1F). Altogether, low expression level of miR-500a-3p predicts poor prognosis in CRC patients and it may suppress CRC by regulating tumor proliferation and glucose metabolism.
miR-500a-3p suppresses CRC growth in vitro and in vivo
HCT116/SW480 cells were transfected with miR-500a-3p mimics to investigate the function of miR-500a-3p in vitro. Transfection efficiencies were confirmed by RT-qPCR (Additional file 2: Fig. S2A). CCK-8 and colony formation assays indicated that CRC cell proliferation and colony formation were significantly inhibited by miR-500a-3p mimics (Fig. 2A–D). Next, fluorescence-activated cell sorting (FACS) analysis revealed an increase in G1-phase cells and a concomitant decrease in S-phase cells by miR-500a-3p (Fig. 2E, F). We also found that miR-500a-3p mimics did not promote apoptosis significantly (Additional file 2: Fig. S2B). Transwell assay showed that the miR-500a-3p mimics markedly inhibited CRC cell invasion (Additional file 2: Fig. S2C). We established xenograft tumor models with SW480 and HCT116 cells stably overexpressing miR-500a-3p (Additional file 2: Fig. S2D, E). The mean volume and weight of tumors were lower in the miR-500a-3p-overexpressing group than the vector group (Fig. 2G–I and Additional file 2: Fig. S2F–H). The protein expression of cell proliferation (Ki-67) in the xenograft tumors was lower in miR-500a-3p overexpression group than vector group (Fig. 2J). Altogether, these results indicated that miR-500a-3p suppresses CRC cell proliferation in vitro and in vivo.
miR-500a-3p suppresses aerobic glycolysis in CRC
The expression of miR-500a-3p in CRC patients is negatively correlated with glucose uptake, indicating that miR-500a-3p may have an inhibitory effect on glycolysis. Cellular glycolytic activity was measured. As a result, overexpression of miR-500a-3p significantly decreased glucose uptake (Fig. 3A), lactate production (Fig. 3B) and intracellular ATP levels (Fig. 3C) in both HCT116 and SW480 cells. The extracellular acidification rate (ECAR), measured by the Seahorse analyzer, was downregulated by miR-500a-3p (Fig. 3D). Furthermore, targeted metabolomics focused on the energy metabolism pathway was performed, and lower level of glycolytic metabolites like glucose 6-phosphate, fructose 6-phosphate, fructose 1,6-bisphosphate, dihydroxyacetone phosphate and 3-phospho-glycerate caused by overexpression of miR-500a-3p was found (Fig. 3E and Additional file 2: Fig. S3A–E). We then measured the expression of major genes involved in glucose transport and glycolysis. miR-500a-3p mimics reduced the RNA level of HK2, PFKP, PGK1 and LDHA (Fig. 3F) and the protein level of GLUT1, HK2 and PKM2 (Fig. 3G). Since both RNA and protein levels of HK2 were downregulated by miR-500a-3p, HK2 was evaluated in the CRC cohort. It was negatively correlated with the miR-500a-3p levels (r = − 0.418, P < 0.0001, Fig. 3H, I), confirming the negative regulation of HK2 by miR-500a-3p in CRC tissues. These data suggest that the overexpression of miR-500a-3p causes a decrease in glucose metabolism via the regulation of multiple genes involved in glycolysis.
CDK6 is a direct target and down regulated by miR-500a-3p
To explore the mechanisms by which miR-500a-3p inhibits the proliferation and glycolysis of CRC, we firstly used the publicly available algorithms TargetScan and miRDB to find the potential binding sites of glycolytic enzymes, like HK2, PFKP, PGK1 and LDHA, that are downregulated by miR-500a-3p. However, their 3′UTR did not have the conserved binding site of miR-500a-3p, indicating it may regulate glycolysis indirectly (Additional file 2: Fig. S4A). CDK6, a novel glycolysis regulator , was predicted as a target of miR-500a-3p by TargetScan and miRDB database (Fig. 4A). qPCR and Western blotting analysis showed that miR-500a-3p reduced both mRNA and the protein level of CDK6, indicating that CDK6 may act as a downstream target of miR-500a-3p (Fig. 4B, C). To further confirm this, we constructed two pairs of luciferase reporter plasmids containing the predicted seed sequence of miR-500a-3p in the 3′-UTR of the CDK6 mRNA and a control reporter containing the mutated sequence of the same fragment (Fig. 4D). miR-500a-3p significantly inhibited the luciferase activity in cells containing the wild-type CDK6 sequence, but not the mutant sequence (Fig. 4E, F). Finally, we evaluated the protein level of CDK6 in clinical CRC samples using IHC (Fig. 4G), and a significant negative correlation was revealed between miR-500a-3p and CDK6 expression (P = 0.0012, r = − 0.285, Fig. 4H). Collectively, these findings suggest that CDK6 is a direct target of miR-500a-3p in CRC.
CDK6 is upregulated and promotes aerobic glycolysis in CRC
We examined CDK6 protein levels in 20 paired human CRC tumors and adjacent nontumor tissues using IHC (Fig. 5A). The IHC staining results revealed that it was significantly upregulated in CRC tissues compared to nontumor tissues (Fig. 5B). Survival analysis of a large CRC cohort showed that high CDK6 expression was significantly correlated with poor overall survival (P = 0.0333, Additional file 2: Fig. S5A) and disease-free survival (P = 0.0559, Additional file 2: Fig. S5B). CCK-8 and cell cycle assays indicated that CDK6 enhanced CRC cell proliferation (Additional file 2: Fig. S5C, D). Then, we explored the influence of CDK6 on the aerobic glycolysis of CRC cells. As exhibited in Fig. 5C–F, CDK6 promoted glucose uptake, the production of lactate and ATP and the ECAR in SW480 cells. The expression of major genes involved in glucose transport and glycolysis were also measured. The RNA level of GLUT1, LDHB, HK1 and HK2 and the protein level of GLUT1, HK2 were elevated after overexpressing CDK6 (Fig. 5G, H). In CRC tissues, the IHC scores of HK2 was positively correlated with that of CDK6 (r = 0.411, P < 0.0001, Fig. 5I). These findings indicated that CDK6 plays an important role in promoting aerobic glycolysis in CRC.
miR-500a-3p mediates the inhibition of proliferation and aerobic glycolysis in CRC cells by targeting CDK6
To address whether miR-500a-3p exerted its role through targeting CDK6, we conducted rescue experiments. The CCK-8 and colony formation assays showed that CDK6 promoted CRC cell proliferation and colony formation and offset the inhibitory effect caused by miR-500a-3p (Fig. 6A, B). As shown in Fig. 6C, CDK6 overexpression counteracted the inhibitory effect of miR-500a-3p on the cell cycle of CRC cells. We explored the influence of miR-500a-3p/CDK6 axis on the aerobic glycolysis of CRC cells by measuring glucose uptake, lactate production, intracellular ATP level and ECAR. The inhibition of glycolysis caused by miR-500a-3p could be partly eliminated by CDK6 overexpression (Fig. 6D–G). Together, these results indicated that miR-500a-3p suppressed CRC proliferation and aerobic glycolysis through down-regulating CDK6.
Materials and methods
From January 2015 to March 2015, 134 continuous CRC specimens and 20 matching adjacent non-tumor specimens obtained from CRC patients at the Zhongshan Hospital, Fudan University were collected. The inclusion criteria were as follows: with histologically confirmed CRC, receiving radical surgery, having complete follow-up. The exclusion criteria were receiving preoperative chemotherapy or radiation, pregnancy or breastfeeding, previous malignancy within 5 years or no written informed consent. The primary tumor ki67 score is the routine pathological evaluation of our center and was evaluated by two experienced pathologists. Patients in the cohort (134 continuous CRC specimens) who received PETCT preoperatively (n = 52) were included to analyze the association between miR-500a-3p expression and the SUV of the primary tumors. The study was authorized by the Ethics Committee of the Zhongshan Hospital, Fudan University. Written informed consent was obtained from each patient.
Collection of TCGA public data and identification of significantly differential expression miRNA
We collected gene expression profiles of 551 COADREAD samples from the TCGA cohort (http://gdac.broadinstitute.org/). Data including mRNA expression level and miRNA expression level. All the data have been standardized. For miRNAs differential analysis, we calculated the mean expression level of these miRNAs and divided them into high and low expression groups.
The R package DESeq2 (version 1.22.2) was used to calculate the differential expressed t statistics for microarray and RNA sequencing data. We used the univariate Cox proportional hazards model to examine the associations between gene expression and overall survival. MiRNAs with P values less than 0.05 were considered to be statistically significant and included in consensus survival analysis.
The whole cohort which contained 134 continuous CRC specimens were stained with the CDK6 antibody (ab124821, Abcam) and HK2 antibody (2867T, Cell Signaling Technology). The 20 matched adjacent non-tumor specimens were stained with CDK6 antibody (ab124821, Abcam). IHC scores were conducted according to the ratio and intensity of positive-staining areas. The staining areas were scored as follow: 0–25%, score 1; 25–50%, score 2; 50–75%, score 3; and beyond 75%, score 4. The signal intensity was scored on a scale of 0–3: 0-negative, 1-weak, 2-moderate and 3-strong. IHC scores were averaged from two experienced pathologists who scored the IHC staining independently.
Cell lines and culture
Human CRC cell lines [HCT116 (RRID: CVCL_0291) and SW480 (RRID: CVCL_0546)] and normal human colon epithelial cell line NCM460 (RRID: CVCL_0460) were obtained from the National Collection of Authenticated Cell Cultures (Shanghai, China). All human cell lines have been authenticated using STR profiling within the last three years and all experiments were performed with mycoplasma-free cells. These cell lines were maintained in Dulbecco’s Modified Eagle Medium (Logan Utah, HyClone, USA) with 10% fetal bovine serum (FBS; Gibco), 1% penicillin (10 U/mL) and 1% streptomycin (10 μg/mL) in an incubator with the environment of 37 °C and 5% CO2.
Quantitative real-time polymerase chain reaction (qRT-PCR)
RNA was isolated with Trizol reagent. The purity and the concentration of different RNA samples were measured using NanoDrop ND-1000 (Thermo Scientific, USA). Reverse transcription reaction was performed using miRcute miRNA cDNA synthesis Kit (kr211, Tiangen, China) for miRNAs and SuperScript IV VILO cDNA synthesis kit (11756050, Thermo Scientific, USA) for mRNAs. Quantitative RT-PCR (qRT-PCR) was performed using SYBR® Premix Ex Taq™ (Takara) containing mRQ 3′ Primer on an ABI 7500 platform (Applied Biosystems, Carlsbad, CA, USA). The relative quantities of miR-500a-3p in cells were normalized to U6. Beta-actin was used as an internal control for mRNA detection. The primer sequences for miRNA and mRNA detection are listed in Additional file 1: Table S3.
Lentivirus and miRNA transfection
The lentivirus particles overexpressing human miR-500a-3p and CDK6 were purchased from Genomeditech (Shanghai, China). miR-500a-3p mimics were synthesized by Genomeditech (Shanghai, China). The sequences of miR-500a-3p mimic were 5′-AUGCACCUGGGCAAGGAUUCUG-3′. The miRNA oligonucleotides were transfected using Lipofectamine® RNAiMAX™ (50 nmol/L; Invitrogen, CA, USA).
Cell viability assay
At 72 h after lentivirus infection or at 24 h after miRNA transfection, HCT-116 (2500/well) and SW480 (2000/well) cells were plated into 96-well culture plates. Cell viability was measured at 24, 48, 72 and 96 h post-seeding, using the Cell Counting Kit-8 (Dojindo, Kumamoto, Japan) according to the manufacturer’s instructions.
HCT-116 (1500/well) and SW480 (2000/well) cells were seeded in six-well plates after 72 h lentivirus infection. The culture medium was changed every 3 days. After incubating for 2 weeks, colonies were fixed with 4% paraformaldehyde for 15 min and stained with 5% Giemsa for 20 min. The colonies containing at least 50 cells were scored.
Cell cycle analysis
Seventy-two hours after infection, cells were harvested, washed with phosphate-buffered saline (PBS), and fixed in 70% ethanol at 4 °C overnight. After fixation, cells were washed with PBS before suspension in RNase A/propidium iodide solutions (100 mg/mL RNase A and 5 μg/mL propidium iodide). Cells were incubated at room temperature for an hour. Stained cells were analyzed by a FACScan flow cytometer (BD Biosciences, Mountain View, CA, USA).
Cell apoptosis analysis
Cell apoptosis was assessed by annexin V/propidium iodide (BD Biosciences, San Jose, CA, USA). Cells were harvested, washed with PBS, and resuspended in 1× binding buffer at a concentration of 1 × 106 cells/mL. A total of 100 μL solution was transferred to a new tube and added with 5 μL of APC Annexin V and 5 μL of propidium iodide. Cells were incubated at room temperature for 15 min in the dark and then analyzed by a FACScan flow cytometer.
Cell invasion assay
The cell invasion assay was performed in 24-well transwell chambers pre-coated with Matrigel (Corning, NY, USA). HCT-116 cells (105/well). Serum-free medium (200 μL) was placed into the upper chamber, and 600 μL complete medium (with 10% FBS) was filled into the lower chamber. Cells on the inner membrane were removed with a cotton swab after 24 h. The outer membrane was fixed in 4% formaldehyde (in PBS) and stained with 0.5% crystal violet. Cell numbers were counted and averaged in five random fields at a magnification of 100×.
Xenograft tumor assay
BALB/c nude mice (6 weeks old) were purchased from Shanghai SLAC Laboratory animal co. ltd (shanghai, China) and randomly divided into 2 groups (n = 7/6). SW480 (3 × 106) or HCT116 (4 × 106) cells that stably transfected with miR-NC or miR-500a-3p were injected at the back region of BALB/c mice subcutaneously. The volume of tumors was measured twice a week using a vernier caliper with the method of volume = (length × width2)/2. The mice were killed after inoculation for 30 d, and the tumors were weighed. Tumor tissues were subjected to the expression analysis of ki67. The procedures in this study were permitted by the Animal Research Committee of the Zhongshan Hospital, Fudan University.
Luciferase reporter assay
The putative miR-500a-3p binding sequence and the matching mutant binding sequence in the 3′untranslated region (3′UTR) of CDK6 were amplified and inserted into pGL3 luciferase reporter vector (Genomeditech, shanghai, China). SW480 cells were transfected with miR-NC or miR-500a-3p and the above-constructed reporter plasmids. After transfection for 48 h, the luciferase activities in different groups were determined by the dual-luciferase reporter assay system (Promega).
Glucose, lactate and ATP measurement
SW480 and HCT116 cells were cultured with an FBS-free medium and culture medium was collected after 24 h. Glucose Uptake Colorimetric Assay kit (Biovision, Milpitas, California, USA), Lactate Assay Kit (Biovision) and ATP Colorimetric Assay kit (Biovision) were utilized to detect the glucose consumption, the lactate production and ATP in CRC cells according to the manufacturer’s instructions. The glucose, lactate and ATP levels were normalized to total cell protein.
Extracellular acidification rate (ECAR) was detected using a Seahorse XF96 analyzer (Seahorse Biosciences, USA). SW480 cells (105/well) were seeded in a 96-well XF96 microplate (Seahorse Biosciences, USA). Before experiments, cell culture medium was replaced and cells were then incubated with assay medium for 1 h at 37 °C in a CO2-free incubator. ECAR was detected using a sequential injection of 10 mM glucose, 2 mM oligomycin (Sigma-Aldrich) and 50 mM 2-deoxyglucose (2-DG, Sigma-Aldrich). Each cycle of measurement involved mixing (3 min), waiting (2 min), and measuring (3 min) cycles.
The targeted metabolomics analyses were performed using an HPLC system (Agilent 1290, Agilent Technologies) and a mass spectrometer (Agilent 5500, Agilent Technologies). A 10-cm dish of cultured tumor cells was collected, adding 1 mL acetonitrile/methanol water (v, 2:2:1) and storing at − 80℃ after quick freezing in liquid nitrogen. Sample preparation processes were performed in accordance with the above method of parallel preparation of QC samples. MRM transitions representing the metabolites were simultaneously monitored, and the positive/negative polarity switching was used. Data analyses were performed as instructed by Shanghai Applied Protein Technology .
Whole protein extracts were lysed by radioimmunoprecipitation assay (RIPA) buffer (Thermo Fisher Scientific) according to the manufacturer’s protocol. At that time, 30 μg proteins were run on a 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) gel at 100 V for 2 h and transferred to a polyvinylidene fluoride membrane at 80 V for 2 h. Membranes were blocked with 5% bovine serum albumin (BSA) in Tris-buffered saline with Tween-20 (TBST) at room temperature for 1 h. They were then incubated overnight at 4 °C with one of the following primary antibodies: anti-CDK6 (ab124821, Abcam) (1:1000), anti-GLUT1 (ab115730, Abcam) (1:1000) and anti-HK2 (2867T, Cell Signaling Technology) (1:1000) and anti-beta actin (1:1000) from Santa Cruz Biotechnology (Dallas, Texas, USA). After washing with TBST, membranes were incubated with secondary antibodies for 1 h at room temperature and signals were developed using an enhanced chemiluminescence kit (Pierce, Waltham, MA, USA).
SPSS (version 22.0, SPSS Inc.) or GraphPad Prism software (version 7.0, USA) were used to analyze the data and generate the graphs. The differences between the miR-500a-3p expression and the clinical characteristics of CRC patients were analyzed using χ2 test. Survival curves were generated using the Kaplan–Meier method and log-rank tests. Univariate and multivariate Cox regression analyses were conducted to identify the independent factors. The correlation among the expression of miR-500a-3p, CDK6 and HK2 was analyzed using Spearman’s correlation test. Student’s t-test or the Mann–Whitney U test was used to calculate the P values between two groups. A value of P < 0.05 was identified as statistically significant. R software 3.5.1 was applied to this study for the statistical analyses. All miRNA expression data were normalized.
Both investigating new prognostic biomarkers and uncovering the mechanisms behind the progression of CRC are critical for developing CRC therapies. In the current study, we focused on the role of miR-500a-3p in CRC progression. We found that patients with low miR-500a-3p expression had a worse prognosis and miR-500a-3p may represent an independent prognostic factor for CRC patients. Further in vitro and in vivo experiments showed that miR-500a-3p inhibits CRC cell proliferation and aerobic glycolysis. In addition, we discovered CDK6 as a novel downstream target of miR-500a-3p. These data suggested that miR-500a-3p might serve as a novel prognostic biomarker or therapeutic target of CRC.
Accumulating evidence indicates that miRNAs contribute to CRC progression [19,20,21]. However, research studies screening miRNAs from the sequencing results of large cohorts were relatively rare. We analyzed miRNAs that are significantly related to the prognosis of CRC from the cBioportal database and found that miR-500 has the most significant correlation with the prognosis of CRC patients. Then we focused on miR-500a-3p which has not been reported in CRC.
Interestingly, both tumor-promoting and—suppressive roles of miR-500a-3p have been revealed. In both hepatocellular and gastric carcinoma, miR-500a-3p served as a tumor promoter via regulating the cancer cell stemness [11, 12]. While in non-small cell lung cancer, miR-500a-3p acts as a tumor suppressor via downregulating LY6K expression . The distinct roles of miR-500a-3p might be due to the heterogeneity of tumors and differences of cancer types. In this study, we provide new evidence that miR-500a-3p serves as a tumor suppressor in CRC and more importantly, for the first time, reveal the role of miR-500a-3p on cancer cell metabolism. Specifically, we found that miR-500a-3p downregulated the expression of glycolytic enzymes and inhibited glycolysis. Decreased glycolysis restrains the cellular buildings for rapid cell proliferation , and ultimately contributes to the tumor-suppressive role of miR-500a-3p in CRC.
Previous studies have shown that CDK6 is regulated by several miRNAs, such as miR-29b, miR-211 and miR-497 [23,24,25]. Our study provided additional evidence to support that CDK6 is regulated by miRNAs. We report that CDK6 is a direct functional target of miR-500a-3p in CRC and found that miR-500a-3p inhibits the transcription of CDK6. In different cancer types, CDK6 may regulate different metabolic enzymes and play distinct roles in glycolysis. Wang et al. found that CDK6 inhibits the glycolytic pathway and re-directs the glycolytic intermediates into the pentose phosphate pathway (PPP) and serine pathways . Xing et al. proved that CDK6 promotes glycolysis via phosphorylation of the fructose bisphosphate PFK2 (PFKFB3) in breast cancer . In this study, we found that CDK6 enhanced glycolysis in CRC and HK2 might be a potential downstream target of CDK6.
The inhibition of glycolysis by miR-500a-3p was only partially restored by CDK6 overexpression, indicating other targets might also be involved in miR-500a-3p-mediated regulation of glycolysis. Some reported targets of miR-500a-3p, such as XBP1 and FBXW7, may regulate glycolysis directly or indirectly [11, 27,28,29]. Further studies may be needed to validate the mechanistic link of miR-500a-3p and these targets in modulating glycolysis in CRC.
Taken together, miR-500a-3p was found to be a tumor suppressor in CRC. miR-500a-3p suppressed the proliferation, cell cycle and glycolysis of CRC cells partially via CDK6. These findings will expand our understanding of the versatile role of miR-500a-3p in CRC progression and support the rationale of miR-500a-3p in CRC treatment.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its Additional files.
Quantitative real-time PCR
Cell counting kit 8
Cyclin-dependent kinases 6
Standard Uptake Value
Fluorescence-activated cell sorting
Extracellular acidification rate
Pentose phosphate pathway
Yin Yang 1
Histone deacetylase 2
Platelet isoform of phosphofructokinase
Phosphoglycerate kinase 1
Glucose transporter 1
Pyruvate kinase M2
Positron emission tomography computed tomography
The Cancer Genome Atlas
Fetal bovine serum
Sodium dodecyl sulfate poly-acrylamide gel electrophoresis
Bovine serum albumin
Tris buffered saline with Tween-20
Lymphocyte antigen 6 complex locus K
X-box binding protein 1
F-box and WD repeat domain containing 7
Sung H, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.
Brenner H, Kloor M, Pox CP. Colorectal cancer. Lancet. 2014;383(9927):1490–502.
He L, Hannon GJ. MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004;5(7):522–31.
Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009;324(5930):1029–33.
Vaupel P, Schmidberger H, Mayer A. The Warburg effect: essential part of metabolic reprogramming and central contributor to cancer progression. Int J Radiat Biol. 2019;95(7):912–9.
Zhang LF, Jiang S, Liu MF. MicroRNA regulation and analytical methods in cancer cell metabolism. Cell Mol Life Sci. 2017;74(16):2929–41.
Mirzaei H, Hamblin MR. Regulation of glycolysis by non-coding RNAs in cancer: switching on the Warburg effect. Mol Ther Oncolytics. 2020;19:218–39.
Li S, et al. lncARSR sponges miR-34a-5p to promote colorectal cancer invasion and metastasis via hexokinase-1-mediated glycolysis. Cancer Sci. 2020;111(10):3938–52.
Wang X, et al. Exosome-delivered circRNA promotes glycolysis to induce chemoresistance through the miR-122-PKM2 axis in colorectal cancer. Mol Oncol. 2020;14(3):539–55.
Cai B, et al. Inhibition of microRNA-500 has anti-cancer effect through its conditional downstream target of TFPI in human prostate cancer. Prostate. 2017;77(10):1057–65.
Lin H, et al. Exosomal MiR-500a-3p promotes cisplatin resistance and stemness via negatively regulating FBXW7 in gastric cancer. J Cell Mol Med. 2020;24(16):8930–41.
Jiang C, et al. miR-500a-3p promotes cancer stem cells properties via STAT3 pathway in human hepatocellular carcinoma. J Exp Clin Cancer Res. 2017;36(1):99.
Jiang M, et al. Down-regulation of miR-500 and miR-628 suppress non-small cell lung cancer proliferation, migration and invasion by targeting ING1. Biomed Pharmacother. 2018;108:1628–39.
Tang W, et al. MicroRNA-500a-5p inhibits colorectal cancer cell invasion and epithelial–mesenchymal transition. Int J Oncol. 2020;56(6):1499–508.
Tang W, et al. The p300/YY1/miR-500a-5p/HDAC2 signalling axis regulates cell proliferation in human colorectal cancer. Nat Commun. 2019;10(1):663.
Nebenfuehr S, Kollmann K, Sexl V. The role of CDK6 in cancer. Int J Cancer. 2020;147(11):2988–95.
Wang H, et al. The metabolic function of cyclin D3-CDK6 kinase in cancer cell survival. Nature. 2017;546(7658):426–30.
Cai Y, et al. An integrated targeted metabolomic platform for high-throughput metabolite profiling and automated data processing. Metabolomics. 2015;11(6):1575–86.
Zhuang M, et al. MALAT1 sponges miR-106b-5p to promote the invasion and metastasis of colorectal cancer via SLAIN2 enhanced microtubules mobility. EBioMedicine. 2019;41:286–98.
Sun M, et al. Integrated analysis identifies microRNA-195 as a suppressor of Hippo-YAP pathway in colorectal cancer. J Hematol Oncol. 2017;10(1):79.
Cao L, et al. MiR-760 suppresses human colorectal cancer growth by targeting BATF3/AP-1/cyclinD1 signaling. J Exp Clin Cancer Res. 2018;37(1):83.
Liao XH, Xie Z, Guan CN. MiRNA-500a-3p inhibits cell proliferation and invasion by targeting lymphocyte antigen 6 complex locus K (LY6K) in human non-small cell lung cancer. Neoplasma. 2018;65(5):673–82.
Hoareau-Aveilla C, et al. miR-497 suppresses cycle progression through an axis involving CDK6 in ALK-positive cells. Haematologica. 2019;104(2):347–59.
Xia B, et al. miR-211 suppresses epithelial ovarian cancer proliferation and cell-cycle progression by targeting Cyclin D1 and CDK6. Mol Cancer. 2015;14:57.
Zhu K, et al. MiR-29b suppresses the proliferation and migration of osteosarcoma cells by targeting CDK6. Protein Cell. 2016;7(6):434–44.
Xing Z, et al. Expression of long noncoding RNA YIYA promotes glycolysis in breast cancer. Cancer Res. 2018;78(16):4524–32.
Zhou Y, et al. lncRNA NEAT1 regulates gastric carcinoma cell proliferation, invasion and apoptosis via the miR-500a-3p/XBP-1 axis. Mol Med Rep. 2021;24(1):1–10.
Wang C, et al. Fulvestrant inhibits the glycolysis of prolactinoma GH3 cells by downregulating IRE1/XBP1 signaling pathway. Eur Rev Med Pharmacol Sci. 2018;22(16):5364–70.
Zhan P, et al. FBXW7 negatively regulates ENO1 expression and function in colorectal cancer. Lab Invest. 2015;95(9):995–1004.
Supported by The National Natural Science Foundation of China (82072653, 82072678, 81602035); Shanghai Science and Technology Committee Project (19511121300); Clinical Research Plan of SHDC (No. SHDC2020CR1033B, No. SHDC2020CR5006, No. SHDC2020CR3037B) Shanghai Engineering Research Center of Colorectal Cancer Minimally Invasive (17DZ2252600); Youth fund of Zhongshan Hospital (2019ZSQN28).
Ethics approval and consent to participate
The study was approved by the ethics committee of Fudan University Zhongshan Hospital and was performed according to Good Clinical Practice guidelines and the Declaration of Helsinki. Written informed consent was obtained from every patient.
Consent for publication
Consent for publication was obtained from every patient.
No conflict of interest exits in the submission of this manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
MiRNAs that significantly affect the prognosis of CRC patients in cBioportal database. Table S2. Univariate and multivariate analysis showed miR-500a-3p expression was an independent prognostic factor in CRC paients. Table S3. Primer sequences for miRNA and mRNA detection.
(A, B) Low miR-500 expression was significantly associated with poor overall survival (A) and progression-free survival (B) of CRC patients in cBioportal database. Figure S2. (A) qRT-PCR was carried out to detect the level of miR-500a-3p in SW480 and HCT116 cells transfected with miR-NC or miR-500a-3p. (B) miR-500a-3p mimics did not increase apoptosis in HCT116 cells. (C) miR-500a-3p mimics inhibited cell migration in HCT116 cells. (D, E) Verification of miR-500a-3p overexpression efficiency in SW480 and HCT116 cells. (F) Representative tumor diagrams in different groups were shown. (G, H) Tumor volume and weight in the miR-500a-3p group were significantly lower than those in the miR-NC group. (I) Ki-67 expression was significantly higher in tumors of miR-500a-3p group than that of miR-NC group. Figure S3. (A–E) Glucose 6-phosphate (A), fructose 6-phosphate (B), fructose 1,6-bisphosphate (C), dihydroxyacetone phosphate (D) and 3-phospho-glycerate (E) in the miR-500a-3p group were significantly lower than those in the miR-NC group. Figure S4. The 3′UTR of glycolysis enzymes did not have conserved binding site of miR-500a-3p. Figure S5. (A, B) CDK6 high expression was significantly associated with poor overall survival (A) and disease-free survival (B) in CRC specimens. (C, D) CDK6 enhanced CRC cell proliferation using CCK-8 (C) and cell cycle assays (D). Figure S6. (A) The correlation between miR-500a-3p and CDK6 in CPTAC_COAD database (n = 105). (B) The correlation between miR-500a-3p and HK2 in CPTAC_COAD database (n = 105). (C) The representative images of RNA electrophoresis.
About this article
Cite this article
Liu, Y., Tang, W., Ren, L. et al. Activation of miR-500a-3p/CDK6 axis suppresses aerobic glycolysis and colorectal cancer progression. J Transl Med 20, 106 (2022). https://doi.org/10.1186/s12967-022-03308-8
- Colorectal cancer