Introduction {#Sec1}
============

Desmoid tumors (aggressive fibromatosis) are locally invasive soft tissue tumors that account for 3% of soft tissue sarcomas^[@CR1]^. Current treatments include, surgery, and radiation^[@CR2]^, but treatments frequently have temporary results.  For instance, between 23--50% patients have recurrence after surgical resection^[@CR3],[@CR4]^. The first line of treatment is often watch-and-wait because of this probability of recurrence and lack of targeted therapeutics^[@CR5],[@CR6]^.

Mutations of the beta-catenin encoding gene (*CTNNB1*) were identified with approximately 85% of sporadic tumors by Sanger sequencing, previously^[@CR7]^, but next-generation sequencing estimates that these mutations may account for 90--95% of desmoid tumor cases^[@CR8]^. Adenomatous polyposis coli (*APC*) gene mutations has been linked to the genetic disorder, familial adenomatous polyposis (FAP), in which at least 25% of patients develop desmoid tumors.  Both of these genes are involved in the Wnt signaling pathway, which directs cell proliferation and adhesion, cell polarity, and cell fate determination^[@CR9]^. Two mutations of the beta-catenin (T41A and S45F) account for 70% of the sporadic tumor formations^[@CR10]^, and the S45F mutation has been correlated to an increase likelihood of local recurrence^[@CR10]^. Additionally, in a small investigation, the S45F genotype has shown poor therapeutic response with the clooxygenase-2 (COX-2) selective inhibitor, Meloxicam^[@CR11]^. Besides the genotyping data, additional immunohistochemical analysis yielded no definitive results among typical tumor and oncology markers^[@CR12]^. Commercially available drug libraries can be used to screen compounds targeting desmoid tumor cells, and an ongoing screening procedure is identifying compounds that inhibited desmoid cellular growth but not normal fibroblast viability^[@CR13]^. All of these compounds are either approved for use in patients or have the potential to be rapidly developed for patient care.  Two compounds found in this ongoing screening, Dasatinib and FAK Inhibitor 14, target proteins that directly interact with the beta-catenin upstream of the Wnt pathway.  Dasatinib is Bcr-Abl and Src family tyrosine kinase inhibitor^[@CR14]^, while FAK inhibitors target the intersecting integrin and receptor tyrosine kinase signal transduction pathways^[@CR15]^.

In recent years, the 'rediscovery' of how critical metabolic processes regulate tumor etiology and progression is identifying new prognostic factors and therapeutic approaches in a variety of tumors.  To the best of our knowledge, this approach has not been applied to beta-catenin driven mesenchymal tumors.  The importance of these findings is the opportunity to leverage the divergences from normal cellular functions into targetable features for novel drug development.  To this end, metabolomics has become a valuable scientific tool in deciphering differential metabolic patterns and the factors that drive them.  Cell line models are essential for biomarker discovery research^[@CR16]^, as are establishing clinically relevant cell-based assays for target identification.  Metabolomics has been instrumental, for instance, in understanding differences in metabolomics profiles and perturbed pathways between primary cell lines associated with triple negative breast cancer^[@CR17]^, ovarian cancer, and colon cancer^[@CR18]^. Additionally, it is possible to utilize metabolomics to determine the mechanism of action of novel therapeutics.  One recent Nature Communications article thoroughly described how broad spectrum metabolomics is utilized when evaluating suspension and adherent cell lines in response to known kinase inhibitors^[@CR19]^.

The goal of this project was to understand the underlying metabolomic differences in paired normal fibroblast and desmoid cells from affected patients and normal fibroblast cells from unaffected patients.  Previous breast cancer research has indicated that it is important to understand the differences between normal tissue from affected and unaffected patients^[@CR20],[@CR21]^. Moreover, the metabolic response to two novel desmoid drugs identified by the Desmoid Collaboration for a Cure was investigated.  Broad spectrum metabolomics has not been applied to desmoid tumor research, and metabolic profiles differences in between desmoid tumor cells, normal fibroblast cells, and unaffected normal fibroblast cells were found readily.  Additionally, there were differences identified in the metabolomics profiles associated with the two beta-catenin mutations, T41A and S45F, and in their response to two desmoid therapeutics, Dasatinib and FAK Inhibitor 14, identified from a drug screen.  These data reflect differences in desmoid tumor biology, and the metabolic response to two chemotherapeutics provides important information for understanding how these effective therapies for desmoid tumors modify biological pathways. 

Results {#Sec2}
=======

Metabolomics profiles distinguishes untreated cell lines {#Sec3}
--------------------------------------------------------

To understand the differences between desmoid tumor, normal fibroblast, and unaffected fibroblast, broad spectrum metabolomics was used to investigate the untreated cell lines.  Unsupervised multivariate analysis (PCA) and supervised analysis (OPLS-DA) differentiated the untreated cell lines.  As shown in Fig.  [1A](#Fig1){ref-type="fig"}, the PCA shows the cell types between unaffected and affected, in addition to tumor and normal.  The unaffected cell line was separated from the affected cell lines on the 1^st^ component of the PCA, whereas the matched desmoid tumor and normal cell lines were along the 2^nd^ PCA component.  Supervised analysis (OPLS-DA) was used to distinguish tumor from normal (Fig.  [1B](#Fig1){ref-type="fig"}), normal from unaffected (Fig.  [1C](#Fig1){ref-type="fig"}), and tumor from unaffected (Fig.  [1D](#Fig1){ref-type="fig"}).  In these comparisons, it is clear that the metabolomic profiles of the two tumor and paired normal cell lines differ from the unaffected.  While the predictive 1st component in these three supervised analysis differentiated the cell type, the orthogonal component is dominated, again, by the genotypes.  The NMR bins associated with defining the differences between the cell lines is listed in Table [1](#Tab1){ref-type="table"}, which included library matched metabolites which were found to be VIP ≥ 1.0, p \< 0.1, or FC \> 2.  It was found that the lipids were greatest in the unaffected cell line compared to the affected.  The desmoid tumor cell lines had higher lipids than their matched normal.  The bins with identified amino acids (alanine, leucine, isoleucine, valine, glutamate, glutamine, serine, threonine, and tyrosine) experienced a higher fold change in the unaffected cell line compared to the affected cell lines.  Additionally, these same amino acids were noted also to be higher fold change in the tumor compared to the matched normal cells.  The exception, however, was aspartate and glycine, where the bins associated with this amino acid had a higher fold change in the normal cells than the tumor cells.  It was also interesting that asparagine and phenylalanine were differentiated in the affected and unaffected, while only methionine differentiated in the normal and unaffected binned NMR profiles.  Dimethylamine was found to be lowest in the tumor cell lines.  Metabolites were semi-quantified, and Table [2](#Tab2){ref-type="table"} shows the metabolites that were found to be p \< 0.1 or FC \> 2.  Interestingly, glucose 1-phosphate was significantly higher in both affected cell lines compared to the unaffected cell line, with the normal and tumor median concentrations 210 and 140 µM, respectively. Figure 1Multivariate analysis comparing untreated desmoid, normal, and unaffected. (** A**) In the unsupervised analysis, the 141 (T41A) and 219 (S45F) cell lines from affected subjects separate from the unaffected cell line from ATCC on the primary axis \[5 components, R2X = 0.79, Q2 (cum) = 0.95\]. (** B**) The supervised analysis of tumor and normal fibroblast \[1 predictive + 1 orthogonal components, R2X = 0.71, R2Y = 0.82, Q2 (cum) = 0.71\]. (** C**) The supervised analysis of normal and unaffected \[1 predictive + 1 orthogonal components R2X = 0.97, R2Y = 1.0, Q2 (cum) = 1.0\]. (** D**) The supervised analysis of tumor and unaffected (D) \[1 predictive + 1 orthogonal components, R2X = 0.89, R2Y = 0.993, Q2 (cum) = 0.95\].Table 1The library matched NMR bins that differentiate the untreated desmoid, tumor, and unaffected cell lines where the p values were calculated with the Wilcoxon Rank Sum Test (normal vs tumor) and exact Wilcoxon Rank Sum Test (tumor vs unaffected and normal vs unaffected) (VIP ≥ 1.0, p \< 0.1, or FC \> 2).Library-Matched MetaboliteChemical ShiftTumor vs NormalNormal vs UnaffectedTumor vs UnaffectedCell Line DifferenceVIPp-value\*FC\*\*VIPp-value\*FC\*\*VIPp-value\*FC\*\*Lipids\[0.50.. 0.52\]0.40.024−2.00.50.024−2.6Unaffected \> Normal \> TumorLipids\[0.52.. 0.55\]0.40.024−1.80.50.048−2.0Unaffected \> Normal \> TumorLipids\[0.55.. 0.57\]0.30.024−1.40.40.048−1.3Unaffected \> Tumor \> NormalLipids\[0.57.. 0.59\]0.40.024−1.50.50.024−1.7Unaffected \> Normal \> TumorLipids\[0.59.. 0.62\]0.40.024−1.50.50.024−1.7Unaffected \> Normal \> TumorLipids\[0.62.. 0.64\]0.40.024−1.30.50.024−1.3Unaffected \> Tumor \> NormalLipids\[0.64.. 0.66\]0.40.024−1.40.50.024−1.5Unaffected \> Normal \> TumorLipids\[0.66.. 0.68\]0.50.024−1.50.60.024−1.5Unaffected \> Normal \> TumorLipids\[0.68.. 0.70\]0.50.024−1.40.60.024−1.4Unaffected \> Normal \> TumorLipids\[0.70.. 0.72\]0.50.024−1.40.60.024−1.4Unaffected \> Normal \> TumorLipids\[0.72.. 0.74\]0.50.024−1.40.60.024−1.4Unaffected \> Tumor \> NormalLipids\[0.74.. 0.77\]0.60.024−1.30.60.024−1.3Unaffected \> Tumor \> NormalLipids\[0.77.. 0.79\]0.60.024−1.30.60.024−1.3Unaffected \> Normal \> TumorLipids\[0.79.. 0.81\]0.50.048−1.30.60.024−1.3Unaffected \> Normal \> TumorLipids\[0.81.. 0.83\]0.50.095−1.20.50.048−1.3Unaffected \> Normal \> TumorLipids\[0.83.. 0.89\]1.10.024−1.21.00.548−1.1Unaffected \> Tumor \> NormalLipids\| Glycocholate\[0.89.. 0.91\]0.70.024−1.2Unaffected \> Normal\[0.91.. 0.97\]2.30.2551.12.20.024−1.42.10.024−1.3Unaffected \> Tumor \> NormalValine \| Isoleucine\[0.97.. 1.02\]1.60.3971.11.40.024−1.31.30.024−1.3Unaffected \> Tumor \> NormalValine\[1.02.. 1.07\]1.50.938−1.01.30.024−1.41.50.024−1.4Unaffected \> Normal \> TumorUnknown\[1.07.. 1.10\]0.40.024−1.20.60.024−1.3Unaffected \> Normal \> TumorUnknown\[1.10.. 1.12\]0.40.024−1.20.50.024−1.3Unaffected \> Normal \> TumorUnknown\[1.12.. 1.14\]0.40.024−1.20.50.024−1.2Unaffected \> Normal \> TumorLipids\[1.14.. 1.16\]0.40.024−1.20.50.024−1.2Unaffected \> Normal \> TumorLipids\[1.16.. 1.18\]0.40.095−1.10.40.024−1.2Unaffected \> Normal \> TumorIsoleucine \| 3-Hydroxyisovalerate\[1.23.. 1.27\]0.80.024−1.3Unaffected \> TumorIsoleucine\[1.27.. 1.30\]0.70.024−1.3Unaffected \> TumorLactate \| Threonine\[1.30.. 1.34\]1.70.3201.11.80.024−1.41.70.024−1.2Unaffected \> Tumor \> NormalLipids\[1.34.. 1.36\]0.50.024−1.20.40.024−1.2Unaffected \> Tumor \> NormalLipids\[1.36.. 1.39\]0.50.0931.10.60.024−1.2Unaffected \> Tumor \> NormalLipids\[1.39.. 1.44\]1.00.024−1.20.80.024−1.1Unaffected \> Tumor \> NormalAlanine\| Lipids\[1.44.. 1.49\]2.00.3971.01.80.024−1.41.80.024−1.3Unaffected \> Tumor \> NormalLipids\[1.49.. 1.55\]1.20.3971.01.20.024−1.41.20.024−1.4Unaffected \> Tumor \> NormalLipids\[1.55.. 1.57\]0.60.024−1.30.60.024−1.2Unaffected \> Tumor \> NormalLipids\[1.57.. 1.59\]0.50.024−1.30.50.024−1.2Unaffected \> Tumor \> NormalLipids\[1.59.. 1.61\]0.50.024−1.20.40.048−1.2Unaffected \> Tumor \> NormalLeucine\[1.61.. 1.64\]0.70.024−1.30.60.024−1.2Unaffected \> Tumor \> NormalLeucine\[1.64.. 1.66\]0.70.024−1.30.70.024−1.2Unaffected \> Tumor \> NormalLeucine\[1.66.. 1.68\]0.70.0541.10.80.024−1.30.70.024−1.2Unaffected \> Tumor \> NormalLeucine \| Lysine\| Lipids\[1.68.. 1.73\]2.10.3201.11.90.024−1.41.90.024−1.4Unaffected \> Tumor \> NormalLeucine \| Lysine\| Lipids\[1.73.. 1.79\]1.40.5861.11.40.024−1.41.40.024−1.3Unaffected \> Tumor \> NormalLipids\[1.79.. 1.84\]1.00.024−1.31.00.024−1.2Unaffected \> Tumor \> NormalLipids\[1.84.. 1.86\]0.70.024−1.30.60.024−1.2Unaffected \> Tumor \> NormalLysine\| Lipids\[1.86.. 1.88\]0.70.024−1.40.70.024−1.2Unaffected \> Tumor \> NormalArginine \| Lysine\[1.88.. 1.94\]1.90.3971.11.70.024−1.41.80.024−1.4Unaffected \> Tumor \> NormalLysine\| Proline\[1.94.. 1.96\]0.80.024−1.40.70.024−1.2Unaffected \> Tumor \> NormalProline\[1.96.. 1.99\]1.20.0931.21.10.024−1.41.00.024−1.2Unaffected \> Tumor \> NormalProline\[1.99.. 2.01\]1.10.0931.11.00.024−1.40.90.024−1.2Unaffected \> Tumor \> NormalPyroglutamate \| Glutamate \| Proline\[2.01.. 2.07\]1.60.0931.11.70.024−1.21.30.548−1.1Unaffected \> Tumor \> NormalGlutamate\[2.07.. 2.09\]0.80.0541.10.90.024−1.30.70.024−1.2Unaffected \> Tumor \> NormalGlutamate \| Glutamine\[2.09.. 2.15\]1.70.9381.0Tumor \> NormalGlutamate \| Glutamine \| Glutathione\[2.15.. 2.19\]0.90.024−1.20.90.024−1.1Unaffected \> Tumor \> NormalUnknown\[2.19.. 2.21\]0.50.024−1.10.50.024−1.1Unaffected \> Tumor \> NormalLipids\[2.21.. 2.24\]0.60.024−1.20.60.024−1.1Unaffected \> Tumor \> NormalValine \| Lipids\[2.24.. 2.30\]1.20.3971.11.20.024−1.31.20.024−1.2Unaffected \> Tumor \> NormalGlutamate\[2.30.. 2.32\]0.70.024−1.20.50.048−1.1Unaffected \> Tumor \> NormalGlutamate \| Proline\[2.32.. 2.37\]2.00.2011.12.00.024−1.21.50.381−1.1Unaffected \> Tumor \> NormalPyroglutamate \| Glutamine\[2.37.. 2.40\]0.80.024−1.30.80.024−1.4Unaffected \> Normal \> TumorGlutamine\| Pyroglutamate\[2.40.. 2.42\]0.40.023−1.20.60.024−1.2Unaffected \> Normal \> TumorGlutamine\[2.42.. 2.48\]1.70.040−1.21.20.0241.2Normal \> Tumor \> UnaffectedPyroglutamate\| Glutathione\[2.48.. 2.51\]0.70.024−1.3Unaffected \> TumorGlutathione\[2.51.. 2.57\]0.90.024−1.31.00.095−1.2Unaffected \> Tumor \> NormalGlutathione\[2.57.. 2.60\]0.60.024−1.30.70.048−1.2Unaffected \> Tumor \> NormalUnknown\[2.60.. 2.62\]0.30.024−1.30.40.095−1.3Unaffected \> Tumor \> NormalMethionine\[2.62.. 2.66\]0.50.0241.1Normal \> UnaffectedAspartate\[2.66.. 2.72\]0.90.017−1.11.00.024−1.31.30.024−1.4Unaffected \> Normal \> TumorDimethylamine\[2.72.. 2.77\]0.50.024−1.10.70.024−1.2Unaffected \> Normal \> TumorAspartate\[2.77.. 2.80\]0.80.054−1.10.80.024−1.41.00.024−1.6Unaffected \> Normal \> TumorAspartate\[2.80.. 2.83\]0.70.071−1.10.70.024−1.30.90.024−1.5Unaffected \> Normal \> TumorAsparagine\[2.83.. 2.87\]0.60.024−1.50.70.095−1.3Unaffected \> Tumor \> NormalAsparagine\[2.87.. 2.90\]0.60.024−1.40.70.024−1.4Unaffected \> Tumor \> NormalTyramine\[2.90.. 2.92\]0.60.024−1.7Unaffected \> NormalAsparagine\[2.92.. 2.95\]0.90.024−1.7Unaffected \> NormalAsparagine\[2.95.. 2.97\]0.80.024−1.40.60.024−1.2Unaffected \> Tumor \> NormalGlutathione \| Lysine\| Lipids\[2.97.. 3.02\]1.00.1211.11.10.024−1.31.00.024−1.2Unaffected \> Tumor \> NormalCreatine \| Creatinine \| Lysine\| Tyrosine\[3.02.. 3.07\]1.70.3971.11.70.024−1.61.70.024−1.4Unaffected \> Tumor \> NormalLipids\[3.07.. 3.09\]0.50.024−1.50.60.095−1.2Unaffected \> Tumor \> NormalPhenylalanine\[3.09.. 3.13\]0.90.024−1.30.90.024−1.2Unaffected \> Tumor \> NormalPhenylalanine\[3.13.. 3.17\]0.70.024−1.30.80.024−1.2Unaffected \> Tumor \> NormalTyrosine\[3.17.. 3.19\]0.60.024−1.40.60.024−1.1Unaffected \> Tumor \> NormalCholine \| O-Phosphocholine \| O-Acetylcholine \| sn-Glycero-3-phosphocholine\[3.19.. 3.23\]2.00.0171.31.80.024−1.30.90.024−1.0Unaffected \> Tumor \> Normalmyo-Inositol \| Glucose \| Taurine\[3.23.. 3.29\]1.40.8151.01.50.0241.21.50.3811.2Tumor \> Normal \> Unaffectedmyo-Inositol\[3.29.. 3.31\]1.20.017−2.11.00.0241.9Normal \> Unaffected \> TumorProline\[3.31.. 3.33\]0.60.024−1.40.80.024−1.6Unaffected \> Normal \> TumorMethanol \| Proline\[3.33.. 3.36\]0.90.024−1.30.90.024−1.2Unaffected \> Tumor \> NormalGlucose\|Glucose-1-phosphate\[3.36.. 3.39\]1.40.815−1.01.50.0245.51.60.0245.3Normal \> Tumor \> UnaffectedGlucose\|Glucose-1-phosphate \| Taurine\[3.39.. 3.44\]2.30.0171.22.80.0242.13.50.0242.5Tumor \> Normal \> UnaffectedGlucose\[3.44.. 3.46\]0.80.0241.8Normal \> UnaffectedGlucose\|Glucose-1-phosphate\[3.46.. 3.51\]3.00.8151.13.30.0244.03.70.0244.3Tumor \> Normal \> UnaffectedGlucose\|myo-Inositol \| Glycine\[3.51.. 3.57\]3.50.040−1.73.10.0242.01.71.0001.2Normal \> Tumor \> Unaffectedmyo-Inositol\| Threonine\[3.57.. 3.62\]2.70.040−1.72.30.0241.81.21.0001.0Normal \> Tumor \> UnaffectedGlycerol \| Unknown\[3.62.. 3.68\]1.40.017−1.31.40.0481.30.50.048−1.0Normal \> Unaffected \> TumorGlucose\[3.68.. 3.70\]0.60.024−1.1Unaffected \> TumorGlucose\[3.70.. 3.72\]1.00.0241.6Normal \> UnaffectedGlucose-1-phosphate \|Glutamate\[3.72.. 3.74\]1.81.0001.02.00.0242.62.20.0242.6Tumor \> Normal \> UnaffectedGlucose-1-phosphate \|Glutamate \|Alanine\[3.74.. 3.80\]3.31.0001.13.40.0241.44.00.0241.5Tumor \> Normal \> UnaffectedGlucose-1-phosphate\| Serine\[3.83.. 3.89\]2.41.0001.02.80.0241.73.20.0241.8Tumor \> Normal \> UnaffectedGlucose-1-phosphate \| Creatine\[3.89.. 3.93\]2.71.0001.12.90.0241.93.30.0242.1Tumor \> Normal \> UnaffectedSerine \| Unknown\[3.93.. 3.97\]1.90.0711.11.70.024−1.51.50.024−1.3Unaffected \> Tumor \> NormalSerine \| Unknown\[3.97.. 4.02\]2.40.5861.11.90.024−1.61.80.167−1.4Unaffected \> Tumor \> Normalmyo-Inositol\| Choline\[4.02.. 4.08\]1.90.040−1.51.60.0241.41.01.000−1.1Normal \> Unaffected \> TumorRibose\| Lactate\[4.08.. 4.11\]0.80.024−1.20.80.024−1.2Unaffected \> Tumor \> NormalLactate\| O-Phosphocholine\| Proline\| Uridine\[4.11.. 4.17\]0.90.0231.11.10.024−1.20.60.095−1.1Unaffected \> Tumor \> NormalO-Phosphocholine \| Pyroglutamate\[4.17.. 4.19\]0.60.024−1.30.60.024−1.4Unaffected \> Normal \> TumorThreonine\[4.22.. 4.28\]0.90.024−1.20.80.048−1.2Unaffected \> Tumor \> NormalLipids\[4.28.. 4.34\]1.20.1561.11.20.024−1.41.10.024−1.3Unaffected \> Tumor \> Normal1-Methylnicotinamide\[4.45.. 4.50\]0.70.024−1.30.60.024−1.2Unaffected \> Tumor \> NormalUnknown\[4.50.. 4.54\]0.50.024−1.50.60.024−1.7Unaffected \> Normal \> TumorGlutathione\[4.54.. 4.59\]0.70.024−1.6Unaffected \> NormalUnknown\[4.59.. 4.63\]0.40.024−1.5Unaffected \> NormalGlucose\[4.63.. 4.68\]0.60.0481.6Normal \> UnaffectedGlucose\[5.20.. 5.26\]0.70.0482.0Normal \> UnaffectedUnknown\[5.32.. 5.37\]0.30.317−2.10.40.095−2.6Unaffected \> Normal \> TumorUnknown\[5.37.. 5.42\]0.50.053−8.90.61.0002.30.40.048−4.0Normal \> Unaffected \> TumorGlucose-1-phosphate\[5.42.. 5.48\]2.70.6971.22.90.02472.03.30.02483.7Tumor \> Normal \> UnaffectedUnknown\[5.48.. 5.53\]0.20.076−6.90.20.095Normal \> Tumor \> UnaffectedUracil\[5.76.. 5.82\]0.60.024−8.70.70.0120.0Unaffected \> Normal \> TumorUridine\[5.88.. 5.94\]1.21.0001.00.90.02422.01.00.02422.4Tumor \> Normal \> UnaffectedTyrosine\[6.87.. 6.93\]0.60.024−1.30.60.024−1.2Unaffected \> Tumor \> NormalUnknown\[7.06.. 7.08\]0.30.048−1.50.30.024−1.6Unaffected \> Normal \> Tumorp-Methylhistidine\[7.08.. 7.14\]0.40.024−1.20.30.095−1.1Unaffected \> Tumor \> NormalTyrosine\[7.16.. 7.22\]0.70.024−1.3Unaffected \> NormalUnknown\[7.22.. 7.27\]0.30.024−1.1Unaffected \> NormalPhenylalanine\[7.27.. 7.29\]0.30.024−1.4Unaffected \> NormalPhenylalanine\[7.29.. 7.34\]0.80.024−1.4Unaffected \> NormalPhenylalanine\[7.34.. 7.39\]0.60.024−1.50.60.095−1.4Unaffected \> Tumor \> NormalUnknown\[7.39.. 7.45\]0.80.024−1.60.70.024−1.4Unaffected \> Tumor \> NormalUracil\[7.45.. 7.51\]0.50.024−3.80.40.155−2.5Unaffected \> Tumor \> NormalUnknown\[7.51.. 7.57\]0.80.024−2.80.90.024−3.7Unaffected \> Normal \> TumorUnknown\[7.71.. 7.76\]0.30.317−3.20.40.048−3.30.50.036−10.6Unaffected \> Normal \> TumorUridine\[7.84.. 7.89\]1.00.697−1.10.70.0243.40.70.0243.2Normal \> Tumor \> Unaffectedp-Methylhistidine\[7.91.. 7.97\]0.30.095−1.40.40.095−1.4Unaffected \> Tumor \> NormalUnknown\[8.33.. 8.38\]0.10.070−2.50.30.155−4.8Unaffected \> Normal \> TumorUnknown\[8.44.. 8.47\]0.10.076−19.50.40.024−20.2Unaffected \> Normal \> Tumor1-Methylnicotinamide\[8.92.. 8.95\]0.20.048−3.40.30.024−6.1Unaffected \> Normal \> Tumor1-Methylnicotinamide\[9.27.. 9.31\]0.30.095−1.5Unaffected \> TumorTable 2Semi-quantitated metabolites that were found to be statistically different based on the median values for untreated tumor, normal and unaffected cells.  The p values were calculated with the Wilcoxon Rank Sum Test (normal vs tumor) and exact Wilcoxon Rank Sum Test (tumor vs unaffected and normal vs unaffected) (p \< 0.1, or FC \> 2).MetaboliteTumor vs NormalNormal vs UnaffectedTumor vs UnaffectedCell Line Differencep-value\*FC\*\*p-value\*FC\*\*p-value\*FC\*\*Aspartate0.048−1.7Unaffected \> TumorDimethylamine0.023−1.30.095−1.4Tumor \< (Normal\| Unaffected)Glucose-1-phosphate0.04−1.50.02455.50.02436Normal \> Tumor \> UnaffectedGlutamine0.017−1.50.0241.4Normal \> (Tumor\| Unaffected)Leucine0.024−2Unaffected \> TumorPhenylalanine0.095−1.7Unaffected \> TumorPyroglutamate0.023−1.40.095−1.30.024−1.9Unaffected \> Normal \> TumorTyrosine0.095−1.4Unaffected \> Tumor

Metabolites identified were incorporated in the pathway enrichment analysis which was performed using the knowledge-based canonical pathways and endogenous metabolic pathways in the MetaCore module of the GeneGo software.  Rankings of the canonical pathways were based on hypergeometric *p*-values and false discovery rate (FDR) being \< 0.05.  The top hits included aminoacyl-tRNA biosynthesis in mitochondria and cytoplasm and signal transduction amino acid-dependent mTORC1 activation (Table [3](#Tab3){ref-type="table"}).Table 3Top ten statistically significant pathway maps created using GeneGo enrichment analysis based on library-matched metabolites the tumor, normal, and unaffected cell lines and are ranked using hypergeometric *p*-values and false discovery rate (FDR).Tumor vs NormalTumor vs UnaffectedNormal vs Unaffectedp-valueFDRp-valueFDRp-valueFDR1Aminoacyl-tRNA biosynthesis in mitochondrion7.62E-191.39E-161.09E-212.04E-194.82E-239.11E-212Aminoacyl-tRNA biosynthesis in cytoplasm1.11E-177.83E-162.47E-201.81E-181.34E-211.01E-193Aminoacyl-tRNA biosynthesis in cytoplasm/ Rodent version1.29E-177.83E-162.90E-201.81E-181.61E-211.01E-194Signal transduction_Amino acid-dependent mTORC1 activation5.18E-162.36E-142.24E-151.05E-135.00E-152.36E-135Nociception_Pro-nociceptive action of Nociceptin in spinal cord at low doses1.66E-065.04E-052.78E-061.04E-043.68E-061.39E-046Glycine, serine, cysteine and threonine metabolism1.74E-065.04E-053.51E-061.05E-045.16E-061.55E-047Glycine, serine, cysteine and threonine metabolism/ Rodent version1.94E-065.04E-053.91E-061.05E-045.75E-061.55E-048Nicotine signaling in cholinergic neurons4.03E-069.16E-056.69E-061.56E-048.85E-062.09E-049Galactose metabolism1.00E-052.02E-041.66E-053.44E-042.18E-054.59E-0410Leucine, isoleucine and valine metabolism5.37E-059.78E-048.02E-051.50E-031.00E-041.89E-03

To further investigate the differences between the two genotypes, a supervised analysis exhibited a well-defined separation (Supplemental Fig.  [1](#MOESM1){ref-type="media"}), and Supplemental Tables [1](#MOESM1){ref-type="media"} and [2](#MOESM1){ref-type="media"} list the metabolites identified and semi-quantified that were significantly different.  The metabolites and pathways that differentiated the normal and tumor cells lines were different from the two tumor cell lines.  Two metabolites were noted to have high fold change differences between the genotypes in the binned data.  The bins associated with 1-methylnicotinamide was between 2.1--4.5 fold higher in S45F; additionally, the bin identified as p-methylhistidine was 4.5 fold higher in the S45F cell line.  1-methylnicotinamide is involved in NAD metabolism, while p-methylhistidine is involved in beta-alanine and histadine metabolism.  The semi-quantified metabolites were all greater in the T219 cell line. 

Treatment of cells with Dasatinib and FAK Inhibitor 14 alters metabotype of all cell lines {#Sec4}
------------------------------------------------------------------------------------------

The differences between the response to Dasatinib and FAK Inhibitor 14 by the desmoid tumor, normal, and unaffected cell lines were investigated with the unsupervised multivariate analysis (PCA), as show in Fig.  [2A--B](#Fig2){ref-type="fig"}.  The differentiation of the metabolomics profiles were affected and unaffected on the first component.  The Dasatinib treated cells separated by tumor and normal on the 2^nd^ component (A); the FAK Inhibitor 14 cells did not (B).  The supervised analysis of the Dasatinib and FAK Inhibitor 14 treated cells showed excellent differentiation between the treated and vehicle (Supplemental Fig.  [2A--F](#MOESM1){ref-type="media"}), and Supplemental Table [3](#MOESM1){ref-type="media"} and [4](#MOESM1){ref-type="media"} list the metabolites identified and semi-quantified metabolites that distinguish the two treatments.  In the Dasatinib treatment, it was observed that the lipid bins in the normal cell lines were lower in the vehicle, while there was no significant difference in the tumor cells lines.  In contrast, there was a lipid increase in the tumor cells and a decrease in the normal cells when treated with the FAK Inhibitor 14.  The bins identified as 1-methylnicotinamide were observed to be significantly different in the FAK Inhibitor 14 treatment in the tumor cell lines, while there was no difference in the Dasatinib treatment. Figure 2The unsupervised multivariate analysis (PCA) of desmoid, normal, and unaffected cell lines in the presence of Dasatinib (**A**) and FAK Inhibitor 14 (**B**).  The PCA in the presence of Dasatinib \[6 components, R2X = 0.99, Q2 (cum) = 0.92\].  The PCA in the presence of FAK Inhibitor 14 \[7 components, R2X = 0.99, Q2 (cum) = 0.94\]. 

As was seen in the untreated comparison, the two genotypes differentiated in the treatment responses, and the supplemental Table [5](#MOESM1){ref-type="media"} highlight the semi-quantitative therapeutic responses.  Generally, the treatments caused the metabolites to increase in concentration in the T219 (S45F) cell line and decrease in the T141 (T41A) cell line.  It was observed, however, that the Dasatinib treatment resulted in a 25.5-fold increase in glucose-1-phosphate in the T219 cell line, but no significant change in either treatment in the T141 cell line.  There were similarities in the response for the normal and the T141 cells; both had decreases in glutathione and pyroglutamate, while there was a slight increase or no change in the T219 cells for these two, respectively.  While there was a slight increase in uridine in the T219, Dasatinib treatment, there was an 8.8-fold decrease when treated with the FAK Inhibitor 14. 

Discussion {#Sec5}
==========

To better understand underlying desmoid tumor cell behavior, this pilot broad spectrum metabolomics project was used to investigate metabolomic differences in matched normal and desmoid tumor cells with unaffected fibroblast cells.  These cells are derived from affected patients and may have been the only matched cell lines that exist to date for this rare tumor type.  While a small study, this provides a unique opportunity to understand the differences between affected and unaffected fibroblast cells.  The cells clearly separated by unsupervised and supervised multivariate analysis and semi-quantitative concentration by affected and unaffected, in addition to tumor, normal, and unaffected.  Thus, the metabolomic profile of desmoid tumors and normal fibroblasts from the same patient differ.  Most studies have evaluated affected and unaffected tissue; it is not unexpected that the normal and unaffected fibroblast cells differentiate^[@CR22]^. Huang, Stern, and Zhao showed that the studying tumors micro-environments may be more informative for markers of genetic susceptibility and disease monitoring in breast cancer, kidney clear cell carcinoma, lung adenocarcinoma, liver cancer, lung squamous cell carcinoma, and head and neck cancer.  These paired samples are rare; even in breast cancer research, only 105 paired samples have been identified and studied.  The micro-environment is particularly relevant for desmoid tumors, which lack the ability to metastasize^[@CR7]^.

The metabolites that were found to be significantly different in the normal, desmoid tumor, and unaffected (Tables [1](#Tab1){ref-type="table"} and [2](#Tab2){ref-type="table"}) were represented in the aminoacyl-tRNA biosynthesis in mitochondria and cytoplasm (Table [3](#Tab3){ref-type="table"}).  In most cases, the amino acids and lipids were greatest in the unaffected cell lines, and the desmoid tumor cell lines had greater concentrations than normal.  Differences in the amino acid concentrations may be an indicator that there is a dysfunction in aminoacyl-tRNA synthetases.  The aminoacyl-tRNA synthetases (ARSs) are responsible for protein translation, catalyse the ligation of amino acids to their cognate tRNAs, and interact with various proteins with tumorigenesis implications^[@CR23],[@CR24]^. Aspartate and glycine were two amino acids that were greater in the normal cells compared to the desmoid tumor.  The aspartyl-tRNA synthetase (AspRD) is known to interact with the multisynthetase complex p43 (MSC^p43^), which in turns antiproliferative and proapoptotic signaling pathways related to FUSEbinding protein (FBP), Myc, TNF receptorassociated factor 2 (TRAF2) and p53^[@CR25]^. Glycyl-tRNA synthetase (GRS) has been previously found to be overexpressed in breast cancer cell lines^[@CR26]^, in addition to GRS being secreted to the serum in defense of ERK-activated tumorigenesis^[@CR27]^. Another amino acid that differentiated the cell types was tyrosine.  In addition to the normal ARS functions, Tyrosyl-tRNA synthetase (TyrRS) contains the tripeptide ELR (GluLeuArg) and the EMAPII domain^[@CR28]^. The ELR peptide has been identified in lung cancer as a chemokine^[@CR29]^. Additionally, the EMAPII domain inhibits integrin-dependent cell adhesion with the interaction of α5β1 integrin receptor on the surface of endothelial cells^[@CR28]^. As previously noted, Wnt/beta-Catenin signaling also play a critical role in cell adhesion, and Yum *et al*.  recently investigated the relationship between AIMP2, such as TyrRS, and Wnt/beta-catenin signaling in a murine model of intestinal homeostasis and tumorigenesis^[@CR29]^. The AIMP2 was thusly found to regulate the Wnt/beta-catenin signaling in the intestinal tract; however, it is not clear from their results if or how the signaling cascade would be altered with the known desmoid tumor mutation of beta-catenin, T41A or S45F. ARS have been identified as drug targets, such as aspartyl-tRNA synthetase for Mycobacterium tuberculosis^[@CR30]^, or glycyl-tRNA synthetase for Charcot-Marie-Tooth peripheral neuropathy^[@CR31]^. Because it is possible to target one ARS verses the others using known binding motifs^[@CR32]^, future screening experiments would target several ARS inhibitors against these cell lines to find most efficacious results. 

Another differentiated pathway was signal transduction amino acid-dependent mTORC1 activation pathways.  The mTOR complex 1(mTORC1) signaling is reduced dramatically in the absence of arginine and the branched chain amino acids^[@CR33]^. The activity is recovered when the amino acids, particularly leucine, are added back into the reaction^[@CR34]^. Leucyl-tRNA synthetase (LARS) was identified as sensor in the mTORC1 signaling, also independent of its translational function^[@CR35]^. Furthermore, there was a clear difference in the glucose and glucose-1-phosphate concentrations in the normal, desmoid tumor, and unaffected cell lines.  In an elegant set of experiments, Pusapati *et al*.  found that the activation of mTORC1 signaling in solid tumors can avoid glycolysis block and ensure adequate energy for cellular survival^[@CR36]^. The metabolomics results and these ARS and mTORC1 results indicate that further investigations are needed to further understand this pathways in desmoid tumorigenesis.  mTORC1 has been shown to play a role in cell proliferation and tumor growth in transgenic animals carrying APC mutations similar to what is seen in desmoid tumor patients^[@CR37]^. Sirolimus, also known as rapamycin, is an immunosuppressive agent and known mTOR inhibitor, and is currently being investigated in a pediatric desmoid tumor clinical trial (NCT01265030).  The primary outcome of this trial is to determine if mTOR inhibition is beneficial for children and young adults with desmoid tumors.  Future metabolomic studies will also investigate the response of the sirolimus in the cell lines. 

The metabolomics profiles for the two genotypes (T41A and S45F) were also distinct and revealed 1-methylnicotinamide and p-methylhistidine as being significantly greater in the S45F cell line.  1-methylnicotinamide is involved in NAD metabolism.  Like mTORC1 signaling, NAD metabolism has also been suggested as a vital metabolic rewiring for energy metabolism^[@CR38]^ P-methylhistidine is involved in histidine metabolism, which has been noted previously as markers of patients with cachexia and pancreatic cancer^[@CR39]^. p-methylhistadine is also known to bind histidine decarboxylase, where Masini *et al*.  found that histamine concentration and histidine decarboxylase activity is correlated with colorectal cancer tumor stages and intratumor microvessel density^[@CR40]^. Surgically resected colorectal adenocarcinomas showed an increase in both histamine concentration and histidine decarboxylase with increased staging.  FAP manifests as intra-abdominal desmoid tumors in 25% of FAP patients^[@CR41]^. The genetic cause are APC mutations, rather than beta-catenin, so these results may indicate a similar dysregulation of histadine metabolism.  The NAD(+) salvage pathway has been identified as important for malignancies, such multiple myeloma^[@CR42]^ and breast cancer^[@CR43]^. In three breast cancer cell lines and 176 primary tumor cell lines, Sharif *et al*.  found these cell lines had high levels of the rate-limiting enzyme of the NAD(+) salvage pathways, nicotinamide phosphoribosyltransferase (NAMPT) and lower levels of p73, both correlating to poor survival rates.  While p53 did not have the same correlation with NAMPT or survival, the p53 pathway was recently identified as regulating Wnt/beta-catenin signaling in mesenchymal progenitor cells^[@CR44]^. There has been conflicting results in the literature evaluating the role of p53 in desmoid tumor patients.  For example, there is data that shows the overexpression of p53 and Ki-67 indicates a high probability of recurrence^[@CR45]^, while high expression of both beta-catenin and p53 are markers of the risk reduction^[@CR46]^. Future studies need to evaluate the relationship between NAMPT, p53, and beta-catenin mutations in these mesenchymal tumors. 

In addition to the untreated cells, the broad spectrum metabolomics identified differences in the metabolic response to two drugs identified by the ongoing screening efforts.  Both Dasatinib and FAK Inhibitor 14 are kinase inhibitors.  Dasatinib is l Bcr-Abl and Src family tyrosine kinase inhibitor^[@CR14]^, while FAK inhibitors target the intersecting integrin and receptor tyrosine kinase signal transduction pathways^[@CR15]^. Both treatments target proteins that directly interact with beta-catenin upstream of the Wnt pathway.  Src and FAK form a dual kinase complex, promote cell motility, cell cycle progression and cell survival^[@CR47]^. One of their many downstream molecules is Akt, which phosphorylates GSK3B, regulating beta-catenin phosphorylation.  It is likely that both inhibitors are affecting these same links, although both Src (and other molecules that Dasatinib potentially inhibits) and FAK have their other independent functions.  While FAK activity in desmoid tumors has not yet been evaluated, others have considered Src activity and found it to be variable^[@CR48]^.

The metabolites that were found to be perturbed in the treated cells were the same as those identified as differentiating the untreated cells.  The response to treatment, however, was markedly different by cell and genotype, as expected.  The lipid bins were lower in the FAK Inhibitor 14 treated tumor cells, while there were no changes associated with the Dasatinib treatments compared to the vehicle.  This was not surprising since tyrosine kinase activity is regulated by lipid membranes and contain a lipid-binding domains^[@CR49]^. Lipid dysregulation is widely studied in cancer and cancer metabolism as a means of signaling^[@CR50],[@CR51]^; therefore, FAK Inhibitor 14 could constrict the tumorigenesis signaling pathways through the lipids^[@CR52]^. Additionally, the 1-methylnicotinamide bins were 1.6--3.7-fold higher in the FAK Inhibitor 14 cells compared to vehicle, while no difference was observed in the normal cell lines. 

Glucose and lactate were decreased in the FAK Inhibitor 14 treated cells binned data, suggesting an altered energetic state through NAD metabolism.  Differences in glucose-1-phosphate and uridine in the T219 treated cells are related to the alterations in energy metabolism.  Uridine has been identified as a biomarker in two colorectal cancer studies^[@CR53],[@CR54]^. Nucleotide production is dependent on glycolysis and the TCA cycle metabolic intermediates^[@CR55]^. Targeting glucose and nucleotide metabolism in cancer research has proved to be difficult but important^[@CR56],[@CR57]^, particularly in the tumor micro-environment^[@CR58]^.

The two treatments also caused in a difference in the branched chain amino acids (BCAAs).  Isoleucine, leucine, and valine have been noted in the literature as markers of cancer as integral players in the mTORC1 pathway^[@CR59]^. As shown in Supplemental Table [5](#MOESM1){ref-type="media"}, the BCAAs all increased, except for the T141 treated with FAK 14 Inhibitor, which decreased 2.3-fold compared to vehicle.  Since concentrations of these three amino acids usually trend together. 

Asparagine was another amino acid that had an equal but opposite change in response to the treatments.  There was a 2.7-fold increase in asparagine for T219 Dasatinib treated cells while decreased in T141 FAK Inhibitor 14 by 2.4-fold.  mTORC1 activity and protein synthesis is regulated by asparagine, and helps regulate uptake of serine, arginine and histidine^[@CR60]^. While these amino acids were not observed as differentiating the cells, the mTORC1 pathway was one of the top pathways that matched the data. 

Metabolomics has been used to better understand the impact of another Bcr-Abl tyrosine kinase inhibitor, Gleevec (imatinib), on a gastrointestinal stromal tumor (GIST) T1 cell line^[@CR61]^. Gleevec has been found to be effective in some desmoid tumor studies^[@CR62],[@CR63]^; therefore, it was expected to identify similar changes in the desmoid study.  The most important metabolites for differentiating the treated and untreated GIST cells were decreases in glutathione, glutamate, aspartate, myo-inositol, glycerophosphocholine, and phosphocholine.  A decrease in the glutathione, glycerophosphocholine, and phosphocholine were attributed to inhibited growth and invasiveness, while a decrease in the glutamate and glutamine could be indication of a change in glycolysis moving from cytosol to mitochondria.  Several differences were observed comparing these metabolites and tyrosine kinase inhibitors.  In the desmoid tumor study, aspartate increased in the Dasatinib T219 cells 2.7-fold.  No significant myo-inositol differences were observed when the cells were treated in this study.  Glutathione metabolism alterations have been previously identified in cervical, colorectal, gastric and esophageal, bone marrow, breast, colon, larynx, and lung cancer^[@CR64],[@CR65]^. FAK 14 Inhibitor caused glutathione to decreased 2-fold in T141 cells, while Dasatinib increased glutathione 1.5-fold in the T219.  Glutathione metabolism has been used as a marker of therapeutic efficacy for some time, as it is a marker of proliferative and cell cycle progression^[@CR66]^. For example, glutathione S-transferase polymorphisms have been evaluated as a chemotherapeutic marker toxicity and reduced survival of childhood acute myeloid leukemia^[@CR67]^. Future desmoid tumor studies could investigate targeting glutathione S-transferase in desmoid tumor samples to better understand the these differences in therapeutic responses. 

While this was a small pilot study, this was the first attempt to utilize the broad spectrum metabolomics to explore matched desmoid tumor and normal fibroblast cell lines in comparison to an unaffected cell line.  Differences in the cell lines were found to be related to aminoacyl-tRNA biosynthesis in mitochondria and cytoplasm and signal transduction amino acid-dependent mTORC1 activation pathways.  Future studies will investigate if the aminoacyl-tRNA synthetases differentiate these cell lines in follow-up studies, to include aspartyl-tRNA synthetase, glycyl-tRNA synthetase, tyrosyl-tRNA synthetase, and leucyl-tRNA synthetase.  Once the pathways and proteins have been verified, there will be a concerted effort to determine the mechanism.  Metabolite profile differences were observed in the two beta-catenin genotypes, T41A and S45F. The novel treatments identified by the Desmoid Collaboration for a Cure, Dasatinib and FAK Inhibitor 14, altered the normal and desmoid tumor genotype profiles differentially, signifying a potential opportunity for directed therapies in the future.  These results will need to be validated in a much larger data set but provide a tangible step in this rare disease understanding. 

Methods {#Sec6}
=======

Cell lines and culturing {#Sec7}
------------------------

Primary desmoid and normal fibroblast cell cultures were established as previously described^[@CR68],[@CR69]^. Tumor cell line 141 (T141) was genotyped to have the T41A beta-catenin mutation, while tumor cell line 219 (T219) was S45F beta-catenin mutation.  Both of the desmoid cell lines had matched normal fibroblast cell lines, which are denoted N141 and N219.  For normal fibroblasts, the cells are derived similarly but from skin tissue from the same desmoid tumor patients.  Obtaining paired samples is extremely difficult, and it is possible that at the time of this experiment that these matched cell lines are the only ones of their kind.  The fibroblast cell line (CRL7518) was purchased from ATCC and culture identically to the primary cell lines; this ATCC cell line served as the unaffected cell line.  Monolayer cell cultures were grown in DMEM supplemented with 10% fetal bovine serum and maintained at 37 °C in 5% CO~2~. Cells were divided when confluent and experiments were performed between the third and sixth passages. 

*In vitro* inhibitors treatment and sample preparation {#Sec8}
------------------------------------------------------

Approximately 10 × 10^6^ cells were treated with 1.0 uM Dasatinib (Selleck, Houston, USA) dissolved in DMSO, or 0.5 µM FAK Inhibitor 14 (Cayman Chemicals Company, Michigan, USA) dissolved in water.  Cells were incubated in fresh media containing the inhibitors, or vehicle (DMSO for Dasatinib; H~2~O for FAK Inhibitor 14), for 24 hours.  Both Dasatinib and FAK Inhibitor 14 were found to inhibit desmoid cellular growth but not normal fibroblast viability in our on-going screening studies, and there were no visible signs of cell death after the 24 hours.  An aliquot of the media following treatment was collected, and the remainder of the media was aspirated.  The cells were washed with PBS twice, and quenched with 8 mL ice-cold isotonic 0.9% (w/v) saline for 2 minutes.  Total cellular content was then extracted with 1.7 mL ice-cold acetonitrile/water (50:50, v/v) solution.  Cell extracts were collected using a cell scraper and quickly transferred to MagNA Lyser Green Beads tubes (Roche, Indianapolis, USA) and stored in −80 °C.  Media was added to empty plates and incubated together with the cells for the duration of the experiment served as a blank. 

Cells were homogenized on the MagNA Lyzer, with two 30-sec cycles at 2000 rpm, resting in a −20 °C chilling block for 1 min in between pulses, and centrifuged samples at 16,000 rcf for 4 min.  The cell lysate was transferred to a new 2 mL Lo-Bind Eppendorf tubes, with the final cell count approximately 10 × 10^6^ cells for each sample.  Of the twenty cell lysate samples, six samples had sufficient volume for study samples and to be included in an analytical quality control (QC) total pool.  Aliquots from these cell lysate samples were combined, divided into three total pool aliquots, and processed identically to the cell lysate study samples.  All study and pool samples were lyophilized to dryness and reconstituted in a 0.2 M phosphate buffer, pH 7.4, in D~2~O with 10% Chenomx ISTD. 

NMR data acquisition and analysis {#Sec9}
---------------------------------

Data acquisition, statistics, and pathway analysis were performed as previously described^[@CR17],[@CR70],[@CR71]^. Three NMR spectra were acquired for each of the individual study samples and the pooled samples.  ^1^H NMR spectra of cell lysate samples were acquired on a Bruker Avance III 700 MHz NMR spectrometer (located at the David H. Murdock Research Institute at Kannapolis, NC, USA) using a NOESY1D (noesypr1d) pulse sequence.  NMR spectra were pre-processed using ACD 1D NMR Processor 12.0 (ACD Labs, Toronto, Canada).  NMR bins (0.50--9.30 ppm) were made after excluding water (4.70--5.20 ppm) and regions with low signal to noise^[@CR72]^ (5.95--6.85, 8.47--8.85, 9.00--9.25 ppm) using intelligent binning width of 0.04 ppm and 50% looseness factor.  Integrals of each of the bins were normalized to total integral of each of the spectrum. 

Descriptive statistics and two-sided t-tests, using the Satterthwaite approximation for unequal variances, were conducted for the tumor and normal binned NMR data (SAS Institute Inc, Cary, NC).  When there were at least 6 samples in each group of a binary comparison, the Wilcoxon rank sum test was used; for sample sizes smaller than this, the exact Wilcoxon rank was employed.  Spectral replicates were treated as independent samples for this pilot study, and p-values \< 0.1 were considered to be statistically significant and were not adjusted for multiple testing^[@CR73],[@CR74]^. Normalized binned NMR data were mean centered and Pareto scaled prior to multivariate analysis.  Multivariate data analysis methods (e.g. principal component analysis \[PCA\], orthogonal partial least squares discriminant analysis \[OPLS-DA\]) were used to reduce the dimensionality and to enable the visualization of the separation of the study groups (SIMCA 14.1, Umetrics, Umeå, Sweden).  The PCA plots were inspected to ensure that the pooled samples were tightly clustered in the center of all of the individual study samples, a quality control method that is widely used in metabolites studies^[@CR75]^. All models used a 7-fold cross-validation to assess the predictive ability of the model (Q^2^).  Loadings plots and variable influence on projections (VIP) plots were inspected, and bins that had a VIP ≥ 1.0 with a jack-knife confidence interval that did not include 0 were determined to be important to differentiating the study groups.  Chenomx NMR Suite 8.2 Professional software (Edmonton, Alberta, Canada), which has a concentration library of approximately 350 compounds, was used to match the signals in the identified bins to metabolites.  Chenomx was also used to semi-quantify metabolites, and all concentrations were adjusted to the cell count for each sample. 

Metabolites identified as important (VIP ≥ 1.0, p \< 0.1, or magnitude of fold change (FC) \> 2) were analyzed for pathway enrichment analysis using the knowledge-based canonical pathways and endogenous metabolic pathways in the MetaCore module in GeneGo software (Chicago, IL).  Ranking of relevant pathways was based on hypergeometric p-values. 

The metabolomics data are available for download at the NIH Common Fund Metabolomics Data Repository and Coordinating Center at the University of California at San Diego (Dr. Shankar Subramaniam, PI, *U01-DK097430*) under study ST000454. 

Electronic supplementary material
=================================

 {#Sec10}

Supplementary Information

**Electronic supplementary material**

**Supplementary information** accompanies this paper at 10.1038/s41598-017-18921-7. 

**Publisher\'s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. 

This work was supported by a New Investigator (Mercier) and Desmoid Collaboration for a Cure (Alman) grants from the Desmoid Tumor Research Foundation.  The authors thank the NIH Eastern Regional Comprehensive Metabolomics Research Core at RTI (RTI RCMRC) for their support.  Appreciation also is given to Dr. Richele Thompson for the careful reading of this manuscript. 

K.M. and B.A. designed this pilot study.  K.M., M.A.J., and R.P. collected the data.  K.M. and Z.A. interpreted the data.  K.M. and M.A.J. drafted the manuscript, and B.A. provided the critical revisions.  All authors approved of the final draft for publication. 

Competing Interests {#FPar1}
===================

The authors declare that they have no competing interests.