ORIGINAL RESEARCH article

Front. Pharmacol., 03 December 2021

Sec. Pharmacology of Anti-Cancer Drugs

Volume 12 - 2021 | https://doi.org/10.3389/fphar.2021.670167

Antiproliferative, Antiangiogenic, and Antimetastatic Therapy Response by Mangiferin in a Syngeneic Immunocompetent Colorectal Cancer Mouse Model Involves Changes in Mitochondrial Energy Metabolism

  • 1. Centro de Investigación y Desarrollo de Medicamentos (CIDEM), La Habana, Cuba

  • 2. Laboratorio de Farmacología, Instituto de Ciencias del Mar (ICIMAR), CITMA, La Habana, Cuba

  • 3. Laboratory of Protein Science, Proteomics and Epigenetic Signaling (PPES) and Integrated Personalized and Precision Oncology Network (IPPON), Department of Biomedical Sciences, University of Antwerp, Campus Drie Eiken, Antwerp, Belgium

  • 4. Department of Food Chemistry, Technology and Biotechnology, Faculty of Chemistry, Gdansk University of Technology, Gdansk, Poland

  • 5. Departamento de Farmacología, Instituto de Ciencias Biológicas (ICB), Universidad Federal de Minas Gerais (UFMG), Horizonte, Brazil

  • 6. Center of Medical Genetics, University of Antwerp and Antwerp University Hospital, Antwerp, Belgium

  • 7. Centro de Estudios para las Investigaciones y Evaluaciones Biológicas (CEIEB), Instituto de Farmacia y Alimentos (IFAL), Universidad de La Habana, La Habana, Cuba

  • 8. Facultad de Ciencias Naturales y Agropecuarias, Universidat de Santander (UDES), Bucaramanga, Colombia

Abstract

In spite of the current advances and achievements in cancer treatments, colorectal cancer (CRC) persists as one of the most prevalent and deadly tumor types in both men and women worldwide. Drug resistance, adverse side effects and high rate of angiogenesis, metastasis and tumor relapse remain one of the greatest challenges in long-term management of CRC and urges need for new leads of anticancer drugs. We demonstrate that CRC treatment with the phytopharmaceutical mangiferin (MGF), a glucosylxanthone present in Mango tree stem bark and leaves (Mangifera Indica L.), induces dose-dependent tumor regression and decreases lung metastasis in a syngeneic immunocompetent allograft mouse model of murine CT26 colon carcinoma, which increases overall survival of mice. Antimetastatic and antiangiogenic MGF effects could be further validated in a wound healing in vitro model in human HT29 cells and in a matrigel plug implant mouse model. Interestingly, transcriptome pathway enrichment analysis demonstrates that MGF inhibits tumor growth, metastasis and angiogenesis by multi-targeting of mitochondrial oxidoreductase and fatty acid β-oxidation metabolism, PPAR, SIRT, NFκB, Stat3, HIF, Wnt and GP6 signaling pathways. MGF effects on fatty acid β-oxidation metabolism and carnitine palmitoyltransferase 1 (CPT1) protein expression could be further confirmed in vitro in human HT29 colon cells. In conclusion, antitumor, antiangiogenic and antimetastatic effects of MGF treatment hold promise to reduce adverse toxicity and to mitigate therapeutic outcome of colorectal cancer treatment by targeting mitochondrial energy metabolism in the tumor microenvironment.

Introduction

The latest GLOBOCAN 2018 database showed that colorectal cancer (CRC) is placed third in incidence and second in mortality worldwide, respectively; being responsible for 10% of all newly diagnosed cancer cases, and 9% of all cancer deaths worldwide (). In Cuba, CRC ranked fifth in incidence and third in mortality in 2013, being more frequently in female and in adult people over 60 years of age in both indicators (). Despite advances in the prevention and diagnosis of CRC, between 60 and 70% of patients are detected at middle- or late-stage CRC, which results in difficult clinical management and less favourable prognosis (). This contributes to 40–50% mortality rate within the first 5 years since the diagnosis, accounting for almost 10% of the cancer mortalities each year (Torre et al., 2016). Most of the CRC patients die from metastases (mCRC), predominantly in the liver, lungs, lymph nodes, and the peritoneum (Vatandoust et al., 2015). The main therapeutic options to mCRC are the surgical removal of liver metastases and/or chemotherapy. The former has increased patient survival in the first 5 years, but only 20–25% of these patients are eligible for surgery (). The effectiveness of current chemotherapies is limited () due to suppression of host immune antitumor responses (Teng et al., 2015; ), therapy-induced tumor resistance, tumor relapse (Shaked, 2016) local/systemic toxicities and increased risk of secondary tumorigenesis (; ). In this scenario more than 80% of patients with mCRC die during the long-term follow-up (Ting et al., 2013). The progression of CRC is also closely related to angiogenesis associated with a poor prognosis and relapse of the disease () and has become an alternative valid therapeutic target in the treatment of colon cancer (Wang et al., 2015; Wang and Zhu, 2016). Since 2004, conventional chemotherapy is combined with targeted monoclonal antibody therapies specifically directed against vascular endothelial growth factor (VEGF) and epidermal growth factor (EGF) and their receptors to block angiogenic and proliferative tumor actions (), (). Unfortunatedly, resistance and adverse effects have also been reported for these targeted treatments (Wang et al., 2015; Wang and Zhu, 2016). This has renewed interest in phytochemicals from medicinal plants for CRC treatment. Herbs with a rich ethnomedicinal history represent a rich resource of bioactive compounds with polypharmaceutical effects with multiple targets in cell growth, cell differentiation, and apoptosis regulation (; ; ). The public WHO global report on traditional and complementary medicine indicates that in developing countries, 80% of population relies on plant-derived medicines for the health care (World Health Organization, 2019). The number of phytochemical products submitted to the FDA is particularly high in the oncological area.

Mangiferin (MGF) is a polyphenolic phytochemical of the C-glycosylated xanthone type (1,3,6,7-tetrahydroxyxanthone-C2-β-d glucoside) present in all parts of Mangifera indica L. (Anacardiaceae), and in some medicinal plants. We and others have previously demonstrated promising antiproliferative, anti-angiogenic, anti-invasive and antimetastatic antitumor effects by poorly characterized mechanisms in different cancer cell types ((; ; ; ; ; ; Zou et al., 2017; Tan et al., 2018; ; )).

Taking into account the interesting pharmacological properties of MGF described above, we applied a systems biology transcriptomics approach to resolve the molecular mechanism of action for the in vivo antitumor effects of MGF in a syngeneic immunocompetent allograft mouse model of murine CT26 colon carcinoma.

Materials and Methods

Reagents

Mangiferin (MGF) was purified and characterized from the aqueous extract of the stem bark and leaves of Mangifera indica L to 95,5% purity at the Departments of Pharmacognosy, Institute of Food and Pharmacy, Havana University and Department of Analytical Chemistry of Center of Drugs Research and Development, CIDEM (; ; ). Cisplatin (cis-diamminedichloroplatinum(II), CDDP) was purchased from Shandong Boyuan Pharmaceutical Co., Ltd, China. The culture medium RPMI-1640, dimethylsulfoxide (DMSO), fetal bovine serum, penicillin-streptomycin (100x), L-Glutamine, trypsin, ethylenediaminetetraacetic acid (EDTA), sodium dodecyl sulphate (SDS), reduced Matrigel (196 USP unit/ mg) and set of Drabkin reagents were purchased from Sigma-Aldrich Inc. (Saint Louis, MO, United States). 3-(4,5-Dimethyl-2-thiazolyl)-2,5-diphenyl-2H-tetrazolium bromide (MTT), 2-propanol, ethanol, hydrochloric acid (HCl), formaldehyde and acetic acid were obtained from Merck KGaA (Germany). Solutions of test compounds were routinely prepared in DMSO or MilliQ water and sterilized using syringe-driven 0.22 µm filters. The final concentration of DMSO in culture medium did never exceed 0.5%.

Tumor Cell Lines and Culture Conditions

CT26. WT is an N-nitroso-N-methylurethane-(NNMU) induced undifferentiated murine adenocarcinoma of the colon, syngeneic with the BALB/c mouse, which was provided by the Molecular Immunology Center, Havana. CT26. WT cells were grown in cell culture as monolayers in T-75 cm2 culture flasks (BD Falcon) in RPMI-1640 with 2 mM of L-Glutamine supplemented with 10% fetal bovine serum, 100 UI/mL of penicillin, and 100 μg/ ml of streptomycin. Only cells of the first two serial passages after cryostorage were used in the syngeneic allograft models (). At the day of implantation, the cells were harvested from subconfluent cultures (70–85%) by trypsinization (0.05% trypsin and 0.02% EDTA), and washed twice in phosphate-buffered saline solution before they were resuspended in non-supplemented RPMI culture medium. Human colon adenocarcinoma cell line HT29 (ATCC HTB-38) was cultured in McCoy’s 5 A medium (Gibco, ThermoFisher, United States) supplemented with 1.5 mM L-glutamine, 2.2 g/L sodium bicarbonate, 100 ml/ L of fetal bovine serum and 100 U/ mL penicillin and 100 g/ L streptomycin. All cell lines were cultured in a humidified atmosphere (37°C, 5% CO2).

Cell Viability Assays

MTT assays were performed to assess cell viability and/or growth inhibition of CT26 and HT29 colon cancer cell lines after exposure to MGF essentially, as previously described (). Briefly, cells were seeded in 96-well tissue culture plates at the concentration of 5,000 cells per well. The next day cells were treated with MGF in a concentration range of 10–400 µM for 24, 48, and 72 h. Control cells were exposed to 0.5% DMSO solvent control. Following treatment, the medium supernatant was aspirated from the wells and replaced with fresh medium. Next, 50 µL of MTT solution was added to all the wells and cells were incubated for another 4 h at 37°C. Subsequently, medium was aspirated and formazan crystals were dissolved in 50 µL of DMSO. The measurement of absorption was performed at 595 nm using (Envision 2103 multilabel reader, Perkin Elmer, United States) plate reader. The experiments were repeated with four technical replicates per treatment.

Cell Migration Wound Scratch Assay

HT29 cells were seeded at a concentration of 35,000 cells per well and left to grow at 37°C under 5% CO2 until they reached 95–100% confluence. Then, cells were treated with 10 μg/ ml of MMC for 4 h to block cell proliferation and distinguish it from cell migration. When the incubation time was finished, a scratch was made in each well with the Incucyte Woundmaker tool. Next, cells were gently washed with PBS to remove any cell debris formed after making a scratch. Afterwards PBS was removed and exchanged with serum-free McCoy’s 5 A medium to prevent cell division. The cells were treated either with MGF alone (at 1, 50 and 400 µM concentrations) or in combination with a proangiogenic cell migration growth factor VEGF (at 10 ng/ ml) for 24 h and incubated at 37°C under 5% CO2 in the Incucyte Live-Cell Analysis system (Sartorius, Germany). Control cells were treated with an appropriate solvent, which in the case of MGF was 0.5% DMSO and in the case of VEGF–MilliQ water. Imaging was performed every 2 h for a total duration of 24 h with a × 10 objective. The area of wound closure was calculated and analysed using the Incucyte Scratch Wound Cell Migration software module of the ZOOM imaging system. The experiments were performed in three biological replicates with three technical replicates for each treatment.

Fatty Acid β-oxidation and Ketone Body Assays

Intracellular lipid accumulation induced by free fatty acid exposure leads to enhanced fatty acid oxidation (FAO) and ketogenesis (beta-hydroxybutyrate production). Protein expression levels of CPT1 (carnitine palmitoyltransferase 1), a key regulatory enzyme of mitochondrial fatty acid β-oxidation were determined by Western analysis. Ketone body (beta-hydroxybutyrate) levels were assessed using a colorimetric assay kit (Cayman Chem). Briefly HT29 cells were seeded at a density of 18*106 cells in a T175 flask. The next day cells were treated with 0.2 mM FFA and 40 µM or 400 µM MGF for 24 h. Stock solutions of 0.66 M oleic acid and (1.32 M) palmitic acid (Sigma-Aldrich, Germany) were prepared in isopropanol. Equal amounts of oleic and palmitic acid were mixed to prepare a FFA stock of 1 mM FFA. Free fatty acid-free bovine serum albumin (BSA) was dissolved in serum-free McCoy’s 5 A or RPMI medium without antibiotics at the final concentration of 1% and then sterilized using syringe-driven 0.22 µm filters. Afterwards, medium was supplemented with the mixture of FFA at a final concentration of 0.2 mM and sonicated for 6–8 h until FFA were completely dissolved using a Branson 3,200 sonication bath. FFA medium was protected from light and stored at 4°C. Following treatment, cell pellets were collected with a cell scraper and used immediately to extract ketone bodies according to the manufacturers protocol. The measurement of absorption was performed at 450 nm using (Envision 2103 multilabel reader, Perkin Elmer, United States) plate reader at room temperature. The experiments were performed in one biological replicate, with three technical replicates per each treatment. Similarly, RNA was extracted using RNeasy Mini Kit (Qiagen, Germany) from the collected cell pellets. RNA quantity was determined using Qubit™ RNA Broad Range Assay Kit with the aid of the Invitrogen Qubit™ 4 Fluorometer (Thermo Fisher Scientific, United States). The extracted RNA was stored at -80°C until further QPCR analysis. For protein analysis, cell pellet was washed once with cold PBS, centrifuged (1,500 rpm, 3 min, 4°C) then cells were lysed using 1x RIPA buffer (Cell Signalling, United States) with a cocktail of protease inhibitors (Sigma-Aldrich, Germany) and left on ice for 15 min. Afterwards, cell lysates were centrifuged (13,000 rpm, 15 min, 4°C) and the supernatants were transferred to the new test tubes. Protein concentration was measured using Pierce™ BCA Protein Assay Kit (Thermo Fisher Scientific, United States) and the measurement of absorption was performed at 560 nm using (Envision 2103 multilabel reader, Perkin Elmer, United States) plate reader. Prior to SDS-PAGE electrophoresis, protein lysates were mixed with Laemmli buffer (Biorad, United States) and 50 mM 1,4-dithiothreitol (DTT) and then heated at 70°C for 10 min. Then, the samples and protein ladder (BenchMark™ Protein Ladder, Thermo Fisher Scientific, United States) were loaded on SDS-PAGE polyacrylamide gels (stacking 6%, resolving 12%) at protein concentration equal to 10 µg/well. The electrophoresis was run in Mini-PROTEAN Tetra Cell System (Biorad, United States) using the high molecular weight buffer (100 mM MOPS, 100 mM Tris, 0.2% SDS, 2 mM EDTA, 5 mM sodium bisulfite) at 70 V until the samples reached the separating gel and then the voltage was increased up to 120 V. Afterwards, the proteins were transferred on pre-wet nitrocellulose membranes (Cytiva, United States) in Power Blotter 1-step transfer buffer (Invitrogen, United States) using Power Blotter XL System (Invitrogen, United States) (7 min at 1.3 A). Nitrocellulose membranes were blocked in blocking buffer (5% milk solution in TBST) for 1 h at room temperature and incubated overnight at 4°C with primary antibodies diluted in blocking buffer as follows: 1:4000 for anti-CPT1A (Proteintech, UK). Next day, membranes were washed three times with TBST and then incubated with HRP-conjugated anti-rabbit secondary antibody diluted in blocking buffer (1:2000). Anti-GAPDH (Cell signaling, United States) antibodies (diluted 1:2000) in blocking buffer were used as loading controls for CPT1A. Protein detection was performed on the Amershan imager 680 (Cytiva, United States) using Pierce™ ECL Western Blotting Substrate (Thermo Fisher Scientific, United States) and quantified using Image-J software.

In vivo Experiments

Housing and all procedures involving animals were performed according to protocols approved by the Drug Research and Development Center Ethics Committee for Animal Research and in compliance with the Cuban National Guidelines for the Care and Use of Laboratory Animals (CECMED Resolution No. 37 2012 for Good Laboratory Practice for Quality Control). Male BALB/c mice were purchased at National Laboratory Animal Production Center, Havana, Cuba. BALB/c mice were used at 6–8 weeks of age with weight between 18 and 22 g, in all experiments. They were supplied with water and food at libitum and were kept in quarantine for 15 days for adaptation to temperature (22°C ± 2°C), humidity (77 ± 3%) and alternating cycles of 12 h light/dark environment. For allograft studies, the animals were inoculated (as specified in more detail for each specific model below) with 1x105 CT26. WT cells/mouse in 100 μL of non-supplemented RPMI-1640 culture medium. The following day, they were weighted, labeled and placed in separate boxes.

Ectopic Subcutaneous Syngeneic CT26. WT Tumor Model

The animals were inoculated subcutaneously with CT26. WT cells on the right dorsal side. Once tumors grew up to 30 mm3, mice were randomly divided at day 14 into five groups (10 animals/group) and treated for 14 days (28 days protocol). The products, doses, and treatment schedules were as follows: Group I, carboxymethylcellulose (CMC) 0.05% as the negative control; Group II, MGF 10 mg/ kg b. w.; Group III, MGF 50 mg/ kg b. w.; Group IV, MGF 100 mg/ kg b. w.; and Group V, CDDP 5 mg/kg b. w. as positive control. MGF doses were prepared in 0.05% CMC and CDDP in sterile saline solution (0,9% NaCl). MGF and CMC were administered daily intragastrically and CDDP every 6 days by intraperitoneal injections in three treatment cycles. The tumor growth was monitored and tumor size was measured by length and width with caliper every 2 days by the same operator, and the tumor volume was calculated (V, mm3 = length x width2 x 0.5). At the end of the experiment, all mice were sacrificed and dissected. The tumors were recovered for photographic imaging and weighed. Three independent experiments were performed.

Lung Metastasis Model

The animals were challenged by an intravenous puncture in the tail veins and distributed in five groups (10/group). The products, doses, and treatment regimens were as follows: Group I, 0.05% CMC; Group II, MGF 10 mg/kg b. w.; Group III, MGF 50 mg/ kg b. w.; Group IV: MGF 100 mg/ kg b. w.; and Group V, CDDP 5 mg/kg b. w. The preparation, route, and frequency of administration were the same as those described in the ectopic/subcutaneous model. The animals were treated for 17 days and sacrificed at day 18. The trachea-lung were extracted, photographed and perfused with 20 Gx1½ tracheal needles with 1 ml 15% Chinese ink in PBS pH 7.2–7.4. The trachea was excised and the lungs were soaked for discoloration and fixation for 5 min in Fekete’s solution (). The lungs were photographed and the superficial metastatic nodules were counted using a stereoscope (Olympus, Japan) and fixed with 4% formaldehyde. Lungs presented with more than 250 nodules, were assigned a cut-off value of 250 (due to the difficulty of counting a larger number because of high nodule density). A survival analysis was performed until the treatments were completed. Three independent experiments were performed.

Angiogenesis Model by Matrigel Implant

The animals were divided into four groups (5/group). The control groups were implanted with 200 μL of a cold mixture of Matrigel, heparin (32 IU) and PBS (negative or basal control: implant A) and another mixture equal to the previous one plus CT26. WT cells (5 x 105) as positive control: implant B. For the treatment groups, two concentration of MGF: 200 μg/ ml (implant C) and 100 μg/ ml(implant D) were included into the implant volume B. The implants were subcutaneously applied in the in the inguinal region. The experimentation time was extended by 12 days. The animals were sacrificed, the plugs were excised, photographed and cleaved with residues of surrounding tissue to facilitate histological orientation. They were weighed and measured with caliper by their length, width, and height, and the volume was calculated using the formula of an ellipsoid: V (mm3) = length x width x height x π/6). The plugs were liquified in 300 μL PBS by incubation at 4°C, overnight (), and the hemoglobin (Hb) released was quantified by the Drabkin method (). For histological analysis, the plugs were fixed in 10% formalin in PBS, embedded in paraffin and cut (5–7 μm). Histological staining was initially performed with hematoxylin-eosin and whereas specific staining for blood vessel counts was done with Masson’s trichrome. A microscopic count of the blood vessels was performed throughout the periphery of the tumor.

Immunohistochemical Studies

Subcutaneous tissue tumors were fixed in formalin (10% w/v in phosphate-buffered saline–PBS pH 7.4) during 24 h and sections (4 µm) were processed for immunohistochemistry analysis. A biotin-peroxidase system was used with identification of the secondary antibody by polymer (ADVANCE HRP). Antigen retrieval was performed in water bath at 98°C, with a citrate buffer solution (pH 6.0) (Target Retrieval Solution) for 20 min. In order to block the endogenous peroxidase activity, slides were incubated with a solution of H2O2 3% in methyl alcohol. Reagents were applied manually with 1 h incubation for the primary antibody and 30 min for the reagents, except for chromogen diaminobenzidine, incubated for 5 min. Minichromosome maintenance complex component 7 (MCM 7 or CDC47), a marker of cell proliferation, was measured using a specific monoclonal antibody CDC47 (clone 47DC141, Neomarkers, 1/300) (Zorde Khvalevsky et al., 2013). Positive and negative controls were included in each bath. The index for CDC47 staining was obtained by counting the percentage of positive cells in 500 tumor cells.

RNA Extraction and Isolation

Subcutaneous and pulmonary tumor tissue samples from each animal (n = 5/group) were individually treated with RNAlater® (Qiagen, Benelux, Belgium) by 24 h at 4°C, and immediately stored at −80°C until the RNA isolation procedure. Up to 30 mg tissue were completely disrupted and homogenized in 700 μL QIAzol® Lysis Reagent (Qiagen, Benelux, Belgium) by the Tissue Lyser II (Qiagen®, Benelux, Belgium) with 5 mm stainless steel beads diameter to extract total RNA. Total RNA was isolated using RNeasy® Microarray Tissue kit (Qiagen, Benelux, Belgium), according to manufacturer’s protocol. Total RNA was purified on RNeasy Mini spin columns according to the manufacturer’s protocol. During the purification process contaminating genomic DNA was removed by means of DNase digestion. Purity and concentration total RNA samples were evaluated by ultraviolet spectroscopy (NanoDrop, Thermo Scientific, Wilmington, United States). The integrity of the RNA samples was evaluated by the Experion Automated Electrophoresis System (BioRad, CA, United States) according to the manufacturer instructions. Good quality of the RNA material samples with A260/280 ratio between 1,9–2,1 were used for further analysis.

RNA Bead Array and Data Processing

Five hundred ng of total RNA from 18 samples (3 replicates by group) subcutaneous and pulmonary tissue tumors were amplified using the Illumina TotalPrep RNA Amplification kit (Life Technologies, Carlsbad, CA, United States). RNA was reverse transcribed using T7 oligo(dT) primers after which biotinylated cRNA was synthesized through in vitro transcription reaction. Seven hundred 50 ng of amplified cRNA was hybridized to a corresponding array on a Mouse WG-6 2v bead chip (Illumina, San Diego, CA, United States). The beadchip array was incubated for 18 h at 58°C in a hybridization oven whilst continuous rocking. After several consecutive washing steps, bead intensities were read on Illumina iScan. Raw data intensities were read in R using the “bead array” package (v2.8.1) (). Intensities were quantile normalized and differential gene expression between samples was estimated using “limma” (v3.14.1) (). Genes were considered to be differentially expressed when they had a p-value < 0.05 and a fold change of at least 1.25 in comparison to controls. Transcriptomic analysis was done by uploading Illumina bead array gene expression data into Ingenuity Pathway Analysis (IPA) software (Ingenuity® Systems, www.ingenuity.com, Redwood City, CA, United States and QIAGEN Inc., https://digitalinsights.qiagen.com/products-overview/discovery-insights-portfolio/analysis-and-visualization/qiagen-ipa/?), and performing a core analysis according to the instructions provided. Each identifier was mapped to its corresponding gene object in the Ingenuity knowledge base. Fisher’s exact test was used to calculate a p-value determining the probability that each biological function and/or disease assigned to that data set is due to chance alone. Metascape systems biology freeware (https://metascape.org/) was used for correlating the transcriptomic profiles of the different in vivo data (Zhou et al., 2019). The Circos plot visualization allows to show how genes from different input gene lists overlap. Heatmaps show Metascape enrichment analysis of all statistically enriched ontology terms (GO/KEGG terms, canonical pathways, hall mark gene sets). Accumulative hypergeometric p-values and enrichment factors are calculated and used for filtering. Remaining significant terms are then hierarchically clustered into a tree dendrogram based on Kappa-statistical similarities among their gene memberships. The term with the best p-value are selected within each cluster as a representative term to be displayed in a hierarchical tree dendrogram. The heatmap cells are colored by their p-values (see color legend). Along the same line, Metascape enrichment analysis of all statistically enriched transcription factor (TF)-target interaction networks is dermined by the TRRUST database (). Protein-protein interactions (PPI) among all input gene lists are extracted from PPI data source to form a PPI network. GO enrichment analysis is applied to the network to assign biological “meanings” of sub-protein networks. GO enrichment analysis is applied to each MCODE network to assign “meanings” to the network component, where top three best p-value terms were retained. MCODE components were identified from the merged network. Each MCODE network is assigned a unique color.

Reverse Transcription-PCR and Real-Time Quantitative PCR Validation

Total RNA (700 ng) isolated by RNeasy Mini Kit (Qiagen), was reverse transcribed to cDNA by the GoScript™ Reverse Transcriptase kit (Promega, United States) according to the manufacturer’s protocol. RT-qPCR was performed using GoTaq® qPCR master mix (Promega, United States) on Rotor-Gene Q (Qiagen, Germany). Gene expression was calculated with the ΔΔCt method after normalization against the glyceraldehyde-3-phosphate dehydrogenase (GAPDH), actinB (ACTB) or beta-2-microtubulin (B2M) housekeeping genes, as further specified. Following gene specific primersets were used for array validation experiments of C-C Motif Chemokine Ligand 3 (CCL3), Collagen Type VI Alpha three chain (COL6A3), C-X-C Motif Chemokine Ligand 9 (CXCL9) and Metalloproteinase 13 (MMP13); and the following gene specific primer sets: CCL3 forward primer: 5′-CAG​CCA​GGT​GTC​ATT​TTC​CTG-3′ and reverse primer: 5′-AGG​TCT​CTT​TGG​AGT​CAG​CG-3´. COL6A3 forward primer: 5′-GGA​ACC​ACG​GAA​GAG​AGC​AA-3′ and reverse primer: 5′-CAGGGAACTGACCCAAGACA-3´.CXCL9 forward primer: 5′- TGT​GGA​GTT​CGA​GGA​ACC​CT-3′ and reverse primer: 5′- AGT​CCG​GAT​CTA​GGC​AGG​TT-3´. MMP13 forward primer: 5′- GGA​GCC​CTG​ATG​TTT​CCC​AT-3′ and reverse primer: 5′- GTC​TTC​ATC​GCC​TGG​ACC​ATA-3´. The experiment was carried out in duplicate in three independent biological replicates. The thermal cycling conditions included initial denaturation at 95°C for 2 min and subsequent 40 cycles of two-step protocol: denaturation at 95°C for 15 s and annealing/extension at 60°C for 60 s.

Statistical Analysis

One-way or two-way ANOVA with Dunnett post hoc test (comparison every mean to control mean) were performed to evaluate statistical significance, as indicated in Figure legends. Survival analysis was done with the log-rank (Mantel-Cox) statistical test. All statistical analyses were performed using GraphPad Prism® (GraphPad Software, La Jolla California, United States) statistical package. P value < 0.05 was considered significant.

Results

Mangiferin Induces CT26.WT Colorectal Tumor Regression in vivo

The CT26 syngeneic colorectal cancer tumor model was used to investigate the effect of MGF on tumor growth in vivo. CT26. WT cells were grown for 14 days after subcutaneous injection. Changes in tumor volume by MGF treatment were measured during 2 weeks of treatment protocol. As shown in Figure 1A, allgrafts treated with MGF show dose dependent tumor regression, comparable to the volume reduction observed upon treatment with the reference chemotherapeutic drug CDDP. Significant treatment antitumor effects start to be observed from day 22, following 8 days of treatment (Figure 1A), reaching a maximal reduction in tumor size of 75% at 100 mg/kg. Overall, CDDP treatment revealed the strongest antitumor effect (approximately 90% tumor size reduction). Along the same line, MGF and CDDP dependent reduction of tumor volume corresponds with decreased immunohistochemical staining intensity of the cell proliferation marker CDC47 reflecting arrest of tumor growth (Figures 1B,C).

FIGURE 1

Mangiferin Inhibition of Tumor Cell Motility (Invasion), Angiogenesis, Vasculogenesis Involves Transcriptional Changes in Fatty Acid Metabolism, PPAR, NFκB, Wnt and GP6 Signaling Pathways

To get more insight in the molecular mechanism of the antitumor of MGF, a gene expression microarray experiment was performed. Gene expression profiles of subcutaneous tumor tissue obtained from mice treated for 2 weeks with 5 mg/ kg CDDP, 10 mg/ kg MGF, 50 mg/ kg MGF and 100 mg /kg were compared to untreated control mice with tumors. In general, MGF showed mild transcriptional changes as compared to CDDP treatment. Using a fold change difference of at least 1.25 and a p-value below 0.05, 454 genes were found to be differentially expressed in the 10 mg/ kg MGF condition, 263 genes in the 50 mg/ kg MGF condition, 812 genes in the 100 mg/ kg MGF condition and 3,340 genes after 5 mg/ kg CDDP treatment (Supplementary Table S1). Highly significant positive correlations could be observed upon comparing the log2 microarray gene expression values and the ΔCt qPCR expression values of CCL3, COL6A3, CXCL9 and MMP13 target genes, supporting the validity our microarray results (Supplementary Figure S1).

Next, bioinformatic analysis of transcriptome changes of the different treatments was performed by correlating the -log10 p-values for every treatment combination (Figure 2A). Whereas strong positive correlations were observed for transcriptional responses for the different MGF doses (r = ≥ 0.5), only weak overlap could be observed for MGF and CDDP treatments (r = 0.13), suggesting a different mechanism of action (Figure 2B). Next, we focused on pathway enrichment analysis of disease and biological pathways of the MGF treatment specific transcriptional changes using Ingenuity Pathway Analysis (IPA®, QIAGEN Redwood City) (Figures 2C,D; Supplementary Table S1). Differentially expressed genes for all MGF doses are enriched for regulation of cancer and colon cancer related disease functions, including cell death and survival and cell motility, gastrointestinal neoplasia, gastrointestinal carcinoma, malignant neoplasm of large intestine, and digestive system cancer (Figure 2C).

FIGURE 2

Along the same line, significant transcriptional changes could be observed at the gene expression level in cytokine signaling (acute phase response, oncostatinM), mitochondrial metabolic stress (oxidative phosphorylation, fatty acid oxidation, Nrf2, Glutathione detoxification, SIRT, and AMPK/mTOR signaling), cell cycle (cyclins, checkpoints, mitotic PLK kinases), metalloproteases (MMPs), cell death (death receptors), iron homeostasis (heme degradation) DNA damage (ATM, BRCA1, GADD45 signaling), cell death (ubiquitination), differentiation (Wnt/catenin) and xenobiotic-metabolizing hormone receptor signaling (Aryl hydrocarbon and aldosterone receptor signaling) (Figure 2D; Supplementary Figure S3; Supplementary Table S1).

Next, we used the IPA regulation z-score algorithm to identify biological functions and pathways that are most significantly activated (Z-score >1) or repressed (Z-score < -1) by MGF based on direction of gene expression changes (increase/decrease). Interestingly, MGF dose dependently suppresses expression of multiple genes involved in tumor cell motility (invasion), angiogenesis, vasculogenesis and development of vasculature functions (Figure 3A). Furthermore, when analyzing possible protein interaction networks of MGF responsive gene targets by the STRING algorithm (https://string-db.org/), we identified various interconnected protein network clusters involved in cell motility, angiogenesis and cell death/autophagy/proliferation, or dysfunctional mitochondrial oxidoreductase metabolism, in line with the results of our IPA analysis (Figures 2D, 3A and Supplementary Figure S2).

FIGURE 3

Mangiferin Suppresses Cell Motility in vitro

To further characterize MGF anticancer effects on cell viability and cell motility in vitro, we next performed MTT and wound scratch assays in mouse CT26 and human HT29 colon cancer in vitro cell models. First, CT26 and HT29 were treated with MGF for 24, 48, and 72 h with a MGF concentration range from 10 to 400 μM, followed by a colorimetric MTT assay, as previously described (). As can be observed from Figure 3B, MGF does not significantly decrease cell viability in the dose range tested and shows no adverse toxicity. Next, MGF effects on prometastic and proangiogenic cell motility were analyzed for 24 h in a wound scratch assay in VEGF treated HT29 cells by realtime Incucyte Live Cell imaging in IncuCyte® ImageLock 96-well microplates, after scratching cells with the Incucyte Woundmaker tool (). Time dependent changes in wound closure of the different treatments were quantified by the Incucyte Scratch Wound Cell Migration software module of the ZOOM imaging system. In line with our gene expression pathway analysis in mouse CT26 cells, MGF treatment was also found to significantly suppress VEGF induced prometastatic and proangiogenic cell motility in human HT29 cells (Figures 3C,D).

Mangiferin Exerts Antiangiogenic Effects in vivo in a Matrigel Plug Angiogenesis Mouse Model

Next, we further verified whether inhibition of vasculogenesis by MGF could also be experimentally confirmed in a CT26. WT cells engrafted matrigel plug assay, allowing to measure effects on neovascularization of tumor cells in vivo (). Angiogenesis is an attractive cancer therapeutic target, not only because it supplies oxygen and nutrients for the survival of tumor cells, but also provides the route for metastatic spread of these cancer cells.

Tumor formation in the subcutaneous matrigel plug tissue, with abundant blood supply, was observed in the CT26. WT engrafted control animals. In contrast, the plug in control animals colors yellow lacking blood irrigation in the matrigel subcutaneous implant (Figure 4A). Clearly, MGF treatment promotes a reduction of peritumoral and intratumoral blood vessel supply and tumor size in a dose-dependent manner (Figure 4A from left to right). Moreover, quantitative analysis of Hb content shows that MGF dose dependently inhibits angiogenesis, i.e. 58% and 81% decrease of Hb levels at respectively 100 and 200 μg/ml MGF treatment (p < 0.05) (Figure 4B). As expected, CT26. WT tumor cells show an adequate proangiogenic response in matrigel (81% increase in Hb content) as compared to the control setup (Figure 4B). Besides, the antiangiogenic effect of MGF also results in dose dependent reduction in tumor volume (84% reduction at 200 μg/ ml MGF, Figure 4C) or tumor weight (73% reduction at 200 μg/ ml MGF, Figure 4D). Besides, a histological study was performed using the hematoxylin-eosin and Masson’s trichrome staining. This allows clear differentiation of the collagen fibers (colored in blue) from the blood vessels (colored in red). As can be observed in the histological staining (Figure 4E) or blood vessel count (Figure 4F), MGF dose dependently inhibits the formation of blood vessels in the tumor periphery of the CT26. WT engrafted Matrigel (p < 0.001). Altogether, these results support our gene expression pathway enrichment analysis which identified various MGF responsive target genes involved in suppression of angiogenesis-vasculogenesis.

FIGURE 4

Mangiferin Inhibits Lung Metastasis of CT26 Cells in an Experimental Metastasis Mouse Model

To further evaluate suppression of cancer cell motility-migration in response to MGF, as predicted in our transcriptome analysis, we next investigated whether MGF could reduce number of lung metastasis nodules upon intravenous puncture of colorectal cells in the mouse tail veins. The number of CT26 metastasis nodules in lungs were quantified following 17 days treatment in the different experimental setups as shown in Figure 5. A significant reduction (p < 0.001) in the number of pulmonary surface metastatic nodules could be observed by increasing doses of MGF, respectively 46, 77, and 88% decrease upon treatment with 10, 50 or 100 mg/ kg b. w. MGF as compared to the vehicle control group (Figure 5A). Similarly, metastatic nodules could be reduced 83% upon treatment with 5 mg/kg b. w. CDDP (83%) (Figure 5A). Treatments with 50 and 100 mg/ kg b. w. MGF, or 5 mg/ kg b. w. CDDP also resulted in a significant improvement of longterm survival rate as compared to the vehicle control group (MGF 50 mg/ kg: p = 0.027; MGF 100 mg/ kg: p = 0.020; CDDP: p = 0.022) (Figure 5B). Remarkably, gene expression analysis of pulmonary metastatic colon tumor tissue biopsies from mice treated with MGF revealed similar gene expression patterns and pathway enrichment as in the primary CT26 tumor model (Figure 5C; Supplementary Table S2). Accordingly, MGF treatment targets mitochondrial energy metabolism (fatty acid oxidation, sirtuin, TCA) NFκB, PPAR, Wnt and GP6 signaling to suppress cancer cell migration-metastasis processes (Figure 5C). Since tumor tissue biopsies are a heterogenous mixture of CT26 tumor metastasis and healthy lung cells, tumor gene expression changes by MGF treatment at 100 mg/kg are less pronounced than at 50 mg/kg since the fraction of CT26 tumor cells is strongly decreased at the highest MGF concentration, which significantly reduces the MGF effect size in the mixed cell population.

FIGURE 5

Mangiferin Regulated Protein Interaction Networks Target Mitochondrial Energy Metabolism and Inflammatory Signaling to Elicit Antitumor and Antimetastatic Therapy Response in vivo

To further characterize molecular mechanism of action of MGF, we further crosscompared MGF responsive genes in the primary CT26 tumor model and the metastatic cancer models by Metascape systems biology freeware (https://metascape.org/) (Zhou et al., 2019) (Supplementary Figure S4; Supplementary Table S3). As can be observed in the circos plot (Figure 6A), we observe significant overlap between gene expression changes as well as biological pathways, as both experiments probably target similar biological processes. This is also clear in the overlapping heatmap representation (Figure 6B), showing highly significant pathway enrichment related to blood vessel development, lipid metabolism, chemotaxis, locomotion, wound healing, generation of energy metabolites, hormone stimulus, growth factor response and leukocyte activation among others (Supplementary Table S3). Further TRRUST analysis of all statistically enriched TF-target interaction networks () identified major regulatory roles for STAT3, PPARA/G, HIF1A and NFκB1/RelA transcription factor families, which are key players in lipid metabolism, hypoxia and inflammatory tumor signaling (Figure 6C). Furthermore, Metascape protein-protein-interaction (PPI) analysis of all MGF responsive target genes reveals a global protein network involved in reprogramming of mitochondrial energy metabolism in CT26 tumor cells (Figures 6D,E). More specifically, 4 major protein interaction subnetworks could be identified (Figure 6F) involved in fatty acid β-oxidation and toll like receptor signaling (MCODE1), generation of energy metabolites (oxidation organic compounds, ribose phosphate metabolism) (MCODE2), mitochondrial respiratory chain complex 1 (MCODE3) and NFκB signaling (MCODE4) (Figure 6G).

FIGURE 6

Mangiferin Targets Fatty Acid β-oxidation Metabolism in Human HT29 Colon Cancer Cells

To verify whether MGF directly targets fatty acid β-oxidation energy metabolism in colon cancer cells, we measured ketone body concentrations (β-hydroxybutyrate) upon exposure of human HT29 colon cancer cells to free fatty acids (to mimick free fatty acids present in the tumor microenvironment in vivo) in presence or absence of MGF. Figure 7A shows that 24 h treatment of HT29 cells with 0.2 mM free fatty acids (FFA) increases intracellular concentrations of β-hydroxybutyrate as compared to untreated HT29 cells. Of special note, cotreatment with 400 µM MGF reverses this effect to levels observed in untreated cells. Along the same line we measured changes in protein expression of carnitine palmitoyltransferase (CPT1), a key metabolic enzyme involved in fatty acid β-oxidation. Similarly, FFA treatment increases CPT1 expression levels whereas, MGF reverses the effect to control levels (Figure 7B).

FIGURE 7

Discussion

Today, various pharmacological anticancer properties of MGF have already been described in vitro, including antiproliferative, pro-apoptotic, anti-invasive, antimetastatic and antiangiogenic bioactivities (; ; ; ; ). Here, we further demonstrate in vivo antitumor efficacy of MGF in different syngeneic immune competent colon cancer (CT26) mouse allograft models and propose a novel “metabolic” anti-cancer mechanism of action.

To further characterize the corresponding molecular mechanisms involved in tumor regression we have applied systems biology pathway enrichment analysis (IPA) of gene expression patterns of tumor samples of the various treatments. This revealed that cell motility, invasion and metastasis are the most significantly suppressed biological pathways upon MGF treatment. This could be confirmed human and mouse colon cancer cells in a wound scratch assay in vitro and a lung metastasis model in vivo. At the gene expression level, MGF suppresses various prometastatic genes, including MMP13, COL4A1, COL4A2, INHBA, KMT2A, ZEB1/2. Tumour invasion and metastasis involves degradation of different components of the extracellular matrix and basement membranes (Type IV collagen), catalysed by proteolytic enzymes such as matrix metalloproteinases (MMPs). Degradation of type IV collagen correlates with metastatic potential, whereas MMP-13 gene expression has been identified as a reliable predictor of liver metastasis in patients with CRC (Yamada et al., 2010).

Serum levels of Collagen IV (Col4A1, Col4A2) abundant in are elevated in patients with colorectal liver metastases and facilitate invasion (Wang et al., 2020; ). The heterodimeric TGF-β family ligand Inhibin (INHA) is a novel paracrine factor involved in cancer angiogenesis and metastasis (Singh et al., 2018a). The histone methyltransferase KMT2A promotes colorectal cancer metastasis via activation of cathepsin Z (). Finally PDPK1 and ZEB1/2 members drive tumor metastasis spread via regulating epithelial–mesenchymal transition (EMT) during hypoxia (; ; Stemmler et al., 2019; Saxena et al., 2020). Since angiogenesis is one route for metastatic spread of cancer cells, it can be noticed that various transcriptional changes blocking angiogenesis, indirectly also reduce metastatic cell migration and both pathways are interconnected (Supplementary Figure S2). Accordingly, MGF decreased expression of various proangiogenic genes such as HMOX1, MMP13, IGF1R, INSR, TGFB3, PDGFR A/B and SEMA5A among others. HMOX1 promotes vascularisation in a hypoxic environment by stimulating paracrine proangiogenic functions of bone marrow derived cells (BMDCs) (). Besides, HMOX1 also promotes proangiogenic effects in endothelial cells by reciprocal regulation of VEGF and SDF1 (). TGFB3 induces angiogenesis by promoting transdifferentiation of mesenchymal stem cells into endothelial cells (). PDGFRA and PDGFRB play different roles in neosvascularisation and promoting vessel stability (Zhang et al., 2009). Matrix metalloprotein (MMP13) by degrading extracellular matrix (ECM) stimulates the release of VEGF from ECM and thereby promoting angiogenesis and metastasis in skin squamous cell carcinoma and colitis induced colon cancer (; ). SEMA5A promotes angiogenesis by increasing endothelial cell proliferation, migration (). Furthermore, oncogenic angiogenic and metastatic functions have been described for INS and IGF1R (). Finally, in line with the bioinformatic analysis, the predicted inhibition of angiogenesis, and suppression of vasculogenesis could also be confirmed in an in vivo matrigel plug implant mouse model.

By crosscomparing MGF responsive genes in the primary CT26 tumor model and the metastatic cancer models by Metascape systems biology freeware (https://metascape.org/), we searched for a common molecular mechanism of action (Figure 6). Interestingly, significant enrichment of various common pathways could be observed related to blood vessel development, lipid metabolism, chemotaxis, locomotion, wound healing, generation of energy metabolites, hormone stimulus, growth factor response and leukocyte activation among others Furthermore, we identified major regulatory roles for STAT3, PPARA/G, SIRT, HIF1A and NFκB1/RelA transcription factor families, which are key players in lipid metabolism, hypoxia and inflammatory tumor signaling. Ultimately, protein-protein-interaction (PPI) analysis of all MGF responsive target genes established a global protein network with 4 major protein interaction subnetworks involved in fatty acid β-oxidation and toll like receptor signaling (MCODE1), generation of energy metabolites (oxidation organic compounds, ribose phosphate metabolism) (MCODE2), mitochondrial respiratory chain complex 1 (MCODE3) and NFκB signaling (MCODE4). MGF regulation of fatty acid β-oxidation metabolism could further be experimentally validated in human HT29 colon cancer cells by colorimetric detection of decrease in β-hydroxybutyrate concentrations and protein expression of the key metabolic enzyme CPT1. Ketone bodies are interwoven with crucial mammalian metabolic pathways such as β-oxidation (FAO), the tricarboxylic acid cycle (TCA), gluconeogenesis, glycolysis and the pentose phosphate pathway (PPP) (), which were identified in Metascape protein-protein interaction analysis as key protein (sub)networks affected by MGF treatment (Figure 6). A common feature of cancer cells is their ability to rewire their metabolism to sustain the production of ATP and macromolecules needed for cell growth, division and survival. These alterations are pivotal to the development and maintenance of the malignant phenotype of cancer cells in unfavorable tumor microenvironment or metastatic sites. In particular, the importance of altered fatty acid metabolism in cancer has received renewed interest as, aside their principal role as structural components of the membrane matrix, they are important secondary messengers, and can also serve as fuel sources for energy production (; ). The ways in which cancer cells utilize lipids are often influenced by the complex interactions within the tumor microenvironment and adjacent stroma (; ). Cancer cells survive in challenging hypoxic or nutrient depleted tumor micro-environments. For example tumor-stroma localized adipocytes can be activated by cancer cells to lipolyze their triglyceride stores, delivering secreted fatty acids to cancer cells for uptake through numerous fatty acid transporters. Cancer-associated fibroblasts are also implicated in lipid secretion for cancer cell catabolism and lipid signaling leading to activation of mitogenic and migratory pathways (; ). As these cancer-stromal interactions are exacerbated during tumor progression, fatty acids secreted into the microenvironment can impact infiltrating immune cell function and phenotype. Regulation of TLR/NFκB/PPAR pathways by (un)saturated fatty acids plays a key role in tumor-immune responses (). Lipid metabolic abnormalities such as increased fatty acid oxidation and de novo lipid synthesis can provide survival advantages for the tumor to resist chemotherapeutic and radiation treatments and alleviate cellular stresses involved in the metastatic cascade. Fatty acid oxidation is also reprogrammed in cancer-associated immune and other host cells, which may contribute to immune suppression and tumorpromoting microenvironment. New evidence indicates that cancer stem cells show high metabolic plasticity and can dynamically transform their metabolic state to favor glycolysis or oxidative metabolism depending on the tumor microenvironment (; Zhu et al., 2020). This illustrates that fatty acids regulate each part of the cancer lifecycle and suggests that therapeutic intervention targeting lipid and fatty acid metabolism signaling pathways by MGF holds promise for colorectal cancer treatment (; ; ). Fatty acid catabolism and anabolism pathway crosstalk are pivotal in cell fate decision during redox regulated ferroptosis cancer therapy and maybe strongly context dependent (; ; Yuan et al., 2021). In clinical trials, MGF treatment was already found to improve metabolic health (diabetes, hyperlipidemia, insulin resistance) by optimization of mitochondrial bioenergetic pathways (fuel switching between fatty acid and carbohydrates) via sirtuin (SIRT) and PPAR (in)dependent metabolic mechanisms (; ; ; ; Singh et al., 2018b; ; ; Zhang et al., 2019). Sirtuin NAD + -dependent histone deacetylase enzymes are central players in the maintenance of cellular energy and metabolic homeostasis, whereas impaired SIRT functions result in several pathological conditions and contribute to the altered metabolic phenotype (metastasis, drug resistance) of malignantly transformed (colorectal) cancer cells, referred to as the Warburg effect (; ). Natural product SIRT modulators (such as MGF) are receiving growing interest to optimize mitochondrial bioenergetic pathways in cancer treatment (; ). Similarly, PPAR transcription factors are fatty acid sensors that enhance mitochondrial fatty acid β-oxidation (; ). Besides their lipid metabolic effects, PPARs also suppress angiogenesis, metastasis by targeting HIF, NFκB and STAT3 TF functions which regulate cancer hypoxia, inflammatory pathways, chemotaxis, cell adhesion (tight junctions), epithelial-mesenchymal transition (wnt signaling) in the tumor-microenvironment ((; Wang et al., 2004; ; ; Zhou et al., 2012; Woolf et al., 2015; ; ; Vu and Datta, 2017; Vallée and Lecarpentier, 2018; ), in line with our systems biology network analysis. Today, various natural and pharmacological PPAR ligands have already been described with anti-metastatic and anti-angiogenic therapeutic properties (; ), which increase the efficacy of chemotherapeutics with minimal toxicity (; ). Similarly, our results further suggest that MGF has potential value to be used as adjuvant phytotherapeutic treatment against CRC to reduce adverse toxicity and to mitigate outcome of colon cancer treatment. To further develop its application in clinical practice, more molecular combination studies with different MGF formulations (; ; ) as well as uniformly controlled clinical trials are needed.

Statements

Data availability statement

All transcriptome data have now been deposited in a public depository (https://zenodo.org/) and can be accessed viahttps://doi.org/10.5281/zenodo.5185362.

Ethics statement

The animal study was reviewed and approved by Drug Research and Development Center Ethics Committee for Animal Research and in compliance with the Cuban National Guidelines for the Care and Use of Laboratory Animals (CECMED Resolution No. 37 2012 for Good Laboratory Practice for Quality Control). No human studies were conducted.

Author contributions

Conceptualization, JR-G, IH-B, WVB, RD-H, and IR-G; methodology, JR-G, IH-B, OQ-M, CP-N, EL, CT, PJ, GVC, DC, and ML; software, KD; validation, JCR-G, IH-B; formal analysis, KD, JR-G, IH-B, and WVB; investigation, JR-G, IH-B, OQ-M, CP-N, EL, GVC, ML, WVB, RD-H, IR-G; resources, WVB, RD-H, IR-G; data curation KD, WVB; writing, review and editing, JR-G, IH-B, WVB, RD-H, and IR-G; supervision, WVB, RD-H, IR-G; project administration, WVB, RD-H, IR-G; funding acquisition, WVB, RD-H, IR-G; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the P211LH005-019 project (CITMA), Cuba; Visiting Research Program, PVE400768/2014–3 project (CNPq), Brazil, and the projects ZEIN2011PR383 and ZEIN2016PR418 (Flemish VLIR collaboration program), Belgium. This work was partially supported by the project PN223LH010-035 from Ministerio de Ciencia, Tecnologia y Medio Ambiente (Cuba), Programa de Ciencias Basicas y Naturales.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.670167/full#supplementary-material

References

  • 1

    AbdullahL. N.ChowE. K. (2013). Mechanisms of Chemoresistance in Cancer Stem Cells. Clin. Transl Med.2 (1), 3. 10.1186/2001-1326-2-3

  • 2

    AiF.ZhangX.LiX.QinZ.YeQ.TianL.et al (2015). Up-regulation of Matrix Metalloproteinases in a Mouse Model of Chemically Induced Colitis-Associated Cancer: the Role of microRNAs. Oncotarget6 (7), 54125425. 10.18632/oncotarget.3027

  • 3

    ApontesP.LiuZ.SuK.BenardO.YounD. Y.LiX.et al (2014). Mangiferin Stimulates Carbohydrate Oxidation and Protects against Metabolic Disorders Induced by High-Fat Diets. Diabetes63 (11), 36263636. 10.2337/db14-0006

  • 4

    AtanasovA. G.ZotchevS. B.DirschV. M.SupuranC. T.International Natural Product Sciences Taskforce (2021). Natural Products in Drug Discovery: Advances and Opportunities. Nat. Rev. Drug Discov.20, 200216. 10.1038/s41573-020-00114-z

  • 5

    AyerA.ZarjouA.AgarwalA.StockerR. (2016). Heme Oxygenases in Cardiovascular Health and Disease. Physiol. Rev.96 (4), 14491508. 10.1152/physrev.00003.2016

  • 6

    BaánA.AdriaensensP.LammensJ.Delgado HernandezR.Vanden BergheW.PietersL.et al (2019). Dry Amorphisation of Mangiferin, a Poorly Water-Soluble Compound, Using Mesoporous Silica. Eur. J. Pharm. Biopharm.141, 172179. 10.1016/j.ejpb.2019.05.026

  • 7

    BaoL.HaqueA.JacksonK.HazariS.MorozK.JetlyR.et al (2011). Increased Expression of P-Glycoprotein Is Associated with Doxorubicin Chemoresistance in the Metastatic 4T1 Breast Cancer Model. Am. J. Pathol.178 (2), 838852. 10.1016/j.ajpath.2010.10.029

  • 8

    BatlleR.AndrésE.GonzalezL.LlonchE.IgeaA.Gutierrez-PratN.et al (2019). Regulation of Tumor Angiogenesis and Mesenchymal-Endothelial Transition by P38α through TGF-β and JNK Signaling. Nat. Commun.10 (1), 3071. 10.1038/s41467-019-10946-y

  • 9

    BeattyG.PatersonY. (2001). IFN-gamma-dependent Inhibition of Tumor Angiogenesis by Tumor-Infiltrating CD4+ T Cells Requires Tumor Responsiveness to IFN-Gamma. J. Immunol.166 (4), 22762282. 10.4049/jimmunol.166.4.2276

  • 10

    BelisarioD. C.KopeckaJ.PasinoM.AkmanM.De SmaeleE.DonadelliM.et al (2020). Hypoxia Dictates Metabolic Rewiring of Tumors: Implications for Chemoresistance. Cells9 (12), 2598. 10.3390/cells9122598

  • 11

    BhatiaA.KumarY. (2016). Cancer Stem Cells and Tumor Immunoediting: Putting Two and Two Together. Expert Rev. Clin. Immunol.12 (6), 605607. 10.1586/1744666X.2016.1159133

  • 12

    BrayF.FerlayJ.SoerjomataramI.SiegelR. L.TorreL. A.JemalA. (2018). Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin.68 (6), 394424. 10.3322/caac.21492

  • 13

    ChalkiadakiA.GuarenteL. (2015). The Multifaceted Functions of Sirtuins in Cancer. Nat. Rev. Cancer15 (10), 608624. 10.1038/nrc3985

  • 14

    ChenS. T.DouJ.TempleR.AgarwalR.WuK. M.WalkerS. (2008). New Therapies from Old Medicines. Nat. Biotechnol.26 (10), 10771083. 10.1038/nbt1008-1077

  • 15

    CornK. C.WindhamM. A.RafatM. (2020). Lipids in the Tumor Microenvironment: From Cancer Progression to Treatment. Prog. Lipid Res.80, 101055. 10.1016/j.plipres.2020.101055

  • 16

    Delgado-HernándezR.Hernández-BalmasedaI.Rodeiro-GuerraI.Cesar Rodriguez GonzalezJ.De WeverO.LogieE.et al (2020). Anti-angiogenic Effects of Mangiferin and Mechanism of Action in Metastatic Melanoma. Melanoma Res.30 (1), 3951. 10.1097/CMR.0000000000000647

  • 17

    DhirM.SassonA. R. (2016). Surgical Management of Liver Metastases from Colorectal Cancer. J. Oncol. Pract.12 (1), 3339. 10.1200/JOP.2015.009407

  • 18

    DrabkinD.AustinJ. (1935). Spectrophotometric Studies. II. Preparations from Washed Blood Cells: Nitric Oxide Hemoglobin and Sulfhemoglobin. J. Biol. Chem.112 (51), 19351936. 10.1016/s0021-9258(18)74968-5

  • 19

    DuJ.YangM.ChenS.LiD.ChangZ.DongZ. (2016). PDK1 Promotes Tumor Growth and Metastasis in a Spontaneous Breast Cancer Model. Oncogene35 (25), 33143323. 10.1038/onc.2015.393

  • 20

    du Plessis-StomanD.du PreezJ.van de VenterM. (2011). Combination Treatment with Oxaliplatin and Mangiferin Causes Increased Apoptosis and Downregulation of NFκB in Cancer Cell Lines. Afr. J. Tradit Complement. Altern. Med.8 (2), 177184. 10.4314/ajtcam.v8i2.63206

  • 21

    DunningM. J.SmithM. L.RitchieM. E.TavaréS. (2007). Beadarray: R Classes and Methods for Illumina Bead-Based Data. Bioinformatics23 (16), 21832184. 10.1093/bioinformatics/btm311

  • 22

    DupuyF.TabarièsS.AndrzejewskiS.DongZ.BlagihJ.AnnisM. G.et al (2015). PDK1-Dependent Metabolic Reprogramming Dictates Metastatic Potential in Breast Cancer. Cell Metab22 (4), 577589. 10.1016/j.cmet.2015.08.007

  • 23

    EvansJ. P.SuttonP. A.WiniarskiB. K.FenwickS. W.MalikH. Z.VimalachandranD.et al (2016). From Mice to Men: Murine Models of Colorectal Cancer for Use in Translational Research. Crit. Rev. Oncol. Hematol.98, 94105. 10.1016/j.critrevonc.2015.10.009

  • 24

    FangY.ZhangD.HuT.ZhaoH.ZhaoX.LouZ.et al (2019). KMT2A Histone Methyltransferase Contributes to Colorectal Cancer Development by Promoting Cathepsin Z Transcriptional Activation. Cancer Med.8 (7), 35443552. 10.1002/cam4.2226

  • 25

    GaálZ.CsernochL. (2020). Impact of Sirtuin Enzymes on the Altered Metabolic Phenotype of Malignantly Transformed Cells. Front. Oncol.10, 45. 10.3389/fonc.2020.00045

  • 26

    García-RiveraD.DelgadoR.BougarneN.HaegemanG.BergheW. V. (2011). Gallic Acid Indanone and Mangiferin Xanthone Are strong Determinants of Immunosuppressive Anti-tumour Effects of Mangifera Indica L. Bark in MDA-Mb231 Breast Cancer Cells. Cancer Lett.305 (1), 2131. 10.1016/j.canlet.2011.02.011

  • 27

    GilE.HernándezR. G.NuevasL.HoogmartensJ. (2014). Determination of Mangiferin in Mangifera Indica L. Stem Bark Extract (Vimang®) and Pharmaceuticals by Liquid Chromatography. Emirates J. Food Agric.26, 592601. 10.9755/ejfa.v26i7.18186

  • 28

    GirnunG. D.ChenL.SilvaggiJ.DrapkinR.ChirieacL. R.PaderaR. F.et al (2008). Regression of Drug-Resistant Lung Cancer by the Combination of Rosiglitazone and Carboplatin. Clin. Cancer Res.14 (20), 64786486. 10.1158/1078-0432.CCR-08-1128

  • 29

    GirnunG. D.NaseriE.VafaiS. B.QuL.SzwayaJ. D.BronsonR.et al (2007). Synergy between PPARgamma Ligands and Platinum-Based Drugs in Cancer. Cancer Cell11 (5), 395406. 10.1016/j.ccr.2007.02.025

  • 30

    GouQ.GongX.JinJ.ShiJ.HouY. (2017). Peroxisome Proliferator-Activated Receptors (PPARs) Are Potential Drug Targets for Cancer Therapy. Oncotarget8 (36), 6070460709. 10.18632/oncotarget.19610

  • 31

    Grochot-PrzeczekA.KotlinowskiJ.KozakowskaM.StarowiczK.JagodzinskaJ.StachurskaA.et al (2014). Heme Oxygenase-1 Is Required for Angiogenic Function of Bone Marrow-Derived Progenitor Cells: Role in Therapeutic Revascularization. Antioxid. Redox Signal.20 (11), 16771692. 10.1089/ars.2013.5426

  • 32

    GuoF.HuangC.LiaoX.WangY.HeY.FengR.et al (2011). Beneficial Effects of Mangiferin on Hyperlipidemia in High-Fat-Fed Hamsters. Mol. Nutr. Food Res.55 (12), 18091818. 10.1002/mnfr.201100392

  • 33

    HanH.ChoJ. W.LeeS.YunA.KimH.BaeD.et al (2018). TRRUST V2: an Expanded Reference Database of Human and Mouse Transcriptional Regulatory Interactions. Nucleic Acids Res.46 (D1), D380D386. 10.1093/nar/gkx1013

  • 34

    HassanniaB.WiernickiB.IngoldI.QuF.Van HerckS.TyurinaY. Y.et al (2018). Nano-targeted Induction of Dual Ferroptotic Mechanisms Eradicates High-Risk Neuroblastoma. J. Clin. Invest.128 (8), 33413355. 10.1172/JCI99032

  • 35

    HeW.YouY.DuS.LeiT.WangH.LiX.et al (2019). Anti-neoplastic Effect of Mangiferin on Human Ovarian Adenocarcinoma OVCAR8 Cells via the Regulation of YAP. Oncol. Lett.17 (1), 10081018. 10.3892/ol.2018.9708

  • 36

    HeideggerI.KernJ.OferP.KlockerH.MassonerP. (2014). Oncogenic Functions of IGF1R and INSR in Prostate Cancer Include Enhanced Tumor Growth, Cell Migration and Angiogenesis. Oncotarget5 (9), 27232735. 10.18632/oncotarget.1884

  • 37

    Hernandez-BalmasedaI.GuerraI. R.DeclerckK.Herrera IsidrónJ. A.Pérez-NovoC.Van CampG.et al (2021). Marine Seagrass Extract of Thalassia Testudinum Suppresses Colorectal Tumor Growth, Motility and Angiogenesis by Autophagic Stress and Immunogenic Cell Death Pathways. Mar. Drugs19 (2), 52. 10.3390/md19020052

  • 38

    HwangD. H.KimJ. A.LeeJ. Y. (2016). Mechanisms for the Activation of Toll-like Receptor 2/4 by Saturated Fatty Acids and Inhibition by Docosahexaenoic Acid. Eur. J. Pharmacol.785, 2435. 10.1016/j.ejphar.2016.04.024

  • 39

    ImranM.ArshadM. S.ButtM. S.KwonJ. H.ArshadM. U.SultanM. T. (2017). Mangiferin: a Natural Miracle Bioactive Compound against Lifestyle Related Disorders. Lipids Health Dis.16 (1), 84. 10.1186/s12944-017-0449-y

  • 40

    IsideC.ScafuroM.NebbiosoA.AltucciL. (2020). SIRT1 Activation by Natural Phytochemicals: An Overview. Front. Pharmacol.11, 1225. 10.3389/fphar.2020.01225

  • 41

    JuK. D.LimJ. W.KimH. (2017). Peroxisome Proliferator-Activated Receptor-Gamma Inhibits the Activation of STAT3 in Cerulein-Stimulated Pancreatic Acinar Cells. J. Cancer Prev.22 (3), 189194. 10.15430/JCP.2017.22.3.189

  • 42

    KadochiY.MoriS.Fujiwara-TaniR.LuoY.NishiguchiY.KishiS.et al (2017). Remodeling of Energy Metabolism by a Ketone Body and Medium-Chain Fatty Acid Suppressed the Proliferation of CT26 Mouse colon Cancer Cells. Oncol. Lett.14 (1), 673680. 10.3892/ol.2017.6195

  • 43

    Karaman MayackB.SipplW.Ntie-KangF. (2020). Natural Products as Modulators of Sirtuins. Molecules25 (14), 3287. 10.3390/molecules25143287

  • 44

    KastanaP.ZahraF. T.NtenekouD.Katraki-PavlouS.BeisD.LionakisM. S.et al (2019). Matrigel Plug Assay for In Vivo Evaluation of Angiogenesis. Methods Mol. Biol.1952, 219232. 10.1007/978-1-4939-9133-4_18

  • 45

    KhuranaR. K.KaurR.LohanS.SinghK. K.SinghB. (2016). Mangiferin: a Promising Anticancer Bioactive. Pharm. Pat Anal.5 (3), 169181. 10.4155/ppa-2016-0003

  • 46

    KobeltD.WaltherW.SteinU. S. (2021). Real-Time Cell Migration Monitoring to Analyze Drug Synergism in the Scratch Assay Using the IncuCyte System. Methods Mol. Biol.2294, 133142. 10.1007/978-1-0716-1350-4_9

  • 47

    KoundourosN.PoulogiannisG. (2020). Reprogramming of Fatty Acid Metabolism in Cancer. Br. J. Cancer122 (1), 422. 10.1038/s41416-019-0650-z

  • 48

    KumarN.Perez-NovoC.ShawP.LogieE.Privat-MaldonadoA.DewildeS.et al (2021). Physical Plasma-Derived Oxidants Sensitize Pancreatic Cancer Cells to Ferroptotic Cell Death. Free Radic. Biol. Med.166, 187200. 10.1016/j.freeradbiomed.2021.02.026

  • 49

    LederleW.HartensteinB.MeidesA.KunzelmannH.WerbZ.AngelP.et al (2010). MMP13 as a Stromal Mediator in Controlling Persistent Angiogenesis in Skin Carcinoma. Carcinogenesis31 (7), 11751184. 10.1093/carcin/bgp248

  • 50

    LeiroJ.ArranzJ. A.YáñezM.UbeiraF. M.SanmartínM. L.OralloF. (2004). Expression Profiles of Genes Involved in the Mouse Nuclear Factor-Kappa B Signal Transduction Pathway Are Modulated by Mangiferin. Int. Immunopharmacol4 (6), 763778. 10.1016/j.intimp.2004.03.002

  • 51

    LiJ.LiuM.YuH.WangW.HanL.ChenQ.et al (2018). Mangiferin Improves Hepatic Lipid Metabolism Mainly through its Metabolite-Norathyriol by Modulating SIRT-1/AMPK/SREBP-1c Signaling. Front. Pharmacol.9, 201. 10.3389/fphar.2018.00201

  • 52

    LiJ.ZhangH.DaiZ. (2021). Cancer Treatment with the Ketogenic Diet: A Systematic Review and Meta-Analysis of Animal Studies. Front. Nutr.8, 594408. 10.3389/fnut.2021.594408

  • 53

    LiuY.ZhangJ.ChenY.SohelH.KeX.ChenJ.et al (2020). The Correlation and Role Analysis of COL4A1 and COL4A2 in Hepatocarcinogenesis. Aging (Albany NY)12 (1), 204223. 10.18632/aging.102610

  • 54

    LiuZ.ApontesP.FomenkoE. V.ChiN.SchusterV. L.KurlandI. J.et al (2018). Mangiferin Accelerates Glycolysis and Enhances Mitochondrial Bioenergetics. Int. J. Mol. Sci.19 (1), 201. 10.3390/ijms19010201

  • 55

    MaY.TemkinS. M.HawkridgeA. M.GuoC.WangW.WangX. Y.et al (2018). Fatty Acid Oxidation: An Emerging Facet of Metabolic Transformation in Cancer. Cancer Lett.435, 92100. 10.1016/j.canlet.2018.08.006

  • 56

    Mammadova-BachE.Gil-PulidoJ.SarukhanyanE.BurkardP.ShityakovS.SchonhartC.et al (2020). Platelet Glycoprotein VI Promotes Metastasis through Interaction with Cancer Cell-Derived Galectin-3. Blood135 (14), 11461160. 10.1182/blood.2019002649

  • 57

    MayerR. J. (2009). Targeted Therapy for Advanced Colorectal Cancer-Mmore Is Not Always Better. N. Engl. J. Med.360 (6), 623625. 10.1056/NEJMe0809343

  • 58

    MichalikL.DesvergneB.WahliW. (2004). Peroxisome-proliferator-activated Receptors and Cancers: Complex Stories. Nat. Rev. Cancer4 (1), 6170. 10.1038/nrc1254

  • 59

    MihalacheA.RogoveanuI. (2014). Angiogenesis Factors Involved in the Pathogenesis of Colorectal Cancer. Curr. Health Sci. J.40 (1), 511. 10.12865/CHSJ.40.01.01

  • 60

    Ministerio de Salud Pública. Dirección de Registros Médicos y Estadísticas de SaludL. H.de CubaRepública. (2017). Anuario Estadístico de Salud 2016.

  • 61

    MoraesL. A.SpyridonM.KaiserW. J.JonesC. I.SageT.AthertonR. E.et al (2010). Non-genomic Effects of PPARgamma Ligands: Inhibition of GPVI-Stimulated Platelet Activation. J. Thromb. Haemost.8 (3), 577587. 10.1111/j.1538-7836.2009.03732.x

  • 62

    MorozkinaS. N.Nhung VuT. H.GeneralovaY. E.SnetkovP. P.UspenskayaM. V. (2021). Mangiferin as New Potential Anti-cancer Agent and Mangiferin-Integrated Polymer Systems-A Novel Research Direction. Biomolecules11 (1), 79. 10.3390/biom11010079

  • 63

    NaL.ZhangQ.JiangS.DuS.ZhangW.LiY.et al (2015). Mangiferin Supplementation Improves Serum Lipid Profiles in Overweight Patients with Hyperlipidemia: a Double-Blind Randomized Controlled Trial. Sci. Rep.5, 10344. 10.1038/srep10344

  • 64

    NagarajanS. R.ButlerL. M.HoyA. J. (2021). The Diversity and Breadth of Cancer Cell Fatty Acid Metabolism. Cancer Metab.9 (1), 2. 10.1186/s40170-020-00237-2

  • 65

    NakamuraK.TonouchiH.SasayamaA.AshidaK. (2018). A Ketogenic Formula Prevents Tumor Progression and Cancer Cachexia by Attenuating Systemic Inflammation in Colon 26 Tumor-Bearing Mice. Nutrients10 (2), 206. 10.3390/nu10020206

  • 66

    NewmanD. J.CraggG. M. (2020). Natural Products as Sources of New Drugs over the Nearly Four Decades from 01/1981 to 09/2019. J. Nat. Prod.83 (3), 770803. 10.1021/acs.jnatprod.9b01285

  • 67

    NiuY.LiS.NaL.FengR.LiuL.LiY.et al (2012). Mangiferin Decreases Plasma Free Fatty Acids through Promoting its Catabolism in Liver by Activation of AMPK. PLoS One7 (1), e30782. 10.1371/journal.pone.0030782

  • 68

    Núñez SellesA. J.DagliaM.RastrelliL. (2016). The Potential Role of Mangiferin in Cancer Treatment through its Immunomodulatory, Anti-angiogenic, Apoptopic, and Gene Regulatory Effects. Biofactors42 (5), 475491. 10.1002/biof.1299

  • 69

    Núñez SellésA. J.Vélez CastroH. T.Agüero-AgüeroJ.González-GonzálezJ.NaddeoF.De SimoneF.et al (2002). Isolation and Quantitative Analysis of Phenolic Antioxidants, Free Sugars, and Polyols from Mango (Mangifera Indica L.) Stem Bark Aqueous Decoction Used in Cuba as a Nutritional Supplement. J. Agric. Food Chem.50 (4), 762766. 10.1021/jf011064b

  • 70

    PanigrahyD.HuangS.KieranM. W.KaipainenA. (2005). PPARgamma as a Therapeutic Target for Tumor Angiogenesis and Metastasis. Cancer Biol. Ther.4 (7), 687693. 10.4161/cbt.4.7.2014

  • 71

    PetrelliF.TomaselloG.BorgonovoK.GhidiniM.TuratiL.DalleraP.et al (2017). Prognostic Survival Associated with Left-Sided vs Right-Sided Colon Cancer: A Systematic Review and Meta-Analysis. JAMA Oncol.3 (2), 211219. 10.1001/jamaoncol.2016.4227

  • 72

    PoulsenLl.SiersbækM.MandrupS. (2012). PPARs: Fatty Acid Sensors Controlling Metabolism. Semin. Cell Dev Biol23 (6), 631639. 10.1016/j.semcdb.2012.01.003

  • 73

    PritchardJ. R.LauffenburgerD. A.HemannM. T. (2012). Understanding Resistance to Combination Chemotherapy. Drug Resist. Updat15 (5-6), 249257. 10.1016/j.drup.2012.10.003

  • 74

    PuchalskaP.CrawfordP. A. (2017). Multi-dimensional Roles of Ketone Bodies in Fuel Metabolism, Signaling, and Therapeutics. Cell Metab25 (2), 262284. 10.1016/j.cmet.2016.12.022

  • 75

    RitchieM. E.PhipsonB.WuD.HuY.LawC. W.ShiW.et al (2015). Limma powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res.43 (7), e47. 10.1093/nar/gkv007

  • 76

    SadanandamA.RosenbaughE. G.SinghS.VarneyM.SinghR. K. (2010). Semaphorin 5A Promotes Angiogenesis by Increasing Endothelial Cell Proliferation, Migration, and Decreasing Apoptosis. Microvasc. Res.79 (1), 19. 10.1016/j.mvr.2009.10.005

  • 77

    SarkarA.SreenivasanY.RameshG. T.MannaS. K. (2004). Beta-D-Glucoside Suppresses Tumor Necrosis Factor-Induced Activation of Nuclear Transcription Factor kappaB but Potentiates Apoptosis. J. Biol. Chem.279 (32), 3376833781. 10.1074/jbc.M403424200

  • 78

    SaxenaK.JollyM. K.BalamuruganK. (2020). Hypoxia, Partial EMT and Collective Migration: Emerging Culprits in Metastasis. Transl Oncol.13 (11), 100845. 10.1016/j.tranon.2020.100845

  • 79

    ShakedY. (2016). Balancing Efficacy of and Host Immune Responses to Cancer Therapy: the Yin and Yang Effects. Nat. Rev. Clin. Oncol.13 (10), 611626. 10.1038/nrclinonc.2016.57

  • 80

    SinghA. K.RajV.KeshariA. K.RaiA.KumarP.RawatA.et al (2018). Isolated Mangiferin and Naringenin Exert Antidiabetic Effect via PPARγ/GLUT4 Dual Agonistic Action with strong Metabolic Regulation. Chem. Biol. Interact280, 3344. 10.1016/j.cbi.2017.12.007

  • 81

    SinghP.JenkinsL. M.HorstB.AlersV.PradhanS.KaurP.et al (2018). Inhibin Is a Novel Paracrine Factor for Tumor Angiogenesis and Metastasis. Cancer Res.78 (11), 29782989. 10.1158/0008-5472.CAN-17-2316

  • 82

    StemmlerM. P.EcclesR. L.BrabletzS.BrabletzT. (2019). Non-redundant Functions of EMT Transcription Factors. Nat. Cell Biol21 (1), 102112. 10.1038/s41556-018-0196-y

  • 83

    TanH. Y.WangN.LiS.HongM.GuoW.ManK.et al (2018). Repression of WT1-Mediated LEF1 Transcription by Mangiferin Governs β-Catenin-Independent Wnt Signalling Inactivation in Hepatocellular Carcinoma. Cell Physiol Biochem47 (5), 18191834. 10.1159/000491063

  • 84

    TengM. W.GalonJ.FridmanW. H.SmythM. J. (2015). From Mice to Humans: Developments in Cancer Immunoediting. J. Clin. Invest.125 (9), 33383346. 10.1172/JCI80004

  • 85

    TingW. C.ChenL. M.PaoJ. B.YangY. P.YouB. J.ChangT. Y.et al (2013). Common Genetic Variants in Wnt Signaling Pathway Genes as Potential Prognostic Biomarkers for Colorectal Cancer. PLoS One8 (2), e56196. 10.1371/journal.pone.0056196

  • 86

    TorreL. A.SiegelR. L.WardE. M.JemalA. (2016). Global Cancer Incidence and Mortality Rates and Trends--An Update. Cancer Epidemiol. Biomarkers Prev.25 (1), 1627. 10.1158/1055-9965.EPI-15-0578

  • 87

    ValléeA.LecarpentierY. (2018). Crosstalk between Peroxisome Proliferator-Activated Receptor Gamma and the Canonical WNT/β-Catenin Pathway in Chronic Inflammation and Oxidative Stress during Carcinogenesis. Front. Immunol.9, 745. 10.3389/fimmu.2018.00745

  • 88

    VatandoustS.PriceT. J.KarapetisC. S. (2015). Colorectal Cancer: Metastases to a Single Organ. World J. Gastroenterol.21 (41), 1176711776. 10.3748/wjg.v21.i41.11767

  • 89

    VuT.DattaP. K. (2017). Regulation of EMT in Colorectal Cancer: A Culprit in Metastasis. Cancers (Basel)9 (12), 171. 10.3390/cancers9120171

  • 90

    WangJ.ZhuC. (2016). Anticoagulation in Combination with Antiangiogenesis and Chemotherapy for Cancer Patients: Evidence and Hypothesis. Onco Targets Ther.9, 47374746. 10.2147/OTT.S103184

  • 91

    WangL. H.YangX. Y.ZhangX.HuangJ.HouJ.LiJ.et al (2004). Transcriptional Inactivation of STAT3 by PPARgamma Suppresses IL-6-responsive Multiple Myeloma Cells. Immunity20 (2), 205218. 10.1016/s1074-7613(04)00030-5

  • 92

    WangT.JinH.HuJ.LiX.RuanH.XuH.et al (2020). COL4A1 Promotes the Growth and Metastasis of Hepatocellular Carcinoma Cells by Activating FAK-Src Signaling. J. Exp. Clin. Cancer Res.39 (1), 148. 10.1186/s13046-020-01650-7

  • 93

    WangZ.DabrosinC.YinX.FusterM. M.ArreolaA.RathmellW. K.et al (2015). Broad Targeting of Angiogenesis for Cancer Prevention and Therapy. Semin. Cancer Biol.35 Suppl (Suppl. l), S224S243. 10.1016/j.semcancer.2015.01.001

  • 94

    WoolfE. C.CurleyK. L.LiuQ.TurnerG. H.CharltonJ. A.PreulM. C.et al (2015). The Ketogenic Diet Alters the Hypoxic Response and Affects Expression of Proteins Associated with Angiogenesis, Invasive Potential and Vascular Permeability in a Mouse Glioma Model. PLoS One10 (6), e0130357. 10.1371/journal.pone.0130357

  • 95

    World Health Organization (2019). WHO Global Report on Traditional and Complementary Medicine 2019. Geneva: World Health Organization.

  • 96

    YamadaT.OshimaT.YoshiharaK.TamuraS.KanazawaA.InagakiD.et al (2010). Overexpression of MMP-13 Gene in Colorectal Cancer with Liver Metastasis. Anticancer Res.30 (7), 26932699.

  • 97

    YuanZ. H.LiuT.WangH.XueL. X.WangJ. J. (2021). Fatty Acids Metabolism: The Bridge between Ferroptosis and Ionizing Radiation. Front Cell Dev Biol9, 675617. 10.3389/fcell.2021.675617

  • 98

    ZhangJ.CaoR.ZhangY.JiaT.CaoY.WahlbergE. (2009). Differential Roles of PDGFR-Alpha and PDGFR-Beta in Angiogenesis and Vessel Stability. FASEB J.23 (1), 153163. 10.1096/fj.08-113860

  • 99

    ZhangQ.KongX.YuanH.GuanH.LiY.NiuY. (2019). Mangiferin Improved Palmitate-Induced-Insulin Resistance by Promoting Free Fatty Acid Metabolism in HepG2 and C2C12 Cells via PPARα: Mangiferin Improved Insulin Resistance. J. Diabetes Res.2019, 2052675. 10.1155/2019/2052675

  • 100

    ZhouJ.ZhangS.XueJ.AveryJ.WuJ.LindS. E.et al (2012). Activation of Peroxisome Proliferator-Activated Receptor α (PPARα) Suppresses Hypoxia-Inducible Factor-1α (HIF-1α) Signaling in Cancer Cells. J. Biol. Chem.287 (42), 3516135169. 10.1074/jbc.m112.367367

  • 101

    ZhouY.ZhouB.PacheL.ChangM.KhodabakhshiA. H.TanaseichukO.et al (2019). Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat. Commun.10 (1), 15231533. 10.1038/s41467-019-09234-6

  • 102

    ZhuX.ChenH. H.ikGaoC. Y.ZhangX. X.JiangJ. X.ZhangY.et al (2020). Energy Metabolism in Cancer Stem Cells. World J. Stem Cell12 (6), 448461. 10.4252/wjsc.v12.i6.448

  • 103

    Zorde KhvalevskyE.GabaiR.RachmutI. H.HorwitzE.BrunschwigZ.OrbachA.et al (2013). Mutant KRAS Is a Druggable Target for Pancreatic Cancer. Proc. Natl. Acad. Sci. U S A.110 (51), 2072320728. 10.1073/pnas.1314307110

  • 104

    ZouB.WangH.LiuY.QiP.LeiT.SunM.et al (2017). Mangiferin Induces Apoptosis in Human Ovarian Adenocarcinoma OVCAR3 Cells via the Regulation of Notch3. Oncol. Rep.38 (3), 14311441. 10.3892/or.2017.5814

Summary

Keywords

Mangiferin (PubChem CID: 5081647), antitumor, antiangiogenesis, antimetastasic, colon carcinoma, mitochondrial metabolism

Citation

Rodriguez-Gonzalez JC, Hernández-Balmaseda I, Declerck K, Pérez-Novo C, Logie E, Theys C, Jakubek P, Quiñones-Maza OL, Dantas-Cassali G, Carlos dos Reis D, Van Camp G, Lopes Paz MT, Rodeiro-Guerra I, Delgado-Hernández R and Vanden Berghe W (2021) Antiproliferative, Antiangiogenic, and Antimetastatic Therapy Response by Mangiferin in a Syngeneic Immunocompetent Colorectal Cancer Mouse Model Involves Changes in Mitochondrial Energy Metabolism. Front. Pharmacol. 12:670167. doi: 10.3389/fphar.2021.670167

Received

23 February 2021

Accepted

03 November 2021

Published

03 December 2021

Volume

12 - 2021

Edited by

Cristiana Perrotta, University of Milan, Italy

Reviewed by

Seung Ho Shin, Gyeongsang National University, South Korea

Olga Krysko, Ghent University, Belgium

Updates

Copyright

*Correspondence: René Delgado-Hernández, , ; Wim Vanden Berghe,

† These authors have contributed equally to this work and share first authorship

This article was submitted to Pharmacology of Anti-Cancer Drugs, a section of the journal Frontiers in Pharmacology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics