Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Gender-Specific Metabolomic Profiling of Obesity in Leptin-Deficient ob/ob Mice by 1H NMR Spectroscopy

  • Eun-Young Won ,

    Contributed equally to this work with: Eun-Young Won, Mi-Kyung Yoon

    Affiliations Medical Proteomics Research Center, KRIBB, Daejeon, Republic of Korea, Department of Biological Sciences, Korea Advanced Institute of Science and Technology, Daejeon, Republic of Korea

  • Mi-Kyung Yoon ,

    Contributed equally to this work with: Eun-Young Won, Mi-Kyung Yoon

    Affiliation Medical Proteomics Research Center, KRIBB, Daejeon, Republic of Korea

  • Sang-Woo Kim,

    Affiliation Laboratory Animal Center, KRIBB, Daejeon, Republic of Korea

  • Youngae Jung,

    Affiliation Integrated Metabolomics Research Group, Seoul Center, Korea Basic Science Institute, Seoul, Republic of Korea

  • Hyun-Whee Bae,

    Affiliations Integrated Metabolomics Research Group, Seoul Center, Korea Basic Science Institute, Seoul, Republic of Korea, Graduate School of Analytical Science and Technology, Chungnam National University, Daejeon, Republic of Korea

  • Daeyoup Lee,

    Affiliation Department of Biological Sciences, Korea Advanced Institute of Science and Technology, Daejeon, Republic of Korea

  • Sung Goo Park,

    Affiliation Medical Proteomics Research Center, KRIBB, Daejeon, Republic of Korea

  • Chul-Ho Lee ,

    chullee@kribb.re.kr (CHL); gshwang@kbsi.re.kr (GSH); swchi@kribb.re.kr (SWC)

    Affiliation Laboratory Animal Center, KRIBB, Daejeon, Republic of Korea

  • Geum-Sook Hwang ,

    chullee@kribb.re.kr (CHL); gshwang@kbsi.re.kr (GSH); swchi@kribb.re.kr (SWC)

    Affiliations Integrated Metabolomics Research Group, Seoul Center, Korea Basic Science Institute, Seoul, Republic of Korea, Graduate School of Analytical Science and Technology, Chungnam National University, Daejeon, Republic of Korea

  • Seung-Wook Chi

    chullee@kribb.re.kr (CHL); gshwang@kbsi.re.kr (GSH); swchi@kribb.re.kr (SWC)

    Affiliation Medical Proteomics Research Center, KRIBB, Daejeon, Republic of Korea

Abstract

Despite the numerous metabolic studies on obesity, gender bias in obesity has rarely been investigated. Here, we report the metabolomic analysis of obesity by using leptin-deficient ob/ob mice based on the gender. Metabolomic analyses of urine and serum from ob/ob mice compared with those from C57BL/6J lean mice, based on the 1H NMR spectroscopy in combination with multivariate statistical analysis, revealed clear metabolic differences between obese and lean mice. We also identified 48 urine and 22 serum metabolites that were statistically significantly altered in obese mice compared to lean controls. These metabolites are involved in amino acid metabolism (leucine, alanine, ariginine, lysine, and methionine), tricarbocylic acid cycle and glucose metabolism (pyruvate, citrate, glycolate, acetoacetate, and acetone), lipid metabolism (cholesterol and carnitine), creatine metabolism (creatine and creatinine), and gut-microbiome-derived metabolism (choline, TMAO, hippurate, p-cresol, isobutyrate, 2-hydroxyisobutyrate, methylamine, and trigonelline). Notably, our metabolomic studies showed distinct gender variations. The obese male mice metabolism was specifically associated with insulin signaling, whereas the obese female mice metabolism was associated with lipid metabolism. Taken together, our study identifies the biomarker signature for obesity in ob/ob mice and provides biochemical insights into the metabolic alteration in obesity based on gender.

Introduction

Once associated only with more developed countries, obesity is now prevalent in less developed countries, reaching global epidemic proportions across all genders and age groups [1]. Obesity results from excessive caloric intake compared with caloric consumption. The risks of obesity are associated with severe diseases such as noninsulin-dependent diabetes mellitus (NIDDM), cardiovascular diseases, hypertension, hypercholesterolemia, arthritis, asthma, and cancers [2]. Despite extensive studies, the physiological and biochemical basis of obesity remains largely unknown, as it is caused by multiple factors, such as genetic predisposition, environment, and lifestyle as well as complex interactions between these factors.

Among the known genetic factors, leptin, a key hormone regulating a balance between food intake and energy expenditure, first identified in 1994 [3], was highlighted as an important genetic component of human obesity [4], [5] Mutations in leptin or leptin receptor genes have been reported to result in human obesity [6][9]. In particular, leptin-deficiency was shown to be associated with severe early-onset obesity in a 9-year-old girl and her cousin [6] and leptin treatment reduced their body weights to the normal range [10]. The C57BL/6J-Lepob/Lepob (hereafter referred to as ob/ob) mouse contains a spontaneous mutation in the ob gene encoding leptin, resulting in profound obesity, NIDDM, and phenotypes that resemble human obesity [3].

Metabolomics, analogous to proteomics and genomics, is now a widely utilized bio-analytical methodology in systems biology which employs nuclear magnetic resonance (NMR) spectroscopy or mass spectrometry to define small molecules present in biological samples in response to genetic modifications or physiological stimuli [11]. The NMR-based analysis platform, which contains a wealth of quantitative chemical information, in combination with multivariate statistical analysis, has been successfully utilized to yield information pertaining to the qualitative and quantitative alterations of small molecules following various stimuli across data sets. To date, metabolomics has been utilized for biomarker discovery [12] and understanding disease processes [13] and gut microbial-host metabolic interactions [14]. Since obesity is a systemic disorder involving metabolic changes, metabolic profilings of obesity have been performed by using various animal models, including db/db mice [12], Zucker rats [15][17], growth hormone mutant mice [18], and high fat diet (HFD)-induced mice [19][21]. Despite its importance in obese animal model, however, metabolomics analysis has rarely been performed on ob/ob mice.

Here, we report the metabolomic profiling of obesity using the ob/ob mice. Our 1H-NMR-based metabolomics analyses, in combination with multivariate statistical analysis, demonstrate not only significantly altered metabolism in obesity but also gender-specific metabolic changes in obesity.

Materials and Methods

Ethics Statement

All experimental procedures were approved by Institutional Animal Care and Use Committee of Korea Research Institute of Bioscience and Biotechnology, and performed in accordance with the NIH Guide for the Care and Use of Laboratory Animals (NIH Publication 1985, revised 1996).

Animals

C57BL/6J (B6) and ob/ob mice were purchased from the Jackson Laboratory (Bar Harbor, ME) and housed in a humidity (60%)-, temperature (25±2°C)-, and light (12 h light/dark cycle)-controlled environment. The mice were fed commercial standard mouse chow diet (Teklad rodent diet 2018S, Harlan, Madison, WI) and water ad libitum. At 10 weeks of age, we collected 24 h-urine from the mice housed in metabolic cages (Tecniplast, Buguggiate, Italy) and serum from the same mice following overnight fasting. The samples were stored at −80°C until NMR analysis.

Sample Preparation

After urine and serum samples were thawed at room temperature, they were centrifuged at 13,000 rpm for 10 min to remove the pellet. The 400 µL of resulting supernatant from urine sample was mixed with 0.2 M phosphate-buffered saline (PBS) (pH 7.0), 0.018% NaN3, 5 mM 4,4-dimethyl-4-silapentane-1-sulfonic acid (DSS), and 10% D2O to the final volume of 600 µL. For serum sample, 200 µL of resulting supernatant was mixed with 0.9% NaCl and 10% D2O to the final volume of 600 µL. The final samples were placed in a 5-mm NMR tube (Wilmad Lab Glass, Buena, NJ).

1H NMR Spectroscopy

NMR data were acquired on a 600 MHz Agilent NMR spectrometer (Agilent, Santa Clara, California, CA) equipped with a triple-resonance 5-mm HCN salt-tolerant cryoprobe at Korea Basic Science Institute. The NOESYPRESAT with the presaturation of water or water suppressed Carr–Purcell–Meiboom–Gill (CPMG) spin-echo pulse sequence was used to obtain 1H NMR spectra for urine or serum samples, respectively. The NOESYPRESAT spectra were collected with 64 transients into 67,568 data points using a spectral width of 10,000 Hz, relaxation delay of 2 s, and a mixing time of 100 ms. The CPMG spectra were collected with 128 transients into 67,568 data points using a spectral width of 10,000 Hz, and relaxation delay of 2 s.

NMR Data Processing

NMR spectra were phased and baseline-corrected using TOPSPIN software (version 3.0, Bruker Biospin, Rheinstetten, Germany). All acquired NMR spectra were referenced to DSS at 0 ppm or lactate at 1.32 ppm for urine or serum sample, respectively. Spectra were then divided into 0.005 ppm bins excluding the spectral regions 4.70–5.05 ppm and 0.00–0.63 ppm, corresponding to water and DSS (for urine sample), respectively. After normalization of the spectra to the total spectral area to account for differences within the samples, the binning data were converted to the ASCII format, which were further imported into MATLAB (R2008a; MathWorks, Inc., Natick, MA). All spectra were aligned using the correlation optimized warping (COW) method [22] to reduce chemical shift variability across the spectra.

Multivariate Statistical Analysis

The resulting binning data were Pareto-scaled and the targeted profiling data (as shown below) were unit variance (UV)-scaled, which were used for multivariate statistical analysis by using SIMCA-P (version 12.0, Umetrics, Umeå, Sweden). Principle components analysis (PCA), a classical unsupervised multivariate pattern recognition method, was employed to examine the intrinsic variation within a group and to assess the clustering behavior between groups. An orthogonal partial least-squares discriminant analysis (OPLS-DA), a supervised pattern recognition method, was further performed to maximize the variation between groups and to determine the variables that contribute to this variation. The quality of models was validated by determining R2 (goodness of fit parameter) and Q2 (goodness of prediction parameter) values [23].

Targeted Metabolite Profiling

For the identification and quantification of metabolites, processed NMR spectra were imported into Chenomx software, a module of NMR Suite (version 7.1, Chenomx, Edmonton, Canada). The assignment of ambiguous peaks due to peaks overlap was confirmed by spiking standard compounds. The 600 MHz library from Chenomx was used for the identification of individual compounds. The internal reference standard DSS of a known concentration was used to determine their concentrations for urine samples. For serum samples, the intensities of the identified metabolites were used to determine their relative concentrations. The resulting concentration data sets of each metabolite were further used for multivariate statistical analyses including PCA and OPLS-DA, in addition to the binning input data, as described above.

Statistical Analysis of Identified Metabolites

For the detection of statistical differences of the metabolite levels between lean and obese groups, nonparametric Mann-Whitney t-tests were performed using PRISM (version 5.0, GraphPad Software, Inc., La Jolla, CA) with 95% probability level (p<0.05).

Pathway Analysis

For the metabolites that significantly changed in obesity with p<0.05, pathway analysis and metabolite set enrichment analysis (MESA) were performed using MetaboAnalyst, a web-based metabolomics data analysis software [24]. The concentrations of statistically significant metabolites that have p<0.05 in t-tests and VIP score >1 from OPLS-DA were used for quantitative pathway enrichment analysis. The “metabolic pathway-associated metabolite sets” was selected as metabolite set library. All input data were autoscaled.

Results

B6 and ob/ob Mice

The leptin-deficient ob/ob mouse was used as a genetically induced obesity mouse model and the B6 mouse, a commonly studied experimental model of metabolic diseases, was employed as a control. As expected, the ob/ob mice weighed two and three times as much as the B6 mice for male and female, respectively. The sample number and average body weight of male and female B6 and ob/ob mice are listed in Table 1.

thumbnail
Table 1. The number and body weight of the mice used in this study.

https://doi.org/10.1371/journal.pone.0075998.t001

1H NMR Spectra of Urine and Serum Samples

The representative 1H NMR spectra of urine and serum samples were overlaid for obese and lean mice (Figure S1). The clear spectral differences between obese and lean samples were mostly observed at the up-field region between 0 and 4 ppm (Figure 1), which contains information about the aliphatic groups of metabolites. Through the visual inspection of the NMR spectra, more intense peaks were observed for obese mice samples than lean ones, which was even more obvious for female than male and for urine than serum. This suggests potential metabolic changes in obese groups, when compared to lean controls, and gender-specific variations in obese groups.

thumbnail
Figure 1. The overlaid representative 1H NMR spectra of urine (A) and serum (B) samples obtained from lean and obese mice.

Black and red spectra are overlaid and shown for B6 mice and ob/ob mice, respectively. Samples from male and female mice are shown on the left and right, respectively. Each peak is labeled with the number that corresponds to the following metabolites. For urine samples (A), 1: pantothenate 2∶2-hydroxyisovalerate 3∶2-hydroxybutyrate 4∶2-hydroxyisobutyrate, lactate 5: alanine 6: acetate 7: N-acetyl aspartate 8: acetone 9: acetoacetate 10: succinate 11∶2-oxoisocaproate 12∶2-oxoglutarate 13: trimethylamine 14: methyl malonate 15: TMAO 16: taurine 17: creatinine 18: glycerol 19: guanidoacetate 20: pantothenate 21∶2-hydroxyvalerate 22: valine 23: isobutyrate 24: threonine 25: leucine 26: glutaric monomethyl ester 27: suberate 28: N-acetyal glutamate 29: methionine 30: pyruvate 31: methylamine 32: citrate 33∶5-aminolevulinate 34: dimethylamine 35: cis-aconitate 36: creatine 37: indole-3-acetate 38: choline 39: isopropanol 40: trigonelline. For serum samples (B), 1: VLDL/LDL cholesterol 2: isoleucine/leucine 3: valine 4: lactate 5: citrate 6: alanine 7: acetate 8: lysine 9: ornithine 10: TMAO 11: acetoacetate 12: glucose 13: ethylene glycol 14: alanine 15: glycolate 16: creatine 17: acetone 18: pyruvate 19: carnithine 20: methionine 21: phenylalanine 22: choline 23∶4-hydroxyphenylacetate 24: arginine 25: glycine 26: creatine 27: serine.

https://doi.org/10.1371/journal.pone.0075998.g001

Next, we identified and quantified 71 urine and 29 serum metabolites using the 600 MHz library within the Chenomx software. For identification of unambiguous peaks, spiking experiments were performed using several standard compounds. The concentrations of the identified metabolites were determined with respect to the internal standard DSS of a known concentration or lactate for urine and serum sample, respectively.

Multivariate Statistical Analysis

In order to investigate whether there are statistically significant metabolic differences between B6 and ob/ob mice, we carried out multivariate statistical analysis using concentrations of the identified metabolites derived from targeted profiling. The PCA, an unsupervised multivariate pattern recognition method, does not show clear discrimination between B6 and ob/ob mice for urine samples, while it does for serum samples (Figure S2, left panels). Based on the NMR spectra showing the gender differences between B6 and ob/ob mice particularly for urine samples (Figure 1), we hypothesized that the separation between B6 and ob/ob mice was hindered by gender variation. As expected, a gender-specific PCA showed a better separation between B6 and ob/ob mice, as indicated by increased R2 (goodness of fit parameter) and Q2 (goodness of prediction parameter) values [23] (Figure S2, middle and right panels). Such effect was more obvious for the urine than serum samples. For urine samples, there is a separation between obese and lean subjects only when analyzed based on gender, whereas for serum samples, the separation that exists for all genders becomes better when analyzed based on gender. Although similar results were obtained when spectral binning data were used as an input (data not shown), artifacts can arise from the integration of input bins. For this reason, we used metabolic concentration data for further analysis.

To maximize the differences between the groups and determine the variables that contribute to discrimination, the OPLS-DA, a supervised pattern recognition method, was further employed for the same set of metabolic concentration data (Figure 2). The OPLS-DA model was validated using R2 and Q2 values [23], as shown in Figure 2. Although the OPLS-DA showed distinct separation between lean and obese mice regardless of gender (Figure 2, left panels), a gender-specific analysis provided an even better clustering between groups (Figure 2, middle and right panels). This suggests that the metabolic changes occur in obese mice compared to lean controls and such metabolic changes are dependent on gender.

thumbnail
Figure 2. Orthogonal partial least-squares discriminant analysis (OPLS-DA) score plots obtained from the 1H NMR spectra of urine (A) and serum (B) samples from lean and obese mice.

Data for all (left panels), male (middle panels), and female (right panels) are shown. Closed and open circles represent obese and lean mice, respectively.

https://doi.org/10.1371/journal.pone.0075998.g002

Furthermore, to confirm the gender-specific metabolic changes in obesity, an OPLS-DA was also performed based on gender (Figure 3). Indeed, OPLS-DA of urine and serum samples from obese groups showed clear differences between males and females (Figure 3, left panels). Again, such differences were more obvious for urine than serum sample. However, gender variation was also observed for urine and serum samples from lean controls (Figure 3, right panels), indicating intrinsic metabolic variations between male and female. Thus, the gender variation observed in obesity is likely due to the intrinsic metabolic differences between males and females. Taken together, our multivariate statistical analyses indicate statistically significant metabolic differences between obese and lean mice, which can be further differentiated by gender.

thumbnail
Figure 3. OPLS-DA score plots derived from the 1H NMR spectra of urine (A) and serum (B) samples from lean and obese mice based on gender.

Data for obese (left panels) and lean (right panels) are demonstrated. Closed and open circles represent male and female mice, respectively.

https://doi.org/10.1371/journal.pone.0075998.g003

Identification of Biomarker for Obesity

We generated variable importance in project (VIP) plots from the OPLS-DA with a threshold of 1.0 (Figure S3) to identify the metabolites that significantly contribute to the clustering between obese and lean groups. The metabolites with VIP score >1, considered to be strong contributors, include 2-oxoisocaproate, creatine, leucine, 2-hydroxybutyrate, methionine, 2-hydroxyvalerate, suberate, alanine, creatinine, methylamine, trigonelline, choline, glutaric acid monomethyl ester (MME), hippurate, p-cresol, valine, N-acetylglutamate, guanidoacetate, phenylalanine, isopropanol, tyrosine, 5-aminolevulinate, allantoin, 3,5-dibromotyrosine, threonine, isobutyrate, 2-hydroxyisovalerate, and nicotinamide N-oxide for urine sample and arginine, creatine, lysine, very low-density lipoprotein (VLDL)/low-density lipoprotein (LDL) cholesterol, ornithine, glycolate, alanine, carnitine, trimethylamine N-oxide (TMAO), and pyruvate for serum sample (listed from high to low VIP scores).

For the concentrations of the identified metabolites, Mann-Whitney t-test was further employed with 95% probability level (p<0.05) to select the metabolites that were significantly altered in ob/ob mice, compared to lean mice. Using this t-test together with VIP score from OPLS-DA, we were able to select 28 out of 71 urine metabolites and 10 out of 29 serum metabolites for all obese mice including males and females, 37 and 22 urine metabolites from males and females, respectively, and 11 and 15 serum metabolites from male and female mice, respectively, all with 95% probability level (p<0.05) and VIP score >1. They are categorized and listed in Tables 2 and 3. Interestingly, most of the significantly altered metabolites are down-regulated in ob/ob mice, compared to B6 mice (Tables 2 and 3). This is probably because ob/ob mice are hypometabolic due to reduced circulating levels of thyrotropin and thyroid hormones. [25] The up-regulated metabolites include acetoacetate, acetone, citrate, fumarate, 2-oxoglutarate, succinate, trimethylamine (TMA), and 3-hydroxybutyrate for urine sample only from obese females and acetoacetate, acetone, succinate, carnitine, VLDL/LDL cholesterol, and TMAO for serum from obese mice, many of which are involved in tricarboxylic acid (TCA) cycle and lipid metabolism (Tables 2 and 3).

thumbnail
Table 2. Summary of significant urinary metabolites differentiating between ob/ob and B6 mice.

https://doi.org/10.1371/journal.pone.0075998.t002

thumbnail
Table 3. Summary of significant serum metabolites differentiating between ob/ob and B6 mice.

https://doi.org/10.1371/journal.pone.0075998.t003

To identify the metabolites that contribute to the clustering between genders, we generated VIP plots from the OPLS-DA with a threshold of 1.0, based on the genders (Figure S4). After removing the metabolites that contribute to the clustering between genders in urine and serum samples from lean mice, the metabolites that have unique association with gender variation in obese mice are identified and listed in Table 4. The VIP score and/or p value of some metabolites were significantly altered in obese mice, while the relatively high levels of other metabolites, based on gender, were reversed in obese male mice compared to lean mice. The latter includes urea, succinate, pyruvate, methyl malonate, 2-hydroxyvalerate, glycerol, nicotinamide N-oxide, and threonine for urine sample and alanine and isoleucine for serum sample and the former includes the remainders. Interestingly, 72% of urine metabolites from obese mice with VIP score >1 have a unique association with gender, i.e. 28% of them already exist in lean mice differentiating the genders, whereas 60% of serum metabolites from obese mice have a unique association with genders. This is consistent with the previous statistical analysis that the gender variations were more distinct for the urine than serum samples (Figures 1, 2, and S2).

thumbnail
Table 4. Summary of significant metabolites differentiating between males and females in ob/ob mice.

https://doi.org/10.1371/journal.pone.0075998.t004

Pathway Analysis

To gain insight into the metabolic mechanism of obesity, metabolic pathways of the significantly altered metabolites were analyzed using the “pathway analysis” module within the MetaboAnalyst software [24]. We identified a total of 27 and 31 distinct metabolic pathways that were significantly altered in the urine (Table S1) and serum (Table S2) samples from obese group, respectively. Next, we performed metabolite set enrichment analysis (MSEA) for the same set of significant metabolic concentration data to test if there arebiologically significant groups of metabolites that are enriched in obesity. Using this method, we identified a number of potentially affected metabolic pathways in urine and serum from obese mice and found that many of them were involved in amino acid metabolism (Figure 4). This was expected because many of the significantly altered metabolites with obesity were amino acids and their degradation products, as listed in Tables 2 and 3. In addition to amino acid metabolism, other affected pathwayswere involved in lipid-related metabolisms, which included betain metabolism, phospholipid biosynthesis in both urine and serum and steroid biosynthesis, steroidogenesis, and sphingolipid metabolism in serum. This is consistent with the fact that obesity is characterized by excessive fat deposits in tissue [26]. The enriched pathways in amino acid metabolism and lipid metabolism are related to biotin metabolism in serum, as biotin is involved in fatty acid synthesis and amino acid catabolism by acting as a prosthetic group for biotin-dependent carboxylases [27]. Urea cycle, another perturbed pathway in both urine and serum, is related to creatine metabolism [28].

thumbnail
Figure 4. Summary of pathway enrichment analysis of urine (A) and serum (B) samples.

Each pathway is ranked from low (red) to high (white) Holm p-value, corresponding to the most and the least significant values, respectively.

https://doi.org/10.1371/journal.pone.0075998.g004

