Next Article in Journal
In Vivo Estimation of Ketogenesis Using Metabolic Flux Analysis—Technical Aspects and Model Interpretation
Next Article in Special Issue
Metabolomics Studies in Psoriatic Disease: A Review
Previous Article in Journal
Application of Lipid Class Ratios for Sample Stability Monitoring—Evaluation of Murine Tissue Homogenates and SDS as a Stabilizer
Previous Article in Special Issue
Alterations of Extracellular Matrix Components in the Course of Juvenile Idiopathic Arthritis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 11
Issue 5
10.3390/metabo11050278
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Metabolomic Profile Predictive of New Osteoporosis or Sarcopenia Development

by
Kana Miyamoto
1,
Akiyoshi Hirayama
2,
Yuiko Sato
3,4,5,
Satsuki Ikeda
2,
Midori Maruyama
2,
Tomoyoshi Soga
2,
Masaru Tomita
2,
Masaya Nakamura
3,
Morio Matsumoto
3,
Noriko Yoshimura
6 and
Takeshi Miyamoto
1,3,4,5,*
1
Department of Orthopedic Surgery, Kumamoto University, 1-1-1 Honjo, Chuo-ku, Kumamoto 860-8556, Japan
2
Institute for Advanced Biosciences, Keio University, 246-2 Mizukami, Kakuganji, Tsuruoka, Yamagata 997-0052, Japan
3
Department of Orthopedic Surgery, School of Medicine, Keio University, 35 Shinano-machi, Shinjuku-ku, Tokyo 160-8582, Japan
4
Department of Advanced Therapy for Musculoskeletal Disorders II, Keio University School of Medicine, 35 Shinano-machi, Shinjuku-ku, Tokyo 160-8582, Japan
5
Department of Musculoskeletal Reconstruction and Regeneration Surgery, Keio University School of Medicine, 35 Shinano-machi, Shinjuku-ku, Tokyo 160-8582, Japan
6
Department of Preventive Medicine for Locomotive Organ Disorders, 22nd Century Medical and Research Center, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-8655, Japan
*
Author to whom correspondence should be addressed.
Metabolites 2021, 11(5), 278; https://doi.org/10.3390/metabo11050278
Submission received: 4 April 2021 / Revised: 26 April 2021 / Accepted: 26 April 2021 / Published: 28 April 2021
(This article belongs to the Special Issue Biomarkers and Human Blood Metabolites)
Review Reports Versions Notes

Abstract

:
The increasing number of patients with osteoporosis and sarcopenia is a global concern among countries with progressively aging societies. The high medical costs of treating those patients suggest that prevention rather than treatment is preferable. We enrolled 729 subjects who attended both the second and third surveys of the Research on Osteoarthritis/Osteoporosis Against Disability (ROAD) study. Blood samples were collected from subjects at the second survey, and then a comprehensive metabolomic analysis was performed. It was found that 35 had newly developed osteoporosis at the third survey performed four years later, and 39 were newly diagnosed with sarcopenia at the third survey. In the second survey, we found that serum Gly levels were significantly higher even after adjustment for age, sex, and BMI in subjects with newly developed osteoporosis relative to those who remained osteoporosis-negative during the four-year follow-up. We also show that serum taurine levels were significantly lower at the second survey, even after adjustment for age, sex, and BMI in subjects with newly developed sarcopenia during the four-year follow-up compared with those not diagnosed with sarcopenia at the second or third surveys. Though our sample size and odds ratios were small, increased Gly and decreased taurine levels were found to be predictive of new development of osteoporosis and sarcopenia, respectively, within four years.

1. Introduction

Osteoporosis and sarcopenia are the most important musculoskeletal diseases of the elderly [1]. Both develop with aging and are risks for loss of activity of daily living (ADL) and quality of life (QOL), as well as mortality, in the elderly [2,3]. As the both the number and rate in the elderly increase in various populations, preventing these diseases before they develop is mandatory. To do this, factors predictive of future disease are needed. However, both osteoporosis and sarcopenia are multifactorial and not attributable to a single factor, making prediction difficult.
Osteoporosis is characterized by reduced bone strength and increased risk for fragility fractures, and morbidity rates reportedly increase with age [4]. Known risks for osteoporosis include aging, menopause, smoking, excessive alcohol consumption, steroid use, genetic factors, and inflammatory or metabolic disease. Metabolites associated with low bone mass have been identified [5,6], but no longitudinal predictors have been reported. The fracture risks assessment tool (FRAX®) can predict fragility fractures based on responses to questionnaires relevant to age, sex, weight, height, previous fracture, parent fractured hip, current smoking and alcohol habits, glucocorticoid use, rheumatoid arthritis, and secondary osteoporosis; however, no biological markers are used in this assessment [7].
Sarcopenia is defined as a progressive loss of muscle mass with reduced muscle strength and/or function with age [8,9]. Known risks for sarcopenia include aging, immobilization, and cachexia [10]. The European Working Group on Sarcopenia in Older People (EWGSOP) developed a practical clinical definition and diagnostic criteria for sarcopenia in 2010 [8]. After that, the Asian Working Group for Sarcopenia (AWGS) provided appropriate diagnostic cutoff values for Asian populations [9]. Osteoporosis and sarcopenia frequently simultaneously develop in the elderly [11], suggesting that common predictors in addition to aging may exist for these diseases.
The present study was performed using Osteoarthritis/Osteoporosis Against Disability (ROAD) cohorts established in 2005. The ROAD study is a national, prospective study of musculoskeletal diseases that consists of population-based cohorts from several communities in Japan. Details of the cohort profiles have been reported elsewhere [4,12]. Briefly, between 2005 and 2007, a baseline database including the clinical and genetic information of 3040 residents (1061 men and 1979 women with a mean age of 70.3 years (standard deviation (SD): 11.0); 71.0 (10.7) years for men and 69.9 (11.2) years for women) was created. The study revealed that sarcopenia prevalence increases by 60 years of age and that four years after sarcopenia development, subjects become at risk for osteoporosis [11].
Metabolites are the end products of various metabolic pathways, and their levels overall reflect individuals’ metabolomic status [13,14]. Metabolomic analysis can detect changes or differences in metabolite levels and is thus useful to diagnose diseases such as cancer and to detect low bone mass or bone status in human and animals [15,16,17,18]. Since osteoporosis and sarcopenia are metabolic diseases that disrupt bone and muscle homeostasis, metabolomic analyses could predict diseases before development.
Here, we undertook the metabolomic analysis of serum samples collected at the second follow-up of ROAD and found that levels of the amino acid glycine (Gly) were significantly higher in subjects not identified as having osteoporosis at that second survey but who had newly developed osteoporosis during the four-year follow-up, compared with osteoporosis-negative subjects at the second and third surveys, after adjustment for age, sex, and BMI. Similarly, levels of taurine were significantly lower in subjects newly defined as exhibiting sarcopenia at the third survey after adjustment for age, sex, and BMI. Thus, Gly and taurine could serve as predictors of the future development of osteoporosis and sarcopenia, respectively.

