- Primary research
- Open Access
rs61991156 in miR-379 is associated with low capability of glycolysis of gastric cancer by enhanced regulation of PKM2
© The Author(s) 2018
Received: 11 April 2018
Accepted: 26 June 2018
Published: 4 July 2018
Glycolysis is an important metabolic oncogenic change also play a pivot role in the Warburg effect. Glycolysis related gene PKM2 expressed differently individually. Presently, we sought to investigate the effect of single nucleotide polymorphism (SNP) at rs61991156 of miR-379 on gastric cancer (GC) proliferation and metabolism.
The genotype of rs61991156 in miR-379 was investigated by using real-time PCR. The glycolysis-related metabolites were determined by using GC–TOF–MS. The biological effects of rs61991156 in miR-379 was explored by in vitro studies.
In this study, we found that rs61991156 in miR-379 was involved in the occurrence of GC by acting on the 3′UTR region of PKM2. The clinical data analysis revealed that A > G in rs187960998 was significantly associated with better differentiation, small tumor size, and non-metastasis. In vitro study further revealed that A > G SNP of miR-379 could decrease GC cell proliferation as well as the promoter activity and expression of PKM2. The glycolysis of the patients with miR-379 GG genotype was significantly lower than AG and AA genotype by metabolomics analysis. The patients with AA genotype have significantly lower PKM2 expression compared to the G carrier, while there is no significant expression difference in miR-379 expression. Patients with AA genotype have significantly shorter survival rate compared to the G carrier.
rs61991156 in miR-379 was highly associated with a decreased risk, well differentiation and better post-surgery survival in Chinese population by inhibiting the expression of PKM2.
The Warburg effect is observed in most cancer cells, reflected as predominantly producing energy at a high rate of glycolysis followed by lactic acid fermentation. However, most normal cells demonstrated a comparatively low rate of glycolysis followed by oxidation of pyruvate in mitochondria. Malignant, rapidly growing tumor cells typically have glycolytic rates up to 200 times higher than those of their normal tissues of origin; this occurs even if oxygen is plentiful [1–3]. In the other point of view, the glycolysis is positively correlated to the tumor growth and malignancy. Today, mutations in oncogenes and tumor suppressor genes are thought to be responsible for malignant transformation, and the Warburg effect is considered to be a result of these mutations rather than a cause .
Since the glycolysis provides most of the building blocks required for cell proliferation, cancer cells (and normal proliferating cells) have been proposed to need to activate glycolysis, despite the presence of oxygen, to proliferate . A set of genes regulated the process of glycolysis, three of them are most important. They are hexokinase (HK), phosphofructokinase-1 (PFK1) and pyruvate kinase (PK), in which PKM2 is the last step within glycolysis and was detected to expressing in most cancers .
Similar to other genes, PKM2 was regulated by miRNA. It has been reported that PKM2 can be targeted by the tumor-suppressive miRNA including miR-326 , miR-122 [7, 8], miR-124 , miR-137 , etc., thus to decrease the Warburg effects. However, although the expression of these miRNAs was detectable within tumors, over-expression of their targeting gene: PKM2 still appeared frequently. There are many explanations, in which SNP in miRNA, especially mature miRNA are one of the convincing reasons.
MicroRNAs (miRNAs) are endogenous 22 nt non-coding RNAs which play important regulatory roles in animals and plants by targeting 3′UTR of mRNAs for cleavage or translational repression [10, 11]. Some SNPs in pre-microRNAs, flanking regions or target sites have been demonstrated to affect certain physiological processes or related to diseases . In the present study, we found one possible valuable SNP in the mature miR-379, in which could potentially affect the binding ability to PKM2. We postulated that this SNP might contribute to the various expression of PKM2 within invidious and further affects the metabolism and growth of the tumor.
Materials and methods
Clinical sample information
Clinical characteristic of gastric cancer patients and cancer-free controls
Cases (n = 871)
Controls (n = 812)
H. pylori infection
Tumor size (cm)
Cell lines and cell culture
Gastric cancer cell lines including MKN-45 and AGS were purchased from American Type Culture Collection (ATCC). All cells were cultured in Dulbecco modified Eagle medium (DMEM) purchased from Gibco (CA, USA) supplemented with 10% fetal bovine serum (Invitrogen, Carlsbad, USA) and maintained in humidified 5% CO2 at 37 °C.
Construction of luciferase-based reporter plasmids
The full length of PKM2 cDNA as well as its 3′UTR were synthesized and sub-cloned into pGL3 plasmids. For the gene promoter activity assays, 3′UTR of PKM2 was synthesized and sub-cloned into a pGL3 plasmid (Promega, WI, USA). The construction containing different genotype of miR-379 was also synthesized and cloned into pSilence 2.1-U6. All the DNA synthesis and clones were performed in Genscript Co. (Nanjing, China).
Dual-luciferase reporter assay
The treated cells harvested 48 h after miRNA treatment, and the firefly luciferase expression was measured and normalized to Renilla activities. Dual-luciferase assays (Promega, Madison, WI) were performed according to the manufacturer’s protocol and detected with a Fluoroskan microplate reader (Thermo Labsystems, Helsinki, Finland). Transfection was repeated three times in triplicate (Additional file 1).
Cell proliferation assays
Cell proliferation was determined by using CCK-8 (Dojin Laboratories, Kumamoto, Japan) according to the manufacturer’s instructions. Briefly, the control and infected cells were seeded at a density of 1 × 103 cells/well in 96-well plates. 10 μL of CCK-8 was added to each well containing 100 µL of the culture medium, and the plate was incubated for 2 h at 37 °C. The viability of cells was evaluated by measuring the absorbance at 450 nm, using a microplate reader (Thermo Labsystems, CA).
Genomic DNA was extracted from peripheral blood by using QIAamp DNA blood mini kits (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Genotyping was performed with the TaqMan SNP Genotyping Assay. The PCR reactions were carried out in a total volume of 5 μL containing TaqMan Universal Master Mix, SNP Genotyping AssayMix, DNase-free water and genomic DNA. The PCR conditions were 2 min at 50 °C, 10 min at 95 °C, followed by 40 cycles at 95 °C for 15 s and 60 °C for 1 min. The 384-well ABI 7900HT real-time PCR system was applied (ABI, CA, USA).
Sections were stained according to the previous publication . The section was incubated within primary mouse anti-human Ab for PKM2 (ab38237), the sections were stained with DAB according to manufacturer’s protocols and mounted and photographed using a digitalized microscope camera (Nikon, Tokyo, Japan).
Plasma samples preparation and analysis by GC–TOF–MS
Metabolites extracted from plasma samples were analyzed using an Agilent 7890N gas chromatograph coupled with a Pegasus HT TOF mass spectrometer (Leco Corporation). Briefly, a 1 μL aliquot of the derivatized solution was injected with the splitless mode. Rxi-5 ms capillary column (30 m × 250 μm I.D., 0.25-μm film thickness; Restek Corporation, Bellefonte, PA, USA) was used for metabolites separation, with helium as the carrier gas at a constant flow rate of 1.0 mL/min. The temperature settings for injection, transfer interface, and ion source were 260, 260, and 210 °C, respectively. The separation was achieved with the following GC temperature program: 80 °C for 2 min, 10 °C/min to 220 °C, 5 °C/min to 240 °C, and 25 °C/min to 290 °C, and kept at 290 °C for 8 min. The data was collected with full scan mode (m/z 40–600), and an acquisition rate of 20 spectra/s. Electron impact ionization (70 eV) was used.
Metabolomics data analysis
The data from GC–TOFMS was processed with ChromaTOF software (v4.22, Leco Co., CA, USA). Compound annotation was performed by comparing the mass fragments with NIST 08 Standard mass spectral databases with a similarity of more than 70% and finally verified by available reference standards. The annotated compounds from GC–TOFMS were imported to SIMCA-P software 12.0.1 (Umetrics, Umeå, Sweden) for multivariate statistical analysis. Supervised orthogonal partial least squares-discriminant analysis (OPLS-DA) was used to compare between groups. Differential metabolites were selected based on the criteria of variable importance in the projection (VIP) > 1 in OPLS-DA model and P value < 0.05 from Student’s t test.
All experiments were performed in triplicate and repeated at least three times. Data were expressed as mean ± SD. The association between rs61991156 genotypes and the risk of GC was evaluated by calculating the odds ratios (ORs) and their 95% confidence intervals (CIs) using univariate and multivariate logistic regression analysis. Differences between two independent groups were tested with Student’s t test. All statistical analyses were carried out using SPSS version 18.0 and presented with Graph-pad prism software. Kaplan–Meier survival curves were plotted, and the log-rank test was done. The significance of various variables for survival was analyzed by the Cox proportional hazards model in a multivariate analysis. The results were considered to be statistically significant at P < 0.05.
Clinical significance of rs61991156 in human gastric cancer
Total 871 GC cases and 812 healthy controls were involved in our study, and their clinical characteristics were listed in Table 1. There is no significant difference in “age” and “Gender” between the case and control, significantly more H. pylori-infected cases were involved compared to the controls. The GC patients were divided into four groups according to the pathological differentiation degree staged from G1 to G4.
Genotype frequencies of the miR-379 at rs61991156 polymorphism among GC cases and controls
Cases (n = 871)
Controls (n = 821)
OR (95% CI)a
Stratified analysis of the correlation between miR-379 SNP and GC
Stratified analysis of SNP rs61991156 with clinicopathological parameters of GC
AG vs. GG
AA vs. GG
H. Pylori infection
Tumor size (cm)
The role of rs61991156 in regulatory effects of miR-379 on PKM2 expression and cell proliferation in vitro
The role of rs61991156 in regulatory glycolysis and other related metabolomics in GC patients
A/G SNP was associated with low expression of PKM2 and longer postoperative survival in gastric cancer
The general point of views on miR-379 were controversial. Some reports regarded it as a tumor suppressor which capable of down-regulating many oncogenes by targeting their 3′UTR region. For example, miR-379 was reported to be boosted by rifampicin and blocking the expression of ABCC2 (multidrug resistance-associated protein 2) thus to sensitize the HCC cells to chemotherapy . Others also indicated that miR-379 was overexpressed in malignancies, and might serve as an indicator for bad prognosis. miR-379 expression was elevated in bone-metastatic prostate cancer cell lines and tissues. The expression of miR-379 was also correlated with shortened progression-free survival of patients with prostate cancer . In our study, the results revealed that miR-379 was a tumor suppressor in human GC. From the clinical investigation, miR-379 GG genotype was associated with small tumor size, well differentiation, and non-metastasis which is related relatively low expression of PKM2 in gastric cancer. In line with the expression level, the glycolysis level within GC patients with GG genotype was also the weakest. We speculated that these results might result from rs61991156 in the mature form of miR-379.
In the present study, we found the SNP rs61991156 located within the mature form, generating three different genotypes of miR-379, among which A > G mutation generate an 8-mer complementary sequence in the 3′UTR of PKM2. The G allele might have significantly stronger the binding affinity to 3′UTR of PKM2 compared to the A allele. We thought this was the reason why miR-379 GG can significantly decrease the expression of PKM2 and in turn to attenuate both the proliferation and glycolysis of GC cells. Similar results concerning miRNA SNP have been reported previously resulting in either “Gain” or “LOSS” regulation of the targeting genes. miR-SNPs in miR-125a and Kaposi’s sarcoma-associated herpes virus-encoded miR-K5 were reported to impair miRNA processing by the Drosha/DGCR8 complex [16, 17]. SNP of miR-196a2 at rs11614913 in the mature miR-196a2 was reported to be associated with a significantly decreased rate of survival in individuals with non-small cell lung cancer, and the same research team of this study also suggested an association of rs11614913 with enhanced processing of mature miR-196a . miR-146a-SNP (rs2910164) within the pre-miR-146a sequence reduced both the amount of pre- and mature miR-146a and apparently affected the Drosha/DGCR8 processing step [18, 19]. miR-196a2-SNP, miR-146a-SNP, miR-149-SNP (rs2292832), and miR-499-SNP (rs3746444) are each associated with increased breast cancer risk .
So far, there is almost no report concerning the SNP of miR-379, we reported firstly that A > G SNP in 12nt of miR-379 might enhance the binding affinity to the 3′UTR of PKM2, thus to might be associated with low glycolysis level, well differentiation, as well as slower tumor growth. And the detection of rs61991156 might be associated with low occurrence and less aggressiveness of gastric cancer in Chinese population due to the enhanced regulations on PKM2.
Conception and design: XWW and NC. Collection and assembly of data: ML and NC; Data analysis and interpretation: YRW and ML. Contribution of reagents, materials, and analysis tools, wrote the paper: XWW and NC. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Morandi A, Taddei ML, Chiarugi P, Giannoni E. Targeting the metabolic reprogramming that controls epithelial-to-mesenchymal transition in aggressive tumors. Front Oncol. 2017;7:40.View ArticlePubMedPubMed CentralGoogle Scholar
- Schwartz L, Supuran CT, Alfarouk KO. The warburg effect and the hallmarks of cancer. Anticancer Agents Med Chem. 2017;17:164–70.View ArticlePubMedGoogle Scholar
- Zhang LF, Jiang S, Liu MF. MicroRNA regulation and analytical methods in cancer cell metabolism. Cell Mol Life Sci. 2017;74:2929–41.View ArticlePubMedGoogle Scholar
- Mitchell MI, Engelbrecht AM. Metabolic hijacking: a survival strategy cancer cells exploit? Crit Rev Oncol Hematol. 2017;109:1–8.View ArticlePubMedGoogle Scholar
- Schwartz L, Seyfried T, Alfarouk KO, Da Veiga Moreira J, Fais S. Out of Warburg effect: an effective cancer treatment targeting the tumor specific metabolism and dysregulated pH. Semin Cancer Biol. 2017;43:134–8.View ArticlePubMedGoogle Scholar
- Wu H, Wang Y, Wu C, Yang P, Li H, Li Z. Resveratrol induces cancer cell apoptosis through MiR-326/PKM2-mediated ER stress and mitochondrial fission. J Agric Food Chem. 2016;64:9356–67.View ArticlePubMedGoogle Scholar
- Liu AM, Xu Z, Shek FH, Wong KF, Lee NP, Poon RT, Chen J, Luk JM. miR-122 targets pyruvate kinase M2 and affects metabolism of hepatocellular carcinoma. PLoS ONE. 2014;9:e86872.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakao K, Miyaaki H, Ichikawa T. Antitumor function of microRNA-122 against hepatocellular carcinoma. J Gastroenterol. 2014;49:589–93.View ArticlePubMedGoogle Scholar
- Sun Y, Zhao X, Zhou Y, Hu Y. miR-124, miR-137 and miR-340 regulate colorectal cancer growth via inhibition of the Warburg effect. Oncol Rep. 2012;28:1346–52.View ArticlePubMedGoogle Scholar
- SanGiovanni JP, SanGiovanni PM, Sapieha P, De Guire V. miRNAs, single nucleotide polymorphisms (SNPs) and age-related macular degeneration (AMD). Clin Chem Lab Med. 2017;55:763–75.View ArticlePubMedGoogle Scholar
- Simchovitz A, Heneka MT, Soreq H. Personalized genetics of the cholinergic blockade of neuroinflammation. J Neurochem. 2017;142(Suppl 2):178–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Rajendran A, Chatterjee A, Pan A. Computational approaches and related tools to identify MicroRNAs in a species: A Bird’s Eye View. Interdiscip Sci. 2017; (Epub ahead of print).Google Scholar
- Cai C, Ashktorab H, Pang X, Zhao Y, Sha W, Liu Y, Gu X. MicroRNA-211 expression promotes colorectal cancer cell growth in vitro and in vivo by targeting tumor suppressor CHD5. PLoS ONE. 2012;7:e29750.View ArticlePubMedPubMed CentralGoogle Scholar
- Rigalli JP, Ciriaci N, Arias A, Ceballos MP, Villanueva SS, Luquita MG, Mottino AD, Ghanem CI, Catania VA, Ruiz ML. Regulation of multidrug resistance proteins by genistein in a hepatocarcinoma cell line: impact on sorafenib cytotoxicity. PLoS ONE. 2015;10:e0119502.View ArticlePubMedPubMed CentralGoogle Scholar
- Gururajan M, Josson S, Chu GC, Lu CL, Lu YT, Haga CL, Zhau HE, Liu C, Lichterman J, Duan P, Posadas EM, Chung LW. miR-154* and miR-379 in the DLK1-DIO3 microRNA mega-cluster regulate epithelial to mesenchymal transition and bone metastasis of prostate cancer. Clin Cancer Res. 2014;20:6559–69.View ArticlePubMedPubMed CentralGoogle Scholar
- Duan R, Pak C, Jin P. Single nucleotide polymorphism associated with mature miR-125a alters the processing of pri-miRNA. Hum Mol Genet. 2007;16:1124–31.View ArticlePubMedGoogle Scholar
- Gottwein E, Cai X, Cullen BR. A novel assay for viral microRNA function identifies a single nucleotide polymorphism that affects Drosha processing. J Virol. 2006;80:5321–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Hu Z, Chen J, Tian T, Zhou X, Gu H, Xu L, Zeng Y, Miao R, Jin G, Ma H, Chen Y, Shen H. Genetic variants of miRNA sequences and non-small cell lung cancer survival. J Clin Invest. 2008;118:2600–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Shen J, Ambrosone CB, DiCioccio RA, Odunsi K, Lele SB, Zhao H. A functional polymorphism in the miR-146a gene and age of familial breast/ovarian cancer diagnosis. Carcinogenesis. 2008;29:1963–6.View ArticlePubMedGoogle Scholar
- Liang J, Chen P, Hu Z, Zhou X, Chen L, Li M, Wang Y, Tang J, Wang H, Shen H. Genetic variants in fibroblast growth factor receptor 2 (FGFR2) contribute to susceptibility of breast cancer in Chinese women. Carcinogenesis. 2008;29:2341–6.View ArticlePubMedGoogle Scholar