Next Article in Journal
Prognostic Significance of Preoperative NLR, MLR, and PLR Values in Predicting the Outcome of Primary Cytoreductive Surgery in Serous Epithelial Ovarian Cancer
Previous Article in Journal
Usefulness of Contrast-Enhanced Endoscopic Ultrasound (CH-EUS) to Guide the Treatment Choice in Superficial Rectal Lesions: A Case Series
Previous Article in Special Issue
Association between the Immunophenotype of Peripheral Blood from mCRPC Patients and the Outcomes of Radium-223 Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Urinary Metabolomic Analysis of Prostate Cancer by UPLC-FTMS and UPLC-Ion Trap MSn

1
Molecular Medicine Research Center, Chang Gung University, Kwei-San, Taoyuan 33302, Taiwan
2
Department of Urology, Chang Gung Memorial Hospital, College of Medicine, Chang Gung University, Kwei-San, Taoyuan 33305, Taiwan
3
Department of Biomedical Sciences, College of Medicine, Chang Gung University, Taoyuan 33302, Taiwan
4
Department of Cell and Molecular Biology, College of Medicine, Chang Gung University, Kwei-San, Taoyuan 33302, Taiwan
*
Author to whom correspondence should be addressed.
Submission received: 31 May 2023 / Revised: 25 June 2023 / Accepted: 29 June 2023 / Published: 4 July 2023
(This article belongs to the Special Issue Biomarkers of Urological Disease)

Abstract

:
Accumulative evidence suggests metabolic disorders correlate with prostate cancer. Metabolic profiling of urine allows the measurement of numerous metabolites simultaneously. This study set up a metabolomic platform consisting of UPLC-FTMS and UPLC-ion trap MSn for urine metabolome analysis. The platform improved retention time, mass accuracy, and signal stability. Additionally, the product ion spectrum obtained from ion trap MSn facilitated structure elucidation of candidate metabolites, especially when authentic standards were not available. Urine samples from six hernia patients and six BPH patients were used for the initial establishment of the analytic platform. This platform was further employed to analyze the urine samples of 27 PCa and 49 BPH patients. Choosing the upper and lower 16% of metabolites, 258 metabolite candidates were selected. Twenty-four of them with AUC values larger than 0.65 were further selected. Eighteen of the twenty-four features can be matched in METLIN and HMDB. Eleven of the eighteen features can be interpreted by MSn experiments. They were used for the combination achieving the best differential power. Finally, four metabolites were combined to reach the AUC value of 0.842 (CI 95, 0.7559 to 0.9279). This study demonstrates the urinary metabolomic analysis of prostate cancer and sheds light on future research.

1. Introduction

According to Globocan 2020, prostate cancer is the most frequently diagnosed cancer and the second leading cause of cancer death in men. There were an estimated 1.4 million new cases of prostate cancer worldwide in 2020 [1]. Prostate cancer tends to develop in men aged 65 and over [2]. It has been shown that metabolic syndrome and obesity correlate with benign prostatic hyperplasia (BPH) and prostate cancer (PCa) [3]. Metabolic analysis of prostate disease is urgently needed.
Prostate cells possess a unique metabolism specialized for producing, storing, and secreting PSA, polyamines, and citrate [4,5,6]. PCa is distinguishable from BPH and normal prostate by its high lactate and low citrate levels; hence, metabolomic analysis is ideal for studying the disease [7,8,9,10,11,12,13,14].
Various analytic platforms have been used for metabolomic analyses [9,10,15,16,17,18,19]. However, nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry (MS) are the most common tools. The advantage of MS methods over NMR is much higher sensitivity and detection of metabolites at much lower concentrations. Moreover, the acquired sample information provided by MS is specific to each metabolite, while NMR provides information about the chemical structure and quantitative information about metabolite concentrations.
MS is generally hyphenated with gas or liquid chromatography (GC, LC). For GC-MS analysis, derivatization is important to improve the volatility of metabolites. Sometimes, derivatization can also result in metabolite degradation [19]. LC-MS does not need derivatization and is widely applied in metabolomic analysis. Due to the sample complexities of metabolomic analysis, analysis by conventional LC was swiftly superseded by ultra-performance liquid chromatography (UPLC). UPLC equipped with a sub-2-μm column is acclaimed for its great resolution, high sensitivity, and rapid analytical speed [20,21]. For metabolomic analysis, Q-Tof or Tof analyzers are commonly utilized to match the fast speed of UPLC [22,23]. Though the accuracy of Q-Tof and Tof MS are listed in the range of 1~5 ppm [24], it was found that errors could be enlarged. Sometimes, the error was as significant as 20 ppm in the metabolomic analysis [25]. However, the first step for chemical structural identification of metabolites using database searching needs accurate results to remove false positive candidates. Fourier transform mass spectrometry (FTMS) is an accurate and stable analyzer. The analyzer is shielded by an isothermal shell and provides high resolution and highly sensitive detection. The accuracy is generally in the range of 1.5 to sub-ppm. These characteristics facilitate the first step of metabolite identification in the metabolomic analysis.
Generally, two independent and orthogonal data related to an authentic compound analyzed under identical experimental conditions are proposed as a necessity to validate non-novel metabolite identifications [26]. In addition to accurate mass, a double check is essential when identifying metabolites. For some commercially available metabolites, identifications rely on matching unknown compounds with authentic standards on the retention and the accuracy mass. However, most of the time, many metabolites on the search list are unavailable. Therefore, another possibility for metabolite identification is performing a product ion experiment to identify the structure information if the authentic compound is unavailable.
If a metabolite has been identified only through a parent ion mass match to a database, its identity should only be considered “very tentative”. To confirm metabolite candidates, researchers are required to purchase a lot of standard compounds. On the other hand, if a compound has been identified by matching the retention time, parent ion mass, and mass fragment patterns of a known standard, the identification of this compound is confident. To avoid the need to purchase many standards, we employed MSn to facilitate our protocol in two aspects. The first is to narrow down the candidate list before matching the endogenous metabolites using standard compounds, and the second is to provide more confidence for the metabolites if the authentic compounds are unavailable [27].
FTMS is not the best for MSn experiments to analyze complex samples with a wide range of concentrations [28]. The fringe field of the high field magnet may reject large ions, and the high radio frequency of small product ions may result in poor trapping efficiency of small ions. These reasons make the sensitivity of the MSn experiment of FTMS not satisfactory, especially for those low abundant and small or high molecular weight metabolites. To overcome these limitations, we employed an ion trap mass spectrometer as the analyzer for MSn analysis. Ion trap MS has the ability of ion storage. Even though the ion is low abundant, it could be collected longer to accumulate enough ion signals [29]. The scanning speed of ion trap MS is faster than FTMS. Thus, the ion detection efficiency is better for ion trap MS in the low mass region. Additionally, the MSn in ion trap MS could provide the structural information of metabolite candidates. Collectively, the sensitivity of the product ion spectrum from ion trap MS is better than those obtained from FTMS.
This study established the UPLC-FTMS and UPLC-ion trap MSn that might better resolve and characterize metabolically related differences between PCa and BPH.

2. Materials and Methods

2.1. Sample Collection and Preparation

All urine samples were collected at Chang Gung Memorial Hospital, Taoyuan, Taiwan. The study protocol was approved by the Medical Ethics and Human Clinical Trial Committee at Chang Gung Memorial Hospital. All urine samples were collected on the first morning after admission and before any operation. As the normal range of urine specific gravity is 1.003~1.030, urine samples were discarded if out of this range. All prostate cancer patients were diagnosed before the commencement of androgen deprivation or any focal therapy such as radiation. For the concern of metabolic interference, any BPH patient exposed to 5- α -reductase inhibitor was also excluded from this study. The collection, processing, and storage of urine samples were described in our previous publication [30]. Briefly, the urine was collected in the presence of a protease inhibitor cocktail tablet (one tablet per 50 mL urine) and sodium azide (1 mM). Removal of cells and debris by centrifugation (5000× g for 30 min at 4 °C) was carried out within one hour after sample collection. The sample was then kept at −80 °C for long-term storage. Before MS analysis, the sample was thawed at 4 °C. Freezing point depression was measured to determine the osmolarities of samples using a model 3320 osmometer (Advanced Instruments, Norwood, MA, USA). In total, 50 μL of urine were diluted four times with methanol and centrifuged at 13,200× g for 15 min at 4 °C. The supernatant was dried under N2. The sample was re-dissolved with 50 μL solvent consisting of MeOH: H2O = 2:1. The sample was centrifuged at 13,200× g for 15 min at 4 °C again; the supernatant was collected and used for LC-MS analysis.
In the initial establishment of the metabolomic analysis, the urine samples from 6 hernia patients (age 50~68, mean 62.7) and 6 BPH patients (age 64~77, mean 68.3) were used. They were all male with a comparable age range. This platform was further employed to differentiate the metabolomic profile of the prostatic disease between BPH and prostate cancer.

2.2. Chemicals

Protease inhibitor cocktail tablets were purchased from Roche (Mannheim, Germany). Solvents for the sample preparation and LC-MS analysis, such as methanol, water, acetonitrile with 0.1% formic acid, and water with 0.1% formic acid, were purchased from Sigma–Aldrich (St. Louis, MO, USA). These solvents were of CHROMASOLV grade. N2, N2-Dimethylguanosine, 7-Ketodeoxycholic acid, and Phenylacetyl L-glutamine were purchased from Toronto Research Chemicals Inc. (North York, ON, Canada). The other chemicals were all purchased from Sigma–Aldrich (St. Louis, MO, USA). The peptide was synthesized by Kelowna International Scientific Inc. (Kaohsiung, Taiwan).

2.3. Instrumentation

A UPLC system (Waters, Milford, MA, USA) equipped with a C18 reverse phase column (2.1 × 100 mm, 1.8 μm, HSS-T3, Waters, Milford, MA, USA) was coupled with either an FT-ICR MS (Apex-Q 94, Bruker Daltonics, Billerica, MA, USA) or an ion trap MS (HCT, Bruker Daltonics, Billerica, MA, USA) in the positive ESI mode. The flow rate was 0.1 mL/min at the beginning (99% solvent A: 0.1% formic acid; 1% solvent B: acetonitrile with 0.1% formic acid). One μL sample was injected. After injection, 1% solvent B was held for 5 min and increased to 50% in 9 min, then to 90% in 6 min, and 99% in 12 min held for 2 min, and the flow rate was increased to 0.5 mL/min, 5 min later. The flow rate was then reduced to 0.1 mL/min and back to 1% in 0.1 min and held for 7 min.
The mass spectra were recorded using the full scan mode, covering a mass range of m/z 100–1000. The sample size was 256 k, the calibration standard was 10 mM sodium formate in isopropanol: H2O = 9:1. The resolution was 17,944 at m/z 430.9138 (the ion formula is Na(NaCOOH)6), the source accumulation time was 0.01 s, and the collision cell time was 0.08 s. The capillary voltage was set to 4200 V, and the spray shield voltage was set to 3500 V. The nebulizing gas and drying gas flow were set to 4.0 and 4.6 L/min, respectively. The drying temperature was 190 °C.
An equal amount of each sample in the analysis was mixed to be the QC sample. Before batch analysis, the QC sample was injected 10 times for the condition of the UPLC column, and the sample injection sequence was randomized to reduce the effect of contamination from the previous injections. The sample injection sequence was randomized according to the suggestions of Wang EJ et al. [31]. In the injection sequence, every 10 analyses, the QC sample was injected to check the system stability throughout the whole analysis.

2.4. Data Analysis

Principal component analysis (PCA) and t-test analyses were performed by the software of ProfileAnalysis 2.0 (Bruker Daltonics, Billerica, MA, USA). The MS data from 0.5 to 17 min were used for analysis. The advanced bucket for time and mass was set to 0.3 s and 5 mDa, respectively. The S/N threshold was set to 5 under the “calculate compound” mode. The raw data of LC-MS was imported into the ProfileAnalysis, and PCA analysis was performed first. Additionally, the “count number” filter was applied. The “count number” means the encounters of a feature detected in the sample batch. It is used for reducing the low confident features. In this study, about half number of one group sample was used. If the sample numbers were different for different sample groups, the smallest sample number was selected to set the criteria. Then, the features were filtered by the “count number”, and the t-test calculation was performed. All tests were two-sided, and a p-value less than 0.05 was considered statistically significant.

2.5. Normalization and Significance Criteria

The numbers of detected features vs. their log (ratio) plot were drawn for the evaluation of the need for normalization. By fitting the curve with Gaussian distribution, the SD and shift (x0) were determined. If the x0 to 0 was less than |±0.04|it meant that the error of ratio was smaller than 10%, and the normalization was not processed. Using 0 ± 1 SD as the criteria, the top and bottom 16% features were selected for further analysis.

2.6. Web Searching

The mass list was searched against two web databases, the METLIN and the HMDB databases [32,33,34]. The mass tolerance is 5 ppm for METLIN and 0.01 Dalton for HMDB. The formulas of candidates obtained from the search were compared with the calculated isotope pattern of the corresponding feature. The matched candidates were further investigated by MS2 and/or MS3 experiments.

2.7. Statistical Analysis

Differences in concentration levels of urinary metabolites between different groups were analyzed using the nonparametric Mann–Whitney U-test. In addition, SPSS software was also used for statistical analyses of receiver operator characteristic (ROC) curves and areas under the curve (AUC) to check the abilities of metabolites to discriminate different groups.

3. Results

3.1. Workflow for Urinary Metabolomic Analysis by Integration of UPLC-FTMS and UPLC-Ion Trap

Figure 1 demonstrates the MS-based metabolic analysis workflow for urine in this study. A UPLC system was employed for the separation of metabolites. The FTMS and ion trap MSn (n = 2~3) were for the exact mass detection and structure determination, respectively.

3.2. Optimization of UPLC-FTMS System

FTMS is chiefly characterized by its high mass accuracy, high resolution, and high sensitivity, and these characteristics make it an excellent tool for small molecule detection. However, FTMS needs a longer time to collect stable signals (cycle time > 5 s) than other types of mass spectrometers; for example, the Q-Tof mass spectrometer can be operated at a speed of 0.1 s/scan [35]. If the signal collection time of FTMS is not long enough, the signal will be poor or fluctuate. When MS was hyphenated with chromatography, the first criterion for detecting a good peak shape is that the eluted peak width should be wide enough, and the MS can detect enough data points to describe the original peak shape. Generally, 10~20 data points for a peak are preferred.
The advantages of UPLC were recognized for its high speed and high resolution; moreover, we found that it also provided high accuracy in the retention time of peaks. Untargeted metabolomic analysis has the advantage of discovering novel metabolites. There is no internal standard added to the samples, and the analyzed results are directly compared to determine the difference. The analytical instruments must be consistent and accurate throughout the duration of the analysis in order to produce results that are comparable between sessions. In this circumstance, the stability of UPLC can help in the untargeted metabolomic analysis.
In connecting UPLC to FTMS, the drawback was the incompatibility of analysis speed in common analytical conditions. For example, the peak width of UPLC is generally about 1.5~2.5 s (FWHM) in a flow rate of 0.6 mL/min on a 2.1 × 100 mm column, and the cycle time of FTMS is about 4.5 s at the resolution of 32,000 (512 k, single spectrum). According to the Savitzky–Golay smoothing filtering method, an LC peak needs at least 11 data points to describe its sharpness for quantitation [36]. Ten of the eleven data points are used for peak description, and one more for baseline determination (Figure 2a).
In the time frame of UPLC peaks, the FTMS should have a cycle time of less than 0.3~0.5 s. Nevertheless, in the time frame of FTMS, the separation speed of UPLC should slow down to make the LC peak wider than 45 s. ESI is known to be a concentration-dependent detector. When we increase the width of an LC peak by reducing the flow rate of LC, the peak height is inversely proportional to the peak width, and the signal of the mass spectrometer decreases with the peak height decreasing.
To resolve this challenge, we reduced the flow rate of UPLC to provide a longer time for FTMS scanning, speed up the scanning of FTMS to obtain more data points, and restrain the peak width from maintaining the sensitivity of ESI detection. The cycle time was reduced from 4.5 to 3 s/spectrum to speed up the scanning speed of FTMS. Under such conditions, 10 data points would take about 30 s to obtain a stable signal in the FTMS system; it took too long for the UPLC separation, and the mass sensitivity dropped significantly. Therefore, the data points were reduced to five to save the separation time. As illustrated in Figure 2b, when the data points of a peak were reduced to five points, the peak shape could be almost maintained, and the sensitivity would not be significantly reduced. If the data points were further reduced to four points, the peak shape would be distorted, and the sensitivity and resolution would be decreased (Figure 2c). The flow rate used for the 2.1 mm × 100 mm UPLC column was generally in the range of 0.4~0.6 mL/min, and the peak width was about 1.5~2.5 s (FWHM). After the flow rate adjustment, the flow rate was reduced to 0.1 mL/min to widen the peak width to 7.5 s (FWHM), and six data points could be detected for an ion with an intensity of 2 × 105 counts in our system (Figure 2d), which could provide a minimum data point for describing the sharpness of LC peak.
As suggested by the Savitzky–Golay smoothing filtering method, when the data point was less than 10, the shoulder peak could not be drawn perfectly, and the quantitation might have problems. However, because of the high resolving power of FTMS, the high mass resolution could help in separating two overlapping peaks in LC (in the time domain) into two distinct peaks in MS (in the mass domain). Only when the co-eluted metabolites are structural isomers, the quantitation will fail. Additionally, UPLC is also a high-resolution instrument. Poor separation for a UPLC is minor compared to that of a conventional LC. All the features detected in the metabolomic analysis are currently checked manually for overlaps.

3.3. Analysis of Urine Samples from BPH and Hernia Patients

Urine samples from six BPH and six hernia patients were used to test the platform. In this analysis batch, 79 analyses (6 BPH samples × 6 injections, 6 hernia samples × 6 injections, and 7 QC runs) were run after the system was conditioned with 10 QC samples.
A typical base peak chromatogram of BPH urine is shown in Figure 3a. Generally, the urine metabolites will be separated into two groups in the LC analysis. One group was eluted at about 1.8~6.0 min; the other group was eluted at 9.5~19.0 min. To evaluate the performance of the UPLC-FTMS platform, seven QC runs in the analysis sequence were used to check the data qualities. In the QC samples, we selected four features at different retention times with different mass values and intensities. They were m/z 205.0685, 334.1126, 551.2124, and 407.2798 ions at retention times of 2.2, 10.3, 12.1, and 16.8 min, respectively (Figure 3a). Their retention times were very accurate, their standard deviations were all 0.0 min, and the standard deviations of detected mass values were less than 0.9 ppm. Moreover, their standard deviations of intensities were all less than 14%. These four features were later identified as D-mannitol, N2, N2-dimethylguanosine, phenylacetyl L-glutamine, and 7-ketodeoxycholic acid (Supplementary Table S3), respectively. Their actual mass errors were all less than 1.3 ppm. The results indicate that our metabolomic analysis system was reliable and accurate within the analysis batch.
The principal component analysis (PCA) plot for the 79 analyses is shown in Figure 3b. Each point in the PCA plot represents one analysis of a sample. It appeared that all the QC samples were clustered together, demonstrating that the analysis conditions were consistent and reproducible throughout the whole analysis time. A clear separation of BPH and Her sample groups could be observed. Additionally, for each sample, although the sample injection sequence had been randomized, the six analytical replicates for each sample were still clustered closely in the plot. The results further demonstrated that our platform is stable and reliable for urinary metabolomic analysis.
After peak alignment, 17,708 features were detected in the PCA analysis. To remove the non-reproducible signals, we set the filter to half of the lower limit of the sample number for feature detection. In this analysis, the total analysis numbers of BPH and Her samples were all 36 (6 samples × 6 injections); therefore, the filter was set to 18. After filtration, 1010 features were detected.
First, the resolution of peak separation was manually checked to determine if the reduction of MS data points would affect the quality of LC-MS separation in the view of extracted ion chromatogram. Among the 1010 features, 56 features (5.5% of features) were found to be not baseline separated in their extracted ion chromatograms. The error rate for the relative quantitation experiment is about 5.5%. To avoid false positive results, we suggested that when potential markers were found, it was necessary to recheck the extracted ion chromatograms of the markers; then, the errors for quantitation could be largely reduced.

3.4. Analysis of Urine Samples from PCa and BPH Patients

Using the platform established as described above, we further analyzed urine samples from 27 PCa and 49 BPH patients for the discovery of differential metabolites. The demographic data of the patients enrolled in this study are shown in Supplementary Table S1a,b. According to the suggestion of the Chemical Analysis Working Group, biological replicates are preferred over analytical replicates, as biological variance consistently exceeds analytical variance [26]. The sample number for each group in the study was more than 27 cases, and our metabolomic platform has been shown to be highly stable; thus, each sample was analyzed once. Figure 4a shows the PCA plot of these analyses. Again, all the QC samples were clustered together in this analysis which helped us to confirm the system stability. However, the separation of BPH and PCa samples in the PCA plot was not differentiated clearly.
A total of 15,662 features were detected in the PCA analysis. Again, the counted number of features was set to half of the PCa sample number. In this study, each sample was evaluated once. Hence, the counted number of features filter was set to 14. Following filtration, 617 features were identified, and a t-test was performed on these features.
Figure 4b is the plot of feature number vs. the log ratioBPH/PCa. The curve was fitted with Gaussian distribution. The x0 was −0.0008, which meant that the shift in ratioBPH/PCa was less than 0.2%. According to the experiment results, data were not necessary to carry out the normalization. In the analysis, the b parameter was 0.10, equal to the standard deviation (SD). ±1 SD was utilized to identify candidates with more significant fold changes. It was equal to the ratios of 1.2 or 0.8. The features with ratioBPH/PCa either larger than 1.20 or smaller than 0.80 were chosen. A total of 258 features were selected. We further applied the ROC calculation to these features. By filtering with the medium differential power AUC value, we found 24 features with an AUC value larger than 0.65.
Even though the analysis was performed at high mass accuracy (<1.3 ppm), many chemically possible formulas could be obtained, especially in higher mass regions. A simple way to identify the candidates of features is to combine high mass accuracy data with database searching in general mass spectral databanks or biochemical databases. We selected two human-related databases, METLIN and HMDB [32,33,34], for database search, producing a shortlist of potential candidates. The list can be shortened further by applying isotope ratios and performing MSn experiments. As shown in Figure 1, the first step of metabolite identification was to check the feature’s exact mass and the isotope pattern. The exact mass can be utilized to search against the database, and the isotope pattern could be used to calculate the metabolite formula. However, in the LC-MS analysis, the metabolites are very likely to form clusters or metal adducts, and different element compositions will affect the calculation results of the isotope pattern. If the database is too complicated, the search results will also be complicated for the high possibility of false positive matching. By choosing a suitable database, candidates for metabolites were selected in advance based on the database’s properties. Thus, the possibility of false positive identification will be reduced. Furthermore, the resulting list of candidates will focus on those highly relevant metabolites and reduce the complexity of the candidate list. After the web search, 18 out of 24 features can be matched to compounds collected in METLIN and/or HMDB databases.
Among the 18 matched features, 11 can be interpreted by their product ion spectra (Supplementary Table S3). For example, a feature was detected at 10.4 min with m/z 334.1135. The mass value was used to search against the METLIN database. The search results showed that there were 10 possible candidates. By checking their formula, five categories of compounds and five element compositions were possible. They were imazethapyr (a pesticide) ([M + 2Na − H]+), 4-nitrophenol-alpha-D-galactopyranoside ([M + CH3OH + H]+), three sialic acid derivatives ([M + H-H2O]+), two dipeptides, and three nucleosides (one [M + ACN + Na]+ and two [M + Na]+). According to these formulas, the elements of C, H, O, N, and Na were used for calculating the possible isotope ratio of m/z 334.1135 ions. By matching the formulas obtained from the database search with the results from the isotope ratio calculation, we found that only three formulas (eight compounds) were possible. They were [M + − H2O]+, [M + ACN + Na]+, and [M + Na]+. We then performed an MS2 experiment to further explore the metabolite. As shown in Figure 5a, the fragment ions were very simple; only m/z 266, 202 ions were detected. Some clues can be found in the product ion spectrum: (1) no fragment ion was observed because of the loss of acetonitrile (Δm = 41); (2) the most abundant fragment ion resulted from the elimination of a fragment of mass 132. According to clue 1, the formula [M + ACN + Na]+, with solvent and sodium adducts, was not a possible formula. Two dipeptides and one nucleoside were excluded. In the product ion spectrum of m/z 344.11, no water loss was observed (Δm = 18). These sialic acid derivatives included five hydroxyl groups and lost water molecules readily during collision-induced dissociation (CID) processes. According to clue 2, we tentatively excluded the sialic acid derivatives.
Finally, two potential nucleosides remained, N2, N2-Dimethylguanosine and 1,7-Dimethylguanosine. Since they were all nucleosides, it made sense that these compounds may lose ribosyl moieties (Δm = 132) in their respective product ion spectra. By examining the fragment ions in the product ion spectrum of m/z 334.1134 in further detail, we identified an ion of mass 266. It was hypothesized that the ion would be generated by the breaking of the C–N bond on the side chain of N2-dimethylguanosine’s nucleobase. It appeared unlikely that the ion would result from the cleavage of 1, 7-dimethylguanosine. The identification of the metabolite as N2, N2-Dimethylguanosine was validated by comparison to its authentic standard (Figure 5b).
Table 1 lists the details of these 11 features (10 metabolites), and Table 2 lists the statistical results of the 10 metabolites. All 10 metabolites have medium differential power for differentiating the PCa and BPH patients. By further checking their p-values obtained from the t-test and fold change, there were seven metabolites with a p-value less than 0.05, and four of them were detected with a fold change of more than 2.0 or less than 0.52. All four metabolites have AUC values of more than 0.658. These four metabolites were used to reconstruct a new PCA plot. Figure 4c demonstrates that the BPH and PCa groups can be separated more clearly than their original PCA plot. The four metabolites were one dipeptide, a Vitamin D3 derivative, 7-Dehydropregnenolone, and L-Urobilinogen.
The dipeptide, Ala-Asp, was synthesized and confirmed by matching the retention time, accurate mass, and product ion spectrum (Supplementary Table S3). The standard compounds of the Vitamin D3 derivative, 7-Dehydropregnenolone and L-Urobilinogen, were unavailable. The Vitamin D3 derivative had four possible candidates. An MS3 experiment was then used to confirm the identity of the Vitamin D3 derivative. As shown in Figure 6a, the product ion spectrum of m/z 459 ions can be interpreted by the fragmentation of four possible Vitamin D3 derivatives (Figure 6b). They had different side chains and chiral centers. Unfortunately, there was no fragment ion with which the side chain structure could be confirmed. The m/z 367 ion was a hydronium adduct ion of the cleavage of Vitamin D3 derivative at the thioether bond. The m/z 367 ion in the product ion spectrum was selected for further fragmentation. The product ion spectrum of m/z 367 showed two fragment ions (Figure 6c); the m/z 313 ion resulted from losing three water molecules, and the m/z 155 ion from losing the fused ring. MS2 and MS3 experiments confirmed the skeleton of Vitamin D3, but the structure of its side chain and chirality could still not be determined by these MSn experiments.

