Next Article in Journal
Ultramicronized Palmitoylethanolamide Inhibits NLRP3 Inflammasome Expression and Pro-Inflammatory Response Activated by SARS-CoV-2 Spike Protein in Cultured Murine Alveolar Macrophages
Next Article in Special Issue
A Longitudinal 1H NMR-Based Metabolic Profile Analysis of Urine from Hospitalized Premature Newborns Receiving Enteral and Parenteral Nutrition
Previous Article in Journal
Digesta and Plasma Metabolomics of Rainbow Trout Strains with Varied Tolerance of Plant-Based Diets Highlights Potential for Non-Lethal Assessments of Enteritis Development
Previous Article in Special Issue
NMR-Based Metabolomics Approach to Investigate the Effects of Fruits of Acanthopanax sessiliflorus in a High-Fat Diet Induced Mouse Model
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 11
Issue 9
10.3390/metabo11090591
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Unique Urinary Metabolic Feature for the Determination of Bladder Cancer, Prostate Cancer, and Renal Cell Carcinoma

by
Sujin Lee
1,
Ja Yoon Ku
2,
Byeong Jin Kang
3,
Kyung Hwan Kim
3,
Hong Koo Ha
4 and
Suhkmann Kim
1,*
1
Department of Chemistry and Chemistry Institute for Functional Materials, Institute for Plastic Information and Energy Materials, Pusan National University, Busandaehak-ro 63, Geumjeong-gu, Busan 46241, Korea
2
Department of Urology, Dongnam Institute of Radiological & Medical Sciences Cancer Center, Busan 46033, Korea
3
Department of Urology, College of Medicine, Pusan National University, Busan 49241, Korea
4
Department of Urology, College of Medicine, Pusan National University and Biomedical Research Institute, Pusan National University Hospital, Busan 49241, Korea
*
Author to whom correspondence should be addressed.
Metabolites 2021, 11(9), 591; https://doi.org/10.3390/metabo11090591
Submission received: 9 August 2021 / Revised: 28 August 2021 / Accepted: 31 August 2021 / Published: 2 September 2021
(This article belongs to the Special Issue NMR-Based Metabolic Profiling with Applications in Disease and Food Control)
Browse Figures
Versions Notes

Abstract

:
Prostate cancer (PCa), bladder cancer (BCa), and renal cell carcinoma (RCC) are the most prevalent cancer among urological cancers. However, there are no cancer-specific symptoms that can differentiate them as well as early clinical signs of urological malignancy. Furthermore, many metabolic studies have been conducted to discover their biomarkers, but the metabolic profiling study to discriminate between these cancers have not yet been described. Therefore, in this study, we aimed to investigate the urinary metabolic differences in male patients with PCa (n = 24), BCa (n = 29), and RCC (n = 12) to find the prominent combination of metabolites between cancers. Based on 1H NMR analysis, orthogonal partial least-squares discriminant analysis was applied to find distinct metabolites among cancers. Moreover, the ranked analysis of covariance by adjusting a potential confounding as age revealed that 4-hydroxybenzoate, N-methylhydantoin, creatinine, glutamine, and acetate had significantly different metabolite levels among groups. The receiver operating characteristic analysis created by prominent five metabolites showed the great discriminatory accuracy with area under the curve (AUC) > 0.7 for BCa vs. RCC, PCa vs. BCa, and RCC vs. PCa. This preliminary study compares the metabolic profiles of BCa, PCa, and RCC, and reinforces the exploratory role of metabolomics in the investigation of human urine.

Graphical Abstract

1. Introduction

Urological cancers include tumors of the kidney, bladder, and in men, also of the prostate and testis. Among these, prostate cancer (PCa), bladder cancer (BCa), and renal cell carcinoma (RCC) are the most frequent cancer types. In the United States of America, an estimated 191,930 new cases and 33,330 deaths from PCa and 62,110 new cases and 13,050 deaths from BCa have been reported [1]. In most cases, the early symptoms of urologic cancer do not appear until the cancer has progressed, and because the incidence is higher in men than women, these cancers are considered to be dangerous cancers for men. In addition, there are no cancer-specific clinical symptoms between PCa, BCa, and RCC. Urological cancers are accompanied by common symptoms such as dysuria, hematuria, and frequent urination, which may be a sign of cancer [2]. However, these symptoms can also appear with urinary tract infection [3]. Thus, diagnostic methods such as cystoscopy of the bladder, digital rectal examination of the prostate including biopsies, and computed tomography (CT) scans are needed to distinguish the three representative urological malignancies. However, these diagnoses are very invasive and require radiation exposure, rendering them unpleasant to patients as well as risking unnecessary complication [4,5,6]. Therefore, although these examinations remain the optimal norm for detecting urological cancers, given their invasiveness and high cost, non-invasive approaches for early discrimination between three urological cancers are essential.
Metabolomics involves the comprehensive systemic profiling of metabolite concentrations and their cellular perturbations that can be influenced by numerous factors, such as genetics, pathological conditions, and external factors (e.g., diet, lifestyle, and environment) [7]. This technique is closely related to the phenotype in biological systems, so is increasingly being applied toward a potential diagnostic tool to monitor disease status [8,9,10]. The main analytical platforms for metabolic profiling are nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry (MS). Both approaches are appropriate for metabolic analysis, but they have different analytical strengths and weaknesses [7]. MS can identify metabolites with high sensitivity. However, the matrix effect and ionization suppression can be influenced by the presence of other chemical compounds, resulting in discrepancies [11]. In addition, this requires extensive sample preparation and consequently results in sample disintegration. By contrast, although NMR has low sensitivity, it provides high reproducibility. The sample preparation is straightforward and non-destructive, enabling further re-analysis of the same sample [11,12,13]. Using these approaches, cancer-related alterations of metabolic pathway as well as molecular markers could be discovered. Bansal et al. explored the metabolic differences between the serum of BCa patients and that of healthy controls using NMR, and they proposed that the combination of dimethylamine, malonate, lactate, glutamine, histidine, and valine could differentiate BCa from the control [14]. Ming Cao et al. showed that the serum metabolic profile of BCa patients was distinctly different from that of the calculi and healthy individuals using NMR spectroscopy. They reported a decrease in the levels of isoleucine/leucine, tyrosine, lactate, glycine, and citrate as well as an increase in those of lipids and glucose. They offered potentially valuable information for the pathogenic process of BCa and a noninvasive method for the diagnosis of the disease [15]. In addition, sarcosine (N-methylglycine) levels, as a well-known biomarker of prostate cancer, was elevated during prostate cancer progression compared to very low levels in a healthy control [16]. In an NMR-based metabolomics study, several metabolites including citrate, myo-inositol, and spermine was identified as tentative biomarkers of prostate cancer [17,18]. Furthermore, in metabolomics study for renal cell carcinoma, Monteiro et al. defined 32 metabolite markers that clearly discriminate RCC patients from healthy controls by principal component analysis [19].
From these many previous studies, it has been proven that NMR-based metabolomics is a superior methodology for clinical diagnosis of various disease. On the other hand, these metabolic studies in urological oncology have been focused on biomarker discovery by comparing the control group of healthy subjects, whereas the studies about metabolic comparison among bladder, prostate cancer, and renal cell carcinoma have not yet been described. In this study, using NMR spectroscopy, we mainly aim to compare the metabolic differences between BCa, PCa, and RCC and to find the discriminant metabolites among cancer groups. Therefore, the purpose of our study is not to determine whether or not cancer is cancerous compared to a healthy control, but to differentiate between urological cancer groups in metabolomics. Furthermore, through this preliminary study, we expect to provide a faster determination of urological cancer types in patients with voiding dysfunction.

