International Journal of Analytical Chemistry

International Journal of Analytical Chemistry / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 8839215 |

Mengchan Fang, Fan Liu, Lingling Huang, Liqing Wu, Lan Guo, Yiqun Wan, "A Urine Metabonomics Study of Rat Bladder Cancer by Combining Gas Chromatography-Mass Spectrometry with Random Forest Algorithm", International Journal of Analytical Chemistry, vol. 2020, Article ID 8839215, 9 pages, 2020.

A Urine Metabonomics Study of Rat Bladder Cancer by Combining Gas Chromatography-Mass Spectrometry with Random Forest Algorithm

Academic Editor: David M. Lubman
Received15 Jun 2020
Revised06 Sep 2020
Accepted09 Sep 2020
Published21 Sep 2020


A urine metabolomics study based on gas chromatography-mass spectrometry (GC-MS) and multivariate statistical analysis was applied to distinguish rat bladder cancer. Urine samples with different stages were collected from animal models, i.e., the early stage, medium stage, and advanced stage of the bladder cancer model group and healthy group. After resolving urea with urease, the urine samples were extracted with methanol and, then, derived with N, O-Bis(trimethylsilyl) trifluoroacetamide and trimethylchlorosilane (BSTFA + TMCS, 99 : 1, v/v), before analyzed by GC-MS. Three classification models, i.e., healthy control vs. early- and middle-stage groups, healthy control vs. advanced-stage group, and early- and middle-stage groups vs. advanced-stage group, were established to analyze these experimental data by using Random Forests (RF) algorithm, respectively. The classification results showed that combining random forest algorithm with metabolites characters, the differences caused by the progress of disease could be effectively exhibited. Our results showed that glyceric acid, 2, 3-dihydroxybutanoic acid, N-(oxohexyl)-glycine, and D-turanose had higher contributions in classification of different groups. The pathway analysis results showed that these metabolites had relationships with starch and sucrose, glycine, serine, threonine, and galactose metabolism. Our study results suggested that urine metabolomics was an effective approach for disease diagnosis.

1. Introduction

Bladder cancer (BC) is a common malignant tumor disease of the urinary tract, and its incidence and mortality have always occupied the first place in the urinary reproductive system tumors. Due to the easy relapse characteristic, BC has been the focus of researchers to search tumor markers for the early diagnosis and postoperative evaluation to improve the survival rate of bladder cancer patients [1].

Metabonomics has been widely used in the research of diseases diagnosis [2], pharmacological [3] and toxicological mechanisms [4], plant and microorganism metabolism [5, 6], and so on [7]. Two analysis methods were widely used in metabonomics. One is the metabonomics technology basing on nuclear magnetic resonance (NMR) [810], and the other is chromatography-mass spectrometry [1113]. NMR technology has the advantages of fast analysis and simple sample preparation, while it also possesses the disadvantage of low sensitivity [14]. Chromatography-mass spectrometry technology mainly contains gas chromatography-mass spectrometry (GC-MS) [15, 16], liquid chromatography-mass spectrometry (LC-MS) [1719], and capillary electrophoresis-mass spectrometry (CE-MS) [2022]. Among these technologies, GC-MS has been widely used in metabonomics studies owing to its high sensitivity, strong analysis ability, and possessing more mature commercial mass spectrum library [23]. The samples used in metabonomics commonly are biological fluids, such as urine [24, 25], serum [26, 27], interstitial fluid [28], and cerebrospinal fluid [29]. Due to the characteristic of weak volatility and poor thermal stability, many analytes in biological fluids, such as amino acids and organic acids, must be derivatized before GC-MS analysis. The development of modern derivatization technology has greatly promoted the application of GC-MS in metabonomics [30, 31].

Some studies were proposed to discover biomarkers for BC [3236] in recent years. Pasikanti et al. [37] have developed a noninvasive method for the diagnosis and surveillance of BC progression by using the GC × GC-MS method. Peng et al. [38] reported a chemical isotope labeling LC-MS metabolomics method based on the universal metabolome-standard method with low CV for all quantified metabolites. This method was used to screen potential biomarkers from urine samples for bladder cancer diagnosis. However, the current study for diagnosis of BC is mostly focusing on high-grade tumors. Discovering biomarkers which could characterize different stages of BC is of more useful in diagnosis and prognostics of bladder cancer [39]. Alberice et al. [40] followed up 48 urothelial bladder cancer patients by using urine metabolomics, and 27 metabolites were highlighted between different BC stages and recurrence. However, the BC patients following up is difficult, especially the early-stage BC samples collection is of a hard work since most BC patients are in high stage when diagnosed.

Therefore, in this study, the rat bladder cancer model was established, and the urine metabonomics was studied with GC-MS technology. The rat urine samples of the advanced stage, medium stage, and early stage of the bladder cancer model group and healthy group had been detected to establish the bladder cancer urine metabolic fingerprint, and then, the experimental data were analyzed with algorithm of random forests [41, 42] and used for the preliminary exploration of the tumor markers of BC.

2. Experiment

2.1. Standards and Reagents