3.5. Clinical Narratives

Vitamin D (Vit D) deficiency has been linked to prostate cancer [37]. Sunlight exposure may boost the skin’s vitamin D production. Geographically, the incidence of prostate cancer is lowest in tropical zones, where exposure to sunshine is highest, compared to the subtropical and temperate zones [1]. Skin pigmentation influences cutaneous Vit D production. Black men in North America have a higher incidence of prostate cancer than white men.
Urobilinogen is a breakdown product of bilirubin, which is formed when red blood cells are broken down in the liver. Urobilinogen is then excreted in the urine and feces. The presence of urobilinogen in the urine can indicate liver or gallbladder problems, while the absence of urobilinogen in the urine can indicate a problem with red blood cell production or destruction. There is no direct relationship between urobilinogen and vitamin D. However, some liver disorders that affect the production of urobilinogen may also affect the absorption or metabolism of vitamin D. Therefore, it is important to monitor the levels of both urobilinogen and vitamin D in individuals with liver disorders. In the current study, it is intriguing to note that the urinary metabolites of vitamin D and urobilinogen are lower in PCa patients than in BPH patients (Table 2, BPH/PCa Ratio).
Cytochrome P450scc, in vitro or ex vivo in adrenal glands, cleaves the side chain of 7-dehydrocholesterol, producing 7-dehydropregnenolone. 7-dehydropregnenolone and its hydroxyderivatives can undergo UVB-induced transformation to corresponding pregnacalciferol and lumisterol-like compounds [38]. 7-dehydropregnenolone is a precursor to several hormones in the body. It can be converted into testosterone, which is important for male sexual development and function.
Combing the four metabolites, the AUC value reached 0.842 (CI 95, 0.7559 to 0.9279, Figure 4d). The value was in the middle–high range of the ROC analysis.

4. Conclusions

A platform comprised of UPLC-FTMS and UPLC-ion trap MSn has been established. The platform is reliable and capable of conducting high-quality urinary metabolomic analysis. The identification of endogenous metabolites was more certain with MSn prior to confirmation with authentic standards. The platform can also offer guidance for synthesizing the standards; not even the authentic standards were accessible. This study detected one peptide, one Vitamine D3 derivative, 7-dehydropregnenolone, and L-urobilinogen as potential biomarkers of prostate cancer. After combing the four metabolites, the AUC value reached 0.842 (CI 95, 0.7559 to 0.9279). The underlying metabolic pathway could elicit future investigation. Instead of a clinical approach, this study demonstrates urine metabolomic analysis of disease that could be an adjunct of disease diagnosis and shed light on future research.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/diagnostics13132270/s1, Table S1a: Clinical data of 27 PCa patients and 49 BPH patients with analysis results of four metabolites; Table S1b: The demography of patients with benign prostatic hyperplasia (BPH) and prostate cancer (PCa); Table S2: The product ion spectra of metabolites in Table 1; Table S3: The product ion spectra of identified metabolites and their authentic standards.

Author Contributions

Conceptualization, C.-L.C. and B.-R.J.; methodology, Y.-T.C., W.-Y.L. and B.-R.J.; software, W.-Y.L. and B.-R.J.; validation, Y.-T.C. and B.-R.J.; formal analysis, B.-R.J.; investigation, B.-R.J.; resources, B.-R.J.; data curation, W.-Y.L. and B.-R.J.; writing—original draft preparation, B.-R.J.; writing—review and editing, C.-L.C.; supervision, Y.-S.C. and J.-S.Y.; visualization, C.-L.C. and B.-R.J.; funding acquisition, C.-L.C. and Y.-T.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Chang Gung Memorial Hospital (CMRPD1L0091), and the Ministry of Science and Technology (MOST), Taiwan (NSC 111-2113-M-182-002).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of Chang Gung Medical Foundation Institutional Review Board (protocol code 1041606C1 and date of approval: 21 January 2010).

Informed Consent Statement

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

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AUCArea Under Curve
BPHBenign Prostatic Hyperplasia
CIDCollision-Induced Dissociation
ESIElectrospray Ionization
FTMSFourier Transform Mass Spectrometry
FWHMFull Width at Half Maximum
HerHernia
HMDB The Human Metabolome Database
LC-MS Liquid Chromatography-Mass Spectrometry
METLIN The METLIN Metabolite and Chemical Entity Database
(https://metlin.scripps.edu/landing_page.php?pgcontent=mainPage (accessed on 29 May 2023))
MSnMass Spectrometry in n (n = 2~3) rounds
NMRNuclear Magnetic Resonance
PCAPrincipal Component Analysis
PCaProstate Cancer
QCQuality Control
Q-Tof Quadruple Time-of-flight
ROC curveReceiver Operating Characteristic curve
UPLCUltra-Performance Liquid Chromatography

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef]
  2. Siegel, R.L.; Miller, K.D.; Fuchs, H.E.; Jemal, A. Cancer statistics, 2022. CA Cancer J. Clin. 2022, 72, 7–33. [Google Scholar] [CrossRef] [PubMed]
  3. De Nunzio, C.; Aronson, W.; Freedland, S.J.; Giovannucci, E.; Parsons, J.K. The correlation between metabolic syndrome and prostatic diseases. Eur. Urol. 2012, 61, 560–570. [Google Scholar] [CrossRef] [PubMed]
  4. Takyi, E.E.; Fuller, D.J.; Donaldson, L.J.; Thomas, G.H. Deoxyribonucleic acid and polyamine synthesis in rat ventral prostrate. Effects of age of the intact rat and androgen stimulation of the castrated rat with testosterone, 5 alpha-dihydrotestosterone and 5 alpha-androstane-3 beta, 17 beta-diol. Biochem. J. 1977, 162, 87–97. [Google Scholar] [CrossRef]
  5. Schick, F.; Bogers, H.; Kurz, S.; Jung, W.I.; Pfeffer, M.; Lutz, O. Localized proton MR spectroscopy of citrate in vitro and of the human prostate in vivo at 1.5 T. Magn. Reson. Med. 1993, 29, 38–43. [Google Scholar] [CrossRef] [PubMed]
  6. Stenman, U.H.; Leinonen, J.; Zhang, W.M.; Finne, P. Prostate-specific antigen. Semin. Cancer Biol. 1999, 9, 83–93. [Google Scholar] [CrossRef]
  7. Trock, B.J. Application of metabolomics to prostate cancer. Urol. Oncol. 2011, 29, 572–581. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Truong, M.; Yang, B.; Jarrard, D.F. Toward the detection of prostate cancer in urine: A critical analysis. J. Urol. 2013, 189, 422–429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. 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] [Green Version]
  10. DeFeo, E.M.; Wu, C.L.; McDougal, W.S.; Cheng, L.L. A decade in prostate cancer: From NMR to metabolomics. Nat. Rev. Urol. 2011, 8, 301–311. [Google Scholar] [CrossRef]
  11. 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] [PubMed] [Green Version]
  12. Hyndman, M.E.; Mullins, J.K.; Bivalacqua, T.J. Metabolomics and bladder cancer. Urol. Oncol. 2011, 29, 558–561. [Google Scholar] [CrossRef] [PubMed]
  13. Ganti, S.; Weiss, R.H. Urine metabolomics for kidney cancer detection and biomarker discovery. Urol. Oncol. 2011, 29, 551–557. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Weiss, R.H.; Kim, K. Metabolomics in the study of kidney diseases. Nat. Rev. Nephrol. 2011, 8, 22–33. [Google Scholar] [CrossRef]
  15. Serkova, N.J.; Spratlin, J.L.; Eckhardt, S.G. NMR-based metabolomics: Translational application and treatment of cancer. Curr. Opin. Mol. Ther. 2007, 9, 572–585. [Google Scholar]
  16. Ryan, D.; Robards, K.; Prenzler, P.D.; Kendall, M. Recent and potential developments in the analysis of urine: A review. Anal. Chim. Acta 2011, 684, 8–20. [Google Scholar] [CrossRef]
  17. Theodoridis, G.A.; Gika, H.G.; Want, E.J.; Wilson, I.D. Liquid chromatography-mass spectrometry based global metabolite profiling: A review. Anal. Chim. Acta 2012, 711, 7–16. [Google Scholar] [CrossRef]
  18. Xiao, J.F.; Zhou, B.; Ressom, H.W. Metabolite identification and quantitation in LC-MS/MS-based metabolomics. Trends Anal. Chem. 2012, 32, 1–14. [Google Scholar] [CrossRef] [Green Version]
  19. Halket, J.M.; Waterman, D.; Przyborowska, A.M.; Patel, R.K.; Fraser, P.D.; Bramley, P.M. Chemical derivatization and mass spectral libraries in metabolic profiling by GC/MS and LC/MS/MS. J. Exp. Bot. 2005, 56, 219–243. [Google Scholar] [CrossRef] [Green Version]
  20. Mazzeo, J.R.; Neue, U.D.; Kele, M.; Plumb, R.S. Advancing LC Performance with Smaller Particles and Higher Pressure. Anal. Chem. 2005, 77, 460–467. [Google Scholar] [CrossRef] [Green Version]
  21. Wu, N.; Clausen, A.M. Fundamental and practical aspects of ultrahigh pressure liquid chromatography for fast separations. J. Sep. Sci. 2007, 30, 1167–1182. [Google Scholar] [CrossRef]
  22. Lenz, E.M.; Williams, R.E.; Sidaway, J.; Smith, B.W.; Plumb, R.S.; Johnson, K.A.; Rainville, P.; Shockcor, J.; Stumpf, C.L.; Granger, J.H.; et al. The application of microbore UPLC/oa-TOF-MS and 1H NMR spectroscopy to the metabonomic analysis of rat urine following the intravenous administration of pravastatin. J. Pharm. Biomed. Anal. 2007, 44, 845–852. [Google Scholar] [CrossRef] [PubMed]
  23. Plumb, R.; Castro-Perez, J.; Granger, J.; Beattie, I.; Joncour, K.; Wright, A. Ultra-performance liquid chromatography coupled to quadrupole-orthogonal time-of-flight mass spectrometry. Rapid Commun. Mass Spectrom. 2004, 18, 2331–2337. [Google Scholar] [CrossRef] [PubMed]
  24. Kind, T.; Fiehn, O. Advances in structure elucidation of small molecules using mass spectrometry. Bioanal. Rev. 2010, 2, 23–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Juo, C.G.; Chiu, D.T.; Shiao, M.S. Liquid chromatography-mass spectrometry in metabolite profiling. Biofactors 2008, 34, 159–169. [Google Scholar] [CrossRef] [PubMed]
  26. Sumner, L.W.; Amberg, A.; Barrett, D.; Beale, M.H.; Beger, R.; Daykin, C.A.; Fan, T.W.; Fiehn, O.; Goodacre, R.; Griffin, J.L.; et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics 2007, 3, 211–221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Shao, C.H.; Chen, C.L.; Lin, J.Y.; Chen, C.J.; Fu, S.H.; Chen, Y.T.; Chang, Y.S.; Yu, J.S.; Tsui, K.H.; Juo, C.G.; et al. Metabolite marker discovery for the detection of bladder cancer by comparative metabolomics. Oncotarget 2017, 8, 38802–38810. [Google Scholar] [CrossRef] [Green Version]
  28. Keller, B.O.; Sui, J.; Young, A.B.; Whittal, R.M. Interferences and contaminants encountered in modern mass spectrometry. Anal. Chim. Acta 2008, 627, 71–81. [Google Scholar] [CrossRef] [PubMed]
  29. Jonscher, K.R.; Yates, J.R., 3rd. The quadrupole ion trap mass spectrometer--a small solution to a big challenge. Anal. Biochem. 1997, 244, 1–15. [Google Scholar] [CrossRef] [PubMed]
  30. Juo, C.G.; Chen, C.L.; Lin, S.T.; Fu, S.H.; Chen, Y.T.; Chang, Y.S.; Yu, J.S. Mass accuracy improvement of reversed-phase liquid chromatography/electrospray ionization mass spectrometry based urinary metabolomic analysis by post-run calibration using sodium formate cluster ions. Rapid Commun. Mass Spectrom. 2014, 28, 1813–1820. [Google Scholar] [CrossRef]
  31. Want, E.J.; Wilson, I.D.; Gika, H.; Theodoridis, G.; Plumb, R.S.; Shockcor, J.; Holmes, E.; Nicholson, J.K. Global metabolic profiling procedures for urine using UPLC-MS. Nat. Protoc. 2010, 5, 1005–1018. [Google Scholar] [CrossRef] [PubMed]
  32. Smith, C.A.; O’Maille, G.; Want, E.J.; Qin, C.; Trauger, S.A.; Brandon, T.R.; Custodio, D.E.; Abagyan, R.; Siuzdak, G. METLIN: A metabolite mass spectral database. Ther. Drug Monit. 2005, 27, 747–751. [Google Scholar] [CrossRef] [PubMed]
  33. Wishart, D.S.; Knox, C.; Guo, A.C.; Eisner, R.; Young, N.; Gautam, B.; Hau, D.D.; Psychogios, N.; Dong, E.; Bouatra, S.; et al. HMDB: A knowledgebase for the human metabolome. Nucleic Acids Res. 2009, 37, D603–D610. [Google Scholar] [CrossRef]
  34. Wishart, D.S.; Tzur, D.; Knox, C.; Eisner, R.; Guo, A.C.; Young, N.; Cheng, D.; Jewell, K.; Arndt, D.; Sawhney, S.; et al. HMDB: The Human Metabolome Database. Nucleic Acids Res. 2007, 35, D521–D526. [Google Scholar] [CrossRef]
  35. Xie, G.; Plumb, R.; Su, M.; Xu, Z.; Zhao, A.; Qiu, M.; Long, X.; Liu, Z.; Jia, W. Ultra-performance LC/TOF MS analysis of medicinal Panax herbs for metabolomic research. J. Sep. Sci. 2008, 31, 1015–1026. [Google Scholar] [CrossRef] [PubMed]
  36. Savitzky, A.; Golay, M.J.E. Smoothing and Differentiation of Data by Simplified Least Squares Procedures. Anal. Chem. 1964, 36, 1627–1639. [Google Scholar] [CrossRef]
  37. Kanaan, Y.; Copeland, R.L. The link between vitamin D and prostate cancer. Nat. Rev. Cancer 2022, 22, 435. [Google Scholar] [CrossRef]
  38. Slominski, A.T.; Janjetovic, Z.; Fuller, B.E.; Zmijewski, M.A.; Tuckey, R.C.; Nguyen, M.N.; Sweatman, T.; Li, W.; Zjawiony, J.; Miller, D.; et al. Products of vitamin D3 or 7-dehydrocholesterol metabolism by cytochrome P450scc show anti-leukemia effects, having low or absent calcemic activity. PLoS ONE 2010, 5, e9907. [Google Scholar] [CrossRef]
Figure 1. Workflow for urinary metabolomic analysis by integration of UPLC-FTMS and UPLC-ion trap. The numbers in the parenthesis indicated the feature numbers after the filtration by different criteria. Colors in the embedded image 2 indicate grouping during analysis.
Figure 1. Workflow for urinary metabolomic analysis by integration of UPLC-FTMS and UPLC-ion trap. The numbers in the parenthesis indicated the feature numbers after the filtration by different criteria. Colors in the embedded image 2 indicate grouping during analysis.
Diagnostics 13 02270 g001
Figure 2. (a) A diagram of 11 data points describing a peak. According to the Savitzky–Golay smoothing filtering method, an LC peak needs at least 11 data points to describe its sharpness for quantitation. Ten data points are used for peak description and one point for the baseline determination. (b) When the data points are only five points for an LC peak, the data points can still keep the peak sharp, but the ability to find the shoulder peak will be lost. (c) If the data points are less than 5 points, for example, 4 data points, the peak will broaden, and the peak sensitivity will decrease after smoothing. The light blue dot line delineates the original sharp peak, and the dark blue solid line denotes the simulated peak after the data points are reduced. (d) An actual peak detected in our UPLC-FT MS system. Six data points were recorded for the LC peak.
Figure 2. (a) A diagram of 11 data points describing a peak. According to the Savitzky–Golay smoothing filtering method, an LC peak needs at least 11 data points to describe its sharpness for quantitation. Ten data points are used for peak description and one point for the baseline determination. (b) When the data points are only five points for an LC peak, the data points can still keep the peak sharp, but the ability to find the shoulder peak will be lost. (c) If the data points are less than 5 points, for example, 4 data points, the peak will broaden, and the peak sensitivity will decrease after smoothing. The light blue dot line delineates the original sharp peak, and the dark blue solid line denotes the simulated peak after the data points are reduced. (d) An actual peak detected in our UPLC-FT MS system. Six data points were recorded for the LC peak.
Diagnostics 13 02270 g002
Figure 3. (a) A typical base peak chromatogram of BPH urine. The four arrows indicate the retention times of four selected features. They were m/z 205.0685, 334.1126, 551.2124, and 407.2798 at retention times of 2.2, 10.3, 12.1, and 16.8 min, respectively. These four features are different in retention times, molecular weights, and ion intensities. They were used to evaluate the stability of the UPLC-FTMS. They were finally identified as D-mannitol, N2, N2-dimethylguanosine, phenylacetyl L-glutamine, and 7-ketodeoxycholic acid, respectively. The detected parameters are listed. (b) The PCA plot of BPH and Her urines. The symbol (Δ) denotes BPH, (x) Her, and (O) QC. The QC was made by mixing an equal amount of each urine sample.
Figure 3. (a) A typical base peak chromatogram of BPH urine. The four arrows indicate the retention times of four selected features. They were m/z 205.0685, 334.1126, 551.2124, and 407.2798 at retention times of 2.2, 10.3, 12.1, and 16.8 min, respectively. These four features are different in retention times, molecular weights, and ion intensities. They were used to evaluate the stability of the UPLC-FTMS. They were finally identified as D-mannitol, N2, N2-dimethylguanosine, phenylacetyl L-glutamine, and 7-ketodeoxycholic acid, respectively. The detected parameters are listed. (b) The PCA plot of BPH and Her urines. The symbol (Δ) denotes BPH, (x) Her, and (O) QC. The QC was made by mixing an equal amount of each urine sample.
Diagnostics 13 02270 g003
Figure 4. (a) The PCA plot of 27 PCa and 49 BPH urine samples. The symbol (Δ) denotes PCa, (O) BPH, and (Χ) QC. (b) The ratio distribution of 441 features. The distribution was fitted by the Gaussian distribution Curve; the three parameters were a: 170.99, b: 0.10, and x0: −0.0008. b was equal to one standard deviation of the curve, and x0 was the shift of the fitted curve to 0 (log1); (c) the reconstructed PCA plot made by four metabolites in PCa and BPH urines. The red spots denote PCa, and the blue spots denote BPH. (d) The ROC curve was calculated by the combination of four metabolites.
Figure 4. (a) The PCA plot of 27 PCa and 49 BPH urine samples. The symbol (Δ) denotes PCa, (O) BPH, and (Χ) QC. (b) The ratio distribution of 441 features. The distribution was fitted by the Gaussian distribution Curve; the three parameters were a: 170.99, b: 0.10, and x0: −0.0008. b was equal to one standard deviation of the curve, and x0 was the shift of the fitted curve to 0 (log1); (c) the reconstructed PCA plot made by four metabolites in PCa and BPH urines. The red spots denote PCa, and the blue spots denote BPH. (d) The ROC curve was calculated by the combination of four metabolites.
Diagnostics 13 02270 g004
Figure 5. (a) The product ion spectrum of m/z 334.1 ion. (b) The product ion spectrum of N2, N2-dimethylganosine.
Figure 5. (a) The product ion spectrum of m/z 334.1 ion. (b) The product ion spectrum of N2, N2-dimethylganosine.
Diagnostics 13 02270 g005
Figure 6. (a) The product ion spectrum of m/z 459.3 ion; (b) the interpretations of the fragmentation of m/z 459.3 ion and the fragmentation of the m/z 367.1 product ion. The chirality of molecular structure was not indicated. (c) The MS3 spectrum of m/z 459.3 to 367 ions. The m/z 459.3 ion was tentatively assigned as the Vitamin D3 derivative, but the side chain cannot be determined by MSn experiments.
Figure 6. (a) The product ion spectrum of m/z 459.3 ion; (b) the interpretations of the fragmentation of m/z 459.3 ion and the fragmentation of the m/z 367.1 product ion. The chirality of molecular structure was not indicated. (c) The MS3 spectrum of m/z 459.3 to 367 ions. The m/z 459.3 ion was tentatively assigned as the Vitamin D3 derivative, but the side chain cannot be determined by MSn experiments.
Diagnostics 13 02270 g006
Table 1. The 10 metabolites identified in the PCa and BPH urines. The features were first selected with fold changes over 1.2 or less than 0.8, and the ratio criteria were determined by the ratio distribution plot. After applying the ROC calculation on the selected features, they were filtered with the medium differential power with an AUC value of 0.65, and then these features were used to search against the METLIN and HMDB databases. After the web searching, 18 out of 24 features can be matched to compounds collected in METLIN and/or HMDB databases. Among the 18 matched features, 11 features (10 metabolites) can be interpreted by their product ion spectra. #: metabolites; a: the number of the product ion spectrum of the metabolite in Supplementary Table S3.
Table 1. The 10 metabolites identified in the PCa and BPH urines. The features were first selected with fold changes over 1.2 or less than 0.8, and the ratio criteria were determined by the ratio distribution plot. After applying the ROC calculation on the selected features, they were filtered with the medium differential power with an AUC value of 0.65, and then these features were used to search against the METLIN and HMDB databases. After the web searching, 18 out of 24 features can be matched to compounds collected in METLIN and/or HMDB databases. Among the 18 matched features, 11 features (10 metabolites) can be interpreted by their product ion spectra. #: metabolites; a: the number of the product ion spectrum of the metabolite in Supplementary Table S3.
#RT
(min)
Exact MassAdduct IonError (ppm)Metlin IDHMDB IDNamea Spectrum
115.4315.2365[M + H ]+14.9xHMDB13121 7-Dehydropregnenolone1
211.7318.1915[M + CAN + H]+3.166,195xp-Coumaroylagmatine2
310.7207.1107
185.1286
[M + CAN + Na] +
[M + CAN + H]+
−1.2
0.9
6589x1-Aminocyclohexanecarboxylic acid3
4
413.1344.2436[M + NH4]+−1.636,202x2,3-dinor-11b-PGF2α5
51.9441.2987[M + Na]+2.841,996
41,998
x
x
(24R)-1α,24-dihydroxy-22-oxacholecalciferol
(24S)-1α,24-dihydroxy-22-oxacholecalciferol
6
610.4334.1135[M + Na]+3.87086HMDB04824N2,N2-Dimethylguanosine7
712.5187.0225[M + Na]+0.744,770xLumazine8
814.9679.4108[M + 2ACN + H ]+4.858,221HMDB04157L-Urobilinogen9
912.7459.3336[M + H − H2O]+−1.342,352
42,353
42,354
42,355
x
x
x
x
1α,25-dihydroxy-24a-homo-26,27-dimethyl-22-thiacholecalciferol
1α,25-dihydroxy-26,27-dimethyl-24a-homo-22-thia-20-epicholecalciferol
1α,25-dihydroxy-24a,24b,24c-trihomo-22-thiacholecalciferol
1α,25-dihydroxy-24a,24b,24c-trihomo-22-thia-20-epicholecalciferol
10
1012.5227.0599[M + Na]+16.923,907xAla Asp11
Table 2. The statistics of the 10 metabolites in Table 1. The most abundant ion, the sodium adduct ion of 1-Aminocyclohexanecarboxylic acid, was used for the ratio calculation to achieve the best sensitivity.
Table 2. The statistics of the 10 metabolites in Table 1. The most abundant ion, the sodium adduct ion of 1-Aminocyclohexanecarboxylic acid, was used for the ratio calculation to achieve the best sensitivity.
#NameRatio (BPH/PCa)pROC Curve
AUCSD95% CI
17-Dehydropregnenolone0.52 0.02070.6850.0670.5526 To 0.8171
2p-Coumaroylagmatine0.620.06620.6770.0650.5488 To 0.8042
31-Aminocyclohexanecarboxylic acid0.660.007730.7190.0660.5902 To 0.8482
42,3-dinor-11b-PGF2α0.740.121470.6710.0660.5413 To 0.8011
5(24R)-1α,24-dihydroxy-22-oxacholecalciferol
(24S)-1α,24-dihydroxy-22-oxacholecalciferol
0.760.003370.6920.0660.5636 To 0.8212
6N2,N2-Dimethylguanosine0.770.031430.6970.0670.5657 To 0.8281
7Lumazine1.540.055920.6520.0680.5186 To 0.7852
8L-Urobilinogen2.530.000060.7160.0570.6034 To 0.8282
91α,25-dihydroxy-24a-homo-26,27-dimethyl-22-thiacholecalciferol
1α,25-dihydroxy-26,27-dimethyl-24a-homo-22-thia-20-epicholecalciferol
1α,25-dihydroxy-24a,24b,24c-trihomo-22-thiacholecalciferol1α,25-dihydroxy-24a,24b,24c-trihomo-22-thia-20-epicholecalciferol
2.90.009830.6580.0610.5373 To 0.7779
10Ala Asp3.960.023710.6970.0590.5810 To 0.8121
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, C.-L.; Chen, Y.-T.; Liao, W.-Y.; Chang, Y.-S.; Yu, J.-S.; Juo, B.-R. Urinary Metabolomic Analysis of Prostate Cancer by UPLC-FTMS and UPLC-Ion Trap MSn. Diagnostics 2023, 13, 2270. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics13132270

AMA Style

Chen C-L, Chen Y-T, Liao W-Y, Chang Y-S, Yu J-S, Juo B-R. Urinary Metabolomic Analysis of Prostate Cancer by UPLC-FTMS and UPLC-Ion Trap MSn. Diagnostics. 2023; 13(13):2270. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics13132270

Chicago/Turabian Style

Chen, Chien-Lun, Yi-Ting Chen, Wen-Yu Liao, Yu-Sun Chang, Jau-Song Yu, and Bao-Rong Juo. 2023. "Urinary Metabolomic Analysis of Prostate Cancer by UPLC-FTMS and UPLC-Ion Trap MSn" Diagnostics 13, no. 13: 2270. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics13132270

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