Skip to main content
Download PDF

Metabolic consequences of perioperative oral carbohydrates in breast cancer patients — an explorative study

BMC Cancer volume 19, Article number: 1183 (2019) Cite this article

Abstract

Background

The metabolic consequences of preoperative carbohydrate load in breast cancer patients are not known. The present explorative study investigated the systemic and tumor metabolic changes after preoperative per-oral carbohydrate load and their influence on tumor characteristics and survival.

Methods

The study setting was on university hospital level with primary and secondary care functions in south-west Norway. Serum and tumor tissue were sampled from a population-based cohort of 60 patients with operable breast cancer who were randomized to either per-oral carbohydrate load (preOp™; n = 25) or standard pre-operative fasting (n = 35) before surgery. Magnetic resonance (MR) metabolomics was performed on serum samples from all patients and high-resolution magic angle spinning (HR-MAS) MR analysis on 13 tumor samples available from the fasting group and 16 tumor samples from the carbohydrate group.

Results

Fourteen of 28 metabolites were differently expressed between fasting and carbohydrate groups. Partial least squares discriminant analysis showed a significant difference in the metabolic profile between the fasting and carbohydrate groups, compatible with the endocrine effects of insulin (i.e., increased serum-lactate and pyruvate and decreased ketone bodies and amino acids in the carbohydrate group). Among ER-positive tumors (n = 18), glutathione was significantly elevated in the carbohydrate group compared to the fasting group (p = 0.002), with a positive correlation between preoperative S-insulin levels and the glutathione content in tumors (r = 0.680; p = 0.002). In all tumors (n = 29), glutamate was increased in tumors with high proliferation (t-test; p = 0.009), independent of intervention group. Moreover, there was a positive correlation between tumor size and proliferation markers in the carbohydrate group only. Patients with ER-positive / T2 tumors and high tumor glutathione (≥1.09), high S-lactate (≥56.9), and high S-pyruvate (≥12.5) had inferior clinical outcomes regarding relapse-free survival, breast cancer-specific survival, and overall survival. Moreover, Integrated Pathway Analysis (IPA) in serum revealed activation of five major anabolic metabolic networks contributing to proliferation and growth.

Conclusions

Preoperative carbohydrate load increases systemic levels of lactate and pyruvate and tumor levels of glutathione and glutamate in ER-positive patients. These biological changes may contribute to the inferior clinical outcomes observed in luminal T2 breast cancer patients.

Trial of registration

ClinicalTrials.gov; NCT03886389. Retrospectively registered March 22, 2019.

Peer Review reports

Background

Breast cancer is the most common female malignancy and one of the most frequent causes of death among women in the Western world [1]. Breast cancer incidence has more than doubled in the last 50 years, probably due to increased estrogen exposure and a change towards high levels of alimentary carbohydrates and fat [2, 3]. Even though breast cancer originates locally in the breast, circulating tumor cells (CTCs) may spread to the systemic circulation before and during surgery [4] and establish distant micrometastases [5]. These CTCs must thrive and survive attacks from the innate and adaptive immune system. Thus, tumor cells have to establish a favorable metabolism that can produce energy, protection mechanisms, and the necessary biomass to survive the journey from the breast tumor to remote locations, including transformation into dormancy [6]. The luminal breast cancer subtype, which express estrogen receptor (ER) and/or progesterone receptor (PR) in the tumor cells, comprise the largest subgroup, accounting for approximately 75% of all breast cancers. Endocrine resistance in this subtype can creates micrometastases that escape anti-estrogen therapy and can hibernate for many years before they become clinically overt [7]. The molecular features underlying these cellular characteristics are driven by hallmarks of cancer [8], including changes in cellular energetics and metabolism, followed by a vast number of necessary metabolic modifications to strengthen the metabolic needs of breast cancer cells [9]. A well-known cellular characteristic of tumor cells is increased glucose consumption and glycolysis towards lactate despite the presence of oxygen, a feature called ‘the Warburg effect’ [10, 11]. This metabolic switch includes the production of ribose for DNA synthesis and allowing amino acids to be a source for ATP production [12]. Furthermore, the Warburg effect extends to increased choline metabolism for cell membrane synthesis and increased amino acid turnover for protein synthesis [10, 13].

Even though much is known about metabolism in breast cancer cells [14], little is known about the influence of carbohydrate loading in the early recovery after surgery (ERAS) program [15] on peri-operative metabolism in the systemic circulation and locally in the breast tumor. We recently conducted a randomized controlled trial (RCT) in which operable breast cancer patients were treated with either two oral loads of enriched carbohydrate solution or a standard fasting procedure comprising free drinking of tap water before surgery [16]. In this study, luminal breast cancer patients, who received oral pre-operative carbohydrates, had a higher tumor proliferation and an adverse survival. The goal of the present paper, using the same patients, was to further explore the metabolic differences in serum and the tumor. Based on our previous findings, we hypothesize that the metabolic changes after carbohydrate loading will correlate with proliferation and outcome in patients with ER positive tumors. Also, we also wanted to study whether such metabolic alterations correlate with other tumor characteristics or translate into differences in clinical outcome.

Methods

Ethics statement