2. Results and Discussion

2.1. Study Population Characteristics

In the present study, a total of 65 patients were classified into the following groups: 29 patients with BCa, 24 patients with PCa, and 12 patients with RCC. The mean age of the patients with BCa was 67.92 ± 10.42 years (range, 58–79 years), and those of the patients with PCa and RCC were 68.21 ± 5.54 years (range, 59–79 years) and 62.17 ± 11.76 years (range, 34–78 years), respectively. All subjects in BCa, PCa, and RCC were matched for sex. On the other hand, serum creatinine (SCr) and glomerular filtration rate (eGFR) of BCa, PCa, and RCC were observed to be less than 1.2 mg/dL and more than 60 mL/min/1.73 m2, respectively, which were normal. These results mean that the renal function of the patients in each group is normal [20]. In the BCa group, 15/29 (48.27%) were classified as low grade by WHO/ISUP classification. The majority of the patients with BCa were found to be stage Ta (44.83%) and T1 (37.93%) at the final pathology, followed by T2 (15.15%) and Tis (6.89%). Among the 24 patients with PCa, 5 (20.9%), 3 (12.5%), 8 (33.3%), and 8 (33.3%) had clinical Gleason scores of 6, 7, 8, and above 8, respectively, upon preoperative biopsy. A total of 8 and 16 urine samples were collected for hormone-sensitive prostate cancer (HSPC) and castration-resistant prostate cancer (CRPC), respectively. The patients with PCa were found to be stage T1 (12.5%), T2 (25.0%), T3 (58.33%), and T4 (4.17%) at the final pathology. In the RCC group, 8/12 (66.67%) patients were classified as low-grade (Fuhrman grade 2). Furthermore, 7/12 (58.33%) patients were pT1 stage, followed by T3 (4/12, 33.33%), and T4 (1/12, 8.33%) stages. The detailed clinical characteristics of the patients are presented in Table 1.

2.2. Metabolic Differences for Urological Cancers by Using Multivariate Pattern Recognition Analysis

Representative 1H NMR spectra with water suppression obtained from the urine of patients with PCa, BCa, and RCC are shown in Figure 1.
A total of 40 metabolites were commonly identified in urinary cancer groups and quantified using Chenomx suite 8.4 NMR 600 MHz library based on human metabolome database (HMDB) database. In addition, to accurately assign the resonances of metabolites, literature data and public databases, such as Pubchem were used [21]. To explore the possible differences in the metabolic profiles, a multivariate statistical analysis for the three groups was performed by applying supervised method, partial least squares discriminant analysis (PLS-DA) and orthogonal partial least squares discriminant analysis (OPLS-DA) score plot, which require prior knowledge of sample clustering for the purpose of discriminant analysis and are used to elucidate the most reliable variables for group separation [22,23].
The PLS-DA model which was built using the two principal components showed the overall discriminatory pattern among BCa, PCa, and RCC. This model was assessed by the goodness of fit R2Y = 0.547 and predictive quality Q2 = 0.281 (Figure 2A). For improved interpretability, we preferred the OPLS-DA model which filters out the structured noise in the data set, uncorrelated to the major discriminating response, thus reducing the complexity of the model [23]. As shown in Figure 2B, we observed a clear separation among the three urological groups compared with the PLS-DA model. This OPLS-DA showed a proper fitting of data (R2Y value = 0.846) and predictability (Q2 value = 0.328). The significance of these classification models was assessed by cross-validated predictive residuals p(CV-ANOVA) (estimated by p-value < 0.05). Furthermore, the differences of metabolic pattern between intergroups were identified. The OPLS-DA models of BCa vs. RCC (Figure S1A) and PCa vs. RCC (Figure S1B) showed good predictability, whereas the model of PCa vs. BCa showed lower predictability compared to other classification models, which indicated a low metabolic difference between PCa and BCa (Figure S1C).
To select significant metabolites that mostly contributed to the differences among the three groups, the variable selection method (VIP) was applied to the OPLS-DA model based on t [2] axis. Generally, a cutoff for VIP higher than 1 is regarded as proper [24], therefore 22 metabolites with VIP > 1 were selected as influential metabolites in separation of groups (Figure 2C). Particularly, creatinine showed the highest value, followed by dimethylamine, taurine, mannitol, myo-inositol, glucose, citrate, and N-methylhydantoin. Studies based on NMR spectroscopy for the detection of metabolite markers of prostate cancer have mainly been performed on prostate biopsy tissue or prostatic secretion. From these studies, myo-inositol, citrate, and choline-containing compounds were reported as potential markers of PCa [18,25,26]. In addition, other urinary metabolomics studies of BCa suggested that carnitine, pyroglutamate, pseudouridine, histidine, and hippurate could be putative biomarkers [27,28,29,30,31]. Similarly, in RCC, Monteiro et al. revealed 32 urinary signature metabolites containing 4-hydroxyphenylacetate and hippurate to be excreted at lower levels in RCC [19]. Furthermore, Kim et al. showed that quinolinate, 4-hydroxybenzoate, and genistate were significantly differentially expressed in RCC patients compared with the control group [32]. Taking into consideration the accumulating evidence regarding the urinary metabolomics of urological cancers, we hypothesized that the 22 metabolites with VIP values higher than 1 identified in this study could be considered as tentative discriminatory variables of BCa, PCa, and RCC.

2.3. Comparison of Differential Metabolites Levels among Urological Cancers Using Univariate Statistical Analysis