2. Results

2.1. High Serum Gly Levels Are Significantly Associated with New Osteoporosis Development

Among study participants, 460 did not meet the osteoporosis criteria at the second survey and were thus judged osteoporosis-negative, but 35 (7 males and 28 females) of those 460 were newly defined as positive for osteoporosis at the third survey. We compared patient characteristics at the second survey between subjects evaluated at the third survey as new osteoporosis (new OP) or non-osteoporosis (non-OP) (Table 1). That analysis further indicated that the proportion of subjects either of relatively advanced age or female gender was significantly higher in new- versus non-OP groups (Table 1). In contrast, values indicative of BMD, grip power, and walking speed were significantly lower in new- versus non-OP groups (Table 1). Grip power and walking speed are criteria used to diagnose sarcopenia, and the percentage of subjects exhibiting sarcopenia newly present at the third survey was significantly higher in new OP than non-OP cases (Table 1).
We then analyzed levels of IGF1 and various metabolites in sera collected at the second survey and compared those levels between new OP and non-OP cases (Table 2). Levels of some metabolites (shown in Table 2) differed significantly between new OP and non-OP subjects.
Age is one of the most common risk factors for osteoporosis development, and age indeed differed significantly between new OP and non-OP cases (Table 1). Thus, the parameters for which levels were significantly different between new OP and non-OP cases were adjusted for age (Table 3). Moreover, female gender and BMI were considered potential risk factors, and indeed, percentages associated with female gender and BMI values differed significantly between new OP and non-OP cases (Table 1). Thus, after adjusting parameters significantly associated with new osteoporosis development for age, we further adjusted them for age, sex, and BMI (Table 3). After adjustment for age, sex, and BMI, we found that relatively high Gly levels remained significantly associated with new OP cases (odds ratio: 1.003; 95% confidence interval: 1.000–1.006; p = 0.021) (Table 3).

2.2. Low Serum Taurine Levels Are Significantly Associated with New Sarcopenia Development

Among participants, 492 did not satisfy the sarcopenia criteria at the second survey, while at the third survey, 39 (14 males and 25 females) of that 492 did. We analyzed patient characteristics at the second survey and compared them between new sarcopenia (new SP) and non-sarcopenia (non-SP) subjects, (Table 4). Age was significantly higher and BMI significantly lower in the new SP group compared to non-SP group (Table 4). Moreover, BMD, grip power, and walking speed were all significantly lower in the new SP group relative to the non-SP group (Table 4). Osteoporosis prevalence at the second and third surveys was significantly higher in new SP subjects (Table 4), suggesting an association with sarcopenia.
Subsequently, we analyzed levels of IGF1 and various metabolites in sera collected at the second survey. Metabolites whose levels significantly differed between new SP and non-SP groups are shown in Table 5. Serum IGF1 levels were significantly lower in the new SP group than the non-SP group (Table 5). Valine, leucine, and isoleucine, all branched-chain amino acids known to be abundant in skeletal muscle, were significantly lower in new SP cases compared to non-SP cases. Similarly, levels of taurine, also a skeletal muscle-related amino acid, were significantly lower in new SP cases compared to non-SP cases (Table 5).
We adjusted metabolites whose levels differed significantly between the new SP and non-SP groups by age (Table 6). We then adjusted factors significantly associated with new SP after adjustment with by age, sex, and BMI (Table 6). A logistic regression analysis after adjustment revealed that relatively reduced taurine levels (odds ratio: 0.991; 95% confidence interval: 0.982–0.999; p = 0.039) and existing osteoporosis (odds ratio: 3.118; 95% confidence interval: 1.215–7.997; p = 0.018) at the second survey predicted new sarcopenia development within the next four years (Table 6).

3. Discussion

