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

Potential Biomarkers of Fatigue Identified by Plasma Metabolome Analysis in Rats

  • Satoshi Kume,

    Affiliation Cellular Function Imaging Team, Division of Bio-function Dynamics Imaging, RIKEN Center for Life Science Technologies, Kobe, Hyogo, Japan

  • Masanori Yamato,

    Affiliation Cellular Function Imaging Team, Division of Bio-function Dynamics Imaging, RIKEN Center for Life Science Technologies, Kobe, Hyogo, Japan

  • Yasuhisa Tamura,

    Affiliation Cellular Function Imaging Team, Division of Bio-function Dynamics Imaging, RIKEN Center for Life Science Technologies, Kobe, Hyogo, Japan

  • Guanghua Jin,

    Affiliation Cellular Function Imaging Team, Division of Bio-function Dynamics Imaging, RIKEN Center for Life Science Technologies, Kobe, Hyogo, Japan

  • Masayuki Nakano,

    Affiliations Cellular Function Imaging Team, Division of Bio-function Dynamics Imaging, RIKEN Center for Life Science Technologies, Kobe, Hyogo, Japan, Department of Physiology, Osaka City University Graduate School of Medicine, Osaka, Japan

  • Yukiharu Miyashige,

    Affiliation Cellular Function Imaging Team, Division of Bio-function Dynamics Imaging, RIKEN Center for Life Science Technologies, Kobe, Hyogo, Japan

  • Asami Eguchi,

    Affiliation Cellular Function Imaging Team, Division of Bio-function Dynamics Imaging, RIKEN Center for Life Science Technologies, Kobe, Hyogo, Japan

  • Yoshiyuki Ogata,

    Affiliation Graduate School of Life and Environmental Sciences, Osaka Prefecture University, Osaka, Japan

  • Nobuhito Goda,

    Affiliation Department of Life Science and Medical Bio-Science, School of Advanced Science and Engineering, Waseda University, Tokyo, Japan

  • Kazuhiro Iwai,

    Affiliation Department of Molecular and Cellular Physiology, Graduate School of Medicine, Kyoto University, Kyoto, Japan

  • Emi Yamano,

    Affiliation Department of Physiology, Osaka City University Graduate School of Medicine, Osaka, Japan

  • Yasuyoshi Watanabe,

    Affiliations Department of Physiology, Osaka City University Graduate School of Medicine, Osaka, Japan, Pathophysiological and Health Science Team, RIKEN Center for Life Science Technologies, Kobe, Hyogo, Japan

  • Tomoyoshi Soga,

    Affiliation Institute for Advanced Biosciences, Keio University, Tsuruoka, Yamagata, Japan

  • Yosky Kataoka

    kataokay@riken.jp

    Affiliations Cellular Function Imaging Team, Division of Bio-function Dynamics Imaging, RIKEN Center for Life Science Technologies, Kobe, Hyogo, Japan, Department of Physiology, Osaka City University Graduate School of Medicine, Osaka, Japan

Abstract

In the present study, prior to the establishment of a method for the clinical diagnosis of chronic fatigue in humans, we validated the utility of plasma metabolomic analysis in a rat model of fatigue using capillary electrophoresis-mass spectrometry (CE-MS). In order to obtain a fatigued animal group, rats were placed in a cage filled with water to a height of 2.2 cm for 5 days. A food-restricted group, in which rats were limited to 10 g/d of food (around 50% of the control group), was also assessed. The food-restricted group exhibited weight reduction similar to that of the fatigued group. CE-MS measurements were performed to evaluate the profile of food intake-dependent metabolic changes, as well as the profile in fatigue loading, resulting in the identification of 48 metabolites in plasma. Multivariate analyses using hierarchical clustering and principal component analysis revealed that the plasma metabolome in the fatigued group showed clear differences from those in the control and food-restricted groups. In the fatigued group, we found distinctive changes in metabolites related to branched-chain amino acid metabolism, urea cycle, and proline metabolism. Specifically, the fatigued group exhibited significant increases in valine, leucine, isoleucine, and 2-oxoisopentanoate, and significant decreases in citrulline and hydroxyproline compared with the control and food-restricted groups. Plasma levels of total nitric oxide were increased in the fatigued group, indicating systemic oxidative stress. Further, plasma metabolites involved in the citrate cycle, such as cis-aconitate and isocitrate, were reduced in the fatigued group. The levels of ATP were significantly decreased in the liver and skeletal muscle, indicative of a deterioration in energy metabolism in these organs. Thus, this comprehensive metabolic analysis furthered our understanding of the pathophysiology of fatigue, and identified potential diagnostic biomarkers based on fatigue pathophysiology.

Introduction

Fatigue is best defined as difficulty in initiating or sustaining voluntary activities [1] and is thought to be accompanied by deterioration of performance [2]. Half of the general population in modern society currently experiences chronic and complex (a combination of physical and mental) fatigue caused by continual stress and prolonged deficiency of rest or sleep [3]. Besides the chronic and complex fatigue seen in our daily life, there are patients with chronic fatigue syndrome (CFS) exhibiting persistent fatigue which is substantially unrelieved by rest, and accompanied by other symptoms such as circadian rhythm sleep disorder and cognitive dysfunction for a minimum of 6 months [4]. CFS is thought to have a worldwide prevalence of 0.4–1.0%. The pathogenesis of CFS remains incompletely understood, and the objective diagnostics based on the pathophysiology remains to be developed.

