Applying 1H NMR Spectroscopy to Detect Changes in the Urinary Metabolite Levels of Chinese Half-Pipe Snowboarders after Different Exercises

Abstract

Monitoring physical training is important for the health and performance of athletes, and real-time assessment of fatigue is crucial to improve training efficiency. The relationship between key biomarkers and exercise has been reported. The aim of this study was to determine the effects of different levels of training exercises on the urine metabolome. 1H NMR-based metabolomics analysis was performed on urine samples from half-pipe snowboarders, and spectral profiles were subjected to PCA and PLS-DA. Our results show that metabolic profiles varied during different stages of exercises. Lactate, alanine, trimethylamine, malonate, taurine, and glycine levels decreased while TMAO and phenylalanine levels increased in the stage with higher amount and intensity of exercise. Although the amount of exercise was reduced in subsequent stage, no significant variations of metabolic profile were found. Metabolic changes induced by training level were analyzed with related metabolic pathway. Studying metabolome changes can provide a better understanding of the physiology of athletes and could aid in adjusting training.

Fuqiu Wang 1 and Jiao Han 2 and Qing He 3 and Zhufeng Geng 2 and Zhiwei Deng 2 and Decai Qiao 1

Academic Editor:Gowda A. Nagana Gowda

1, College of P.E and Sports, Beijing Normal University, No. 19 Xinjiekouwai Street, Haidian District, Beijing 100875, China
2, Center of Analysis and Test, Beijing Normal University, No. 19 Xinjiekouwai Street, Haidian District, Beijing 100875, China
3, School of Chemical Engineering and Technology, Tianjin University, 92 Weijin Road, Nankai District, Tianjin 300072, China

Received 11 March 2015; Revised 10 May 2015; Accepted 11 May 2015; 25 May 2015

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.

1. Introduction

The aim of physical training is to enhance athletic performance and prevent serious injury. Normally, physical training is often associated with a type of tissue trauma called "adaptive microtrauma" [1]. With adequate recovery, the body adapts to the training. Helping the body adapt to higher-intensity stimuli, therefore, is important for physical training [2, 3]. However, overtraining can lead to overtraining syndrome, which is often accompanied by physiological, biochemical, immunological, and psychological symptoms. Margonis et al. [4] found that 37% of professional athletes in all types of sports have suffered from overtraining syndrome at least once in their careers. Hence, monitoring physical training is essential. Many studies indicate that exercises stimulate the body producing changes in the metabolite levels of athletes, regardless of whether the exercise follow acute or chronic, planned or unplanned protocols [5, 6], even in strength-endurance sports [3]. Blood and urine were usually selected as indicators in monitor procedure. Previous studies utilized invasive methods and focused on single diagnostic metabolites, like lactate, phosphocreatine, or certain hormones [7]. However, variations in some of these indicators were not observed through detection of target biochemicals such as serum cortisol, testosterone, or creatine kinase [3].

Following the development of "omics" techniques, the use of metabonomics has increased rapidly in the past decade. As Nicholson et al. [8] stated, metabonomics is a measurement of the dynamic metabolic response of living systems to stimuli or modification. Metabonomics approaches can detect changes in the levels of many small molecules that arise as a response to cells or an organism to endogenous or exogenous stimuli, such as disease, development, environmental changes, or nutrient availability [9].

Metabonomics technologies have already been applied in the physiological monitoring of human movement [10]. However, studies applying metabonomics in sports are not numerous [7]. The reason is that there are many technical barriers, such as identifying plenty of metabolites and interpreting acquired data through advanced multivariate statistical analysis [9]. Hence, metabolomics is conducted by using sensitive analytical techniques, such as liquid chromatography-mass spectrometry (LC-MS), gas chromatography-mass spectrometry (GC-MS), and nuclear magnetic resonance (NMR) spectroscopy [11]. To date, many studies have analyzed metabolite profiles or conducted metabolite fingerprinting. Recent metabolomics studies of physical training have focused on blood serum or plasma in human and animal models. These studies have studied the effects of strenuous endurance exercises [12], strength-endurance training [3], specific diets [13], or ingestion of different beverages after exercise [14-17]. The effects of acute and chronic exercise and different types of exercise on the human urine metabolome have also been studied [6, 18]. Several studies have used saliva to assess the levels of metabolites in response to exercise and training or to evaluate physical performance [7, 19]. Enea et al. [18] detected changes in the concentration of 11 metabolites after exercise by using ¹ H NMR analysis of urine from female athletes and nonathletes.

Generally, previous studies have shown that exercise with a variety of forms and levels could cause metabolic changes in the body's biological fluids, such as blood, urine, and saliva. The test subjects in majority of studies were under strict requirements (diets, fasting, specific training plans with rest intervals, etc.) to obtain variations of key metabolites, aiming for hypotheses testing regarding specific biochemical pathways [6]. However, under actual circumstances, athletic training follows schedules, which was a combination of exercises in various forms. It is important to control the amount and intensity of physical exercise so as to obtain positive effects without fatigue. Thus, the present study uses traditional monitoring methods and NMR-based metabonomics to determine changes in the levels of small-molecule metabolites in urine during training process.

We describe nontargeted ¹ H NMR-based metabonomics analyses of urine samples obtained from professional snowboarders after weekly physical training. The aim of this work was to explore how physical exercise broadly impacts and changes metabolite profiles. Furthermore, we determined the feasibility of using metabonomics to assess and monitor periodic physical exercise.

2. Experimental Section

2.1. Sample Collection

Twelve professional half-pipe snowboarders between the ages of 21-25 years with an average Body Mass Index (BMI) of [figure omitted; refer to PDF] kg/m² participated in this study. All athletes had a standard diet and none used pharmaceutical products or tobacco during the training period. Physical training consisted of typical waist and abdominal strength and endurance and trampoline exercises, which were divided into three stages (stages I, II, and III; see Figure 1). Before the three-week training, subjects underwent a one-week adaptive training period. One training stage consisted of 5.5 days of training and 1.5 days of rest. All urine and blood samples were collected before breakfast on Monday morning. At the beginning of formal training, subjects completed strength, endurance, and trampoline exercises of a predefined level for one week. The next week, the amount and intensity of exercise were both increased by about 10% (warm up and stretch are not included in intensity calculation). Since there is an individual difference in the athletes, especially the gender difference, the changes of training level were quantitatively estimated according to the heart rate of athletes and the professional coach's personalized training. During the third week, the training intensity is the same as the second week, but the amount of exercise was reduced by about 15%. Nine urine samples were collected for stage I, and ten samples were collected for stage II. Only eight urine samples were collected for stage III. Blood samples were collected for all three stages. Urine samples were stored at -20°C with sodium azide until analysis.

Figure 1: The flowchart for illustration of three-stage/level training exercises.

[figure omitted; refer to PDF]

2.2. Sample Preparation and Clinical Data Collection

Urine samples were thawed and centrifuged at 13,000 ×g for 10 min. Supernatant aliquots of 300 μ L were transferred to standard 5 mm NMR tubes and mixed with 300 μ L phosphate buffer (pH 7.4) and 50 μ L D₂ O containing sodium 3-(trimethylsilyl)-1-propanesulfonate (DSS, 0.1% w/v in D₂ O).

Serum cortisol and testosterone concentrations were measured by automatic biochemical analyzer. CK concentration was measured using diagnostic kits. BUN was measured by the urease UV method. Routine blood tests were used to analyze WBC, LYM, MXD, NEUT, RBC, HGB, HCT, MCV, and PLT.

2.3. ¹ H NMR Spectroscopy of Urine Samples

¹ H NMR acquisition was performed using a Bruker DRX-500 spectrometer equipped with a BBFO probe with [figure omitted; refer to PDF] -axis gradients maintained at 298 K. One-dimensional NMR spectra were recorded using a standard 1D NOESY pulse sequence with water suppression. A total of 256 scans were collected in 32 K data points with a spectral width of 8012.8 Hz. The acquisition time was 2.04 s, and relaxation delay was 2 s. All spectra were processed using MestReNova software (v. 8.1.2; Mestrelab Research SL, Santiago de Compostela, Spain). Spectra were manually referenced to the chemical shift of DSS before manual phase and baseline correction. On the purpose of identification, DQF-COSY and HSQC spectrum of representative sample were acquired (see Figures 1S, 2S, and 3S in Supplementary Material available online at http://dx.doi.org/10.1155/2015/315217).

2.4. Statistical Analysis

The spectral region of δ 0.5-10 ppm of all spectra was segmented into 0.02 ppm portions, excluding the 4.2-6.2 ppm region. Signal intensity was then normalized to the total spectral area and imported into Matlab (R2010a; Mathworks, Inc., Natick, MA, USA).

Multivariate data analysis of NMR data was performed using PLS Toolbox (Eigenvector Research, Inc., Manson, WA, USA). A principal component analysis (PCA) model was applied to determine the distribution and separation trends. To enlarge the difference and determine which variables contribute to separation, a partial least squares projection to latent structures and discriminant analysis (PLS-DA) model were then constructed. Cross-validation was conducted to ascertain the predictive ability of the model using the leave-one-out method. [figure omitted; refer to PDF] represents the variance captured by the model and [figure omitted; refer to PDF] represents the prediction ability of the model. Furthermore, a permutation test was applied to validate whether models were overfit.

To identify important biomarkers, the clinical parameters of three groups were analyzed by Student's [figure omitted; refer to PDF] -test. A nonparametric version of one-way ANOVA test [20] was also performed on NMR data of identified metabolites. The statistical significance level was set at [figure omitted; refer to PDF] .

Hierarchical clustering analysis [21] was used to assemble all variables into a single hierarchical tree and cluster relatively homogeneous elements together. In this study, we generated dendrograms and heat maps in Matlab based on measured data.

3. Results

3.1. Clinical Data

There were no significant differences in the concentrations of testosterone, cortisol, creatine kinase, BUN, WBC, LYM, or NEUT between the first two stages. The concentrations of MXD, RBC, MCV, and PLT between the first two stages were significantly different, as were the HGB values in samples from male individuals. However, most of the clinical parameters, such as MXD and RBC did not significantly differ between the second and third groups. MCV, PLT, and HGB concentrations from male individuals reached levels of statistical significance in the last two groups. Comparing the levels of MCV, PLT, and HGB (M) in these three groups, the HGB (M) level rose in the second stage and decreased in the third stage, whereas MCV and PLT values showed decreased levels in the second stage with higher levels approaching initial values in the third stage (Table 1).

Table 1: Clinical indicators and blood cell levels in different stages.

	Stage I	Stage II	Stage III	[figure omitted; refer to PDF] value	[figure omitted; refer to PDF] value
T (ng/dL)
M	571.4 ± 173.0	557.8 ± 176.5	600.0 ± 148.8	0.68	0.41
F	48.6 ± 13.9	45.3 ± 6.9	44.1 ± 11.1	0.48	0.75
C (μ g/dl)	17.8 ± 2.5	17.01 ± 2.98	17.94 ± 4.5	0.11	0.50
CK (U/L)	620.3 ± 1007.7	163.6 ± 111.0	153.6 ± 131.3	0.14	0.84
BUN (mmol/L)	3.8 ± 0.8	4.0 ± 1.0	3.6 ± 0.7	0.49	0.08
WBC ( [figure omitted; refer to PDF] 1 [figure omitted; refer to PDF] /μ L)	5.3 ± 0.9	5.2 ± 1.0	5.4 ± 1.2	0.71	0.60
LYM ( [figure omitted; refer to PDF] 1 [figure omitted; refer to PDF] /μ L)	2.3 ± 0.5	2.3 ± 0.7	2.2 ± 0.7	0.83	0.73
MXD ( [figure omitted; refer to PDF] 1 [figure omitted; refer to PDF] /μ L)	0.5 ± 0.1	0.3 ± 0.1	0.4 ± 0.2	0.01	0.20
NEUT ( [figure omitted; refer to PDF] 1 [figure omitted; refer to PDF] /μ L)	2.6 ± 0.5	2.7 ± 0.7	2.8 ± 0.8	0.90	0.72
RBC ( [figure omitted; refer to PDF] 1 [figure omitted; refer to PDF] /μ L)	4.7 ± 0.5	4.9 ± 0.7	4.8 ± 0.5	0.02	0.31
HGB (g/L)
M	152.4 ± 6.3	160.8 ± 9.1	149.0 ± 6.6	0.03	0.03
F	129.6 ± 7.7	129.0 ± 8.8	132.6 ± 12.1	0.77	0.13
HCT	0.42 ± 0.04	0.43 ± 0.05	0.43 ± 0.04	0.07	0.57
MCV (μ m3 )	89.3 ± 2.4	88.3 ± 2.6	89.3 ± 2.8	2.6 [figure omitted; refer to PDF] 10-3	5 [figure omitted; refer to PDF] 10-3
PLT ( [figure omitted; refer to PDF] 1 [figure omitted; refer to PDF] /μ L)	226.6 ± 47.5	165.0 ± 34.2	223.8 ± 38.4	1.5 [figure omitted; refer to PDF] 10-4	3.65 [figure omitted; refer to PDF] 10-5

Values are given as the means ± SD.

T: testosterone; C: cortisol; CK: creatine kinase; BUN: blood urea nitrogen; WBC: white blood cells; LYM: lymphocyte count; MXD: mononuclear cell count; NEUT: neutrophil count; RBC: red blood cells; HGB: hemoglobin; HCT: hematocrit; MCV: mean corpuscular volume; PLT: platelet; M: male; F: female.

3.2. ¹ H NMR Spectra and Metabolic Changes

A representative ¹ H NMR spectrum of a urine sample from an athlete is shown in Figure 2. Peaks of major metabolites could be observed in the spectra and were identified by comparing with a database [22] and data reported in the literature [6, 23]. With the urine samples, the PCA model was used to observe the distribution and the separation trends (Figure 3). A clear separation was found in the score plot of the PCA model based on the first two stages. To enlarge the difference and detect important variables, the PLS-DA model was established using one predictive component with [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , as shown in Figure 4(a). The model was validated by a permutation test (Figure 5). As Figure 4 showed, [figure omitted; refer to PDF] value of model was much higher than any other [figure omitted; refer to PDF] from permutated model. This result indicated that the model with satisfied predictability was not overfitted. Since the PCA model of NMR data from stages two and three did not show clear separation (Figure 3), PLS was applied with [figure omitted; refer to PDF] , which indicates an unsatisfied model with low predictive ability. Figure 4(b) displays a loading plot of the PLS-DA model from the first two stages. The loading plot illustrates the trends of major metabolite levels between the two groups. Concentrations of major metabolites, including TMAO, taurine, and hippurate, were higher in urine samples from the relatively higher training intensity and amount group, whereas creatinine, glycine, and an unknown signal at 3.68 ppm, which falls in the spectral region of sugars, polyols, and amino acid CH-groups, were obviously higher in the first stage group. To validate this trend and explore other metabolites that contribute to variations between the two groups, an integral value of identified metabolites was compared using the significance test. Decreased urinary lactate, alanine, trimethylamine, malonate, taurine, and glycine levels, and increased urinary TMAO and phenylalanine levels, were observed in samples from the higher training intensity and amount stage (Table 2). These urinary metabolites with statistically significant levels could be used as potential biomarkers contributing to metabolic changes in athletes in different physical states. In addition to those mentioned above, other metabolites, such as succinate, dimethylamine, creatinine, and an unknown signal in the region of 3.5-3.7 ppm, were close to the margin of statistical significance. The levels of these metabolites were decreased in urine samples from stages one to two. The observed trends from the significance test were in agreement with those from the PLS loading plot.

Table 2: List of peaks intensity with [figure omitted; refer to PDF] from urinary spectra of first two stages.

Number	[figure omitted; refer to PDF] (ppm)	Metabolite	Trend	Significance
1	1.32	Lactate	[arrow down]	[figure omitted; refer to PDF]
2	1.47	Alanine	[arrow down]	[figure omitted; refer to PDF]
3	2.91	Trimethylamine	[arrow down]	[figure omitted; refer to PDF]
4	3.14	Malonate	[arrow down]	[figure omitted; refer to PDF]
5	3.26	TMAO	[arrow up]	[figure omitted; refer to PDF]
6	3.42	Taurine	[arrow down]	[figure omitted; refer to PDF]
7	3.55	Glycine	[arrow down]	[figure omitted; refer to PDF]
8	7.42	Phenylalanine	[arrow up]	[figure omitted; refer to PDF]

The potential biomarkers were labeled with downregulation ([arrow down]) and upregulation ([arrow up]). [figure omitted; refer to PDF] and [figure omitted; refer to PDF] .

Figure 2: Typical 500 MHz 1D-NOESY NMR spectra of the urine sample from an athlete. (1) Lactate (d, 1.32 ppm); (2) alanine (d, 1.47 ppm); (3) acetate (s, 1.91 ppm); (4) succinate (s, 2.39 ppm); (5) citrate (dd, 2.51, 2.54, 2.65, and 2.68 ppm); (6) dimethylamine (s, 2.71 ppm); (7) creatinine (s, 3.03 ppm, s, 4.04 ppm); (8) TMAO (s, 3.26 ppm); (9) taurine (t, 3.26 ppm, t, 3.42 ppm); (10) glycine (s, 3.55 ppm); (11) phenylalanine (m, 7.42 ppm); (12) hippurate (d, 7.82 ppm).

[figure omitted; refer to PDF]

Figure 3: PCA score plot of NMR data from three stages (stage I: red dots; stage II: blue dots; stage III: green dots).

[figure omitted; refer to PDF]

Figure 4: (a) Score plot of PLS model based on NMR data from first two stages (red dots represent samples from stage I; blue dots represent samples from stage II). (b) PLS loading plot. Positive value corresponds to metabolites at higher concentration in the urine sample of stage I and vice versa.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

Figure 5: Permutation test of PLS model based on NMR data of first two stages.

[figure omitted; refer to PDF]

3.3. Hierarchical Cluster Analysis

To determine the relationships among major metabolites, correlation analysis and hierarchical clustering analysis were performed on NMR data and illustrated by a dendrogram and heat map (Figure 6). The correlation structure was conveniently shown through different colors in correlation map. Metabolites that highly correlate should cluster together in the dendrogram. Three clusters were clearly observed: (1) dimethylamine (δ 2.71) and creatinine (δ 3.03); (2) acetate (δ 1.91), alanine (δ 1.47), lactate (δ 1.32), citrate (δ 2.53), succinate (δ 2.39), and glycine (δ 3.55); and (3) taurine (δ 3.42), TMAO (δ 3.26), phenylalanine (δ 7.42), and hippurate (δ 7.84).

Figure 6: Correlation map (a) of the identified metabolites in stage I and stage II. Heat map (b) for identified metabolites in samples from first two stages (group A: first training stage; group B: second training stage). The upper dendrogram showed hierarchical clustering of variances corresponding to column data in the heat map. Rows: samples; columns: chemical shift of significant signal of metabolites.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

4. Discussion

As a routine examination, blood test was accomplished each week through the entire training program. During this period, hormone concentrations in serum, such as testosterone and cortisol, did not show significant variance because of one-day rest before sampling. This result is in accordance with a previous study [21]. As a sensitive indicator of training, the average value of creatine kinase concentration of first-stage samples was [figure omitted; refer to PDF] U/L, which is much higher than the other two sample sets. By inspecting the original data (298, 228, 332, 125, 338, 171, 60, 745, 1011, 83, 3842, and 210 U/L), three values are out of the normal range (<400 U/L). However, the abnormal creatine kinase values returned to normal in the next training stage, which illustrates adaptation to training from the so-called repeated-bout effect [24]. Mononuclear cell counts decreased markedly at first but then increased along with the amount reduction of training, whereas other leukocytes (lymphocytes and neutrophils) did not show significant changes during training. The normal range of hemoglobin values was different between male and female samples, indicating that it depends on gender. The significant variations in the hemoglobin values in male samples suggest a reasonable training program for the athletes.

In this study, the urinary metabolic profiles of professional half-pipe snowboarders were investigated and multivariate analysis was performed on the datasets to detect variations in metabolites that are caused by changes in the amount and intensity of training. After performing high-throughput metabolomics analyses, the dynamic metabolic profiles of urine samples from the athletes were recorded. Variation between the first two training stages was clearly observed, which reflects changes of physical states of the athletes. In contrast, differences in the metabolic profiles between the last two training stages were not significant based on the PLS model.

By analyzing the differential metabolites caused by training level, they were related to each other through metabolic pathways and are depicted in Figure 7. The tricarboxylic acid (TCA) cycle was considered an important pathway related to energy generation in exercise science research [25]. As an intermediate in the TCA cycle, the decreased levels of succinate between stages one and two indicate that the TCA cycle could be affected by long-term training based on different exercise loads.

Figure 7: Summary of the metabolic pathways with indication of the effect of first two weeks of training on urine metabolite levels. Color coding: black, no change; red, upregulation; green, downregulation.

[figure omitted; refer to PDF]

Based on the heat map in Figure 6, lactate (peak at 1.32 ppm) exhibits a positive correlation with alanine (peak at 1.47 ppm). This result is in accordance with expectation because both of them served as substrates in gluconeogenesis [26, 27]. The concentration of creatinine also showed a distinct trend toward significance in the PLS model. It displayed a decreased level after one week training of higher level compared to the initial week. From here we could presume that the storage ability of phosphagen in muscle was improved, and consequently the level of excreted creatinine decreased after the second stage.

Trimethylamine levels in urine samples decreased in the second training stage. In contrast, the level of trimethylamine-N -oxide (TMAO) which is synthesized from endogenous trimethylamine increased. Because TMAO level could be elevated by consuming seafood which was included in the daily diet, the increased level of TMAO in urine samples may be due to increased food intake with high physical consumption during the second training stage. TMAO concentration also displayed remarkable variations in urine samples collected before and after exercise in some studies [6, 17, 18]. Thus, regulation of one-week training could be another reason why concentration level changed.

In addition, hippurate and its precursor phenylalanine [28] displayed a high correlation in hierarchical cluster analysis. As an essential amino acid, phenylalanine plays an important role in the synthesis of neurotransmitters. The increased level indicates metabolic adaption in the athlete as a response to the increased training load.

5. Conclusion

The present study evaluated changes in the urinary metabolic profiles of half-pipe snowboarders as a result of training, using a metabolomics strategy. The PCA score plot showed the tendency to separate athletes with different training amounts and intensities. Furthermore, PLS model was established to explore the differential metabolites caused by internal metabolic changes. These results proved the perturbations of related pathway and demonstrated that organisms reach a relatively stable physical state to adapt to the training load after long-term training. The study of metabolic adaption to different training level could help reveal the influence of long-term exercise on physiological status and provide useful information for formulating training programs.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (31171138) and Beijing Natural Science Foundation (5142012). Firstly, special thanks are due to Yang Li minister of Chinese Ski Association and medical staff Yang Cheng from Shenyang Sport University. Secondly, thanks are due to Yanchun Li from the Center of Science and Research of Beijing Sport University.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] L. L. Smith, "Cytokine hypothesis of overtraining: a physiological adaptation to excessive stress?," Medicine and Science in Sports and Exercise , vol. 32, no. 2, pp. 317-331, 2000.

[2] D. J. Smith, D. Roberts, "Effects of high volume and/or intense exercise on selected blood chemistry parameters," Clinical Biochemistry , vol. 27, no. 6, pp. 435-440, 1994.

[3] B. Yan, A. Jiye, G. Wang, H. Lu, X. Huang, Y. Liu, W. Zha, H. Hao, Y. Zhang, L. Liu, S. Gu, Q. Huang, Y. Zheng, J. Sun, "Metabolomic investigation into variation of endogenous metabolites in professional athletes subject to strength-endurance training," Journal of Applied Physiology , vol. 106, no. 2, pp. 531-538, 2009.

[4] K. Margonis, I. G. Fatouros, A. Z. Jamurtas, "Oxidative stress biomarkers responses to physical overtraining: implications for diagnosis," Free Radical Biology and Medicine , vol. 43, no. 6, pp. 901-910, 2007.

[5] J. K. Nicholson, J. C. Lindon, "Systems biology: metabonomics," Nature , vol. 455, no. 7216, pp. 1054-1056, 2008.

[6] A. Pechlivanis, S. Kostidis, P. Saraslanidis, A. Petridou, G. Tsalis, V. Mougios, H. G. Gika, E. Mikros, G. A. Theodoridis, "H-1 NMR-based metabonomic investigation of the effect of two different exercise sessions on the metabolic fingerprint of human urine," Journal of Proteome Research , vol. 9, no. 12, pp. 6405-6416, 2010.

[7] C. Santone, V. Dinallod, M. Paci, S. D'Ottavio, G. Barbato, S. Bernardini, "Saliva metabolomics by NMR for the evaluation of sport performance," Journal of Pharmaceutical and Biomedical Analysis , vol. 88, pp. 441-446, 2014.

[8] J. K. Nicholson, J. C. Lindon, E. Holmes, "'Metabonomics': understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological NMR spectroscopic data," Xenobiotica , vol. 29, no. 11, pp. 1181-1189, 1999.

[9] A. M. Weljie, F. R. Jirik, "Hypoxia-induced metabolic shifts in cancer cells: moving beyond the Warburg effect," International Journal of Biochemistry and Cell Biology , vol. 43, no. 7, pp. 981-989, 2011.

[10] J. C. Lindon, E. Holmes, M. E. Bollard, E. G. Stanley, J. K. Nicholson, "Metabonomics technologies and their applications in physiological monitoring, drug safety assessment and disease diagnosis," Biomarkers , vol. 9, no. 1, pp. 1-31, 2004.

[11] M.-E. Dumas, E. C. Maibaum, C. Teague, H. Ueshima, B. Zhou, J. C. Lindon, J. K. Nicholson, J. Stamler, P. Elliott, Q. Chan, E. Holmes, "Assessment of analytical reproducibility of 1H NMR spectroscopy based metabonomics for large-scale epidemiological research: the INTERMAP study," Analytical Chemistry , vol. 78, no. 7, pp. 2199-2208, 2006.

[12] E. Pohjanen, E. Thysell, P. Jonsson, C. Eklund, A. Silfver, I.-B. Carlsson, K. Lundgren, T. Moritz, M. B. Svensson, H. Antti, "A multivariate screening strategy for investigating metabolic effects of strenuous physical exercise in human serum," Journal of Proteome Research , vol. 6, no. 6, pp. 2113-2120, 2007.

[13] G. M. Kirwan, V. G. Coffey, J. O. Niere, J. A. Hawley, M. J. Adams, "Spectroscopic correlation analysis of NMR-based metabonomics in exercise science," Analytica Chimica Acta , vol. 652, no. 1-2, pp. 173-179, 2009.

[14] A. B. Hodgson, R. K. Randell, N. Boon, U. Garczarek, D. J. Mela, A. E. Jeukendrup, D. M. Jacobs, "Metabolic response to green tea extract during rest and moderate-intensity exercise," Journal of Nutritional Biochemistry , vol. 24, no. 1, pp. 325-334, 2013.

[15] E. Chorell, T. Moritz, S. Branth, H. Antti, M. B. Svensson, "Predictive metabolomics evaluation of nutrition-modulated metabolic stress responses in human blood serum during the early recovery phase of strenuous physical exercise," Journal of Proteome Research , vol. 8, no. 6, pp. 2966-2977, 2009.

[16] S. J. Bruce, I. Breton, J. Decombaz, C. Boesch, E. Scheurer, I. Montoliu, S. Rezzi, S. Kochhar, P. A. Guy, "A plasma global metabolic profiling approach applied to an exercise study monitoring the effects of glucose, galactose and fructose drinks during post-exercise recovery," Journal of Chromatography B: Analytical Technologies in the Biomedical and Life Sciences , vol. 878, no. 29, pp. 3015-3023, 2010.

[17] A. Miccheli, F. Marini, G. Capuani, A. T. Miccheli, M. Delfini, M. E. Di Cocco, C. Puccetti, M. Paci, M. Rizzo, A. Spataro, "The influence of a sports drink on the postexercise metabolism of elite athletes as investigated by nmr-based metabolomics," Journal of the American College of Nutrition , vol. 28, no. 5, pp. 553-564, 2009.

[18] C. Enea, F. Seguin, J. Petitpas-Mulliez, N. Boildieu, N. Boisseau, N. Delpech, V. Diaz, M. Eugène, B. Dugue, " ¹ H NMR-based metabolomics approach for exploring urinary metabolome modifications after acute and chronic physical exercise," Analytical and Bioanalytical Chemistry , vol. 396, no. 3, pp. 1167-1176, 2010.

[19] E. Papacosta, G. P. Nassis, "Saliva as a tool for monitoring steroid, peptide and immune markers in sport and exercise science," Journal of Science and Medicine in Sport , vol. 14, no. 5, pp. 424-434, 2011.

[20] R. Howard, A. L. Carriquiry, W. D. Beavis, "Parametric and nonparametric statistical methods for genomic selection of traits with additive and epistatic genetic architectures," G3-Genes Genomes Genetics , vol. 4, no. 6, pp. 1027-1046, 2014.

[21] M. B. Eisen, P. T. Spellman, P. O. Brown, D. Botstein, "Cluster analysis and display of genome-wide expression patterns," Proceedings of the National Academy of Sciences of the United States of America , vol. 95, no. 25, pp. 14863-14868, 1998.

[22] D. S. Wishart, D. Tzur, C. Knox, R. Eisner, A. C. Guo, N. Young, D. Cheng, K. Jewell, D. Arndt, S. Sawhney, C. Fung, L. Nikolai, M. Lewis, M.-A. Coutouly, I. Forsythe, P. Tang, S. Shrivastava, K. Jeroncic, P. Stothard, G. Amegbey, D. Block, D. D. Hau, J. Wagner, J. Miniaci, M. Clements, M. Gebremedhin, N. Guo, Y. Zhang, G. E. Duggan, G. D. MacInnis, A. M. Weljie, R. Dowlatabadi, F. Bamforth, D. Clive, R. Greiner, L. Li, T. Marrie, B. D. Sykes, H. J. Vogel, L. Querengesser, "HMDB: the human metabolome database," Nucleic Acids Research , vol. 35, no. 1, pp. D521-D526, 2007.

[23] L. G. Rasmussen, F. Savorani, T. M. Larsen, L. O. Dragsted, A. Astrup, S. B. Engelsen, "Standardization of factors that influence human urine metabolomics," Metabolomics , vol. 7, no. 1, pp. 71-83, 2011.

[24] M. P. McHugh, "Recent advances in the understanding of the repeated bout effect: the protective effect against muscle damage from a single bout of eccentric exercise," Scandinavian Journal of Medicine & Science in Sports , vol. 13, no. 2, pp. 88-97, 2003.

[25] K. R. Howarth, P. J. LeBlanc, G. J. F. Heigenhauser, M. J. Gibala, "Effect of endurance training on muscle TCA cycle metabolism during exercise in humans," Journal of Applied Physiology , vol. 97, no. 2, pp. 579-584, 2004.

[26] G. Perriello, R. Jorde, N. Nurjhan, M. Stumvoll, G. Dailey, T. Jenssen, D. M. Bier, J. E. Gerich, "Estimation of glucose-alanine-lactate-glutamine cycles in postabsorptive humans: role of skeletal muscle," The American Journal of Physiology-Endocrinology and Metabolism , vol. 269, no. 3, pp. E443-E450, 1995.

[27] F.-C. Tang, "Influence of branched-chain amino acid supplementation on urinary protein metabolite concentrations after swimming," Journal of the American College of Nutrition , vol. 25, no. 3, pp. 188-194, 2006.

[28] H. J. Lees, J. R. Swann, I. D. Wilson, J. K. Nicholson, E. Holmes, "Hippurate: the natural history of a mammalian-microbial cometabolite," Journal of Proteome Research , vol. 12, no. 4, pp. 1527-1546, 2013.