Distinct Metabolomic Profiles of Papillary Thyroid Carcinoma and

Jul 1, 2015 - Papillary thyroid carcinoma (PTC) and benign thyroid adenoma (BTA) are the most common head and neck tumors. However, the metabolic ...
1 downloads 0 Views 2MB Size
Article pubs.acs.org/jpr

Distinct Metabolomic Profiles of Papillary Thyroid Carcinoma and Benign Thyroid Adenoma Yanan Xu,†,‡,∥ Xiaojiao Zheng,§ Yunping Qiu,§ Wei Jia,§ Jiadong Wang,*,‡,∥,⊥ and Shankai Yin*,†,∥,⊥ †

Department of Otolaryngology Head and Neck Surgery, Shanghai Jiao Tong University Affiliated Sixth People’s Hospital, Yishan Road 600, 200233 Shanghai, China ‡ Department of Head and Neck Surgery, Renji Hospital of Shanghai Jiao Tong University School of Medicine, Shandongzhong Road 145,200001 Shanghai, China § Center for Translational Medicine, Shanghai Jiao Tong University Affiliated Sixth People’s Hospital, Yishan Road 600, 200233 Shanghai, China ∥ Otolaryngological Institute of Shanghai Jiao Tong University, Yishan Road 600, 200233 Shanghai, China S Supporting Information *

ABSTRACT: Papillary thyroid carcinoma (PTC) and benign thyroid adenoma (BTA) are the most common head and neck tumors. However, the metabolic differences between PTC and BTA have not been characterized. The aim of this study was to identify the metabolic profiles of these two types of tumors using a metabolomics approach. Tumors and adjacent nontumor specimens collected from 57 patients with PTC and 48 patients with BTA were profiled using gas chromatography−time-of-flight mass spectrometry and ultraperformance liquid chromatography-quadrupole time-of-flight mass spectrometry. A panel of 46 and 44 differentially expressed metabolites were identified in the PTC and BTA specimens, respetively, and compared with nontumor tissues. Common metabolic signatures, as characterized by increased glycolysis, amino acid metabolism, one carbon metabolism and tryptophan metabolism, were found in both types of tumors. Purine and pyrimidine metabolism was significantly elevated in the PTC specimens, and taurine and hypotaurine levels were also higher in the PTC tissues. Increased fatty acid and bile acid levels were found, especially in the BTA tissues. The metabolic profiles of the PTC and BTA tissues include both similar and remarkably different metabolites, suggesting the presence of common and unique mechanistic pathways in these types of tumors during tumorigenesis. KEYWORDS: metabolomics, papillary thyroid carcinoma, benign thyroid adenoma, gas chromatography−time-of-flight mass spectrometry, ultraperformance liquid chromatography-quadrupole time-of-flight mass spectrometry



INTRODUCTION

malignant tumor in the population, as evidenced by cadaveric studies from Finland, where one-third of patients who died from non-thyroid-related conditions were found to have thyroid cancer.4 Thyroid cancer is currently the fifth most common

The incidence of thyroid nodules has dramatically increased in the past decade.1 Depending on the populations evaluated and the methods of detection used, the prevalence of thyroid nodules varies from 5% by palpation to 30−67% by ultrasound evaluation.2 Although most of these thyroid nodules are benign, 5−20% are malignant.3 As a result, thyroid cancer could be a common © 2015 American Chemical Society

Received: April 24, 2015 Published: July 1, 2015 3315

DOI: 10.1021/acs.jproteome.5b00351 J. Proteome Res. 2015, 14, 3315−3321

Article