Osteoporosis and sarcopenia are age-related diseases that are increasing in the number of cases and prevalence in developed countries. Both worsen ADL and QOL in the elderly. These conditions may be correlated and influence each other. In support of this idea, we found that among patients with newly developed osteoporosis, 15.2% also had newly developed sarcopenia. In contrast, 5.8% of subjects enrolled in the study who did not develop osteoporosis newly developed sarcopenia (Table 1). By contrast, 17.2% of patients with newly developed sarcopenia also newly developed osteoporosis, while 6.7% of subjects who did not newly develop sarcopenia newly developed osteoporosis (Table 4). Moreover, we found that elevated serum Gly levels predicted new osteoporosis development within four years but not the development of sarcopenia, while decreased serum taurine levels were associated with the new development of sarcopenia but not osteoporosis in a four-year period. Gly levels were higher in subjects with newly developed sarcopenia compared with those who did not, but differences were not statistically significant (non-SP vs. new SP: 494.3 ± 118.0 vs. 499.8 ± 138.6, respectively; p > 0.05). Moreover, taurine levels were significantly higher in subjects that newly developed osteoporosis relative to those who did not (non-OP vs. new OP: 252.0 ± 80.0 vs. 283.1 ± 126.1, respectively; p = 0.036). However, that difference also was not statistically significant after adjustment for age and BMI. Thus, although osteoporosis and sarcopenia are closely related diseases, methods to predict them likely differ overall.
Five subjects developed both osteoporosis and sarcopenia (the “both” group), but their Gly levels were comparable to those seen in individuals who developed osteoporosis but not sarcopenia (new OP alone) (both vs. new OP alone: 484.3 ± 93.7 vs. 561.0 ± 145.4, respectively; p > 0.05). Likewise, taurine levels were similar between the “both” group and subjects who developed sarcopenia but not osteoporosis (new SP alone) (both vs. new SP alone: 242.5 ± 49.7 vs. 221.3 ± 47.0, respectively; p > 0.05). In contrast, levels of the amino acid citrulline were significantly higher in the “both” group relative to new SP alone (57.8 ± 22.7 vs. 41.2 ± 10.3, respectively; p < 0.05). Citrulline supplementation reportedly promotes increased skeletal muscle mass [19], but at present, the mechanisms underlying increased citrulline levels in subjects who develop both osteoporosis and sarcopenia remain unclear.
Bone and muscle interact [20,21,22,23], and myokines produced from muscle regulate bone homeostasis [24]. For example, irisin, a myokine secreted from muscle with exercise, controls bone metabolism [25,26]. Extracellular matrix vesicles produced in and secreted by muscles reportedly regulate bone homeostasis [27]. Insulin-like growth factor 1 (IGF1), which is produced by osteoblasts and stored in extracellular matrix proteins in bones, is a bone-remodeling factor [28] that is released upon osteoclastic bone resorption and stimulates subsequent osteoblastic bone-forming activities. IGF1 also stimulates anabolic pathways and inhibits catabolic pathways in muscles [29]. Previously, we demonstrated that reduced IGF1 levels promoted decreased muscle mass in adult mice [30]. Indeed, here we showed lower serum IGF1 levels in the new- versus non-SP cases (Table 5). Transforming growth factor beta 1 (TGFβ1) is also a bone-remodeling factor released from the bone matrix [31]. TGFβ1 is activated by osteoclastic activity and stimulates bone formation by osteoblasts [31]. On the other hand, TGFβ1 activity promotes muscle atrophy [32]. Thus, the activities of bone and muscle are tightly linked, and the disruption of that interaction could promote the development of both osteoporosis and sarcopenia. Indeed, some subjects assessed here developed osteoporosis and sarcopenia concomitantly, although we showed that factors predictive of these conditions differ.
Several studies combining cohort and metabolomic analysis have been reported. For example, gut microbiome α-diversity is reportedly explained by a subset of 40 plasma metabolites in two independent human cohort blood samples [33]. Metabolite profiles are also reported to predict future diabetes [34]. Low plasma lysophosphatidylcholine (LPC) 18:2 is reportedly a predictor of decline in gait speed in older adults [35]. Here, although small sample size was a limitation of our study, we combined comprehensive metabolome analysis with a longitudinal cohort study to show for the first time that Gly and taurine levels in sera are predictors of osteoporosis and sarcopenia, respectively.
Glycine (Gly) is a non-essential amino acid, and increased glycine levels are reportedly seen in male idiopathic osteoporosis patients [36]. Moreover, in humans, supplementation with six amino acids including glycine is associated with a higher BMD in the spine and forearm [37]. We found that Gly levels were significantly and negatively correlated with bone mineral density (BMD) at either the lumbar spine L2–4 levels or the femoral neck (r = −0.1818 (p < 0.05) or −0.2134 (p < 0.05), respectively) at the second survey, supporting the idea that increased Gly levels are associated with decreased BMD. Gly is the most abundant amino acid residue in collagen [38], and collagen degradation products such as collagen cross-linked N-telopeptide (NTx) are frequently detected at high levels in the blood and urine of osteoporosis patients. Taurine is a non-essential amino acid synthesized from methionine and cysteine by an cysteine sulfinic acid decarboxylase (CSD). Older adults with sarcopenia reportedly exhibit higher levels of several amino acids in serum, including taurine [39]. Sarcopenia is known to be associated with undernutrition [40], and nutritional supplementation with taurine reportedly counteracts the development and progression of sarcopenia in human subjects [41]. We found that taurine levels were positively correlated with walk speed, one of the criteria used to assess sarcopenia, but that association was not statistically significant (r = 0.0378 and p > 0.05). The mechanisms underlying the changes in the serum levels of these amino acid levels before osteoporosis or sarcopenia development remain unclear; however, we conclude that the monitoring of these levels in individuals could be predictive of musculoskeletal disorders.

4. Materials and Methods

4.1. Subjects

This study was conducted using ROAD study cohorts. Subjects were enrolled from resident registration listings in three communities with different characteristics: 1350 subjects were from an urban region in Itabashi, Tokyo; 864 were from a mountainous region in Hidakagawa, Wakayama; and 826 were from a coastal region in Taiji, Wakayama. After the baseline study, the second survey was conducted in the same communities from 2008 to 2010 [42], and the third survey was performed from 2012 to 2013. In the second and third surveys, all examinations for osteoporosis (OP), sarcopenia (SP), and frailty were performed only in cohorts from mountainous and coastal regions. Therefore, here, 1551 participants (521 men and 1030 women) were enrolled from those two cohorts. Blood samples were collected from these subjects in the second survey and stocked. Additionally, from among the 1551 participants, those aged ≥60 years were selected based on the AWGS criteria for sarcopenia [9]. Consequently, 1083 (372 men and 711 women; mean age: 72.1 (SD: 7.4) years (72.7 (SD: 7.5) years for men; 71.7 (SD: 7.3) years for women)) participants were recruited as eligible subjects for the follow-up of sarcopenia. The third survey was carried out in the same subjects between 2012 and 2013, and 767 individuals (253 men and 514 women) participated in both the second and third surveys of the ROAD study among 1083 participants and completed assessment for both OP and SP. Both osteoporosis and sarcopenia were evaluated at the second and third surveys, and blood samples were collected and stored from participants from the mountainous region of Hidakagawa. We analyzed 729 participants (254 men and 475 women) aged 68.7 ± 11.1 years (69.3 ± 10.7 in men and 68.4 ± 11.3 in women) from Hidakagawa, and all completed the second and third follow-up surveys and provided blood samples. Written informed consent were provided by all participants, and the study was conducted with the approval of the ethics committees of the University of Tokyo (Nos. 1264 and 1326) and the University of Wakayama Medical University (No. 373). All procedures were conducted in accordance with the ethical standards as described in the 1964 Declaration of Helsinki and its later amendments.
Osteoporosis and sarcopenia were defined as described in [11]. Briefly, osteoporosis was defined based on bone mineral density (BMD) values, and subjects who exhibited a BMD T score lower than −2.5 SD were judged to be osteoporotic. Sarcopenia was defined based on the criteria of the AWGS [9,11]. New cases of either condition were defined as appearing in subjects who did not meet the osteoporosis or sarcopenia criteria at the second survey but satisfied those criteria at the third survey.