Twenty amino acid and nine carbohydrate standards, including isoleucine, lysine hydrochloride, cysteine, glycine, tryptophan, glutamic acid, valine, threonine, proline, glutamate, leucine, methionine, phenylalanine, cystine, aspartic acid, asparagine, histidine, serine, 4-hydroxy proline, arginine, galactose, fructose, sucrose, glucose, mannose, ribose, xylose, arabinose, and glucuronic acid, were purchased from Sigma-Aldrich (St. Louis, MO, USA). Urease was of analytical grade and purchased from Sigma-Aldrich (St. Louis, MO, USA). Methanol, acetonitrile, and acetone were of chromatographical grade and purchased from Merck Company (Darmstadt, Germany); n-docosane, pyridine, and n-heptane were of analytical grade and purchased from Sinopharm Chemical Reagent Co. Ltd. (Shanghai, China). N,O-Bis(trimethylsilyl) trifluoroacetamide + trimethylchlorosilane (BSTFA + TMCS, 99 : 1, v/v) and methoxyamine hydrochloride were purchased from Supelco (Bellefonte, PA, USA).

2.2. Rat Bladder Cancer Model

All animal experiments were conducted according to the institutional guidelines of Medical College of Nanchang University. One hundred and eighty male SD rats (6 weeks old) weighing 150–160 g were purchased from Hunan Experimental Animal Co. LTD. (Hunan, China), with the license number of SCXK (Hunan) 2009–0004 and the qualified number of HNASLKJ20102985. The rats were randomly assigned to two groups, forty-five rats were in the control group and one hundred and thirty-five rats were in the model group. Bladder tumors were induced by adding 0.05% BBN (N-butyl-N-(4-hydroxybutyl) nitrosamine) (Tokyo, Japan) to freely available drinking water after a week of adaptive period in the experimental animal room. BBN was continuously administered for 35 weeks in the experiment.

2.3. Histologic Examination

After induced by BBN, three rats of the BC model group were randomly selected for histological examination every five weeks. The rats were anesthetized by intraperitoneal injection of ketamine (0.6 mL/50 g), and the rats were killed under deep anesthesia. After death, rats were catheterized and 0.2 mL buffered formalin was instilled into the bladder. The urethras were, then, ligated, and the bladders were removed. Fixed specimens were embedded in paraffin, and 4–5 μm thick horizontal slices in each rat were prepared at 2 mm intervals followed by routine hematoxylin and eosin staining. Then, the prepared pathological sections were detected with optical microscopy using ordinary white light by the same pathologists.

2.4. Sample Preparation

The rat urine samples of the control group and model group were collected with a metabolism cage for 24 h every five weeks. The urine samples were centrifuged immediately for 20 min at 4000 r/min to remove protein with a TDL-5-A low-speed tabletop centrifuge (Shanghai Anting Scientific Instrument Factory, China). Then, 150 μL supernatant liquid was mixed with 100 μL of urease solution (2 mg/mL) in a 1 mL centrifuge tube and heated at 37°C in a DGG-9140BD constant temperature oven (Shanghai Senxin Experimental Instrument Co., Ltd., China) for 30 min to decompose urea. After adding 800 μL methanol, the mixture was homogenized for 1 min at 1800 r/min with a MS2 mini shaker (Guangzhou Yike Lab Technology LTM Co., China), followed by ultrasonic processing in an ice bath for 10 min with a KQ-100DE ultrasonic cleaner (Kunshan Ultrasonic Instruments Co., Ltd., China) and centrifuging for 10 min at 12000 r/min with a TDL-16G high-speed tabletop centrifuge (Shanghai Anting Scientific Instrument Factory, China). Then, 500 μl supernatant liquid was transferred into a 2 mL centrifuge tube and evaporated to dryness under a gentle nitrogen stream, and then, 75 μL methoxyamine hydrochloride solution (20 mg/mL in pyridine) was added into the tube to react for 1 h at 70°C. After the reaction, the mixture was cooled to room temperature and reacted with 75 μl BSTFA + TMCS (99 : 1, v/v) for 1 h at room temperature to form trimethylsilyl (TMS) derivatives. Finally, the reaction was terminated by the addition of 150 μL n-heptane (containing 0.1 g/L n-docosane, which was used as an internal standard substance), and the products of the derivative reaction were analyzed by GC-MS.

2.5. GC-MS Analysis

GC-MS analysis was carried out using an Agilent 6890N Gas Chromatograph (Agilent Technologies, Palo Alto, CA, USA) integrated with an Agilent 7683 series autosampler and a 5973 I mass selective detector (MSD). The analytes were separated on a 30 m × 0.25 mm i.d. × 0.25 μm film thickness DB-5MS fused-silica capillary column (Agilent Technologies). The injector was set at 270°C, and the carrier gas was UHP helium at a flow rate of 1.0 mL/min. The samples were injected in a splitless mode. The oven temperature was initially at 70°C for 5 min; increased at a rate of 20°C/min up to 160°C, held for 4 min; increased at a rate of 10°C/min up to 300°C; and held for 1.5 min. The ion source, quadrupole, and transfer line temperature were set at 230, 150°C, and 280°C, respectively. The detector was operated at 70 ev in an electron impact (EI) mode with full scan (60∼600 amu). The solvent delay and injection volume were set as 5 min and 2 μL, respectively. All the data was collected and analyzed with MSD ChemStation D.01.02 software (Agilent Technologies), with the NIST02 mass database.

3. Results and Discussion

3.1. Optimization of the Urine Sample Preparation Method
3.1.1. Selection of Extraction Solvent

There are many kinds of endogenous metabolites with high polarity difference and wide concentration range in rat urine. To obtain more information of endogenous metabolites, the extraction efficiencies of different organic solvents were investigated. It was found that methanol and acetonitrile performed better than others. Finally, methanol was chosen as the extraction solvent in this study because of the lower toxicity of methanol than acetonitrile.

3.1.2. Optimization of the Urea Decomposition Conditions

The urine should be treated with urease to decompose the urea before GC-MS analysis, since the urea with high concentration in urine may affect other compounds analysis. Effects of several factors including the dosage of urease, decomposition temperature, and decomposition time were investigated. 50 µL, 100 µL, and 200 µL urease solution (2 mg/mL) was added to a 150 µL urine sample and reacted at 60°C, 37°C, and 20°C for 15 min, 30 min, and 60 min. Then, its decomposition effect was examined. The results are listed in Supplementary Materials. Finally, the optimum conditions for urea decomposing in the 150 μL urine sample were set as dosage of urease solution (2 mg/mL), 100 μL; decomposition temperature, 37°C; and decomposition time, 30 min.

3.1.3. Optimization of the Multistep Temperature Program

To obtain better separation and get more information of endogenous metabolites, four multistep temperature programs as shown in below were carried out on a DB-5MS capillary column in this study. Considering the number of chromatographic peaks and their separation efficiency, the third temperature-rising program was selected as the experimental condition. The total ion chromatogram (TIC) of the actual sample is shown in Figure 1, and there were more than 40 peaks obtained in 29 min with good separation efficiency.(1)The temperature was initially at 85°C for 5 min, increased at a rate of 10°C/min up to 300°C, and held for 10 min(2)The temperature was initially at 100°C for 3 min, increased at a rate of 8°C/min, and held for 2 min(3)The temperature was initially at 70°C for 5 min, increased at a rate of 20°C/min up to 160°C, held for 4 min, increased at a rate of 10°C/min up to 300°C, and held for 1.5 min(4)The temperature was initially at 85°C for 5 min, increased at a rate of 8°C/min up to 205°C, held for 5 min, increased at a rate of 8°C/min up to 300°C, and held for 5 min

3.2. Urine Sample Analysis

The different stages of bladder cancer rats were confirmed by the histopathology analysis. Because the rat bladder of fifteenth, twenty-fifth, and thirty-fifth week after BBN inducing displayed the typical characteristics of early, medium, and advanced stage tumors (shown in Figure 2), the urine samples collected in these three weeks were chosen as early-, medium-, and advanced-stage samples, respectively. The early, medium, and advanced stage of the bladder cancer model group and healthy group include 45 samples each, and all the samples were analyzed by GC-MS under the optimized conditions with 41 common peaks obtained. To identify the complex metabolites in urine, the metabolites were divided into three categories: amino acids, carbohydrates, and fatty acids. Because the NIST database has rich information of fatty acid derivatives but poor information of amino acid derivatives and carbohydrate derivatives, the fatty acid derivatives were identified directly with the NIST database, while the amino acids and carbohydrate derivatives were identified by standard compounds. The quantitative results of metabolites were given as the peak area ratio of the analyte to the internal standard of n-docosane. Finally, all the 41 common peaks were identified successfully, and the qualitative and quantitative results are listed in Table 1.

No. (min)Endogenous metabolitesHealthy group (Ax/Ai, n = 45)Early stage (Ax/Ai, n = 45)Medium stage (Ax/Ai, n = 45)Advanced stage (Ax/Ai, n = 45)