For the analysis of gender-specific metabolic pathway changes, we performed MSEA based on gender (Figure S5). The input data for this analysis were concentrations of 37 and 22 urine metabolites from male and female mice, respectively, and 11 and 15 serum metabolites from male and female mice, which were significantly altered in obesity with VIP score >1 and p<0.05 (Tables 2 and 3). We computed the normalized log values of – (Holm p-value) (PN) to enable direct comparison of the relative impact of individual pathways within the same gender with high PN having the high impact. The comparison of PN values of MSEA for each sample (Figure 5) demonstrated that amino acid metabolisms that showed distinct differentiation between obese and lean mice (Figure 4, Tables 2 and 3), did not significantly differ between male and female mice. This suggests that amino acid metabolites are significant biomarkers of obesity regardless of gender. The metabolic pathways enriched only in male were associated with diabetes-related pathways, including insulin signaling, glycolysis and galactose metabolism. In contrast, female-specific metabolic features included lipid-related metabolisms, such as betaine metabolism, sphingolipid metabolism, and phospholipid biosynthesis. Relatively large PN values were observed for steroid biosynthesis, steroidogenesis, and sphingolipid metabolism in the females.

thumbnail
Figure 5. Summary of the relatively enriched pathways urine (A) and serum (B) samples based on gender.

The log values of – (Holm p-value) are normalized (PN) within the same gender to enable the direct comparison of the relatively enriched metabolic pathways between male and female. The pathways for male and female are shown as black and gray bars, respectively.

https://doi.org/10.1371/journal.pone.0075998.g005

Discussion

Our metabolomic profiling of obesity in ob/ob mice reveals significant gender variations in obesity as well as distinct metabolic differences between obese and lean mice, as monitored by 1H NMR spectroscopy. Despite extensive metabolomics studies on obesity, gender has been rarely considered to be an important factor in analyzing metabolomics data in obesity to date. Here, we report the systemic metabolomic profiling of obesity based on gender. Despite its importance in the obesity animal model, ob/ob mice have rarely been utilized for the metabolomics study of obesity. To our knowledge, only one report using this mouse model has been reported, but was only limited to analysis of fat content of fatty pancreas [29]. Similar to ob/ob mice, leptin receptor-deficient db/db mice have been used in several metabolomic studies as a model of NIDDM. [12], [29], [30] Although the several phenotypes are shared by both mouse models (e.g. obesity and insulin resistance), ob/ob mouse model differs from db/db model such that leptin-deficiency was shown to be specifically associated with severe early-onset obesity in obese humans [6]. The metabolic alterations and gender variations observed in ob/ob mice, compared with B6 mice, are discussed below.

Metabolic Alterations in Obesity

Amino acid metabolism.

The most significant metabolites differentiating between obese and lean mice are related to amino acid metabolism (Figure 4 and Tables 2 and 3). In particular, we observed significantly depleted level of leucine (VIP score of 1.75 for all; p<0.001) in urine samples from obese mice compared with those from lean controls (Table 2). The depletion of leucine is consistent throughout our data. The depletion of 2-oxoisocaproate, a degradation product from leucine by leucine transaminase, and pyruvate and alanine, the metabolites related to leucine catabolism and elevation of acetoacetate, an end product of leucine in obese mice, compared to lean mice (Tables 2 and 3), also indicate the depletion of leucine. Leucine plays key roles in modulation of insulin signaling [31], protein synthesis in skeletal muscle [32] and production of alanine/glutamine in skeletal muscle [33], [34] via mammalian target of rapamycin (mTOR) pathway [35], [36]. The reduced levels of leucine were observed in HFD-induced obese Sprague-Dawley rat [37], HFD-B6 mice [19], [20], and NIDDM individuals [38], which are consistent with our results (Table 2). Importantly, the roles of leucine in the regulation of muscle protein and insulin signaling seem to be consistent features throughout our data. Since leucine stimulates muscle protein synthesis and inhibits protein degradation in skeletal muscle and liver [32], [39], significant depletion of leucine in obesity is related to skeletal muscle atrophy. Indeed, the urinary excretion of creatinine, degraded from creatine in muscle in proportion to muscle mass [40], was significantly reduced in obese mice compared to lean mice (VIP score of 1.59; p<0.001, Table 2), suggesting skeletal muscle atrophy. Consistent with this, obesity has been shown to be associated with skeletal muscle atrophy [41]. Regarding the role of leucine in insulin signaling, down-regulated glycolysis as indicated by depletion of pyruvate and lactate, a feature of insulin resistance, was observed in obese mice (Table 3). In addition, significantly reduced levels of alanine in both urine and serum from obese mice (Tables 2 and 3) are related to depletion of leucine because leucine is the major nitrogen source for alanine synthesis in muscle. [33], [34] Previously, leucine has been suggested to be used to treat obesity [35]. Even leucine alone has been shown to reduce insulin resistance, improve glucose tolerance, and restore metabolic abnormalities in HFD-induced mice [42], and reduce insulin resistance in db/db mice [43]. Taken together, leucine can be used as a significant biomarker for obesity.

Other significantly altered amino acids in urine and/or serum samples from ob/ob mice include arginine, lysine, methionine, glycine, serine, threonine, and tryptophan (Tables 2 and 3). In particular, the significantly lowered level of arginine in serum sample from obese mice (VIP score of 1.86; p<0.001, Table 3) is in agreement with the anti-obesity effect of arginine [44], [45]. Several lines of evidences have shown that dietary arginine supplementation reduced obesity in genetically obese rats, HFD-induced rats, and obese humans by promoting the oxidation of glucose and fatty acids as well as by decreasing de novo synthesis of glucose and triacylglycerols [45][47]. In addition, the depletion of lysine and methionine suggests increased insulin resistance. Consistent with previous metabolomic profiling of HFD-mice [20], we observed significantly lowered levels of lysine in serum from obese mice compared to lean controls (VIP score of 1.59; p<0.001, Table 3). Lysine is involved in energy metabolism, implying the perturbations in insulin resistance [20]. The concentration of methionine, a key metabolite in lipid homeostasis and insulin sensitivity [19], was significantly reduced in ob/ob mice compared to lean controls (VIP score of 1.65; p<0.001 for urine sample, Table 2), which is consistent with the previous report on HFD-B6 mice [19]. Another altered amino acid glycine was also depleted in obese samples (Table 3), which is consistent with previous reports [16], [21], [48], [49]. The depletion of glycine was likely due to the decreased levels of serine and pyruvate in serum sample from obese mice (Table 3), because glycine can be generated from serine, which is, in turn, derived from pyruvate [40]. In addition, it is related to the lowered level of urea in urine sample from female obese mice (Table 2) because glycine is an important precursor for urea biosynthesis [40].

TCA cycle and glucose metabolism.

The TCA cycle is a series of biochemical reactions that are utilized by all aerobic organisms to produce energy. Acetate, in the form of acetyl coenzyme A (acetyl-CoA), derived from pyruvate via glycolysis, enters the TCA cycle to produce citrate and other metabolites under aerobic conditions. Our metabolomics analysis of urine demonstrates elevated levels of several TCA cycle metabolites including acetoacetate, acetone, citrate, fumarate, 2-oxoglutarate, and succinate particularly for obese females (Tables 2), suggesting the up-regulation of TCA cycle in obese female mice. In particular, the level of citrate was elevated in the urine from the obese female mice, whereas depleted in the serum from the obese female mice (Tables 2 and 3). The citrate is endogenously generated from fatty acid and glucose metabolism and regulated by insulin and glucose levels [50]. The increased levels of citrate were observed in alloxan diabetic rats [51] and HFD-induced obese animals [20], [21] due to hyperglycaemia and insulin resistance [20], while the lowered levels were seen in other obese animals [17], [18], [21]. Despite this discrepancy regarding the levels of citrate, alterations of citrate levels suggest the disturbances in insulin or glucose levels in obesity. In contrast to the TCA cycle, depleted levels of pyruvate and acetate (Tables 2 and 3), corresponding to the entry points into the TCA cycle, were observed, suggesting down-regulated glycolysis in obesity. In addition, the level of lactate, a product of pyruvate under anaerobic condition, was lowered in serum from obese mice (Table 2), again suggesting down-regulated glycolysis under anaerobic condition.

Lipid metabolism.

Fat deposits in tissue are a significant feature of obesity and are often associated with elevated levels of free fatty acids or cholesterol in the blood [26]. We observed significant elevation of VLDL/LDL cholesterol in serum from obese mice (VIP score of 1.50; p<0.001, Table 3), which is consistent with the previous studies [15], [21], [26]. Another key metabolite identified in our study is carnitine. Carnitine plays key roles in fatty acid metabolism by participating in the transport of fatty acids from the cytosol into the mitochondria to produce energy via β-oxidation [52], [53]. Due to the direct association of carnitine levels with β-oxidation, depleted levels of carnitine was previously observed in liver tissue of HFD-induced mice [54] and blood of obese humans [55]. However, elevated carnitine levels were also observed in the serum of HFD-induced mice [26], [54]. In this study, we observed significantly high levels of carnitine in both obese male and female serum, compared to lean controls. Taken together, whether increased or decreased, the significant changes in carnitine level indicate the disturbance in energy metabolism with obesity.

Creatine metabolism.

Creatine and creatine phosphate serve as a rapid source of high-energy phosphate by the reaction of creatine phosphate with ADP to produce ATP. They are broken down into creatinine in muscle [56]. Since creatinine is usually produced at levels proportional to muscle mass [40] and obesity is associated with skeletal muscle atrophy [41], depletion of creatine, creatine phosphate, and creatinine is expected. Indeed, our metabolomic analysis showed a significant depletion of urinary excretion of creatine, creatine phosphate, and creatinine and serum concentration of creatine in obese mice compared to lean controls (Tables 2 and 3). In addition, depletion of creatine is also associated with the reduced levels of arginine and glycine in obese serum (Table 3), because creatine is biosynthesized from these amino acids [56]. Major metabolic pathway networks, including amino acid metabolism, TCA cycle, glucose metabolism, and creatine metabolism, are summarized and illustrated in Figure 6.

thumbnail
Figure 6. Major metabolic pathways for the metabolites that are significantly altered in urine and serum from ob/ob mice compared to lean controls, observed by 1H NMR spectroscopy.

The pathways were analyzed on the basis of the KEGG pathway database (http://www.genome.jp/kegg/pathway.html). The levels of the identified metabolites that were elevated and depleted in obese mice, compared to lean mice, are color coded as red and blue, respectively, and the metabolites in gray were not identified in our study. Direction of change is shown for all obese mice including males and females. For the direction of change based on gender, see Tables 2 and 3. The metabolites with asterisk are shown several times in this diagram to avoid the complexity of the pathways. The metabolites with sharp indicate the significant ones in male or female group only. Solid and dashed arrows indicate the single and multiple steps involved, respectively, between two metabolites.

https://doi.org/10.1371/journal.pone.0075998.g006

Gut microbiome-derived metabolism.

Gut microbiome has been reported to be able to influence susceptibility to obesity by altering the efficiency of energy harvest from the diet [57]. Although the two dominant bacterial divisions, the Bacteroidetes and the Firmicutes, are found in a majority of human and mouse gut, their relative abundance is altered by obesity [58]. In obese mice and humans, the number of Firmicutes increased, whereas the number of Bacteroidetes decreased compared to lean individuals. Thus, the identified metabolites and their related pathways, originating from gut microflora, are often altered in obese human and mice. In particular, the pathways from choline to TMA and TMAO are indicative of the gut microbiota status [40], [59]. Our metabolomics analysis showed that the levels of choline and methylamine were significantly decreased in obese urine and serum, whereas the level of TMAO was significantly increased in obese serum, compared to lean controls (Tables 2 and 3). Consistent with our results, methylamine levels were reported to be attenuated in obese phenotype of db/db mice [12]. The elevation of serum TMAO was also observed in GHR mutant obese mice [18], whereas decreased TMAO was reported in serum and urine from HFD-induced obese mice or Zucker rats [16], [19], [20]. An inverse relationship of microbiome-originated metabolites might be due to variations in the mouse strain utilized [60]. Recently, Koeth et al. reported that TMAO can be produced from dietary carnitine by intestinal microbiota, resulting in the acceleration of atherosclerosis in mice [61]. Consistent with this finding, our analysis shows the elevation of carnitine and TMAO in serum from obese mice compared to lean controls (Table 3), suggesting the gut-microbiota-mediated metabolism of TMAO.

Another microbiota-derived metabolite, hippurate, is generated from the conjugation of glycine and benzoic acid that is produced from aromatic amino acids by gut microbes and excreted into the urine [62]. The significance of the urinary excretion of hippurate has been considered for obesity and diabetes [63]. The urinary level of hippurate in our study was significantly decreased, compared to lean controls (p<0.005; VIP score of 1.38, Table 2), which is consistent with previous studies using other obese animal models [16], [20], [64], [65] as well as obese humans [66]. A significantly decreased urinary excretion of p-cresol, generated from tyrosine by gut microflora [67], was observed in the obese mice (p<0.005; VIP score of 1.37, Table 2). In addition, depletion of gut microflora-originated short-chain fatty acid isobutyrate in obese urine (Table 2) suggests decreased bacterial fermentation of leucine and valine, resulting in the alteration of host energy metabolism [68]. Urinary excretion of 2-hydroxyisobutyrate, a compound derived from the microbial degradation of dietary protein, has been shown to be associated with some members of the gut microflora, for example, F. prausnitzii [69]. The urinary level of trigonelline, produced via the conversion of S-adenosylmethionine (SAM) to S-adenosylhomocysteine (SAH) by gut microflora [70], was also significantly decreased in obese mice, when compared to lean controls (p<0.005; VIP score of 1.49, Table 2). This suggested that SAM was consumed and thus depleted in the trans-sulfuration pathway to produce glutathione that was depleted in obesity [66], [70].

Gender Variations in Obesity

Our metabolomics analysis showed significant gender variations in obesity. Although our OPLS-DA analyses showed distinct separation between lean and obese subjects, such discrimination was even more clear when analyzed based on gender (Figure 2). In addition, this gender variation was more obvious for urine than serum samples and for female than male, as observed in the NMR spectra (Figure 1) and OPLS-DA analyses based on obesity (Figure 2) and gender (Figure 3).

Pathway enrichment analysis based on gender (Figure 5) showed that female-specific pathways are mainly related to lipid metabolism. One of the metabolites that significantly contribute to the metabolic variations in female obese group is taurine. Taurine plays key roles in conjugation of cholesterol and bile acids, antioxidation, osmoregulation, and calcium signaling [71]. Previous reports showed that the levels of taurine were decreased [16], [19] and that taurine supplementation decreased serum lipids and body weight in obese human [72]. In our analysis, taurine was significantly depleted in urine from obese females, whereas no significant alteration was observed in obese males (Table 2). Female-specific lipid metabolism may arise from the female sex hormone, estrogen. Previous studies reported the strong positive correlation of estrogen with lipid synthesis and the excretion of citrate [73], [74] because estrogen increases the availability of fatty acids [75] and urinary citrate excretion [74]. In addition to lipid-related metabolism, the urinary excretion of citrate was significantly increased in obese female urine with high VIP score (1.59), whereas only marginal increase was observed in obese male urine with low VIP score (0.22) (Table 2). Other female-specific metabolites are urea and glycine. The urinary level of urea was significantly lowered in obese females (p<0.05; VIP score of 1.31, Table 2), whereas only marginally lowered in obese males (p>0.05; VIP score of 0.35, Table 2). In addition, urea had the highest VIP score (1.80) from the OPLS-DA based on gender (Table 4). The depletion of urea in obese females is related to significant decrease in the level of glycine in serum from obese females (Table 3) because glycine is an important precursor for urea biosynthesis [40]. Contrary to the level of glycine in obese males, it was significantly decreased in obese females (p<0.01; VIP score of 1.24, Table 3).

In contrast to females, males are specifically associated with diabetes-related pathways, including insulin signaling, glycolysis and galactose metabolism (Figure 5). Some of the metabolites that significantly contributed to the metabolic variations in male obese group include pyruvate in urine (Table 2) and glucose in serum (Table 3) from obese male mice. These metabolites were significantly decreased compared to obese females (Tables 2 and 3) and contribute to the separation between males and females in obese mice (Table 4). The depletion of urinary pyruvate in obese males suggests down-regulation of glycolysis and thus an insulin resistance. Although some previous studies indicated that there were no gender difference in diabetes and obesity, other reports suggested that the male sex hormone was a risk factor in the development of diabetes due to the changes in the hormone levels in diabetes [76][79]. In particular, significantly low levels of testosterone have been observed in obese humans [80] and diabetic men [76], [78]. Furthermore, hypoandrogenism has been proposed as an early biomarker for the disturbances in insulin and glucose metabolism [77]. Thus, the alterations in insulin signaling unique to obese males in our study may arise from such depletion of testosterone.

These gender-specific metabolic differences are thought to originate from intrinsic metabolic differences between males and females because there are clear metabolic differences between genders even in lean controls (Figure 3). The intrinsic gender-specific metabolic differences have been also reported in normal mice [81] and humans [82], suggesting the fundamental metabolic differences between genders. Thus, it is likely that the fundamental metabolic differences between genders, regardless of metabolic abnormality, led to the gender-specific metabolic variation in obesity. Despite the significant effect of gender on the metabolism of obesity, only one report has been previously published to our knowledge, which showed the gender variation in obese Zucker rats [64]. Based on our results, we suggest that the gender be considered in metabolomics profiling of obesity.

The risks and damages caused by obesity is not only limited to individual health problems; it also leads to social economic problems by increasing health care cost [83]. Given that obesity is a multifactorial problem, early diagnosis and treatment are critical to improve personal health and prevent the resulting economic burden. In this study, our 1H-NMR-based metabolomic profiling of obesity using leptin-deficient ob/ob mouse has shown metabolic changes and gender variations in obesity. Future work is needed in human samples to validate the biomarker signatures of obesity identified in our study and to characterize leptin-deficient obesity, which can be further used for anti-obesity drug development.

Supporting Information

Figure S1.

Full overlaid spectra of representative 1H NMR spectra of urine sample from lean and obese mice (Figure 1A). Samples from male and female mice are shown on the top and bottom, respectively. Black and red spectra are overlaid and shown for lean mice and ob/ob mice, respectively. Since there are almost no peaks in the upfield region of the spectra for serum samples and all the peaks are shown in Figure 1B, full view of the spectra are not shown here.

https://doi.org/10.1371/journal.pone.0075998.s001

(TIF)

Figure S2.

PCA score plots obtained from the 1H NMR spectra of urine (A) and serum (B) samples from lean and obese mice. All the mice including males and females, male mice, and female mice are shown on the left, middle, and right, respectively.

https://doi.org/10.1371/journal.pone.0075998.s002

(TIF)

Figure S3.

Variable importance in the projection (VIP) plots of lean vs. obese mice obtained from OPLS-DA with a threshold of 1.0. All the mice including males and females, males, and females are shown on the top, middle, and down, respectively, for urine (A) and serum (B) samples.

https://doi.org/10.1371/journal.pone.0075998.s003

(TIF)

Figure S4.

Variable importance in the projection (VIP) plots of male vs. female mice obtained from OPLS-DA with a threshold of 1.0. The urine (A) and serum (B) samples from the obese (top) and lean (down) mice are shown.

https://doi.org/10.1371/journal.pone.0075998.s004

(TIF)

Figure S5.

Pathway enrichment analysis of urine (A) and serum (B) samples. The 47 and 23 urine metabolites from male and female mice, respectively, and 13 and 19 serum metabolites from male and female mice, respectively, all with p<0.05 in Mann-Whitney t-tests, were used for quantitative pathway enrichment analysis.

https://doi.org/10.1371/journal.pone.0075998.s005

(TIF)

Table S1.

Summary of significantly altered urine metabolites and the related metabolic pathways in ob/ob mice in detail.

https://doi.org/10.1371/journal.pone.0075998.s006

(DOCX)

Table S2.

Summary of significantly altered serum metabolites and the related metabolic pathways in ob/ob mice in detail.

https://doi.org/10.1371/journal.pone.0075998.s007

(DOCX)

Acknowledgments

We thank D.H. Choi for technical assistance of animal.

Author Contributions

Conceived and designed the experiments: CHL GSH SWC. Performed the experiments: EYW SWK. Analyzed the data: EYW MKY YJ HWB SWC. Contributed reagents/materials/analysis tools: DL SGP CHL GSH SWC. Wrote the paper: EYW MKY CHL GSH SWC.

References

  1. 1. WHO website, Available: http://www.who.int/features/factfiles/obesity/en/.Accessed 2013 Aug 27.
  2. 2. Conway B, Rene A (2004) Obesity as a disease: no lightweight matter. Obes Rev 5: 145–151.
  3. 3. Zhang Y, Proenca R, Maffei M, Barone M, Leopold L, et al. (1994) Positional cloning of the mouse obese gene and its human homologue. Nature 372: 425–432.
  4. 4. Friedman JM, Halaas JL (1998) Leptin and the regulation of body weight in mammals. Nature 395: 763–770.
  5. 5. Zhang F, Basinski MB, Beals JM, Briggs SL, Churgay LM, et al. (1997) Crystal structure of the obese protein leptin-E100. Nature 387: 206–209.
  6. 6. Montague CT, Farooqi IS, Whitehead JP, Soos MA, Rau H, et al. (1997) Congenital leptin deficiency is associated with severe early-onset obesity in humans. Nature 387: 903–908.
  7. 7. Strobel A, Issad T, Camoin L, Ozata M, Strosberg AD (1998) A leptin missense mutation associated with hypogonadism and morbid obesity. Nat Genet 18: 213–215.
  8. 8. Fischer-Posovszky P, von Schnurbein J, Moepps B, Lahr G, Strauss G, et al. (2010) A new missense mutation in the leptin gene causes mild obesity and hypogonadism without affecting T cell responsiveness. J Clin Endocrinol Metab 95: 2836–2840.
  9. 9. Clement K, Vaisse C, Lahlou N, Cabrol S, Pelloux V, et al. (1998) A mutation in the human leptin receptor gene causes obesity and pituitary dysfunction. Nature 392: 398–401.
  10. 10. Farooqi IS, Matarese G, Lord GM, Keogh JM, Lawrence E, et al. (2002) Beneficial effects of leptin on obesity, T cell hyporesponsiveness, and neuroendocrine/metabolic dysfunction of human congenital leptin deficiency. J Clin Invest 110: 1093–1103.
  11. 11. Nicholson JK, Lindon JC (2008) Systems biology: Metabonomics. Nature 455: 1054–1056.
  12. 12. Connor SC, Hansen MK, Corner A, Smith RF, Ryan TE (2010) Integration of metabolomics and transcriptomics data to aid biomarker discovery in type 2 diabetes. Mol Biosyst 6: 909–921.
  13. 13. Spratlin JL, Serkova NJ, Eckhardt SG (2009) Clinical applications of metabolomics in oncology: a review. Clin Cancer Res 15: 431–440.
  14. 14. Martin FP, Wang Y, Sprenger N, Yap IK, Lundstedt T, et al. (2008) Probiotic modulation of symbiotic gut microbial-host metabolic interactions in a humanized microbiome mouse model. Mol Syst Biol 4: 157.
  15. 15. Serkova NJ, Jackman M, Brown JL, Liu T, Hirose R, et al. (2006) Metabolic profiling of livers and blood from obese Zucker rats. J Hepatol 44: 956–962.
  16. 16. Waldram A, Holmes E, Wang Y, Rantalainen M, Wilson ID, et al. (2009) Top-down systems biology modeling of host metabotype-microbiome associations in obese rodents. J Proteome Res 8: 2361–2375.
  17. 17. Zhao LC, Zhang XD, Liao SX, Gao HC, Wang HY, et al. (2010) A metabonomic comparison of urinary changes in Zucker and GK rats. J Biomed Biotechnol 2010: 431894.
  18. 18. Schirra HJ, Anderson CG, Wilson WJ, Kerr L, Craik DJ, et al. (2008) Altered metabolism of growth hormone receptor mutant mice: a combined NMR metabonomics and microarray study. PLoS One 3: e2764.
  19. 19. Duggan GE, Hittel DS, Sensen CW, Weljie AM, Vogel HJ, et al. (2011) Metabolomic response to exercise training in lean and diet-induced obese mice. J Appl Physiol 110: 1311–1318.
  20. 20. Shearer J, Duggan G, Weljie A, Hittel DS, Wasserman DH, et al. (2008) Metabolomic profiling of dietary-induced insulin resistance in the high fat-fed C57BL/6J mouse. Diabetes Obes Metab 10: 950–958.
  21. 21. Li H, Xie Z, Lin J, Song H, Wang Q, et al. (2008) Transcriptomic and metabonomic profiling of obesity-prone and obesity-resistant rats under high fat diet. J Proteome Res 7: 4775–4783.
  22. 22. Larsen FH, van den Berg F, Engelsen SB (2006) An exploratory chemometric study of 1H NMR spectra of table wines. Journal of Chemometrics 20: 198–208.
  23. 23. Eriksson L, Johansson E, Kettaneh-Wold N, Trygg J, Wikström C, et al.. (2006) Multi- and megavariate data analysis. Basic Principles and Applications.: Umetrics AB.
  24. 24. Xia J, Wishart DS (2010) MSEA: a web-based tool to identify biologically meaningful patterns in quantitative metabolomic data. Nucleic Acids Res 38: W71–77.
  25. 25. Pelleymounter MA, Cullen MJ, Baker MB, Hecht R, Winters D, et al. (1995) Effects of the obese gene product on body weight regulation in ob/ob mice. Science 269: 540–543.
  26. 26. Golay A, Swislocki AL, Chen YD, Jaspan JB, Reaven GM (1986) Effect of obesity on ambient plasma glucose, free fatty acid, insulin, growth hormone, and glucagon concentrations. J Clin Endocrinol Metab 63: 481–484.
  27. 27. Dakshinamurti K, Chauhan J (1989) Biotin. Vitam Horm 45: 337–384.
  28. 28. Brosnan JT, Brosnan ME (2010) Creatine metabolism and the urea cycle. Mol Genet Metab 100 Suppl 1S49–52.
  29. 29. Zyromski NJ, Mathur A, Gowda GA, Murphy C, Swartz-Basile DA, et al. (2009) Nuclear magnetic resonance spectroscopy-based metabolomics of the fatty pancreas: implicating fat in pancreatic pathology. Pancreatology 9: 410–419.
  30. 30. Gipson GT, Tatsuoka KS, Ball RJ, Sokhansanj BA, Hansen MK, et al. (2008) Multi-platform investigation of the metabolome in a leptin receptor defective murine model of type 2 diabetes. Mol Biosyst 4: 1015–1023.
  31. 31. Hay N, Sonenberg N (2004) Upstream and downstream of mTOR. Genes Dev 18: 1926–1945.
  32. 32. Kimball SR, Jefferson LS (2001) Regulation of protein synthesis by branched-chain amino acids. Curr Opin Clin Nutr Metab Care 4: 39–43.
  33. 33. Ruderman NB (1975) Muscle amino acid metabolism and gluconeogenesis. Annu Rev Med 26: 245–258.
  34. 34. Holecek M (2001) The BCAA-BCKA cycle: its relation to alanine and glutamine synthesis and protein balance. Nutrition 17: 70.
  35. 35. Layman DK, Walker DA (2006) Potential importance of leucine in treatment of obesity and the metabolic syndrome. J Nutr 136: 319S–323S.
  36. 36. Lynch CJ (2001) Role of leucine in the regulation of mTOR by amino acids: revelations from structure-activity studies. J Nutr 131: 861S–865S.
  37. 37. Calles-Escandon J, Cunningham J, Felig P (1984) The plasma amino acid response to cafeteria feeding in the rat: influence of hyperphagia, sucrose intake, and exercise. Metabolism 33: 364–368.
  38. 38. van Doorn M, Vogels J, Tas A, van Hoogdalem EJ, Burggraaf J, et al. (2007) Evaluation of metabolite profiles as biomarkers for the pharmacological effects of thiazolidinediones in Type 2 diabetes mellitus patients and healthy volunteers. Br J Clin Pharmacol 63: 562–574.
  39. 39. Garlick PJ (2005) The role of leucine in the regulation of protein metabolism. J Nutr 135: 1553S–1556S.
  40. 40. Xie B, Waters MJ, Schirra HJ (2012) Investigating potential mechanisms of obesity by metabolomics. J Biomed Biotechnol 2012: 805683.
  41. 41. Sishi B, Loos B, Ellis B, Smith W, du Toit EF, et al. (2011) Diet-induced obesity alters signalling pathways and induces atrophy and apoptosis in skeletal muscle in a prediabetic rat model. Exp Physiol 96: 179–193.
  42. 42. Macotela Y, Emanuelli B, Bang AM, Espinoza DO, Boucher J, et al. (2011) Dietary leucine–an environmental modifier of insulin resistance acting on multiple levels of metabolism. PLoS One 6: e21187.
  43. 43. Hinault C, Mothe-Satney I, Gautier N, Lawrence JC Jr, Van Obberghen E (2004) Amino acids and leucine allow insulin activation of the PKB/mTOR pathway in normal adipocytes treated with wortmannin and in adipocytes from db/db mice. FASEB J 18: 1894–1896.
  44. 44. Fu WJ, Haynes TE, Kohli R, Hu J, Shi W, et al. (2005) Dietary L-arginine supplementation reduces fat mass in Zucker diabetic fatty rats. J Nutr 135: 714–721.
  45. 45. McKnight JR, Satterfield MC, Jobgen WS, Smith SB, Spencer TE, et al. (2010) Beneficial effects of L-arginine on reducing obesity: potential mechanisms and important implications for human health. Amino Acids 39: 349–357.
  46. 46. Jobgen W, Meininger CJ, Jobgen SC, Li P, Lee MJ, et al. (2009) Dietary L-arginine supplementation reduces white fat gain and enhances skeletal muscle and brown fat masses in diet-induced obese rats. J Nutr 139: 230–237.
  47. 47. Wu G, Bazer FW, Davis TA, Kim SW, Li P, et al. (2009) Arginine metabolism and nutrition in growth, health and disease. Amino Acids 37: 153–168.
  48. 48. Oberbach A, Bluher M, Wirth H, Till H, Kovacs P, et al. (2011) Combined proteomic and metabolomic profiling of serum reveals association of the complement system with obesity and identifies novel markers of body fat mass changes. J Proteome Res 10: 4769–4788.
  49. 49. Backman L, Hallberg D, Kallner A (1975) Amino acid pattern in plasma before and after jejuno-ileal shunt operation for obesity. Scand J Gastroenterol 10: 811–816.
  50. 50. Piloquet H, Ferchaud-Roucher V, Duengler F, Zair Y, Maugere P, et al. (2003) Insulin effects on acetate metabolism. Am J Physiol Endocrinol Metab 285: E561–565.
  51. 51. DeVilliers DC Jr, Dixit PK, Lazarow A (1966) Citrate metabolism in diabetes. I. Plasma citrate in alloxan-diabetic rats and in clinical diabetes. Metabolism 15: 458–465.
  52. 52. Randle PJ (1998) Regulatory interactions between lipids and carbohydrates: the glucose fatty acid cycle after 35 years. Diabetes Metab Rev 14: 263–283.
  53. 53. Cha YS (2008) Effects of L-carnitine on obesity, diabetes, and as an ergogenic aid. Asia Pac J Clin Nutr 17 Suppl 1306–308.
  54. 54. Kim HJ, Kim JH, Noh S, Hur HJ, Sung MJ, et al. (2011) Metabolomic analysis of livers and serum from high-fat diet induced obese mice. J Proteome Res 10: 722–731.
  55. 55. Kim JY, Park JY, Kim OY, Ham BM, Kim HJ, et al. (2010) Metabolic profiling of plasma in overweight/obese and lean men using ultra performance liquid chromatography and Q-TOF mass spectrometry (UPLC-Q-TOF MS). J Proteome Res 9: 4368–4375.
  56. 56. Walker JB (1961) Metabolic control of creatine biosynthesis. II. Restoration of transamidinase activity following creatine repression. J Biol Chem 236: 493–498.
  57. 57. Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, et al. (2006) An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 444: 1027–1031.
  58. 58. Delzenne NM, Neyrinck AM, Backhed F, Cani PD (2011) Targeting gut microbiota in obesity: effects of prebiotics and probiotics. Nat Rev Endocrinol 7: 639–646.
  59. 59. Jung JY, Kim IY, Kim YN, Kim JS, Shin JH, et al. (2012) 1H NMR-based metabolite profiling of diet-induced obesity in a mouse mode. BMB Rep 45: 419–424.
  60. 60. Dumas ME, Barton RH, Toye A, Cloarec O, Blancher C, et al. (2006) Metabolic profiling reveals a contribution of gut microbiota to fatty liver phenotype in insulin-resistant mice. Proc Natl Acad Sci U S A 103: 12511–12516.
  61. 61. Koeth RA, Wang Z, Levison BS, Buffa JA, Org E, et al. (2013) Intestinal microbiota metabolism of l-carnitine, a nutrient in red meat, promotes atherosclerosis. Nat Med 19: 576–585.
  62. 62. Nicholson JK, Holmes E, Wilson ID (2005) Gut microorganisms, mammalian metabolism and personalized health care. Nat Rev Microbiol 3: 431–438.
  63. 63. Lees HJ, Swann JR, Wilson ID, Nicholson JK, Holmes E (2013) Hippurate: The Natural History of a Mammalian-Microbial Cometabolite. J Proteome Res.
  64. 64. Williams RE, Lenz EM, Evans JA, Wilson ID, Granger JH, et al. (2005) A combined (1)H NMR and HPLC-MS-based metabonomic study of urine from obese (fa/fa) Zucker and normal Wistar-derived rats. J Pharm Biomed Anal 38: 465–471.
  65. 65. Salek RM, Maguire ML, Bentley E, Rubtsov DV, Hough T, et al. (2007) A metabolomic comparison of urinary changes in type 2 diabetes in mouse, rat, and human. Physiol Genomics 29: 99–108.
  66. 66. Calvani R, Miccheli A, Capuani G, Tomassini Miccheli A, Puccetti C, et al. (2010) Gut microbiome-derived metabolites characterize a peculiar obese urinary metabotype. Int J Obes (Lond) 34: 1095–1098.
  67. 67. Smith EA, Macfarlane GT (1996) Enumeration of human colonic bacteria producing phenolic and indolic compounds: effects of pH, carbohydrate availability and retention time on dissimilatory aromatic amino acid metabolism. J Appl Bacteriol 81: 288–302.
  68. 68. Macfarlane GT, Gibson GR, Beatty E, Cummings JH (1992) Estimation of short-chain fatty acid production from protein by human intestinal bacteria based on branched-chain fatty acid measurements. FEMS Microbiology Letters 101: 81–88.
  69. 69. Li M, Wang B, Zhang M, Rantalainen M, Wang S, et al. (2008) Symbiotic gut microbes modulate human metabolic phenotypes. Proc Natl Acad Sci U S A 105: 2117–2122.
  70. 70. Sun J, Schnackenberg LK, Holland RD, Schmitt TC, Cantor GH, et al. (2008) Metabonomics evaluation of urine from rats given acute and chronic doses of acetaminophen using NMR and UPLC/MS. J Chromatogr B Analyt Technol Biomed Life Sci 871: 328–340.
  71. 71. Satoh H (1994) Cardioprotective actions of taurine against intracellular and extracellular calcium-induced effects. Adv Exp Med Biol 359: 181–196.
  72. 72. Zhang M, Bi LF, Fang JH, Su XL, Da GL, et al. (2004) Beneficial effects of taurine on serum lipids in overweight or obese non-diabetic subjects. Amino Acids 26: 267–271.
  73. 73. Kochhar S, Jacobs DM, Ramadan Z, Berruex F, Fuerholz A, et al. (2006) Probing gender-specific metabolism differences in humans by nuclear magnetic resonance-based metabonomics. Anal Biochem 352: 274–281.
  74. 74. Dey J, Creighton A, Lindberg JS, Fuselier HA, Kok DJ, et al. (2002) Estrogen replacement increased the citrate and calcium excretion rates in postmenopausal women with recurrent urolithiasis. J Urol 167: 169–171.
  75. 75. D’Eon T, Braun B (2002) The roles of estrogen and progesterone in regulating carbohydrate and fat utilization at rest and during exercise. J Womens Health Gend Based Med 11: 225–237.
  76. 76. Laaksonen DE, Niskanen L, Punnonen K, Nyyssonen K, Tuomainen TP, et al. (2004) Testosterone and sex hormone-binding globulin predict the metabolic syndrome and diabetes in middle-aged men. Diabetes Care 27: 1036–1041.
  77. 77. Stellato RK, Feldman HA, Hamdy O, Horton ES, McKinlay JB (2000) Testosterone, sex hormone-binding globulin, and the development of type 2 diabetes in middle-aged men: prospective results from the Massachusetts male aging study. Diabetes Care 23: 490–494.
  78. 78. Grossmann M, Thomas MC, Panagiotopoulos S, Sharpe K, Macisaac RJ, et al. (2008) Low testosterone levels are common and associated with insulin resistance in men with diabetes. J Clin Endocrinol Metab 93: 1834–1840.
  79. 79. Gale EA, Gillespie KM (2001) Diabetes and gender. Diabetologia 44: 3–15.
  80. 80. Cao J, Chen TM, Hao WJ, Li J, Liu L, et al. (2012) Correlation between sex hormone levels and obesity in the elderly male. Aging Male 15: 85–89.
  81. 81. Plumb R, Granger J, Stumpf C, Wilson ID, Evans JA, et al. (2003) Metabonomic analysis of mouse urine by liquid-chromatography-time of flight mass spectrometry (LC-TOFMS): detection of strain, diurnal and gender differences. Analyst 128: 819–823.
  82. 82. Slupsky CM, Rankin KN, Wagner J, Fu H, Chang D, et al. (2007) Investigations of the effects of gender, diurnal variation, and age in human urinary metabolomic profiles. Anal Chem 79: 6995–7004.
  83. 83. Bae KH, Kim WK, Lee SC (2012) Involvement of protein tyrosine phosphatases in adipogenesis: new anti-obesity targets? BMB Rep 45: 700–706.