Abstract
Purpose: Mitochondrial glycerophosphate dehydrogenase (MGPDH) is the key enzyme connecting oxidative phosphorylation (OXPHOS) and glycolysis as well as a target of the antidiabetic drug metformin in the liver. There are no data on the expression and role of MGPDH as a metformin target in cancer. In this study, we evaluated MGPDH as a potential target of metformin in thyroid cancer and investigated its contribution in thyroid cancer metabolism.
Experimental Design: We analyzed MGPDH expression in 253 thyroid cancer and normal tissues by immunostaining and examined its expression and localization in thyroid cancer–derived cell lines (FTC133, BCPAP) by confocal microscopy. The effects of metformin on MGPDH expression were determined by qRT-PCR and Western blot analysis. Seahorse analyzer was utilized to assess the effects of metformin on OXPHOS and glycolysis in thyroid cancer cells. We analyzed the effects of metformin on tumor growth and MGPDH expression in metastatic thyroid cancer mouse models.
Results: We show for the first time that MGPDH is overexpressed in thyroid cancer compared with normal thyroid. We demonstrate that MGPDH regulates human thyroid cancer cell growth and OXPHOS rate in vitro. Metformin treatment is associated with downregulation of MGPDH expression and inhibition of OXPHOS in thyroid cancer in vitro. Cells characterized by high MGPDH expression are more sensitive to OXPHOS-inhibitory effects of metformin in vitro and growth-inhibitory effects of metformin in vitro and in vivo.
Conclusions: Our study established MGPDH as a novel regulator of thyroid cancer growth and metabolism that can be effectively targeted by metformin. Clin Cancer Res; 24(16); 4030–43. ©2018 AACR.
The antidiabetic drug metformin has been shown to be an effective anticancer agent in numerous studies. We previously demonstrated that metformin-treated patients with diabetes were characterized by smaller thyroid cancer size, higher complete remission rate, and longer progression-free survival than diabetics not treated with metformin. In this study, we provide evidence that MGPDH is a target of metformin in thyroid cancer. We demonstrate for the first time that not only does MGPDH regulate human thyroid cancer cell growth and metabolism, but the susceptibility to the growth-inhibitory effects of metformin in thyroid cancer depends on the level of MGPDH expression. Our observations may have a significant impact for an individualized approach to cancer therapy as high tumor tissue MGPDH expression may predict susceptibility to the growth-inhibitory effects of metformin in vivo. The study also creates a scope to test metformin efficacy in other cancers characterized by high MGPDH expression.
Introduction
Epidemiologic evidence suggests that therapy with the antidiabetic drug metformin is associated with decreased incidence of certain cancers, such as colon, liver, or lung, as well as improved response to cancer therapy and decreased cancer mortality (1). However, there is a significant heterogenicity and discrepancy between these studies; some document an association between metformin use and beneficial cancer treatment outcome with reduced mortality (2–5), while others fail to document such beneficial effects (6, 7). These observations suggest the presence of a specific molecular signature of cancer that increases its susceptibility to the antineoplastic effects of metformin.
Metformin is a hydrophilic, positively charged drug under physiologic conditions and, thus, requires active intracellular transport through organic cation transporters (OCT) 1, 2, or 3. OCT1 is highly expressed in the major target tissue of metformin action—the liver (8). Mitochondrial glycerophosphate dehydrogenase (MGPDH) has been found to be a metformin target responsible for hepatic gluconeogenesis inhibition (9). MGPDH is a flavin-linked respiratory chain dehydrogenase which, via the glycerol-3-phosphate shuttle, connects glycolysis with oxidative phosphorylation (OXPHOS). Mammalian MGPDH is encoded by a single, highly conserved GPD2 gene, located on human chromosome 2q24.1 (10).
Although glycolysis and OXPHOS are the two major metabolic adaptation pathways in cancer (11), there are no data on the role of MGPDH as a metformin target in cancer or its contribution in cancer cell metabolism. To analyze the role of MGPDH in cancer metabolism, we utilized thyroid cancer as a model system. Currently, thyroid cancer is the most common endocrine malignancy, with an incidence increasing faster than any other cancer type (12). We used two human thyroid cancer cell line models derived from follicular and papillary thyroid cancer tissues (13).
We previously documented that thyroid cancer in metformin-treated patients with diabetes is characterized by smaller tumor size, higher complete remission rate, and longer progression-free survival than in patients with diabetes not treated with metformin (14). We investigated the pathophysiology of this association by studying in vitro models of human thyroid cancer and documented that the growth-inhibitory effects of metformin were due to downregulation of the mTOR signaling pathway (14). Interestingly, we observed a differential susceptibility of different thyroid cancer cell lines to the antiproliferative effects of metformin, and showed that the availability of metabolic substrates (i.e., glucose) modifies the response to metformin in vitro (15). This observation formed the rationale to test the role of MGPDH in growth and metabolism of thyroid cancer cell lines in vitro, metastatic mouse models in vivo, and in a transgenic mouse model that spontaneously develops thyroid cancer.
In this study, we document for the first time that MGPDH is overexpressed in thyroid cancer compared with normal thyroid tissue. We show that MGPDH regulates thyroid cancer cell growth and mitochondrial metabolism—with MGPDH overexpression associated with increased growth and OXPHOS rate, and, conversely, decreased proliferation and mitochondrial respiration with MGPDH downregulation. Furthermore, we provide evidence that MGPDH is a metformin target in thyroid cancer.
Methods
Cell lines and culture conditions
Thyroid cancer cell lines FTC133 [male-derived, follicular thyroid cancer (FTC) with a PTEN R130STOP and TP53 R273H mutation] and BCPAP [female-derived, papillary thyroid cancer (PTC) with a BRAFV600E and TP53 D259Y mutation] were utilized (9, 13). The cell lines were authenticated by short tandem repeat: >80% FTC133; 100% BCPAP (Appendix 2). The cells were grown in DMEM-high glucose medium (Gibco) supplemented with 10% FBS (Thermo Fisher Scientific), 2 μg/mL insulin (Thermo Fisher Scientific), 1 IU/100 mL thyrotropic hormone (Sigma Aldrich), 10 U/mL penicillin/streptomycin (Gibco), and 0.25 μg/mL Amphotericin B (Gibco; ref. 16).
Cells were treated with 1 mmol/L and 5 mmol/L metformin (Sigma Aldrich) for 24 and 48 hours, and 50, 100, and 200 nmol/L concentrations of T3 (Sigma Aldrich) for 48 and 72 hours, and combined therapy with metformin 5 mmol/48 hours and T3 100 nmol/L/72 hours.
Luciferase-transfected FTC133 and BCPAP cells were used for in vivo studies (17). Cells were transfected with a linearized pGL4.51[luc2/CMV/Neo] vector (Promega) encoding the luciferase reporter gene luc2 (Photinus pyralis) using Lipofectamine 2000 reagent (Invitrogen) and selected with Geneticin (600 μg/mL)-containing culture media and the G418-resistant cells were amplified in the same medium. IMPACT II PCR analysis was performed on the cells to check for any contamination and the test results were reported to be negative (Appendix 3).
siRNA transfection
FTC133 and BCPAP cells were transfected with commercially available MGPDH siRNA (hs.Ri.MGPDH.13.2, Integrated DNA Technologies) or negative control (NC) siRNA (51-01-14-04, Integrated DNA Technologies) using Lipofectamine RNAiMAX (13778075, Invitrogen) as the transfection agent. Cells were transfected with 100 pmoles si-MGPDH or si-NC. qRT-PCR and Western blot analysis demonstrated successful silencing at 48 hours posttransfection.
For Seahorse assay, cells were transfected using a reverse transfection protocol. To transfect all the wells, a complex of siRNA and Lipofectamine RNAiMAX was prepared in a optiMEM I medium (31985088, Gibco). Cells were treated with metformin for cellular energetic studies.
CRISPR/Cas9 gene editing
MGPDH gene knockdown in FTC133 and BCPAP cells was accomplished utilizing commercially available pCas guide vector and donor template DNA containing homologous arms and functional cassette (KN213341, OriGene). OriGene protocol was followed to transfect cells with guide RNA (1 μg) and donor template (1 μg) using turbofectin 8.0 (TF81001, OriGene). The puromycin-resistant cells were analyzed for MGPDH (monoallelic deletion) knockdown by Western blot analysis.
Stable overexpressing cell lines
To generate MGPDH-overexpressing FTC133 and BCPAP cell lines, commercially available MGPDH (EX-A3080-M03)-expressing plasmid along with an empty control vector (EX-NEG-M03) was purchased from GeneCopoeia. Cells were transfected with 1 μg of MGPDH-expressing plasmid DNA or empty control using Turbofectin 8.0 as the transfection agent. Plasmid DNA–Turbofectin complex was prepared in optiMEM I medium. Cells were subjected to G418 selection to obtain resistant cell colonies. qRT-PCR and Western blot analysis confirmed MGPDH overexpression.
Cell proliferation assay
Cell proliferation assay was performed according to the Vybrant MTT Cell Proliferation Assay Kit (V13154, Invitrogen) protocol.
RNA extraction and quantitative real-time PCR
Total RNA was isolated using RNeasy Mini Kit (Qiagen). cDNA was synthesized using iScript cDNA Synthesis Kit (Bio-Rad). Gene expression was analyzed using iQ SYBR Green Supermix (Bio-Rad) on a StepOnePlus Real-Time PCR Detection System (Applied Biosystems) using commercially available primers for MGPDH (Hs_MGPDH_1_SG QuantiTect Primer Assay) and β-actin (Hs_ACTB_1_SG QuantiTect Primer Assay). The custom-made primers (Integrated DNA technologies) for RT-PCR analysis of other genes were as follows: OCT1 [forward primer (FP)-GGATGAAGATGGAGTTATC, reverse primer (RP)- GTAGGCAAGTATGAGGAA]; COX6B2 (FP-CTAACGCCCTCACCATTT, RP-TCGCATCATTGTACTCCAG); NDUFA6 (FP- TTCCTATCCAAGTTCTATGTT, RP-AAGTGACCATTGTATAGTGA); and NDUFC2 (FP-ATTGAGTAGTTGTGATAGA, RP- ACAGAATCCTCCATAATC).
Mitochondrial metabolism PCR array
The PCR array was performed using an RT² Profiler PCR Array Human Mitochondrial Energy Metabolism Kit (PAHS-008Z; Qiagen) as per manufacturer instructions.
Mitochondrial DNA copy-number analysis
The average mitochondrial DNA (mtDNA) copy-number was determined using a NovaQUANT Human Mitochondrial to Nuclear DNA Ratio Kit (72620; Millipore Sigma) as per the manufacturer's instructions.
Protein extraction and immunoblot analysis
Cells were lysed using M-PER (Mammalian Protein Extraction Reagent; Thermo Scientific) as per manufacturer's instructions. Amicon Centrifugal filters (Millipore) were used for protein concentration. Following protein estimation by DC Protein Assay Kit II, samples were loaded onto 4% to 12% SDS-PAGE gel (Novex) and transferred to a membrane utilizing a Trans-blot turbo system (Bio-Rad). The membranes were blocked with 3% BSA for 1 hour and incubated overnight with a primary antibody against MGPDH (1:700; Abcam; ab188585), OCT1 (SLC22A1; 1:1,000; Abcam; ab181022) and total OXPHOS (1:500; Abcam; ab110411). GAPDH (1:1,000; Invitrogen; MA5-15738) and Porin (1:1,000; Abcam; ab14734) were used as endogenous controls. After 3× washing with TBST, the membrane was incubated with the secondary antibody (1:10,000) and washed three times. The blot was visualized using ECL (Thermo Scientific) and imaging system (Bio-Rad). The relative protein expression was quantified by densitometry using Fiji ImageJ software.
Mitochondrial stress assay
Mitochondrial function was determined by measuring oxygen consumption rate (OCR) of each cell line using XF Cell Mito Stress Test Kit (Agilent Technologies). FTC133 and BCPAP cells were seeded in an XF96 cell culture microplate and treated with metformin for 48 hours. The sensor cartridge and base medium were prepared by adding 1 mmol/L pyruvate, 2 mmol/L glutamine, and 10 mmol/L glucose and stored as per the manufacturer's instructions. Seahorse assay was run in XF96 Extracellular Flux Analyzer (Agilent Technologies). Following three baseline OCR measurements, cells were exposed sequentially to Oligomycin (0.5 μmol/L), Carbonyl cyanide-4 (trifluoromethoxy) phenylhydrazone (FCCP; 1 μmol/L), and rotenone/antimycin A (0.5 μmol/L). Oligomycin inhibits ATP synthase (complex V), and the decrease in OCR following injection of oligomycin correlates to the mitochondrial respiration associated with cellular ATP production. FCCP is an uncoupling agent that collapses the proton gradient and disrupts the mitochondrial membrane potential, allowing cells to achieve maximal OCR. As a result, electron flow through the electron transfer chain (ETC) is uninhibited and oxygen is maximally consumed by complex IV. The FCCP-stimulated OCR can then be used to calculate spare respiratory capacity, defined as the difference between maximal respiration and basal respiration. Spare respiratory capacity is a measure of the ability of the cell to respond to increased energy demand. The third injection is a mix of rotenone, a complex I inhibitor, and antimycin A, a complex III inhibitor. This combination shuts down mitochondrial respiration and enables the calculation of nonmitochondrial respiration driven by processes outside the mitochondria. Three measurements were recorded after every injection. To normalize the results, cell number/well was quantified using the Celígo Imaging Cytometer (Nexcelom). The assay results were analyzed using Wave program 2.3.0 (Seahorse Bioscience).
Glycolysis stress assay
The conversion of glucose to pyruvate, and subsequently lactate, results in a net production and extrusion of protons into the extracellular medium and consequently acidification of the medium. The Seahorse XF Glycolysis Stress Test Kit (Agilent Technologies) was utilized to measure the extracellular acidification rate (ECAR) because of glycolytic metabolism. Cells were seeded in a XF96 cell culture microplate and treated with metformin for 48 hours. The sensor cartridge and assay medium preparation was performed as per manufacturer's instructions. Following three baseline ECAR measurements, cells were exposed sequentially to glucose (10 mmol/L), oligomycin (1.0 μmol/L), and 2-deoxy-glucose (50 mmol/L). The first injection is a saturating concentration of glucose. The cells utilize the glucose injection and catabolize it through the glycolytic pathway to pyruvate, producing ATP, NADH, water, and protons. The extrusion of protons into the surrounding medium causes a rapid increase in ECAR. This glucose-induced response is reported as the rate of glycolysis under basal conditions. The second injection is oligomycin, an ATP synthase inhibitor. Oligomycin inhibits mitochondrial ATP production, and shifts the energy production to glycolysis, with the subsequent increase in ECAR revealing the cellular maximum glycolytic capacity. The final injection is 2-deoxy-glucose (2-DG), a glucose analogue, which inhibits glycolysis through competitive binding to glucose hexokinase, the first enzyme in the glycolytic pathway. The resulting decrease in ECAR confirms that the ECAR produced in the experiment is due to glycolysis. ECAR, prior to glucose injection, is referred to as nonglycolytic acidification. Three measurements were recorded after every injection. The results were analyzed and normalized as described above.
MGPDH enzyme assay
Mitochondria were isolated from cultured FTC133 and BCPAP cells using a commercially available Mitochondrial Isolation Kit (Thermo Fisher Scientific), according to the manufacturer's instructions. MGPDH enzymatic activity was measured as rate of reduction of cytochrome c following the addition of glycerophosphate. The assay was performed in incubation medium containing 50 mmol/L KCl, 1 mg/mL BSA, 10 mmol/L Tris-HCl, 1 mmol/L EDTA, 1 mmol/L KCN, 50 μmol/L cytochrome c, pH 7.4. Prior to the assay, isolated mitochondria were incubated with either 1 mmol/L or 5 mmol/L metformin for 10 minutes. The reaction was initiated by the addition of 20 mmol/L glycerophosphate and gain of absorbance was measured at 550 nm for 5 minutes at 30°C. The enzyme activity was measured as rate of change of absorbance/minute and presented as activity percentage relative to control.
Complex I activity assay
FTC133 and BCPAP cells were treated with either 1 or 5 mmol/L metformin for 4 hours. The enzyme activity of complex I was measured according to the Complex I Enzyme Activity Microplate Assay Kit protocol (Abcam; ab109721).
In vivo studies
The animal protocol for the experiments was approved by NCI-Bethesda Animal Care and Use Committee. Xenograft NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ mice, purchased from Jackson Laboratory were used for in vivo studies. At 8 to 10 weeks of age, mice received tail vein injection of either FTC133 (1.6 × 105 cells) or BCPAP (5 × 105 cells) cell lines containing a linearized pGL4.51[luc2/CMV/Neo] vector.
Mice were treated with metformin (12.5 mg/250 μL) or water via oral gavage daily for 4 weeks. Tumor burden was checked weekly through Xenogen IVIS (In Vivo Imaging System; Caliper Life Sciences Inc.). Mice were injected with luciferin approximately 15 minutes prior to imaging. After the last whole-body imaging, the mice were euthanized by CO2 inhalation, and the organs were isolated and fixed in formalin for histologic analysis. All in vivo experiments were performed twice.
Immunostaining
Immunostaining was performed on paraffin-embedded tissue sections. Sections were dewaxed, soaked in alcohol, and, after microwave treatment in antigen unmasking solution (Vector Laboratories) for 10 minutes, incubated in 3% hydrogen peroxide for 15 minutes to inactivate endogenous peroxidase activity. Sections were incubated for 10 minutes in working solution of blocking serum and incubated at 4°C overnight with anti-MGPDH (2 mg/mL; Abcam; catalog No. ab188585) or with anti-OCT1 (2 mg/mL; Abcam; catalog No. 182144). Immunostaining was performed using the Vectastain Universal Quick Kit according to manufacturer's instructions. Sections were incubated with biotinylated secondary antibody for 10 minutes and in streptavidin/peroxidase complex working solution for 10 minutes. Peroxidase staining was revealed with ImmPact DAB, peroxidase substrate Kit (Vector Laboratories). Sections were counterstained with hematoxylin and mounted. Antiserum was omitted in the negative control.
The results of staining were interpreted and scored independently by two investigators (V.V. Vasko and J. Klubo-Gwiezdzinska).
Immunofluorescence staining and confocal microscopy
The 2.5 × 105 cells were plated in glass coverslips in 6-well plates and allowed to attach for 24 hours at 37°C (5% CO2). The cells were incubated in 1:5,000 concentration of Mitotracker Deep Red (Invitrogen) for 15 minutes, washed with prewarmed media, then fixed and permeabilized in 100% methanol for 5 minutes at 4°C. Next, the cells were blocked with 5% BSA in PBS-TT (0.5% Tween and 0.1% Triton) for 1 hour and immunostained for mGPDH using mAb anti-MGPDH (Abcam) for 2 hours then secondary antibodies conjugated with Alexa Fluor 488 or Alexa Fluor 569 for 1 hour (Life Technologies, Thermo Fisher Scientific). DNA was stained with DAPI (Vector Laboratories). Fluorescence images were detected by confocal microscope NLO 710 Zeiss, and collected using Carl Zeiss Zen Software (Zeiss).
Widefield fluorescence microscopy and image analysis
Images were acquired on a Keyence BZ-9000 digital microscope (Keyence Corporation of America), using the equipped Nikon 10× objective lens (Nikon Instruments, Inc.). Filter sets used for DAPI and Alexa 488 were 49028 and 49002, respectively (Chroma Technology Corporation).
Quantification of green and red signal in the images was performed using a custom macro written for Fiji/ImageJ software (18). Cell counts/image, based on DAPI nuclei identification, ranged from 141 to 977. Intensity measurements were divided by the number of nuclei to get per-cell readings. For each experimental group, a background/cell was determined from the negative control and subtracted from the other per-cell measurements in the group.
Statistical analysis
Data were analyzed using GraphPad Prism 7.0 software. For noncategorical data, a two-tailed unpaired Student t test was used. Data are presented as mean ± SD. Seahorse graphs are presented with SE bars. For categorical data, a two-tailed Fisher exact test was used. P value of <0.05 was considered significant. Bonferroni correction was used for multiple comparisons. All in vitro experiments were performed in two to 18 technical replicates (n represents technical replicates) and repeated with at least two biological replicates while in vivo experiments involved eight to 14 biological replicates and were performed twice.
Results
Metformin intracellular transporter, OCT1 (SLC22A1), is overexpressed in thyroid cancer compared with normal thyroid tissue
Metformin transport into cells involves an active uptake process mediated by organic cation transporter (OCT) 1, 2, or 3 (8). To determine the level of OCT1 expression in thyroid cancer, we analyzed RNA-Seq data derived from The Cancer Genome Atlas (TCGA) database incorporating 507 patients with thyroid cancer. We examined the OCT1 expression in different subtypes of papillary thyroid cancer (Supplementary Fig. S1A). TCGA database includes 65 paired samples of tumor versus normal tissue derived from the same patients. Analysis of these paired samples demonstrated significantly higher OCT1 expression in thyroid cancer compared with normal thyroid (log fold change 0.3; P = 0.04; raw data—Appendix 1).
We next performed IHC staining for OCT1 in 200 tissue samples—27 normal thyroid, 69 follicular thyroid cancer (FTC), 68 papillary thyroid cancer (PTC), and 36 poorly differentiated thyroid cancer (PDTC). OCT1 protein level was significantly higher in thyroid cancer compared with normal thyroid (P < 0.001), with the highest level observed in PTC (Fig. 1A).
Expression of OCT1 (SLC22A1) and MGPDH in human thyroid cancer tissues and cell lines. A, Significantly higher expression of OCT1 in thyroid cancer compared with normal thyroid tissue as documented by IHC staining (n = 200, P < 0.001). The expression of OCT1 is significantly higher in papillary thyroid cancer compared with follicular (P < 0.001) and poorly differentiated thyroid cancer (P = 0.004), which are characterized by similar OCT1 expression. NT, normal thyroid; FTC, follicular thyroid cancer; PTC, papillary thyroid cancer; PDTC, poorly differentiated thyroid cancer. B, The scheme representing the MGPDH role in glycerol-3-phosphate shuttle. MGPDH transfers reducing equivalents from cytosolic NADH into the mitochondrial electron transport chain connecting glycolysis with OXPHOS. C, Significantly higher expression of MGPDH in thyroid cancer compared with normal thyroid tissue as documented by IHC staining (n = 207, P < 0.001). D, BCPAP cells express higher OCT1 mRNA (n = 4; P < 0.01) and protein than FTC133 cells. No effect of metformin 5 mmol/L/48 hours on OCT1 expression. E, Significantly higher mRNA expression of MGPDH in FTC133 compared with BCPAP human thyroid cancer cell line (n = 4, P < 0.0001). Significantly higher protein expression of MGPDH in FTC133 compared with BCPAP human thyroid cancer cell line (n = 4, P < 0.001). F, Intracellular mitochondria-associated localization of MGPDH in FTC133 cells documented by IF staining and confocal microscopy imaging, (n = 4-6). Blue, DAPI nuclear staining; purple, Mitotracker mitochondrial staining; green, MGPDH staining. G, Intracellular mitochondria-associated localization of MGPDH in BCPAP cells documented by IF staining and confocal microscopy imaging (n = 4–6). Blue, DAPI nuclear staining; purple, Mitotracker mitochondrial staining; green, MGPDH staining.
Expression of OCT1 (SLC22A1) and MGPDH in human thyroid cancer tissues and cell lines. A, Significantly higher expression of OCT1 in thyroid cancer compared with normal thyroid tissue as documented by IHC staining (n = 200, P < 0.001). The expression of OCT1 is significantly higher in papillary thyroid cancer compared with follicular (P < 0.001) and poorly differentiated thyroid cancer (P = 0.004), which are characterized by similar OCT1 expression. NT, normal thyroid; FTC, follicular thyroid cancer; PTC, papillary thyroid cancer; PDTC, poorly differentiated thyroid cancer. B, The scheme representing the MGPDH role in glycerol-3-phosphate shuttle. MGPDH transfers reducing equivalents from cytosolic NADH into the mitochondrial electron transport chain connecting glycolysis with OXPHOS. C, Significantly higher expression of MGPDH in thyroid cancer compared with normal thyroid tissue as documented by IHC staining (n = 207, P < 0.001). D, BCPAP cells express higher OCT1 mRNA (n = 4; P < 0.01) and protein than FTC133 cells. No effect of metformin 5 mmol/L/48 hours on OCT1 expression. E, Significantly higher mRNA expression of MGPDH in FTC133 compared with BCPAP human thyroid cancer cell line (n = 4, P < 0.0001). Significantly higher protein expression of MGPDH in FTC133 compared with BCPAP human thyroid cancer cell line (n = 4, P < 0.001). F, Intracellular mitochondria-associated localization of MGPDH in FTC133 cells documented by IF staining and confocal microscopy imaging, (n = 4-6). Blue, DAPI nuclear staining; purple, Mitotracker mitochondrial staining; green, MGPDH staining. G, Intracellular mitochondria-associated localization of MGPDH in BCPAP cells documented by IF staining and confocal microscopy imaging (n = 4–6). Blue, DAPI nuclear staining; purple, Mitotracker mitochondrial staining; green, MGPDH staining.
Two thyroid cancer cell lines, FTC133 (FTC) and BCPAP (PTC), served as a model for investigating thyroid cancer. We found moderate OCT1 expression at the mRNA (Fig. 1D) and protein level (Fig. 1D), with BCPAP cells characterized by higher OCT1 expression than FTC133 cells. There was no change in OCT1 expression after metformin treatment (Fig. 1D). OCT2 and OCT3 expression was not detected in the thyroid cancer cell lines (data not shown).
We show that the level of OCT1 is elevated in human thyroid cancer tissue samples. OCT1 is expressed in human FTC and PTC cell lines. These observations established a robust model for functional in vitro and in vivo studies to test the metabolic and growth-inhibitory effects of metformin in thyroid cancer.
MGPDH is overexpressed in thyroid cancer compared with normal thyroid tissue
While MGPDH has been found to be a target of metformin in the liver (9), there are no data on its expression and role in cancer cell metabolism and response to metformin. Using the RNA-Seq data from TCGA database consisting of 507 patients with PTC, we examined the MGPDH expression in all PTC subtypes (Supplementary Fig. S1B). Among the 65 paired samples of tumor versus normal tissue, MGPDH showed significantly higher expression in thyroid cancer compared with normal thyroid (log fold change 0.28, P < 0.001; raw data—Appendix 1). IHC staining in 207 tissue samples derived from 117 patients (72 FTC, 72 PTC, 36 PDTC and 27 normal thyroid tissues) showed significantly higher expression of MGPDH protein in thyroid cancer compared with normal thyroid tissue (P < 0.001; Fig. 1C). We next evaluated MGPDH expression in thyroid cancer cell lines by qRT-PCR, Western blot (WB), immunofluorescence staining, and its intracellular localization by confocal microscopy. We found mitochondrial expression of MGPDH based on costaining with mitochondrial marker, Mitotracker, in both thyroid cancer cell lines (Fig. 1F and G). FTC133 cells are characterized by significantly higher MGPDH expression than BCPAP cells (P < 0.001), (Fig. 1E; Supplementary Fig. S1C). FTC133 cells are also characterized by higher mitochondrial staining as documented by quantification of immunofluorescence staining with Mitotracker (Supplementary Fig. S1D).
In summary, we demonstrated that MGPDH is overexpressed in thyroid cancer, thus establishing in vitro models characterized by high (FTC133) and low (BCPAP) MGPDH expression to test the effect of metformin on MGPDH expression and cancer cell metabolism.
Metformin treatment decreases MGPDH expression in thyroid cancer cell lines in vitro
We previously demonstrated that metformin inhibits proliferation of thyroid cancer cell lines in vitro via activation of the intracellular energy sensor, AMPK (AMP-activated protein kinase), and downstream inhibition of the mTOR/pS6/cyclin D1 signaling pathway (14, 19). In this study, we tested whether the same growth-inhibitory and metabolic stress–inducing concentrations of metformin affect MGPDH expression, the key enzyme connecting glycolysis with OXPHOS via the glycerol-3-phosphate shuttle (Fig. 1B). Treatment with 5 mmol/L metformin for 48 hours resulted in significant MGPDH downregulation in the two examined thyroid cancer cell lines at the mRNA (P = 0.001; Fig. 2A and B), and protein level (Fig. 2C). We also assessed the proliferation rate after 5 days of exposure to metformin 5 mmol/L. We found a significant inhibition of proliferation by 58.01% ± 7.273% in FTC133 cells and 61.96% ± 9.345% in BCPAP cells (Supplementary Fig. S2G). However, longer incubation with metformin was not associated with a more significant decrease in MGPDH protein expression in comparison with 48-hour exposure (Supplementary Fig. S2H).
FTC133 cells are characterized by higher level of MGPDH expression and higher mitochondrial respiration rate than BCPAP cells. Metformin treatment is associated with reduction in MGPDH expression in both cell lines. A, Decrease in mRNA expression of MGPDH in FTC133 cells treated with 5 mmol metformin for 48 hours. M1, metformin 1 mmol/48 hours; M5, metformin 5 mmol/48 hours (n = 4, P = 0.001). B, Decrease in mRNA expression of MGPDH in BCPAP cells treated with 5 mmol metformin for 48 hours. M1, metformin 1 mmol/48 hours; M5, metformin 5 mmol/48 hours (n = 4, P = 0.002). C, Reduction of MGPDH protein expression after treatment with metformin (1 and 5 mmol/48 hours) in FTC133 and BCPAP cells documented by Western blot analysis. M1, treatment with metformin 1 mmol/48 hours; M5, treatment with metformin 5 mmol/48 hours; P, positive control recombinant human MGPDH. D, Higher baseline MGPDH activity in FTC133 cell compared with BCPAP cells (n = 3, P < 0.001). Significant reduction of MGPDH activity after treatment with metformin (n = 3). M1, treatment with metformin 1 mmol/10 min (FTC133-P < 0.05; BCPAP-P < 0.01); M5, treatment with metformin 5 mmol/10 min (FTC133-P < 0.01; BCPAP-P < 0.001). E, Significantly higher baseline mitochondrial respiration, ATP production, maximal respiration and spare respiratory capacity in FTC133 compared with BCPAP cells documented by Seahorse Mitostress Kit (n = 6–9). Basal OCR – FTC133 49.2 ± 4.6 pmol/min/cell versus BCPAP 7.4 ± 1.3 pmol/min/cell (P < 0.0001). ATP production – FTC133 40.6 ± 3.6 pmol/min/cell versus BCPAP 3.5 ± 1.6 pmol/min/cell (P < 0.0001). Maximum respiration – FTC133 132 ± 14.2 pmol/min/cell versus BCPAP 7.8 ± 2.2 pmol/min/cell (P < 0.0001). Spare respiratory capacity – FTC133 82.8 ± 12.2 pmol/min/cell versus BCPAP 2 ± 1.3 pmol/min/cell (P < 0.0001). F, Comparable level of glycolysis, glycolytic reserve and capacity of FTC133 and BCPAP cells measured by Seahorse Glycolysis Stress Kit, (n = 6–9). Glycolysis rate – FTC133 26.8 ± 4.6 mpH/min/cell versus BCPAP 29.9 ± 3 mpH/min/cell (P = 0.18). Glycolytic capacity – FTC133 51.8 ± 7.5 mpH/min/cell versus BCPAP 46.6 ± 5.7 mpH/min/cell (P = 0.21). Glycolytic reserve – FTC133 25.7 ± 3.4 mpH/min/cell versus BCPAP 16.7 ± 3.4 mpH/min/cell (P = 0.1); ECAR, extracellular acidification rate mpH/min/cell.
FTC133 cells are characterized by higher level of MGPDH expression and higher mitochondrial respiration rate than BCPAP cells. Metformin treatment is associated with reduction in MGPDH expression in both cell lines. A, Decrease in mRNA expression of MGPDH in FTC133 cells treated with 5 mmol metformin for 48 hours. M1, metformin 1 mmol/48 hours; M5, metformin 5 mmol/48 hours (n = 4, P = 0.001). B, Decrease in mRNA expression of MGPDH in BCPAP cells treated with 5 mmol metformin for 48 hours. M1, metformin 1 mmol/48 hours; M5, metformin 5 mmol/48 hours (n = 4, P = 0.002). C, Reduction of MGPDH protein expression after treatment with metformin (1 and 5 mmol/48 hours) in FTC133 and BCPAP cells documented by Western blot analysis. M1, treatment with metformin 1 mmol/48 hours; M5, treatment with metformin 5 mmol/48 hours; P, positive control recombinant human MGPDH. D, Higher baseline MGPDH activity in FTC133 cell compared with BCPAP cells (n = 3, P < 0.001). Significant reduction of MGPDH activity after treatment with metformin (n = 3). M1, treatment with metformin 1 mmol/10 min (FTC133-P < 0.05; BCPAP-P < 0.01); M5, treatment with metformin 5 mmol/10 min (FTC133-P < 0.01; BCPAP-P < 0.001). E, Significantly higher baseline mitochondrial respiration, ATP production, maximal respiration and spare respiratory capacity in FTC133 compared with BCPAP cells documented by Seahorse Mitostress Kit (n = 6–9). Basal OCR – FTC133 49.2 ± 4.6 pmol/min/cell versus BCPAP 7.4 ± 1.3 pmol/min/cell (P < 0.0001). ATP production – FTC133 40.6 ± 3.6 pmol/min/cell versus BCPAP 3.5 ± 1.6 pmol/min/cell (P < 0.0001). Maximum respiration – FTC133 132 ± 14.2 pmol/min/cell versus BCPAP 7.8 ± 2.2 pmol/min/cell (P < 0.0001). Spare respiratory capacity – FTC133 82.8 ± 12.2 pmol/min/cell versus BCPAP 2 ± 1.3 pmol/min/cell (P < 0.0001). F, Comparable level of glycolysis, glycolytic reserve and capacity of FTC133 and BCPAP cells measured by Seahorse Glycolysis Stress Kit, (n = 6–9). Glycolysis rate – FTC133 26.8 ± 4.6 mpH/min/cell versus BCPAP 29.9 ± 3 mpH/min/cell (P = 0.18). Glycolytic capacity – FTC133 51.8 ± 7.5 mpH/min/cell versus BCPAP 46.6 ± 5.7 mpH/min/cell (P = 0.21). Glycolytic reserve – FTC133 25.7 ± 3.4 mpH/min/cell versus BCPAP 16.7 ± 3.4 mpH/min/cell (P = 0.1); ECAR, extracellular acidification rate mpH/min/cell.
Because we found mitochondrial expression of MGPDH, we next examined whether MGPDH expression inhibition could be due to decreased intracellular mitochondrial content. Metformin therapy did not affect mitochondrial staining in any of the examined cell lines, as documented by the quantification of immunofluorescence staining with Mitotracker (Supplementary Fig. S1D) and by mtDNA copy-number analysis (Supplementary Fig. S1E).
We also performed a PCR expression profile array of other genes involved in mitochondrial metabolism (Supplementary Table S4). The only genes with mRNA expression, which was nearly twofold downregulated in both cell lines after metformin treatment (5 mmol/48 hours), were ATP4B (subunit of ATPase synthase Complex), COX412 (subunit of complex IV), and COX6B2 (subunit of complex IV). Among these three, ATP4B and COX4I2 had very low expression in both cell lines with Ct values greater than 35. In FTC133 cells, we observed more than twofold downregulation in NDUFC2 and NDUFA6 genes, which are complex I subunits. We performed validation of COX6B2, NDUFC2, and NDUFA6 with targeted PCR, and found that none of the targets were consistently downregulated in both cell lines (Supplementary Fig. S2A). We also performed Western blot analysis using total OXPHOS Human WB Antibody Cocktail enabling concomitant analysis of main proteins of each complex in the electron transport chain—complex I subunit NDUFB8, complex II subunit 30kDa, complex III subunit Core 2, complex IV subunit II, and ATP synthase subunit alpha. There was no difference in protein expression after treatment with metformin 5 mmol/L for 48 hours (Supplementary Fig. S2B). We also analyzed the enzymatic activity of MGPDH in mitochondria extracted from FTC133 and BCPAP cells. We found a significant decrease in MGPDH activity after 10 minutes of exposure to metformin (Fig. 2D). We did not find any difference in mitochondrial complex I activity in cancer cells after up to 4-hour exposure to metformin (Supplementary Fig. S2F). These results do not rule out the involvement of complex I activity in the intracellular effects of metformin in our model if it binds at or near the ubiquinone binding site. To the best of our knowledge, the binding site of metformin into the mitochondrial complex I is not elucidated and is beyond the scope of this article. The activity assay utilized in our experiments measured diaphorase-type activity of complex I. This activity is not dependent on the presence of ubiquinone and, therefore, inhibitors that bind at or near the ubiquinone binding site do not inhibit this assay.
Taken together, our data suggest that MGPDH may serve as a metformin target in thyroid cancer cells and formed the rationale to test the role of MGPDH in thyroid cancer metabolism and its response to metformin.
MGPDH-rich thyroid cancer cell line FTC133 is characterized by higher baseline OXPHOS rate compared with BCPAP cells with low MGPDH expression
To investigate thyroid cancer cell metabolism, we next analyzed OXPHOS and glycolysis rate in FTC133 and BCPAP cells using the Seahorse technique. The well-established Seahorse Mito Stress Kit analyzes key parameters of mitochondrial function by directly measuring the oxygen consumption rate (OCR) of cells utilizing sequential compound injections (Fig. 2E). We documented significantly higher baseline OCR, ATP production, maximal respiration, and spare respiratory capacity in MGPDH-rich FTC133 compared with BCPAP cells (P < 0.001; Fig. 2E).
We next used a Seahorse Glycolysis Stress Kit to compare the glycolysis rate, ECAR, among thyroid cancer cell lines utilizing sequential compound injections (Fig. 2F). No significant difference was observed between the baseline glycolysis, glycolytic capacity, and glycolytic reserve in the thyroid cancer cell lines (Fig. 2F).
These results show that MGPDH-rich FTC133 cells are characterized by a different metabolic phenotype of higher mitochondrial respiration rate than MGPDH-poor BCPAP cells, which are predominantly glycolysis-dependent. Of note, MGPDH-rich FTC133 cells are characterized by a significantly higher proliferation rate (doubling time, 27 hours) compared with BCPAP cells (doubling time, 30–36 hours; refs. 13, 20).
Metformin inhibits OXPHOS more significantly in FTC133 cells compared with BCPAP cells and switches cancer cell metabolism to glycolysis
We determined the metformin effect on mitochondrial respiration in the two-human thyroid cancer cell lines. After treatment with 5 mmol/L concentration of metformin for 48 hours (M5), we observed significant inhibition of baseline respiration, ATP production, maximum respiration, and spare respiratory capacity in FTC133 cells, as well as minimal inhibition of baseline and maximal respiration and ATP production in BCPAP cells (Fig. 3A and B). The metformin-induced OXPHOS inhibition resulted in the switching of glucose metabolism to be more anaerobic glycolysis-dependent in FTC133 and BCPAP cells, leading to a significant inhibition of glycolytic reserve in both cell lines (Fig. 3C and D).
Significantly higher reduction in mitochondrial respiration in FTC133 than in BCPAP cells after metformin therapy leading to shifting cancer cell metabolism to glycolysis. A, Significant inhibition of basal and maximal mitochondrial respiration, ATP production, and spare respiratory capacity in FTC133 cells pretreated with metformin (5 mmol/48 hours; n = 6–9). Basal OCR – FTC133 49.2 ± 1.5 pmol/min/cell versus FTC133 M5 12.05 ± 0.5 pmol/min/cell (P < 0.0001). ATP production – FTC133 40.6 ± 1.2 pmol/min/cell versus FTC133 M5 3.4 ± 0.4 pmol/min/cell (P = 0.004). Maximum respiration – FTC133 132 ± 4.7 pmol/min/cell versus FTC133 M5 46.2 ± 2.3 pmol/min/cell (P < 0.0001). Spare respiratory capacity – FTC133 82.8 ± 4.1 pmol/min/cell versus FTC133 M5 34.2 ± 1.9 pmol/min/cell (P < 0.0001). B, Significant inhibition of basal and maximal mitochondrial respiration and ATP production in BCPAP cells pretreated with metformin (5 mmol/48 hours; n = 6–9). Basal OCR – BCPAP 5.8 ± 1.1 pmol/min/cell versus BCPAP M5 1.5 ± 1.1 pmol/min/cell (P = 0.006). ATP production – BCPAP 3.5 ± 0.5 pmol/min/cell versus BCPAP M5 0.3 ± 0.1 pmol/min/cell (P = 0.0004). Maximum respiration – BCPAP 7.8 ± 0.7 pmol/min/cell versus BCPAP M5 2.2 ± 0.6 pmol/min/cell (P < 0.0001). Spare respiratory capacity – BCPAP 2 ± 0.4 pmol/min/cell versus BCPAP M5 3.1 ± 0.7 pmol/min/cell, (P = 0.22). C, Increased glucose utilization via glycolysis leading to decreased glycolytic reserve in FTC133 cells after metformin treatment (5 mmol/48 hours; n = 6–9). Glycolysis rate – FTC133 26.1 ± 12.2 mpH/min/cell versus FTC133 M5 56.2 ± 19.7 mpH/min/cell (P = 0.01.) Glycolytic capacity – FTC133 51.8 ± 22.4 mpH/min/cell versus FTC133 M5 44.4 ± 15.8 mpH/min/cell (P = 0.18). Glycolytic reserve – FTC133 25.7 ± 10.3 mpH/min/cell versus FTC133 M5 0.0 ± 0.0 mpH/min/cell (P = 0.001). D, Increased glucose utilization via glycolysis leading to decreased glycolytic capacity and glycolytic reserve in BCPAP cells after metformin treatment (5 mmol/48 hours; n = 6–9). Glycolysis rate – BCPAP 29.9 ± 9.0 mpH/min/cell versus BCPAP M5 45.6 ± 15.1 mpH/min/cell (P = 0.05). Glycolytic capacity – BCPAP 46.6 ± 17.0 mpH/min/cell versus BCPAP M5 31.1 ± 11.6 mpH/min/cell (P = 0.05). Glycolytic reserve – BCPAP 16.7 ± 10.2 mpH/min/cell versus BCPAP M5 0.0 ± 0.0 mpH/min/cell (P = 0.001). M5, metformin treatment (5 mmol/48 hours); OCR, oxygen consumption rate pmol/min/cell; ECAR, extracellular acidification rate mpH/min/cell.
Significantly higher reduction in mitochondrial respiration in FTC133 than in BCPAP cells after metformin therapy leading to shifting cancer cell metabolism to glycolysis. A, Significant inhibition of basal and maximal mitochondrial respiration, ATP production, and spare respiratory capacity in FTC133 cells pretreated with metformin (5 mmol/48 hours; n = 6–9). Basal OCR – FTC133 49.2 ± 1.5 pmol/min/cell versus FTC133 M5 12.05 ± 0.5 pmol/min/cell (P < 0.0001). ATP production – FTC133 40.6 ± 1.2 pmol/min/cell versus FTC133 M5 3.4 ± 0.4 pmol/min/cell (P = 0.004). Maximum respiration – FTC133 132 ± 4.7 pmol/min/cell versus FTC133 M5 46.2 ± 2.3 pmol/min/cell (P < 0.0001). Spare respiratory capacity – FTC133 82.8 ± 4.1 pmol/min/cell versus FTC133 M5 34.2 ± 1.9 pmol/min/cell (P < 0.0001). B, Significant inhibition of basal and maximal mitochondrial respiration and ATP production in BCPAP cells pretreated with metformin (5 mmol/48 hours; n = 6–9). Basal OCR – BCPAP 5.8 ± 1.1 pmol/min/cell versus BCPAP M5 1.5 ± 1.1 pmol/min/cell (P = 0.006). ATP production – BCPAP 3.5 ± 0.5 pmol/min/cell versus BCPAP M5 0.3 ± 0.1 pmol/min/cell (P = 0.0004). Maximum respiration – BCPAP 7.8 ± 0.7 pmol/min/cell versus BCPAP M5 2.2 ± 0.6 pmol/min/cell (P < 0.0001). Spare respiratory capacity – BCPAP 2 ± 0.4 pmol/min/cell versus BCPAP M5 3.1 ± 0.7 pmol/min/cell, (P = 0.22). C, Increased glucose utilization via glycolysis leading to decreased glycolytic reserve in FTC133 cells after metformin treatment (5 mmol/48 hours; n = 6–9). Glycolysis rate – FTC133 26.1 ± 12.2 mpH/min/cell versus FTC133 M5 56.2 ± 19.7 mpH/min/cell (P = 0.01.) Glycolytic capacity – FTC133 51.8 ± 22.4 mpH/min/cell versus FTC133 M5 44.4 ± 15.8 mpH/min/cell (P = 0.18). Glycolytic reserve – FTC133 25.7 ± 10.3 mpH/min/cell versus FTC133 M5 0.0 ± 0.0 mpH/min/cell (P = 0.001). D, Increased glucose utilization via glycolysis leading to decreased glycolytic capacity and glycolytic reserve in BCPAP cells after metformin treatment (5 mmol/48 hours; n = 6–9). Glycolysis rate – BCPAP 29.9 ± 9.0 mpH/min/cell versus BCPAP M5 45.6 ± 15.1 mpH/min/cell (P = 0.05). Glycolytic capacity – BCPAP 46.6 ± 17.0 mpH/min/cell versus BCPAP M5 31.1 ± 11.6 mpH/min/cell (P = 0.05). Glycolytic reserve – BCPAP 16.7 ± 10.2 mpH/min/cell versus BCPAP M5 0.0 ± 0.0 mpH/min/cell (P = 0.001). M5, metformin treatment (5 mmol/48 hours); OCR, oxygen consumption rate pmol/min/cell; ECAR, extracellular acidification rate mpH/min/cell.
We also analyzed the acute effects of metformin on thyroid cancer metabolism. We performed Seahorse analysis after 0.5, 1, 2, and 4 hours of exposure to metformin. We found a significant decrease in baseline and maximal OCR and ATP production in metformin-treated cells compared with control cells at 2 and 4 hours (Supplementary Fig. S2C and S2D). However, there was no difference in MGPDH expression after acute exposure to metformin at these time points (Supplementary Fig. S2E). Metformin inhibited the enzymatic activity of MGPDH as early as 10 minutes after exposure (Fig. 2D). These data suggest that the acute effects of metformin on cancer cells’ OXPHOS are due to the inhibition of the enzymatic activity of respiratory chain enzymes, targeting MGPDH enzymatic activity.
In summary, we demonstrate that FTC133 cells, characterized by higher MGPDH expression and mitochondrial respiration rate, are more sensitive to OXPHOS-inhibitory effects of metformin in vitro compared with BCPAP cells. Metformin therapy results in switching cancer cell metabolism toward glycolysis.
Metformin-induced inhibition of growth and mitochondrial respiration is dependent on MGPDH in thyroid cancer cells
To determine the direct effect of MGPDH on growth and metabolism of thyroid cancer cells, we utilized two methods of MGPDH silencing: transient silencing using siRNA and a knockdown model using CRISPR/Cas9 gene editing. Silencing efficacy was documented by qRT-PCR and Western blot analysis in the siRNA model (Fig. 4A and B) and by Western blot analysis in the knockdown model (Fig. 4B).
MGPDH silencing leads to decreased cancer cell growth and mitochondrial respiration, and causes resistance to growth and OXPHOS inhibition by metformin. A, Decreased mRNA expression of MGPDH after siRNA silencing in FTC133 (P < 0.0001) and BCPAP (P = 0.007) cells (n = 2–5). B, Decreased protein expression of MGPDH after siRNA silencing in FTC133 and BCPAP cells. Decreased protein expression of MGPDH after monoallelic knockdown in FTC133 and BCPAP cells. V, empty vector; KD, MGPDH knockdown. C, Decreased cell growth after silencing of MGPDH in FTC133 and BCPAP cells, (n = 18, P < 0.0001). D, Decreased cell growth after MGPDH knockdown in FTC133 (P < 0.0001) and BCPAP (P = 0.0001) cells. Reduced growth-inhibitory metformin action in MGPDH-knockdown cells compared with control cells (n = 10–11). E, Decreased baseline and maximal mitochondrial respiration, ATP production, and spare respiratory capacity in FTC133 cells after MGPDH knockdown. Reduced mitochondrial respiration inhibitory effects of metformin in knockdown cells compared with control cells, (n = 6–9). Basal OCR vector (V) versus knockdown (KD) – FTC133 V 106 ± 7 pmol/min/cell versus FTC133 KD 14.8 ± 2.7 pmol/min/cell, P < 0.0001. ATP production – FTC133 V 84.8 ± 4.4 pmol/min/cell versus FTC133 KD 9.02 ± 1.9 pmol/min/cell, P < 0.0001. Maximum respiration – FTC133 V 164.8 ± 14.8 pmol/min/cell versus FTC133 KD 27.1 ± 2.8 pmol/min/cell, P < 0.0001. Spare respiratory capacity – FTC133 V 57.9 ± 7.6 pmol/min/cell versus FTC133 KD 12.3 ± 0.52 pmol/min/cell, P < 0.0001. F, Decreased baseline and maximal mitochondrial respiration, ATP production, and proton leak in BCPAP cells after MGPDH knockdown. Reduced mitochondrial respiration–inhibitory effects of metformin in knockdown cells compared with control cells, (n = 6–9). Basal OCR vector (V) versus knockdown (KD)–BCPAP V 12.1 ± 2 pmol/min/cell versus BCPAP KD 4.05 ± 0.51 pmol/min/cell, P = 0.0001. ATP production – BCPAP V 9 ± 1.6 pmol/min/cell versus BCPAP KD 2.8 ± 0.3 pmol/min/cell, P = 0.0002. Maximum respiration – BCPAP V 9.4 ± 2.7 pmol/min/cell versus BCPAP KD 6.7 ± 0.3 pmol/min/cell, P = 0.01. Spare respiratory capacity – BCPAP V 0 pmol/min/cell versus BCPAP KD 2.7 ± 0.2 pmol/min/cell, P < 0.0001. NC, control cells; si, cell transfected with MGPDH siRNA; V, control cells transfected with pCAS empty vector; KD, cells with monoallelic MGPDH knockdown; M5, metformin 5 mmol/48 hours.
MGPDH silencing leads to decreased cancer cell growth and mitochondrial respiration, and causes resistance to growth and OXPHOS inhibition by metformin. A, Decreased mRNA expression of MGPDH after siRNA silencing in FTC133 (P < 0.0001) and BCPAP (P = 0.007) cells (n = 2–5). B, Decreased protein expression of MGPDH after siRNA silencing in FTC133 and BCPAP cells. Decreased protein expression of MGPDH after monoallelic knockdown in FTC133 and BCPAP cells. V, empty vector; KD, MGPDH knockdown. C, Decreased cell growth after silencing of MGPDH in FTC133 and BCPAP cells, (n = 18, P < 0.0001). D, Decreased cell growth after MGPDH knockdown in FTC133 (P < 0.0001) and BCPAP (P = 0.0001) cells. Reduced growth-inhibitory metformin action in MGPDH-knockdown cells compared with control cells (n = 10–11). E, Decreased baseline and maximal mitochondrial respiration, ATP production, and spare respiratory capacity in FTC133 cells after MGPDH knockdown. Reduced mitochondrial respiration inhibitory effects of metformin in knockdown cells compared with control cells, (n = 6–9). Basal OCR vector (V) versus knockdown (KD) – FTC133 V 106 ± 7 pmol/min/cell versus FTC133 KD 14.8 ± 2.7 pmol/min/cell, P < 0.0001. ATP production – FTC133 V 84.8 ± 4.4 pmol/min/cell versus FTC133 KD 9.02 ± 1.9 pmol/min/cell, P < 0.0001. Maximum respiration – FTC133 V 164.8 ± 14.8 pmol/min/cell versus FTC133 KD 27.1 ± 2.8 pmol/min/cell, P < 0.0001. Spare respiratory capacity – FTC133 V 57.9 ± 7.6 pmol/min/cell versus FTC133 KD 12.3 ± 0.52 pmol/min/cell, P < 0.0001. F, Decreased baseline and maximal mitochondrial respiration, ATP production, and proton leak in BCPAP cells after MGPDH knockdown. Reduced mitochondrial respiration–inhibitory effects of metformin in knockdown cells compared with control cells, (n = 6–9). Basal OCR vector (V) versus knockdown (KD)–BCPAP V 12.1 ± 2 pmol/min/cell versus BCPAP KD 4.05 ± 0.51 pmol/min/cell, P = 0.0001. ATP production – BCPAP V 9 ± 1.6 pmol/min/cell versus BCPAP KD 2.8 ± 0.3 pmol/min/cell, P = 0.0002. Maximum respiration – BCPAP V 9.4 ± 2.7 pmol/min/cell versus BCPAP KD 6.7 ± 0.3 pmol/min/cell, P = 0.01. Spare respiratory capacity – BCPAP V 0 pmol/min/cell versus BCPAP KD 2.7 ± 0.2 pmol/min/cell, P < 0.0001. NC, control cells; si, cell transfected with MGPDH siRNA; V, control cells transfected with pCAS empty vector; KD, cells with monoallelic MGPDH knockdown; M5, metformin 5 mmol/48 hours.
MGPDH silencing by siRNA resulted in cell proliferation inhibition by 33.6% ± 5% in FTC133 cells and 30.7% ± 6.7% in BCPAP cells (Fig. 4C), while monoallelic knockdown resulted in inhibition of proliferation by 38.6% ± 2.6% in FTC133 cells and 25.2% ± 5.2% in BCPAP cells (Fig. 4D). MGPDH knockdown was associated with mitochondrial respiration inhibition in FTC133 and BCPAP cells (Fig. 4E and F), with a more pronounced effect in FTC133 cells. Similar effects on mitochondrial respiration were observed in both cell lines in a transient silencing model (Supplementary Fig. S3A and S3B).
We next examined the effect of metformin with MGPDH-silencing conditions. While metformin inhibited growth of control cells, its growth-inhibitory effect in MGPDH knockdown model was significantly lowered with growth inhibition by 19% ± 3.1% in FTC133 cells and 17.1% ± 5.1% in BCPAP cells (Fig. 4D). Metformin significantly inhibited basal OCR, ATP production, maximum respiration, and spare respiratory capacity in control cells, while the net OCR-inhibitory effect of metformin in knockdown cells was significantly reduced in both cell lines (Fig. 4E and F; Supplementary Table S1). The quantification of a net difference (delta) in the OCR-inhibitory effect of metformin in control cells and MGPDH-knockdown cells is depicted in the Supplementary Fig. S3C and S3D.
Taken together, we showed that MGPDH silencing leads to inhibition of growth and mitochondrial respiration in thyroid cancer cells and decreased efficacy of the antiproliferative and OXPHOS-inhibitory effects of metformin in vitro.
MGPDH overexpression leads to increased thyroid cancer cell growth and mitochondrial respiration and sensitizes thyroid cancer cells to metformin
MGPDH transcription is activated by TR/RXR binding to a TRE (thyroid hormone response element) sequence in the promoter region of the MGPDH gene (21). We found that treatment with thyroid hormone (T3) resulted in significant MGPDH upregulation in both thyroid cancer cell lines, with the maximum effect observed at 100 nmol/L concentration for 72 hours (Fig. 5A and B). We also successfully induced MGPDH overexpression in thyroid cancer cell lines by an MGPDH-containing plasmid (Fig. 5C and D).
MGPDH overexpression results in increased cancer cell growth and mitochondrial respiration and sensitizes cancer cells to OCR-inhibitory effects of metformin. A, 100 nmol/L/72-hour triiodothyronine (T3) treatment results in increased MGPDH expression at the mRNA level in FTC133 (P < 0.0001) and BCPAP (P = 0.002) cells (n = 4–8). B, 100 nmol/L/72-hour T3 treatment results in increased MGPDH expression at the protein level in FTC133 and BCPAP cells. C, Plasmid induced overexpression of MGPDH at the mRNA level in FTC133 and BCPAP cells (n = 3–4, P = 0.0005). D, Plasmid induced overexpression of MGPDH at the protein level in FTC133 and BCPAP cells. E, Metformin treatment (5 mmol/48 hours) is associated with decreased expression of MGPDH mRNA in control FTC133 (n = 4, P = 0.001) and MGPDH-overexpressing FTC133 cells, (n = 4; P = 0.03). F, MGPDH overexpression is associated with increased FTC133 thyroid cancer cell growth and sensitizes FTC133 cells to metformin growth–inhibitory effects (n = 10–11; P = 0.003). G, Significant increase in basal and maximum respiration and ATP production in MGPDH-overexpressing FTC133 compared with control cells and increased sensitivity to OXPHOS-inhibitory effects of metformin (5 mmol/48 hours; n = 6–9; Supplementary Table S2). H, Metformin treatment (5 mmol/48 hours) is associated with decreased expression of MGPDH mRNA in control BCPAP cells (n = 4; P = 0.002), but not in plasmid induced MGPDH-overexpressing BCPAP cells (n = 4; P = 0.27). I, MGPDH-overexpression is associated with increased BCPAP cell growth and sensitizes BCPAP cells to metformin growth-inhibitory effects (n = 10–11, P < 0.0001). J, Significant increase in basal respiration and ATP production in MGPDH-overexpressing BCPAP compared with control cells and increased sensitivity to OXPHOS-inhibitory effects of metformin (5 mmol/48 hours). The basal respiration of MGPDH-overexpressing BCPAP cells operates on maximum capacity, as adding the uncoupler–FCCP–does not stimulate mitochondrial respiration any further, (n = 6–9), (Supplementary Table S2). T3, triiodothyronine 100 nmol/L/72 hours; V, control cells transfected with empty vector; OE, cells with MGPDH-overexpression; M5, metformin 5 mmol/48 hours.
MGPDH overexpression results in increased cancer cell growth and mitochondrial respiration and sensitizes cancer cells to OCR-inhibitory effects of metformin. A, 100 nmol/L/72-hour triiodothyronine (T3) treatment results in increased MGPDH expression at the mRNA level in FTC133 (P < 0.0001) and BCPAP (P = 0.002) cells (n = 4–8). B, 100 nmol/L/72-hour T3 treatment results in increased MGPDH expression at the protein level in FTC133 and BCPAP cells. C, Plasmid induced overexpression of MGPDH at the mRNA level in FTC133 and BCPAP cells (n = 3–4, P = 0.0005). D, Plasmid induced overexpression of MGPDH at the protein level in FTC133 and BCPAP cells. E, Metformin treatment (5 mmol/48 hours) is associated with decreased expression of MGPDH mRNA in control FTC133 (n = 4, P = 0.001) and MGPDH-overexpressing FTC133 cells, (n = 4; P = 0.03). F, MGPDH overexpression is associated with increased FTC133 thyroid cancer cell growth and sensitizes FTC133 cells to metformin growth–inhibitory effects (n = 10–11; P = 0.003). G, Significant increase in basal and maximum respiration and ATP production in MGPDH-overexpressing FTC133 compared with control cells and increased sensitivity to OXPHOS-inhibitory effects of metformin (5 mmol/48 hours; n = 6–9; Supplementary Table S2). H, Metformin treatment (5 mmol/48 hours) is associated with decreased expression of MGPDH mRNA in control BCPAP cells (n = 4; P = 0.002), but not in plasmid induced MGPDH-overexpressing BCPAP cells (n = 4; P = 0.27). I, MGPDH-overexpression is associated with increased BCPAP cell growth and sensitizes BCPAP cells to metformin growth-inhibitory effects (n = 10–11, P < 0.0001). J, Significant increase in basal respiration and ATP production in MGPDH-overexpressing BCPAP compared with control cells and increased sensitivity to OXPHOS-inhibitory effects of metformin (5 mmol/48 hours). The basal respiration of MGPDH-overexpressing BCPAP cells operates on maximum capacity, as adding the uncoupler–FCCP–does not stimulate mitochondrial respiration any further, (n = 6–9), (Supplementary Table S2). T3, triiodothyronine 100 nmol/L/72 hours; V, control cells transfected with empty vector; OE, cells with MGPDH-overexpression; M5, metformin 5 mmol/48 hours.
MGPDH overexpression was associated with increased thyroid cancer cell proliferation by 13% ± 3.5% in FTC133 cells and 37.4% ± 5.8% in BCPAP cells (Fig. 5F and I). MGPDH overexpression leads to significant activation of baseline and maximal mitochondrial respiration as well as ATP production in FTC133 and BCPAP cells, with the net effect larger in BCPAP cells (Supplementary Fig. S4A and S4B). Of note, BCPAP cells with MGPDH overexpression were operating on maximal respiratory rate at baseline because FCCP addition did not further stimulate maximal respiration (Supplementary Fig. S4B). Similar effects of increased cancer cell mitochondrial respiration were observed after induction of MGPDH expression by T3 treatment (100 nmol/L/72 hours; Supplementary Fig. S4C and S4D).
In plasmid-induced MGPDH-overexpressing thyroid cancer cell lines, metformin treatment inhibited MGPDH expression only in FTC133 cells characterized by high endogenous MGPDH level, but not in BCPAP cells (Fig. 5E and H). Metformin inhibited cell proliferation by 42.2% ± 3.5% in MGPDH-overexpressing FTC133 cells and by 36.5% ± 3.5% in control cells (Fig. 5F). The growth-inhibitory effect of metformin in mGPDH-overexpressing cells was higher in the BCPAP model, with growth inhibition of 49.6% ± 5.5% in MGPDH-overexpressing cells compared with 32.3% ± 5.5% in control cells (Fig. 5I). FTC133 and BCPAP cells with stable MGPDH overexpression were characterized by higher metformin-induced inhibition of basal and maximal respiration as well as ATP production compared with control cells (Fig. 5G and J; Supplementary Table S2). The quantification of a net difference (delta) in the OCR-inhibitory effect of metformin in control and MGPDH-overexpressing cells is depicted in Supplementary Fig. S4E and S4F.
Similar effects were observed in FTC133 and BCPAP cells with T3-induced MGPDH overexpression (Supplementary Fig. S5A and S5B; Supplementary Table S3). The net difference in mitochondrial respiration–inhibitory effects of metformin in control and T3-treated cells are summarized in Supplementary Fig. S5C and S5D.
Taken together, these data showed that overexpression of MGPDH in thyroid cancer cells shifts their metabolism to increased mitochondrial respiration and stimulates thyroid cancer cell proliferation. Therefore, thyroid cancer cells characterized by MGPDH overexpression are more sensitive to metformin-induced inhibition of growth and OXPHOS.
Metformin treatment is associated with cancer growth inhibition in an FTC133 metastatic mouse model, but not in a BCPAP metastatic mouse model
We established a metastatic mouse model by tail vein injection of luciferase-labeled FTC133 and BCPAP thyroid cancer cells into immunocompromised mice (17). This model results in the development of lungs metastases in 100% of the animals, and liver and bone metastases in a subset of mice. We began metformin therapy in these models once IVIS (In Vivo Imaging System) confirmed the development of lung metastases (2 weeks postinjection; Fig. 6C and D). Metformin treatment led to a significant reduction in thyroid cancer metastatic growth in the FTC133 model, but not in the BCPAP model (Fig. 6A, B, E, and G). A significantly lower proportion of metformin-treated FTC133 mice developed liver metastases, while there was no difference in the rate of liver metastases in the BCPAP model (Fig. 6F and H). Most importantly, the growth-inhibitory effects of metformin in vivo in the FTC133 model were achieved utilizing a metformin concentration reflecting a therapeutic concentration in humans (7.9 ± 3.9 μmol/L). Metformin-treated animals were characterized by significantly lower serum glucose levels than placebo-treated animals (FTC133 model – MF glucose 132.2 ± 44.9 mg/dL vs. placebo glucose 207.3 ± 74.8 mg/dL, P < 0.001; BCPAP model – MF glucose 126.3 ± 67.9 mg/dL vs. placebo glucose 171.5 ± 85.6 mg/dL, P = 0.049). There was no difference in the insulin level in control and metformin-treated animals (FTC133 model – MF insulin 0.266 ± 0.21 ng/mL vs. placebo insulin 0.399 ± 0.3 ng/mL, P = 0.1; BCPAP model – MF insulin 0.53 ± 0.49 ng/mL vs. placebo insulin 0.57 ± 0.31 ng/mL, P = 0.7).
Metformin inhibits growth of thyroid cancer metastases in FTC133 metastatic mouse model, but not in BCPAP metastatic mouse model. A, Significantly decreased growth rate of metastases in metformin-treated mice (n = 14) compared with placebo-treated animals (n = 13) documented by weekly bioluminescence measurement in FTC133 metastatic mouse model (P = 0.046). B, No difference in the growth rate of metastases in metformin-treated mice (n = 12) compared with placebo-treated animals (n = 13) documented by weekly bioluminescence measurement in BCPAP metastatic mouse model (P = 0.6). C, The representative images of the metastatic disease burden at the baseline and the end of the in vivo experiment in control and metformin-treated animals in metastatic FTC133 mouse model. D, The representative images of the metastatic disease burden at the baseline and the end of the in vivo experiment in control and metformin-treated animals in metastatic BCPAP mouse model. E, Representative hematoxylin and eosin (H&E)-stained slides of lung and liver metastases derived from FTC133 metastatic mouse model in control, placebo-treated, and metformin-treated animals. F, Significantly lower rate of liver metastases in metformin-treated mice compared with placebo-treated animals in FTC133 metastatic mouse model (P = 0.005). G, Representative H&E-stained slide of the lung metastases derived from BCPAP metastatic mouse model. H, No difference in the liver metastases rate in metformin-treated mice compared with placebo (P = 0.8). C, control placebo-treated animals; MF, metformin-treated animals.
Metformin inhibits growth of thyroid cancer metastases in FTC133 metastatic mouse model, but not in BCPAP metastatic mouse model. A, Significantly decreased growth rate of metastases in metformin-treated mice (n = 14) compared with placebo-treated animals (n = 13) documented by weekly bioluminescence measurement in FTC133 metastatic mouse model (P = 0.046). B, No difference in the growth rate of metastases in metformin-treated mice (n = 12) compared with placebo-treated animals (n = 13) documented by weekly bioluminescence measurement in BCPAP metastatic mouse model (P = 0.6). C, The representative images of the metastatic disease burden at the baseline and the end of the in vivo experiment in control and metformin-treated animals in metastatic FTC133 mouse model. D, The representative images of the metastatic disease burden at the baseline and the end of the in vivo experiment in control and metformin-treated animals in metastatic BCPAP mouse model. E, Representative hematoxylin and eosin (H&E)-stained slides of lung and liver metastases derived from FTC133 metastatic mouse model in control, placebo-treated, and metformin-treated animals. F, Significantly lower rate of liver metastases in metformin-treated mice compared with placebo-treated animals in FTC133 metastatic mouse model (P = 0.005). G, Representative H&E-stained slide of the lung metastases derived from BCPAP metastatic mouse model. H, No difference in the liver metastases rate in metformin-treated mice compared with placebo (P = 0.8). C, control placebo-treated animals; MF, metformin-treated animals.
We analyzed MGPDH expression in the lung and liver metastases derived from the metastatic mouse models. In the FTC133 model, we found significantly lower MGPDH expression in tumors derived from metformin-treated animals compared with robust expression of MGPDH in the placebo group (P = 0.03; Fig. 7A). The MGPDH expression was very low (score 0 and 1) in metastases derived from the BCPAP model, without a significant difference in control and metformin-treated groups (P = 0.51; Fig. 7B).
Significantly lower expression of MGPDH in metastatic lesions derived from FTC133 metastatic mouse model and thyroid tumor tissue derived from transgenic TRβPV/PVPten ± mouse model with documented cancer growth- and progression-inhibitory effects of metformin. A, FTC133-metastatic mouse model—significantly lower expression of MGPDH in lung and liver metastases derived from metformin-treated (n = 14) compared with placebo-treated mice (n = 13; P = 0.03). B, BCPAP-metastatic mouse model—very low and comparable expression of MGPDH in lung metastases derived from metformin-treated (n = 12) and placebo-treated mice (n = 13; P = 0.51). C, Transgenic TRβPV/PVPten ± mouse model – significantly lower expression of MGPDH in thyroid tumor tissue derived from metformin-treated (n = 8) compared with placebo-treated mice (n = 8; P = 0.001).
Significantly lower expression of MGPDH in metastatic lesions derived from FTC133 metastatic mouse model and thyroid tumor tissue derived from transgenic TRβPV/PVPten ± mouse model with documented cancer growth- and progression-inhibitory effects of metformin. A, FTC133-metastatic mouse model—significantly lower expression of MGPDH in lung and liver metastases derived from metformin-treated (n = 14) compared with placebo-treated mice (n = 13; P = 0.03). B, BCPAP-metastatic mouse model—very low and comparable expression of MGPDH in lung metastases derived from metformin-treated (n = 12) and placebo-treated mice (n = 13; P = 0.51). C, Transgenic TRβPV/PVPten ± mouse model – significantly lower expression of MGPDH in thyroid tumor tissue derived from metformin-treated (n = 8) compared with placebo-treated mice (n = 8; P = 0.001).
We next examined MGPDH expression in a transgenic TRβPV/PVPten ± mouse model, which spontaneously develops FTC, characterized by increased tumor aggressiveness when subjected to high-fat diet–induced obesity. Metformin was found to block the progression of thyroid cancer in this model (22). We obtained thyroid cancer tissue samples from TRβPV/PVPten ± mice derived from control and metformin-treated animals. We found a significantly higher expression of MGPDH in the control group compared with the metformin-treated animals (P = 0.001; Fig. 7C).
In summary, we demonstrated the inhibitory effects of metformin on both tumor growth and progression in two thyroid cancer in vivo models characterized by high tissue MGPDH expression, and lack of growth-inhibitory effects of metformin in a metastatic mouse model characterized by low tissue MGPDH expression. Therefore, our in vitro and in vivo data suggest that MGPDH is a target of metformin in thyroid cancer and may serve as a potential biomarker of the response to metformin treatment (Supplementary Fig. S6).
Discussion
We demonstrate for the first time that MGPDH is overexpressed in thyroid cancer. We show that MGPDH overexpression leads to increased thyroid cancer cell growth and stimulation of mitochondrial respiration to meet the metabolic demands of increased proliferation. Conversely, we demonstrate MGPDH knockdown is associated with growth inhibition and reduction of OXPHOS. Our data show that MGPDH is a metformin target in thyroid cancer. Metformin treatment is associated with decreased MGPDH expression, inhibition of growth and mitochondrial respiration, as well as the shifting of cancer cell metabolism toward glycolysis. MGPDH silencing decreases sensitivity to the inhibition of OXPHOS and growth by metformin, while MGPDH overexpression sensitizes thyroid cancer cells to mitochondrial respiration and the growth-inhibitory effects of metformin. We also show that in vivo metastatic and transgenic thyroid cancer models characterized by high tissue MGPDH expression are responsive to cancer progression–inhibitory effects of therapeutic concentrations of metformin, while metformin does not inhibit the growth of metastatic lesions in a metastatic thyroid cancer model characterized by low tumor tissue MGPDH expression.
To study the direct effects of metformin on potential intracellular targets in cancer cells, it is necessary to prove that metformin intracellular transporters are expressed. We found that thyroid cancer is characterized by a significant overexpression of OCT1 transporter. OCT1 has also been found to be overexpressed in other malignancies such as leukemia, prostate cancer, hepatocellular carcinoma, colorectal cancer, or renal cell carcinoma (23–27). OCT1 downregulation limits the antineoplastic effects of metformin, further supporting the significant role of OCT1 expression in cancer response to metformin (28). The same metformin transporter plays a major role in metformin action in the liver (8).
MGPDH has been found to be a major metformin target in the liver, acting by shifting the NADH/NAD+ ratio via inhibition of glycerol-3-phosphate shuttle, directly blocking gluconeogenesis from glycerol and preventing clearance of cytosolic NADH, impairing glucose production from lactate (9). Interestingly, an MGPDH knockdown in vivo model recapitulated the effects of metformin treatment, and metformin had no additional gluconeogenesis-inhibitory function in these animals (9). Of note, this landmark study documented metformin-induced noncompetitive inhibition of the enzymatic activity of MGPDH, but did not investigate the role of metformin in regulating MGPDH expression (9). One of the reasons for noncompetitive inhibition is decreased expression of the enzyme. We now demonstrate that metformin inhibits MGPDH activity and expression at the mRNA and protein level. The proposed mechanism of inhibition of MGPDH enzymatic activity is via conformational change at the FAD/FADH molecule of the MGPDH enzyme complex as suggested by Madiraju and colleagues (9).
It was not known whether MGPDH could act as a molecular target of metformin in cancer. In fact, several other mechanisms of metformin action in cancer have been identified, including activation of the intracellular energy sensor AMPK and LKB1-AMPK–induced inhibition of Raptor–mTOR (mTORC1) complex, leading to suppression of cellular protein synthesis and growth. Several studies suggest that the major mechanism of induction of energetic stress in cancer cells leading to activation of AMPK is metformin-induced inhibition of mitochondrial complex I (29, 30). AMPK-independent mechanisms of metformin action have also been described, such as the inhibitory effect on AKT phosphorylation in breast and lung cancers, suppression of HER2 oncoprotein, suppression of NFκB and CREB signaling, direct inhibition of hexokinase 2 and indirect inhibition of IGFR signaling due to decreased plasma concentrations of insulin (31, 32).
Our study, demonstrating the role of MGPDH as a metformin target in thyroid cancer cells, facilitates our understanding of a cross-talk between major signaling pathways and cancer metabolism. Acknowledging the role of MGPDH in the glycerol-3-phosphate shuttle connecting the major cancer metabolic pathways, glycolysis with OXPHOS, underscores its potentially important role in shifting the metabolic phenotype of cancer cells. There have been major discoveries in cancer metabolism recently, including identification of the “reverse Warburg effect”—the cross-talk between glycolysis-dependent cancer cells and a population of cancer cells characterized by a high level of OXPHOS (33, 34). Cancer cells behave as metabolic parasites with glycolysis-dependent cells shifting the glycolysis endproduct l-Lactate to another subpopulation of cancer cells characterized by active mitochondrial respiration. l-Lactate is then utilized to fuel OXPHOS-based energy production required for cancer growth (33, 34). Birsoy and colleagues documented a significant role of OXPHOS as a metabolic pathway warranting optimal cancer cell proliferation specifically in low glucose conditions (35). The authors postulated that cell lines characterized by impaired glucose utilization or mtDNA mutation in complex I genes are the most sensitive to biguanide action, such as phenformin, in in vitro and in vivo xenograft models. We found that the cells characterized by increased OXPHOS were more vulnerable to growth and OXPHOS-inhibitory effects of metformin. These data are consistent with the observation by Javeshghani and colleagues, who showed that the cells that are forced to rely on OXPHOS are more affected by the actions of metformin (36). We demonstrate that one of the mechanisms of sensitivity to metformin is modulated via MGPDH expression and enzymatic activity. Therefore, targeting MGPDH as an important factor that is involved in the regulation of OXPHOS and shifting cancer metabolism toward glycolysis might be an effective strategy against cancer.
Our study revealed significant MGPDH overexpression in a large number of human thyroid cancer tissue samples. There are very limited published data on MGPDH expression in cancer. Studies from the early 70s have documented significantly higher expression of MGPDH in Novikoff hepatoma compared with the normal liver (37), while more recent studies document MGPDH overexpression in prostate cancer and prostate cancer cell lines (38, 39).
With MGPDH expression regulated by thyroid hormones, predominantly T3, it is worthwhile to speculate that our observation of significant MGPDH overexpression in thyroid cancer tissue samples might be partially due to the fact that the majority of intermediate and high-risk patients with thyroid cancer are treated with the supraphysiologic doses of thyroid hormones to achieve TSH suppression. Unfortunately, we did not have data on the thyroid hormone levels available to correlate MGPDH tissue expression with free T4 and T3 serum levels in a cohort of patients whose tissue samples were analyzed. We demonstrate that T3 treatment results in increased MGPDH expression in thyroid cancer cells and leads to OXPHOS activation in thyroid cancer cell lines, particularly those characterized by low baseline MGPDH and OXPHOS levels. A similar observation documenting T3-induced MGPDH overexpression in Moris hepatoma, characterized by low baseline MGPDH level, but not in Novikoff hepatoma, characterized by baseline high MGPDH expression, was reported in the early 1970s by Hunt and colleagues (37).
Our in vivo models further support the role of MGPDH as a metformin target as evidenced by cancer growth and progression-inhibitory effects of metformin in metastatic and transgenic mouse models characterized by high tumor tissue MGPDH expression, and lack thereof in a metastatic mouse model characterized by low MGPDH expression. Moreover, metformin treatment was associated with a significant reduction in tumor tissue MGPDH expression in MGPDH-rich metastatic and transgenic mouse models. Metformin-treated mice were characterized by significantly lower serum glucose levels than control animals. Our data suggest that glycemia did not play a major role in modulating cancer growth-inhibitory effects of metformin, as lower metformin-induced glucose levels were observed in both metastatic mouse models, one characterized by a significant inhibition of cancer progression (FTC133) and the other without the growth-inhibitory effects of metformin (BCPAP). Shen and colleagues documented that metformin inhibits glucose uptake by cancer cells as shown by 18FDG-PET imaging in PTC subcutaneous xenograft models (KTC1 and BCPAP cells; ref. 40). In this study, the authors did not provide information on cancer growth in vivo nor information about the serum glucose levels, which could have significantly affected 18FDG-PET imaging (41).
Several studies have suggested indirect cancer progression–inhibitory effects of metformin in vivo via reduction of serum insulin levels leading to decreased IGFR signaling (31, 32). Our study suggests that this mechanism of metformin action may not play a key role in thyroid cancer metastatic mouse models because, although we observed lower insulin levels in metformin-treated animals compared with placebo, the difference was not statistically significant. To the best of our knowledge, there are no other studies testing the role of metformin in thyroid cancer metastatic mouse models.
In conclusion, based on in vitro and in vivo functional studies, we demonstrated that MGPDH is overexpressed in thyroid cancer and serves as a metformin target. Metformin-induced inhibition of MGPDH leads to a shift in thyroid cancer metabolism with significant inhibition of OXPHOS and ATP production. High tumor tissue MGPDH expression may predict susceptibility to the growth-inhibitory effects of metformin in vivo. Our observations may have a significant impact for an individualized approach to cancer therapy. Clinical trials in humans will be necessary to confirm the role of MGPDH as a predictor of response to metformin therapy.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: S. Thakur, E. Kebebew, J. Klubo-Gwiezdzinska
Development of methodology: S. Thakur, J. Klubo-Gwiezdzinska
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): K. Gaskins, V.V. Vasko, M. Boufraqech, D. Patel, C. Sourbier, J.M. Reece, E. Kebebew, J. Klubo-Gwiezdzinska
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): S. Thakur, B. Daley, V.V. Vasko, D. Patel, C. Sourbier, J.M. Reece, S.-Y. Cheng, E. Kebebew, S.K. Agarwal, J. Klubo-Gwiezdzinska
Writing, review, and/or revision of the manuscript: S. Thakur, B. Daley, M. Boufraqech, D. Patel, S.-Y. Cheng, E. Kebebew, S.K. Agarwal, J. Klubo-Gwiezdzinska
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): S. Thakur, B. Daley, K. Gaskins, S.-Y. Cheng, E. Kebebew
Study supervision: E. Kebebew, J. Klubo-Gwiezdzinska
Acknowledgments
We thank all co-workers who supported this work at our institution, with special thanks to Oksana Gavrilova—NIDDK Animal Core Facility, Kristopher Krausz – NCI Mass Spectrometry Core Facility, Matthew Breymaier—assistance with TCGA database, Jeongwon Park and Sunmi Park—preparation of tissue samples derived from the TRβPV/PVPten ± mouse model, and Lee Weinstein—critical analysis of the manuscript. This work was supported by an NIH intramural grant DK075140.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
References
Supplementary data
TCGA normal vs tumor raw data
STR PROFILE
IMPACT II