15.922Aminoethane0.2456 ± 0.07050.1782 ± 0.03520.1905 ± 0.05670.1958 ± 0.0551
26.593Ethylene glycol0.0182 ± 0.00200.0373 ± 0.01070.0530 ± 0.04280.0746 ± 0.0626
36.840N, N-diethylacetamide0.0657 ± 0.00870.0584 ± 0.00990.0476 ± 0.02020.0557 ± 0.0107
47.716Lactic acid0.0872 ± 0.03740.0765 ± 0.029340.0952 ± 0.05920.1482 ± 0.2155
57.934Acetic acid0.0856 ± 0.03330.0341 ± 0.02540.0229 ± 0.01400.0412 ± 0.0203
610.010Phosphate2.1278 ± 0.91732.7785 ± 1.40371.4730 ± 0.73811.3767 ± 0.9361
710.200l-threonine0.0173 ± 0.00980.0052 ± 0.00290.0108 ± 0.00680.0096 ± 0.0065
810.297Phenylacetic acid0.0047 ± 0.00230.0183 ± 0.00540.0159 ± 0.01030.0147 ± 0.0097
910.382Succinic acid0.0311 ± 0.01290.0117 ± 0.00840.0098 ± 0.00310.0119 ± 0.0086
1010.4471, 2-Hydroquinone0.0120 ± 0.00720.0103 ± 0.00730.0078 ± 0.00470.0067 ± 0.0039
1110.503Glyceric acid0.0961 ± 0.02660.0250 ± 0.01200.0400 ± 0.02320.0183 ± 0.0087
1210.723(R, R)-2, 3-dihydroxybutanoic acid0.0167 ± 0.00530.0045 ± 0.00130.0037 ± 0.00140.0053 ± 0.0029
1311.3572, 4-Dihyoxybutanoic acid0.0147 ± 0.00510.0131 ± 0.00410.0155 ± 0.00800.0166 ± 0.0047
1411.583(R, S)-3, 4-dihydroxybutanoic acid0.0304 ± 0.00980.0161 ± 0.00530.0132 ± 0.00640.0178 ± 0.0107
1511.797N-(1-oxobutyl)- glycine0.0653 ± 0.02440.0191 ± 0.00960.0319 ± 0.01860.0274 ± 0.0151
1612.341Isovaleroglycine0.0356 ± 0.01340.0134 ± 0.00410.0160 ± 0.00790.0107 ± 0.0073
1712.483D-threitol0.0714 ± 0.02730.0665 ± 0.02600.0290 ± 0.01300.0251 ± 0.0151
1812.645N-crotonyl glycine0.0240 ± 0.01460.328 ± 0.00740.0207 ± 0.01290.0148 ± 0.0099
1912.973, 13.2032, 3, 4-Trihydroxybutyrate0.1276 ± 0.01620.0475 ± 0.02440.0631 ± 0.03430.0412 ± 0.0250
2014.530N-(1-oxohexyl)-glycine0.0960 ± 0.03190.0112 ± 0.0390.0421 ± 0.02730.0232 ± 0.0081
2114.5803-Hydroxyphenylacetic acid0.0326 ± 0.01000.0122 ± 0.09100.0140 ± 0.00810.0134 ± 0.0088
2214.713D-xylose0.0408 ± 0.01500.0225 ± 0.02060.0182 ± 0.00440.0193 ± 0.0053
2314.823, 15.057D-ribose0.0926 ± 0.03700.0682 ± 0.03400.0252 ± 0.01420.0250 ± 0.0179
2415.509, 15.733Arabitol0.0287 ± 0.01640.0252 ± 0.01680.0283 ± 0.01790.0278 ± 0.0215
2516.0236-Deoxy-D-galactose0.0336 ± 0.00830.0370 ± 0.00640.0177 ± 0.01000.0149 ± 0.0104
2616.087Mannonic acid0.0505 ± 0.01770.0419 ± 0.01680.0211 ± 0.01430.0168 ± 0.0138
2716.200cis-Aconitic acid0.0535 ± 0.02880.0524 ± 0.02090.0105 ± 0.00790.0168 ± 0.0147
2816.357Phosphoric acid0.0414 ± 0.02020.0463 ± 0.01770.0230 ± 0.01410.0212 ± 0.0168
2917.177Isocitric acid0.0348 ± 0.01210.0410 ± 0.01570.0140 ± 0.00930.0248 ± 0.0138
3017.563Hippuric acid0.0470 ± 0.01260.0201 ± 0.01110.0180 ± 0.00740.0156 ± 0.0096
3117.850, 17.960D-fructose0.0512 ± 0.02860.0728 ± 0.01060.0371 ± 0.01450.0480 ± 0.0131
3218.087N-phenyl glycine0.0596 ± 0.02140.0363 ± 0.01940.0455 ± 0.02720.0389 ± 0.0287
3318.197, 18.147D-glucose0.3785 ± 0.16180.2254 ± 0.13280.1741 ± 0.06540.1859 ± 0.0736
3418.507Altronic acid0.0302 ± 0.00690.0164 ± 0.01060.0185 ± 0.01000.0102 ± 0.0074
3518.577, 18.650D-sorbitol0.0896 ± 0.02690.0623 ± 0.01710.0254 ± 0.01870.0300 ± 0.0275
3618.983, 19.533Galactonic acid0.0613 ± 0.02820.0387 ± 0.01860.0617 ± 0.03280.0441 ± 0.0351
3719.990Palmitic acid0.0084 ± 0.00090.0085 ± 0.00470.0067 ± 0.00170.0071 ± 0.0025
3820.403Myo-inositol0.0347 ± 0.02280.0089 ± 0.00330.0097 ± 0.00370.0134 ± 0.0129
3925.465D-turanose0.0216 ± 0.01380.0145 ± 0.01180.0197 ± 0.00900.0510 ± 0.0099
4025.653, 25.783D- (+)-lactose monohydrate1.0400 ± 0.33490.9997 ± 0.23850.7475 ± 0.23660.6559 ± 0.3286
4125.927Lactose0.0142 ± 0.00430.0231 ± 0.00390.0143 ± 0.00750.0190 ± 0.0163

No: the serial number of the common peak; : retention time. Ax/Ai: the ratio of the peak area of the analyte to that of the internal standard.
3.3. Experimental Data Analysis

The classifiers commonly used in metabolomics include partial least squares discriminant analysis (PLS-DA), support vector machine (SVM), and random forest (RF) [43]. PLS-DA can reduce the impact of multiple correlations between variables, but it is easy to overfit the data, and the selected biomarkers are not robust enough. SVM can solve small sample classification, high-dimensional data classification, and nonlinear problems; however, it is more difficult to train large-scale samples and deal with multiclassification problems. Random forests (RF) algorithm was first proposed by Breiman in 2001 [44] and widely used [45, 46] since it can distinguish the differences between different group samples effectively. RF is a supervised machine learning classifier, including a collection of classification and regression trees. It consists of many different decision trees, which are grown based on various guide samples. Each tree voted for the sample for classification, and RF chose the majority vote to determine the final classification result. It has good performance and has great advantages compared with other algorithms. RF can handle high-dimensional data (many feature data) and does not need to make feature selection. After training, random forest can screen out more important features. Compared with other classification models, its biased estimation of classification results is low, which makes random forests applicable to many research fields. RF was adopted to classify metabolites among the four groups of rat urines, and the obtained multidimensional scaling (MDS) figure is shown in Figure 3. As shown in Figure 3, the differences among different groups were obviously in the classification plot. The cancer groups could be effectively distinguished from the healthy groups; moreover, the advanced-stage groups could also be distinguished from the early and middle-stage groups. However, the sample points of the early stage and middle stage are somewhat overlapped. The results suggested that the metabolic pathways of rat bladder cancer were obviously different from that of healthy rats, which were similar in the early and middle stage, but in the advanced stage, the metabolic pathways had changed significantly due to tumor deterioration and excessive nutrient consumption.

Thus, we further established three classification models, i.e., healthy control vs. early- and middle-stage groups, healthy control vs. advanced-stage group, and early- and middle-stage groups vs. advanced stage group, respectively. The classification results for three models are listed in Table 2. It is obvious that the healthy control vs. advanced-stage group had the best classification accuracy, indicating the signification differences of metabolic features between the healthy control and advanced-stage group; and the healthy control vs. early- and middle-stage groups also had a classification accuracy of 96.06%. These results suggested that the proposed metabolomics approach can reflect the differences among different groups, with the progress of the disease.

ResultsSensitivity (%)Specificity (%)Accuracy (%)MCCAUC

Healthy control vs. early- and middle-stage group83.3398.8096.060.86090.9775
Healthy control vs. advanced-stage group97.2297.6097.540.91940.9868
Early- and middle-stage groups vs. advanced stage group80.5697.0194.090.79090.9546

MCC: Matthews correlation coefficient, AUC: area under curve.

During the classification model establishing, the importance of the metabolites was calculated. The metabolites with higher importance values in the classification have more contributions to clinical diagnosis of disease, which means these metabolites can be used as potential biomarkers for disease diagnosis, especially in the early diagnosis. The variable importance of metabolites for each classification model is shown in Figure 4. Glyceric acid, 2, 3-dihydroxybutanoic acid, N-(oxohexyl)-glycine, and D-turanose showed the highest variable importance (higher than 0.45, Figure 4), which are more likely to be useful markers for BC diagnosis.

Glyceric acid is an important intermediate in the lipid metabolism, which can be produced by the oxidation of fatty acids and the hydrolysis of phosphoglyceric acids (such as 2-phosphoglyceric acid, 3-phosphoglyceric acid, and 1, 3-bisphosphoglyceric acid). Phosphoglyceric acids are important intermediate products of the tricarboxylic acid cycle in the organisms and directly involved in the metabolism and transformation of energy, such as 1, 3-bisphosphoglyceric acid, which is a high-energy phosphate compound in vivo and can produce one molecule of ATP to the living body under the catalysis of the phosphoglycerate enzyme.

2, 3-Dihydroxybutanoic acid is related to the metabolic pathway of L-threonine and generated by the metabolites of L-threonine, which is a ketogenic amino acid, and its metabolites can directly join in the energy metabolism.

N-(oxohexyl)-glycine is one of the acyl amino acids, exactly acyl glycine, in organisms and usually produced in the metabolic process of fatty acids with very small quantity. Acyl glycine is usually produced under the catalysis of acyltransferase. The reaction is as follows: glycine + acyl-coenzyme A ⟶ acyl glycine + coenzyme A. Acetyl coenzyme A can directly provide a molecular dicarbon compound for the tricarboxylic acid cycle. Furthermore, the combination of oxaloacetate and acetyl coenzyme A is believed to be the initial step in the citric acid cycle. As a result, the abnormalities of the acyl glycine metabolism may affect the energy metabolism of cells and form the specific metabolic pathways of tumor.

All the abovementioned results indicated that the metabolic pathways of lipid and some amino acids changed significantly as the bladder tumor grew. Thus, based on these metabolites, the pathway analysis was implemented by using Metaboanalyte software. The pathway analysis results showed that starch and sucrose metabolism, glycine, serine, and threonine metabolism, and galactose metabolism have strong relationships with selected metabolites (shown in Figure 5).

4. Conclusions

Metabolomics is an effective approach to discover biomarkers by analyzing global changes in the metabolic profiles. To collect early-stage bladder cancer (BC) samples and follow-up the BC progress, the rat bladder cancer model was established by BBN inducing in this study. The metabolites in rat urine were detected with GC-MS and analyzed with random forests algorithm to distinguish the early, middle, and advanced stage of the bladder cancer group and healthy group. The results showed that urinary levels of some metabolites had a significant difference between the cancer group and the healthy group and advanced-stage group and the other two stage groups, which suggested that the growth of bladder tumor might result in the abnormality of the metabolism of lipids and some amino acids. Furthermore, glyceric acid, 2, 3-dihydroxybutanoic acid, N-(oxohexyl)-glycine, and D-turanose with the highest variable importance might be the potential markers of bladder cancer, and their metabolic pathways were studied. Yet, the data reported here are preliminary and need to be confirmed by large scale of samples. Further studies should be required to value the significance of the four compounds as the potential marker in human urine for bladder cancer.

Data Availability

The figure data and related data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest.


This work was financially supported by the Natural Science Foundation of China (No. 21465017) and the Science and Technology Innovation Platform Project of Jiangxi Province (No. 20192BCD40001). The authors gratefully acknowledge the financial supports.

Supplementary Materials

Figure S1: the total ion chromatograms (TIC) of samples with different extraction solvents. Figure S2: the total ion chromatograms (TIC) of samples with different dosages of urease. Figure S3: the total ion chromatograms (TIC) of samples with different decomposition temperature. Figure S4: the total ion chromatograms (TIC) of samples with different decomposition time. (Supplementary Materials)


  1. K. K. Pasikanti, K. Esuvaranathan, P. C. Ho et al., “Noninvasive urinary metabonomic diagnosis of human bladder cancer,” Journal of Proteome Research, vol. 9, no. 6, pp. 2988–2995, 2010. View at: Publisher Site | Google Scholar
  2. N. Mahendran, G.-X. xie, and W. Jia, “Metabonomics of human colorectal cancer: new approaches for early diagnosis and biomarker discovery,” Journal of Proteome Research, vol. 13, no. 9, pp. 3857–3870, 2014. View at: Google Scholar
  3. D. G. Robertson, M. D. Reily, and J. D. Baker, “Metabonomics in pharmaceutical discovery and development,” Journal of Proteome Research, vol. 6, no. 2, pp. 526–539, 2007. View at: Publisher Site | Google Scholar
  4. W. Xie, W. Zhang, J. Ren et al., “Metabonomics indicates inhibition of fatty acid synthesis, β-oxidation, and tricarboxylic acid cycle in triclocarban-induced cardiac metabolic alterations in male mice,” Journal of Agricultural and Food Chemistry, vol. 66, no. 6, pp. 1533–1542, 2018. View at: Publisher Site | Google Scholar
  5. C. Li and S. Sturm, “Analytical aspects of plant metabolite profiling platforms: current standings and future aims,” Journal of Proteome Research, vol. 6, no. 2, pp. 480–497, 2007. View at: Google Scholar
  6. K. Eloh, N. Sasanelli, and A. Maxia, “Untargeted metabolomics of tomato plants after root-knot nematode infestation,” Journal of Agricultural and Food Chemistry, vol. 64, no. 29, pp. 5963–5968, 2016. View at: Publisher Site | Google Scholar
  7. J. T. Caboni, J. K. Nicholson, P. M. Schofield, D. J. Grainger, and E. Holmes, “Application of chemometrics to 1H NMR spectroscopic data to investigate a relationship between human serum metabolic profiles and hypertension,” The Analyst, vol. 128, no. 1, pp. 32–36, 2003. View at: Publisher Site | Google Scholar
  8. A.-H. M. Grainger, J. S. Merzaban, and S. M. H. Serrai, “Theory and applications of NMR-based metabolomics in human disease diagnosis,” Applications of NMR Spectroscopy, vol. 1, pp. 93–130, 2015. View at: Publisher Site | Google Scholar
  9. H. M. Sonawat and S. Sharma, “Host responses in malaria disease evaluated through nuclear magnetic resonance-based metabonomics,” Clinics in Laboratory Medicine, vol. 32, no. 2, pp. 129–142, 2012. View at: Publisher Site | Google Scholar
  10. M. Cao, L. Zhao, H. Chen, and D. Lin, “NMR-based metabolomic analysis of human bladder cancer,” Analytical Sciences, vol. 28, no. 5, pp. 451–456, 2012. View at: Publisher Site | Google Scholar
  11. Y. Xue, K. K. Pasikanti, and C. W. Yap, “Comparative evaluation of software for retention time alignment of gas chromatography/time-of-flight mass spectrometry-based metabonomic data,” Journal of Chromatography A, vol. 1217, no. 52, pp. 8308–8316, 2010. View at: Publisher Site | Google Scholar
  12. J. M. Chan, D. Waterman, A. M. Przyborowska, P. D. Fraser, and P. M. Bramley, “Chemical derivatization and mass spectral libraries in metabolic profiling by GC/MS and LC/MS/MS,” Journal of Experimental Botany, vol. 56, no. 410, pp. 219–243, 2004. View at: Publisher Site | Google Scholar
  13. E. M. Patel and I. D. Wilson, “Analytical strategies in metabonomics,” Journal of Proteome Research, vol. 6, no. 2, pp. 443–458, 2007. View at: Publisher Site | Google Scholar
  14. M. B. Lauridsen, H. Bliddal, R. Christensen et al., “1H NMR spectroscopy-based interventional metabolic phenotyping: a cohort study of rheumatoid arthritis patients,” Journal of Proteome Research, vol. 9, no. 9, pp. 4545–4553, 2010. View at: Publisher Site | Google Scholar
  15. S. Danneskiold-Samsøe, E. Johansson, L. Sjöström et al., “Visualization of GC/TOF-MS-based metabolomics data for identification of biochemically interesting compounds using OPLS class models,” Analytical Chemistry, vol. 80, no. 1, pp. 115–122, 2008. View at: Publisher Site | Google Scholar
  16. S. M. M. Mellerowicz, C. Nebel, N. C. Nielsen et al., “A GC-MS-based metabonomic investigation of blood serum from irritable bowel syndrome patients undergoing intervention with acidified milk products,” European Food Research and Technology, vol. 233, no. 6, pp. 1013–1021, 2011. View at: Publisher Site | Google Scholar
  17. S. Andersen, K. Scholz, M. Donegan, J. Wingate, and W. Völkel, “Metabonomics and biomarker discovery: LC−MS metabolic profiling and constant neutral loss scanning combined with multivariate data analysis for mercapturic acid analysis,” Analytical Chemistry, vol. 78, no. 4, pp. 1296–1305, 2006. View at: Publisher Site | Google Scholar
  18. M. Burton, J. C. Serrano, N. Ortega et al., “Multicompartmental LC-Q-TOF-based metabonomics as an exploratory tool to identify novel pathways affected by polyphenol-rich diets in mice,” Journal of Proteome Research, vol. 10, no. 8, pp. 3501–3512, 2011. View at: Google Scholar
  19. A. M. Portero-Otín, L. G. Mullin, I. D. Wilson et al., “Development of a rapid profiling method for the analysis of polar analytes in urine using HILIC–MS and ion mobility enabled HILIC–MS,” Metabolomics, vol. 15, no. 2, 2019. View at: Publisher Site | Google Scholar
  20. M. Cauchi, C. M. Weber, B. J. Bolt et al., “Evaluation of gas chromatography mass spectrometry and pattern recognition for the identification of bladder cancer from urine headspace,” Analytical Methods, vol. 8, no. 20, pp. 4037–4046, 2016. View at: Publisher Site | Google Scholar
  21. R. Ramautar, G. W. Somsen, and G. J. de Jong, “CE-MS in metabolomics,” Electrophoresis, vol. 30, no. 1, pp. 276–291, 2009. View at: Publisher Site | Google Scholar
  22. A. Hirayama, K. Kami, M. Sugimoto et al., “Quantitative metabolome profiling of colon and stomach cancer microenvironment by capillary electrophoresis time-of-flight mass spectrometry,” Cancer Research, vol. 69, no. 11, pp. 4918–4925, 2009. View at: Publisher Site | Google Scholar
  23. K. K. Sugawara, P. C. Ho, and E. C. Y. Chan, “Gas chromatography/mass spectrometry in metabolic profiling of biological fluids,” Journal of Chromatography B, vol. 871, no. 2, pp. 202–211, 2008. View at: Publisher Site | Google Scholar
  24. J. L. McClay, D. E. Adkins, N. G. Isern et al., “1H nuclear magnetic resonance metabolomics analysis identifies novel urinary biomarkers for lung function,” Journal of Proteome Research, vol. 9, no. 6, pp. 3083–3090, 2010. View at: Publisher Site | Google Scholar
  25. M. O’Connell, S. H. Hansen, and J. W. Jaroszewski, “Human urine as test material in1H NMR-based metabonomics: recommendations for sample preparation and storage,” Analytical Chemistry, vol. 79, no. 3, pp. 1181–1186, 2007. View at: Publisher Site | Google Scholar
  26. Y. Cornett, R. Song, Z. Zhang et al., “The development of plasma pseudotargeted GC-MS metabolic profiling and its application in bladder cancer,” Analytical and Bioanalytical Chemistry, vol. 408, no. 24, pp. 6741–6749, 2016. View at: Publisher Site | Google Scholar
  27. S. J. Lu, I. Tavazzi, V. Parisod, S. Kochhar, and P. A. Guy, “Investigation of human blood plasma sample preparation for performing metabolomics using ultrahigh performance liquid chromatography/mass spectrometry,” Analytical Chemistry, vol. 81, no. 9, pp. 3285–3296, 2009. View at: Publisher Site | Google Scholar
  28. T. Rezzi, D. A. Broszczak, J. A. Broadbent, H. Lu, and T. J. Parker, “The biochemistry of blister fluid from pediatric burn injuries: proteomics and metabolomics aspects,” Expert Review of Proteomics, vol. 13, no. 1, pp. 35–53, 2016. View at: Publisher Site | Google Scholar
  29. Y. Cuttle, Y. Wang, M. Li et al., “1H NMR-based metabolomics exploring biomarkers in rat cerebrospinal fluid after cerebral ischemia/reperfusion,” Molecular BioSystems, vol. 9, no. 3, pp. 431–439, 2013. View at: Publisher Site | Google Scholar
  30. S. H. Xu, H. M. Woo, B. H. Jung et al., “Metabolomic approach to evaluate the toxicological effects of nonylphenol with rat urine,” Analytical Chemistry, vol. 79, no. 16, pp. 6102–6110, 2007. View at: Publisher Site | Google Scholar
  31. J. L. Lee, “Artifacts in trimethylsilyl derivatization reactions and ways to avoid them,” Journal of Chromatography A, vol. 844, no. 1-2, pp. 1–22, 1999. View at: Publisher Site | Google Scholar
  32. H. Shi, X. Li, Q. Zhang, and X. Zhang, “Discovery of urine biomarkers for bladder cancer via global metabolomics,” Biomarkers, vol. 21, no. 7, pp. 578–588, 2016. View at: Publisher Site | Google Scholar
  33. J. Yang, G. Xie, R. Che et al., “Distinct metabolic signature of human bladder cancer cells carrying an impaired Fanconi anemia tumor-suppressor signaling pathway,” Journal of Proteome Research, vol. 15, no. 4, pp. 1333–1341, 2016. View at: Publisher Site | Google Scholar
  34. Y. Su, X. Yang, X. Deng et al., “Metabolomics in bladder cancer: a systematic review,” International Journal of Clinical and Experimental Medicine, vol. 8, no. 7, pp. 11052–11063, 2015. View at: Google Scholar
  35. Y. Lu, K. Esuvaranathan, and K. K. Pasikanti, “Metabonomic diagnosis of bladder cancer: what are the advantages and considerations?” Bioanalysis, vol. 7, no. 21, pp. 2727–2729, 2015. View at: Publisher Site | Google Scholar
  36. E. C. Y. Yong Chan, K. K. Pasikanti, Y. Hong et al., “Metabonomic profiling of bladder cancer,” Journal of Proteome Research, vol. 14, no. 2, pp. 587–602, 2014. View at: Publisher Site | Google Scholar
  37. K. K. Ho, K. Esuvaranathan, Y. Hong et al., “Urinary metabotyping of bladder cancer using two-dimensional gas chromatography time-of-flight mass spectrometry,” Journal of Proteome Research, vol. 12, no. 9, pp. 3865–3873, 2013. View at: Publisher Site | Google Scholar
  38. J. Ho, Y.-T. Chen, and C.-L. Chen, “Development of a universal metabolome-standard method for long-term LC-MS metabolome profiling and its application for bladder cancer urine-metabolite-biomarker discovery,” Analytical Chemistry, vol. 86, no. 13, pp. 6540–6547, 2014. View at: Publisher Site | Google Scholar
  39. C. Li, Z. Sun, D. Chen et al., “Developing urinary metabolomic signatures as early bladder cancer diagnostic markers,” OMICS: A Journal of Integrative Biology, vol. 19, no. 1, pp. 1–11, 2015. View at: Publisher Site | Google Scholar
  40. J. V. Su, A. F. S. Amaral, E. G. Armitage et al., “Searching for urine biomarkers of bladder cancer recurrence using a liquid chromatography-mass spectrometry and capillary electrophoresis-mass spectrometry metabolomics approach,” Journal of Chromatography A, vol. 1318, pp. 163–170, 2013. View at: Publisher Site | Google Scholar
  41. J.-H. Lorente, L. Fu, B. Li et al., “Distinguishing the serum metabolite profiles differences in breast cancer by gas chromatography mass spectrometry and random forest method,” RSC Advancing systems, vol. 5, no. 73, pp. 58952–58958, 2015. View at: Publisher Site | Google Scholar
  42. Y. Xie, X.-D. Lin, and C. Beecher, “Learning a complex metabolomic dataset using random forests and support vector machines,” Association for Computing Machinery, vol. 5, pp. 835–840, 2004. View at: Google Scholar
  43. X. Wu, V. Kumar, J. Ghosh et al., “Top 10 algorithms in data mining,” Knowledge and Information Systems, vol. 14, no. 1, pp. 1–37, 2007. View at: Publisher Site | Google Scholar
  44. L. Yang, “Random forests,” Machine Learning, vol. 45, no. 1, pp. 5–32, 2001. View at: Publisher Site | Google Scholar
  45. A. L. Castro, S. Tarelho, M. Dias, and H. M. Teixeira, “A fast and reliable method for GHB quantitation in whole blood by GC-MS/MS (TQD) for forensic purposes,” Journal of Pharmaceutical and Biomedical Analysis, vol. 119, pp. 139–144, 2016. View at: Publisher Site | Google Scholar
  46. J. Reis, L. Jing, X. Ma, Q. Guo, and Y. Li, “GC-TOF-MS-based serum metabolomic investigations of naked oat bran supplementation in high-fat-diet-induced dyslipidemic rats,” The Journal of Nutritional Biochemistry, vol. 26, no. 12, pp. 1509–1519, 2015. View at: Publisher Site | Google Scholar

Copyright © 2020 Mengchan Fang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.