Quantified metabolites were normalized by total area to minimize the dilution effect of urine among groups. The median value and interquartile range (IQR) of the normalized metabolite concentration for BCa, PCa, and RCC are presented in Table S1.
Before identifying the metabolites, which distinguish the three urological cancers, the normal distribution of the urinary metabolic profile dataset was examined using the Kolmogorov-Smirnov test, which showed that the samples under analysis did not exhibit normal distribution (p < 0.05) within each group (Table S1). Then, an age-adjusted non-parametric analysis of covariance (Ranked ANCOVA) was used to determine significant metabolites among the groups with a Bonferroni correction method (adjusted p-value < 0.05) (Table S1) [33]. A total of five metabolites were differentially expressed in urological cancer groups: 4-hydroxybenzoate (4-HBA), N-methylhydantoin (N-MH), creatinine, glutamine, and acetate (Figure 3A–E). Urinary creatinine levels among these metabolites may be affected by renal function, but as described above, SCr and eGFR levels between cancer groups were confirmed to be normal, and this means that there were no metabolic differences according to renal function among groups.
4-HBA is a phenolic derivative of benzoate with anti-cancer properties [34]. Seidel et al. found that 4-HBA, a histone deacetylase-specific inhibitor, reduced the proliferation and viability of prostate cancer cells, and decreased androgen receptor levels, which play a role in the development and progression of prostate cancer [35]. Recently, in plasma metabolomic profiling by GC/MS, 4-HBA was found to discriminate between the control, prostatic intraepithelial neoplasia, and prostate cancer, and its levels increased with disease progression. In our study, 4-HBA had the highest concentration in PCa (Figure 3A), and based on previous studies, we explain that the levels of 4-HBA could be useful in determining the abnormal prostate metabolic profile [36].
N-MH is derived from creatinine-by-creatinine deaminase through an alternative creatinine-creatine pathway [37,38,39,40]. N-MH is further degraded to N-carbamoylsarcosine and sarcosine, and consequently to glycine. Several studies have explored the relation between the urinary levels of sarcosine and prostate cancer, which indicated increased sarcosine levels as a potential marker of prostate cancer [16,41,42,43]. Although other intermediates in the pathway for creatinine degradation to N-MH (e.g., N-carbamoylsarcosine, sarcosine) could not be detected in this study, we confirmed the highest median levels of N-methylhydantoin and glycine and lower levels of creatinine in PCa, suggesting that the metabolic pathway for creatinine degradation to N-MH could be a specific active pathway in prostate cancer cells (Figure 3B,C and Table S1).
Glutamine is converted to glutamate and pyroglutamate, which is produced by the cycling of glutamine involved in GSH metabolism [44,45]. Whenever there is an increased demand for GSH to defend against stress responses, this cyclic form of pyroglutamate is opened by 5-oxoprolinase to regenerate glutamate, which then produces GSH [46]. Ahmad et al. found an increased rate of synthesis of GSH from glutamine in papillary renal cell carcinoma tissues by proteomic and metabolomic analysis [47]. Similarly, we identified significantly lower levels of glutamine in RCC (Figure 3D and Table S1), suggesting that GSH synthesis from glutamine is promoted in renal cell carcinoma.
As a metabolite relevant to major energy metabolism pathways, acetate, which is produced by glucose-derived pyruvate, supports acetyl-CoA metabolism and plays an important role in bioenergetics [48,49]. In tumor cells, acetate is mainly used as fuel for lipid synthesis by formed to acetyl-CoA [50]. Gao et al. reported that acetate, which is a crucial product of fatty acid β-oxidation, contributed to the discrimination of paired adjacent tissue from advanced RCC [51]. Furthermore, many previous studies have investigated an association between renal cell carcinoma and lipid metabolism, revealing that the lipid accumulation in tumor cell is the histological and morphological characteristics of RCC [52,53,54]. There is also evidence to suggest that lipogenesis is a feature of the malignant phase of RCC [55]. Based on the previous studies, we suggest that significantly high acetate levels in RCC could be supported to previous studies related to lipid metabolism in RCC (Figure 3E).

2.4. Evaluation of the Predictive Effectiveness of Significant Metabolites

To confirm the discriminative accuracy of individual metabolites between groups, we performed the univariate receiver operating characteristic (ROC) curve analysis using the prominent five metabolites (Figure 3). The criteria for assessing the accuracy of the signature based on the area under the curve (AUC), summarized into a single metric of ROC curve, were as follows: 0.9 to 1.0 = excellent; 0.8 to 0.9 = good to very good; 0.7 to 0.8 = fair; 0.6 to 0.7 = poor; 0.5 to 0.6 = fail [56]. Therefore, we found significant metabolites with AUC > 0.7 (p < 0.05) in each univariate ROC curve analyses (Table 2), and the multivariate ROC curves for intergroup were created by combining the significant metabolites in univariate ROC analysis and assessed for their predictability.
The multivariate ROC curve of BCa vs. RCC combined with creatinine and glutamine showed an AUC value of 0.778 with sensitivity = 69.7% and specificity = 75% (Figure 4A). The multivariate ROC curve with 4-HBA, N-MH, and creatinine showed with an AUC 0.898, sensitivity = 93.1%, specificity = 79.2% in classification between PCa and BCa (Figure 4B). The composition of 4-HBA, N-MH, glutamine, and acetate was able to accurately differentiate 95.1% of PCa from RCC, with high sensitivity (83.3%) and specificity (83.3%) (Figure 4C). The all-multivariate ROC curves between BCa, PCa and RCC showed high performance with empirical p-values less than 0.05 (Figure S2A–C), suggesting that these multivariate ROC findings can be used to determine which cancer is under BCa, PCa, and RCC.

3. Materials and Methods

3.1. Reagents

Deuterium oxide (D2O, 99% in D) was purchased from Cambridge Isotope Laboratories, Inc., (Andover, MA, USA). Potassium phosphate dibasic (K2HPO4, ACS grade reagent), sodium phosphate monobasic (anhydrous), hydrochloric acid (HCl, ACS grade reagent), sodium hydroxide solution (NaOH, HPLC grade), and 3-trimethylsilyl propionic-2,2,3,3-d4 acid sodium salt (TSP-d4) were obtained from Sigma-Aldrich (St. Louis, MO, USA).

3.2. Sample Collection

A total of 65 patients were recruited from Pusan National University Hospital (PNUH), including 29 with bladder cancer, 24 with prostate cancer, and 12 with renal cancer. The study was approved by the Institutional Review Board (IRB) of PNUH (IRB No 1802-004-063). Urine samples were collected during outpatient follow-up or prior to urological surgery. At the outpatient department, patients were asked to take mid-stream urine in an aseptic container. In the operating room, urine samples were acquired in an aseptic container via catheterization using 7 Fr. Nelaton catheter under general anesthesia. Before catheterization, preoperative skin preparation was conducted on the external genital organs using povidone-iodine solution. All the samples were over 50 mL in volume and immediately stored in a −80 ℃ freezer.

3.3. Preparation of Urine Sample for NMR Analysis

Urine samples for NMR analysis were prepared according to criterion procedures [57,58]. Samples were thawed at 25 ℃ (Room temperature), shaken before use, and 700 μL of each urine aliquot were centrifuged at 12,000× g for 15 min at 4 ℃ to precipitate protein and other macromolecules [45]. To minimize chemical shift variation described by Xiao et al. [59], 70 μL of K2HPO4 and NaH2PO4 buffer (pH 7.4, 1.5 M in D2O) containing 20 mM TSP-d4 was added to 630 μL of urine sample supernatant. Then, 700 μL of the mixture were transferred into a 5 mm NMR tube for NMR measurement.

3.4. NMR Spectra Acquisition and Spectral Processing

NMR measurements were conducted using an Agilent NMR spectrometer (Agilent, Santa Clara, CA, USA) equipped with a 600 MHz magnet at 298 K. All 1H NMR data collection and pre-processing were performed using VnmrJ Software (Agilent Technologies, Santa Clara, CA, USA). All NMR data of each sample were measured by a pulse sequence PRESAT in all cases to suppress the water signal during a relaxation time of 2 s and 90-degree pulse (pw) of 11.4 μs at 128 scans into 32 K data points. NMR data were obtained in the form of FID files, and the spectrum was acquired by Fourier transformation. The baseline and phase of the spectrum were manually corrected, and all spectra were calibrated with reference to the TSP-d4 signal at δ 0.00 ppm. The NMR spectrum segmented into buckets of 0.005 ppm with a range of 0.5 to 8.5 ppm was analyzed after excluding the water and urea peak region (4.2 to 6.1 ppm) and normalized to total area to minimize the urinary dilution effect among groups, then the peaks were aligned. NMR data were processed using Mnova (Mestrelab Research, Santiago de Compostela, Spain) because of the software suite’s useful functions such as peak binning and alignment, spectral normalization, and spectral superimposing. The peak signals of metabolites were identified and quantified using 600 MHz library from Chenomx NMR suite 8.4 professional (Chenomx Inc., Edmonton, AB, Canada), which could be obtained by comparing a known reference peak signal (TSP) with signals from a library of metabolites. In addition, to accurately confirm the identified compound signal, we also double-checked using data from public databases such as the human metabolome database (HMDB, www.hmdb.ca, accessed on 10 June 2021) [60] and Pubchem, and the literature [21]. The relative concentration of metabolites was analyzed to compare the significant differences among groups.

3.5. Statistical Analysis

To confirm the differences in the metabolic profile among urological cancers, multivariate pattern recognition analysis was performed for all processed NMR spectrum data using the SIMCA-P 12.0.1 software (Umetrics, Umeå, Sweden). The PLS/OPLS-DA analysis was performed with pareto scaling to examine the possible minor metabolic differences among groups, and the VIP analysis was used to identify the metabolites contributing for group separation in OPLS-DA model. The quality of the model was assessed using the cross-validation parameters R2Y(cum), Q2 (cum) which is reported that Q2 > 0.5 is admitted for good predictability in the model [61]. Furtherly, to assess the reliability of these classification models, analysis of variance testing of p(CV-ANOVA) was performed and p-value under 0.05 was determined to be significant. Before examining the significantly differential levels of metabolites among groups, the distribution of normalized concentration was estimated by the Kolmogorov-Smirnov test, and the metabolites with abnormal distribution were conducted by a non-parametric analysis of covariance (Ranked-ANCOVA) by adjusting the covariate (age) as a potential confounding factor in metabolites data set. The post-hoc tests of a ranked ANCOVA were executed using with a Bonferroni correction method by adjusted p-value (p < 0.05). Statistical analyses for the quantities of metabolites were performed using IBM SPSS Statistics 25 (SPSS, Inc., Chicago, IL, USA). ROC analysis between intergroups was performed using a web server for metabolomics data analysis (MetaboAnalyst 5.0, https://www.metaboanalyst.ca, accessed on 12 July 2021) [62]. Firstly, univariate ROC analysis was examined, and AUC and p-values of each ROC curves were used to evaluate the predictability. Then to improve the discriminatory accuracy, multivariate ROC curves were plotted with false positive rate and true positive rate using combination of significant metabolites with AUC > 0.7 (p-value < 0.05).

4. Conclusions

In summary: the present study revealed the different urinary metabolic profiles of bladder cancer (BCa), prostate cancer (PCa), and renal cell carcinoma (RCC) based on the 1H NMR spectroscopy with multivariate pattern recognition analysis and quantitative analysis. For seeking the prominent metabolites, several criteria were applied: (i) the metabolites with VIP value > 1.0 in the OPLS-DA model; (ii) ranked-ANCOVA with a Bonferroni-adjusted p-value < 0.05. Consequently, a total of five metabolites were identified as significant. These metabolites were used in univariate and multivariate ROC curve analysis to examine their sensitivity and specificity in discriminating among cancer groups. Multivariate ROC curves created by combining each significant metabolites for BCa vs. RCC (including of creatinine, and glutamine), PCa vs. BCa (including 4-hydroxybenzoate, N-methylhydantoin, and creatinine), and RCC vs. PCa (including 4-hydroxybenzoate, N-methylhydantoin, glutamine, and acetate) showed better discriminatory accuracy than univariate ROC model, suggesting that these multiple metabolites panel could be determinant of the urological cancers. Although this preliminary study demonstrated a comparison of holistic metabolic profiles between the urological cancers and supports the explorative role of metabolomics in the investigation of human urine, it is limited in that the sample size is small. Therefore, further experiments should be performed with more subjects with matched clinical characteristics (e.g., stage of cancer, age, etc.), including healthy controls. Through further study, we expect to provide a faster determination of urological malignancy types in patients with voiding dysfunction.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/metabo11090591/s1, Figure S1: Orthogonal partial least-squares discriminant analysis (OPLS-DA) score plot between inter cancer groups, Figure S2: Empirical p-value plots of multivariate ROC curve analysis. Table S1: Identified and quantified urinary metabolites from 1H NMR spectra.

Author Contributions

All authors have made substantial contribution in this work. Conceptualization, H.K.H. and S.K.; Data curation, J.Y.K. and B.J.K.; Formal analysis, S.L.; Investigation, S.L.; Project administration, S.K.; Resources J.Y.K., B.J.K. and H.K.H.; Supervision, S.K.; Writing—original draft, S.L. and K.H.K.; Writing—review & editing, K.H.K. and H.K.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board (IRB) of PNUH (IRB No 1802-004-063).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data are contained within the article or supplementary.

Acknowledgments

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. NRF-2020R1I1A2075016).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2020. CA Cancer J. Clin. 2020, 70, 7–30. [Google Scholar] [CrossRef]
  2. Bruyninckx, R.; Buntinx, F.; Aertgeerts, B.; Van Casteren, V. The diagnostic value of macroscopic haematuria for the diagnosis of urological cancer in general practice. Br. J. Gen. Pract. 2003, 53, 31–35. [Google Scholar]
  3. Gadler, T.; Keedy, M.; Rivas, N. A case of hematuria. Adv. Emerg. Nurs. J. 2010, 32, 30–41. [Google Scholar] [CrossRef]
  4. Urquidi, V.; J Rosser, C.; Goodison, S. Molecular diagnostic trends in urological cancer: Biomarkers for non-invasive diagnosis. Curr. Med. Chem. 2012, 19, 3653–3663. [Google Scholar] [CrossRef] [Green Version]
  5. Badalament, R.A.; Hermansen, D.K.; Kimmel, M.; Gay, H.; Herr, H.W.; Fair, W.R.; Whitmore, W.F., Jr.; Melamed, M.R. The sensitivity of bladder wash flow cytometry, bladder wash cytology, and voided cytology in the detection of bladder carcinoma. Cancer 1987, 60, 1423–1427. [Google Scholar] [CrossRef]
  6. Catalona, W.J.; Richie, J.P.; Ahmann, F.R.; Hudson, M.L.A.; Scardino, P.T.; Flanigan, R.C.; Dekernion, J.B.; Ratliff, T.L.; Kavoussi, L.R.; Dalkin, B.L.; et al. Comparison of Digital Rectal Examination and Serum Prostate Specific Antigen in the Early Detection of Prostate Cancer: Results of a Multicenter Clinical Trial of 6,630 Men. J. Urol. 1994, 151, 1283–1290. [Google Scholar] [CrossRef]
  7. Beckonert, O.; Keun, H.C.; Ebbels, T.M.; Bundy, J.; Holmes, E.; Lindon, J.C.; Nicholson, J.K. Metabolic profiling, metabolomic and metabonomic procedures for NMR spectroscopy of urine, plasma, serum and tissue extracts. Nat. Protoc. 2007, 2, 2692–2703. [Google Scholar] [CrossRef]
  8. Capati, A.; Ijare, O.B.; Bezabeh, T. Diagnostic applications of nuclear magnetic resonance–based urinary metabolomics. Magn. Reson. Insights 2017, 10, 1178623X17694346. [Google Scholar] [CrossRef] [Green Version]
  9. Guijas, C.; Montenegro-Burke, J.R.; Warth, B.; Spilker, M.E.; Siuzdak, G. Metabolomics activity screening for identifying metabolites that modulate phenotype. Nat. Biotechnol. 2018, 36, 316–320. [Google Scholar] [CrossRef]
  10. Liang, Q.; Liu, H.; Zhang, T.; Jiang, Y.; Xing, H.; Zhang, A.-h. Potential urine biomarkers from a high throughput metabolomics study of severe sepsis in a large Asian cohort. RSC Adv. 2015, 5, 102204–102209. [Google Scholar] [CrossRef]
  11. Dunn, W.B.; Bailey, N.J.; Johnson, H.E. Measuring the metabolome: Current analytical technologies. Analyst 2005, 130, 606–625. [Google Scholar] [CrossRef]
  12. Psihogios, N.G.; Gazi, I.F.; Elisaf, M.S.; Seferiadis, K.I.; Bairaktari, E.T. Gender-related and age-related urinalysis of healthy subjects by NMR-based metabonomics. NMR Biomed. 2008, 21, 195–207. [Google Scholar] [CrossRef]
  13. Lehman-Mckeeman, L.D.; Car, B.D. Metabonomics: Application in Predictive and Mechanistic Toxicology. Toxicol. Pathol. 2016, 32, 94–95. [Google Scholar] [CrossRef] [Green Version]
  14. Bansal, N.; Gupta, A.; Mitash, N.; Shakya, P.S.; Mandhani, A.; Mahdi, A.A.; Sankhwar, S.N.; Mandal, S.K. Low- and high-grade bladder cancer determination via human serum-based metabolomics approach. J. Proteome Res. 2013, 12, 5839–5850. [Google Scholar] [CrossRef]
  15. Cao, M.; Zhao, L.; Chen, H.; Xue, W.; Lin, D. NMR-based metabolomic analysis of human bladder cancer. Anal. Sci. 2012, 28, 451–456. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Sreekumar, A.; Poisson, L.M.; Rajendiran, T.M.; Khan, A.P.; Cao, Q.; Yu, J.; Laxman, B.; Mehra, R.; Lonigro, R.J.; Li, Y.; et al. Metabolomic profiles delineate potential role for sarcosine in prostate cancer progression. Nature 2009, 457, 910–914. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Kurhanewicz, J.; Vigneron, D.B.; Nelson, S.J.; Hricak, H.; MacDonald, J.M.; Konety, B.; Narayan, P. Citrate as an in vivo marker to discriminate prostate cancer from benign prostatic hyperplasia and normal prostate peripheral zone: Detection via localized proton spectroscopy. Urology 1995, 45, 459–466. [Google Scholar] [CrossRef]
  18. Lynch, M.J.; Nicholson, J.K. Proton MRS of human prostatic fluid: Correlations between citrate, spermine, and myo-inositol levels and changes with disease. Prostate 1997, 30, 248–255. [Google Scholar] [CrossRef]
  19. Monteiro, M.S.; Barros, A.S.; Pinto, J.; Carvalho, M.; Pires-Luis, A.S.; Henrique, R.; Jeronimo, C.; Bastos, M.L.; Gil, A.M.; Guedes de Pinho, P. Nuclear Magnetic Resonance metabolomics reveals an excretory metabolic signature of renal cell carcinoma. Sci. Rep. 2016, 6, 37275. [Google Scholar] [CrossRef] [Green Version]
  20. Edwards, K.D.; Whyte, H.M. Plasma creatinine level and creatinine clearance as tests of renal function. Australas Ann. Med. 1959, 8, 218–224. [Google Scholar] [CrossRef]
  21. Bouatra, S.; Aziat, F.; Mandal, R.; Guo, A.C.; Wilson, M.R.; Knox, C.; Bjorndahl, T.C.; Krishnamurthy, R.; Saleem, F.; Liu, P.; et al. The human urine metabolome. PLoS ONE 2013, 8, e73076. [Google Scholar] [CrossRef] [Green Version]
  22. Lever, J.; Krzywinski, M.; Altman, N. Points of significance: Principal component analysis. Nat. Methods 2017, 14, 641–643. [Google Scholar] [CrossRef] [Green Version]
  23. Bylesjö, M.; Rantalainen, M.; Cloarec, O.; Nicholson, J.K.; Holmes, E.; Trygg, J. OPLS discriminant analysis: Combining the strengths of PLS-DA and SIMCA classification. J. Chemom. 2006, 20, 341–351. [Google Scholar] [CrossRef]
  24. Chong, I.-G.; Jun, C.-H. Performance of some variable selection methods when multicollinearity is present. Chemom. Intell. Lab. Syst. 2005, 78, 103–112. [Google Scholar] [CrossRef]
  25. Swanson, M.G.; Keshari, K.R.; Tabatabai, Z.L.; Simko, J.P.; Shinohara, K.; Carroll, P.R.; Zektzer, A.S.; Kurhanewicz, J. Quantification of choline- and ethanolamine-containing metabolites in human prostate tissues using 1H HR-MAS total correlation spectroscopy. Magn. Reson Med. 2008, 60, 33–40. [Google Scholar] [CrossRef] [Green Version]
  26. Serkova, N.J.; Gamito, E.J.; Jones, R.H.; O’Donnell, C.; Brown, J.L.; Green, S.; Sullivan, H.; Hedlund, T.; Crawford, E.D. The metabolites citrate, myo-inositol, and spermine are potential age-independent markers of prostate cancer in human expressed prostatic secretions. Prostate 2008, 68, 620–628. [Google Scholar] [CrossRef]
  27. Zhang, W.-T.; Zhang, Z.-W.; Guo, Y.-D.; Wang, L.-S.; Mao, S.-Y.; Zhang, J.-F.; Liu, M.-N.; Yao, X.-D. Discovering biomarkers in bladder cancer by metabolomics. Biomark. Med. 2018, 12, 1347–1359. [Google Scholar] [CrossRef]
  28. Kim, W.T.; Yun, S.J.; Yan, C.; Jeong, P.; Kim, Y.H.; Lee, I.S.; Kang, H.W.; Park, S.; Moon, S.K.; Choi, Y.H.; et al. Metabolic Pathway Signatures Associated with Urinary Metabolite Biomarkers Differentiate Bladder Cancer Patients from Healthy Controls. Yonsei Med. J. 2016, 57, 865–871. [Google Scholar] [CrossRef]
  29. Pasikanti, K.K.; Esuvaranathan, K.; Hong, Y.; Ho, P.C.; Mahendran, R.; Raman Nee Mani, L.; Chiong, E.; Chan, E.C. Urinary metabotyping of bladder cancer using two-dimensional gas chromatography time-of-flight mass spectrometry. J. Proteome Res. 2013, 12, 3865–3873. [Google Scholar] [CrossRef]
  30. Huang, Z.; Lin, L.; Gao, Y.; Chen, Y.; Yan, X.; Xing, J.; Hang, W. Bladder cancer determination via two urinary metabolites: A biomarker pattern approach. Mol. Cell Proteom. 2011, 10, M111.007922. [Google Scholar] [CrossRef] [Green Version]
  31. Alberice, J.V.; Amaral, A.F.; Armitage, E.G.; Lorente, J.A.; Algaba, F.; Carrilho, E.; Marquez, M.; Garcia, A.; Malats, N.; Barbas, C. Searching for urine biomarkers of bladder cancer recurrence using a liquid chromatography-mass spectrometry and capillary electrophoresis-mass spectrometry metabolomics approach. J. Chromatogr. A 2013, 1318, 163–170. [Google Scholar] [CrossRef]
  32. Kim, K.; Taylor, S.L.; Ganti, S.; Guo, L.; Osier, M.V.; Weiss, R.H. Urine metabolomic analysis identifies potential biomarkers and pathogenic pathways in kidney cancer. Omics J. Integr. Biol. 2011, 15, 293–303. [Google Scholar] [CrossRef] [Green Version]
  33. Cangür, Ş.; Sungur, M.A.; Ankarali, H. The Methods Used in Nonparametric Covariance Analysis. Düzce Tıp Fakültesi Dergrisi 2018, 20, 1–6. [Google Scholar] [CrossRef]
  34. Manuja, R.; Sachdeva, S.; Jain, A.; Chaudhary, J. A comprehensive review on biological activities of p-hydroxy benzoic acid and its derivatives. Int. J. Pharm. Sci. Rev. Res. 2013, 22, 109–115. [Google Scholar]
  35. Seidel, C.; Schnekenburger, M.; Mazumder, A.; Teiten, M.H.; Kirsch, G.; Dicato, M.; Diederich, M. 4-Hydroxybenzoic acid derivatives as HDAC6-specific inhibitors modulating microtubular structure and HSP90alpha chaperone activity against prostate cancer. Biochem. Pharm. 2016, 99, 31–52. [Google Scholar] [CrossRef]
  36. Markin, P.A.; Brito, A.; Moskaleva, N.; Lartsova, E.V.; Shpot, Y.V.; Lerner, Y.V.; Mikhajlov, V.Y.; Potoldykova, N.V.; Enikeev, D.V.; La Frano, M.R.; et al. Plasma metabolomic profile in prostatic intraepithelial neoplasia and prostate cancer and associations with the prostate-specific antigen and the Gleason score. Metabolomics 2020, 16, 74. [Google Scholar] [CrossRef]
  37. Wyss, M.; Kaddurah-Daouk, R. Creatine and creatinine metabolism. Physiol. Rev. 2000, 80, 1107–1213. [Google Scholar] [CrossRef]
  38. Szulmajster, J. Bacterial Fermentation of Creatinine I.: Isolation of N-Methyl-Hydantoin. J. Bacteriol. 1958, 75, 633–639. [Google Scholar] [CrossRef] [Green Version]
  39. Yamada, H.; Shimizu, S.; Kim, J.M.; Shinmen, Y.; Sakai, T. A novel metabolic pathway for creatinine degradation in Pseudomonas putida 77. FEMS Microbiol. Lett. 1985, 30, 337–340. [Google Scholar] [CrossRef]
  40. Shimizu, S.; Kim, J.M.; Shinmen, Y.; Yamada, H. Evaluation of two alternative metabolic pathways for creatinine degradation in microorganisms. Arch. Microbiol 1986, 145, 322–328. [Google Scholar] [CrossRef]
  41. Gkotsos, G.; Virgiliou, C.; Lagoudaki, I.; Sardeli, C.; Raikos, N.; Theodoridis, G.; Dimitriadis, G. The Role of Sarcosine, Uracil, and Kynurenic Acid Metabolism in Urine for Diagnosis and Progression Monitoring of Prostate Cancer. Metabolites 2017, 7, 9. [Google Scholar] [CrossRef] [Green Version]
  42. Cernei, N.; Heger, Z.; Gumulec, J.; Zitka, O.; Masarik, M.; Babula, P.; Eckschlager, T.; Stiborova, M.; Kizek, R.; Adam, V. Sarcosine as a potential prostate cancer biomarker--a review. Int. J. Mol. Sci. 2013, 14, 13893–13908. [Google Scholar] [CrossRef] [Green Version]
  43. Cavaliere, B.; Macchione, B.; Monteleone, M.; Naccarato, A.; Sindona, G.; Tagarelli, A. Sarcosine as a marker in prostate cancer progression: A rapid and simple method for its quantification in human urine by solid-phase microextraction-gas chromatography-triple quadrupole mass spectrometry. Anal. Bioanal. Chem. 2011, 400, 2903–2912. [Google Scholar] [CrossRef]
  44. Welbourne, T.C. Ammonia production and glutamine incorporation into glutathione in the functioning rat kidney. Can. J. Biochem. 1979, 57, 5. [Google Scholar] [CrossRef] [PubMed]
  45. Yeh, S.-L.; Shih, Y.-M.; Lin, M.-T. Glutamine and its antioxidative potentials in diabetes. In Diabetes; Academic Press: Cambridge, MA, USA, 2020; pp. 255–264. [Google Scholar]
  46. Lord, R.S.; Bralley, J.A. Clinical applications of urinary organic acids. Part I: Detoxification markers. Altern Med. Rev. 2008, 13, 205–215. [Google Scholar]
  47. Al Ahmad, A.; Paffrath, V.; Clima, R.; Busch, J.F.; Rabien, A.; Kilic, E.; Villegas, S.; Timmermann, B.; Attimonelli, M.; Jung, K.; et al. Papillary Renal Cell Carcinomas Rewire Glutathione Metabolism and Are Deficient in Both Anabolic Glucose Synthesis and Oxidative Phosphorylation. Cancers 2019, 11, 1298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Comerford, S.A.; Huang, Z.; Du, X.; Wang, Y.; Cai, L.; Witkiewicz, A.K.; Walters, H.; Tantawy, M.N.; Fu, A.; Manning, H.C.; et al. Acetate dependence of tumors. Cell 2014, 159, 1591–1602. [Google Scholar] [CrossRef] [Green Version]
  49. Liu, X.; Cooper, D.E.; Cluntun, A.A.; Warmoes, M.O.; Zhao, S.; Reid, M.A.; Liu, J.; Lund, P.J.; Lopes, M.; Garcia, B.A.; et al. Acetate Production from Glucose and Coupling to Mitochondrial Metabolism in Mammals. Cell 2018, 175, 502–513. [Google Scholar] [CrossRef] [Green Version]
  50. Schug, Z.T.; Vande Voorde, J.; Gottlieb, E. The metabolic fate of acetate in cancer. Nat. Rev. Cancer 2016, 16, 708–717. [Google Scholar] [CrossRef]
  51. Gao, H.; Dong, B.; Jia, J.; Zhu, H.; Diao, C.; Yan, Z.; Huang, Y.; Li, X. Application of ex vivo (1)H NMR metabonomics to the characterization and possible detection of renal cell carcinoma metastases. J. Cancer Res. Clin. Oncol. 2012, 138, 753–761. [Google Scholar] [CrossRef]
  52. Pinthus, J.H.; Whelan, K.F.; Gallino, D.; Lu, J.P.; Rothschild, N. Metabolic features of clear-cell renal cell carcinoma: Mechanisms and clinical implications. Can. Urol. Assoc. J. 2011, 5, 274–282. [Google Scholar] [CrossRef]
  53. Rezende, R.B.; Drachenberg, C.B.; Kumar, D.; Blanchaert, R.; Ord, R.A.; Ioffe, O.B.; Papadimitriou, J.C. Differential diagnosis between monomorphic clear cell adenocarcinoma of salivary glands and renal (clear) cell carcinoma. Am. J. Surg. Pathol. 1999, 23, 1532–1538. [Google Scholar] [CrossRef]
  54. Sanchez, D.J.; Simon, M.C. Genetic and metabolic hallmarks of clear cell renal cell carcinoma. Biochim. Biophys. Acta Rev. Cancer 2018, 1870, 23–31. [Google Scholar] [CrossRef]
  55. Tun, H.W.; Marlow, L.A.; von Roemeling, C.A.; Cooper, S.J.; Kreinest, P.; Wu, K.; Luxon, B.A.; Sinha, M.; Anastasiadis, P.Z.; Copland, J.A. Pathway signature and cellular differentiation in clear cell renal cell carcinoma. PLoS ONE 2010, 5, e10696. [Google Scholar] [CrossRef]
  56. Emwas, A.H.; Roy, R.; McKay, R.T.; Ryan, D.; Brennan, L.; Tenori, L.; Luchinat, C.; Gao, X.; Zeri, A.C.; Gowda, G.A.; et al. Recommendations and Standardization of Biomarker Quantification Using NMR-Based Metabolomics with Particular Focus on Urinary Analysis. J. Proteome Res. 2016, 15, 360–373. [Google Scholar] [CrossRef] [Green Version]
  57. Vignoli, A.; Ghini, V.; Meoni, G.; Licari, C.; Takis, P.G.; Tenori, L.; Turano, P.; Luchinat, C. High-Throughput Metabolomics by 1D NMR. Angew Chem. Int. Ed. Engl. 2019, 58, 968–994. [Google Scholar] [CrossRef]
  58. Ghini, V.; Quaglio, D.; Luchinat, C.; Turano, P. NMR for sample quality assessment in metabolomics. N Biotechnol. 2019, 52, 25–34. [Google Scholar] [CrossRef]
  59. Xiao, C.; Hao, F.; Qin, X.; Wang, Y.; Tang, H. An optimized buffer system for NMR-based urinary metabonomics with effective pH control, chemical shift consistency and dilution minimization. Analyst 2009, 134, 916–925. [Google Scholar] [CrossRef]
  60. Wishart, D.S.; Feunang, Y.D.; Marcu, A.; Guo, A.C.; Liang, K.; Vazquez-Fresno, R.; Sajed, T.; Johnson, D.; Li, C.; Karu, N.; et al. HMDB 4.0: The human metabolome database for 2018. Nucleic Acids Res. 2018, 46, D608–D617. [Google Scholar] [CrossRef]
  61. Umetrics, A. User Guide to SIMCA-P 13.0; Umetrics Inc.: Kinnelon, NJ, USA, 2012. [Google Scholar]
  62. Chong, J.; Wishart, D.S.; Xia, J. Using MetaboAnalyst 4.0 for Comprehensive and Integrative Metabolomics Data Analysis. Curr. Protoc. Bioinform. 2019, 68, e86. [Google Scholar] [CrossRef]