This paper is an explorative study based upon a recently published randomized controlled trial (RCT) approved by the Regional Ethics Committee in Western Norway (#2015/1445) and was retrospectively registered at Clinicaltrials.gov (NCT03886389).

Patients

Details on these patients have been described previously [16]. In short, between 12 May 2009 and 23 June 2010 a population-based cohort of 61 operable breast cancer patients (Stage I and II) were randomized into an intervention group receiving preoperative per-oral carbohydrate loading (n = 26) or a control group (n = 35) receiving the standard preoperative fasting protocol.

The patients in the carbohydrate group drank 200 mL pre-Op™ (Nutricia, the Netherlands). This non-carbonated carbohydrate enriched drink contained 100 kCal per bottle containing 4.2 g (2.1%) glucose and 20 g (10%) polysaccharides. A loading dose of two bottles pre-Op™ were given 18 h before surgery (i.e. the evening before surgery) and another 2 bottles were administered 2–4 h before surgery (i.e. the morning of the operation day). In contrast, the control group practiced the standard fasting procedure with free intake of tap water 12–14 h before surgery. From this cohort, patients with available fresh frozen tissue and serum samples were included in the present study (Fig. 1). The patient characteristics are given in Table 1.

Fig. 1
figure 1

Flowchart of study participants

Full size image
Table 1 Clinical variables in the randomized groups
Full size table

Blood sampling

Blood samples were drawn immediately before surgery. In total three serum gel tubes and one EDTA plasma tube were drawn in this study. One serum gel tube and one EDTA plasma tube were delivered within an hour to the department of medical biochemistry for standard analysis. For metabolomics analyses, two serum gel tubes were centrifuged within one hour at 4 °C, 2500 x g in 10 min. After centrifugation, the serum of the two tubes were mixed and a minimum of 1.1 mL serum were sent for analyses in Haukeland University Hospital, Bergen, Norway, the rest of the serum were stored in 1 mL cryotubes at − 80 °C in the biobank at Stavanger University Hospital, Stavanger, Norway.

Tumor tissue sampling

Immediately after removal of the surgical specimen from the systemic circulation, it was transported to the Department of Pathology for further sampling. To avoid necrotic areas, cancerous tissue from the invasive front of the tumor (i.e. tumor periphery) was immediately snap-frozen in liquid nitrogen and stored at − 80 °C until assayed for tissue metabolomics. Before HR-MAS analysis, tissues from all of the patients were analyzed consecutively for histopathology and immunohistochemistry as described preciously [16].

Serum hormone and protein analyses

Serum was transported to the Hormone Laboratory, Haukeland University Hospital, Bergen, Norway. Insulin, insulin c-peptide, insulin growth factor 1 (IGF-1), and insulin growth factor binding protein 3 (IGFBP-3) were measured by the IMMULITE 2000 two-site chemiluminescent immunometric assay (Siemens Medical Solutions Diagnostics).

Serum metabolomics analyses

A separate aliquot of serum was transported to the MR Core Facility at NTNU, Trondheim, Norway for metabolomics analyses. Thawed samples (100 μL) were mixed with bacteriostatic buffer (100 μL; pH 7.4, 0.075 mM Na2HPO4, 5 mM NaN3, 5 mM TSP), transferred to 3-mm NMR tubes, and stored at 5 °C until analysis (< 15 h). The MR analysis was performed using a Bruker Avance III Ultrashielded Plus 600 MHz spectrometer (Bruker Biospin GmbH, Germany) equipped with a 5 mm QCI Cryoprobe with integrated, cooled pre-amplifiers for 1H, 2H, and 13C. Experiments were fully automated using the SampleJet™ in combination with Icon-NMR in TopSpin 3.1 software (Bruker Biospin). One-dimensional 1H Nuclear Overhauser effect spectroscopy (NOESY) and Carr–Purcell–Meiboom–Gill (CPMG) spectra with water presaturation were acquired at 310.15 K. The spectra were Fourier transformed to 128 K after 0.3 Hz exponential line broadening and automatically phased and baseline-corrected. Spectra were further processed in Matlab 2013b (The Mathworks Inc., Natick, MA, USA). The CPMG spectral region between 0.1 and 4.2 ppm was selected for further processing. Chemical shifts were referenced to the left alanine peak at 1.47 ppm. Metabolites were identified based on previous assignment [17, 18]. Twenty-eight metabolites were identified as measurable and their areas calculated by integrating the area under the signal curve.

Breast tumor tissue metabolomics analyses

In the 29 patients with available tissue, the tumors were larger (45% vs. 9% pT2/3/4, p = 0.003), had a higher histological grade (52% vs. 18% grade 3, p = 0.022), were more often ER-negative (35% vs. 3%, p = 0.002), and had higher proliferation (59% vs. 27% PPH3-positive, p = 0.002) than those without tissue. Thus, we had a selection bias of larger, non-luminal and a more proliferative tumors into the present study compared to the original study [16]. Tissue was transported on dry ice to the MR Core Facility at NTNU, Trondheim, Norway, for metabolomics analyses. Tissue samples were prepared frozen on a metal plate bathed in liquid nitrogen to minimize tissue degradation. Biopsies (11.0 ± 2.3 mg) were cut to fit 30 μL disposable inserts (Bruker Biospin Corp, USA) filled with 3 μL D2O containing 25 mM formate. The insert containing the frozen sample was placed in a 4-mm diameter zirconium rotor (Bruker, Biospin GmbH, Germany) and kept at − 20 °C until analysis (< 8 h). Spin-echo spectra were acquired on a Bruker Avance DRX600 spectrometer with a 1H/13C magic angle spinning (MAS) probe with gradient (Bruker Biospin GmbH, Germany) using the following parameters: 5 KHz spin rate, 5 °C probe temperature, 5-min temperature acclimatization before shimming and spectral acquisition, CPMG pulse sequence (cpmgpr1d; Bruker) with 4 s water suppression prior to a 90° excitation pulse, total echo time 77 ms, 256 scans, and spectral width 20 ppm. Spectra were Fourier transformed into 64 K following 0.3 Hz line broadening. Phase correction was performed automatically for each spectrum using TopSpin 3.1.

Spectra were preprocessed in Matlab 2013b as follows [19]. The spectral region between 1.4–4.70 ppm, which contained the majority of the metabolite signals, was selected for further processing. Chemical shifts were referenced to the creatine peak at 3.03 ppm. The spectra were baseline-corrected using asymmetric least squares [20] with parameters λ = 1e7 and p = 0.0001, setting the lowest point in each spectrum to zero. Lipid peaks at 4.34–4.27, 4.19–4.14, 2.90–2.7, 2.31–2.18, 2.11–1.92, and 1.68–1.5, and ethanol at 3.67–3.62, were excluded. The resulting spectra were normalized to the total area to correct for differences in sample size and tumor cell content. Metabolite peak assignment was based on previous identification [21]. Twenty metabolites were identified as measurable, and the area under the signal curve in the preprocessed spectra was used to calculate their relative intensities. The metabolite integrals were log10 transformed to satisfy prerequisite assumptions of normality.

Endpoints

Proliferation differences between the carbohydrate and fasting groups were evaluated by Ki67 (< 15% or ≥ 15 and < 30% or ≥ 30%), mitotic activity index (MAI; < 10 or ≥ 10), and PPH3 (< 13 or ≥ 13). The metabolic response to preoperative oral carbohydrate loading was evaluated in serum (preoperative) by 1H NMR and in tumor tissue by HR-MAS MRS.

Univariate analysis

Metabolite differences between groups were assessed by student T-tests. Correlations between continuous variables were assessed by Pearson correlation. Categorical variables were compared by Chi square tests. P-values were considered significant when p < 0.05. When multiple variables were compared, the resulting p-value tables were corrected for multiple testing by the Benjamini-Hochberg method [22].

Multivariate analyses (serum and tissue)

Multivariate analyses were performed in R V.3.5 [23] using the package PLS [24] and MetaboAnalyst [25]. Metabolite values were auto-scaled (mean-centered and divided by variance) before multivariate analysis. Principal component analysis (PCA) was performed to evaluate the data sets for outliers. Partial least squares discriminant analysis (PLS-DA) was performed to explore differences in serum and tissue metabolic profiles between categories: carbohydrate loading vs fasting. Partial least squares (PLS) was used to find correlations between the tissue metabolic profile and variables (MAI, PPH3, Ki67, serum (S)-glucose, S-insulin, S-insulin c-peptide, S-IGFR, S-IGFPB3, S-estradiol). Metabolites were evaluated by Variable Importance in Projection (VIP) score. The VIP score is a measure of how important each variable was for creating the discrimination model. It is calculated as a weighted sum of squares of the PLS loadings, where the weights are based on the amount of y-variance explained in each dimension [26]. PLS and PLS-DA classification parameters were evaluated by ‘leave-one-out’ cross validation due to the limited sample numbers. Permutation testing was carried out as an additional model validation; sample classes or responses were shuffled, and the model rebuilt with the same numbers of latent variables as the original model. One thousand permutations were performed, and models were considered significant if the final accuracy (of classification models) or R2 (of regression models) were > 95% of the permuted accuracy values (p < 0.05).

Thresholds in survival analyses

Relapse-free survival (RFS) was defined as the time from surgery until a relapse from any site. Breast cancer-specific survival (BCSS) was defined as the time from surgery until death from breast cancer, whereas overall survival (OS) was until death from any cause. Receiver-operator characteristic (ROC) analysis identified optimal thresholds for the various continuous metabolite variables using relapse ‘Yes/No’ as the categorical variable (Table 8 in Appendix). The cut-off values obtained in RFS analysis were also used in the BCSS and OS analyses. In ER-negative patients, none of the explanatory variables with ROC-derived thresholds were significant for analysis of RFS, BCSS, or OS. Therefore, further analyses were limited to ER-positive patients. The ROC-obtained thresholds were confirmed with the minimal p-value/maximal Wald-value in a Cox model. In the multivariabel Cox analyses the ‘Forward Wald’ method was primarily used. In cases of an unstable model, a stepwise backward analysis was performed.

Metabolite set enrichment analysis and ingenuity pathway analysis (IPA)

Serum metabolite levels were uploaded to the Enrichment module of MetaboAnalyst to explore the pathways affected by the carbohydrate intervention. Pathway-associated metabolite sets with sets containing at least two metabolites were used. Pathways with p-values ≤0.05 (after FDR correction) were interpreted as significant. Serum metabolites with significantly different expression (p = 0.05) and their corresponding fold changes were imported into the Ingenuity Pathway Analysis (IPA) software (Ingenuity, Redwood City, USA) to explore which biological and molecular functions these metabolites were involved in and how these and their direct and indirect target molecules were connected, using the network function in IPA. Additionally, we examined if there were a direct or indirect connection between the top network and seven microRNAs related to tamoxifen resistance from our previous paper [27], using the grow function with a moderate or experimentally observed confidence level.

Results

Systemic metabolism

The results of the quantification of serum metabolites in the carbohydrate and fasting groups are given in Table 2. Fourteen out of 28 metabolites were significantly altered between the groups. PLS-DA revealed a significant difference in metabolic profiles between the two groups.; (one component, classification accuracy = 0.85; p < 0.001; Fig. 2a). The main increased markers were increased serum (S) lactate and S-pyruvate in the carbohydrate group (p < 0.0001; Fig. 2a and b). Among the patients in the fasting group, the levels of ketone bodies, such as S-acetate, S-acetoacetate, and S-3-hydroxybutyrate, were increased (Table 2). In addition, we observed increased S-N-acetylated groups, S-leucine, S-valine and S-isoleucine in the fasting group (all p < 0.05; Fig. 2b). We found positive correlations between tumor size and S-lactate (r = 0.344; p = 0.016) and tumor size and S-pyruvate (r = 0.370; p = 0.009).

Table 2 Serum metabolites with p-values from t-tests, fasting group versus carbohydrate (CH), for all patients and for the ER positive subset
Full size table
Fig. 2
figure 2

Partial Least Square Discriminant Analysis (PLS-DA) in serum. a Scores plot showing serum samples from the fasting group (green) and carbohydrate group (red). The carbohydrate and fasting groups have significantly different metabolic profiles as evidenced by permutation testing. b Variable Importance in Projection (VIP) scores showing the top 14 metabolites contributing to differences between the groups. The right column indicates increased (red) or decreased (green) metabolite in the indicated group

Full size image

In the carbohydrate group, there was a positive linear correlation between proliferation (Ki-67) and tumor size (r = 0.782, p = 0.038). When Ki-67, PPH3 and MAI were included in a forward and backward stepwise linear regression MAI was the only independent factor explaining increment in tumor size with a Beta = 0.530 (95%CI, 0.201 to 0.875) P = 0.009. In the fasting group, there was no correlation between tumor size and proliferation.’

Serum glucose and insulin responses

The mean fasting glucose and insulin values at admission were 5.4 mmol/L (95% CI 5.1 to 10.0) and 9.4 mIU (95% CI 6.8 to 32.5), respectively (normal ranges: glucose, 4.0 to 6.0 mmol/L; insulin, 6.0 to 27.0 mIU; c-peptide, 0.3 to 2.4 nmol/L). In the carbohydrate group, the mean preoperative insulin value was 35.6 mIU (26.7 to 106 mIU), compared to 9.1 (8.6 to 22 mIU) in the fasting group (student’s t-test p < 0.001). For C-peptide, the mean values in the carbohydrate and fasting groups were 2.10 nmol/L and 0.76 nmol/L, respectively (p < 0.001). We found significant univariate correlations between the serum concentrations of preoperative insulin (Table 3), Insulin C peptide (Table 9 in Appendix) IGFBP3 (Table 10 in Appendix), but not to IGF1 (Table 11 in Appendix). Multivariate analysis with leave-one-out cross-validation showed significant correlations between the serum metabolic profile and insulin (Cross-validated (CV) (R2 = 0.33, p < 0.001; Fig. 3a+b), Insulin C-peptide (CV R2 = 0.35, p < 0.001; Fig. 3c+d), IGFBP3 (CV R2 = 0.11, p < 0.001; Fig. 3e+f), but not IGF-1. For both insulin and insulin C-peptide, the most important metabolites for predictions were increased S-glucose, S-lactate and decreased S-Leucine. For IGFBP3, the most important metabolites were increased S-Acetone, S-Glycoprotein, and S-Leucine. We also found positive correlations between S-lactate and the preoperative increase in S-insulin and S-insulin / c-peptide (r = 0.57; p < 0.001 and r = 0.61; p < 0.0001), and between S-pyruvate and the increase in preoperative S-insulin and S-insulin c-peptide (r = 0.54; p < 0.001 and r = 0.60; p < 0.001).

Table 3 Serum metabolite values correlated to insulin (Pearson’s correlation) for the total study population, and the carbohydrate and fasting groups separately
Full size table
Fig. 3
figure 3

Correlation between serum metabolic profile and serum insulin, insulin C-peptide, and IGFBP3. Samples from carbohydrate-fed patients are shown in red, while samples from fasting patients are shown in blue. Metabolites are colored according to their variable importance in projection (VIP) score and labeled when VIP≥1. a Measured insulin vs. predicted insulin levels based on metabolic profile (cross-validated measurements). b Metabolites versus regression coefficient for insulin. Increased S-glucose, S-lactate, and decreased S-Leucine are important to prediction of serum insulin from the metabolic profile. c Measured insulin C peptide vs. predicted insulin C-peptide levels. d Regression weight plot showing metabolites versus the regression coefficient for insulin C-peptide. Increased S-Glucose, S-Lactate, and decreased S-Leucine are important to prediction of serum insulin C-peptide from the metabolic profile. e Measured Insulin Growth Factor Binding Protein 3 (IGFBP3) vs. predicted IGFBP3 based on metabolic profile. f Regression weight plot showing metabolites versus the regression coefficient for IGFBP3. Increased S-Acetone, S-Glycoproteins, and S-Leucine are important to prediction of serum IGFBP3 from the metabolic profile

Full size image

Tumor metabolism

Metabolites included in the analysis are presented in Table 4. PLS-DA did not result in a significant model discriminating between fasting and carbohydrate-fed patients, and no metabolites were significantly different in univariate testing when all tumors were analyzed (Fig. 4a). However, for ER-positive tumors (n = 18), glutathione was significantly elevated in the carbohydrate group compared to the fasting group (p = 0.002; Fig. 4b), even after adjusting for tumor size. In the ROC analysis, we found an area under the curve (AUC) of 0.894 (95%CI = 0.687–1.000, p = 0.0015) for glutathione in discriminating between fasting and carbohydrate-fed patients with ER-positive tumors (Fig. 4c). The difference was also significant in the ER-positive tumors with low proliferation (MAI < 10; n = 7). Moreover, we found a positive correlation between preoperative S-insulin levels and the glutathione content in tumor tissue (r = 0.680; p = 0.002). Furthermore, we observed a higher level of tissue glutamate in tumors with a high proliferation as measured by Ki67</≥ 15% (p = 0.004). This association remained significant when adjusted for intervention group using a general linear model with intervention status as fixed factor, Ki67</≥ 15% as random factor, and tissue Glutamate as dependent variable (p = 0.009). Also, choline (p = 0.002) and phosphoetanolamine (p = 0.019) were increased in T2 tumors compared to T1 tumors.

Table 4 Tumor metabolites with fold changes and t-test p-values in carbohydrate vs fasting groups
Full size table
Fig. 4
figure 4

a Principal Component Analysis (PCA) of tumor metabolites. No grouping of fasting vs carbohydrate groups observed. b Glutathione levels in ER positive tumors. c ROC curve for classification into carbohydrate or fasting group by glutathione concentration in ER-positive tumors. AUC= 0.894; 95%CI=0.0.687-1.000, P=0.002

Full size image

Survival analysis

First, we used S-lactate, S-pyruvate, and tissue (T) glutathione as continuous variables in a univariate Cox model for RFS, BCSS and OS. Both S-pyruvate and S-lactate, but not T-glutathione reached significance with a hazard ratio (HR) for RFS of 1.53 (95% CI, 1.11 to 2.11; p = 0.009) and 1.08 (95% CI, 1.01 to 1.17; p = 0.029), respectively. For BCSS the HR for the continuous variables of S-pyruvate and S-Lactate were 1.85 (95%CI, 1.15 to 2.97; p = 0.011) and 1.13 (95%CI, 1.01 1.26; p = 0.028) respectively. The corresponding observations for OS were 1.63 (95%CI, 1.11 to 2.40; p = 0.014) for lactate and 1.10 (95%CI, 1.002 to 1.20; p = 0.045) for pyruvate. Thereafter, the following independent variables were dichotomized according to the optimal ROC-derived thresholds: S-lactate, S-pyruvate, preoperative S-insulin, preoperative S-insulin-c-peptide, and tissue glutathione. In addition, the well-established prognostic factors tumor size, nodal status, histological grade, MAI 10, Ki-67-30 and PPH3–13 were deemed clinically relevant and included as explanatory variables in the multivariable analyses. The results of the univariate RFS, BCSS, and OS analyses are given in Tables 5, 6 and 7, respectively. Patients with a high glutathione content in the tumor (≥1.09) had a 37% risk of experiencing a relapse and 37% risk of dying of breast cancer compared to no relapses and no deaths in patients with a low glutathione content in the tumor (both comparisons: p = 0.038; HR = Inf.; Fig. 5a and d). Patients with high S-lactate (≥56.9) had RFS of 71% compared to 97% for those with lower S-lactate (p = 0.002, HR = 7.47; 95% CI 1.66–33.6; Fig. 5b). Patients with S-pyruvate ≥12.5 had an adverse RFS of 50% compared to 95% for the patients with S-pyruvate < 12.5 (p < 0.0001; HR = 13.6; 95% CI 2.61–70.6; Fig. 5c). The same pattern was observed in the BCSS and OS analyses for these three prognostic variables (Fig. 5e-i). Notably, only one contralateral relapse occurred in the fasting group – all others were in the carbohydrate group. Even though the relapses were restricted to patients with T2 tumors, tumor category was not an independent prognostic factor in the multivariable analyses. In the multivariable analysis for RFS, S-pyruvate was the only factor left in the final model (HR = 12.8; 95% CI, 2.47 to 66.8), and only S-lactate remained in the final multivariable model for BCSS (HR = 14.8; 95% CI 1.54 to 142). Furthermore, S-pyruvate was the sole factor to reach significance in the multivariable model of the OS analysis (HR = 18.2; 95% CI 2.03 to 164).

Table 5 Univariate analysis of Relapse Free Survival in ER+ patients
Full size table
Table 6 Univariate analysis of Breast Cancer Specific Survival in ER+ patients
Full size table
Table 7 Univariate analysis of Overall Survival in ER+ patients
Full size table
Fig. 5
figure 5

Survival analyses for Tumor-Glutathione, Serum-lactate and Serum-pyruvate. a-c Relapse Free Survival (RFS); d-f Breast Cancer Survival (BCSS); g-i Overall Survival (OS)

Full size image

Pathway analyses

In the Pathway analyses, MetaboAnalyst and IPA showed complimentary information. Quantitative metabolite set enrichment analysis (MSEA) identified biologically meaningful patterns in serum metabolite concentration changes (Fig. 6a and Table 12 in Appendix). Significantly enriched pathways included energy associated metabolic pathways (amino sugar metabolism and pyruvate metabolism which links to glutamate metabolism, the citric acid cycle, gluconeogenesis and the Warburg effect). IPA showed the main functions of the involved metabolites as cellular growth and proliferation, molecular transport, small molecule biochemistry, carbohydrate metabolism and amino acid metabolism (Fig. 6b). Interestingly, the metabolites showed a pattern congruent with growth of organism (Fig. 6c) with metabolites increased in carbohydrate-fed patients activating growth pathways, and downregulation of metabolites acting as inhibitors of growth. Finally, four (miR-26a-5p, miR-30c-5p, miR126-3p and miR-210-3p) out of the seven microRNAs found to be involved in resistance to tamoxifen in our previous review [27] could indirectly be associated with the metabolic network through insulin signaling pathways (Fig. 6d). The same metabolic pathways were evident when only ER positive patients were considered.

Fig. 6
figure 6

Pathway analyses in serum metabolites. a Metabolite Set Enrichment Analysis of serum metabolism. Significantly enriched pathways are annotated in the pathway network. The circle size denotes significance of the pathway, and lines denote at least 25% shared metabolites in the pathways. b Ingenuity pathway analysis (IPA) bar chart showing the top 5 functions enriched in the dataset. c IPA pathway network showing the metabolites connected to four microRNAs found to be involved in tamoxifen resistance. Metabolites in green are downregulated in carbohydrate-fed patients, while metabolites in red are upregulated. MicroRNAs are colored purple. d IPA Function plot showing metabolites involved in organismal growth. Orange arrows indicate activation, while blue arrows indicate inhibition

Full size image

Discussion

We present the first study to examine the effect of per-oral preoperative carbohydrate load on perioperative metabolism in operable breast cancer patients. Among the 15 different serum metabolites that distinguished fasting from the per-oral carbohydrate load, we observed increased systemic lactate and pyruvate, decreased ketone bodies, increased glycerol, and reduced amino acids in the patients who received the carbohydrate load. Moreover, we found highly significant positive correlations between S-insulin and S-lactate and S-pyruvate. Thus, changes in these 15 key metabolites are consistent with increased glycolysis, increased ketolytic activity, reduced lipolysis, and reduced proteolysis, which are exactly the same metabolic modifications seen after carbohydrate challenge in healthy persons [28]. Being able to capture these well-known metabolic effects of insulin increases the reliability of our model to detect other changes that may follow a carbohydrate load.

It may be considered that 18 h is too short to expect effect of the carbohydrate load on tumor cell proliferation and metabolism. However, in vitro studies show that glucose fed MCF-7 cells increase their proliferation after 12–24 h [3] .Others found the same pattern in three different breast cancer cell lines [29]. As the cell lines lack the in vivo endocrine response to glucose the increased proliferation was based on GTP-ase driven phosphorylation of EGFR with increased activity and longevity of this receptor as a consequence. Also, animals fed with a diet containing increased glucose show an increased epithelial mesenchymal transition (EMT) [30].

The increased S-lactate and S-pyruvate in the carbohydrate patients stems primarily from two sources. Firstly, lactate is the product of glycolysis, especially in muscle cells, and is transported to the liver for conversion back to glucose, known as the Cori cycle [31]. The intended effect of preOP is to contra act and reduce insulin resistance that follows surgical stress [32]. This stressor leads to reduced mitochondrial ATP production and lactate formation [33, 34]. In healthy individuals, an oral glucose tolerance test (OGTT) showed a negative correlation between differences in S-glucose concentrations and differences in S-lactate levels (i.e. a rise in S-glucose leads to a reduction in S-lactate) [35]. Moreover, during 180 min after an OGTT among non-insulin dependent diabetic mellitus (NIDDM) patients there was no significant alteration in S-lactate levels [36]. Thus, it is unlikely that preOp itself creates a systemic lactate production. Therefore, S-lactate in our patients may come from excretion of intracellular lactate and pyruvate produced in the breast cancer cells. Consequently, lactate and pyruvate in the present study are probably translocated into the systemic circulation via mono carboxylate transporter type 4 (MCT-4), which is a known part of the Warburg effect [13]. Despite the fact that systemic metabolite concentrations are functional read outs of the numerous homeostatic reactions in the body, which will blur the contribution from the cancer cell metabolism to the serum levels [18], our present observation of positive correlation between larger tumor size and increasing S-lactate is supported by Hui S et al. [37].. Also, the positive correlation between proliferation and tumor size solely occurs in the carbohydrate group this suggests that carbohydrate exposure to larger tumors (i.e.T2 tumors) increases both proliferation and S-lactate. Thus, this indicates that lactate from the Warburg effect in the tumor cells may have a substantial contribution to the systemic lactate and pyruvate levels. This observation also adheres to the lack of correlation between intra tumor lactate/pyruvate and fasting/carbohydrate status in the present study, as the former are probably excreted from the cells into the systemic environment.

Moreover, tumor cells not only produce lactate for excretion through MCT-4, they are also able to take up systemic circulating lactate and pyruvate via the MCT-1 transporters [38]. Regardless of the source, systemic lactate and pyruvate will certainly benefit the free CTCs shed from the tumor during surgery that are on their way to distant tissue to form micrometastases [39] but may also benefit the preoperatively established occult micrometastases [40, 41]. Lactate and pyruvate are the most preferred substrates for lactate/pyruvate dehydrogenase (LDH/PDH), ensuring a 1:1 ratio between lactate and pyruvate when equilibrium is reached. Thus, LDH provides substrate for both the production of ATP via the tricarboxylic acid (TCA) cycle [37] and also increased gluconeogenesis for the production of ribose for nucleotide synthesis via the pentose phosphate pathway (PPP) [38]. Notably, increased levels of serum LDH [42] and increased expression of LDH in breast cancer tissue [43] and lung tumors [44] are associated with an inferior prognosis.

In line with other studies [45], we observed a positive correlation between higher proliferation and increased glutamate content in tumor tissue. Glutamate is a metabolic product of glutaminolysis, which drives membrane trafficking to promote breast cancer cell invasiveness [46]. In addition, the expression of glutaminase genes GLS and GLS2 correlates with increased tumor growth rates [47]. Many tumors become glutamine-dependent, as it serves as a direct route into the TCA cycle at the alpha-ketoglutaric acid level with consequential ATP production. Together with glycine and cysteine, glutamate is a precursor to the tripeptide glutathione, which is an antioxidant molecule that serves to ‘buffer’ superoxide insults encountered in the tumor microenvironment [45]. Glutathione is the major thiol-containing endogenous antioxidant and serves as a redox buffer against various sources of oxidative stress. In tumors, maintaining a supply of glutathione is critical for cellular survival because it allows cells to resist the oxidative stress associated with rapid metabolism, DNA-damaging agents, and inflammation, among others [48, 49]. Glucose metabolism and biosynthesis of glutathione are often modulated by the PI3K/Akt pathway, which is often dysregulated in breast cancer tumors [50, 51]. Importantly, one of the effects of targeting the PI3K/Akt-pathway upstream [52] and downstream [19] is reduced glutathione content in tumor cells. In the PPP-pathway, NAD+ and NADP are converted into NADH and NADPH, respectively, which contribute to maintaining glutathione (GSSG) in the reduced state (GSH) [53]. Thus, the PPP-pathway in the Warburg effect secures a high intracellular level of glutathione, which is regarded as the most important cellular protection system against attack from reactive oxygen species (ROS) in both dividing and hibernating luminal cells [10], and also in cancer stem cells [54]. Thus, preoperative carbohydrate loading seems to create a doubly favorable environment that will probably serve the CTCs liberated during surgery [39] more than the already established micrometastases [55]. First, CTCs have a surplus of cellular fuel via lactate and pyruvate available systemically. Second, they benefit from an increased level of intracellular protection systems against ROS via increased tumor glutathione. Both effects will increase the probability of CTCs thriving and surviving as micrometastases, which then may erupt as clinical relapse years later, compatible with the tumor biology of luminal breast cancers. However, our observed clinical endpoint between 3 to 7 years must be regarded as ‘early relapses’ when coming to luminal cancers [56]. Thus, we need a much longer follow up to capture the late recurrences in order to get the correct picture of the clinical outcome of the present study.

Several attempts have been made to reverse the above-mentioned metabolic pathways for treatment purposes. The first attempt was to reverse the Warburg effect with the polyphenol resveratrol, which blocks PDH/LDH. In colon cancer cells, resveratrol inhibits proliferation, gluconeogenesis, and PPP [57]. By blocking PDH, resveratrol promotes mitochondrial electron transport chain overload with increased ROS production, ultimately resulting in apoptosis [58]. Secondly, a ketogenic diet has been shown to be effective in preclinical studies [59]. A ketogenic diet produces a large amount of intracellular ketone bodies that have a direct cytotoxic effect. Furthermore, the ketogenic state inhibits insulin/IGF signaling and downstream signaling pathways, such as PI3K/Akt/mTOR [60]. Interestingly, in the present study, the patients in the fasting group reached a ketogenic state with increased ketone bodies, which may have created an unfavorable environment for the cancer cells in the tumor and for the liberated CTCs. This is in line with a recent RCT of using ketogenic diet as adjuvant treatment in one of the study arms. They observed a better overall survival in the group that received ketogenic diet [61]. Others have recently shown a profound effect of ketogenic diet in a xenografted breast cancer mouse model with increased ketone bodies and increased aminoacidic [62], which is in line with our observations. The authors hypothesize that the anti-cancer effect may be mediated through immunological mechanisms [62]. Thus, use of a ketogenic diet as adjuvants to conventional therapy is rooted in several studies [63].

Likewise, physical activity is known to prevent and improve survival in several cancer forms and is thus recommended as a measure to both prevent and treat breast cancer [64, 65]. One of the mechanisms behind these observation is a change in the estrogen metabolism after 180 min exercise pr. week. They found an increased 2 hydroxy-estrone level known to antagonize the estradiol action [66] This observation is important for both in the preventive setting as breast cancer risk is correlated to total life exposure of estrogens [67]. Also, changes in diet affect the cancer incidence [68], and also prognosis in breast cancer patients [69].

A combination of calorie restriction and physical exercise in postmenopausal women did also reduce insulin levels [70]. In our patients, we found that metabolic changes after the carbohydrate load affected the ER-positive breast cancer patients. Thus, ketogenic diet combined with physical exercise would probably be beneficial for our patients as this approach will affect both the ER and insulin signaling pathways.

Interestingly, intermittent fasting (i.e. caloric restriction for 16–48 h [71] has been proven to affect the metabolism and disease process in a beneficial manner. Notably, intermittent fasting in animal studies have demonstrated reduction of tumor size [72]. In humans, intermittent fasting improves insulin sensitivity and thus reduces insulin and IGF-1 related signaling in over weighted individuals [72, 73]. Preclinical studies show that intermittent fasting more than 2 days is as effective as chemotherapy to reduce cancer load [74]. Thus, the ketones derived from intermittent fasting decreases cancer cell viability by attacking several hallmarks of cancer [75].

The IPA-analyses confirmed that the systemic response to the carbohydrate load converge towards pathways involved in proliferation and growth of the organism. Moreover, other pathways related to the Warburg effect were also involved. Thus, peroral preoperative carbohydrate load shifts the systemic metabolism towards a very fortunate and beneficial environment for CTCs liberated from the tumor under the operation. Interestingly, four out of seven microRNAs related to endocrine resistance [27] also regulate the same metabolic pathways through insulin signaling pathways, which are known to be involved in endocrine resistance with reduced effect of tamoxifen and aromatase inhibitors. Thus, it seems plausible to introduce metformin early on as adjuvant treatment to regain the endocrine sensitivity. Intriguingly, circulating microRNAs from the tumor in exosomes [76] can perform cell-independent microRNA biogenesis and promote tumorigenesis away from the primary tumor [77]. Thus, we may speculate that one of the steps in the metastatic process is to control the systemic metabolic pathways to ensure a beneficial environment and survival of the liberated cancer cells [54]. Moreover, increased cellular uptake of glucose via the Warburg effect [10] favor differentiating glycosylation of intracellular proteins included paucimannosylation [78]. Intriguingly, the metastatic Epithelial-Mesenchymal-Transition (EMT) process is regulated through glycosylation of key regulator proteins, that are frequently modulated via the insulin /IGF signaling [79]. Thus, glycosylation opens up a connection between the glucose/insulin signaling and increased survival of CTCs trough enhancement of the EMT-processes.

Taken together, this explorative study indicates that the carbohydrate loading state and fasting state have opposite systemic and micro-environmental effects, which may explain why the relapses in the present study were skewed towards the carbohydrate group, with an inferior RFS, BCSS, and OS in patients with high tissue glutathione, high S-lactate, and high S-pyruvate. The favorable macro- and micro-environmental changes for the tumor that come from carbohydrate loading reflect the Warburg effect, which serves the CTCs and micrometastases more than the patient [80]. In luminal cancers, the Warburg pathway enzyme PFKFB4 acts as a molecular fulcrum that couples sugar metabolism to transcriptional activation by stimulating the ER co-activator SRC-3 to promote aggressive metastatic tumors [81].

The present study has several weak points. First, it is a post hoc explorative analysis of an RCT. Therefore, the various analyses are not sufficiently powered regarding the various endpoints. In addition, tissue samples were not available for all patients, which reduces the number of patients in the various analyses. Thus, this creates a greater risk of a type II error than a type I error. Furthermore, the tissue analyses were skewed towards patients with larger tumors. This could introduce systematic error in the analysis. However, tumor size was not included in the final Cox models in any survival analysis, indicating that this error was not strong enough to blur the effects of the metabolites. Also, including diet recalls and demographic data of the patients would have strengthened the study. Detecting the well-known endocrine metabolic fingerprint of insulin strengthens the method and the reliability of the various findings in this study. However, the study is too small to conclude on preoperative preparation guidelines; fasting or carbohydrate loading. Moreover, the pilot nature of the present study calls for validation in a larger study with a long-term follow-up. Introducing a ketogenic diet as a third study arm may test out whether ketone bodies could wipe out the liberated CTCs and thus improve survival.

Conclusion

Preoperative oral glucose loading increases systemic levels of lactate and pyruvate, and tumor levels of glutathione and glutamate in luminal breast cancer patients. In fasting patients, the proapoptotic ketone bodies are increased. These biological changes may contribute to the survival differences observed between these two study groups. Integrated Pathway Analysis (IPA) in serum revealed activation of five major anabolic metabolic networks contributing to proliferation and growth mainly through insulin signaling pathways.

Availability of data and materials

The data that support the findings of this study are available from Stavanger Breast Cancer Research Group, but restrictions apply to the availability of these data, which were used under license for the current study and as such are not publicly available. However, data are available from the authors upon reasonable request and with permission from Stavanger Breast Cancer Research Group.

Abbreviations

ATP:

Adenosine triphosphate

AUC:

Area under the curve

BCSS:

Breast cancer specific survival

CI:

Confidence interval

CPMG:

Carr–Purcell–Meiboom–Gill spectra

CTC:

Circulating tumor cell

EMT:

Epithelial mesenchymal transition

ER:

Estrogen receptor

ERAS:

Enhanced recovery after surgery

GSH:

Glutathione, reduced form

GSSG:

Glutathione, oxidized form

HER2:

Human epithelial growth factor receptor 2

HR:

Hazard ratio

HR-MAS-MR:

High-resolution magic angle spinning - magnetic resonance

IGF1:

Insulin-like growth factor 1

IGF1R:

Insulin-like growth factor 1 receptor

IGFBP3:

Insulin-like growth factor 1 binding protein 3

IHC:

Immunohistochemistry

IPA:

Ingenuity Pathway Analysis

IR:

Insulin receptor

LDH:

Lactate dehydrogenase

MAI:

Mitotic activity index

MCT-1:

Monocarboxylate transporter 1

MCT-4:

Monocarboxylate transporter 4

MR:

Magnetic resonance

MRI:

Magnetic resonance imaging

MSEA:

Metabolite Set Enrichment Analysis

NAD:

Nicotinamide adenine dinucleotide

NADP:

Nicotinamide adenine dinucleotide phosphate

NOESY:

One-dimensional 1H Nuclear Overhauser effect spectroscopy

NSD:

Norwegian Center for Research Data

OS:

Overall survival

PHD:

Pyruvate dehydrogenase

PLS:

Partial least square

PLS-DA:

Partial least square discriminant analysis

PPH3:

Phosphorylated phosphohistone 3

PPP:

Pentose phosphate pathway

PR:

Progesterone receptor

RCT:

Randomized controlled trial

RFS:

Relapse free survival

ROS:

Reactive oxygen species

S:

Serum

SRC-3:

Steroid receptor co-activator 3

TCA-cycle:

Tri carboxyl acid cycle

References

  1. Lundqvist A, Andersson E, Ahlberg I, Nilbert M, Gerdtham U. Socioeconomic inequalities in breast cancer incidence and mortality in Europe-a systematic review and meta-analysis. Eur J Pub Health. 2016;26(5):804–13.

    Article  Google Scholar 

  2. Yap YS, Lu YS, Tamura K, Lee JE, Ko EY, Park YH, Cao AY, Lin CH, Toi M, Wu J, Lee SC. Insights into breast Cancer in the east vs the west: a review. JAMA Oncol. 2019. https://doi.org/10.1001/jamaoncol.2019.0620.

    Article  PubMed  Google Scholar 

  3. Sun YS, Zhao Z, Yang ZN, Xu F, Lu HJ, Zhu ZY, Shi W, Jiang J, Yao PP, Zhu HP. Risk factors and preventions of breast Cancer. Int J Biol Sci. 2017;13(11):1387–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Baum M, Demicheli R, Hrushesky W, Retsky M. Does surgery unfavourably perturb the "natural history" of early breast cancer by accelerating the appearance of distant metastases? Eur J Cancer. 2005;41(4):508–15.

    Article  PubMed  Google Scholar 

  5. Pukazhendhi G, Gluck S. Circulating tumor cells in breast cancer. J Carcinog. 2014;13:8. https://doi.org/10.4103/1477-3163.135578.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Price TT, Burness ML, Sivan A, Warner MJ, Cheng R, Lee CH, Olivere L, Comatas K, Magnani J, Kim Lyerly H, et al. Dormant breast cancer micrometastases reside in specific bone marrow niches that regulate their transit to and from bone. Sci Transl Med. 2016;8(340):340–73.

    Article  Google Scholar 

  7. Wangchinda P, Ithimakin S. Factors that predict recurrence later than 5 years after initial treatment in operable breast cancer. World J Surg Oncol. 2016;14(1):223.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.

    Article  CAS  PubMed  Google Scholar 

  9. Zhao H, Shen J, Moore SC, Ye Y, Wu X, Esteva FJ, Tripathy D, Chow WH. Breast cancer risk in relation to plasma metabolites among Hispanic and African American women. Breast Cancer Res Treat. 2019. https://doi.org/10.1007/s10549-019-05165-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009;324(5930):1029–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Warburg OPK, Negelein E. Ueber den stoffwechhsel der tumoren. Biochem Z. 1924;152(1):319–44.

    Google Scholar 

  12. Fernandez-de-Cossio-Diaz J, Vazquez A. Limits of aerobic metabolism in cancer cells. Sci Rep. 2017;7(1):13488.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Tran Q, Lee H, Park J, Kim SH, Park J. Targeting Cancer metabolism - revisiting the Warburg effects. Toxicol Res. 2016;32(3):177–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Hart CD, Tenori L, Luchinat C, Di Leo A. Metabolomics in breast Cancer: current status and perspectives. Adv Exp Med Biol. 2016;882:217–34.

    Article  CAS  PubMed  Google Scholar 

  15. Ljungqvist O. ERAS--enhanced recovery after surgery: moving evidence-based perioperative care to practice. JPEN J Parenter Enteral Nutr. 2014;38(5):559–66.

    Article  PubMed  Google Scholar 

  16. Lende TH, Austdal M, Varhaugvik AE, Skaland I, Gudlaugsson E, Kvaløy JT, Akslen LA, Søiland H, Janssen EAM, Baak JPA. Influence of pre-operative oral carbohydrate loading vs. standard fasting procedure on tumor proliferation and clinical outcome in breast cancer patients — a randomized trial. BMC Cancer. 2019;19(1):1076.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Austdal M, Tangeras LH, Skrastad RB, Salvesen K, Austgulen R, Iversen AC, Bathen TF. First trimester urine and serum metabolomics for prediction of preeclampsia and gestational hypertension: a prospective screening study. Int J Mol Sci. 2015;16(9):21520–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Giskeodegard GF, Madssen TS, Euceda LR, Tessem MB, Moestue SA, Bathen TF. NMR-based metabolomics of biofluids in cancer. NMR Biomed. 2018:e3927. https://doi.org/10.1002/nbm.3927.

  19. Euceda LR, Hill DK, Stokke E, Hatem R, El Botty R, Bièche I, Marangoni E, Bathen TF, Moestue SA. Metabolic response to Everolimus in patient-derived triple-negative breast Cancer Xenografts. J Proteome Res. 2017;16(5):1868–79.

    Article  CAS  PubMed  Google Scholar 

  20. Eilers PH. Parametric time warping. Anal Chem. 2004;76(2):404–11.

    Article  CAS  PubMed  Google Scholar 

  21. Sitter B, Sonnewald U, Spraul M, Fjosne HE, Gribbestad IS. High-resolution magic angle spinning MRS of breast cancer tissue. NMR Biomed. 2002;15(5):327–37.

    Article  CAS  PubMed  Google Scholar 

  22. Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001;29(4):1165–88.

    Article  Google Scholar 

  23. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/.

  24. pls: Partial Least Squares and Principal Component Regression. R package version 2.7–1. https://cran.r-project.org/web/packages/pls/index.html

  25. Xia JMR, Sinelnikov IV, Broadhurst D, Wishart DS. MetaboAnalyst 2.0—a comprehensive server for metabolomic data analysis. Nucleic Acids Res. 2012;40(web server issue). https://doi.org/10.1093/nar/gks374.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Mehmood TLK, Snipen L, Sæbø S. A review of variable selection methods in partial least squares regression. Chemom Intell Lab Syst. 2012;118:62–9.

    Article  CAS  Google Scholar 

  27. Egeland NG, Lunde S, Jonsdottir K, Lende TH, Cronin-Fenton D, Gilje B, Janssen EA, Søiland H. The role of MicroRNAs as predictors of response to Tamoxifen treatment in breast Cancer patients. Int J Mol Sci. 2015;16(10):24243–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Shaham O, Wei R, Wang TJ, Ricciardi C, Lewis GD, Vasan RS, Carr SA, Thadhani R, Gerszten RE, Mootha VK. Metabolic profiling of the human response to a glucose challenge reveals distinct axes of insulin sensitivity. Mol Syst Biol. 2008;4:214.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Hou Y, Zhou M, Xie J, Chao P, Feng Q, Wu J. High glucose levels promote the proliferation of breast cancer cells through GTPases. Breast Cancer (Dove Med Press). 2017;9:429–36.

    CAS  Google Scholar 

  30. Tang FY, Pai MH, Chiang EP. Consumption of high-fat diet induces tumor progression and epithelial-mesenchymal transition of colorectal cancer in a mouse xenograft model. J Nutr Biochem. 2012;23(10):1302–13.

    Article  CAS  PubMed  Google Scholar 

  31. Passarella S. Schurr a: l-lactate transport and metabolism in mitochondria of Hep G2 cells-the Cori cycle revisited. Front Oncol. 2018;8:120.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Scott MJ, Fawcett WJ. Oral carbohydrate preload drink for major surgery - the first steps from famine to feast. Anaesthesia. 2014;69(12):1308–13.

    Article  CAS  PubMed  Google Scholar 

  33. Atkins RPVK, Constantin-Teodosiu D, Lobo DN, Greenhaff PL. Rates of skeletal muscle mitochondrial ATP production are reduced during elective abdominal surgery in humans. J Am Coll Surg. 2011;2013:S59.

    Article  Google Scholar 

  34. Varadan KKAR, Dimitru CT, Blackshaw E, Perkins AC, Greenhaff PL, Lobo DN. Gastrointestinal surgery mediated increases in gut permeability and expression of IL6 and PDK4 mRNAs in quadriceps muscle may underpin the post-operative increase in whole-body insulin resistance in humans. J Am Coll Surg. 2011;2013:S53.

    Article  Google Scholar 

  35. Prando R, Cheli V, Buzzo P, Melga P, Ansaldi E, Accoto S. Blood lactate behavior after glucose load in diabetes mellitus. Acta Diabetol Lat. 1988;25(3):247–56.

    Article  CAS  PubMed  Google Scholar 

  36. Moore MC, Davis SN, Mann SL, Cherrington AD. Acute fructose administration improves oral glucose tolerance in adults with type 2 diabetes. Diabetes Care. 2001;24(11):1882–7.

    Article  CAS  PubMed  Google Scholar 

  37. Hui S, Ghergurovich JM, Morscher RJ, Jang C, Teng X, Lu W, Esparza LA, Reya T, Le Z, Yanxiang Guo J, et al. Glucose feeds the TCA cycle via circulating lactate. Nature. 2017;551(7678):115–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Faubert B, Li KY, Cai L, Hensley CT, Kim J, Zacharias LG, Yang C, Do QN, Doucette S, Burguete D, et al. Lactate metabolism in human lung tumors. Cell. 2017;171(2):358–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Papavasiliou P, Fisher T, Kuhn J, Nemunaitis J, Lamont J. Circulating tumor cells in patients undergoing surgery for hepatic metastases from colorectal cancer. Proc (Baylor Univ Med Cent). 2010;23(1):11–4.

    Article  Google Scholar 

  40. Braun S, Vogl FD, Naume B, Janni W, Osborne MP, Coombes RC, Schlimok G, Diel IJ, Gerber B, Gebauer G, et al. A pooled analysis of bone marrow micrometastasis in breast cancer. N Engl J Med. 2005;353(8):793–802.

    Article  CAS  PubMed  Google Scholar 

  41. de Boer M, van Dijck JA, Bult P, Borm GF, Tjan-Heijnen VC. Breast cancer prognosis and occult lymph node metastases, isolated tumor cells, and micrometastases. J Natl Cancer Inst. 2010;102(6):410–25.

    Article  PubMed  Google Scholar 

  42. Bidard FC, Hajage D, Bachelot T, Delaloge S, Brain E, Campone M, Cottu P, Beuzeboc P, Rolland E, Mathiot C, Pierga JY. Assessment of circulating tumor cells and serum markers for progression-free survival prediction in metastatic breast cancer: a prospective observational study. Breast Cancer Res. 2012;14(1):R29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Liu D, Wang D, Wu C, Zhang L, Mei Q, Hu G, Long G, Sun W. Prognostic significance of serum lactate dehydrogenase in patients with breast cancer: a meta-analysis. Cancer Manag Res. 2019;11:3611–9.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Liu L, He Y, Ge G, Li L, Zhou P, Zhu Y, Tang H, Huang Y, Li W, Zhang L. Lactate dehydrogenase and creatine kinase as poor prognostic factors in lung cancer: a retrospective observational study. PLoS One. 2017;12(8):e0182168.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. Choi JS, Yoon D, Koo JS, Kim S, Park VY, Kim EK, Kim S, Kim MJ. Magnetic resonance metabolic profiling of estrogen receptor-positive breast cancer: correlation with currently used molecular markers. Oncotarget. 2017;8(38):63405–16.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Dornier E, Rabas N, Mitchell L, Novo D, Dhayade S, Marco S, Mackay G, Sumpton D, Pallares M, Nixon C, et al. Glutaminolysis drives membrane trafficking to promote invasiveness of breast cancer cells. Nat Commun. 2017;8(1):2255.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Knox WE, Horowitz ML, Friedell GH. The proportionality of glutaminase content to growth rate and morphology of rat neoplasms. Cancer Res. 1969;29(3):669–80.

    CAS  PubMed  Google Scholar 

  48. Estrela JM, Ortega A, Obrador E. Glutathione in cancer biology and therapy. Crit Rev Clin Lab Sci. 2006;43(2):143–81.

    Article  CAS  PubMed  Google Scholar 

  49. Franco R, Cidlowski JA. Apoptosis and glutathione: beyond an antioxidant. Cell Death Differ. 2009;16(10):1303–14.

    Article  CAS  PubMed  Google Scholar 

  50. Lien EC, Lyssiotis CA, Juvekar A, Hu H, Asara JM, Cantley LC, Toker A. Glutathione biosynthesis is a metabolic vulnerability in PI (3) K/Akt-driven breast cancer. Nat Cell Biol. 2016;18(5):572–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. DeBerardinis RJ, Cheng T. Q's next: the diverse functions of glutamine in metabolism, cell biology and cancer. Oncogene. 2010;29(3):313–24.

    Article  CAS  PubMed  Google Scholar 

  52. Fack F, Espedal H, Keunen O, Golebiewska A, Obad N, Harter PN, Mittelbronn M, Bahr O, Weyerbrock A, Stuhr L, et al. Bevacizumab treatment induces metabolic adaptation toward anaerobic metabolism in glioblastomas. Acta Neuropathol. 2015;129(1):115–31.

    Article  CAS  PubMed  Google Scholar 

  53. Liberti MV, Locasale JW. The Warburg effect: how does it benefit Cancer cells? Trends Biochem Sci. 2016;41(3):211–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Lee SY, Ju MK, Jeon HM, Lee YJ, Kim CH, Park HG, Han SI, Kang HS. Oncogenic metabolism acts as a prerequisite step for induction of Cancer metastasis and Cancer stem cell phenotype. Oxidative Med Cell Longev. 2018;2018:1027453.

    Google Scholar 

  55. Klauber-DeMore N, Van Zee KJ, Linkov I, Borgen PI, Gerald WL. Biological behavior of human breast cancer micrometastases. Clin Cancer Res. 2001;7(8):2434–9.

    CAS  PubMed  Google Scholar 

  56. Engstrom MJ, Opdahl S, Hagen AI, Romundstad PR, Akslen LA, Haugen OA, Vatten LJ, Bofin AM. Molecular subtypes, histopathological grade and survival in a historic cohort of breast cancer patients. Breast Cancer Res Treat. 2013;140(3):463–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Saunier E, Antonio S, Regazzetti A, Auzeil N, Laprevote O, Shay JW, Coumoul X, Barouki R, Benelli C, Huc L, Bortoli S. Resveratrol reverses the Warburg effect by targeting the pyruvate dehydrogenase complex in colon cancer cells. Sci Rep. 2017;7(1):6945.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. Blanquer-Rossello MD, Hernandez-Lopez R, Roca P, Oliver J, Valle A. Resveratrol induces mitochondrial respiration and apoptosis in SW620 colon cancer cells. Biochim Biophys Acta Gen Subj. 2017;1861(2):431–40.

    Article  CAS  PubMed  Google Scholar 

  59. Weber DD, Aminazdeh-Gohari S, Kofler B. Ketogenic diet in cancer therapy. Aging (Albany NY). 2018;10(2):164–5.

    Article  Google Scholar 

  60. Tan-Shalaby J. Ketogenic diets and Cancer: emerging evidence. Fed Pract. 2017;34(Suppl 1):37S–42S.

    PubMed  PubMed Central  Google Scholar 

  61. Khodabakhshi A, Akbari ME, Mirzaei HR, Mehrad-Majd H, Kalamian M, Davoodi SH. Feasibility, safety, and beneficial effects of MCT-based Ketogenic diet for breast Cancer treatment: a randomized controlled trial study. Nutr Cancer. 2019:1–8. [Epub ahead of print].

  62. Licha D, Vidali S, Aminzadeh-Gohari S, Alka O, Breitkreuz L, Kohlbacher O, Reischl RJ, Feichtinger RG, Kofler B, Huber CG. Untargeted metabolomics reveals molecular effects of Ketogenic diet on healthy and tumor Xenograft mouse models. Int J Mol Sci. 2019;20(16):E3873.

    Article  PubMed  Google Scholar 

  63. Paoli A, Rubini A, Volek JS, Grimaldi KA. Beyond weight loss: a review of the therapeutic uses of very-low-carbohydrate (ketogenic) diets. Eur J Clin Nutr. 2013;67(8):789–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. McTiernan A, Friedenreich CM, Katzmarzyk PT, Powell KE, Macko R, Buchner D, Pescatello LS, Bloodgood B, Tennant B, Vaux-Bjerke A, et al. Physical activity in Cancer prevention and survival: a systematic review. Med Sci Sports Exerc. 2019;51(6):1252–61.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Irwin ML, Smith AW, McTiernan A, Ballard-Barbash R, Cronin K, Gilliland FD, Baumgartner RN, Baumgartner KB, Bernstein L. Influence of pre- and postdiagnosis physical activity on mortality in breast cancer survivors: the health, eating, activity, and lifestyle study. J Clin Oncol. 2008;26(24):3958–64.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Smith AJ, Phipps WR, Thomas W, Schmitz KH, Kurzer MS. The effects of aerobic exercise on estrogen metabolism in healthy premenopausal women. Cancer Epidemiol Biomark Prev. 2013;22(5):756–64.

    Article  CAS  Google Scholar 

  67. Yager JD, Davidson NE. Estrogen carcinogenesis in breast Cancer. N Engl J Med. 2006;354(3):270–82.

    Article  CAS  PubMed  Google Scholar 

  68. Donaldson MS. Nutrition and cancer: a review of the evidence for an anti-cancer diet. Nutr J. 2004;3:19.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  69. Thomas F, Rome S, Mery F, Dawson E, Montagne J, Biro PA, Beckmann C, Renaud F, Poulin R, Raymond M, Ujvari B. Changes in diet associated with cancer: an evolutionary perspective. Evol Appl. 2017;10(7):651–7.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Chlebowski RT, Pettinger M, Stefanick ML, Howard BV, Mossavar-Rahmani Y, McTiernan A. Insulin, physical activity, and caloric intake in postmenopausal women: breast cancer implications. J Clin Oncol. 2004;22(22):4507–13.

    Article  PubMed  Google Scholar 

  71. Mattson MP, Longo VD, Harvie M. Impact of intermittent fasting on health and disease processes. Ageing Res Rev. 2017;39:46–58.

    Article  PubMed  Google Scholar 

  72. Harvie MN, Howell T. Could intermittent energy restriction and intermittent fasting reduce rates of Cancer in obese, overweight, and Normal-weight subjects? A summary of evidence. Adv Nutr. 2016;7(4):690–705.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Longo VD, Panda S. Fasting, circadian rhythms, and time-restricted feeding in healthy lifespan. Cell Metab. 2016;23(6):1048–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Lee C, Raffaghello L, Brandhorst S, Safdie FM, Bianchi G, Martin-Montalvo A, Pistoia V, Wei M, Hwang S, Merlino A, et al. Fasting cycles retard growth of tumors and sensitize a range of cancer cell types to chemotherapy. Sci Transl Med. 2012;4(124):ra127.

    Article  Google Scholar 

  75. Poff AM, Ari C, Arnold P, Seyfried TN, D'Agostino DP. Ketone supplementation decreases tumor cell viability and prolongs survival of mice with metastatic cancer. Int J Cancer. 2014;135(7):1711–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Xu R, Rai A, Chen M, Suwakulsiri W, Greening DW, Simpson RJ. Extracellular vesicles in cancer - implications for future improvements in cancer care. Nat Rev Clin Oncol. 2018;15(10):617–38.

    Article  CAS  PubMed  Google Scholar 

  77. Melo SA, Sugimoto H, O'Connell JT, Kato N, Villanueva A, Vidal A, Qiu L, Vitkin E, Perelman LT, Melo CA, et al. Cancer exosomes perform cell-independent microRNA biogenesis and promote tumorigenesis. Cancer Cell. 2014;26(5):707–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Chatterjee S, Lee LY, Kawahara R, Abrahams JL, Adamczyk B, Anugraham M, Ashwood C, Sumer-Bayraktar Z, Briggs MT, Chik JHL, et al. Protein Paucimannosylation is an Enriched N-glycosylation Signature of Human Cancers. Proteomics. 2019:e1900010. https://doi.org/10.1002/pmic.201900010.

    Article  CAS  Google Scholar 

  79. de-Freitas-Junior JCM, Andrade-da-Costa J, Silva MC, Pinho SS. Glycans as Regulatory Elements of the Insulin/IGF System: Impact in Cancer Progression. Int J Mol Sci. 2017;18(9):E1921.

    Article  PubMed  CAS  Google Scholar 

  80. Lu J. The Warburg metabolism fuels tumor metastasis. Cancer Metastasis Rev. 2019;38(1–2):157–64.

    Article  PubMed  Google Scholar 

  81. Dasgupta S, Rajapakshe K, Zhu B, Nikolai BC, Yi P, Putluri N, Choi JM, Jung SY, Coarfa C, Westbrook TF, et al. Metabolic enzyme PFKFB4 activates transcriptional coactivator SRC-3 to drive breast cancer. Nature. 2018;556(7700):249–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

In memory of our late and beloved co-author Bianca van Diermen Hidle, who all too early became a victim of cancer. We are very much in debt to her legacy, and for her excellent work on the quantitative pathology analysis — provided also for this paper.

Advices on the various statistical analyses given by Prof. Jan Terje Kvaløy at University of Stavanger, Stavanger, Norway are very much appreciated.

We would also like to thank the former department heads Dr. Ottar Bjerkeset and Dr. Kjell H. Kjellevold, who facilitated the study flow in the Department of Surgery and the Department of Pathology, respectively.

Funding

The present study was funded by Marathon Oil, the Folke Hermannsen Foundation, and the Inge Steenslands Foundation, Stavanger, Norway. The funding covered the cost of preOp™, blood chemistry, and hormone analysis. The funding bodies were not involved in the design of the study, data collection, analysis, interpretation of the data, or in writing the manuscript.

Author information

Author notes

  1. Tone Hoel Lende and Marie Austdal shared First authorship.

  2. Emiel A. M. Janssen, Håvard Søiland and Jan PA Baak are equal senior contributor.

Authors and Affiliations

  1. Department of Breast & Endocrine Surgery, Stavanger University Hospital, Helse Stavanger HF, P.O. Box 8100, N-4068, Stavanger, Norway

    Tone Hoel Lende, Siri Lunde & Håvard Søiland

  2. Centre for Cancer Biomarkers CCBIO, Department of Clinical Medicine, Faculty of Medicine and Dentistry, University of Bergen, Jonas Lies vei 87, N-5012, Bergen, Norway

    Tone Hoel Lende & Lars A. Akslen

  3. Department of Research, Stavanger University Hospital, Helse Stavanger HF, P.O. Box 8100, N-4068, Stavanger, Norway

    Marie Austdal & Kristin Jonsdottir

  4. Department of Pathology, Stavanger University Hospital, Helse Stavanger HF, P.O. Box 8100, N-4068, Stavanger, Norway

    Marie Austdal, Anne Elin Varhaugvik, Ivar Skaland, Einar Gudlaugsson, Nina G. Egeland, Emiel A. M. Janssen & Jan P. A. Baak

  5. Department of Circulation and Medical Imaging, Norwegian University of Science and Technology, Trondheim, Norway

    Tone Frost Bathen

  6. Department of Pathology, Helse Møre og Romsdal, Ålesund, Norway

    Anne Elin Varhaugvik

  7. Department of Chemistry, Bioscience and Environmental Technology, University of Stavanger, P.O. Box 8600 Forus, N-4036, Stavanger, Norway

    Nina G. Egeland & Emiel A. M. Janssen

  8. Department of Clinical Science, University of Bergen, Jonas Lies vei 87, N-5012, Bergen, Norway

    Håvard Søiland

  9. Dr. Med. Jan Baak AS, Risavegen 66, N-4056, Tananger, Norway

    Jan P. A. Baak

Authors
  1. Tone Hoel Lende
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Marie Austdal
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Tone Frost Bathen
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Anne Elin Varhaugvik
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Ivar Skaland
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Einar Gudlaugsson
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Nina G. Egeland
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Siri Lunde
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Lars A. Akslen
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Kristin Jonsdottir
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Emiel A. M. Janssen
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Håvard Søiland
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Jan P. A. Baak
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

THL included and operated on all of the patients, built the database, and contributed to statistical analyses and interpretation of data. MA contributed to analyses of the metabolite spectra of all serum and tissue samples and also to the IPA analyses. TFB provided all the metabolomic analyses in her MR-imaging lab and interpreted the results of the metabolic spectra. AEV obtained all of the blood samples and performed the laboratory analyses. IS contributed to the laboratory analyses. EG performed surgical pathological analysis with histological grading and morphological analysis of the tumor. NGE and SL provided detailed information of the micro RNA in the IPA analyses. LAA contributed with scientific support and advice. KJ performed the IPA analyses. EAMJ contributed to the concept of the study, assessed the pathological parameters and scorings, and participated in interpreting the data. HS contributed to the concept of the study, statistical analyses and interpretation of data. JPAB contributed to the concept of the study and the analysis and interpretation of the data. All co-authors contributed to writing the manuscript and gave their final approval of the last version to be published.

Authors’ information

N.A.

Corresponding author

Correspondence to Tone Hoel Lende.

Ethics declarations

Ethics approval and consent to participate

The RCT, which this study is based upon, was approved by the Regional Ethics Committee (Accession number 2015/1445), Norwegian Centre for Research Data (#20984), and The Norwegian Biobank Registry (#2239). An informed consent form was signed by each patient. The trial was retrospectively registered at Clinicaltrials.gov (NCT03886389).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Appendix

Appendix

Table 8 ROC – analysis with ‘Relapse / No relapse’ as dichotomous variable in ER+ patients
Full size table
Table 9 Metabolites correlated to serum insulin C peptide (Pearson’s correlation) for all patients, carbohydrate group, and fasting groups
Full size table
Table 10 Metabolites correlated to serum Insulin Growth Factor Binding Protein 3 (IGFBP3) (Pearson’s correlation) for all patients, carbohydrate group, and fasting groups
Full size table
Table 11 Metabolites correlated to serum Insulin Growth Factor 1 (IGF1) (Pearson’s correlation) for all patients, carbohydrate group, and fasting groups
Full size table
Table 12 Results from Quantitative Metabolite Set Enrichment Analysis
Full size table

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lende, T.H., Austdal, M., Bathen, T.F. et al. Metabolic consequences of perioperative oral carbohydrates in breast cancer patients — an explorative study. BMC Cancer 19, 1183 (2019). https://doi.org/10.1186/s12885-019-6393-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-019-6393-7

Keywords

BMC Cancer

ISSN: 1471-2407

Contact us