Metabolic alterations induced by fatigue have been previously reported. After exhaustive or sustained exercise in humans, blood levels of branched-chain amino acids (BCAAs) are decreased compared to pre-exercise levels, whereas free tryptophan levels are increased [5,6]. In addition, mental fatigue loading in humans leads to a decrease in blood levels of the BCAAs, tyrosine, cysteine, methionine, lysine, and arginine [7].

In order to investigate the pathophysiology of fatigue in multiple organs, we have established a complex fatigue animal model where rats are exposed to relatively long-lasting stress and partial sleep deprivation, which humans often experience in their daily lives [8]. In this model, rats were housed in a cage filled with water to a depth of 2.2 cm for 5 days. The fatigued rats showed a significant decrease in swimming time in a weight-loaded forced swimming test, indicative of fatigue [2,8]. The fatigued rats exhibited dysfunction in multiple systems, including the neuro-immuno-endocrine system; specifically, they exhibited increased turnover of 5-hydroxytryptamine in the central nervous system [8], chronic up-regulation of gene expression for the precursor peptides of adrenocorticotropin, β-endorphin, and α-melanocyte stimulating hormone (α-MSH) in the pituitary gland [9], an increase in the peripheral blood level of α-MSH released from the intermediate lobe of hypophysis [10], and an accumulation and an activation of microglial cells in the dorsal horn of spinal cord [11]. Furthermore, levels of BCAAs, such as valine, leucine, and isoleucine, were increased in the plasma, skeletal muscle, liver, and brain, whereas levels of other amino acids, such as glutamine, phenylalanine, tyrosine, arginine, and threonine, were decreased in fatigued rats [2]. The decrease in plasma arginine suggested the facilitation of nitric oxide (NO) production by NO synthase in fatigued rats. Rats experiencing food restriction, which show a weight reduction comparable to fatigued rats, exhibited a different amino acid profile from that of fatigued rats [2]. These findings indicated that the pathophysiology of fatigue includes abnormalities in the neuro-immuno-endocrine system, as well as in metabolism. Although those reports have implied that blood amino acids levels could serve as biomarkers for the diagnosis of physical and/or mental fatigue, more comprehensive metabolite profiling is expected to result in the identification of optimal biomarkers based on the pathophysiology of fatigue.

Recently, metabolome analysis has become a rapidly evolving analytical technology for the comprehensive identification and quantification of endogenous metabolites in samples collected from blood or tissues [12]. This method enables the identification of metabolites or metabolic pathways involved in disease pathophysiology. The efficacy of this experimental approach was demonstrated by the discovery of biomarkers for aging [13], drug–induced gastric injury [14], dilated cardiomyopathy [15], and type 2 diabetes [16].

In the present study, prior to engaging in clinical trials using metabolomic analysis to identify diagnostic biomarkers of chronic fatigue in humans, we conducted a plasma metabolomic analysis of a rat model of fatigue using capillary electrophoresis-mass spectrometry (CE-MS). The objective of this study was to evaluate the utility of a comprehensive analysis of metabolites for identifying potential diagnostic biomarkers of chronic fatigue. In order to identify food intake-independent metabolic changes, plasma samples from fatigued animals were compared to those obtained from a food-restricted group, which exhibited a weight reduction similar to that of the fatigued group. Metabolites potently influenced by food restriction were eliminated from the biomarker candidates for chronic fatigue. Multivariate analyses were employed for determination of potential biomarkers reflecting the pathophysiology of fatigue. Furthermore, we investigated adenosine triphosphate (ATP) levels in organs and erythrocytes, as well as the levels of plasma NO, which are closely associated with metabolic profiles.

Materials and Methods

Study design for a fatigue-loaded animal model

Eight-week old male Sprague–Dawley rats (Shizuoka Laboratory Animal Cooperative; Shizuoka, Japan) were used in this study. The animals were housed in a cage with a raised mesh base under constant environmental conditions (room temperature, 22–23°C; relative humidity, 50%-60%) and a 12-h light-dark cycle (08:00/20:00). Prior to the experiment, food and water were provided ad libitum. All experimental protocols were approved by the Ethics Committee on Animal Care and Use of the RIKEN Center for Life Science Technologies (MAH19–01–13) and were performed in accordance with the Principles of Laboratory Animal Care (NIH publication No. 85–23, revised 1985). All efforts were made to minimize animal suffering and the number of animals used for the studies.

Rats were randomly divided into three groups; control, food-restriction, and fatigued groups. In order to produce an animal model of relatively long-lasting fatigue, rats underwent deprivation of rest and sleep by being housed in a cage filled with water (23±1°C) to a height of 2.2 cm for 5 days, as previously reported by Jin et al [2]. Because the rats exhibited a reduction in body weight during the fatigue period, we included a food-restricted group that produced a degree of weight reduction similar to the fatigued rats. The food-restricted group was housed in normal cages for 5 days and their food intake was restricted to 10 g/day, as described previously [2]. A non-treated control group was also included.

Plasma and tissue preparations

Rats were deeply anaesthetized with intraperitoneal injections of sodium pentobarbital (100 mg/kg), and whole blood samples were obtained from the abdominal aorta. The brain (cerebral cortex), liver, and skeletal muscle were removed after perfusion with ice-cold phosphate buffered saline. The blood samples (1.5 ml) were collected into tubes containing 35 μl of 0.06 g/ml (w/v) EDTA and plasma was obtained following centrifugation at 1000 × g at 4°C for 10 min.

Capillary electrophoresis time-of-flight mass spectrometry (CE-TOFMS) analysis

For extraction of plasma metabolites, 50 μl plasma samples were vortexed after adding 450 μl of methanol containing internal standards (20 μM each of L-methionine sulfone, D-Camphol-10-sulfonic acid and 2-morpholinoethanesulfonic acid). The samples were then mixed with 200 μl Milli-Q water and 500 μl chloroform and centrifuged at 4,600 × g at 4°C for 5 min. Subsequently, the aqueous solution was centrifugally filtered through a 5-kDa cutoff filter (Millipore) at 9,100 × g at 4°C for 3h. The filtered samples (300 μl) were centrifugally concentrated at 35°C for 2h. The dried samples were dissolved in 25 μl Milli-Q water containing reference compounds (200 μM each of 3-aminopyrrolidine and trimesate) for capillary electrophoresis time-of-flight mass spectrometry (CE-TOFMS) analysis [17,18].

CE-TOFMS analysis was performed using an Agilent CE capillary electrophoresis system (Agilent Technologies, Waldbronn, Germany), an Agilent G3250AA liquid chromatography/mass selective detector time-of-flight system (LC/MSD TOF; Agilent Technologies, Palo Alto, CA, USA), an Agilent 1100 series binary high-performance liquid chromatography (HPLC) pump with a G1603A Agilent CE-MS adapter and a G1607A Agilent capillary electrophoresis electrospray ionization mass spectrometry (CE-ESI-MS) sprayer kit. For anion analysis, an Agilent G710060041 platinum ESI needle was used. G2201AA Agilent ChemStation software for CE and Analyst QS Agilent TOFMS software were used for system control and data acquisition. Capillary electrophoresis tandem mass spectrometry (CE-MS/MS) analyses for compound identification were performed on a QStar XL Hybrid liquid chromatography tandem mass spectrometry system (LC-MS/MS; Applied Biosystems, Foster City, CA. USA) connected to an Agilent CE instrument.

CE-TOFMS conditions for cationic metabolite analysis

Separations were carried out in a fused silica capillary (50 mm internal diameter and 100 cm total length) filled with 1 M formic acid as the electrolyte [1921]. Approximately 3 nl of sample solution were injected at 50 mbar for 3 sec, and a 30 kV voltage was applied. The capillary temperature was maintained at 20°C and the sample tray was cooled below 5°C. Methanol/water (50% v/v) containing 0.1 μM Hexakis (2,2-difluoroethoxy)phosphazene was delivered as the sheath liquid at 10 μl/min. CE-ESI-TOFMS was operated in the positive ion mode and the capillary voltage was set at 4 kV. A flow rate of heated dry nitrogen gas (heater temperature, 300°C) was maintained at 10 psig. In TOFMS, the fragmentor, skimmer, and octapole radio frequency voltages (Oct RFV) were set at 75V, 50V, and 125V, respectively. Automatic recalibration of each acquired spectrum was achieved using the masses of reference standards ([13C isotopic ion of a protonated methanol dimer (2MeOH+H)]+, m/z 66.0631; and [Hexakis (2,2-difluoroethoxy)phosphazene + H]+, m/z 622.0290). Exact mass data were acquired at a rate of 1.5 spectra/sec over a 50–1,000 m/z range.

CE-TOFMS conditions for anionic metabolite analysis

A commercially available COSMO(+), chemically coated with cationic polymer, capillary (50 mm internal diameter and 110 cm total length; Nacalai Tesque, Kyoto, Japan) was used as the separation capillary [22]. A 50 mM ammonium acetate solution (pH 8.5) was used as the electrolyte solution for CE separation. Sample solutions (30 nl) were injected at 50 mbar for 30 sec, and a 30 kV voltage was applied. Methanol/water (50% v/v) containing 5 mM ammonium acetate and 0.1 μM Hexakis was delivered as the sheath liquid at 10 μl/min. ESI-TOFMS was conducted in the negative ion mode with the capillary voltage set to 3.5 kV. For TOFMS, the fragmentor, skimmer, and Oct RFV voltages were set at 100V, 50V, and 200V, respectively. Automatic recalibration of each acquired spectrum was performed using the masses of reference standards ([13C isotopic ion of deprotonated acetic acid dimer (2CH3COOH-H)], m/z 120.0383; and [Hexakis + deprotonated acetic acid (CH3COOH-H)], m/z 680.0355). Exact mass data were acquired at a rate of 1.5 spectra/sec over a 50–1,000 m/z range. Other conditions were identical to those used in cationic metabolite analysis.

Measurement of adenosine triphosphate (ATP)

ATP levels of samples were measured using a firefly bioluminescence assay kit (AMERIC-ATP kit; Wako Pure Chemical Industries, Osaka, Japan) [23]. The tissues (brain, liver, and skeletal muscle) and erythrocytes were homogenized immediately in phenol (3 ml) after isolation. The extracted sample (1 ml) was shaken with chloroform (200 μl) and centrifuged at 13,800 × g at 4°C for 5 min. The upper aqueous phase was diluted 2,000-fold for the brain and liver and 10,000-fold for erythrocytes and skeletal muscle with distilled water. The diluted extract (10 μl) was then injected into 90 μl of a luciferin/luciferase mixture and the bioluminescence product was immediately measured using a luminometer (GloMax-96 Microplate Luminometer; Promega, Tokyo, Japan). ATP levels in each sample were calculated from a calibration curve made with ATP standards (0–10–7 M).

Measurement of Nitric Oxide (NO) in Plasma

Nitric oxide production (total concentrations of nitrate and nitrite) was measured using a QuantiChrom Nitric Oxide Assay Kit (Bioassay Systems, Hayward, CA, USA), according to the Griess method [24,25]. A sodium nitrite standard curve was prepared (0–150 μM). Deproteination of the plasma was accomplished by mixing 150 μl of plasma with 8 μl ZnSO4, followed by mixing with 8 μl NaOH, vortexing, and centrifuging at 12,000 × g at 4°C for 10 min. Subsequently, 100 μl of each supernatant and the standard were added to 200 μl of the working reagent and incubated at 37°C for 60 min. After centrifugation of the reacted samples at 6,200 × g at 4°C for 5 min, the supernatant was diluted 2-fold with distilled water. Optical density was measured at 550 nm using an Eppendorf BioPhotometer Plus (Eppendorf, Hamburg, Germany). The total nitrite concentrations (μM) were determined according to the sodium nitrite standard curve.

Data processing, bioinformatics, and Statistical analysis

For data processing of the metabolomic data, missing values were imputed with the lower limit of detection for a given metabolite. The data were normalized using the auto-scaling method (mean-centering and variance-scaling).

For data resampling, outliers in each group were identified using a one-sample t-test. The threshold of data removal was set to P-values less than 0.01. The metabolomic data were analyzed with principal component (PC) analysis, and the scores of the primary component of variance (PC1) in each group were tested using the one sample t-test. In this analysis, no outliers were detected. In the ATP and NO data, outliers in the data were identified using a one-sample t-test and removed prior to statistical processing.

Multivariate analyses based on the hierarchical clustering heatmap analysis, PC analysis (PCA), and random forest analysis were performed using the R Packages for Metabolomics Univariate and Multivariate Statistical Analysis [26], randomForest, and gplots on the R statistics platform (R Foundation for Statistical Computing, http://www.r-project.org) [27]. Visualizations of metabolic profiles were performed using VANTED version 2.0.1 [28].

Statistically significant differences were evaluated using a one-way analysis of variance (ANOVA) and Tukey’s honestly significant difference (HSD) post hoc test for comparisons of multiple groups using the R statistics package [27] and MetaboAnalyst 2.0 web portal (http://www.metaboanalyst.ca/MetaboAnalyst/) [29,30]. P < 0.05 and a false discovery rate (FDR) < 0.1 were considered statistically significant.

Results

Profiling of plasma metabolites

The plasma samples in control, food-restricted, and fatigued rats were analyzed with the CE-TOFMS system, and 48 different cationic and anionic metabolites were identified and quantified (Tables 1 and S1). Hierarchical clustering heatmap analysis and principal component (PC) analysis showed characteristic patterns of expression for the metabolites (Figs. 1 and 2). The hierarchical clustering analysis identified three clusters corresponding to the control, food-restricted, and fatigued animal groups. The two-dimensional PCA plot (PC1 vs. PC2) accounted for 45.4% of the total variance (Fig. 2A). The PCA plot revealed that a clear difference was observed in the metabolome profile, indicating that the PC1 axis mainly comprises the variance information of the fatigued group against both the control and food-restricted groups (Fig. 2A). The PC2 axis contained the variance information of the control against the food-restricted group (Fig. 2A). In the PCA loading plot, phenylalanine, valine, leucine, and isoleucine positively contributed to the loadings of PC1, and hydroxyproline, threonine, glycine, and citrulline negatively contributed (Fig. 2B); while ornithine, proline, and tyrosine positively contributed to the loadings of PC2, and inosine, guanosine, and creatine negatively contributed (Fig. 2B).

thumbnail
Fig 1. Two-way hierarchical clustering heatmap of plasma metabolome data.

Each column shows the metabolic pattern of individual animals in the control group (n = 5), food-restricted group (n = 5), and fatigued group (n = 6). The amount of each metabolite in individual samples is expressed as relative value obtained by the auto-scaling method and is represented by the color scheme in which red and blue indicate high and low concentrations of metabolites, respectively.

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

thumbnail
Fig 2. Score and loading plots of principal component analysis of plasma metabolome CE-MS data.

(A) PCA score plot of PC2 versus PC1. The control group (n = 5), food-restricted group (n = 5), and fatigued group (n = 5) are shown as black, green, and red circles, respectively. Black ellipses represent the 90% confidence intervals for each group. (B) PCA loading plot of PC2 versus PC1. The data were analyzed after being mean-centered and variance scaled.

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

Additionally, one-way ANOVA comparison (P < 0.05 and FDR < 0.1) identified 25 metabolites that were differentially expressed among the control, food-restricted, and fatigued groups. Hierarchical clustering analysis was performed on the statistically different metabolites, resulting in the demonstration of clear differences among the three treatment groups (Fig. 3). Further, the metabolites were ranked using the mean decrease in accuracy from the random forest analysis (Fig. 4). The most important metabolites for making a discrimination of the fatigued group among three groups (control, food-restricted, and fatigued) contained phenylalanine, valine, hypotaurine, serine, and hydroxyproline (Fig. 4A), while those between two groups (fatigued and a combined control and food-restricted group) contained hydroxyproline, serine, glycine, phenylalanine, and leucine (Fig. 4B). There was a high degree of consistency between the top-scoring metabolites identified by the random forest analysis (Fig. 4) and those identified by hierarchical clustering (Fig. 3). Most of the identified metabolites were closely related to the TCA cycle (succinate), BCAA metabolism (valine, leucine, isoleucine, and 2-oxoisopentanoate), urea cycle (ornithine, citrulline, and arginine), proline metabolism (hydroxyproline), and other amino acid metabolism. These results indicate that fatigue-loading potently influences metabolic profiles in rat plasma.

thumbnail
Fig 3. Supervised hierarchical clustered heatmap of 25 metabolites identified by one-way ANOVA.

Each column shows the metabolic pattern of individual animals in the control, food-restricted, and fatigued groups. The amount of each metabolite in individual samples is expressed as relative value obtained by the auto-scaling method and is represented by the color scheme in which red and blue indicate high and low concentrations of metabolites, respectively.

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

thumbnail
Fig 4. Random forest analysis of the plasma metabolome.

The mean decrease in accuracy from the random forest analysis was used to rank metabolites according to their prognostic importance for fatigue status. The 15 most important metabolites among three groups (control, food-restricted, and fatigued) or two groups (fatigued and a combined control and food-restricted group) are shown in A and B, respectively.

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

Potential biomarkers of fatigue in plasma metabolites

In the first steps of the TCA cycle, the levels of organic acids in the fatigued group showed a slight reduction in comparison with those in the control group (citrate decreased by 12.1%, cis-aconitate by 7.7%, and isocitrate by 4.5%), while a subset of metabolites in the TCA cycle showed a trend toward overcompensation (2-oxoglutarate increased by 11.2%, succinate by 15.7%, fumarate by 27.0%, and malate by 24.6%) (Fig. 5). The decreases in citrate, cis-aconitate and isocitrate in the fatigued group were confirmed by the LC/MS measurements (significant reduction in cis-aconitate and isocitrate; S1 Fig and S2 Table). In the food-restricted group, the organic acids of the TCA cycle were slightly changed: cis-aconitate decreased by 10.3%, isocitrate by 17.9%, and succinate by 24.7%, while 2-oxoglutarate increased by 16.2%, fumarate by 18.1%, and malate by 11.0% (Fig. 5). Plasma ATP levels in the three groups were below the detection limit of the assay.

thumbnail
Fig 5. Relative concentrations of plasma metabolites related to the TCA cycle, urea cycle and proline metabolism, and BCAA metabolism.

Relative amounts of each metabolite in the control (C, n = 5), food-restricted (R, n = 5), and fatigued (F, n = 6) groups are shown and expressed as a percent of the control group. Data are presented as mean ± S.D. P < 0.05, significantly different than the control group. §P < 0.05, significantly different between the food-restricted and fatigued group.

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

The urea cycle is known to counteract the toxic effects of ammonia in the liver. Plasma levels of metabolites related to the urea cycle were reduced in the fatigued group: ornithine decreased by 14.7%, citrulline by 33.1% (P < 0.01, compared with the control), and arginine by 31.3% (P < 0.01, compared with the control) (Fig. 5). The citrulline levels in the fatigued group were also lower compared to the food-restricted group (Fig. 5). In addition, plasma levels of aspartate and proline in the fatigued group were slightly decreased compared to the control group, while hydroxyproline in the fatigue group was significantly lower than in the control and food-restricted groups (P < 0.01) (Fig. 5). Conversely, the fatigued group showed significant increases in the BCAAs valine, leucine, and isoleucine, as well as 2-oxoisopentanoate, compared with those in both the control and food-restricted groups. Thus, the changes in citrulline, hydroxyproline, and BCAAs in the fatigued group were attributed to fatigue-loading, but not to food restriction. In addition, the purine metabolites inosine and guanosine were significantly increased only in the food-restricted group (Table 1).

Tissue ATP contents

The changes in plasma metabolites involved in the TCA cycle in fatigued rats suggested abnormal energy metabolism during fatigue loading. In order to evaluate energy production and the demand in organs during food restriction or fatigue, we examined tissue ATP levels in the liver, skeletal muscle, brain, and erythrocytes. In the liver, ATP levels in the fatigued group were significantly deceased compared to the control group (P = 0.013) (Fig. 6). In the skeletal muscle, ATP levels in the fatigued group showed a significant decrease compared with those in the control (P < 0.01) and food restricted groups (P < 0.01) (Fig. 6). In the brain and erythrocytes, ATP levels in the fatigued group and food restricted group showed slight reductions compared with that of the control group (Fig. 6). These data indicate that fatigue loading decreases energy production and/or increases energy demand, particularly in the liver and skeletal muscle.

thumbnail
Fig 6. ATP levels in the liver, brain, skeletal muscle, and erythrocytes.

Relative levels of ATP in each tissue in the control (C, n = 5–6), food-restricted (R, n = 5–6), and fatigued (F, n = 5–7) groups are shown and expressed as a percent of the control group. Data are presented as mean ± S.D. P < 0.05, significantly different from the control group. §P < 0.05, significantly different between the food-restricted and fatigued group.

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

Fatigue-induced changes in plasma NOx levels

In order to assess the relationships between nitric oxide (NO) levels and arginine-citrulline balance under the fatigue condition, we measured total concentrations of nitrate and nitrite (NOx) in plasma as NO levels. The total NOx in the plasma of the fatigued group was significantly increased compared with the control (P = 0.01) and food-restricted groups (P < 0.01, Fig. 7). Food-restriction induced a significant reduction in NOx levels compared to control (P < 0.01, Fig. 7). These results indicate that plasma NO levels are increased under the fatigue condition.

thumbnail
Fig 7. NOx content in the plasma.

Plasma NOx content in the control (C, n = 5), food-restricted (R, n = 6), and fatigued (F, n = 6) groups are shown. Data are presented as mean ± S.D. P < 0.05, significantly different from the control group. §P < 0.05, significantly different between the food-restricted and fatigued group.

https://doi.org/10.1371/journal.pone.0120106.g007

Discussion

In this study, we performed a comprehensive determination of metabolic alterations in the plasma of fatigued rats. The candidate fatigue biomarkers were differentiated from those affected by alterations in food-intake, and were subjected to multivariate analysis. The metabolomic analysis revealed that fatigue-induced metabolic changes in the plasma were clearly discriminated from those in the non-treated controls and the food-restricted group (Fig. 2). We previously reported fatigue-induced changes in the levels of several amino acids in the plasma of fatigued rats, including increases in BCAAs (valine, leucine, and isoleucine) and phenylalanine, as well as decreases in tryptophan, methionine, lysine, arginine, histidine, serine, threonine, asparagine, aspartate, glutamine, glutamate, glycine, alanine, and proline [2]. As reported by Lehmann et al. [5], plasma levels of many amino acids such as threonine, serine, arginine, proline, glycine, alanine, valine, ornithine, lysine, histidine, and taurine were significantly decreased after the human extreme exercise in the 1993 Colmar ultra triathlon. These systematical decreases in the amino acids are considered to be due to the high consumption of amino acids for compensating a decrease of energy metabolism within a relatively short time. The differences in changing profiles of metabolites between the human extreme exercise and the present rat fatigue model may reflect different in duration and strength of fatigue-load. Recently, Matsui et al. reported that prolonged exhaustive exercise decreased glycogen levels in brain, skeletal muscle, and liver to less than 50% of pre-exercise levels [31]. We identified changes in the plasma metabolites related to the TCA cycle in the fatigued model that trended toward significance (Fig. 5, S1 Fig, and S2 Table), indicating abnormal energy metabolism in several organs. The levels of ATP in liver and skeletal muscle were significantly decreased in the fatigued group (Fig. 6). These results suggested that fatigue influences multiple pathways in primary metabolism, resulting from increased energy demand and decreased supply, along with decreased levels of glycogen and amino acids. The increased energy demand may be correlated to up-regulation of BCAAs (leucine, isoleucine, and valine) due to increased proteolysis in skeletal muscle, as discussed previously [2].

Multivariate analyses using PCA and random forest analysis identified significant decreases in the levels of citrulline and hydroxyproline in the fatigued group, compared with the control and food-restricted groups (Figs. 2 and 4). Citrulline is a component of the urea cycle, is involved in major liver function detoxification, and is produced from ornithine and carbamoyl phosphate by ornithine transcarbamylase [32]. Citrulline is also produced from arginine by the NO production pathway and is catalysed by nitric oxide synthase, which is reported to be up-regulated by AMP-activated protein kinase in response to ADP/ATP imbalance [33], as well as other signalling inputs [34]. The metabolic balance of citrulline and arginine is known to influence intracellular and extracellular lipid peroxidation levels [35,36]. In fatigued animals, plasma NOx levels were observed to be increased (Fig. 7), and plasma oxidative stress levels were also found to be elevated. Indeed, thiobarbituric acid-reactive lipoperoxide was reported to be increased in the liver tissue of fatigued animals [37]. Furthermore, oxidative stress is known to stimulate modifications of amino acids, such as the cleavage and oxidation of proline [38]. The decrease in plasma hydroxyproline observed in the present study might be related to oxidative modifications of proline in the liver. Taken together, chronic fatigue is thought to disturb the balance of citrulline and arginine due to metabolic stress (ADP/ATP imbalance), which leads to oxidative stress, including peroxidation in the liver.

We identified fatigue-induced changes in primary metabolism, including urea cycle, proline metabolism, and BCAA metabolism, in rat plasma. The influence of food-intake and increased proteolysis was taken into consideration, and citrulline and hydroxyproline were newly identified as metabolic signs reflecting complex fatigue in the rat fatigue model. These metabolites may be promising diagnostic biomarkers for human chronic fatigue and/or CFS. As described in Geigy Scientific Tables edited by Lentner, C. [39], some of plasma metabolites such as hydroxyproline are known to show sexual differences in human. If we apply such metabolites to the diagnostic biomarkers in humans, sexually-segregated plasma levels of the metabolites will be needed to be taken into consideration. Taken together, our results may be helpful for realizing the effective prevention and treatment of chronic fatigue by controlling metabolism in humans.

Supporting Information

S1 Fig. Relative concentrations of metabolites related to the TCA cycle in rat plasma.

Plasma samples were analyzed using a LC-10A series HPLC (Shimadzu, Kyoto, Japan) equipped with an API 5000 triple quadrupole mass spectrometer (AB Sciex, Foster City, CA, USA). Relative concentrations of TCA cycle metabolites (citrate, cis-aconitate, isocitrate, succinate, and malate) in the control group (C, n = 3) and fatigued group (F, n = 3) are shown and expressed as a percent of the control group. Data are presented as mean ± S.D. P < 0.05, significantly different from the control group. The LC/MS measurements of rat plasma samples revealed a trend toward decreased citrate (P = 0.053) and significant decreases in cis-aconitate (P = 0.014) and isocitrate (P = 0.036) levels.

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

(TIF)

S1 Table. Quantitative metabolome data of rat plasma obtained by CE-MS study.

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

(XLS)

S2 Table. Quantitative metabolome data of rat plasma obtained by LC-MS study.

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

(XLS)

Acknowledgments

We thank Mr Tamitoshi Kobayashi (Center for Life Science Technologies, Kobe, Japan.) for the technical assistance. The Sumika Chemical Analysis Service is acknowledged for performing the LC-MS measurements.

Author Contributions

Conceived and designed the experiments: SK YK. Performed the experiments: SK MY YT GJ MN YM AE NG TS YK. Analyzed the data: SK YO NG KI EY YW. Contributed reagents/materials/analysis tools: TS YW YK. Wrote the paper: SK MY YO NG KI EY YW TS YK.

References

  1. 1. Chaudhuri A, Behan PO. Fatigue in neurological disorders. Lancet. 2004; 363: 978–988. pmid:15043967
  2. 2. Jin G, Kataoka Y, Tanaka M, Mizuma H, Nozaki S, Tahara T, et al. Changes in plasma and tissue amino acid levels in an animal model of complex fatigue. Nutrition. 2009; 25: 597–607. pmid:19216057
  3. 3. Walker EA, Katon WJ, Jemelka RP. Psychiatric disorders and medical care utilization among people in the general population who report fatigue. J Gen Intern Med. 1993; 8: 436–440. pmid:8410409
  4. 4. Fukuda K, Straus SE, Hickie I, Sharpe MC, Dobbins JG, Komaroff A. The chronic fatigue syndrome: a comprehensive approach to its definition and study. International Chronic Fatigue Syndrome Study Group. Ann Intern Med. 1994; 121: 953–959. pmid:7978722
  5. 5. Lehmann M, Huonker M, Dimeo F, Heinz N, Gastmann U, Treis N, et al. Serum amino acid concentrations in nine athletes before and after the 1993 Colmar ultra triathlon. Int J Sports Med. 1995; 16: 155–159. pmid:7649705
  6. 6. Blomstrand E, Hassmen P, Ekblom B, Newsholme EA. Administration of branched-chain amino acids during sustained exercise—effects on performance and on plasma concentration of some amino acids. Eur J Appl Physiol Occup Physiol. 1991; 63: 83–88. pmid:1748109
  7. 7. Mizuno K, Tanaka M, Nozaki S, Yamaguti K, Mizuma H, Sasabe T, et al. Mental fatigue-induced decrease in levels of several plasma amino acids. J Neural Transm. 2007; 114: 555–561. pmid:17160370
  8. 8. Tanaka M, Nakamura F, Mizokawa S, Matsumura A, Nozaki S, Watanabe Y. Establishment and assessment of a rat model of fatigue. Neurosci Lett. 2003; 352: 159–162. pmid:14625009
  9. 9. Ogawa T, Kiryu-Seo S, Tanaka M, Konishi H, Iwata N, Saido T, et al. Altered expression of neprilysin family members in the pituitary gland of sleep-disturbed rats, an animal model of severe fatigue. J Neurochem. 2005; 95: 1156–1166. pmid:16271049
  10. 10. Ogawa T, Shishioh-Ikejima N, Konishi H, Makino T, Sei H, Kiryu-Seo S, et al. Chronic stress elicits prolonged activation of alpha-MSH secretion and subsequent degeneration of melanotroph. J Neurochem. 2009; 109: 1389–1399. pmid:19476550
  11. 11. Yasui M, Yoshimura T, Takeuchi S, Tokizane K, Tsuda M, Inoue K, et al. A Chronic fatigue syndrome model demonstrates mechanical allodynia and muscular hyperalgesia via spinal microglial activation. Glia. 2014.
  12. 12. Villas-Boas SG, Mas S, Akesson M, Smedsgaard J, Nielsen J. Mass spectrometry in metabolome analysis. Mass Spectrom Rev. 2005; 24: 613–646. pmid:15389842
  13. 13. Houtkooper RH, Argmann C, Houten SM, Canto C, Jeninga EH, Andreux PA, et al. The metabolic footprint of aging in mice. Sci Rep. 2011; 1: 134. pmid:22355651
  14. 14. Takeuchi K, Ohishi M, Ota S, Suzumura K, Naraoka H, Ohata T, et al. Metabolic profiling to identify potential serum biomarkers for gastric ulceration induced by nonsteroid anti-inflammatory drugs. J Proteome Res. 2013; 12: 1399–1407. pmid:23336624
  15. 15. Maekawa K, Hirayama A, Iwata Y, Tajima Y, Nishimaki-Mogami T, Sugawara S, et al. Global metabolomic analysis of heart tissue in a hamster model for dilated cardiomyopathy. J Mol Cell Cardiol. 2013; 59: 76–85. pmid:23454301
  16. 16. Bain JR, Stevens RD, Wenner BR, Ilkayeva O, Muoio DM, Newgard CB. Metabolomics applied to diabetes research: moving from information to knowledge. Diabetes. 2009; 58: 2429–2443. pmid:19875619
  17. 17. Soga T, Sugimoto M, Honma M, Mori M, Igarashi K, Kashikura K, et al. Serum metabolomics reveals gamma-glutamyl dipeptides as biomarkers for discrimination among different forms of liver disease. J Hepatol. 2011; 55: 896–905. pmid:21334394
  18. 18. Hirayama A, Kami K, Sugimoto M, Sugawara M, Toki N, Onozuka H, et al. Quantitative metabolome profiling of colon and stomach cancer microenvironment by capillary electrophoresis time-of-flight mass spectrometry. Cancer Res. 2009; 69: 4918–4925. pmid:19458066
  19. 19. Soga T, Heiger DN. Amino acid analysis by capillary electrophoresis electrospray ionization mass spectrometry. Anal Chem. 2000; 72: 1236–1241. pmid:10740865
  20. 20. Ishii N, Nakahigashi K, Baba T, Robert M, Soga T, Kanai A, et al. Multiple high-throughput analyses monitor the response of E. coli to perturbations. Science. 2007; 316: 593–597. pmid:17379776
  21. 21. Soga T, Baran R, Suematsu M, Ueno Y, Ikeda S, Sakurakawa T, et al. Differential metabolomics reveals ophthalmic acid as an oxidative stress biomarker indicating hepatic glutathione consumption. J Biol Chem. 2006; 281: 16768–16776. pmid:16608839
  22. 22. Soga T, Igarashi K, Ito C, Mizobuchi K, Zimmermann HP, Tomita M. Metabolomic profiling of anionic metabolites by capillary electrophoresis mass spectrometry. Anal Chem. 2009; 81: 6165–6174. pmid:19522513
  23. 23. Chida J, Yamane K, Takei T, Kido H. An efficient extraction method for quantitation of adenosine triphosphate in mammalian tissues and cells. Anal Chim Acta. 2012; 727: 8–12. pmid:22541816
  24. 24. Hussein SZ, Mohd Yusoff K, Makpol S, Mohd Yusof YA. Gelam Honey Inhibits the Production of Proinflammatory, Mediators NO, PGE(2), TNF-alpha, and IL-6 in Carrageenan-Induced Acute Paw Edema in Rats. Evid Based Complement Alternat Med. 2012; 2012: 109636. pmid:22919407
  25. 25. Bulau P, Zakrzewicz D, Kitowska K, Leiper J, Gunther A, Grimminger F, et al. Analysis of methylarginine metabolism in the cardiovascular system identifies the lung as a major source of ADMA. Am J Physiol Lung Cell Mol Physiol. 2007; 292: L18–24. pmid:16891395
  26. 26. Edoardo Gaude FC, Spiliotopoulos D, Spitaleri A, Ghitti M, Garcìa-Manteiga JM, Mari S, et al. An R Package for Metabolomics Univariate and MultivariateStatistical Analysis. Current Metabolomics. 2013; 1: 180–189.
  27. 27. Team RC. R: A language and environment for statistical computing. 2014.
  28. 28. Rohn H, Junker A, Hartmann A, Grafahrend-Belau E, Treutler H, Klapperstuck M, et al. VANTED v2: a framework for systems biology applications. BMC Syst Biol. 2012; 6: 139. pmid:23140568
  29. 29. Xia J, Mandal R, Sinelnikov IV, Broadhurst D, Wishart DS. MetaboAnalyst 2.0—a comprehensive server for metabolomic data analysis. Nucleic Acids Res. 2012; 40: W127–133. pmid:22553367
  30. 30. Xia J, Wishart DS. Web-based inference of biological patterns, functions and pathways from metabolomic data using MetaboAnalyst. Nat Protoc. 2011; 6: 743–760. pmid:21637195
  31. 31. Matsui T, Ishikawa T, Ito H, Okamoto M, Inoue K, Lee MC, et al. Brain glycogen supercompensation following exhaustive exercise. J Physiol. 2012; 590: 607–616. pmid:22063629
  32. 32. Mori M, Miura S, Tatibana M, Cohen PP. Cell-free translation of carbamyl phosphate synthetase I and ornithine transcarbamylase messenger RNAs of rat liver. Effect of dietary protein and fasting on translatable mRNA levels. J Biol Chem. 1981; 256: 4127–4132. pmid:7217073
  33. 33. Chen ZP, Mitchelhill KI, Michell BJ, Stapleton D, Rodriguez-Crespo I, Witters LA, et al. AMP-activated protein kinase phosphorylation of endothelial NO synthase. FEBS Lett. 1999; 443: 285–289. pmid:10025949
  34. 34. Sessa WC. eNOS at a glance. J Cell Sci. 2004; 117: 2427–2429. pmid:15159447
  35. 35. Lubec B, Hayn M, Kitzmuller E, Vierhapper H, Lubec G. L-Arginine reduces lipid peroxidation in patients with diabetes mellitus. Free Radic Biol Med. 1997; 22: 355–357. pmid:8958162
  36. 36. Kuhlencordt PJ, Chen J, Han F, Astern J, Huang PL. Genetic deficiency of inducible nitric oxide synthase reduces atherosclerosis and lowers plasma lipid peroxides in apolipoprotein E-knockout mice. Circulation. 2001; 103: 3099–3104. pmid:11425775
  37. 37. Tanaka M, Baba Y, Kataoka Y, Kinbara N, Sagesaka YM, Kakuda T, et al. Effects of (-)-epigallocatechin gallate in liver of an animal model of combined (physical and mental) fatigue. Nutrition. 2008; 24: 599–603. pmid:18455659
  38. 38. Berlett BS, Stadtman ER. Protein oxidation in aging, disease, and oxidative stress. J Biol Chem. 1997; 272: 20313–20316. pmid:9252331
  39. 39. Lentner C. Geigy Scientific Tables, 8th Rev edition. Ciba-Geigy Corp, Basel, Switzerland. 1981–1992.