- Research article
- Open Access
- Open Peer Review
Metabolomics approaches in pancreatic adenocarcinoma: tumor metabolism profiling predicts clinical outcome of patients
BMC Medicinevolume 15, Article number: 56 (2017)
Pancreatic adenocarcinomas (PAs) have very poor prognoses even when surgery is possible. Currently, there are no tissular biomarkers to predict long-term survival in patients with PA. The aims of this study were to (1) describe the metabolome of pancreatic parenchyma (PP) and PA, (2) determine the impact of neoadjuvant chemotherapy on PP and PA, and (3) find tissue metabolic biomarkers associated with long-term survivors, using metabolomics analysis.
1H high-resolution magic angle spinning (HRMAS) nuclear magnetic resonance (NMR) spectroscopy using intact tissues was applied to analyze metabolites in PP tissue samples (n = 17) and intact tumor samples (n = 106), obtained from 106 patients undergoing surgical resection for PA.
An orthogonal partial least square-discriminant analysis (OPLS-DA) showed a clear distinction between PP and PA. Higher concentrations of myo-inositol and glycerol were shown in PP, whereas higher levels of glucose, ascorbate, ethanolamine, lactate, and taurine were revealed in PA. Among those metabolites, one of them was particularly obvious in the distinction between long-term and short-term survivors. A high ethanolamine level was associated with worse survival. The impact of neoadjuvant chemotherapy was higher on PA than on PP.
This study shows that HRMAS NMR spectroscopy using intact tissue provides important and solid information in the characterization of PA. Metabolomics profiling can also predict long-term survival: the assessment of ethanolamine concentration can be clinically relevant as a single metabolic biomarker. This information can be obtained in 20 min, during surgery, to distinguish long-term from short-term survival.
Pancreatic adenocarcinomas (PAs) are extremely aggressive cancers and have one of the poorest prognoses among all cancers . With an estimated 48,960 new cases in 2015 in the USA, pancreatic cancer is the twelfth most common cancer, representing 3.0% of all the new cancers diagnosed in the USA . The majority of pancreatic cancers are pancreatic ductal adenocarcinomas and are localized in the head of the pancreas [3, 4].
Surgery is the only potentially curative treatment for PA. Pancreatic surgery is associated with significant morbidity and mortality. Indeed the mortality rate, even in highly specialized centers, ranges from 2% to 5%, and morbidity can be as high as 70%, especially in left pancreatic resection [5, 6]. The extension of lymphadenectomy, with its high morbidity and low evidence-based data, is actually a matter for debate [7, 8]. In recent years, a significant shift towards targeted surgical interventions has been proposed, relying on accurate characterization of PA, made possible by preoperative imaging to help evaluate its resectability .
Even though computed tomography (CT) and magnetic resonance imaging (MRI) studies now precisely differentiate resectable or locally advanced versus unresectable PA, there is currently a debate on whether survival can be predicted in resected patients on the basis of imaging. Notably, the vascular contacts may be complex to assess, especially after endoprosthesis or after radio- and/or chemotherapy, where the prognostic value of imaging is still debated . The main prognostic factors for survival after cephalic duodeno-pancreatectomy (CDP) are histological parameters, namely R0 margins, nodal status, and differentiation [11, 12]. Although technical breakthroughs have been achieved in the field of pancreatic surgery, it has nevertheless been shown that the rate of R0 resection is rarely more than 20% and that it consistently impacts survival [13, 14]. Nodal extension depends on the extension of lymphadenectomy, as shown by the impact of invaded node-to-total examined node ratio . However, the benefit of extended lymphadenectomy in PA has not been demonstrated. Therefore, there is a lack of accurate prognostic factors, and currently no tissue biomarkers have been identified to predict long-term survival in patients with pancreatic cancer.
Gross examination and intraoperative extemporaneous microscopic examination are reliable for diagnosis, but cannot predict overall survival.
In recent years, metabolomics, or global metabolite profiling, has been used to investigate the metabolite changes associated with pancreatic cancers [16–21]. Metabolomics is the latest stage of the multi-omics approaches. After genomics, transcriptomics, and proteomics, metabolomics has been generating increasing interest in scientific and medical communities in the past few years, particularly in oncology  and more precisely in pancreatic cancers. Currently, well-recognized tools for metabolomics are nuclear magnetic resonance (NMR) spectroscopy and gas (GC-MS) or liquid chromatography-mass spectrometry (LC-MS). 1H high-resolution magic angle spinning (HRMAS) NMR spectroscopy technology is particularly suitable for the analysis of small samples of intact tissue. This technique avoids the need for chemical extraction procedures or for handling the samples, both of which are required by MS and liquid-state NMR. HRMAS NMR spectroscopy enables identification and quantification of several metabolites from spectra with excellent resolution and signal-to-noise ratio.
Beyond serum markers [16, 19], better characterization of pancreatic tissue would be of particular interest in PA. Consequently, there is a need for accurate tissue biomarkers that could help surgeons distinguish between long-term and short-term survivors. The aims of this study were thus to (1) define the metabolome of pancreatic parenchyma (PP, healthy tissue) and PA, (2) determine the impact of neoadjuvant chemotherapy on healthy tissue (PP) and PA, and (3) by using metabolomics analysis, find metabolic biomarkers associated with long-term survival in patients with PA.
This study included 123 samples obtained from 106 patients retrospectively selected after they had undergone PA resection, from May 2000 to March 2011, in the Department of Visceral Surgery and Transplantation (University Hospitals of Strasbourg, Hautepierre Hospital, Strasbourg, France). These patients fulfilled the following criteria: (1) histological diagnosis of PA, (2) all follow-up patients, (3) only patients with tumor-related deaths, (4) homogeneous adjuvant treatment using the same chemotherapy (gemcitabine) for all patients (and no radiotherapy), and (5) samples of pancreatic tissue collected just after resection and then snap-frozen in liquid nitrogen before storage.
Among the 106 samples obtained from patients with PA, there were:
44 samples from patients who did not receive any neoadjuvant chemotherapy
62 samples from patients who did
Among the 17 samples obtained from PP (healthy tissue), there were:
9 samples from patients who did not receive any neoadjuvant chemotherapy
8 samples from patients who did
Finally, two groups with extremely different prognoses (PA samples) were compared:
Long-term survivors (survival >3 years), 8 samples, no neoadjuvant chemotherapy
Short-term survivors (survival <1 year), 9 samples, no neoadjuvant chemotherapy
For this investigation, the tissue samples were obtained from the tumor bio-bank of Strasbourg University Hospitals. A written informed consent was given by all the included patients.
Tissue sample preparation for HRMAS NMR spectroscopy
All tissue specimens were collected during surgery just after tumor removal and were snap-frozen in liquid nitrogen for –80 °C temperature storage. Then, the sample preparation was performed at a temperature of –20 °C. The amount of tissue used for the HRMAS analysis ranged from 15 mg to 20 mg. Each tissue sample was placed in a 30-μL disposable insert. Next, 8 μL of deuterium oxide with 0.75 weight percent 2,2,3,3-D4-3-(trimethylsilyl) propionic acid was added in every biopsy’s insert in order to get a chemical shift reference for the NMR spectrometer. Finally, inserts were kept at –80 °C until the HRMAS analysis was performed. The insert was placed in a 4-mm ZrO2 rotor just before the HRMAS analysis.
HRMAS NMR data acquisition
All HRMAS NMR spectra were obtained on a Bruker Avance III 500 spectrometer (installed at Hautepierre Hospital, Strasbourg) operating at a proton frequency of 500.13 MHz and equipped with a 4 mm triple resonance gradient HRMAS probe (1H, 13C, and 31P). The temperature was maintained at 277.15 K throughout the acquisition time in order to reduce the effects of tissue degradation during the spectra acquisition. A one-dimensional (1D) proton spectrum using a Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence was acquired with an interpulse delay of 285 μs and an acquisition time of 10 min for each tissue sample (Bruker GmbH, Germany). The number of loops was set to 328, giving the CPMG pulse train a total length of 93 ms. The chemical shift was calibrated to the peak of the methyl proton of l-lactate at 1.33 ppm. In order to confirm resonance assignments in a few representative samples, two-dimensional (2D) heteronuclear experiments (1H – 13C) were also recorded immediately after ending the 1D spectra acquisition. Metabolites were assigned using standard metabolite chemical shift tables available in the literature (Table 1) .
HRMAS NMR data processing and statistical analyses
The HRMAS NMR data processing and the metabolites’ quantification have been previously detailed . Briefly, the region between 7.50 and 0.70 ppm of each 1D HRMAS NMR spectrum was automatically bucketed into integral regions of 0.01 ppm, using AMIX 3.9.14 software (Bruker GmbH, Germany). Once the data set was obtained, it was then exported and analyzed into SIMCA P (version 13.0.3, Umetrics AB, Umeå, Sweden). An orthogonal partial least square-discriminant analysis (OPLS-DA) was performed to analyze the data. The following OPLS-DA model was considered: PP versus PA (both without neoadjuvant chemotherapy). Two measurements of model quality were reported for OPLS-DA: R 2Y and Q2. R2Y > 0.7 and Q2 >0.5 can be considered as a good predictor.
When the population is small, instead of applying OPLS-DA analysis, network analyses using the “algorithm to determine expected metabolite level alterations” using mutual information (ADEMA) are justified . This is why ADEMA has been applied to metabolite quantification values. ADEMA evaluates the changes in groups of metabolites between the case and the control instead of analyzing metabolites one by one. ADEMA includes the metabolic network topology and uses mutual information to find out if those metabolites are biomarkers when considered together, and it can predict the expected change in direction per metabolite when the metabolic network topology is considered. The network was constructed using the Kyoto Encyclopedia of Genes and Genomes [26, 27] and Selway’s work .
The following groups of metabolites were compared related to involved metabolic pathways:
Choline, phosphorylcholine, glycerophosphocholine, ethanolamine
Glucose, ascorbate, glycine, glutamate
Glutamate, glutamine, glycine
Glutamate, glycine, creatine
The metabolites were quantified using the PULCON method, which is a very accurate quantification method (with a very limited percentage of error ).
The repetition time used for this study (2 s) and the total acquisition time (10 min) are hence a good compromise in order to quantify the metabolites by HRMAS NMR spectroscopy of intact tissue. However, under these conditions, we cannot see the whole of the metabolites, but only the freer parts of them. Furthermore, the latter are underestimated about 20% in comparison with studies performing tissue extractions (data not shown). This method is widely used in the literature.
Metabolite quantification was performed using an external reference standard of lactate (3 μmol), scanned under the same analytical conditions as the tissue samples. Spectra were normalized according to sample weight. Peaks of interest were automatically defined by an in-house program using MATLAB 7.0 (MathWorks, Natick, MA, USA). Peak integration was then compared to the one obtained with the lactate reference and was corrected according to the number of protons. Only well-defined peaks with no overlapping in the 1D CPMG spectra were selected for quantification (Tables 2, 3, 4, and 5). Quantification results were expressed as nanomoles per milligram of tissue.
Continuous variables are expressed as mean ± standard deviation (SD). The Mann-Whitney U test was performed to compare the metabolites’ concentrations of (1) PP and PA (both without neoadjuvant chemotherapy), (2) PP without and PP with neoadjuvant chemotherapy, (3) PA without and PA with neoadjuvant chemotherapy, and (4) long-term and short-term survival in patients with PA (Tables 2, 3, 4, and 5). The Mann-Whitney U tests were performed using R software .
Receiver operating characteristic (ROC) curves and Kaplan-Meier curves were used to perform a survival analysis and therefore to evaluate the clinical utility of metabolite quantification in the long-term survival characterization (R software ).
All the spectra obtained from the 123 analyzed specimens were of high quality without any obvious evidence of tissue necrosis. A total of 31 metabolites were identified within the range of 7.50–0.70 ppm from the spectra obtained from all pancreatic tissue samples (Table 1). Among the 31 identified metabolites, only 18 metabolites were quantified: only well-defined peaks with no overlapping in the 1D CPMG spectra were selected for quantification.
The representative 1D HRMAS NMR CPMG spectra of PP (healthy tissue) and PA samples (both without neoadjuvant chemotherapy) are shown in Figs. 1a and 2a. Some discriminant metabolites were highlighted using the Mann-Whitney U test. Choline (p = 0.0014), ethanolamine (p = 0.0226), glycerol (p = 0.0037), glycine (p = 0.0005), lactate (p = 0.0006), and taurine (p = 0.0021) were statistically significant between PP and PA (both without any neoadjuvant chemotherapy) (Table 2). Metabolomic profiles of PP and PA were clearly separated by a bi-component OPLS-DA (R2Y = 0.82; Q2 = 0.69) (Fig. 3). A higher concentration of myo-inositol and glycerol was shown in PP tissue samples. By contrast, a higher level of glucose, ascorbate, ethanolamine, lactate, and taurine was revealed in PA tissue samples.
Impact of neoadjuvant chemotherapy on PP and PA
Seventeen samples from PP (healthy tissue) were included in this model: 8 samples received a neoadjuvant chemotherapy, while the 9 others did not (Fig. 1). No discriminant metabolites were found using the Mann-Whitney U test (Table 3). The network analysis showed that a decreased level of succinate, aspartate, taurine, phosphorylcholine, glucose, tyrosine, lactate, and glutamine was predicted in PP samples from patients with neoadjuvant chemotherapy. Moreover, a higher level of threonine and glycine was predicted in PP tissue samples from patients with neoadjuvant chemotherapy. The other metabolites were predicted to be equivalent between the two groups (Fig. 4).
Among the 106 samples of PA, 62 received neoadjuvant chemotherapy, while the 44 others did not (Fig. 2). Some discriminant metabolites were highlighted using the Mann-Whitney U test. Aspartate (p = 0.0017) was statistically significant between PA samples from patients who received neoadjuvant chemotherapy and those who did not (Table 4). As shown in Fig. 5, the network analysis revealed a decreased level of lactate, which was predicted in PA samples from patients with neoadjuvant chemotherapy. Moreover, higher levels of ethanolamine, tyrosine, phenylalanine, aspartate, glucose, glycerol, succinate, glycine, glutamate, glutamine, and creatine were predicted in PA tissue samples from patients with neoadjuvant chemotherapy. The other metabolites were predicted to be equivalent between the two groups.
For some of these analyses, the Mann-Whitney U test did not yield any significant differences, whereas ADEMA reports some changes for those metabolites. Thus, ADEMA appears more effective, as it compares different groups of metabolites, uses mutual information, and does not require a large population of samples.
Long-term versus short-term survival in patients with PA
Patients’ characteristics are described in Table 6. The representative 1D HRMAS NMR CPMG spectra of long-term and short-term survivors are shown, respectively, in Fig. 6a and b. Seventeen samples were studied: 8 long-term survivors and 9 short-term survivors (both without any neoadjuvant chemotherapy). Some discriminant metabolites were highlighted using the Mann-Whitney U test. Choline (p = 0.0150), ethanolamine (p = 0.0078), lactate (p = 0.0360), and phenylalanine (p = 0.0055) were statistically significant between long-term and short-term survival in patients with PA (Table 5).
The network analysis showed higher levels of glucose, ascorbate, and taurine; this was predicted in long-term survivors. Moreover, decreased levels of choline, ethanolamine, glycerophosphocholine, phenylalanine, tyrosine, aspartate, threonine, succinate, glycerol, lactate, glycine, glutamate, glutamine, and creatine were predicted in long-term survivors. Phosphorylcholine was predicted to be equivalent between the two groups (Fig. 7).
Statistical significant differences were highlighted when we focused on the comparison between long-term and short-term survival patients. Only the use of ethanolamine as a single screening test showed a higher accuracy in distinguishing long-term from short-term survivors. As shown in Fig. 8, samples with a decreased level of ethanolamine had a high probability of being assigned to long-term survivors. All the patients showing a low level of ethanolamine were long-term survival patients. The area under the curve (AUC) was 0.861 ± 0.101. The optimal ethanolamine threshold was 0.740 nmol/mg when distinguishing long-term from short-term survivors. With this threshold point, sensitivity and specificity were, respectively, 77.80% and 75.00% (Fig. 8a). The predictive positive value was 47.00% and the negative predictive value was 52.90%. A Kaplan-Meier analysis showed that the overall survival probabilities were significantly higher in patients with low tumor ethanolamine concentrations compared to those with high tumor ethanolamine concentrations (Fig. 8b).
To the best of our knowledge, this is the first study that has evaluated the metabolome of intact tissues for PP and PA.
Preoperative assessment of resectability of PA is not yet an adequate way of predicting survival. However, given the morbidity and mortality of CDP, a better evaluation of the balance between risks and benefits is a complex quest. Indeed, serum markers such as carcinoembryonic antigen (CEA) or carbohydrate antigen (CA 19-9) are poorly correlated to long-term results and should not be used for contraindicating the only potentially curative treatment of this often fatal disease. Up to now only histological parameters have enabled surgeons to retrospectively evaluate the potential benefit of CDP, but these parameters are available only after surgery. Moreover, none of these parameters is really specific for predicting the necessity of pancreatic resection. In order to evaluate the benefit of pancreatic resection, long-surviving patients were compared to short-term survivors. Of note, among the long-term survivors, two of them showed an infra-millimetric margin, and three had histological invasion of the portal vein (Table 6). Thus, the use of currently described predictive factors should not lead to contraindicating resection in patients who can otherwise tolerate surgery . The use of a more specific marker tightly linked to the pathophysiology of pancreatic cancer would be of particular interest.
HRMAS NMR provides new insights into the relationships between metabolic pathways and pancreatic cancers. This method allows for the identification of cell membranes and phospholipid metabolism, cellular energy production via neoglucogenesis, the tricarboxylic acid (TCA) cycle, and oxidative stress. The majority of PAs have KRAS mutation (90% with activating mutations in this oncogene) [32–34].
Our results about PP’s metabolome show a higher level of lactate, glucose, phosphorylcholine, taurine, aspartate, lactate, glutamine, and succinate in patients who received neoadjuvant chemotherapy (gemcitabine). Mutations in TCA cycle enzymes are known to promote cancer development and growth. For example, mutations in succinate dehydrogenase (SDH), fumarate hydratase (FH), and isocitrate dehydrogenase 1 and 2 (IDH1, IDH2) can be cited [35, 36]. Moreover, these mutations have been associated with distinct cancer subsets and different patient prognoses [37, 38]. Mutations in TCA cycle enzymes have also been associated with measurable changes in the levels of metabolites. Nevertheless, these mutations have not been a major feature of PA. Reprogrammed cellular metabolism has increasingly become an obvious field of research for PA. Research has also focused on the alterations that involve the TCA cycle and mutant KRAS-induced tumor cell dependencies for glucose, glutamine, and extracellular protein.
Our results about PA’s metabolome show a higher level of lactate in patients who did not receive any neoadjuvant chemotherapy. Elevated expression of lactate dehydrogenase A (LDHA) is highlighted in breast and colorectal cancers [39, 40] and has been recently identified in pancreatic cancers . Some studies have already shown the different steps where LDHA is involved (tumor initiation, maintenance, and progression) [42, 43]. Given the elevated glycolysis in PA, leading to increased lactate production, the ability to utilize lactate could provide an additional advantage to PA cells . Indeed, the inhibition of glycolysis through suppressing LDHA expression by RNA interference decreased the growth of PA cells . Our study may provide new targets to treat PA. When patients who received neoadjuvant chemotherapy and patients who did not were compared, neoadjuvant chemotherapy seemed to have an inhibiting effect on LDHA expression. A higher level of lactate in patients who did not receive any neoadjuvant chemotherapy showed that survival of PA was highly dependent on LDHA activity in a hypoxic environment. Then, this higher level of lactate could be linked to improved responses to neoadjuvant chemotherapy. Indeed, lactate might be a predictive marker for assessing the response of tumor cells to neoadjuvant chemotherapy. Finally, lactate could be correlated with tumor neoadjuvant chemotherapy in predicting responses to this therapy. Elevated levels of lactate are prognostic biomarkers for poor survival in several cancers. Our findings suggest that lactate could be an important marker for screening the efficiency of cancer treatments.
A decreased level of glucose has also been shown in PA without neoadjuvant chemotherapy. As described in the work of Koong et al.  and Guillaumont et al. , areas within PA tumors are hypoxic, and this has implications for the cellular metabolism. Hypoxia promotes tumor growth by particularly stimulating glycolysis.
When PP and PA, both without chemotherapy, were compared, some discriminant metabolites were highlighted, particularly lactate and taurine. Our study shows a higher level of lactate and taurine in patients with PA. These results are consistent with the work of Wang et al. . The work of Nishiumi et al. has also underlined elevated levels of lactate in the patients’ serum . As described in the work of Wang et al., higher levels of taurine have also been reported in the literature, for several types of cancer. As previously explained, this increased lactate level in PA could come from hypoxia. Our study is in accordance with other studies that have found higher levels of taurine in different types of cancer, probably due to apoptosis [49, 50]. Even so, further studies will be needed to understand the complex biology of that particular type of cancer. If confirmed by other studies, this could deepen our knowledge of pancreatic cancer pathogenesis and might also lead to the identification of new targets for diagnosis, early detection, imaging, or even future therapeutic options.
Our results about PA’s metabolome in short-term survivors are in accordance with other studies that found a higher level of choline in several cancers. The role of choline kinase α (Chk-α) in malignant transformation and progression in several cancers has been well described in the literature. Increased levels and activity of Chk-α have been observed in human breast , colorectal , lung [52, 53], prostate , ovarian , and more recently in endometrial  and pancreatic cancers .
The altered choline metabolism we observed is possibly due to an overexpression of Chk-α. Increased levels of choline could be due to an increase in the membrane activity, particularly due to enhanced cellular proliferation. This activity is increased in inflammatory or tumoral conditions. In the most malignant cases, this membrane activity is much higher still. Increased levels of choline are also due to malignant transformation, and particularly due to an overexpression of Chk-α. Hypoxia may also affect choline phosphorylation through effects on Chk-α and transporters ; in turn, choline phosphorylation may be affected by hypoxia through effects on Chk-α and transporters. In some studies, particularly in prostate cancer , increased levels in choline have been described, leading to the development of 18F-fluorocholine positron emission tomography (PET) imaging in order to detect the increased uptake and the phosphorylation of the tracer. Consequently, our results are in agreement with observations in the literature and could represent the biological substrate and justification to the use of 18F-fluorocholine PET imaging in pancreatic cancers [59–61]. Metabolic targets in choline phospholipid metabolism may provide new therapeutic options for PAs that have severely limited options. Moreover, a decreased amount of succinate was shown by the network analysis in long-term survivors compared to short-term survivors, suggesting a decreased activity of the TCA cycle. PAs from long-term survivors also show decreased levels of threonine, aspartate, glycerol, glutamate, and glutamine compared to those from short-term survivors. Glutamine has been particularly studied for its role in cancer metabolism because it appears to be required for the growth of many types of tumors . Decreased levels of creatine and lactate were shown in long-term survival patients too. Other studies showed the relationship between lower levels of creatine, lactate, and choline and overall survival . Understanding the different metabolic links within pancreatic cancer is a promising approach to identifying novel prognostic markers (long-term survival) and therapeutic programs in patient care .
Short-term survivors showed higher levels of choline, glycerophosphocholine, ethanolamine, and fatty acids, as depicted in Fig. 6b (respectively, numbers 19, 25 and 21, 22). Fatty acids are a major factor in the growth of tumor cells. Some explorations have begun in order to study the PA metabolism and more precisely the role of individual complex fatty acids. Indeed, as detailed in the work of Guo et al.  and Zadra et al. , fatty acids can be both pro-tumorigenic and anti-tumorigenic, as described for various cancers, making their biology difficult to explain. Reducing the levels of certain fatty acids seems to be important for PA. But, we have to determine first which fatty acids are cytotoxic for tumor cells and which fatty acids provide the tumor with metabolic substrates . Focusing on ethanolamine, the survival analysis showed that survival was longer for patients with lower tumoral ethanolamine concentrations. The assessment of ethanolamine concentration can be clinically relevant as a single metabolic biomarker for distinguishing long-term survivors from short-term survivors in patients with PA.
The current study demonstrates that metabolomics profiling may provide prognostic information in patients with pancreatic cancer. Research has usually focused on enzymatic steps within the TCA cycle, as it potentially influences the progression of disease, as well as on alterations of the phospholipid metabolism within the choline/ethanolamine membrane.
Only very few patients are needed to build a model that can predict oncological outcome in pancreatic cancer very accurately; this fact alone should help in promoting our technique.
We acknowledge some limitations to the present study. First, the number of patients is limited, particularly regarding the long-term survivors included into our study. Pancreatic cancers are extremely aggressive and have one of the poorest prognoses among all cancers. Thus, very few patients get a chance of long-term survival, and even less so without any neoadjuvant chemotherapy. Second, further studies should take into account other elements that affect the survival of patients. We believe these data are preliminary and should be validated in further series. We encourage others to validate the findings and to carry out multicenter studies. Indeed, in order to include co-factors, the patient population needs to be significantly enhanced (because PAs have one of the poorest prognoses among all cancers). Third, we acknowledge that a comparative effectiveness study should be performed in real time in order to evaluate usual serum markers versus HRMAS NMR spectroscopy before making any definitive conclusions. Fourth, due to the limited number of patients and the preliminary character of our results, it is difficult to draw a conclusion about the level of ethanolamine. For the moment, there is no specific concept. Fifth, although we acknowledge that the assessment of ethanolamine concentration can be clinically relevant as a single metabolic biomarker for distinguishing long-term survivors from short-term survivors in patients with PA, this result should be put into perspective. Indeed, sensitivity and specificity are respectively 77.80% and 75%. Lastly, this study was retrospective and may involve some bias that could have been unaccounted for.
Metabolomics analysis could be validated as an intraoperative discriminant method for distinguishing healthy tissues from PA tissues. This could deepen our knowledge of PA metabolism and may also lead to the identification of new targets for diagnosis, imaging, or future therapeutic options. If these results are confirmed in further studies, it is expected that the role of intraoperative HRMAS NMR spectroscopy could then be evaluated in the setting of PA. This approach, which we call “metabolomics-guided surgery,” could help surgeons to extend the excision if necessary. Since HRMAS NMR spectroscopy enables rapid characterization of intact tissue, it could also be used as an intraoperative method. HRMAS analysis only takes 20 min. Data analysis is also very quick (<10 min). Furthermore, the cost is < $50 per sample.
In conclusion, as we are able to distinguish PP from PA, we could imagine using this technique to analyze several samples collected from the excision cavity in addition to tissue tumor samples. This last step could help surgeons in the detection of residual tumor cells in the excision cavity and the control of margins. This approach could be used in clinical routine for prediction of long-term survival in patients with PA; indeed, this information can be obtained in 20 min during surgery. Finally, due to the high morbidity and mortality during surgery, we could also imagine using this technique before surgery, with the use of endoscopic or percutaneous biopsy, although these remain invasive techniques. On the whole, our own technique could prove useful and have a positive impact on patient care.
High-resolution magic angle spinning
Nuclear magnetic resonance
Orthogonal partial least square-discriminant analysis
Principal component analysis
Sousa CM, Kimmelman AC. The complex landscape of pancreatic cancer metabolism. Carcinogenesis. 2014;35:1441–50.
National Cancer Institute. SEER Cancer Statistics Factsheets: Pancreas Cancer. National Cancer Institute Bethesda, MD. http://seer.cancer.gov/statfacts/html/pancreas.html. Accessed 13 May 2016.
Vincent A, Herman J, Schulick R, Hruban RH, Goggins M. Pancreatic cancer. Lancet. 2011;378:607–20.
Li D, Xie K, Wolff R, Abbruzzese JL. Pancreatic cancer. Lancet. 2004;363:1049–57.
McMillan MT, Christein JD, Callery MP, Behrman SW, Drebin JA, Hollis RH, et al. Comparing the burden of pancreatic fistulas after pancreatoduodenectomy and distal pancreatectomy. Surgery. 2016;159:1013–22.
Vollmer Jr CM, Lewis RS, Hall BL, Allendorf JD, Beane JD, Behrrman SW, et al. Establishing a quantitative benchmark for morbidity in pancreatoduodenectomy. Ann Surg. 2015;261:527–36.
Jang JY, Kang MJ, Heo JS, Choi SH, Choi DW, Park SJ, et al. A prospective randomized controlled study comparing outcomes of standard resection and extended resection, including dissection of the nerve plexus and various lymph nodes, in patients with pancreatic head cancer. Ann Surg. 2014;259:656–64.
Pedrazzoli S, DiCarlo V, Dionigi R, Mosca F, Pederzoli P, Pasquali C, et al. Standard versus extended lymphadenectomy associated with pancreatoduodenectomy in the surgical treatment of adenocarcinoma of the head of the pancreas: a multicenter, prospective, randomized study. Lymphadenectomy Study Group. Ann Surg. 1998;228:508–17.
Andersson R, Vagianos CE, Williamson RCN. Preoperative staging and evaluation of resectability in pancreatic ductal adenocarcinoma. HPB (Oxford). 2004;6:5–12.
Katz MH, Fleming JB, Bhosale P, Varadhachary G, Lee JE, Wolff R, et al. Response of borderline resectable pancreatic cancer to neoadjuvant therapy is not reflected by radiographic indicators. Cancer. 2012;118:5749–56.
Wang F, Gill AJ, Neale M, Puttaswamy V, Ganadha S, Pavlakis N, et al. Adverse tumor biology associated with mesenterico-portal vein resection influences survival in patients with pancreatic ductal adenocarcinoma. Ann Surg Oncol. 2014;21:1937–47.
Schnelldorfer T, Ware AL, Sarr MG, Smyrk TC, Zhang L, Qin R, et al. Long-term survival after pancreatoduodenectomy for pancreatic adenocarcinoma: is cure possible? Ann Surg. 2008;247:456–62.
Yekebas EF, Bogoevski D, Cataldegirmen G, Kunze C, Marx A, Vashist YK, et al. En bloc vascular resection for locally advanced pancreatic malignancies infiltrating major blood vessels: perioperative outcome and long-term survival in 136 patients. Ann Surg. 2008;247:300–9.
Delperro JR, Bachellier P, Regenet N, Le Treut YP, Paye F, Carrere N, et al. Pancreaticoduodenectomy for pancreatic ductal adenocarcinoma: a French multicentre prospective evaluation of resection margin in 150 evaluable specimens. HPB (Oxford). 2014;16:20–33.
Vuarnesson H, Lupinacci RM, Semoun O, Svrcek M, Julié C, Balladur P, et al. Number of examined lymph nodes and nodal status assessment in pancreaticoduodenectomy for pancreatic adenocarcinoma. Eur J Surg Oncol. 2013;39:1116–21.
Ritchie SA, Akita H, Takemasa I, Eguchi H, Pastural E, Nagano H. Metabolic system alterations in pancreatic cancer patient serum: potential for early detection. BMC Cancer. 2013;13:416.
Urayama S. Pancreatic cancer early detection: expanding higher-risk group with clinical and metabolomics parameters. World J Gastroenterol. 2015;21:1707–17.
Fontana A, Mathur A, Copetti M, Di Gangi IM, Mazza T, Tavano F, et al. Development of a metabolites risk score for one-year mortality risk prediction in pancreatic adenocarcinoma patients. Oncotarget. 2016;7:8968–78.
LaConti JJ, Laiakis EC, Mays AD, Peran I, Kim SE, Shay JW, et al. Distinct serum metabolomics profiles associated with malignant progression in the KrasG12D mouse model of pancreatic ducta adenocarcinoma. BMC Genomics. 2015;16 Suppl 1:S1.
Zhang G, He P, Tan H, Budhu A, Gaedcke J, Ghadimi BM, et al. Integration of metabolomics and transcriptomics revealed a fatty acid network exerting growth inhibitory effectsin human pancreatic cancer. Clin Cancer Res. 2013;19:4983–93.
He XH, Li WT, Gu YJ, Yang B, Deng Hw YY, et al. Metabonomic studies of pancreatic cancer response to radiotherapy in a mouse xenograft model using magnetic resonance spectroscopy and principal components analysis. World J Gastroenterol. 2013;19:4200–8.
Marengo E, Robotti E. Biomarkers for pancreatic cancer: recent achievements in proteomics and genomics through classical and multivariate statistical methods. World J Gastroenterol. 2014;20:13325–42.
Wishart DS, Jewison T, Guo AC, Wilson M, Knox C, Liu Y, et al. HMDB 3.0—The Human Metabolome Database. Nucleic Acids Res. 2013;41:D801–7.
Imperiale A, Elbayed K, Moussallieh FM, Reix N, Piotto M, Bellocq JP, et al. Metabolomic profile of the adrenal gland: from physiology to pathological. Endocr Relat Cancer. 2013;20:705–16.
Cicek AE, Bederman I, Henderson L, Drumm ML, Ozsoyoglu G. ADEMA: an algorithm to determine expected metabolite level alterations using mutual information. PLoS Comput Biol. 2013;9:e1002859.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42:D199–205.
Selway ZZ. Metabolism at a glance. 3rd ed. Malden (MI): Blackwell Publishing; 2014.
Wider G, Dreier L. Measuring protein concentrations by NMR spectroscopy. J Am Chem Soc. 2006;128:2571–6.
R Development Core Team. R: A language and environement for statistical computing. 2016. http://www.R-project.org.
Mierke F, Hempel S, Distler M, Aust DE, Saeger HD, Weitz J, et al. Impact of portal vein involvement from pancreatic cancer on metastatic pattern after surgical resection. Ann Surg Oncol. 2016;23:730–6.
Almoquera C, Shibata D, Forrester K, Martin J, Arnheim N, Perucho M. Most human carcinomas of the exocrine pancreas contain mutant c-K-ras genes. Cell. 1988;53:549–54.
Uemura T, Hibi K, Kaneko T, Takeda S, Inoue S, Okochi O, et al. Detection of K-ras mutations in the plasma DNA of pancreatic cancer patients. J Gastroenterol. 2004;39:56–60.
Löhr M, Klöppel G, Maisonneuve P, Lowenfels AB, Lüttges J. Frequency of K-ras mutations in pancreatic intraductal neoplasias associated with pancreatic ductal adenocarcinoma and chronic pancreatitis: a meta-analysis. Neoplasia. 2005;1:17–23.
Desideri E, Vegliante R, Cririolo MR. Mitochondrial dysfunctions in cancer: genetic defects and oncogenic signaling impinging on TCA cycle activity. Cancer Lett. 2015;356:217–23.
Morin A, Letouze E, Gimenez-Roqueplo AP, Favier J. Oncometabolites-driven tumorigenesis: from genetics to targeted therapy. Int J Cancer. 2014;135:2237–48.
Srinivasan R, Ricketts CJ, Sourbier C, Linehan WM, et al. New strategies in renal cell carcinoma: targeting the genetic and metabolic basis of disease. Clin Cancer Res. 2015;21:10–7.
Hai Y, Parsons DW, Jin G, McLendon R, Rasheed A, Yuan W, et al. IDH1 and IDH2 mutations in gliomas. N Engl J Med. 2009;360:765–73.
Thornburg JM, Nelson KK, Clem BF, Lane AN, Arumugam S, Simmons A, et al. Targeting aspartate aminotransferase in breast cancer. Breast Cancer Res. 2008;10:R84.
Sanchez-Tena S, Alcarraz-Vizan G, Marin S, Torres JL, Cascante M. Epicatechin gallate impairs colon cancer cell metabolic productivity. J Agric Food Chem. 2013;61:4310–7.
Lu QY, Zhang L, Yee JK, Go VL, Lee WN. Metabolic consequences of LDHA inhibition by epigallocatechin gallate and oxamate in MIA PaCa-2 pancreatic cancer cells. Metabolomics. 2015;11:71–80.
Le A, Cooper CR, Gouw AM, Dinavahi R, Maitra A, Deck LM, et al. Inhibition of lactate dehydrogenase A indices oxidative stress and inhibits tumor progression. Proc Natl Acad Sci USA. 2010;107:2037–42.
Fantin VR, St-Pierre J, Leder P. Attenuation of LDH-A expression uncovers a link between glycolysis, mitochondrial physiology, and tumor maintenance. Cancer cells. 2006;9:425–34.
Zhao D, Zou SW, Liu Y, Zhou X, Mo Y, Wang P, et al. Lysine-5 acetylation negatively regulates lactate dehydrogenase A and is decreased in pancreatic cancer. Cancer cell. 2006;23:464–76.
Koong AC, Mehta VK, Le QT, Fisher GA, Terris DJ, Brown JM, et al. Pancreatic tumors show high levels of hypoxia. Int J Radiat Oncol Biol Phys. 2000;48:919–22.
Guillaumont F, Leca J, Olivares O, Lavaut MN, Vidal N, Berthezene P, et al. Strengthened glycolysis under hypoxia supports tumor symbosis and hexosamine biosynthesis in pancreatic adenocarcinoma. Proc Natl Acad Sci USA. 2013;110:3919–24.
Wang AS, Lodi A, Rivera LB, Izquierdo-Garcia JL, Firpo MA, Mulvihill SJ, et al. HR-MAS MRS of the pancreas associated with early pancreatic cancer. NMR Biomed. 2014;27:1361–70.
Nishiumi S, Shinohara M, Ikeda A, Yoshie T, Hatano N, Kakuyama S, et al. Serum metabolomics as a novel diagnostic approach for pancreatic cancer. Metabolomics. 2010;6:518–28.
Opstad KS, Bell BA, Griffiths JR, Howe FA. Taurine: a potential marker of apoptosis in gliomas. Br J Cancer. 2009;100:789–94.
Zhang X, Tu S, Wang Y, Xu B, Wan F. Mechanism of taurine-induced apoptosis in humas colon cancer cells. Acta Biochim Biophys Sin. 2014;46:261–72.
Ramirez de Molina A, Gutierrez R, Ramos MA, Silva JM, Bonilla F, Sanchez JJ, et al. Increased choline kinase activity in human breast carcinomas: clinical evidence for a potential novel antitumor strategy. Oncogene. 2002;21:4317–22.
Ramirez de Molina A, Rodriguez-Gonzalez A, Gutierrez R, Martinez-Pineiro L, Sanchez J, Bonilla F, et al. Overexpression of choline kinase is a frequent feature in human tumor-derived cell lines and in lung, prostate and colorectal human cancers. Biochem Biophys Res Commun. 2002;296:580–3.
Ramirez de Molina A, Sarmentero-estrada J, Belda-Iniesta C, Taron M, Ramirez de Molina V, Cejas P, et al. Expression of choline kinase alpha to predict outcome in patients with early-stage non-small-cell lung cancer: a retrospective study. Lancet Oncol. 2007;8:889–97.
Iorio E, Mezzanzanica D, Alberti P, Spadaro F, Ramoni C, D’Ascenzo S, et al. Alterations of choline phospholipid metabolism in ovarian tumor progression. Cancer Res. 2005;65:9369–76.
Trousil S, Lee P, Pinato DJ, Ellis JK, Dina R, Aboagye EO, et al. Alterations of choline phospholipid metabolism in endometrial cancer are caused by choline kinase alpha overexpression and a hyperactivated deacylation pathway. Cancer Res. 2014;74:6867–77.
Penet MF, Shah T, Bharti S, Krishnamachary B, Artemov D, Mironchik Y, et al. Metabolic imaging of pancreatic ductal adenocarcinoma detects altered choline metabolism. Clin Cancer Res. 2015;21:386–95.
Glunde K, Shah T, Winnard Jr PT, Raman V, Takagi T, Vesuna F, et al. Hypoxia regulates choline kinase expression through hypoxia-inducible factor-1 alpha signaling in a human prostate cancer model. Cancer Res. 2008;68:172–80.
Vali R, Loidl W, Pirich C, Langesteger W, Beheshti M. Imaging of prostate cancer with PET/CT using 18F-fluorocholine. Am J Nucl Med Mol Imaging. 2015;5:96–108.
DeGrado TR, Reiman RE, Price DT, Wang S, Coleman RE. Pharnacokinetics and radiation dosimetry of 18 F-fluorocholine. J Nucl Med. 2002;43:92–6.
Witney TH, Alam IS, Turnton DR, Smith G, Carroll L, Brickute D, et al. Evaluation of deuterated 18 F- and 11C-labeled choline analogs for cancer detection by positron emission tomography. Clin Cancer Res. 2012;18:1063–72.
Challapalli A, Sharma R, Hallet WA, Kozlowski K, Carroll L, Brickute D, et al. Biodistribution and radiation dosimetry of deuterium-substituted 18 F-fluoromethyl-[1,2-2H4] choline in healthy volunteers. J Nucl Med. 2014;55:256–63.
Dr W, Thompson CB. Glutamine addiction: a new therapeutic target in cancer. Trends Biochem Sci. 2010;35:427–33.
Yokota H, Guo J, Matoba M, Higashi K, Tonami H, Nagao Y. Lactate, choline, and creatine levels measured by vitro 1H-MRS as prognostic parameters in patients with non-small-cell lung cancer. J Magn Reson Imaging. 2007;25:992–9.
Guo JY, Karsil-Uzunbas G, Mathew R, Aisner SC, Khamphorst JJ, Strhecker AM, et al. Autophagy suppresses progression of K-ras-induced lung tumors to oncocytomas and maintains lipid homeostasis. Genes Dev. 2013;27:1447–61.
Zadra G, Photopoulos C, Loda M. The fat side of prostate cancer. Biochim Biophys Acta. 2013;1831:1518–32.
The laboratory staff of the Tumor Bio-bank (Centre de Ressources Biologiques) of the Strasbourg University Hospitals are gratefully acknowledged for their technical assistance. We gratefully acknowledge Dr. Karim Elbayed for the MATLAB scripts that he developed allowing for the quantification of metabolites.
Availability of data and materials
Our data sets are now presented within Additional file 1 (Values of metabolites'quantification) (XLSX format).
SB wrote the manuscript, prepared the tissue samples, and performed HRMAS-NMR data acquisition, spectra analysis, and statistical analysis. FF wrote the manuscript and provided patient surgical management. AI reviewed the manuscript. AEC performed network analysis, CH performed statistical analyses, and GA performed histopathological analyses. PB performed patient surgical management and reviewed the manuscript. IJN provided experimental design, data analysis, and manuscript review. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
A written informed consent was obtained from all the included patients. For this investigation, the tissue samples were obtained from the Tumor Bio-bank (Centre de Ressources Biologiques) of the Strasbourg University Hospitals.
Ethics approval and consent to participate
The Ethics Committee of Strasbourg (Comité de Protection des Personnes “Est IV”) approved the study (Registration number: 09/39a). The declaration at the Ministerial level and the authorization from the Tumor Bio-Bank (Centre de Ressources Biologiques) of the Strasbourg University Hospitals correspond to the following number: AC 2008-438/DC 2009-1016.
Written informed consent was obtained from all the included patients. For this investigation, the tissue samples were obtained from the Tumor Bio-bank (Centre de Ressources Biologiques) of the Strasbourg University Hospitals.
Values of metabolites'quantification. (XLSX 80 kb)