Metabolites 11 00591 g001 550
Figure 1. Representative 1H NMR Spectrum with water suppression of urological cancer. (1) Valine, (2) 3-Hydroxyisobutyrate, (3) 3-Hydroxyisovalerate, (4) 2-Hydroxyisobutyrate, (5) Lactate, (6) Alanine, (7) Lysine, (8) Acetate, (9) Methionine, (10) O-Acetylcarnitine, (11) Pyroglutamate, (12) Succinate, (13) Carnitine, (14) Glutamine, (15) Citrate, (16) Dimethylamine, (17) N,N-Dimethylglycine, (18) N-Methylhydantoin, (19) Creatinine, (20) cis-Aconitate, (21) Choline, (22) Trimethylamine N-oxide, (23) Taurine, (24) 4-Hydroxyphenylacetate, (25) Glucose, (26) myo-Inositol, (27) Glycine, (28) Arabinitol, (29) Gluconate, (30) Mannitol, (31) Betaine, (32) Creatine, (33) Glycolate, (34) Hippurate, (35) Tyrosine, (36) 4-Hydroxybenzoate, (37) Histidine, (38) N-Phenylacetylglycine, (39) Pseudouridine, (40) Formate. The region of water and urea was excluded.
Figure 1. Representative 1H NMR Spectrum with water suppression of urological cancer. (1) Valine, (2) 3-Hydroxyisobutyrate, (3) 3-Hydroxyisovalerate, (4) 2-Hydroxyisobutyrate, (5) Lactate, (6) Alanine, (7) Lysine, (8) Acetate, (9) Methionine, (10) O-Acetylcarnitine, (11) Pyroglutamate, (12) Succinate, (13) Carnitine, (14) Glutamine, (15) Citrate, (16) Dimethylamine, (17) N,N-Dimethylglycine, (18) N-Methylhydantoin, (19) Creatinine, (20) cis-Aconitate, (21) Choline, (22) Trimethylamine N-oxide, (23) Taurine, (24) 4-Hydroxyphenylacetate, (25) Glucose, (26) myo-Inositol, (27) Glycine, (28) Arabinitol, (29) Gluconate, (30) Mannitol, (31) Betaine, (32) Creatine, (33) Glycolate, (34) Hippurate, (35) Tyrosine, (36) 4-Hydroxybenzoate, (37) Histidine, (38) N-Phenylacetylglycine, (39) Pseudouridine, (40) Formate. The region of water and urea was excluded.
Metabolites 11 00591 g001
Metabolites 11 00591 g002 550
Figure 2. Metabolic multivariate pattern recognition analysis. (A) PLS-DA score plots from the 1H NMR spectra with water suppression from patients with urological cancers: BCa (n = 29, ▲ red); PCa (n = 24, ◆, blue); RCC (n = 12, ●, green); (B) OPLS-DA score plots based 1H NMR spectra of urinary metabolites among cancer groups. All the score plots were assessed by R2Y, Q2, and pCV-ANOVA (p < 0.05); (C) VIP score plot represents the variables with VIP value higher than 1 selected in the OPLS-DA model. Var ID indicates each contributing metabolite.
Figure 2. Metabolic multivariate pattern recognition analysis. (A) PLS-DA score plots from the 1H NMR spectra with water suppression from patients with urological cancers: BCa (n = 29, ▲ red); PCa (n = 24, ◆, blue); RCC (n = 12, ●, green); (B) OPLS-DA score plots based 1H NMR spectra of urinary metabolites among cancer groups. All the score plots were assessed by R2Y, Q2, and pCV-ANOVA (p < 0.05); (C) VIP score plot represents the variables with VIP value higher than 1 selected in the OPLS-DA model. Var ID indicates each contributing metabolite.
Metabolites 11 00591 g002
Metabolites 11 00591 g003 550
Figure 3. Boxplots of metabolites concentrations that were significantly different among BCa, PCa, and RCC from non-parametric analysis of covariance (Ranked ANCOVA). (A) 4-hydroxybenzoate (4-HBA), (B) N-methylhydantoin (N-MH), (C) creatinine, (D) glutamine, and (E) acetate. Each dot in the boxes represents an individual sample concentration. * <0.05, ** <0.01, *** <0.001.
Figure 3. Boxplots of metabolites concentrations that were significantly different among BCa, PCa, and RCC from non-parametric analysis of covariance (Ranked ANCOVA). (A) 4-hydroxybenzoate (4-HBA), (B) N-methylhydantoin (N-MH), (C) creatinine, (D) glutamine, and (E) acetate. Each dot in the boxes represents an individual sample concentration. * <0.05, ** <0.01, *** <0.001.
Metabolites 11 00591 g003
Metabolites 11 00591 g004 550
Figure 4. Combined multivariate ROC curves between cancer groups. (A) The AUC for creatinine and glutamine to discriminate BCa from RCC was 0.778 (95% CI, 0.522–0.972). (B) The AUC for 4-HBA, N-MH, and creatinine to discriminate PCa from BCa was 0.898 (95% CI, 0.779–1.0). (C) The AUC for 4-HBA, N-MH, glutamine, and acetate to discriminate PCa from RCC was 0.951 (95% CI, 0.816–1.0).
Figure 4. Combined multivariate ROC curves between cancer groups. (A) The AUC for creatinine and glutamine to discriminate BCa from RCC was 0.778 (95% CI, 0.522–0.972). (B) The AUC for 4-HBA, N-MH, and creatinine to discriminate PCa from BCa was 0.898 (95% CI, 0.779–1.0). (C) The AUC for 4-HBA, N-MH, glutamine, and acetate to discriminate PCa from RCC was 0.951 (95% CI, 0.816–1.0).
Metabolites 11 00591 g004
Table
Table 1. Summary of the clinical characteristics of the patients with BCa, PCa and RCC.
Table 1. Summary of the clinical characteristics of the patients with BCa, PCa and RCC.
CharacteristicsVariablesBCa
(n = 29)		PCa
(n = 24)		RCC
(n = 12)
Age (years) 67.92 ± 10.4268.21 ± 5.5462.17 ± 11.76
Sex (%)Male29 (100)24 (100)12 (100)
SCr(mg/dL) 0.98 ± 0.280.79 ± 0.130.98 ± 0.19
eGFR (mL/min/1.73 m2) 81.81 ± 21.11101.18 ±19.4978.30 ±14.98
pT stage (%)Ta13 (44.83)--
Tis2 (6.89)--
T111 (37.93)3 (12.5)7 (58.33)
T23 (15.15)6 (25.0)-
T3-14 (58.33)4 (33.33)
T4-1 (4.17)1 (8.33)
WHO/ISUP
classification (%)		Low15 (48.27)--
High14 (51.72)--
GS (%)6-5 (20.83)-
7-3 (12.50)-
8-8 (33.33)-
Above 8-8 (33.33)-
HSPC/CRPC (%) -8 (33.3)/16 (66.7)-
Fuhrman grade (%)G2--8 (66.67)
G3--2 (16.67)
G4 2 (16.67)
Abbreviations: SCr, Serum creatinine; eGFR, estimated glomerular filtration rate; pT, Pathologic Stage; GS, Gleason Score; HSPC, Hormone Sensitivity Prostate Cancer; CRPC, Castration Resistant Prostate Cancer. Note: Data are presented as mean ± standard deviation. eGFR was calculated using by MDRD formula.
Table
Table 2. Significance assessment of univariate receiver operating characteristic (ROC) curve analysis for metabolites potentially distinguishing between BCa, PCa, and RCC.
Table 2. Significance assessment of univariate receiver operating characteristic (ROC) curve analysis for metabolites potentially distinguishing between BCa, PCa, and RCC.
MetabolitesBCa vs. RCCPCa vs. BCaPCa vs. RCC
AUCp-Value95% CIAUCp-Value95% CIAUCp-Value95% CI
4-HBA0.5830.6660.401–0.7490.7310.0020.718–0.9400.7080.0410.517–0.892
N-MH0.6860.0440.500–0.8720.8408.875 × 10−60.717–0.9460.9171.667 × 10−50.824–1.0
Creatinine0.7140.0070.486–0.9140.7400.0020.587–0.8560.5490.3380.22–0.714
Glutamine0.8152.054 × 10−40.655–0.9470.6720.0240.522–0.810.9241.364 × 10−50.82–0.989
Acetate0.7390.0540.574–0.9040.6170.3450.454–0.7740.8480.0100.691–0.962
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, S.; Ku, J.Y.; Kang, B.J.; Kim, K.H.; Ha, H.K.; Kim, S. A Unique Urinary Metabolic Feature for the Determination of Bladder Cancer, Prostate Cancer, and Renal Cell Carcinoma. Metabolites 2021, 11, 591. https://doi.org/10.3390/metabo11090591

AMA Style

Lee S, Ku JY, Kang BJ, Kim KH, Ha HK, Kim S. A Unique Urinary Metabolic Feature for the Determination of Bladder Cancer, Prostate Cancer, and Renal Cell Carcinoma. Metabolites. 2021; 11(9):591. https://doi.org/10.3390/metabo11090591

Chicago/Turabian Style

Lee, Sujin, Ja Yoon Ku, Byeong Jin Kang, Kyung Hwan Kim, Hong Koo Ha, and Suhkmann Kim. 2021. "A Unique Urinary Metabolic Feature for the Determination of Bladder Cancer, Prostate Cancer, and Renal Cell Carcinoma" Metabolites 11, no. 9: 591. https://doi.org/10.3390/metabo11090591

APA Style

Lee, S., Ku, J. Y., Kang, B. J., Kim, K. H., Ha, H. K., & Kim, S. (2021). A Unique Urinary Metabolic Feature for the Determination of Bladder Cancer, Prostate Cancer, and Renal Cell Carcinoma. Metabolites, 11(9), 591. https://doi.org/10.3390/metabo11090591

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop