Metabolome analysis reveals a diversity of cancer tissues in advanced epithelial ovarian cancer
Cancer Cell International volume 21, Article number: 314 (2021)
Epithelial ovarian cancer remains one of the leading causes of cancer deaths among women worldwide, and advanced epithelial ovarian cancer frequently metastasizes to the omentum. The characteristics of metastatic cancer may differ from those of primary ovarian cancer and reflect the unique omental microenvironment. This study investigated metabolomic differences in epithelial ovarian cancers.
Patients with advanced epithelial ovarian cancer were eligible for this study. Five patients underwent surgery and resection of paired primary ovarian and omental metastatic cancer at Nagoya University. Metabolome analysis was performed in these paired cancer and metastatic cancer tissues through a facility service (C-SCOPE) at Human Metabolome Technologies, Inc. The concentrations of 116 compounds were measured by CE-TOFMS and CE-QqQMS, and 30 metabolic parameters were calculated. For statistical analyses, Welch’s t-test was used for comparisons between two independent groups.
Metabolite profiles were all different, which reflects diversity among these cancer tissues. Of the measured compounds, urea was the only metabolite that was significantly decreased in omental metastatic cancers compared with the primary cancers (p = 0.031). Moreover, in omental metastatic cancers, the pentose phosphate pathway was more dominant than glycolysis. Furthermore, in some cases, lactic acids in omental metastatic cancers were markedly decreased compared with primary cancers. With regard to histological subtype, the total levels of amino acids, especially the percentage of glutamine, were significantly enriched in serous carcinomas compared with nonserous carcinomas (p = 0.004 and p = 0.001). Moreover, the reduced forms of glutathione and polyamines were also more abundant in serous carcinomas than in nonserous carcinomas (p = 0.025 and 0.048).
The metabolite profiles differed depending on tumor location and histological subtype. Metabolome analysis may be a useful tool for identifying cancer diagnostic and prognostic markers.
Epithelial ovarian cancer (EOC) is one of the leading causes of cancer-related deaths among women worldwide with an estimated 184,799 deaths in 2018 . EOC comprises four major histological subtypes: high-grade serous ovarian carcinoma (HGSOC), clear cell carcinoma (CCC), endometrioid carcinoma (EC), and mucinous carcinoma, and the clinical and molecular biological features are different for each histological subtype [2,3,4]. The tumor microenvironment plays a critical role in cancer progression, and EOC typically spreads to the peritoneal cavity and frequently metastasizes to the omentum, which is characterized by an adipose-rich environment [4,5,6,7,8]. Several studies have indicated that adipocytes in the omentum contribute to cancer development by producing cytokines and chemokines [8,9,10]. This implies that the molecular biological characteristics of the cancer cells may change, thus reflecting their unique microenvironment. Therefore, understanding the differences is essential for elucidating EOC pathology.
Since cancer cell-specific metabolism, which is known as the Warburg effect, was discovered, energy metabolism in cancer tissue has been considered to be important [11,12,13,14]. Metabolome analysis represents a robust method that is used to identify and quantitate the levels of small molecules related to cellular metabolic activity, such as sugars, amino acids, and other organic molecules . In molecular biology, the gene expression level is usually regarded as a surrogate for predicting gene function. However, in the case of enzymatic reactions that are part of a metabolic pathway, the enzyme expression level and its function are not necessarily correlated [16, 17]. Therefore, metabolome analysis is suitable for evaluating the functional differences in the microenvironment because it enables analysis of the metabolic changes that reflect enzyme function. Several studies have indicated that serum and tissue metabolites may be potential diagnostic and prognostic biomarkers of EOC [18,19,20,21]. However, metabolomic differences based on the spatial diversity of the tumor microenvironment remain undefined.
In this study, we investigated the metabolic changes in primary ovarian cancer (PC) and omental metastatic cancer (OMC) and the metabolic differences associated with the various histological subtypes of EOC.
Patients and samples
Archival fresh frozen EOC tissue samples were collected from five patients who underwent simultaneous PC and OMC tissue resection at the Nagoya University Hospital (Aichi, Japan) since 2018. Patient characteristics are provided in Table 1. Histologically, Cases 1 and 2 were HGSOC, Cases 3 and 4 were CCC, and Case 5 was EC. Case 3 underwent interval debulking surgery after eight courses of carboplatin and paclitaxel combination chemotherapy, whereas the remaining four patients underwent primary debulking surgery. The study protocol was approved by the Ethics Committee of our institution (approval No. 2017-0497), and written informed consent was obtained from each patient. Immediately following resection, cancer tissues were macroscopically sectioned into small pieces and placed into 1.5-mL tubes. They were immediately frozen in liquid nitrogen and stored at − 80 ℃. Approximately 40–50 mg of tissues from each case was used for the metabolome analysis.
Metabolome analysis was performed using a facility service (C-SCOPE) at Human Metabolome Technologies, Inc. (Yamagata, Japan). Briefly, the tissues were homogenized in a 50% acetonitrile aqueous solution. The samples were then centrifuged at 2,300 × g at 4˚C for 5 min and the upper aqueous layer was centrifugally filtered through a 5 kDa cutoff filter at 9,100 × g for 120 min at 4˚C. The filtrates were resuspended in 50 μl of Milli-Q water and the metabolome analysis was performed using CE-TOFMS and CE-QqQMS. The data were analyzed by MasterHands ver.220.127.116.11 (Keio University, Tsuruoka, Japan) and MassHunter Quantitative Analysis B.06.00 (Agilent Technologies, Santa Clara, CA), and the peak area of each metabolite was calculated. The peak area was then normalized to an internal standard, and the metabolite concentrations were calculated using standard curves. Overall, 116 compounds were measured and 30 metabolic parameters were calculated.
Data analysis and statistical analysis
RStudio (RStudio, Boston, MA) and R software (ver. 3.5.0) were used for the analysis. For heatmap and principal component analyses, values below the detection sensitivity were designated as zero, and the measurements of substances for which the concentration was zero were excluded from all samples. The data was converted to base 10 logarithms and z-scores. The heatmap.2 function of the gplots package (ver. 3.0.1) was used to generate a heatmap. The distance was calculated as “1–Spearman correlation coefficient” and the “ward.D2” clustering method was used. To calculate and visualize the principal component analysis, the prcomp and plot3d functions of the rgl package (ver. 0.99.16) were used. For statistical analyses, Welch’s t-test was used for comparisons between PC and OMC. In addition, metabolite changes between HGSOC and non-HGSOC were compared using Welch’s t-test. P-values less than 0.05 was considered statistically significant.
We performed a metabolome analysis using 10 fresh-frozen EOC tissues from five patients in which 116 compounds were measured (Additional file 1). A heatmap analysis revealed that metabolites from the paired PC and OMC tissues were similar and that metabolites in HGSOC subtype (Cases 1 and 2) were quite different from those in the CCC and EC subtypes (Cases 3–5, Fig. 1A). However, no significant differences were observed in the metabolite profiles between the PC and OMC groups. Likewise, the principal component analysis revealed a diversity in the cancer tissues based on histological subtypes and individual differences (Fig. 1B).
The urea cycle
We investigated the detailed differences in multiple metabolic pathways between the PC and OMC groups. The concentrations of most of the compounds were not significantly different between the two groups (Additional file 1); however, the urea concentration in the OMC group was significantly decreased compared with the PC group (p = 0.031). Notably, urea was decreased by approximately half in Cases 1, 4, and 5 (Fig. 2A). Therefore, other compounds in the urea cycle were subsequently evaluated. The concentration of ornithine in the OMC group was also decreased in all cases as both urea and ornithine are derived from arginine (Fig. 2B). In addition to arginase, nitric oxide synthase (NOS) is a metabolic enzyme of arginine, and the citrulline/ornithine ratio is considered to be an index of the enzyme activities of arginase and NOS (Fig. 2C) . The citrulline/ornithine ratio was elevated in all cases in the OMC group, although citrulline alone was not significantly altered (Fig. 2B, 2D). Therefore, NOS activity was increased in the OMC group, but its levels were not significantly different from the other groups (p = 0.352, Additional file 2). No significant differences were observed in other compounds (Fig. 2B).
Glycolysis and the pentose phosphate pathway (PPP)
Next, we focused on glycometabolism. The glucose-6-phosphate (G6P)/ribose-5-phosphate (R5P) ratio was decreased in the OMC group compared with the paired PC group, but the difference was not significant (p = 0.204, Fig. 3A and Additional file 2). G6P was somewhat decreased in the OMC group, especially in Case 3, whereas R5P was increased in the OMC group, particularly in Case 5 (Fig. 3B). This suggested that the PPP was dominant in the OMC group. However, the concentrations of NADP + and NADPH were almost identical in all the samples (Additional file 1).
Pyruvic acid is produced as a product of glycolysis and is metabolized to lactic acid and acetyl-CoA. Lactic acid was remarkably decreased in the OMC group in Cases 1, 3, and 5 (Fig. 3B and Additional file 1). Moreover, the lactate/pyruvate ratio was decreased in the OMC group compared with the paired PC group in all cases except Case 4 (Fig. 3C). However, due to individual differences, these differences were not statistically significant (p = 0.248). This suggested that the oxygen environment in OMC is different than that in PC.
Amino acid composition
In addition to the PC and OMC groups, we evaluated metabolite differences according to histological subtype. When HGSOC (Cases 1 and 2) was compared with non-HGSOC (Cases 3–5), HGSOC samples contained significantly more amino acids (p = 0.004, Fig. 4A). In all samples, glutamic acid was the major component, whereas glycine and alanine were also abundant (Fig. 4B). Notably, the proportion of glutamine was significantly higher in the HGSOC than in non-HGSOC samples (p = 0.001, Fig. 4B). As previously described, HGSOC samples contained more total amino acids, and thus, the absolute concentration of glutamine was also significantly higher (HGSOC vs. non-HGSOC, p = 0.004). Therefore, glutathione metabolism-related amino acids were enriched in HGSOC compared with non-HGSOC. (Fig. 4C). Glutathione has two forms; a reduced form (GSH) and an oxidized disulfide form (GSSG). GSH is converted into GSSG during oxidative stress . The samples HGSOC contained significantly more GSH than the non-HGSOC (p = 0.025), whereas no differences in GSSG were observed (Fig. 4D). Ornithine, a metabolite of the urea cycle, is metabolized to polyamines including putrescine, spermidine, and spermine (Fig. 4E). The sum of these three polyamine forms was significantly higher in HGSOC samples than in non-HGSOC samples (p = 0.048, Fig. 4F).
Metabolism represents the set of life-sustaining chemical reactions in cells and is essential for cancer progression. In general, gene expression is regulated by various factors, such as microRNAs and 14–3-3 proteins [24,25,26]. Therefore, enzyme activity cannot be determined from the gene expression level alone, and analysis of metabolites is required [15, 17].
First, we showed that the urea concentration was significantly decreased in the OMC group compared with the PC group and that alterations in the urea cycle were observed. The urea cycle is an essential pathway for the disposal of ammonia, and thus, the relationship between alterations in the urea cycle and cancer progression has been presented [27, 28]. Our results indicated that the citrulline/ornithine ratio was increased in OMC and is considered to reflect decreased arginase activity and increased NOS activity. Consequently, NO was also considered to be increased. The function of NO is somewhat ambivalent, but NO primarily promotes cancer-related events, such as angiogenesis, apoptosis, and cell cycle . Unfortunately, in this analysis, NO was not a target compound and further studies are needed to evaluate urea cycle dysfunction and NO during EOC progression.
Next, we focused on glycometabolism, which is essential for cell survival and unique to cancer tissues according to the Warburg effect [11,12,13,14]. In addition to glycolysis, the PPP is also a major pathway for glucose catabolism and plays an important role in cancer cell survival and growth [14, 30,31,32]. The PPP is divided into two biochemical branches known as oxidative and non-oxidative PPP. In the first, the chemical reaction is irreversible and two molecules of NADPH are generated from one molecule of glucose, whereas, in the other phase, which is reversible, the five-carbon sugar is generated from a six-carbon glucose molecule. Therefore, depending on which phase of the PPP is activated, the production of either NADPH or R5P is regulated [30,31,32]. In our results, the low G6P/R5P ratio in OMC indicated the predominance of the PPP. However, NADPH concentrations were similar in all the samples. Therefore, NAPDH consumption may have been increased in OMC due to oxidative stress. Alternatively, the OMC cells may produce more R5P than NADPH by activating the non-oxidative phase of the PPP. In addition, enhanced PPP metabolism might be due to the increased proportion of cancer stem cells in OMC tissues [30, 33]. Next, we evaluated glycolysis as the Warburg effect predicts that cancer cells use glycolysis even under aerobic conditions [11, 14]. We discovered that the lactate/pyruvate ratio was decreased in OMC, which suggested that the oxygen environment is different at metastatic sites and that cancer cells may change energy metabolism process based on their environment. Alterations of glycolysis between PC and OMC were also previously reported, although those detailed metabolic changes were not consistent with our results . Therefore, the heterogeneity of cancer cells and individual differences may be significant and further larger-scale studies are needed.
We demonstrated that HGSOCs contained significantly more glutamine and other amino acids than non-serous carcinomas. Glutamine was the most prominent amino acid, and since the 1950s, it has been reported that proliferating cells require glutamine [14, 35]. In addition, highly proliferative cells produce high levels of reactive oxygen species (ROS) . NADPH and GSH play important roles in the maintenance of intracellular redox balance, and elevated GSH levels are observed in various tumor types [23, 37]. We observed high glutamine and GSH levels in HGSOC, which may reflect a high proliferative potential. Moreover, GSH is known to mediate resistance to platinum analogs through several mechanisms including increased DNA repair and the inhibition of apoptosis [37, 38]. In addition, polyamines were also abundant in HGSOC. Polyamines are involved in many fundamental processes of cell growth and survival, and polyamine metabolism is associated with cancer-driving pathways . In addition, a previous study reported that cyclin E1-driven EOC was associated with activated polyamine synthesis and decreased cancer immunity . Therefore, GSH and polyamines are potential therapeutic targets in EOC, especially HGSOC.
This metabolome-based study has several limitations. First, we included only five cases of EOCs. Hence, although it is difficult to make any significant conclusion, the results presented in this study are very interesting. Therefore, further evaluation is essential to elucidate the metabolic changes in EOCs. Second, detailed molecular biological changes such as genetic mutations, mRNA and protein expression, and enzyme activities were not evaluated. Third, a lipidomics analysis was not performed, although this would be an interesting approach because the omentum is characterized by an adipose-rich microenvironment. However, lipid metabolites are not measurable with of C-SCOPE. Therefore, any metabolomic changes related to those molecules are worth evaluating in the future.
Metabolome analysis revealed different metabolite profiles between PC and OMC, particularly in the urea cycle, glycolysis, and the pentose phosphate pathway. Moreover, HGSOC samples contained more amino acids and GSH than non-serous carcinoma samples. These differences may reflect the diversity of biological processes in cancer cells. Metabolome analysis is expected to be a powerful tool in cancer biology.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its Additional files.
Epithelial ovarian cancer
High-grade serous ovarian carcinoma
Clear cell carcinoma
Primary ovarian cancer
Omental metastatic cancer
Nitric oxide synthase
Pentose phosphate pathway
Reactive oxygen species
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Nagase S, Ohta T, Takahashi F, Enomoto T. Annual report of the committee on gynecologic oncology, the Japan Society of Obstetrics and Gynecology: annual patients report for 2015 and annual treatment report for 2010. J Obstet Gynaecol Res. 2019;45(2):289–98.
Jayson GC, Kohn EC, Kitchener HC, Ledermann JA. Ovarian cancer. Lancet (London, England). 2014;384(9951):1376–88.
Cannistra SA. Cancer of the ovary. N Engl J Med. 2004;351(24):2519–29.
Luo Z, Wang Q, Lau WB, Lau B, Xu L, Zhao L, Yang H, Feng M, Xuan Y, Yang Y, et al. Tumor microenvironment: the culprit for ovarian cancer metastasis? Cancer Lett. 2016;377(2):174–82.
Lengyel E. Ovarian cancer development and metastasis. Am J Pathol. 2010;177(3):1053–64.
Motohara T, Masuda K, Morotti M, Zheng Y, El-Sahhar S, Chong KY, Wietek N, Alsaadi A, Karaminejadranjbar M, Hu Z, et al. An evolving story of the metastatic voyage of ovarian cancer cells: cellular and molecular orchestration of the adipose-rich metastatic microenvironment. Oncogene. 2019;38(16):2885–98.
Rodriguez GM, Galpin KJC, McCloskey CW, Vanderhyden BC. The tumor microenvironment of epithelial ovarian cancer and its influence on response to immunotherapy. Cancers. 2018;10(8):242.
Thibault B, Castells M, Delord JP, Couderc B. Ovarian cancer microenvironment: implications for cancer dissemination and chemoresistance acquisition. Cancer Metastasis Rev. 2014;33(1):17–39.
Nieman KM, Kenny HA, Penicka CV, Ladanyi A, Buell-Gutbrod R, Zillhardt MR, Romero IL, Carey MS, Mills GB, Hotamisligil GS, et al. Adipocytes promote ovarian cancer metastasis and provide energy for rapid tumor growth. Nat Med. 2011;17(11):1498–503.
Warburg O. On the origin of cancer cells. Science (New York, NY). 1956;123(3191):309–14.
Lunt SY, Vander Heiden MG. Aerobic glycolysis: meeting the metabolic requirements of cell proliferation. Annu Rev Cell Dev Biol. 2011;27:441–64.
Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science (New York, NY). 2009;324(5930):1029–33.
Counihan JL, Grossman EA, Nomura DK. Cancer metabolism: current understanding and therapies. Chem Rev. 2018;118(14):6893–923.
Beger RD. A review of applications of metabolomics in cancer. Metabolites. 2013;3(3):552–74.
Ter Kuile BH, Westerhoff HV. Transcriptome meets metabolome: hierarchical and metabolic regulation of the glycolytic pathway. FEBS Lett. 2001;500(3):169–71.
Griffin JL, Shockcor JP. Metabolic profiles of cancer cells. Nat Rev Cancer. 2004;4(7):551–61.
Plewa S, Horala A, Derezinski P, Klupczynska A, Nowak-Markwitz E, Matysiak J, Kokot ZJ. Usefulness of amino acid profiling in ovarian cancer screening with special emphasis on their role in cancerogenesis. Int J Mol Sci. 2017;18(12):2727.
Bachmayr-Heyda A, Aust S, Auer K, Meier SM, Schmetterer KG, Dekan S, Gerner C, Pils D. Integrative systemic and local metabolomics with impact on survival in high-grade serous ovarian cancer. Clin Cancer Res. 2017;23(8):2081–92.
Hilvo M, de Santiago I, Gopalacharyulu P, Schmitt WD, Budczies J, Kuhberg M, Dietel M, Aittokallio T, Markowetz F, Denkert C, et al. Accumulated metabolites of hydroxybutyric acid serve as diagnostic and prognostic biomarkers of ovarian high-grade serous carcinomas. Can Res. 2016;76(4):796–804.
Zhang F, Zhang Y, Ke C, Li A, Wang W, Yang K, Liu H, Xie H, Deng K, Zhao W, et al. Predicting ovarian cancer recurrence by plasma metabolic profiles before and after surgery. Metabolomics. 2018;14(5):65.
Jung C, Gonon AT, Sjoquist PO, Lundberg JO, Pernow J. Arginase inhibition mediates cardioprotection during ischaemia-reperfusion. Cardiovasc Res. 2010;85(1):147–54.
Bansal A, Simon MC. Glutathione metabolism in cancer progression and treatment resistance. J Cell Biol. 2018;217(7):2291–8.
Yoshida K, Yokoi A, Kato T, Ochiya T, Yamamoto Y. The clinical impact of intra- and extracellular miRNAs in ovarian cancer. Cancer Sci. 2020;111(10):3435–44.
De S, Marcinkiewicz JL, Vijayaraghavan S, Kline D. Expression of 14-3-3 protein isoforms in mouse oocytes, eggs and ovarian follicular development. BMC Res Notes. 2012;5(1):57.
Covington TR, De S. Binding partners of 14-3-3 (YWHA) protein isoforms among mammalian species, tissues, and developmental stages. Adv J Grad Res. 2021;10(1):16–22.
Morris SM Jr. Regulation of enzymes of the urea cycle and arginine metabolism. Annu Rev Nutr. 2002;22:87–105.
Keshet R, Szlosarek P, Carracedo A, Erez A. Rewiring urea cycle metabolism in cancer to support anabolism. Nat Rev Cancer. 2018;18(10):634–45.
Choudhari SK, Chaudhary M, Bagde S, Gadbail AR, Joshi V. Nitric oxide and cancer: a review. World J Surg Oncol. 2013;11:118.
Ghanbari Movahed Z, Rastegari-Pouyani M, Mohammadi MH, Mansouri K. Cancer cells change their glucose metabolism to overcome increased ROS: one step from cancer cell to cancer stem cell? Biomed Pharmacother. 2019;112:108690.
Patra KC, Hay N. The pentose phosphate pathway and cancer. Trends Biochem Sci. 2014;39(8):347–54.
Jiang P, Du W, Wu M. Regulation of the pentose phosphate pathway in cancer. Protein Cell. 2014;5(8):592–602.
Debeb BG, Lacerda L, Larson R, Wolfe AR, Krishnamurthy S, Reuben JM, Ueno NT, Gilcrease M, Woodward WA. Histone deacetylase inhibitor-induced cancer stem cells exhibit high pentose phosphate pathway metabolism. Oncotarget. 2016;7(19):28329–39.
Fong MY, McDunn J, Kakar SS. Identification of metabolites in the normal ovary and their transformation in primary and metastatic ovarian cancer. PLoS ONE. 2011;6(5):e19963.
Eagle H, Oyama VI, Levy M, Horton CL, Fleischman R. The growth response of mammalian cells in tissue culture to l-glutamine and l-glutamic acid. J Biol Chem. 1956;218(2):607–16.
Lieu EL, Nguyen T, Rhyne S, Kim J. Amino acids in cancer. Exp Mol Med. 2020;52(1):15–30.
Traverso N, Ricciarelli R, Nitti M, Marengo B, Furfaro AL, Pronzato MA, Marinari UM, Domenicotti C. Role of glutathione in cancer progression and chemoresistance. Oxid Med Cell Longevity. 2013;2013:972913.
Nunes SC, Serpa J. Glutathione in ovarian cancer: a double-edged sword. Int J Mol Sci. 2018;19(7):1882.
Casero RA Jr, Murray Stewart T, Pegg AE. Polyamine metabolism and cancer: treatments, challenges and opportunities. Nat Rev Cancer. 2018;18(11):681–95.
Guo T, Li B, Gu C, Chen X, Han M, Liu X, Xu C. PGC-1alpha inhibits polyamine metabolism in Cyclin E1-driven ovarian cancer. Cancer Med. 2019;8:7754–61.
This study was supported by the Program for Promoting the Enhancement of Research Universities as young researcher units for the advancement of new and undeveloped fields at Nagoya University. We thank Enago (https://www.enago.jp/) for English editing of this manuscript.
This study was supported by a Japan Society for the Promotion of Science Grant-in-Aid for Scientific Research: Grant Number 19K18693.
Ethics approval and consent to participate
The study protocol was approved by the Ethics Committee of our institution (approval No. 2017-0497). We obtained written informed consent from each patient.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yoshida, K., Yoshikawa, N., Kitami, K. et al. Metabolome analysis reveals a diversity of cancer tissues in advanced epithelial ovarian cancer. Cancer Cell Int 21, 314 (2021). https://doi.org/10.1186/s12935-021-02014-7
- Epithelial ovarian cancer
- Omental metastasis