Journal of Proteome Research cancer in women, accounting for 51.3% of head and neck malignancies and 3.8% of all cancers in the United States.1 Papillary thyroid carcinoma (PTC) accounts for ∼80% of all types of thyroid cancer.5 Laryngeal squamous cell carcinonma (LSCC) was the second most frequently occurring head and neck cancer.6 According to GLOBOCAN, there was an estimate of 156 877 new cases and 83 376 deaths of laryngeal cancer in the world in 2012, and the adjusted incidence and mortality rates were 2.1/100 000 and 1.1/100 000, respectively. To date, very few metabolomics studies have been conducted on head and neck cancers.7−10 Deja found that metabolic changes in thyroid cancer were related mainly to osmotic regulators (taurine and scyllo- and myo-inositol), citrate, and amino acids supplying the tricarboxylic acid (TCA) cycle.7 Yao found that the major metabolic differences between benign and malignant thyroid nodules in serum are related to lipid metabolism.8 Tripathi9 and Somashekar10 reported that LSCC tissue exhibited highly active glycolysis with increased amino acids influx (anaplerosis) into the TCA cycle, altered membrane choline phospholipid metabolism, and oxidative and osmotic defense mechanisms. These findings suggest that benign and malignant thyroid tumors have different metabolic phenotypes, and different types of head and neck cancers may also exhibit distinct metabolite profiles, resulting from varied upstream genetic regulations. In this study, metabolomic profiling was applied to PTC, BTA, and LSCC tissues using a combined gas chromatography−timeof-flight mass spectrometry (GC−TOFMS) and ultraperformance liquid chromatography−quadrupole time-of-flight mass spectrometry (UPLC−QTOFMS), which have not been reported in the field of head and neck cancer research. The use of the two mass spectrometry platforms offers high-level sensitivity and selectivity, providing remarkable metabolic information regarding the differences between PTC and BTA. LSCC specimens were used as a positive control, to ensure that the metabolic differences identified between PTC and BTA were specific to the thyroid condition.



Table 1. Demographic and Clinical Chemistry Characteristics of the Human Subjects number of patients age (mean) male/female ratio T stage T1 T2 T3 T4 N stage N0 N1 N2 N3 a

papillary thyroid cancer

benign thyroid adenoma

57 44.5 ± 14.4 14/43

48 48.0 ± 13.2 5/43

46 9 2 0

a

33 24

Without classification.

added, and the sample was homogenized for 4 min, followed by centrifugation at 12 000 rpm for 10 min. A 150 μL aliquot of the supernatant was transferred into a new tube. The residue was extracted with 250 μL of methanol using a homogenizer, followed by centrifugation. A 150 μL aliquot of the supernatant was combined with the previous extraction sample. Pooled quality-control samples were prepared by mixing 20 μL of each sample extraction. After the addition of the internal standard (20 μL of L-2-chlorophenylalanine in water, 0.03 mg/mL), the mixture was vortexed and aliquoted into a 1.5 mL tube (100 μL) and a GC vial (100 μL) for UPLC−QTOFMS and GC−TOFMS analysis, respectively. These samples were then vacuum-dried at room temperature. For UPLC−QTOFMS analysis, 350 μL of a mixture of methanol and acetonitrile (1:9 v/v) was added to the residue, followed by vortexing for 2 min and ultrasonication for 1 min. After centrifugation, the supernatant was transferred to an UPLC−QTOFMS vial. For GC−TOFMS analysis, the residue was chemically derivatized with a two-step procedure following the protocols previously published.12 In brief, the samples were subjected to a derivatization procedure using 80 μL of methoxyamine (15 mg/mL in pyridine) for 90 min at 30 °C, followed by 80 μL of N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA) (1% trimethylchlorosilane (TMCS)) for 60 min at 70 °C.

MATERIALS AND METHODS

Sample Information

Surgical specimens (n = 276) were collected from both tumor and adjacent nontumor tissues from 57 patients with PTC and 48 patients with BTA treated at Shanghai Renji Hospital and 33 patients with LSCC treated at Shanghai Sixth People’s Hospital. Tumor tissues from the center of the lesions and corresponding adjacent nontumor tissues in the same patient were obtained. All tumor tissue samples were carefully dissected to exclude surrounding normal tissue. All samples were collected within 15 min of obtaining the surgical specimens and immediately frozen at −80 °C until metabolomic analysis. The age, gender, and tumor stage of all patients are provided in Table 1. None of the patients were on neoadjuvant chemotherapy before surgical treatment. The protocol was approved by Shanghai Sixth People’s Hospital and Shanghai Renji Hospital, and all participants in this study signed informed consent before the study.

Instrumental Analysis and Metabolite Identification

Metabolomic profiling via UPLC-QTOFMS and GC-TOFMS was performed following our previous publication with minor modifications.12 The samples were run in the order of PTC − PTC adjacent tissue − BTA − BTA adjacent tissue − LSCC − LSCC adjacent tissue. One quality-control sample and one blank vial were run after each of the 10 samples. A Waters ACQUITY ultraperformance liquid chromatography (UPLC) system equipped with a binary solvent delivery manager and a sample manager (Waters Corporation, Milford, MA), coupled to a tandem quadrupole-time-of-flight (Q-TOF) mass spectrometer equipped with an electrospray interface (Waters Corporation) was used as UPLC-QTOFMS platform. An Agilent 6890N gas chromatographer coupled to a Pegasus HT time-of-flight (TOF) mass spectrometry system (Leco Corporation, St. Joseph, MI) was used as the GC-TOFMS platform. The MS data were analyzed to identify potential discriminate variables.12 For UPLC−QTOFMS, the raw data were analyzed

Sample Preparation

The metabolite extraction procedure followed that from our previous publication with minor modifications.11 Each tissue sample weighed ∼50 mg and was extracted via two steps. A 250 μL mixture of chloroform, methanol, and water (2:5:2 v/v/v) was 3316

DOI: 10.1021/acs.jproteome.5b00351 J. Proteome Res. 2015, 14, 3315−3321

Article

Journal of Proteome Research

Figure 1. OPLS-DA scores plot for tumor tissues relative to adjacent nontumor tissues. (A) OPLS-DA scores plot for papillary thyroid cancer (PTC) samples. (B) OPLS-DA scores plot for benign thyroid adenoma (BTA) samples.



using the MarkerLynx applications manager version 4.1 (Waters, Manchester, U.K.). The resulting matrix contained the arbitrarily assigned peak index (retention time (Rt)-m/z pairs), sample names, and peak intensity. After data processing, a list of the intensities of the peaks detected was generated for the first sample, using Rt and m/z data pairs as the identifiers for each peak. The acquired MS files from GC−TOFMS analysis were exported in the NetCDF format using ChromaTOF software (v3.30; Leco, Los Angeles, CA). CDF files were extracted using custom scripts (revised Matlab toolbox HDA) in the MATLAB 7.1 (The MathWorks, Natick, MA) for data pretreatment procedures such as baseline correction, denoising, smoothing, and alignment, time-window splitting, and peak feature extraction (based on multivariate curve resolution algorithm). The resulting 3D data set included sample information, peak retention time, and peak intensities. Compound annotation was performed using our in-house library containing ∼800 mammalian metabolite standards. For UPLC−QTOFMS generated data, identification was performed by comparison with the accurate mass and RT of the reference standards in our in-house library and with the accurate mass of the compounds obtained from web-based resources such as the Human Metabolome Database (http:// www.hmdb.ca/) and METLIN (http://metlin.scripps.edu). For GC−TOFMS generated data, identification was processed by comparison with the mass fragments and RT of the reference standards in our in-house library or with mass fragments from NIST 05 Standard mass spectral databases in NIST MS search 2.0 software (NIST, Gaithersburg, MD) with a similarity of >70%.

RESULTS

In total, 195 metabolites were detected among the study specimens. In total, 69 were detected by GC−TOFMS, 51 were detected by UPLC-QTOFMS in positive ion mode, and 75 were detected by UPLC−QTOFMS in negative ion mode. A plot of the PCA scores from the entire metabolite data set did not show separations, but instead a trend of separation between the PTC and BTA tissues compared with their adjacent nontumor tissues. The application of a supervised multivariate statistical model, the OPLS−DA model, demonstrated clear separations between each tumor and the adjacent nontumor tissue pair (Figure 1, parameters for the OPLS−DA model). Supplementary Tables 1 and 2 in the SI show two panels of 54 and 49 differentially expressed metabolites in PTC and BTA tissues compared with adjacent nontumor tissues (VIP > 1 and P < 0.05, respectively). To evaluate the impact of gender imbalance on the identified discriminant metabolites, we used a smaller, age-matched sample set (14 males vs 14 females) of adjacent nontumor tissues from PTC patients and constructed a new OPLS-DA model. Twentyone metabolites with VIP > 1 and P < 0.05 (Supplementary Table 3 in the SI) were differentially expressed in male PTC nontumor tissues compared with the female nontumor tissues. After eliminating the identified discriminating metabolites associated with gender difference, we obtained a panel of 46 and 44 differentially expressed metabolites in PTC (Supplementary Table 1 in the SI) and BTA (Supplementary Table 2 in the SI) with VIP > 1 and P < 0.05, respectively. Common metabolic signatures, as characterized by increased glycolysis, amino acid metabolism, one-carbon metabolism, and tryptophan metabolism, were found in both PTC and BTA tumor tissues. Citric acid levels were decreased in both PTC and BTA tumor tissues, while the downstream products/intermediates generated in the citric acid cycle (also known as the tricarboxylic acid cycle [TCA] cycle), including succinic acid, fumaric acid, and malic acid, were increased or unchanged in both types of tumor tissues. Figure 2 shows a heat map of the fold changes of these metabolites. The two different types of tumor tissues also showed distinct metabolite profiles. Purine and pyrimidine metabolism were significantly elevated in the PTC tumor tissues, and taurine and hypotaurine were also elevated in the PTC tumor specimens. Increased fatty acid and bile acid levels were found in BTA

Statistical Analysis

Statistical analysis was performed mainly using SPSS (SPSS, Chicago, IL) and SIMCA-P 13.0 (Umetrics, Umea, Sweden). Principal component analysis (PCA) as well as orthogonal partial least-squares discriminant analysis (OPLS−DA) were performed to visualize the metabolic differences among the PTC, BTA, and LSCC groups. Significantly altered metabolites with the variable importance in the projection (VIP) threshold (VIP > 1) in the previously mentioned OPLS-DA model, as well as the Student’s t test (P < 0.05), were selected in the three groups. 3317

DOI: 10.1021/acs.jproteome.5b00351 J. Proteome Res. 2015, 14, 3315−3321

Article

Journal of Proteome Research

Figure 2. Heat map of the fold changes (FCs) of the common metabolites and pathways in the two types of tumor tissues compared with adjacent nontumor tissues, including glycolysis, amino acid metabolism, one-carbon metabolism, and tryptophan metabolism.

Figure 3. Heat map of the fold changes (FCs) of the differential metabolites and pathways in the two types of tumor tissues compared with adjacent nontumor tissues, including purine and pyrimidine, taurine and hypotaurine, fatty acid, and bile acid metabolism.

43 females, which is a reasonable ratio. We investigated the impact of gender using a small age-matched sample and identified a panel of 21 metabolites that were differentially expressed between PTC nontumor tissues from male and female patients. These metabolites were therefore removed from the key discriminating metabolites identified for PTC and BTA. The tumor development process requires an altered cell metabolic program to support its bioenergetics and biosynthetic needs for uncontrolled cellular proliferation.13 In this study, we found that PTC and BTA tumor tissues shared common metabolic signatures, as characterized by increased glycolytic

specimens. Figure 3 shows a heat map of the fold changes of these metabolites. A comparison of the metabolite profiles between the PTC and LSCC specimens showed that they shared common metabolic signatures. The only distinct metabolites between the PTC and LSCC tissues were taurine and hypotaurine, which were elevated in the PTC tissues but unchanged in the LSCC tissues.



DISCUSSION The prevalence of PTC is approximately 3- to 4-fold higher in women than in men. Our PTC group comprised 14 males and 3318

DOI: 10.1021/acs.jproteome.5b00351 J. Proteome Res. 2015, 14, 3315−3321

Article

Journal of Proteome Research

identified as anticarcinogens that characterizing and optimizing such defense systems may be an important part of a strategy of minimizing cancer and other age-related diseases.29 Taurine and its derivatives such as taurolidine and taurochloramin were found to display antineoplastic effect both in vitro and in vivo, through suppressing cell proliferation, enhancement of tumor cell apoptosis,30,31 and through an antiangiogenic effects,32 while enhancing the therapeutic index of some antitumor agent.30 Because of the multiple functions of taurine in tumor cells, its exact role in head and neck cancer is unclear. In our result, the only distinct metabolites between the PTC and LSCC tissues were taurine and hypotaurine, which were elevated in the PTC tissues but unchanged in the LSCC tissues. Yao8 reported major metabolic differences between PTC and BTA in serum with regard to lipid metabolism, where the level of 3-hydroxybutyric acid, an intermediate product of fatty acid metabolism, was much higher in PTC than in BTA. Circulating free fatty acids are important for energy supplementation, particularly when the glucose supply is insufficient. Interestingly, higher levels of six free fatty acids (docosahexaenoic acid, 9-octadecenic acid, adrenic acid, docosapentaenoic acid, palmitoleic acid, and hexadecanoic acid) were observed in BTA tissues in our study, while only eicosadienoic acid was detected at high levels in PTC tissues. These results suggest that benign tumor tissue has higher impact on fatty acid metabolism compared with malignant ones. Oxidation of fatty acids may deliver bioenergy for rapid tumor cell proliferation and growth.33 Therefore, the lower lipid levels of fatty acids in cancer tissues could be associated with increased demand of lipids for membrane biosynthesis of tumor cells, leading to a higher utilization rate of lipids. In addition to altered fatty acid levels, an increased level of bile acids was observed in BTA, a result that may be correlated with higher fatty acid metabolism in BTA. Bile acids contribute to the absorption of dietary fats and induce cellular signaling pathways that affect cell proliferation.34

activity, amino acid metabolism, one-carbon metabolism and tryptophan metabolism, resulting from increased production of metabolic substrates and macromolecular precursors to meet the increased energy and protein synthesis needs for cell proliferation. The Warburg effect is a characteristic feature of tumor metabolism and characterized by an increased rate of glycolysis, followed by lactic acid fermentation during tumor growth.14 Coinciding with the Warburg effect, we found higher levels of lactic acid and TCA cycle intermediates in both PTC and BTA. Citric acid is a TCA cycle product, which was reduced significantly in all thyroid lesions in Deja’s study.7 We also found that citric acid was reduced in PTC and BTA tissues, while the downstream products that are generated in TCA cycle were increased or unchanged in both types of tumor tissues. Reduced citric acid could be explained by an increased conversion of citric acid into acetyl-CoA and further fatty acid biosynthesis. Amino acids are critical substrates needed for mitochondrial metabolism and were usually found to be elevated in tumor tissues. Proline levels are increased in multiple types of cancer, including prostate, lymphoma, and others.15−17 Aspartate is important for the growth of human pancreatic adenocarcinoma.18 Secretion of alanine is higher in melanoma cell lines than in normal melanocytes and is also significant in human colon carcinoma.19,20 Amino acids were increased in both types of tumor tissues in the present study, which may be a substrate for protein synthesis, an essential starting point for nucleotide biosynthesis and replenish TCA cycle metabolite levels to enable their use in anabolic processes. One-carbon metabolism involves the transfer of a carbon unit from methyl donor nutrients to molecules involved in the synthesis and methylation of DNA. Cysteine, an important precursor in glutathione production, has been shown to be elevated in various types of cancers including breast, ovarian, oral, brain, and lung cancers.21 In our study, we also detected an increased level of one-carbon metabolism in both types of tumor tissues. Tryptophan, an essential amino acid, is metabolized in the local microenvironment of tumors through the kynurenine pathway, which is first catalyzed by indoleamine 2,3-dioxygenase (IDO).22 The elevated expression of IDO was suspected to be a mediator of tumor immune tolerance, which may help tumor cells avoid immune attack.23 The increased expression of IDO may result in a higher level of kynurenine, which would ultimately generate more nicotinamide adenosine dinucleotide from tryptophan to meet accelerated growth of tumor cells. Very recently kynurenine was shown in gliomas to be an endogenous ligand of the human aryl hydrocarbon receptor and to suppress antitumor immune responses and also to promote tumor cell survival and motility through the aryl hydrocarbon receptor in an autocrine/paracrine fashion.24 Significant increase in IDO gene expression was observed in the colorectal cancer11 and thyroid cancer.25 In our results, tryptophan metabolism was disrupted not only in malignant tumor (PTC) but also in benign tumor (BTA). Another important amino acid, taurine, acts as an osmoregulator, is similar to scyllo- and myo-inositol. It was reported to be increased in breast cancer,26 colon cancer,27 and head and neck cancer cells.9 This suggests that changes in osmoregulation may be a common metabolic feature for many cancers. Taurine was found elevated in PTC and unchanged in BTA and LSCC in our result. In addition to osmoregulation, taurine is also found to display antioxidant properties.28 Many antioxidants are being



CONCLUSIONS We identified the metabolic profiles of PTC and BTA using a mass-spectrometry-based metabolomics approach, which demonstrated important similarities and differences in specific metabolic pathways. All tumor tissues showed increased glycolysis, amino acid metabolism, one-carbon metabolism, and tryptophan metabolism, suggesting that these two types of tumors share common metabolic pathways resulting from the increased need for energy, macromolecular precursors, and proteins and lipid synthesis. Purine and pyrimidine metabolism was higher particularly in the PTC specimens, and taurine and hypotaurine levels were also higher in PTC tissues. Increased fatty acid and bile acids levels were found in BTA. Taken together, the metabolic profiles of PTC and BTA tissues are composed of similar and remarkably different metabolites, suggesting the presence of common and unique mechanistic pathways in these types of tumors during tumorigenesis.



ASSOCIATED CONTENT

S Supporting Information *

Supplementary Table 1. Identities of 54 differential metabolites between PTC and adjacent normal controls from each of the tissue samples. Supplementary Table 2. Identities of 49 differential metabolites between BTA and adjacent normal controls from each of the tissue samples. Supplementary Table 3. Identities of 21 differential metabolites between male and female 3319

DOI: 10.1021/acs.jproteome.5b00351 J. Proteome Res. 2015, 14, 3315−3321

Article

Journal of Proteome Research

(12) Zheng, X.; Xie, G.; Zhao, A.; Zhao, L.; Yao, C.; Chiu, N. H.; Zhou, Z.; Bao, Y.; Jia, W.; Nicholson, J. K.; Jia, W. The footprints of gut microbial-mammalian co-metabolism. J. Proteome Res. 2011, 10, 5512− 22. (13) Erickson, J. W.; Cerione, R. A. Glutaminase: a hot spot for regulation of cancer cell metabolism? Oncotarget 2010, 1, 734−740. (14) Gatenby, R. A.; Gillies, R. J. Why do cancers have high aerobic glycolysis? Nat. Rev. Cancer 2004, 4, 891−9. (15) Liu, W.; Le, A.; Hancock, C.; Lane, A. N.; Dang, C. V.; Fan, T. W.; Phang, J. M. Reprogramming of proline and glutamine metabolism contributes to the proliferative and metabolic responses regulated by oncogenic transcription factor c-MYC. Proc. Natl. Acad. Sci. U. S. A. 2012, 109, 8983−8. (16) De Ingeniis, J.; Ratnikov, B.; Richardson, A. D.; Scott, D. A.; AzaBlanc, P.; De, S. K.; Kazanov, M.; Pellecchia, M.; Ronai, Z.; Osterman, A. L.; Smith, J. W. Functional specialization in proline biosynthesis of melanoma. PLoS One 2012, 7, e45190. (17) Nilsson, R.; Jain, M.; Madhusudhan, N.; Sheppard, N. G.; Strittmatter, L.; Kampf, C.; Huang, J.; Asplund, A.; Mootha, V. K. Metabolic enzyme expression highlights a key role for MTHFD2 and the mitochondrial folate pathway in cancer. Nat. Commun. 2014, 5, 3128. (18) Son, J.; Lyssiotis, C. A.; Ying, H.; Wang, X.; Hua, S.; Ligorio, M.; Perera, R. M.; Ferrone, C. R.; Mullarky, E.; Shyh-Chang, N.; Kang, Y.; Fleming, J. B.; Bardeesy, N.; Asara, J. M.; Haigis, M. C.; DePinho, R. A.; Cantley, L. C.; Kimmelman, A. C. Glutamine supports pancreatic cancer growth through a KRAS-regulated metabolic pathway. Nature 2013, 496, 101−5. (19) Weinberg, F.; Hamanaka, R.; Wheaton, W. W.; Weinberg, S.; Joseph, J.; Lopez, M.; Kalyanaraman, B.; Mutlu, G. M.; Budinger, G. R.; Chandel, N. S. Mitochondrial metabolism and ROS generation are essential for Kras-mediated tumorigenicity. Proc. Natl. Acad. Sci. U. S. A. 2010, 107, 8788−93. (20) Holm, E.; Hagmuller, E.; Staedt, U.; Schlickeiser, G.; Gunther, H. J.; Leweling, H.; Tokus, M.; Kollmar, H. B. Substrate balances across colonic carcinomas in humans. Cancer Res. 1995, 55, 1373−1378. (21) Gamcsik, M. P.; Kasibhatla, M. S.; Teeter, S. D.; Colvin, O. M. Glutathione levels in human tumors. Biomarkers 2012, 17, 671−91. (22) Terentis, A. C.; Thomas, S. R.; Takikawa, O.; Littlejohn, T. K.; Truscott, R. J.; Armstrong, R. S.; Yeh, S. R.; Stocker, R. The heme environment of recombinant human indoleamine 2,3-dioxygenase. Structural properties and substrate-ligand interactions. J. Biol. Chem. 2002, 277, 15788−94. (23) Zamanakou, M.; Germenis, A. E.; Karanikas, V. Tumor immune escape mediated by indoleamine 2,3-dioxygenase. Immunol. Lett. 2007, 111, 69−75. (24) Opitz, C. A.; Litzenburger, U. M.; Sahm, F.; Ott, M.; Tritschler, I.; Trump, S.; Schumacher, T.; Jestaedt, L.; Schrenk, D.; Weller, M.; Jugold, M.; Guillemin, G. J.; Miller, C. L.; Lutz, C.; Radlwimmer, B.; Lehmann, I.; von Deimling, A.; Wick, W.; Platten, M. An endogenous tumourpromoting ligand of the human aryl hydrocarbon receptor. Nature 2011, 478, 197−203. (25) Moretti, S.; Menicali, E.; Voce, P.; Morelli, S.; Cantarelli, S.; Sponziello, M.; Colella, R.; Fallarino, F.; Orabona, C.; Alunno, A.; de Biase, D.; Bini, V.; Mameli, M. G.; Filetti, S.; Gerli, R.; Macchiarulo, A.; Melillo, R. M.; Tallini, G.; Santoro, M.; Puccetti, P.; Avenia, N.; Puxeddu, E. Indoleamine 2,3-dioxygenase 1 (IDO1) is up-regulated in thyroid carcinoma and drives the development of an immunosuppressant tumor microenvironment. J. Clin. Endocrinol. Metab. 2014, 99, E832−40. (26) Bathen, T. F.; Jensen, L. R.; Sitter, B.; Fjosne, H. E.; Halgunset, J.; Axelson, D. E.; Gribbestad, I. S.; Lundgren, S. MR-determined metabolic phenotype of breast cancer in prediction of lymphatic spread, grade, and hormone status. Breast Cancer Res. Treat. 2007, 104, 181−9. (27) Tessem, M. B.; Selnaes, K. M.; Sjursen, W.; Trano, G.; Giskeodegard, G. F.; Bathen, T. F.; Gribbestad, I. S.; Hofsli, E. Discrimination of patients with microsatellite instability colon cancer using 1H HR MAS MR spectroscopy and chemometric analysis. J. Proteome Res. 2010, 9, 3664−70.

adjacent nontumor tissues from PTC samples.The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jproteome.5b00351.



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. Fax: 086-21-53882205. Tel: 086-21-53882271 (J.W.) *E-mail: [email protected]. Fax: 086-21-64834143. Tel.: 086-21-24058706 (S.Y.) Author Contributions ⊥

J.W. and S.Y. contributed equally to this work.

Notes

The authors declare no competing financial interest.



ABBREVIATIONS PTC, papillary thyroid carcinoma; BTA, benign thyroid adenoma; LSCC, laryngeal squamous cell carcinoma; TCA, tricarboxylic acid; GC−TOFMS, gas chromatography−time-offlight mass spectrometry; UPLC−QTOFMS, ultraperformance liquid chromatography−quadrupole time-of-flight mass spectrometry; PCA, principal component analysis; OPLS−DA, orthogonal partial least-squares−discriminant analysis; VIP, variable importance in the projection



REFERENCES

(1) Siegel, R. L.; Miller, K. D.; Jemal, A. Cancer statistics, 2015. CaCancer J. Clin. 2015, 65, 5−29. (2) Reiners, C.; Wegscheider, K.; Schicha, H.; Theissen, P.; Vaupel, R.; Wrbitzky, R.; Schumm-Draeger, P. M. Prevalence of thyroid disorders in the working population of Germany: ultrasonography screening in 96,278 unselected employees. Thyroid 2004, 14, 926−32. (3) Brito, J. P.; Yarur, A. J.; Prokop, L. J.; McIver, B.; Murad, M. H.; Montori, V. M. Prevalence of thyroid cancer in multinodular goiter versus single nodule: a systematic review and meta-analysis. Thyroid 2013, 23, 449−55. (4) Harach, H. R.; Franssila, K. O.; Wasenius, V. M. Occult papillary carcinoma of the thyroid. A ″normal″ finding in Finland. A systematic autopsy study. Cancer 1985, 56, 531−8. (5) Udelsman, R.; Chen, H. The current management of thyroid cancer. Adv. Surg 1999, 33, 1−27. (6) van Dijk, B. A.; Karim-Kos, H. E.; Coebergh, J. W.; Marres, H. A.; de Vries, E. Progress against laryngeal cancer in The Netherlands between 1989 and 2010. Int. J. Cancer 2014, 134, 674−81. (7) Deja, S.; Dawiskiba, T.; Balcerzak, W.; Orczyk-Pawilowicz, M.; Glod, M.; Pawelka, D.; Mlynarz, P. Follicular adenomas exhibit a unique metabolic profile. (1)H NMR studies of thyroid lesions. PLoS One 2013, 8, e84637. (8) Yao, Z.; Yin, P.; Su, D.; Peng, Z.; Zhou, L.; Ma, L.; Guo, W.; Ma, L.; Xu, G.; Shi, J.; Jiao, B. Serum metabolic profiling and features of papillary thyroid carcinoma and nodular goiter. Mol. BioSyst. 2011, 7, 2608−14. (9) Tripathi, P.; Kamarajan, P.; Somashekar, B. S.; MacKinnon, N.; Chinnaiyan, A. M.; Kapila, Y. L.; Rajendiran, T. M.; Ramamoorthy, A. Delineating metabolic signatures of head and neck squamous cell carcinoma: phospholipase A2, a potential therapeutic target. Int. J. Biochem. Cell Biol. 2012, 44, 1852−61. (10) Somashekar, B. S.; Kamarajan, P.; Danciu, T.; Kapila, Y. L.; Chinnaiyan, A. M.; Rajendiran, T. M.; Ramamoorthy, A. Magic angle spinning NMR-based metabolic profiling of head and neck squamous cell carcinoma tissues. J. Proteome Res. 2011, 10, 5232−41. (11) Qiu, Y.; Cai, G.; Zhou, B.; Li, D.; Zhao, A.; Xie, G.; Li, H.; Cai, S.; Xie, D.; Huang, C.; Ge, W.; Zhou, Z.; Xu, L. X.; Jia, W.; Zheng, S.; Yen, Y.; Jia, W. A distinct metabolic signature of human colorectal cancer with prognostic potential. Clin. Cancer Res. 2014, 20, 2136−46. 3320

DOI: 10.1021/acs.jproteome.5b00351 J. Proteome Res. 2015, 14, 3315−3321

Article

Journal of Proteome Research (28) El Agouza, I. M.; Eissa, S. S.; El Houseini, M. M.; El-Nashar, D. E.; Abd El Hameed, O. M. Taurine: a novel tumor marker for enhanced detection of breast cancer among female patients. Angiogenesis 2011, 14, 321−30. (29) Lin, J.; Manson, J. E.; Selhub, J.; Buring, J. E.; Zhang, S. M. Plasma cysteinylglycine levels and breast cancer risk in women. Cancer Res. 2007, 67, 11123−7. (30) Sadzuka, Y.; Matsuura, M.; Sonobe, T. The effect of taurine, a novel biochemical modulator, on the antitumor activity of doxorubicin. Biol. Pharm. Bull. 2009, 32, 1584−7. (31) Daigeler, A.; Chromik, A. M.; Geisler, A.; Bulut, D.; Hilgert, C.; Krieg, A.; Klein-Hitpass, L.; Lehnhardt, M.; Uhl, W.; Mittelkotter, U. Synergistic apoptotic effects of taurolidine and TRAIL on squamous carcinoma cells of the esophagus. Int. J. Oncol. 2008, 32, 1205−20. (32) Rodak, R.; Kubota, H.; Ishihara, H.; Eugster, H. P.; Konu, D.; Mohler, H.; Yonekawa, Y.; Frei, K. Induction of reactive oxygen intermediates-dependent programmed cell death in human malignant ex vivo glioma cells and inhibition of the vascular endothelial growth factor production by taurolidine. J. Neurosurg. 2005, 102, 1055−68. (33) Kuhajda, F. P. Fatty-acid synthase and human cancer: new perspectives on its role in tumor biology. Nutrition 2000, 16, 202−8. (34) Cronin, J.; Williams, L.; McAdam, E.; Eltahir, Z.; Griffiths, P.; Baxter, J.; Jenkins, G. The role of secondary bile acids in neoplastic development in the oesophagus. Biochem. Soc. Trans. 2010, 38, 337−42.

3321

DOI: 10.1021/acs.jproteome.5b00351 J. Proteome Res. 2015, 14, 3315−3321