|Home | About | Journals | Submit | Contact Us | Français|
Breast cancer is a clinically heterogeneous disease, which necessitates a variety of treatments and leads to different outcomes. As an example, only some women will benefit from chemotherapy. Identifying patients who will respond to chemotherapy and thereby improve their long‐term survival has important implications to treatment protocols and outcomes, while identifying non responders may enable these patients to avail themselves of other investigational approaches or other potentially effective treatments. In this study, serum metabolite profiling was performed to identify potential biomarker candidates that can predict response to neoadjuvant chemotherapy for breast cancer. Metabolic profiles of serum from patients with complete (n = 8), partial (n = 14) and no response (n = 6) to chemotherapy were studied using a combination of nuclear magnetic resonance (NMR) spectroscopy, liquid chromatography–mass spectrometry (LC–MS) and statistical analysis methods. The concentrations of four metabolites, three (threonine, isoleucine, glutamine) from NMR and one (linolenic acid) from LC–MS were significantly different when comparing response to chemotherapy. A prediction model developed by combining NMR and MS derived metabolites correctly identified 80% of the patients whose tumors did not show complete response to chemotherapy. These results show promise for larger studies that could result in more personalized treatment protocols for breast cancer patients.
Breast cancer, although histologically similar, is clinically a very heterogeneous disease, which results in a range of treatment effectiveness and outcomes (Paik et al., 2004). Neoadjuvant chemotherapy can significantly benefit breast cancer patients; however, the varied response to such therapy means that a significant proportion of the patient population is subjected to ineffective treatment while at the same time being exposed to the therapy's toxicities (Wolmark et al., 2001). Pathologic complete response (pCR), which is defined as the disappearance of the invasive cancer cells in the breast after chemotherapy, is used to evaluate patient response and is strongly associated with improved long‐term survival rates (Bear et al., 2004; Fisher et al., 1998; Kuerer et al., 1999). Unfortunately, less than 30% of patients overall show complete response to neoadjuvant chemotherapy (Jones and Smith, 2006). An ability to predict response to chemotherapeutic agents should enable development of personalized treatment protocols, improving survival rates and reducing unnecessary exposure of patients to toxic drugs.
Research focused on finding useful molecular or clinical predictors of pCR to neoadjuvant chemotherapy in breast cancer is relatively sparse. Imaging studies, such as magnetic resonance imaging (MRI) (Padhani et al., 2006) and scintimammography (Marshall et al., 2005; Sciuto et al., 2002) were proposed to predict pathological responses to neoadjuvant chemotherapy, but they are somewhat limited by low sensitivity combined with high costs. High levels of MUC‐1 antigen (CA 15.3) in pre‐treatment serum and its fall after chemotherapy can predict responses as well (Al‐Azawi et al., 2006) but many patients do not exhibit elevation of this marker before treatment and hence it is not helpful for such patients (Kurebayashi et al., 2003). Approaches using genomics and immunohistochemistry have been explored to find serum and tissue biomarkers (Bear et al., 2003; Guarneri et al., 2006; Rouzier et al., 2005; van't Veer et al., 2002). It has been shown that gene signatures such as HER2 overexpression/amplification and lack of ER expression were associated with pCR and certain neoadjuvant chemotherapy regimens (Chen et al., 2011; Gianni et al., 2005; Thuerigen et al., 2006). Other molecular markers such as tumor RNA (Parissenti et al., 2010), glucose‐regulated protein (GRP78) (Lee et al., 2011) and hormone receptors (Van Poznak et al., 2002; Wang et al., 2009) have also been identified as potential predictors of pCR. However, suboptimal performance is a major issue that limits their wide applicability. Circulating tumor cells (CTC) have also been established as providing outcome predictions from particular therapies; however CTCs can be detected in less than 30% of early stage breast cancer patients, which limits their clinical applicability (Hayes and Smerage, 2008).
As an alternative approach for biomarker discovery, metabolomics (or metabolite profiling) enables identification of small‐molecule metabolites in biofluids and tissues that are sensitive to altered pathology (Lindon et al., 2004; Nicholson et al., 1999; Nicholson and Wilson, 2003). High‐throughput analytical techniques of nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry (MS) combined with multivariate statistical analyses provide information on a large number of metabolites, including those that have altered levels between healthy subjects and patients with various diseases including cancer (Lanza et al., 2010; Pan and Raftery, 2007; Zhang et al., 2012). So far, the metabolomics based approaches have been used in a large variety of applications, including early disease detection, drug response, toxicity and nutritional studies, and basic systems biology (Clayton et al., 2006; Gowda et al., 2008; Griffin and Kauppinen, 2007; Serkova and Niemann, 2006; Sreekumar et al., 2009; Zhang et al., 2008). Compared with other biomarker discovery approaches for breast cancer, metabolomics provides a strong link between genotype and phenotype, and may provide some insight into oncogenesis. Also, once established, tests based on metabolic profiles are relatively inexpensive, rapid and can be automated (Spratlin et al., 2009).
A growing number of metabolomics studies are contributing toward an improved understanding of breast cancer, and these advances have been reviewed (Claudino and Quattrone, 2007; Gowda et al., 2008; Oakman et al., 2011a). Many of the studies have focused on identifying altered metabolic levels in breast cancer cells or tissues, and relating these changes to their associated metabolic pathways (Cheng et al., 1998; Sitter et al., 2002; Whitehead and Kieber‐Emmons, 2005; Bathen et al., 2007; Sitter et al., 2006; Yang et al., 2007). A very recent study using metabolic profiling of numerous human cancer cell lines found a high correlation between breast cancer (and other cancer) proliferation and the glycine biosynthetic pathway (Jain et al., 2012). Previously, differences between normal and metastatic mammary epithelial cell lines, including up‐regulation of fatty acid synthesis and alterations in glycolysis, the TCA cycle and others, were detected using 13C stable isotopic label tracing by 2D NMR and GC–MS methods (Yang et al., 2007). Breast cancer tumors could be separated from non‐involved tissues based on intensities from spectra generated by high‐resolution magic angle spinning (HR‐MAS) NMR spectroscopy with a sensitivity of 83% and a specificity of 100%. Some metabolites such as choline and glycine were found to be significantly upregulated in tumors larger than 2 cm (Sitter et al., 2006). In another NMR study, a multivariate statistical model based on 67 urinary metabolites successfully identified all the breast cancer patients with high specificity (93%) (Slupsky et al., 2010). Breast cancer prognostic factors, such as estrogen and progesterone receptor status, could be predicted by HR‐MAS NMR based metabolomics on tissue samples (Giskeodegard et al., 2010). Metastatic breast cancer patients could be differentiated from early stage patients with 72% prediction accuracy using serum samples detected by NMR based metabolomics (Oakman et al., 2011b). For identifying breast cancer recurrence, a predictive model built on 11 biomarkers detected by combining NMR and two‐dimensional gas chromatography–mass spectrometry (GC × GC–MS), provided 86% sensitivity and 84% specificity (Asiago et al., 2010).
In this study, we use a metabolomics approach to predict the response to chemotherapy in the neoadjuvant setting. Serum samples from 28 patients obtained before preoperative chemotherapy have been studied using a combination of NMR, liquid chromatography–mass spectrometry (LC–MS) and multivariate statistics methods. Four metabolites that were identified from NMR and MS methods are well correlated with pCR. A statistical model built based on these metabolites predicts pCR with high sensitivity and specificity.
For NMR experiments, deuterium oxide (D2O, 99.9% D) was purchased from Cambridge Isotope Laboratories, Inc. (Andover, MA). Sodium azide (NaN3) and the sodium salt of trimethylsilylpropionic acid‐d4 (TSP) were purchased from Sigma–Aldrich (Milwaukee, WI). For LC–MS experiments, two internal standards (tridecanoic acid and chlorophenylalanine) and linolenic acid were purchased from Sigma–Aldrich (analytical grade, St. Louis, MO). HPLC‐grade methanol and acetic acid were purchased from Fisher Scientific (Pittsburgh, PA).
The patients were enrolled into the study between 2005 and 2008 (Table 1; see also Supplementary Table S1), and were treated according to current guidelines for neoadjuvant therapy. Those patients with locally advance breast cancer and eligible for neoadjuvant chemotherapy were enrolled in this study. All patients were recruited and treated at the Department of Gynecology and Obstetrics, University of Tuebingen. The criteria for patient selection included: 1) signed informed consent and willingness to comply with requirements; 2) at least 18 years old; 3) histopathologically verified lesions of locally advanced breast cancer; 4) indication for neoadjuvant chemotherapy; 5) consent to neoadjuvant chemotherapy; 6) negative result for a serum human chorionic gonadotropin pregnancy test before start of the study for women of childbearing potential; 7) ability to consent to study participation. Patient's study data were collected in case report forms (CRF). For each included patient a CRF was completed. Documents which identify the patient (e.g. the patient identification log and the signed informed consent) were maintained in confidence by the investigator. More detailed information about the sample collection can be found in the supplementary materials.
Summary of clinicopathological characteristics of the breast cancer patients studied.
The response to neoadjuvant chemotherapy in patients with locally advanced breast cancer was based on the correlation of magnetic resonance imaging (MRI), two dimensional (2D)/three dimensional (3D) ultrasound (US) and mammography (MG) with histopathology (study code BCD‐001, funded by the BMBF, Germany). Based on their response, patients were categorized as having pathologic complete response (pCR), partial response (PR) or stable disease (SD). Here, pCR is defined as the disappearance of all tumor deposits; PR indicates a reduction of tumor volume exceeding 50%; while tumor reduction less than 50% or increase of volume up to 25% is scored as stable disease (SD) (Neubauer et al., 2008). The response of the tumors to the neoadjuvant chemotherapy was evaluated pathologically by classifying the regressive changes using a semiquantitative scoring system from 0 to 4 (0 = no effect, 1 = resorption and tumor sclerosis, 2 = minimal residual invasive tumor [<0.5 cm], 3 = residual noninvasive tumor only, 4 = no tumor detectable) according to the tumor regression grading described by Sinn et al. (1994). A consultant pathologist (U. Vogel) blinded to clinical outcome reviewed all paired biopsy and surgical specimens. Samples were de‐identified in compliance with a protocol approved by the IRB committees at both Eberhard‐Karls‐University of Tuebingen and Purdue University.
Serum samples from all 28 breast cancer patients with baseline samples available before neoadjuvant chemotherapy were obtained. Clinical and histopathological data are listed in Table 1 and Supplementary Table S1. Of these patients, 8 belonged to pCR group, 6 to SD group and 14 to PR group (Table 1 and S1). Blood (2 × 8 mL) was drawn 45 min before the start of the chemotherapy. The collected blood was allowed to clot for 45 min at room temperature and centrifuged for 10 min at 3000 rpm at room temperature; the upper serum phase was then isolated, aliquoted and frozen at −80 °C until further use.
Patients were treated according to the design of GeparQuattro including 4 three weekly EC cycles consisting of Epirubicin (90 mg/m2), with simultaneous administration of Cyclophosphamide (600 mg/m2). Subsequently, patients received Docetaxel (100 mg/m2) in three weekly intervals. Patients with Her2/neu positive tumors were also treated with Trastuzumab (6 mg/kg, i.v.) every three weeks, starting with a loading dose of 8 mg/kg i.v. on day 1 of the first EC‐cycle.
Frozen serum samples were thawed and 530 μL was mixed with 5 μL 0.1% NaN3 solution. A 60 μL 0.5 mM TSP solution was utilized as an internal standard. 1H 1D NMR experiments were performed using the CPMG (Carr–Purcell–Meiboom–Gill) pulse sequence coupled with water presaturation on a Bruker Avance‐500 spectrometer equipped with a TXI gradient cryoprobe. 128 scans with 16 k time domain data points were collected using a spectral width of 6000 Hz. An exponential weighting function corresponding to 1.0 Hz line broadening was applied to the time domain data before Fourier transformation. The spectra were then phased, baseline corrected and referenced to alanine (δ = 1.48 ppm) using Bruker Topspin 3.0 software.
Frozen serum samples (100 μL) were thawed, and protein was precipitated by adding 200 μL methanol containing two internal standards, tridecanoic acid and chlorophenylalanine, used to monitor extraction efficiency. The solutions were centrifuged at 13,200 rpm for 30 min; the resulting supernatants were dried under vacuum and reconstituted in 25 μL methanol/water (1:1). A pooled sample, which was a mixture of small random volumes from all 28 samples, was extracted using the same procedure as above. This sample was used as a quality control (QC) and was analyzed after every 10 patient samples. LC–MS analysis was performed using an Agilent 6520 HPLC‐QTOF system (Agilent Technologies, Santa Clara, CA) consisting of an Agilent 1200 SL liquid chromatograph coupled with a time‐of‐flight (TOF) mass spectrometer. The reconstituted serum samples were gradient‐eluted at 600 μL/min using (A) 0.2% acetic acid in water and (B) 0.2% acetic acid in methanol (2% B to 98% B in 13 min, 98% B for 6 min). A 3 μL sample aliquot was injected onto a 2.1 × 50 mm Agilent Zorbax Extend‐C18 1.8 μm particle column with a 2.1 × 30 mm Agilent Zorbax SB‐C 8 3.5 μm particle guard column heated to 60 °C. Electrospray ionization (ESI) was used in positive mode. The MS interface capillary was maintained at 325 °C, with a sheath gas flow of 9 L/min. The spray voltage for positive ion injection was 4.0 kV. The instrument was scanned over a range of 50–1000 m/z. Agilent MassHunter Workstation LC–TOF and QTOF Acquisition software (B.02.01) was used for automatic peak detection and mass spectrum deconvolution.
The three groups of patients have similar age distributions. Although the pCR cohort has the highest average age (49.5) and the SD cohort has the lowest (43.7), both groups have almost the same age range (37–60 for pCR and 37–59 for SD). In terms of cancer stage and grade, most cancer patients in all three cohorts were in the T1–T2 stage and G2 grade, which represent 75% of pCR patients (n = 6), 71% of PR (n = 10) and 83% of SD patients (n = 5), respectively. All three cohorts contain both pre and post‐menopausal patients. In the pCR and PR groups, pre and post‐menopausal patients are evenly distributed (pre:post = 4:4 and 7:6 for pCR and PR, respectively), while SD cohorts have more pre‐menopausal patients (pre:post = 4:1 for SD). All three groups have diverse distributions of ER/PR/Her2 statuses. Most of the samples (n = 24) are from ER or PR positive women. And almost half of the samples (n = 13) are Her2 positive.
From the 1H NMR spectra, 27 spectral regions with metabolites peaks were identified and integrated after local baseline correction, and normalized to the total sum of all the metabolites. In order to identify distinguishing metabolite biomarkers, the integral for each metabolite was statistically compared using the Student's t‐test between different groups of patients.
The LC–MS data were processed using Agilent's MassHunter Qualitative Analysis software (version B.03.01). A list of ion intensities for each peak detected was generated, matching m/z and retention time (RT) for each ion. Agilent MassHunter Workstation Mass Profiler Professional software (version B.02.00) was then used for compound peak alignment and removal of any peaks with missing values (ion intensity = 1) in more than one sample from any group; 115 metabolites passed this filter. Internal standard peaks were also removed. Finally, the Agilent Formula Database (Agilent 2010) was used for metabolite identification by matching the accurate mass spectrum to a database of compounds. The Student's t‐test was performed between pCR and SD samples, and 9 metabolites with low p‐values (<0.1) (Supplementary Table S2) were selected as biomarker candidates for further statistical analysis.
The statistical analysis for this metabolomics study was performed retrospectively after sample collection. We performed comparisons by computing p‐values between all 3 pairs of the responding groups (SD, PR, pCR) to find any difference or separation between these groups.
All 27 metabolite peak integrals from NMR and 9 LC–MS metabolites with p < 0.1 were imported into Matlab software (Mathworks) installed with the PLS toolbox (Eigenvector Research, Inc., version 4.0) for partial least squares (PLS) analysis to obtain clustering information and identify outliers. The same software was also used for partial least squares discriminant analysis (PLS‐DA) modeling. Leave‐one‐out cross validation (CV) was chosen, and the number of latent variables were selected according to the root mean square error of CV procedure. The R statistical package (version 2.8.0) was used to generate receiver operating characteristics (ROC) curves and box‐and‐whisker plots.
In order to visualize the intrinsic grouping of samples and identify outliers, PLS analysis was performed using the 27 NMR derived metabolite regions, and the results of this analysis are shown in the Supplementary Figure S1. The results showed that the SD patient group has one outlier, which was omitted from further analysis; the new PLS score plot thus obtained is depicted in Figure 1(a). Both pCR and SD groups were separated along the first two latent variables (LV1 and LV2), while the PR group lies between pCR and SD. Similar analysis performed using the 9 metabolites with p < 0.1 obtained from LC–MS data was also performed, and the results are shown in Figure 1(b). As observed for the NMR data, metabolites derived from LC–MS also separated pCR from the SD group of patients. However, again the clustering for the PR group was not as well defined.
PLS score plot based on (a) 27 metabolites detected by NMR and (b) 9 metabolites detected by LC–MS for all the patients. In (a), one outlier was deleted. The LV label on the two axes stands for latent variable. pCR: pathologic complete response; ...
Comparison of the NMR data between different groups of patients using the Student's t‐test showed four metabolites to be statistically significant (p < 0.05) (Table 2). One‐way ANOVA was not used here since the objective is to see the difference between each pair of groups. These p‐values indicate that levels of three metabolites, isoleucine, threonine and glutamine, were significantly different between pCR and SD groups and the levels of two metabolites, threonine and glutamine, were different between PR and SD. Only one metabolite, histidine, differed significantly between pCR and PR. The LC–MS data showed that the most statistically differentiating compounds found were long‐chain lipids or fatty acids. The most interesting of these, linolenic acid, was validated using a pure, commercially obtained compound. This metabolite separated pCR from SD samples perfectly, which makes combining linolenic acid with other LC–MS detected markers unnecessary for model building. Statistical analysis shows linolenic acid to be significantly different between pCR and SD groups (p < 0.01). Performance of these metabolites was evaluated individually using receiver operating characteristic (ROC) curves and box‐and‐whisker plots. Figure 2 shows box‐and‐whisker plots for the five metabolites (4 from NMR and 1 from LC–MS) and Supplementary Figure S2 shows ROC curves along with area under the ROC (AUROC) values. The concentration distribution for all the metabolites except histidine showed a consistent trend from pCR to PR to SD; while threonine, glutamine and linolenic acid increased, isoleucine decreased. Further, four metabolites, threonine, glutamine, isoleucine and linolenic acid each showed a minimum AUROC of 0.80 in classifying pCR and SD patients, while only linolenic acid distinguished pCR from both SD and PR patients with an AUROC of greater than 0.80 (Supplementary Figure S2). As mentioned above, linolenic acid distinctly separated pCR from SD with an AUROC = 1 (Table 3).
Individual metabolite box‐and‐whisker plots for the different groups of patients: (a) threonine; (b) glutamine; (c) isoleucine; (d) histidine; and (e) linolenic acid.
Summary of NMR metabolites having low p‐values.
Summary of data for linolenic acid detected by LC–MS.
Further analysis focused on evaluating the performance of the metabolites in combination. First, we combined the 3 NMR markers, threonine, glutamine and isoleucine that distinguished pCR and SD, performed PLS‐DA and built a statistical classification model. Leave‐one‐out cross validation was used to reduce overfitting and estimate model accuracy. Figure 3 shows prediction results using this model for different patient groups. The AUROC for distinguishing pCR and SD was 0.81. However, as with the individual metabolites, the performance for distinguishing pCR from the other two groups, PR and SD, was relatively poor (AUROC = 0.72).
(a) Predictions results for the PLS‐DA model based on the 3 metabolite markers (threonine, isoleucine and glutamine) detected by NMR; (b) ROC curve for pCR vs. SD using the cross‐validated predicted class values (AUROC = 0.81); ...
Analysis was also performed combining three NMR derived markers (threonine, glutamine and isoleucine) with LC–MS detected linolenic acid. The results from the statistical analysis are shown in Figure 4. The model provides 100% selectivity and 80% sensitivity for the prediction of pCR vs. SD with an AUROC of 0.95, which is significantly better compared to the classification offered by the 3 marker model (Figure 3). PR samples were not well classified by either model.
(a) Prediction results for the PLS‐DA model based on combining isoleucine, glutamine, and threonine detected by NMR and linolenic acid detected by LC–MS; (b) ROC curve for pCR vs. SD using the cross‐validated predicted class values ...
In this study, we present a new metabolomics approach for predicting the outcome of breast cancer neoadjuvant chemotherapy. By combining two different analytical platforms, NMR and LC–MS, three amino acids and one fatty acid are shown to be highly correlated with pCR. Although the clinical and histopathological parameters of the patients and tumors in each response group are quite heterogeneous, the serum samples still fall into distinct clusters. In order to investigate the impact of different subtypes of breast cancer on the performance of four markers we used for model building, p‐values from the Student's t‐test were obtained by comparing samples from different ER/PR/Her2 status for pCR, PR, SD. All of these sample data are summarized in Supplementary Table S3. No p‐value less than 0.05 can be found, indicating that none of the 4 markers is significantly altered by any of these parameters in the data set.
The resulting prediction model has high sensitivity and specificity. Considering the challenges involved in predicting the outcome of almost any cancer treatment, this promising method might open a new approach for selecting the use of certain cancer treatments or even provide better personalized treatment options. Interestingly, however, as seen from the comparison of the clinical data and clustering of the patient samples in the statistical analysis based on blood metabolites (Table 1, Supplementary Table S1 and Figure 4), the derived metabolite biomarkers provide a strong response to chemotherapy and do not indicate an association with any sub types such as ER, PR or Her2 status. These results suggest that altered metabolic pathways associated with the derived biomarkers are unconnected to ER, PR or Her2 status. However, studies on larger patient cohorts are required to substantiate these findings and identify metabolic link between the pathologically different subtypes, if any.
The performance of linolenic acid detected by LC–MS analysis was particularly striking, and showed an AUROC of 0.86 for the classification of pCR from PR and SD, and 1.0 for the classification of pCR from SD (Figure S2). In addition, the combination of MS and NMR detected metabolites provides a number of advantages including enabling identification of additional distinguishing metabolites that provide better insights into the cellular metabolism and provides potentially a more robust model that can perform effectively when tested with a larger cohort of patients. Four metabolites derived from the combination of NMR and MS distinguished the three breast cancer patient groups, pCR, PR and SD with good performance. The predictive model thus developed by combining the four metabolites correctly separated pCR from SD patients with 100% specificity and 80% sensitivity (Figure 4). Such a high degree of sample classification achieved through the serum derived metabolites is particularly noteworthy at a time when the prediction of the response to breast cancer neoadjuvant chemotherapy continues to be challenging. We found that PR samples show large variance in the values for most of the markers we identified. Some PR patients tend to be predicted as pCR; while others cluster with SD samples. This effect may be due to the heterogeneous results of chemotherapy for these patients who are classified between SD and pCR. Pathologic response rates are surrogates of patient outcome. It will be important in a larger cohort of patients to determine if the PRs who were predicted as pCR have a better clinical outcome than those predicted to be SD.
To the best of our knowledge, this is the first study to combine NMR and MS to identify altered metabolites in the serum for use as predictors of response to neoadjuvant chemotherapy. Altered levels of the four metabolites identified in this study, threonine, glutamine, isoleucine and linolenic acid that predict response to neoadjuvant chemotherapy represent changes in the metabolic activity of several pathways. Information about the association of specific metabolic pathways with response prediction remains to be clearly understood. However, numerous studies have reported altered levels of blood serum/plasma amino acids in breast cancer patients and those results match with the association of the amino acids found in the present study (Asiago et al., 2010; Lai et al., 2005; Zhang et al., 2011). For example, a study using blood plasma showed increased concentration of threonine in breast cancer (Kubota et al., 1992). In addition, investigations of breast cancer cells showed increased levels of threonine, under hypoxic conditions, and such an increase is attributed to impaired protein synthesis due to hypoxia (Weljie et al., 2011). Increased glutamine levels in blood plasma have been reported for breast cancer, in addition to esophageal and liver cancers (Gao et al., 2009; Proenza et al., 2003; Zhang et al., 2011). Glutamine is an essential nutrient for the growth of cancer cells, and the catabolism of glutamine is a major route for the proliferation of mammalian cells. Many types of cancer cells are shown to depend on high rates of glutamine uptake and metabolism (Coles and Johnstone, 1962; Eagle, 1955; Wise et al., 2008). An association of isoleucine that we find in the present study has not been reported previously for breast cancer. However, its altered levels in blood plasma are known in liver and colorectal cancers (Georgiannos et al., 1995; Watanabe et al., 1984). Identification of the fatty acid, linolenic acid, is also consistent with several earlier investigations, which showed significant upregulation of fatty acids in breast cancer (Baron et al., 2004; Boros et al., 2002; Yang et al., 2007).
Thus, although, three of the four distinguishing metabolites, threonine, glutamine and linolenic acid, have previously been shown to have links with breast cancer, the mechanism of their down‐regulation in patients that responded to chemotherapy (pCR group) compared to those that do not (SD group) is not well understood. The results, however, suggest that metabolites in the serum of breast cancer patients are indicators of tumor/host metabolism and that they can predict both sensitivity and resistance to chemotherapy a priori. This is consistent with previous reports that indicate some hypoxic cancer cells exhibit higher levels of threonine and that high levels of fatty acids found in malignant breast cancer cells are also associated with clinically aggressive tumor resistance behavior (Kuhajda, 2000; Weljie et al., 2011; Zeng et al., 2010). Moreover, it is now known that some cancer cells display addiction to glutamine. Glutamine is a major oxidizable substrate for tumor cells, which require glutamine for survival (Fuchs and Bode, 2006; Wise and Thompson, 2010; Yuneva et al., 2007). Hence, the higher concentration of glutamine observed for patients (SD group) that did not respond to chemotherapy appears to favor tumor survival and growth, and offer resistance to chemotherapy.
We present a prediction model for the outcome of breast cancer neoadjuvant chemotherapy based on metabolic profiling studies that combine NMR and LC–MS methods. A combination of four metabolites, three detected by NMR (threonine, glutamine and isoleucine), and one by MS (linolenic acid) distinguish groups of patients with no, partial or complete response. While this pilot study was focused on a small patient cohort, it clearly indicates that several blood‐based metabolite markers are sensitive to response, and that the approach is promising for predicting the response to chemotherapy. Studies using larger patient cohorts will be needed to substantiate the results and determine any association of the derived metabolite markers with pathologically different sub types. Validation of these results using a larger and independent sample cohort, and identification of additional metabolite markers, will provide better insights into the pathology at the molecular level that lead to different response outcomes for the three groups of patients. In addition, considering its strong performance as a biomarker in the present study, linolenic acid and possibly other fatty acids might be of particular interest for further validation studies. This approach, which clearly differentiates patients that respond to drugs from those that do not, may provide a useful tool for a non‐invasive prognosis of the treatment outcome.
Daniel Raftery has an executive role at Matrix‐Bio, Inc.
The following is the supplementary material related to this article:
This work was supported by Purdue University and the Purdue Research Foundation. SC is a grantee of the Breast Cancer Research Foundation and her efforts are dedicated to the memory of Evelyn H. Lauder. The authors also thank Dr. Shanaiah Murthy at Matrix‐Bio, Inc, and the Purdue Interdepartmental NMR Facility (PINMRF) for their assistance with NMR experiments.
Supplementary material related to this article can be found at http://dx.doi.org/10.1016/j.molonc.2012.10.003.
Wei Siwei, Liu Lingyan, Zhang Jian, Bowers Jeremiah, Gowda G.A. Nagana, Seeger Harald, Fehm Tanja, Neubauer Hans J., Vogel Ulrich, Clare Susan E., Raftery Daniel, (2013), Metabolomics approach for predicting response to neoadjuvant chemotherapy for breast cancer, Molecular Oncology, 7, doi: 10.1016/j.molonc.2012.10.003.