4.2. IGF1 ELISA

Serum IGF1 levels were measured by ELISA according to the manufacturer’s protocol (MG100, R&D Systems, Minneapolis, MN, USA) using a multiple plate analyzer (Powerscan HT, DS Pharma Biomedical, Osaka, Japan).

4.3. Metabolite Extraction from Serum

Frozen serum samples were thawed, and 40 µL aliquots were added to 360 µL of methanol containing internal standards (20 µmol/L each of methionine sulfone and D-camphor-10-sulfonic acid). Solutions were mixed well, and then both 400 µL of chloroform and 160 µL of Milli-Q water were added, followed by centrifugation at 10,000× g for 3 min at 4 °C. Then, 300 µL of the upper layer were transferred to a 5-kDa-cutoff filter (Human Metabolome Technologies, Tsuruoka, Japan) to remove proteins. The filtrate was dried using a centrifuge concentrator and reconstituted with 40 µL of Milli-Q water containing reference compounds (200 µmol/L each of 3-aminopyrrolidine and trimesic acid) prior to CE-TOFMS analysis.

4.4. Metabolome Analysis by CE-TOFMS

All CE-TOFMS analyses were performed using an Agilent 1600 Capillary Electrophoresis system (Agilent technologies, Waldbronn, Germany) equipped with an Agilent 6210 TOF LC/MS system (Agilent technologies, Santa Clara, CA, USA), as described in [5]. For anionic metabolite analysis, the ESI sprayer was replaced with a platinum needle rather than the original stainless steel needle [43]. Cationic metabolites were separated through a fused-silica capillary (50 µm i.d. × 100 cm) filled with 1 mol/L formic acid as an electrolyte [44] and methanol/water (50%, v/v) containing 0.1 µmol/L hexakis (2,2-difluoroethoxy) phosphazene was delivered as a sheath liquid at a rate of 10 µL/min. Capillary temperature was maintained at 20 °C.
Sample solution was injected at 5 kPa for 3 sec, and the separation voltage was set at 30 kV. ESI-TOFMS was conducted in positive ion mode, and capillary, fragmentor, skimmer, and Oct RF voltages were set at 4000, 75, 50, and 125 V, respectively. Nebulizer gas pressure was configured at 7 psig, and heated nitrogen gas (300 °C) was supplied at a rate of 10 L/min. Anionic metabolites were separated through a COSMO (+) capillary (50 µm i.d. × 105 cm, Nacalai Tesque, Kyoto, Japan) filled with 50 mmol/L ammonium acetate (pH 8.5) as an electrolyte [43], and ammonium acetate (5 mmol/L) in 50% (v/v) methanol/water containing 0.1 µmol/L hexakis (2,2-difluoroethoxy) phosphazene was delivered as a sheath liquid at a rate of 10 µL/min. A sample solution was injected at 5 kPa for 30 s, and −30 kV was applied to promote sample separation. ESI-TOFMS was conducted in negative ion mode, and capillary, fragmentor, skimmer, and Oct RF voltages were set at 3500, 100, 50, and 200 V, respectively. Other conditions were identical for cationic metabolite analysis.
In both modes, an automatic recalibration function was used to correct for the analytical variation of exact masses for each run, as described in [45]. Exact mass data were acquired at a rate of 1.5 cycles/s over a 50–1000 m/z range. Raw data were processed using our proprietary automatic integration software (MasterHands) [46,47]. Each peak was identified by matching m/z values and the normalized migration times of corresponding authentic standard compounds. We assessed 238 anionic and 272 cationic metabolites, and we detected a total 80 (as listed in Supplementary Table S1). Metabolite levels are reported as μM.

4.5. Statistical Analysis

The statistical significance of indicated parameters was analyzed between subjects who were not diagnosed with osteoporosis at both the second and third surveys (non-osteoporosis) and those who were not diagnosed with osteoporosis at the second survey but newly developed osteoporosis after the second survey (new osteoporosis). Similarly, the statistical significance of various parameters was analyzed between subjects who were not diagnosed with sarcopenia at both second and third surveys (non-sarcopenia) and those who were not diagnosed with sarcopenia at the second survey but newly developed sarcopenia after that survey (new sarcopenia). All statistical analysis was performed using the STATA software (STATA Corp., College Station, TX, USA).

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/metabo11050278/s1, Table S1: Detected metabolites.

Author Contributions

T.M.; validation, A.H. and T.M.; formal analysis, K.M.; investigation, K.M., Y.S., S.I. and M.M. (Midori Maruyama); resources, N.Y.; data curation, A.H., N.Y. and T.M.; writing—original draft preparation, T.M.; writing—review and editing, A.H. and N.Y.; supervision, T.S., M.T., M.N. and M.M. (Morio Matsumoto); funding acquisition, A.H. and T.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by JOA-Subsidized Science Project Research, grant number 2020-3” and grants from the Yamagata prefectural government and the city of Tsuruoka.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Ethics Committee of the University of Tokyo (No. 1264).

Informed Consent Statement

Written informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All data are contained within the article.

Acknowledgments

We thank Masaru Uragami for statistical analysis. This study was supported in part by Daiichi Sankyo Co., Ltd.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Edwards, M.H.; Dennison, E.M.; Aihie Sayer, A.; Fielding, R.; Cooper, C. Osteoporosis and sarcopenia in older age. Bone 2015, 80, 126–130. [Google Scholar] [CrossRef] [Green Version]
  2. Reginster, J.Y.; Beaudart, C.; Buckinx, F.; Bruyère, O. Osteoporosis and sarcopenia: Two diseases or one? Curr. Opin. Clin. Nutr. Metab. Care 2016, 19, 31–36. [Google Scholar] [CrossRef] [Green Version]
  3. Curtis, E.; Litwic, A.; Cooper, C.; Dennison, E. Determinants of muscle and bone aging. J. Cell Physiol. 2015, 230, 2618–2625. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Yoshimura, N.; Muraki, S.; Oka, H.; Mabuchi, A.; En-Yo, Y.; Yoshida, M.; Saika, A.; Yoshida, H.; Suzuki, T.; Yamamoto, S.; et al. Prevalence of knee osteoarthritis, lumbar spondylosis, and osteoporosis in Japanese men and women: The research on osteoarthritis/osteoporosis against disability study. J. Bone Miner. Metab. 2009, 27, 620–628. [Google Scholar] [CrossRef] [PubMed]
  5. Miyamoto, T.; Hirayama, A.; Sato, Y.; Koboyashi, T.; Katsuyama, E.; Kanagawa, H.; Miyamoto, H.; Mori, T.; Yoshida, S.; Fujie, A.; et al. A serum metabolomics-based profile in low bone mineral density postmenopausal women. Bone 2017, 95, 1–4. [Google Scholar] [CrossRef] [PubMed]
  6. Miyamoto, T.; Hirayama, A.; Sato, Y.; Koboyashi, T.; Katsuyama, E.; Kanagawa, H.; Fujie, A.; Morita, M.; Watanabe, R.; Tando, T.; et al. Metabolomics-based profiles predictive of low bone mass in menopausal women. Bone Rep. 2018, 9, 11–18. [Google Scholar] [CrossRef]
  7. Kanis, J.A.; Oden, A.; Johnell, O.; Johansson, H.; De Laet, C.; Brown, J.; Burckhardt, P.; Cooper, C.; Christiansen, C.; Cummings, S.; et al. The use of clinical risk factors enhances the performance of BMD in the prediction of hip and osteoporotic fractures in men and women. Osteoporos. Int. 2007, 18, 1033–1046. [Google Scholar] [CrossRef]
  8. Cruz-Jentoft, A.J.; Baeyens, J.P.; Bauer, J.M.; Boirie, Y.; Cederholm, T.; Landi, F.; Martin, F.C.; Michel, J.P.; Rolland, Y.; Schneider, S.M.; et al. 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–423. [Google Scholar] [CrossRef] [Green Version]
  9. Chen, L.K.; Liu, L.K.; Woo, J.; Assantachai, P.; Auyeung, T.W.; Bahyah, K.S.; Chou, M.Y.; Chen, L.Y.; Hsu, P.S.; Krairit, O.; et al. Sarcopenia in Asia: Consensus report of the Asian Working Group for Sarcopenia. J. Am. Med. Dir. Assoc. 2014, 15, 95–101. [Google Scholar] [CrossRef]
  10. Bowden Davies, K.A.; Pickles, S.; Sprung, V.S.; Kemp, G.J.; Alam, U.; Moore, D.R.; Tahrani, A.A.; Cuthbertson, D.J. Reduced physical activity in young and older adults: Metabolic and musculoskeletal implications. Ther. Adv. Endocrinol. Metab. 2019, 10. [Google Scholar] [CrossRef]
  11. Yoshimura, N.; Muraki, S.; Oka, H.; Iidaka, T.; Kodama, R.; Kawaguchi, H.; Nakamura, K.; Tanaka, S.; Akune, T. Is osteoporosis a predictor for future sarcopenia or vice versa? Four-year observations between the second and third ROAD study surveys. Osteoporos. Int. 2017, 28, 189–199. [Google Scholar] [CrossRef]
  12. Yoshimura, N.; Muraki, S.; Oka, H.; Kawaguchi, H.; Nakamura, K.; Akune, T. Cohort profile: Research on Osteoarthritis/Osteoporosis Against Disability study. Int. J. Epidemiol. 2010, 39, 988–995. [Google Scholar] [CrossRef] [Green Version]
  13. Dettmer, K.; Aronov, P.A.; Hammock, B.D. Mass spectrometry-based metabolomics. Mass Spectrom. Rev. 2007, 26, 51–78. [Google Scholar] [CrossRef] [PubMed]
  14. Cuevas-Córdoba, B.; Santiago-García, J. Saliva: A fluid of study for OMICS. OMICS 2014, 18, 87–97. [Google Scholar] [CrossRef] [PubMed]
  15. You, Y.S.; Lin, C.Y.; Liang, H.J.; Lee, S.H.; Tsai, K.S.; Chiou, J.M.; Chen, Y.C.; Tsao, C.K.; Chen, J.H. Association between the metabolome and low bone mineral density in Taiwanese women determined by 1H NMR spectroscopy. J. Bone Miner. Res. 2014, 29, 212–222. [Google Scholar] [CrossRef] [PubMed]
  16. Liu, X.; Liu, Y.; Cheng, M.; Zhang, X.; Xiao, H. A metabolomics study of the inhibitory effect of 17-beta-estradiol on osteoclast proliferation and differentiation. Mol. Biosyst. 2015, 11, 635–646. [Google Scholar] [CrossRef]
  17. Qi, H.; Bao, J.; An, G.; Ouyang, G.; Zhang, P.; Wang, C.; Ying, H.; Ouyang, P.; Ma, B.; Zhang, Q. Association between the metabolome and bone mineral density in pre- and post-menopausal Chinese women using GC-MS. Mol. Biosyst. 2016, 12, 2265–2275. [Google Scholar] [CrossRef]
  18. Lee, M.Y.; Kim, H.Y.; Singh, D.; Yeo, S.H.; Baek, S.Y.; Park, Y.K.; Lee, C.H. Metabolite Profiling Reveals the Effect of Dietary Rubus coreanus Vinegar on Ovariectomy-Induced Osteoporosis in a Rat Model. Molecules 2016, 21, 149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Villareal, M.O.; Matsukawa, T.; Isoda, H. L-Citrulline Supplementation-Increased Skeletal Muscle PGC-1α Expression is Associated With Exercise Performance and Increased Skeletal Muscle Weight. Mol. Nutr. Food Res. 2018, 62, e1701043. [Google Scholar] [CrossRef]
  20. Bettis, T.; Kim, B.J.; Hamrick, M.W. Impact of muscle atrophy on bone metabolism and bone strength: Implications for muscle-bone crosstalk with aging and disuse. Osteoporos. Int. 2018, 29, 1713–1720. [Google Scholar] [CrossRef]
  21. Bonewald, L. Use it or lose it to age: A review of bone and muscle communication. Bone 2019, 120, 212–218. [Google Scholar] [CrossRef]
  22. Colaianni, G.; Storlino, G.; Sanesi, L.; Colucci, S.; Grano, M. Myokines and Osteokines in the Pathogenesis of Muscle and Bone Diseases. Curr. Osteoporos. Rep. 2020, 18, 401–407. [Google Scholar] [CrossRef]
  23. Anastasilakis, A.D.; Polyzos, S.A.; Makras, P.; Gkiomisi, A.; Bisbinas, I.; Katsarou, A.; Filippaios, A.; Mantzoros, C.S. Circulating irisin is associated with osteoporotic fractures in postmenopausal women with low bone mass but is not affected by either teriparatide or denosumab treatment for 3 months. Osteoporos. Int. 2014, 25, 1633–1642. [Google Scholar] [CrossRef] [PubMed]
  24. Schnyder, S.; Handschin, C. Skeletal muscle as an endocrine organ: PGC-1α, myokines and exercise. Bone 2015, 80, 115–125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Colaianni, G.; Cuscito, C.; Mongelli, T.; Pignataro, P.; Buccoliero, C.; Liu, P.; Lu, P.; Sartini, L.; Di Comite, M.; Mori, G.; et al. The myokine irisin increases cortical bone mass. Proc. Natl. Acad. Sci. USA 2015, 112, 12157–12162. [Google Scholar] [CrossRef] [Green Version]
  26. Kim, H.; Wrann, C.D.; Jedrychowski, M.; Vidoni, S.; Kitase, Y.; Nagano, K.; Zhou, C.; Chou, J.; Parkman, V.A.; Novick, S.J.; et al. Irisin Mediates Effects on Bone and Fat via αV Integrin Receptors. Cell 2018, 175, 1756–1768. [Google Scholar] [CrossRef] [Green Version]
  27. Takafuji, Y.; Tatsumi, K.; Ishida, M.; Kawao, N.; Okada, K.; Kaji, H. Extracellular vesicles secreted from mouse muscle cells suppress osteoclast formation: Roles of mitochondrial energy metabolism. Bone 2020, 134, 115298. [Google Scholar] [CrossRef]
  28. Xian, L.; Wu, X.; Pang, L.; Lou, M.; Rosen, C.J.; Qiu, T.; Crane, J.; Frassica, F.; Zhang, L.; Rodriguez, J.P.; et al. Matrix IGF-1 maintains bone mass by activation of mTOR in mesenchymal stem cells. Nat. Med. 2012, 18, 1095–1101. [Google Scholar] [CrossRef] [Green Version]
  29. Tando, T.; Hirayama, A.; Furukawa, M.; Sato, Y.; Kobayashi, T.; Funayama, A.; Kanaji, A.; Hao, W.; Watanabe, R.; Morita, M.; et al. Smad2/3 Proteins Are Required for Immobilization-induced Skeletal Muscle Atrophy. J. Biol. Chem. 2016, 291, 12184–12194. [Google Scholar] [CrossRef] [Green Version]
  30. Nakamura, S.; Sato, Y.; Kobayashi, T.; Oike, T.; Kaneko, Y.; Miyamoto, K.; Funayama, A.; Oya, A.; Nishiwaki, T.; Matsumoto, M.; et al. Insulin-like growth factor-I is required to maintain muscle volume in adult mice. J. Bone Miner. Metab. 2019, 37, 627–635. [Google Scholar] [CrossRef] [PubMed]
  31. Tang, Y.; Wu, X.; Lei, W.; Pang, L.; Wan, C.; Shi, Z.; Zhao, L.; Nagy, T.R.; Peng, X.; Hu, J.; et al. TGF-beta1-induced migration of bone mesenchymal stem cells couples bone resorption with formation. Nat. Med. 2009, 15, 757–765. [Google Scholar] [CrossRef] [Green Version]
  32. Waning, D.L.; Mohammad, K.S.; Reiken, S.; Xie, W.; Andersson, D.C.; John, S.; Chiechi, A.; Wright, L.E.; Umanskaya, A.; Niewolna, M.; et al. Excess TGF-β mediates muscle weakness associated with bone metastases in mice. Nat. Med. 2015, 21, 1262–1271. [Google Scholar] [CrossRef]
  33. Wilmanski, T.; Rappaport, N.; Earls, J.C.; Magis, A.T.; Manor, O.; Lovejoy, J.; Omenn, G.S.; Hood, L.; Gibbons, S.M.; Price, N.D. Blood metabolome predicts gut microbiome α-diversity in humans. Nat. Biotechnol. 2019, 37, 1217–1228. [Google Scholar] [CrossRef]
  34. Wang, T.J.; Larson, M.G.; Vasan, R.S.; Cheng, S.; Rhee, E.P.; McCabe, E.; Lewis, G.D.; Fox, C.S.; Jacques, P.F.; Fernandez, C.; et al. Metabolite profiles and the risk of developing diabetes. Nat. Med. 2011, 17, 448–453. [Google Scholar] [CrossRef] [PubMed]
  35. Gonzalez-Freire, M.; Moaddel, R.; Sun, K.; Fabbri, E.; Zhang, P.; Khadeer, M.; Salem, N., Jr.; Ferrucci, L.; Semba, R.D. Targeted Metabolomics Shows Low Plasma Lysophosphatidylcholine 18:2 Predicts Greater Decline of Gait Speed in Older Adults: The Baltimore Longitudinal Study of Aging. J. Gerontol. A Biol. Sci. Med. Sci. 2019, 74, 62–67. [Google Scholar] [CrossRef]
  36. Pernow, Y.; Thorén, M.; Sääf, M.; Fernholm, R.; Anderstam, B.; Hauge, E.M.; Hall, K. Associations between amino acids and bone mineral density in men with idiopathic osteoporosis. Bone 2010, 47, 959–965. [Google Scholar] [CrossRef] [PubMed]
  37. Jennings, A.; MacGregor, A.; Spector, T.; Cassidy, A. Amino Acid Intakes Are Associated With Bone Mineral Density and Prevalence of Low Bone Mass in Women: Evidence From Discordant Monozygotic Twins. J. Bone Miner. Res. 2016, 31, 326–335. [Google Scholar] [CrossRef] [Green Version]
  38. de Paz-Lugo, P.; Lupiáñez, J.A.; Meléndez-Hevia, E. High glycine concentration increases collagen synthesis by articular chondrocytes in vitro: Acute glycine deficiency could be an important cause of osteoarthritis. Amino Acids 2018, 50, 1357–1365. [Google Scholar] [CrossRef] [Green Version]
  39. Calvani, R.; Picca, A.; Marini, F.; Biancolillo, A.; Gervasoni, J.; Persichilli, S.; Primiano, A.; Coelho-Junior, H.J.; Bossola, M.; Urbani, A.; et al. A Distinct Pattern of Circulating Amino Acids Characterizes Older Persons with Physical Frailty and Sarcopenia: Results from the BIOSPHERE Study. Nutrients 2018, 10, 1691. [Google Scholar] [CrossRef] [Green Version]
  40. Sieber, C.C. Malnutrition and sarcopenia. Aging Clin. Exp. Res. 2019, 31, 793–798. [Google Scholar] [CrossRef]
  41. Scicchitano, B.M.; Sica, G. The Beneficial Effects of Taurine to Counteract Sarcopenia. Curr. Protein Pept. Sci. 2018, 19, 673–680. [Google Scholar] [CrossRef]
  42. Yoshimura, N.; Oka, H.; Muraki, S.; Akune, T.; Hirabayashi, N.; Matsuda, S.; Nojiri, T.; Hatanaka, K.; Ishimoto, Y.; Nagata, K.; et al. Reference values for hand grip strength, muscle mass, walking time, and one-leg standing time as indices for locomotive syndrome and associated disability: The second survey of the ROAD study. J. Orthop. Sci. 2011, 16, 768–777. [Google Scholar] [CrossRef] [PubMed]
  43. Soga, T.; Igarashi, K.; Ito, C.; Mizobuchi, K.; Zimmermann, H.P.; Tomita, M. Metabolomic profiling of anionic metabolites by capillary electrophoresis mass spectrometry. Anal. Chem. 2009, 81, 6165–6174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Soga, T.; Heiger, D.N. Amino acid analysis by capillary electrophoresis electrospray ionization mass spectrometry. Anal. Chem. 2000, 72, 1236–1241. [Google Scholar] [CrossRef] [PubMed]
  45. Hirayama, A.; Sugimoto, M.; Suzuki, A.; Hatakeyama, Y.; Enomoto, A.; Harada, S.; Soga, T.; Tomita, M.; Takebayashi, T. Effects of processing and storage conditions on charged metabolomics profilies in blood. Electrophoresis 2015, 36, 2148–2155. [Google Scholar] [CrossRef]
  46. Sugimoto, M.; Wong, D.T.; Hirayama, A.; Soga, T.; Tomita, M. Capillary electrophoresis mass spectrometry-based saliva metabolomics identified oral, breast and pancreatic cancer-specific profiles. Metabolomics 2010, 6, 78–95. [Google Scholar] [CrossRef] [Green Version]
  47. Sugimoto, M.; Kawakami, M.; Robert, M.; Soga, T.; Tomita, M. Bioinformatics Tools for Mass Spectroscopy-Based Metabolomic Data Processing and Analysis. Curr. Bioinform. 2012, 7, 96–108. [Google Scholar] [CrossRef]
Table
Table 1. Comparison of characteristics between non- and new-osteoporosis subjects at the second survey.
Table 1. Comparison of characteristics between non- and new-osteoporosis subjects at the second survey.
Non-OsteoporosisNew Osteoporosisp Value
Age (years)65.2 ± 10.869.8 ± 8.00.014
Female sex (%)60.080.00.0181
Height (cm)156.0 ± 8.8150.6 ± 8.00.0004
Weight (kg)57.1 ± 9.450.0 ± 6.6<0.0001
BMI (kg/cm2)23.4 ± 3.022.0 ± 2.50.0076
Waist (cm)82.6 ± 8.679.2 ± 9.90.0295
BMD (L2–4) (g/cm2)0.993 ± 0.1720.818 ± 0.147<0.0001
BMD (femoral neck) (g/cm2)0.701 ± 0.1210.562 ± 0.355<0.0001
Grip strength (maximum) (kg)32.1 ± 9.226.4 ± 6.80.0004
Usual walking speed (m/s)1.13 ± 0.241.09 ± 0.230.3843
Sarcopenia at the second survey (%)2.45.70.2314
New sarcopenia at the third survey (%)5.815.20.0353
Non-osteoporosis: subjects who were not diagnosed with osteoporosis at both second and third surveys; new osteoporosis: subjects who were not diagnosed with osteoporosis at the second survey but newly developed osteoporosis after the second survey. Statistical significance was evaluated by a t test.
Table
Table 2. Levels of IGF1 and metabolites in non- and new-osteoporosis subjects at the second survey.
Table 2. Levels of IGF1 and metabolites in non- and new-osteoporosis subjects at the second survey.
Non-OsteoporosisNew Osteoporosisp Value
IGF1 (ng/mL)86.0 ± 31.071.2 ± 12.80.0052
5-Oxoproline81.6 ± 39.796.0 ± 47.70.0422
2-Oxoglutarate13.1 ± 8.68.6 ± 6.50.0027
cis-Aconitate3.61 ± 0.873.98 ± 0.780.0154
Gly486.0 ± 113.6550.1 ± 140.80.0017
3-Aminoisobutyrate2.10 ± 2.163.07 ± 2.020.0109
Taurine252.0 ± 80.1283.1 ± 126.10.0364
Pipecolate2.60 ± 3.094.99 ± 15.850.0099
Non-osteoporosis: subjects who were not diagnosed with osteoporosis at both second and third surveys; new osteoporosis: subjects who were not diagnosed with osteoporosis at the second survey but newly developed osteoporosis after the second survey. Statistical significance was evaluated by a t test. Statistical significance was evaluated by t test.
Table
Table 3. Logistic regression model for new osteoporosis development.
Table 3. Logistic regression model for new osteoporosis development.
p Value after Adjustment for Agep Value after Adjustment for Age, Sex, and BMIOdds Ratio95% Confidence Interval
BMD (L2–4) (g/cm2)<0.001<0.0010.00049.72 × 10−6–0.161
BMD (femoral neck) (g/cm2)<0.001<0.0013.53 × 10−85.47 × 10−11–0.00002
Grip strength (maximum) (kg)0.0030.174
IGF10.088
5-Oxoproline0.07
2-Oxoglutarate0.0030.055
cis-Aconitate0.088
Gly0.0010.0211.0031.000–1.006
3-Aminoisobutyrate0.066
Taurine0.0380.127
Pipecolate0.095
Statistical significance was evaluated by logistic regression analysis.
Table
Table 4. Comparison of characteristics between non- and new-sarcopenia subjects at the second survey.
Table 4. Comparison of characteristics between non- and new-sarcopenia subjects at the second survey.
Non-SarcopeniaNew Sarcopeniap Value
Age (years)64.9 ± 10.575.9 ± 5.7<0.0001
Female sex (%)65.164.10.8983
Height (cm)155.2 ± 8.9151.7 ± 8.80.0197
Weight (kg)56.5 ± 9.548.6 ± 7.1<0.0001
BMI (kg/cm2)23.4 ± 3.021.1 ± 2.4<0.0001
Waist (cm)82.5 ± 8.577.4 ± 8.90.0004
BMD (L2–4) (g/cm2)0.959 ± 0.1890.877 ± 0.2280.0106
BMD (femoral neck) (g/cm2)0.685 ± 0.1250.605 ± 0.1340.0002
Grip strength (maximum) (kg)31.6 ± 9.226.1 ± 7.00.0003
Usual walking speed (m/s)1.13 ± 0.230.98 ± 0.15<0.0001
Osteoporosis at the second survey (%)10.646.2<0.0001
New osteoporosis at the third survey (%)6.717.20.0353
Non-sarcopenia: subjects who were not diagnosed with sarcopenia at both second and third surveys; new sarcopenia: subjects who were not diagnosed with sarcopenia at the second survey but newly developed sarcopenia after the second survey. Statistical significance was evaluated by a t test.
Table
Table 5. Levels of IGF1 and metabolites in non- and new-sarcopenia subjects at the second survey.
Table 5. Levels of IGF1 and metabolites in non- and new-sarcopenia subjects at the second survey.
Non-SarcopeniaNew Sarcopeniap Value
IGF1 (ng/mL)85.0 ± 30.871.5 ± 24.50.0077
2-Hydroxybutyrate27.5 ± 11.623.5 ± 9.20.0376
4-Methyl2oxopentanoate50.3 ± 14.044.0 ± 13.80.0071
2-AB17.1 ± 5.614.9 ± 5.30.0221
Val350.7 ± 72.1326.0 ± 70.40.0040
Ile90.9 ± 26.982.0 ± 21.40.0448
Leu196.5 ± 49.6175.4 ± 46.00.0107
Taurine257.6 ± 86.5227.5 ± 46.30.0328
Trp72.9 ± 13.666.9 ± 14.10.0079
Cystine42.1 ± 8.638.5 ± 7.80.0122
Non-sarcopenia: subjects who were not diagnosed with sarcopenia at both second and third surveys; new sarcopenia, subjects who were not diagnosed with sarcopenia at the second survey but newly developed sarcopenia after the second survey. Statistical significance was evaluated by a t test.
Table
Table 6. Logistic regression model for new sarcopenia development.
Table 6. Logistic regression model for new sarcopenia development.
p Value after Adjustment for Agep Value after Adjustment for Age, Sex, and BMIOdds Ratio95% Confidence Interval
BMD (L2–4) (g/cm2)0.076
BMD (femoral neck) (g/cm2)0.066
Grip strength (maximum) (kg)0.071
Usual walking speed (m/s)0.263
Osteoporosis at the second survey (%)<0.0010.0183.1181.215–7.997
IGF10.915
2-Hydroxybutyrate0.093
4-Methyl2oxopentanoate0.093
2-AB0.077
Val0.06
Ile0.06
Leu0.074
Taurine0.0380.0390.9910.982–0.999
Trp0.102
Cystine0.030.108
Statistical significance was evaluated by logistic regression analysis.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Miyamoto, K.; Hirayama, A.; Sato, Y.; Ikeda, S.; Maruyama, M.; Soga, T.; Tomita, M.; Nakamura, M.; Matsumoto, M.; Yoshimura, N.; et al. A Metabolomic Profile Predictive of New Osteoporosis or Sarcopenia Development. Metabolites 2021, 11, 278. https://doi.org/10.3390/metabo11050278

AMA Style

Miyamoto K, Hirayama A, Sato Y, Ikeda S, Maruyama M, Soga T, Tomita M, Nakamura M, Matsumoto M, Yoshimura N, et al. A Metabolomic Profile Predictive of New Osteoporosis or Sarcopenia Development. Metabolites. 2021; 11(5):278. https://doi.org/10.3390/metabo11050278

Chicago/Turabian Style

Miyamoto, Kana, Akiyoshi Hirayama, Yuiko Sato, Satsuki Ikeda, Midori Maruyama, Tomoyoshi Soga, Masaru Tomita, Masaya Nakamura, Morio Matsumoto, Noriko Yoshimura, and et al. 2021. "A Metabolomic Profile Predictive of New Osteoporosis or Sarcopenia Development" Metabolites 11, no. 5: 278. https://doi.org/10.3390/metabo11050278

APA Style

Miyamoto, K., Hirayama, A., Sato, Y., Ikeda, S., Maruyama, M., Soga, T., Tomita, M., Nakamura, M., Matsumoto, M., Yoshimura, N., & Miyamoto, T. (2021). A Metabolomic Profile Predictive of New Osteoporosis or Sarcopenia Development. Metabolites, 11(5), 278. https://doi.org/10.3390/metabo11050278

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop