Research Paper Volume 10, Issue 10 pp 2624—2635
A joint analysis of metabolomic profiles associated with muscle mass and strength in Caucasian women
- 1 Department of Preventive Medicine, College of Medicine, University of Tennessee Health Science Center, Memphis, TN 38163, USA
- 2 Tulane Center of Bioinformatics and Genomics, Department of Global Biostatistics and Data Science, Tulane University School of Public Health and Tropical Medicine, New Orleans, LA 70112, USA
- 3 Southeast Center for Integrated Metabolomics Core, University of Florida, Gainesville, FL 32610, USA
- 4 Medicinal Chemistry Core, Office of Research, University of Tennessee Health Science Center, Memphis, TN 38163, USA
- 5 Department of Pharmaceutical Sciences, College of Pharmacy, University of Tennessee Health Science Center, Memphis, TN 38163, USA
- 6 Proteomics and Metabolomics Core, Office of Research, University of Tennessee Health Science Center, Memphis, TN 38163, USA
- 7 School of Basic Medical Science, Central South University, Changsha, Hunan 410013, China
Received: July 17, 2018 Accepted: September 24, 2018 Published: October 14, 2018
https://doi.org/10.18632/aging.101574How to Cite
Abstract
Both loss of muscle mass and strength are important sarcopenia-related traits. In this study, we investigated both specific and shared serum metabolites associated with these two traits in 136 Caucasian women using a liquid chromatography-mass spectrometry method. A joint analysis of multivariate traits was used to examine the associations of individual metabolites with muscle mass measured by the body mass index-adjusted appendicular lean mass (ALM/BMI) and muscle strength measured by hand grip strength (HGS). After adjusting for multiple testing, nine metabolites including two amino acids (aspartic acid and glutamic acid) and an amino acid derive (pipecolic acid), one peptide (phenylalanyl-threonine), one carbohydrate (methyl beta-D-glucopyranoside), and four lipids (12S-HETRE, arachidonic acid, 12S-HETE, and glycerophosphocholine) were significant in the joint analysis. Of them, the two amino acids (aspartic acid and glutamic acid) and two lipids (12S-HETRE and 12S-HETE) were associated with both ALM/BMI and HGS, and the other five were only associated with ALM/BMI. The pathway analysis showed the amino acid metabolism pathways (aspartic acid and glutamic acid) might play important roles in the regulation of muscle mass and strength. In conclusion, our study identified novel metabolites associated with sarcopenia-related traits, suggesting novel metabolic pathways for muscle regulation.
Introduction
Age-related sarcopenia has become a major contributor to the risk of physical frailty, functional impairment, poor health-related quality of life, and premature death among the elderly [1–3]. It has emerged as a major health concern and resulted in a huge economic burden as the older population rapidly grows in the US [1]. In 2000, the direct medical costs attributable to sarcopenia was estimated to be about $18.5 billion in the US alone [1]. However, the etiology of sarcopenia remains largely unknown, and there are no approved medications for treating this disorder. Therefore, there is an urgent need for understanding the mechanisms of sarcopenia with the aim of improving its prevention, diagnosis, treatment, and prognosis. Although sarcopenia was originally considered to be solely loss of muscle mass, recently developed definitions of sarcopenia have incorporated both loss of muscle mass and muscle function [4–8]. It is supported by the facts that age-related decline in muscle strength does not parallel with the decline in muscle mass [9] and muscle weakness is a more consistent predictor for impaired physical performance, functional limitation, and physical disability compared to the loss in muscle mass [10]. This consensus highlights the necessity and importance of considering different physiological components of sarcopenia for obtaining a fundamental understanding of the complex biological mechanism underlying the decline in muscle mass and function.
Metabolomics is an emerging approach to systematically profile small molecules in biofluids, cells, and tissues [11]. Because metabolites represent the downstream expression of genome, transcriptome, and proteome, their study is hence most powerful to reveal inherent omics variation closest to the disease risk/phenotype [11]. Some metabolomics studies of muscle mass have been reported, recently [12–17]. Some amino acids (e.g., leucine, isoleucine, and glutamic acid) [12,16], lipids [14], and metabolites related to gut bacterial metabolism have been linked with muscle mass. However, very few metabolomic studies of muscle function have been conducted [13], and no studies have considered both muscle mass and muscle function at the same time. In fact, loss of muscle mass and strength starts at approximately 30 years of age. However, most of the previous metabolomic studies were conducted in the older population aged 60 years or older, which might fail to provide early metabolic markers for the loss of muscle mass and function.
The overall goal of this exploratory study was to identify novel metabolites and metabolic pathways associated with muscle mass and strength, providing metabolic markers for the early prediction of sarcopenia. We used an untargeted metabolomics approach and jointly analyzed the muscle mass and muscle strength to identify specific and shared metabolic factors for these two key sarcopenia traits in a sample of Caucasian women of a relatively young age.
Results
Characteristics of the sample of 136 women are shown in Table 1. The study participants were relatively young with an average age (SD) of 31.5 (5.1) years. The mean BMI, 26.0 kg/m2, was overweight. Over 35% of the participants were current smokers. On average, the participants had physical exercise about three times a week.
Table 1. Characteristics of the study subjects (N=136).
Variable | Mean (SD) or Percentage | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Age, year | 31.5 (5.1) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Weight, kg | 70.3 (21.4) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Height, cm | 164.6 (6.4) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
BMI, kg/m2 | 26.0 (7.5) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Current smoking, % | 35.3 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Alcohol drinking, gram/day | 36.5 (52.4) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Physical activity, times/week | 3.1 (2.2) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Dairy intake, cups/day | 1.6 (1.3) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ALM/BMI index | 0.78 (0.12) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Hand grip strength, kg | 27.0 (8.4) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ALM, appendicular lean mass; BMI, body mass index; SD, standard deviation. |
We identified a total of 206 metabolites with known identities and passed quality control using the LC-MS metabolomics platform. In the joint analysis, nine metabolites were significantly associated with muscle mass and muscle strength after adjusting for multiple testing (FDR q < 0.05) (Table 2). They belong to the chemical classes of organic acids, amino acids, peptides, carbohydrates, and lipids, respectively. All the nine metabolites were significantly associated with ALM/BMI. Among them, two amino acids (aspartate and glutamic acid) and two lipids [12(S)-HETRE and 12(S)-HETE] were also associated with HGS. The directions of the associations of these four metabolites with ALM/BMI and HGS were consistent (Figure 1). Another two lipids, arachidonic acid and glycerophosphocholine, were negatively associated with ALM/BMI as well. In contrast, the carbohydrate, methyl β-D-galactoside, was positively associated with ALM/BMI. Figure 2 shows the correlations among the significant metabolites. Three metabolites, 12(S)-HETRE, 12(S)-HETE, and arachidonic acid, were highly correlated each other with the correlation coefficients over 0.7. We didn’t observe any significant pairwise interactions among these significant metabolites in the relation to ALM/BMI or HGS after adjusting for multiple testing (Supplementary Table 1).
Table 2. The metabolites significantly associated with ALM/BMI or HGS.
Metabolite | Class | m/z | RT | P value | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Joint analysis | ALM/BMI | HGS | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Pipecolic acid | Organic acid | 130.0862 | 1.25 | 7.69×10-5 | 4.32×10-5 | 9.86×10-1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
12(S)-HETRE | Lipid | 321.2438 | 14.17 | 9.34×10-5 | 5.18×10-5 | 2.56×10-2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Aspartic acid | Amino Acid | 134.0446 | 0.73 | 2.25×10-4 | 5.33×10-3 | 4.29×10-4 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Phenylalanyl-threonine | Peptide | 267.1334 | 6.43 | 4.85×10-4 | 1.43×10-4 | 1.35×10-1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Methyl β-D-galactoside | Carbohydrate | 217.0679 | 1.02 | 5.32×10-4 | 1.60×10-4 | 5.67×10-1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Arachidonic acid | Lipid | 305.2467 | 14.8 | 7.86×10-4 | 2.87×10-4 | 8.30×10-2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
12(S)-HETE | Lipid | 319.2275 | 13.8 | 1.00×10-3 | 6.23×10-4 | 3.11×10-2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Glutamic acid | Amino Acid | 146.0458 | 0.76 | 1.02×10-3 | 7.63×10-4 | 2.32×10-2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Glycerophosphocholine | Lipid | 258.1094 | 0.76 | 1.71×10-3 | 4.91×10-4 | 1.81×10-1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ALM, appendicular lean mass; BMI, body mass index; HGS, hand grip strength; RT, retention time. |
Figure 1. Effect sizes of the significant metabolites on ALM/BMI and HGS. The effects were associated with one standard deviation increase in the relative abundance of a metabolite. The bars show the 95% confidence internal for the effect of each metabolite. The metabolites which were significantly associated with both traits are lighted in red. ALM, appendicular lean mass; BMI, body mass index; HGS, hand grip strength.
Figure 2. Pairwise correlation coefficients among the metabolites significantly associated with sarcopenia traits.
In the pathway analysis, eight metabolic pathways were significant in the enrichment analysis (P values < 0.05) (Table 3). Of them, seven pathways include glutamic acid and/or aspartic acid. The impact values of the two amino acids on these pathways ranged widely (from almost 0 to 0.44). In addition, we found only pipecolic acid and glutamic acid significantly associated with age after adjusting for multiple covariables (P values < 0.05) (Table 4).
Table 3. Pathway analysis results of trait-related metabolites.
Pathways | Matched Metabolites | P value a | Impact value b | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Alanine, aspartate and glutamate metabolism | Glutamic acid; Aspartic acid | 9.36×10-4 | 4.4×10-1 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Nitrogen metabolism | Glutamic acid; Aspartic acid | 2.48×10-3 | 6.7×10-4 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Histidine metabolism | Glutamic acid; Aspartic acid | 3.15×10-3 | 5.1×10-4 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Aminoacyl-tRNA biosynthesis | Glutamic acid; Aspartic acid | 9.01×10-3 | 1.1×10-2 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Arginine and proline metabolism | Glutamic acid; Aspartic acid | 9.49×10-3 | 3.6×10-2 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
D-Glutamine and D-glutamate metabolism | Glutamic acid | 2.27×10-2 | 1.1×10-1 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Cyanoamino acid metabolism | Aspartic acid | 3.28×10-2 | 0 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Ether lipid metabolism | Glycerophosphocholine | 4.69×10-2 | 0 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
a.P values from the pathway enrichment analysis. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
b. Impact values from the pathway topology analysis. |
Table 4. Correlations between trait-related metabolites and age.
Metabolite | Correlation coefficient | P value | Partial correlation coefficient a | P value | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Pipecolic acid | -0.28 | 1.24×10-3 | -0.21 | 1.77×10-2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
12(S)-HETRE | 0.13 | 1.23×10-1 | 0.09 | 3.33×10-1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Aspartic acid | -0.07 | 4.55×10-1 | -0.08 | 3.50×10-1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Phenylalanyl-threonine | 0.08 | 3.84×10-1 | 0.01 | 8.67×10-1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Methyl β-D-galactoside | -0.19 | 3.21×10-2 | -0.13 | 1.57×10-1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Arachidonic acid | 0.11 | 2.17×10-1 | 0.08 | 3.43×10-1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
12(S)-HETE | 0.14 | 1.10×10-1 | 0.10 | 2.38×10-1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Glutamic acid | -0.11 | 2.07×10-1 | -0.17 | 4.98×10-2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Glycerophosphocholine | -0.10 | 2.34×10-1 | -0.17 | 5.23×10-2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
a. The covariables including body mass index, smoking, alcohol drinking, physical activities, and diary intakes were adjusted. |
Discussion
The present study, for the first time, jointly investigated the metabolites significantly associated with two important sarcopenia traits, muscle mass and strength. We identified nine metabolites, mainly amino acids and lipids, associated with muscle mass, four of which were also associated with muscle strength. Most of the identified metabolites represented novel findings compared to previous reports, highlighting the importance of amino acids and lipids in the regulation of muscle mass and function.
Previous studies have shown that amino acids are potent stimulators of muscle protein synthesis in both the young and elderly [18–20]. Circulating (plasma or serum) amino acids, especially branched chain amino acids (such as leucine, isoleucine, and valine), have been associated with muscle mass in previous studies [12,17]. We identified two non-branched chain amino acids, aspartic acid and glutamic acid, were significantly associated with both muscle mass and muscle strength. The finding is consistent with prior evidence regarding the association of glutamic acid with lean mass identified in a sample of UK women [16]. However, this is the first time to report that glutamic acid was also associated with muscle strength in human studies. Glutamate (the anion of glutamic acid) in skeletal muscle participates in various metabolic pathways, such as glutathione synthesis, insulin production, tricarboxylic acid cycle, and purine nucleotide cycle. Muscle glutamate metabolism is disturbed in many conditions, such as hypoxia and oxidative stress [21]. Aspartate (the anion of aspartic acid) has been shown to inhibit inflammation-induced muscle loss through regulating phosphorylation of Akt, AMPKα, and FOXO1 in a piglet model [22]. Our pathway analysis further highlights potential important roles of several amino acid metabolic pathways involving both glutamic acid and aspartic acid, such as the pathways of alanine, aspartate, and glutamate metabolism as well as arginine and proline metabolism, in muscle regulation.
Pipecolic acid, the most significant metabolite associated with muscle mass in this study, is also an amino acid-related molecule. It is a metabolite of lysine which is an essential amino acid and very abundant in muscle tissue. It has been recently reported that oral administration of lysine could suppress proteolysis in skeletal muscles of fasted rats [23] and continuous feeding of a lysine-rich diet might increase muscle mass [24]. Supplementation of lysine in rats could reduce the upregulated autophagy activity of skeletal muscle caused by a low-protein diet and activate the Akt/mTOR signaling pathways [24]. As a metabolite of lysine, pipecolic acid also could stimulate the rates of protein synthesis, contributing to the effect of lysine on protein turn over in skeletal muscle [25]. We found that serum pipecolic acid was positively associated with muscle mass and muscle strength, but negatively associated with age. These findings may suggest pipecolic acid is a potential early marker for sarcopenia.
There is increasing evidence that supports a role of lipids and their derives in the regulation of skeletal muscle mass and function through their ability to modulate muscle cell growth, proliferation, and/or differentiation [26]. In particular, a recent study has observed significant changes in lipid mediators, a class of bioactive metabolites of the essential polyunsaturated fatty acids (PUFA), in skeletal muscles with aging in mice [27]. We identified one fatty acid (arachidonic acid) and two fatty acid derivatives [12(S)-HETE and 12(S)-HETRE] that were associated with both muscle mass and strength. Arachidonic acid, a PUFA, is the precursor of prostaglandins, leukotrienes, and related compounds, which have important roles in inflammation and in the regulation of immunity [28]. 12(S)-HETE and 12(S)-HETRE are derivatives of arachidonic acid and dihomo-λ-linolenic acid metabolism through the 12-lipoxygenase pathway, respectively [29]. We found the levels of the three molecules were highly correlated with each other in the serum and negatively associated with muscle mass and strength, suggesting they might be risk factors for the loss of muscle mass and function. However, previous studies showed that arachidonic acid could enhance in vitro skeletal muscle cell growth via a COX-2-dependent pathway [30]. Also, arachidonic acid supplementation could augment strength-training-induced adaptations in resistance-trained males [31]. These inconsistent findings might indicate the complexity of the role of arachidonic acid in skeletal muscle regulation. Their effects might vary by sex, age, and different pathological conditions. Therefore, it will be worth conducting metabolomics studies which target the whole pathway of arachidonic acid metabolism and even related pathways of fatty acid metabolism to investigate their comprehensive role in regulating muscle mass and function.
Arachidonic acid and its lipid mediators, including prostaglandin F2α (PGF2α), PGE2, and PGA2, decrease with aging in mice gastrocnemius muscles [27]. However, just like most of other trait-related metabolites identified in this study, we did not observe significant correlation between serum arachidonic acid and age. This might be because of differential changing of these metabolites in the muscle tissue and blood with aging. The relatively limited age range of the study population (from 20 to 40 years) might also undermine the study power to detect the change in the metabolites with aging. Longitudinal studies are necessary to clarify the changes in these metabolites with aging in the future. Furthermore, we identified another lipid, glycerophosphocholine, might be a risk factor for the loss of muscle mass. Very few studies have investigated the role of glycerophosphocholine in muscle regulation. Only a rat study showed an increased level of glycerophosphocholine in aged gastrocnemius [32].
We identified methyl β-D-galactoside, a member of carbohydrates, was positively associated with muscle mass. Methyl β-D-galactoside has been found in cereals and cereal products. Previous studies showed that carbohydrate consumption could stimulate muscle protein synthesis after resistance exercise [33]. On the contrary, other studies augured that the muscle protein increase during carbohydrate supplementation alone following resistance exercise is mainly through reduction at protein degradation rates rather than through increasing protein synthesis. Little is known about the mechanisms through which methyl β-D-galactoside regulates muscle mass. Further investigations are needed to clarify its role in muscle mass regulation.
Our study has several strengths. First, we considered both muscle mass and muscle strength which reflect important pathophysiological aspects of sarcopenia, in examining the associations with blood metabolites. Second, we used DXA to measure body compositions, which is widely used and well accepted to assess muscle mass [34]. Third, the non-hypothesis-driven untargeted metabolomics approach enabled the study to discover novel metabolites and pathways involved in muscle regulation. Fourth, adjustment for potential confounding factors and multiple testing were conducted to minimize false-positive findings. Finally, we identified both specific and shared metabolites for muscle mass and muscle strengths, providing important knowledge for better understanding the pathogenesis of sarcopenia. However, some limitations should be considered in interpreting the study findings. Specifically, the causal relationship between identified metabolites and sarcopenia traits could not be inferred because of the cross-sectional study design. Also, the relatively small sample size may limit the statistical power to identify sarcopenia-trait related metabolites with small effects. However, this pilot study was developed to generate hypotheses for specific/shared metabolites for muscle mass and strength. Further studies are needed to confirm our findings and clarify the potential molecular mechanisms of the blood metabolites associated with sarcopenia traits.
In summary, we identified trait-specific or shared metabolites associated with muscle mass and strength. Our results highlight the importance of multiple metabolic pathways of amino acids and lipids in the regulation of muscle mass and function. Future studies targeting these pathways may lead to a better understanding of the pathogenesis of sarcopenia and identification of biomarkers for early prediction and prevention of this disorder.
Materials and Methods
Study subjects
A total of 136 Caucasian women participants were from the ongoing Louisiana Osteoporosis Study (LOS) (from 2011 to current), which aims to build a large sample pool and database for investigating genetic and environmental factors for osteoporosis in Southern Louisiana. The inclusion and exclusion criteria of LOS have been described in our previous study [35]. All the study subjects included in this study were aged between 20 and 40 years old. Individuals who were pregnant, had a bilateral oophorectomy, or had any chronic conditions (such as diabetes mellitus, renal failure, liver failure, lung disease, gastrointestinal disease, and inherited bone disease) were excluded from the current study. The study was approved by the institutional review board, and a written consent form was signed by each participant before any data and bio-sample collection.
Clinical measurements
All participants completed an interviewer-assisted comprehensive questionnaire to collect demographic information, lifestyle (including smoking, drinking, and physical activity), dietary factors (including dairy consumption), reproductive and medical history [35]. Weight was measured in light indoor clothing using a calibrated balance beam scale, and height was measured using a calibrated stadiometer without shoes. Body mass index (BMI) was calculated as weight (kg) divided by height squared (m2).
Total body and regional measures of lean mass, fat mass, and bone mineral content were acquired using a dual-energy X-ray absorptiometry (DXA) machine (Hologic Inc., Bedford, MA) by trained and certified research staff. The machine was calibrated daily, and software and hardware were kept up-to-date during the data collection process. More details on data quality control including the usual covariation for repeated measures have been reported previously [35]. Appendicular lean mass (ALM) was calculated as the sum of lean mass in the arms and legs. The BMI-adjusted ALM (ALM/BMI) was used to assess individuals’ muscle mass in the study [8]. Hand grip strength (HGS) was measured using the Jamar 1 hand-held dynamometer (TEC Inc., Clifton, NJ). Two measurements of maximum strength were taken at both hands. The maximum grip strength value was used to assess an individual’s muscle strength.
Metabolomic analyses
The liquid chromatography-mass spectrometry (LC-MS) based metabolomics platform developed by Dr. Garrett’s lab in the Southeast Center for Integrated Metabolomics at University of Florida was used to perform the metabolomic analyses. The laboratory protocols have been previously described [36]. Briefly, frozen serum samples (-80 °C) were thawed at room temperature. Each serum sample (100 μL) was mixed with 20 μL internal standard mix followed by vortex mixing for 20 s. Next, 800 μL of acetonitrile: acetone:methanol (8:1:1, v: v: v) was added and centrifuged at 20,000×g for 10 min at <10°C to remove proteins. The supernatant (250 μL) was transferred to a new 1 mL Eppendorf tube and dried under a gentle stream of nitrogen (Organomation Associates, Berlin, MA, USA). The dried sample was reconstituted in 100 μL of 0.1% formic acid in water containing 4 injection standards (T-Boc amino acids) and placed in an ice bath for 10-15 min followed by centrifugation at 20,000×g for 5 min at <10°C to remove any debris.
A Thermo Q-Exactive High-Resolution Mass Spectrometer (Thermo Fisher Scientific, Fremont, CA) coupled with a Dionex UHPLC (Dionex Corporation, Sunnyvale, CA) was used to conduct the metabolomic analysis. All the samples were analyzed in both positive and negative ion modes with heated electrospray ionization. The mass resolution was 35,000 at m/z 200 with a mass accuracy of less than 5 ppm in positive mode and less than 10 ppm in negative mode. Separation was achieved on an ACE C18-PFP column (100x2.1mm, 2µm) with 0.1% formic acid in water as mobile phase A and acetonitrile as mobile phase B with a column temperature of 25 °C. The flow rate was 350 µL/min with a total run time of 20 min.
Alignment and feature finding was performed using the open source software MZmine [37]. Metabolite identification was performed by searching an internal retention time library of over 600 compounds. Peaks in the MS were quantified using integrated peak height. Raw area counts for each metabolite in every sample was first normalized to the sum of all injection internal standards to correct for subtle injection differences, and then normalized to the total signal of each sample. The batch correction was performed using well-established Bayes method for microarray data by continually running the neat and pooled quality control (QC) samples. Metabolites with missing rates > 20% or coefficients of variation > 20% were excluded from further analyses. Imputation of missing data was performed using the R package ‘missForest’ [38]. The relative abundance data of metabolites were further log transformed and autoscaled to have zero mean and unit variance (z scores) using the R package ‘specmine’.
Statistical analyses
To increase the study power and decrease the multiple testing burden [39], we conducted a joint analysis of multivariate phenotypes to examine the associations of metabolites with sarcopenia traits (ALM/BMI and HGS). The null hypothesis (H0) of the joint analysis was that none of the traits was associated with the tested metabolite. At least one trait associated with the metabolite would reject the null hypothesis. The R package ‘systemfit’ was used to conduct the joint analysis [40]. The seemingly unrelated regression of ‘systemfit’ was conducted to estimate the coefficients of ALM/BMI and HGS with the metabolite, respectively, and the F-statistics were used to test the null hypothesis. Potential confounding factors, including age, BMI (for HGS only), smoking, alcohol drinking, physical activities, and diary intakes, were adjusted when examining the associations of ALM/BMI and HGS with metabolites. Considering the substantial correlations among metabolites, we used the false discovery rate (FDR) method to adjust for multiple testing. A metabolite with an FDR q-value < 0.05 was considered significant [41]. To evaluate the independence of significant metabolic associations, we estimated the pairwise Pearson correlations between all metabolites significantly associated with these traits. In addition, we examined the pairwise interaction effects of significant metabolites on ALM/BMI and HGS. To identify potential pathways represented by the significant metabolites, we conducted a pathway analysis which integrated enrichment analysis and pathway topology analysis using the web tool MetaboAnalyst [42]. The topology analysis takes into account the position of a metabolite in the pathway structure and provides an impact value, ranging from 0 (minimum impact) to 1 (maximum impact), for each tested metabolite. The correlation coefficients between significant metabolites and age without or with adjusting for covariables (including BMI, smoking, alcohol drinking, physical activities, and diary intakes) were also calculated to examine whether these metabolites changed with aging.
Supplementary Materials
Supplementary Table 1
Abbreviations
ALM/BMI: appendicular lean mass/body mass index; DXA: dual-energy X-ray absorptiometry; FDR: false discovery rate; HGS: hand grip strength; LC-MS: liquid chromatography-mass spectrometry; LOS: Louisiana Osteoporosis Study; SD: standard deviation; 12(S)-HETE: 12-hydroxyeicosatetraenoic acid; 12(S)-HETRE: (12S)-hydroxy-(8Z,10E,14Z)-eicosatrienoic acid.
Author Contributions
QZ, HS, and HWD designed the study; HS, LJZ, and CQ conducted data collection; TJG conducted the metabolomics analysis; QZ, HS, KJS, and QT analyzed the data; QZ, HS, and HWD drafted the manuscript; JL and DK provided critical review and contributed to the interpretation of the study findings.
Acknowledgments
We thank the LOS participants for their effort and time to participate in this research.
Conflicts of Interest
The authors have no conflicts of interest.
Funding
The investigators of this work were partially supported by grants from the National Institute of Health (NIH) (R01 AR069055, U19 AG055373, R01 MH104680, R01AR059781 and P20GM109036), and the Edward G. Schlieder Endowment as well as the Drs. W. C. Tsai and P. T. Kung Professorship in Biostatistics from Tulane University. Dr. Qi Zhao was also partially supported by the Scientist Development Grant from the American Heart Association (13SDG17020012). The Southeast Center for Integrated Metabolomics at University of Florida is supported by the NIH grant U24 DK097209.
References
- 1. Janssen I, Shepard DS, Katzmarzyk PT, Roubenoff R. The healthcare costs of sarcopenia in the United States. J Am Geriatr Soc. 2004; 52:80–85. https://doi.org/10.1111/j.1532-5415.2004.52014.x [PubMed]
- 2. Cruz-Jentoft AJ, Landi F, Schneider SM, Zúñiga C, Arai H, Boirie Y, Chen LK, Fielding RA, Martin FC, Michel JP, Sieber C, Stout JR, Studenski SA, et al. Prevalence of and interventions for sarcopenia in ageing adults: a systematic review. Report of the International Sarcopenia Initiative (EWGSOP and IWGS). Age Ageing. 2014; 43:748–59. https://doi.org/10.1093/ageing/afu115 [PubMed]
- 3. Beaudart C, Zaaria M, Pasleau F, Reginster JY, Bruyère O. Health outcomes of sarcopenia: a systematic review and meta-analysis. PLoS One. 2017; 12:e0169548. https://doi.org/10.1371/journal.pone.0169548 [PubMed]
- 4. Cruz-Jentoft AJ, Baeyens JP, Bauer JM, Boirie Y, Cederholm T, Landi F, Martin FC, Michel JP, Rolland Y, Schneider SM, Topinková E, Vandewoude M, Zamboni M, and European Working Group on Sarcopenia in Older People. Sarcopenia: European consensus on definition and diagnosis: Report of the European Working Group on Sarcopenia in Older People. Age Ageing. 2010; 39:412–23. https://doi.org/10.1093/ageing/afq034 [PubMed]
- 5. Fielding RA, Vellas B, Evans WJ, Bhasin S, Morley JE, Newman AB, Abellan van Kan G, Andrieu S, Bauer J, Breuille D, Cederholm T, Chandler J, De Meynard C, et al, and International Working Group on Sarcopenia. Sarcopenia: an undiagnosed condition in older adults. Current consensus definition: prevalence, etiology, and consequences. International working group on sarcopenia. J Am Med Dir Assoc. 2011; 12:249–56. https://doi.org/10.1016/j.jamda.2011.01.003 [PubMed]
- 6. Morley JE, Abbatecola AM, Argiles JM, Baracos V, Bauer J, Bhasin S, Cederholm T, Coats AJ, Cummings SR, Evans WJ, Fearon K, Ferrucci L, Fielding RA, et al, and Society on Sarcopenia, Cachexia and Wasting Disorders Trialist Workshop. Sarcopenia with limited mobility: an international consensus. J Am Med Dir Assoc. 2011; 12:403–09. https://doi.org/10.1016/j.jamda.2011.04.014 [PubMed]
- 7. Chen LK, Liu LK, Woo J, Assantachai P, Auyeung TW, Bahyah KS, Chou MY, Chen LY, Hsu PS, Krairit O, Lee JS, Lee WJ, Lee Y, et al. Sarcopenia in Asia: consensus report of the Asian Working Group for Sarcopenia. J Am Med Dir Assoc. 2014; 15:95–101. https://doi.org/10.1016/j.jamda.2013.11.025 [PubMed]
- 8. Studenski SA, Peters KW, Alley DE, Cawthon PM, McLean RR, Harris TB, Ferrucci L, Guralnik JM, Fragala MS, Kenny AM, Kiel DP, Kritchevsky SB, Shardell MD, et al. The FNIH sarcopenia project: rationale, study description, conference recommendations, and final estimates. J Gerontol A Biol Sci Med Sci. 2014; 69:547–58. https://doi.org/10.1093/gerona/glu010 [PubMed]
- 9. Delmonico MJ, Harris TB, Visser M, Park SW, Conroy MB, Velasquez-Mieyer P, Boudreau R, Manini TM, Nevitt M, Newman AB, Goodpaster BH, and Health, Aging, and Body. Longitudinal study of muscle strength, quality, and adipose tissue infiltration. Am J Clin Nutr. 2009; 90:1579–85. https://doi.org/10.3945/ajcn.2009.28047 [PubMed]
- 10. Manini TM, Clark BC. Dynapenia and aging: an update. J Gerontol A Biol Sci Med Sci. 2012; 67:28–40. https://doi.org/10.1093/gerona/glr010 [PubMed]
- 11. Johnson CH, Ivanisevic J, Siuzdak G. Metabolomics: beyond biomarkers and towards mechanisms. Nat Rev Mol Cell Biol. 2016; 17:451–59. https://doi.org/10.1038/nrm.2016.25 [PubMed]
- 12. Lustgarten MS, Price LL, Chale A, Phillips EM, Fielding RA. Branched chain amino acids are associated with muscle mass in functionally limited older adults. J Gerontol A Biol Sci Med Sci. 2014; 69:717–24. https://doi.org/10.1093/gerona/glt152 [PubMed]
- 13. Lustgarten MS, Price LL, Chalé A, Fielding RA. Metabolites related to gut bacterial metabolism, peroxisome proliferator-activated receptor-alpha activation, and insulin sensitivity are associated with physical function in functionally-limited older adults. Aging Cell. 2014; 13:918–25. https://doi.org/10.1111/acel.12251 [PubMed]
- 14. Moaddel R, Fabbri E, Khadeer MA, Carlson OD, Gonzalez-Freire M, Zhang P, Semba RD, Ferrucci L. Plasma biomarkers of poor muscle quality in older men and women from the Baltimore Longitudinal Study of Aging. J Gerontol A Biol Sci Med Sci. 2016; 71:1266–72. https://doi.org/10.1093/gerona/glw046 [PubMed]
- 15. Lustgarten MS, Price LL, Phillips EM, Kirn DR, Mills J, Fielding RA. Serum predictors of percent lean mass in young adults. J Strength Cond Res. 2016; 30:2194–201. https://doi.org/10.1519/JSC.0b013e31829eef24 [PubMed]
- 16. Korostishevsky M, Steves CJ, Malkin I, Spector T, Williams FM, Livshits G. Genomics and metabolomics of muscular mass in a community-based sample of UK females. Eur J Hum Genet. 2016; 24:277–83. https://doi.org/10.1038/ejhg.2015.85 [PubMed]
- 17. Murphy RA, Moore SC, Playdon M, Meirelles O, Newman AB, Milijkovic I, Kritchevsky SB, Schwartz A, Goodpaster BH, Sampson J, Cawthon P, Simonsick EM, Gerszten RE, et al, and Health ABC Study. Metabolites associated with lean mass and adiposity in older black men. J Gerontol A Biol Sci Med Sci. 2017; 72:1352–59. [PubMed]
- 18. Bohé J, Low A, Wolfe RR, Rennie MJ. Human muscle protein synthesis is modulated by extracellular, not intramuscular amino acid availability: a dose-response study. J Physiol. 2003; 552:315–24. https://doi.org/10.1113/jphysiol.2003.050674 [PubMed]
- 19. Volpi E, Mittendorfer B, Wolf SE, Wolfe RR. Oral amino acids stimulate muscle protein anabolism in the elderly despite higher first-pass splanchnic extraction. Am J Physiol. 1999; 277:E513–20. [PubMed]
- 20. Børsheim E, Bui QU, Tissier S, Kobayashi H, Ferrando AA, Wolfe RR. Effect of amino acid supplementation on muscle mass, strength and physical function in elderly. Clin Nutr. 2008; 27:189–95. https://doi.org/10.1016/j.clnu.2008.01.001 [PubMed]
- 21. Rutten EP, Engelen MP, Schols AM, Deutz NE. Skeletal muscle glutamate metabolism in health and disease: state of the art. Curr Opin Clin Nutr Metab Care. 2005; 8:41–51. https://doi.org/10.1097/00075197-200501000-00007 [PubMed]
- 22. Liu Y, Wang X, Leng W, Pi D, Tu Z, Zhu H, Shi H, Li S, Hou Y, Hu CA. Aspartate inhibits LPS-induced MAFbx and MuRF1 expression in skeletal muscle in weaned pigs by regulating Akt, AMPKα and FOXO1. Innate Immun. 2017; 23:34–43. https://doi.org/10.1177/1753425916673443 [PubMed]
- 23. Sato T, Ito Y, Nagasawa T. Regulation of skeletal muscle protein degradation and synthesis by oral administration of lysine in rats. J Nutr Sci Vitaminol (Tokyo). 2013; 59:412–19. https://doi.org/10.3177/jnsv.59.412 [PubMed]
- 24. Sato T, Ito Y, Nagasawa T. Dietary L-Lysine suppresses autophagic proteolysis and stimulates Akt/mTOR signaling in the skeletal muscle of rats fed a low-protein diet. J Agric Food Chem. 2015; 63:8192–98. https://doi.org/10.1021/acs.jafc.5b03811 [PubMed]
- 25. Sato T, Ito Y, Nagasawa T. Regulatory effects of the L-lysine metabolites, L-2-aminoadipic acid and L-pipecolic acid, on protein turnover in C2C12 myotubes. Biosci Biotechnol Biochem. 2016; 80:2168–75. https://doi.org/10.1080/09168451.2016.1210499 [PubMed]
- 26. Lipina C, Hundal HS. Lipid modulation of skeletal muscle mass and function. J Cachexia Sarcopenia Muscle. 2017; 8:190–201. https://doi.org/10.1002/jcsm.12144 [PubMed]
- 27. Wang Z, Bian L, Mo C, Kukula M, Schug KA, Brotto M. Targeted quantification of lipid mediators in skeletal muscles using restricted access media-based trap-and-elute liquid chromatography-mass spectrometry. Anal Chim Acta. 2017; 984:151–61. https://doi.org/10.1016/j.aca.2017.07.024 [PubMed]
- 28. Calder PC. Polyunsaturated fatty acids, inflammation, and immunity. Lipids. 2001; 36:1007–24. https://doi.org/10.1007/s11745-001-0812-7 [PubMed]
- 29. Ikei KN, Yeung J, Apopa PL, Ceja J, Vesci J, Holman TR, Holinstat M. Investigations of human platelet-type 12-lipoxygenase: role of lipoxygenase products in platelet activation. J Lipid Res. 2012; 53:2546–59. https://doi.org/10.1194/jlr.M026385 [PubMed]
- 30. Markworth JF, Cameron-Smith D. Arachidonic acid supplementation enhances in vitro skeletal muscle cell growth via a COX-2-dependent pathway. Am J Physiol Cell Physiol. 2013; 304:C56–67. https://doi.org/10.1152/ajpcell.00038.2012 [PubMed]
- 31. De Souza EO, Lowery RP, Wilson JM, Sharp MH, Mobley CB, Fox CD, Lopez HL, Shields KA, Rauch JT, Healy JC, Thompson RM, Ormes JA, Joy JM, Roberts MD. Effects of Arachidonic acid supplementation on acute anabolic signaling and chronic functional performance and body composition adaptations. PLoS One. 2016; 11:e0155153. https://doi.org/10.1371/journal.pone.0155153 [PubMed]
- 32. Garvey SM, Dugle JE, Kennedy AD, McDunn JE, Kline W, Guo L, Guttridge DC, Pereira SL, Edens NK. Metabolomic profiling reveals severe skeletal muscle group-specific perturbations of metabolism in aged FBN rats. Biogerontology. 2014; 15:217–32. https://doi.org/10.1007/s10522-014-9492-5 [PubMed]
- 33. Børsheim E, Cree MG, Tipton KD, Elliott TA, Aarsland A, Wolfe RR. Effect of carbohydrate intake on net muscle protein synthesis during recovery from resistance exercise. J Appl Physiol (1985). 2004; 96:674–78. https://doi.org/10.1152/japplphysiol.00333.2003 [PubMed]
- 34. Edwards MH, Buehring B. Novel approaches to the diagnosis of sarcopenia. J Clin Densitom. 2015; 18:472–77. https://doi.org/10.1016/j.jocd.2015.04.010 [PubMed]
- 35. He H, Liu Y, Tian Q, Papasian CJ, Hu T, Deng HW. Relationship of sarcopenia and body composition with osteoporosis. Osteoporos Int. 2016; 27:473–82. https://doi.org/10.1007/s00198-015-3241-8 [PubMed]
- 36. Liu H, Garrett TJ, Su Z, Khoo C, Gu L. UHPLC-Q-Orbitrap-HRMS-based global metabolomics reveal metabolome modifications in plasma of young women after cranberry juice consumption. J Nutr Biochem. 2017; 45:67–76. https://doi.org/10.1016/j.jnutbio.2017.03.007 [PubMed]
- 37. Pluskal T, Castillo S, Villar-Briones A, Oresic M. MZmine 2: modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinformatics. 2010; 11:395. https://doi.org/10.1186/1471-2105-11-395 [PubMed]
- 38. Stekhoven DJ, Bühlmann P. MissForest--non-parametric missing value imputation for mixed-type data. Bioinformatics. 2012; 28:112–18. https://doi.org/10.1093/bioinformatics/btr597 [PubMed]
- 39. Suo C, Toulopoulou T, Bramon E, Walshe M, Picchioni M, Murray R, Ott J. Analysis of multiple phenotypes in genome-wide genetic mapping studies. BMC Bioinformatics. 2013; 14:151. https://doi.org/10.1186/1471-2105-14-151 [PubMed]
- 40. Henningsen A, Hamann JD. systemfit: A package for estimating systems of simultaneous equations in R. J Stat Softw. 2007; 23:1–40. https://doi.org/10.18637/jss.v023.i04
- 41. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. JR Stat Soc. 1995; 57:289–300.
- 42. Xia J, Wishart DS. Using MetaboAnalyst 3.0 for Comprehensive Metabolomics Data Analysis. Curr Protoc Bioinformatics. 2016; 55:14.10.1–14.10.91. https://doi.org/10.1002/cpbi.11 [PubMed]