Muscle plays a more superior role than fat in bone homeostasis: A cross-sectional study of old Asian people

Objectives: The aim of this study was to discover the role of fat and muscle in bone structures, as well as the relationship between obesity and sarcopenia on age-related osteoporosis.

Methods: A total of 400 participants (65.0 ± 8.2 years old, 42.3% women) were recruited. Fat, muscle, bone parameters, basic demographics, medical history, physical performance and activity, and calcium intake of participants were obtained from datasets. The diagnosis of osteoporosis, sarcopenia, and obesity was based on current recommendations. Pearson correlation, non-linear regression models, and decision tree analyses were performed to study the relationship between fat, muscle, and bone. Logistic regression analyses were used to explore the risk of osteoporosis in old people with obesity or sarcopenia via Model 1 (unadjusted) and Model 2 (adjusted by age, physical activity, and calcium intake).

Results: Correlation analysis showed that limb muscle mass and index, and age were best related to bone mineral density (BMD) (|r| = 0.386–0.632, p < 0.001). On the contrary, body mass index (BMI) and increased body fat percentage (BF%) were harmful for bone health. An increase of BMI and fat mass index slowed the increase of BMD in the spine, while skeletal muscle mass index accelerated the increase. People with sarcopenia had low muscle mass and strength. When separating subjects into sarcopenia and non-sarcopenia status, sarcopenia was independently related to higher risks of osteoporosis in both models (OR > 1, p < 0.05). BMI-defined obesity in Model 1 as well as BF%-defined obesity in both models did not reduce the risk of osteoporosis in both models (p > 0.05). The decision tree classification (85% accuracy) showed that greater body weight and larger lower limb muscle performance were negatively related to osteoporosis, while fat mass and percentage did not play roles in this prediction.

Conclusion: Low muscle mass and function were harmful to bone health. Obesity defined by both BMI and BF% had limited protective roles in osteoporosis. The benefits for bone from increased muscle mass and function play a more superior role than increased fat mass in old people. Sarcopenia prevention and treatment instead of controlling obesity should be recommended as an approach to reduce the risks of age-related osteoporosis and fragility fracture for elderly people.

1. Introduction

Osteoporosis is characterized by low bone mineral density (BMD) and leads to a high risk of fracture in the elderly (1, 2). It is estimated that patients with fragility fractures will reach 6.3 million by 2050, posing a huge public health burden (3). The annual costs for fragility fracture patients will also increase 27% from 2017 to 2030 in Europe (4). In the United States (US), an average of US$3,000 per patient is required in the first year after a fracture (5).

The “obesity paradox” in osteoporosis has been reported previously, as people usually have a higher BMD due to increased mechanical loading (6, 7). In fact, studies have shown that for individuals with overweight or obesity, body mass index (BMI) reduces the risk of a fracture (8, 9). However, when body composition of fat and lean mass is taken into consideration, body fat mass is negatively associated with BMD in the younger population, while lean mass plays a positive role (10). In animal models, diet-induced obesity impairs bone turnover through cytokines and osteopontin secreted by adipose tissues, which activate osteoclasts and accelerate bone degradation (11). On the other hand, strain induced by muscle contraction stimulates bone growth as osteoblasts are mechanosensitive (6, 12). The skeletal muscle is one of the major components of lean mass. Muscle-secreted myokines including interleukins, irisin, myostatin, and growth factors can regulate bone metabolism (13). The composition of muscle and fat, as well as mechanical loadings, contributes to the bidirectional relationship of bone homeostasis.

The coexistence of osteoporosis and sarcopenia is prevalent in elderly people (14). Sarcopenia is an age-related muscle disorder, which is strongly related to fragility fractures in the elderly (15, 16). More importantly, it is currently unclear whether people with obesity have better bone microarchitecture compared to those with normal BMI after excluding individuals with sarcopenia. Dividing subjects into normal, obesity, sarcopenia, and sarcopenic obesity phenotypes can help to explore if the “obesity paradox” still exists in people with obesity when compared to healthy controls (17).

The aims of this study were to (1) investigate the correlation between muscle/fat and bone parameters (2); assess the osteoporosis risk according to normal, obesity, sarcopenia, and sarcopenic obesity groups; and (3) identify whether fat or muscle has a greater role in osteoporosis prediction.

2. Materials and methods

2.1. Study design and population

This study included the population aged 50 years or older from the established Normal Reference Study (NRS) cohort in May 2015 to April 2017. The aim of the cohort establishment was to study the reference centile curves of high-resolution peripheral quantitative computed tomography (HR-pQCT) parameters in Chinese adults in Hong Kong. The inclusion and exclusion criteria have been described previously (18). In brief, participants were physically healthy independent Chinese individuals, without medical conditions that affect bone metabolism, history of fragility fracture, or genetic skeletal disorders. Basic demographics, calcium intake, Baecke physical activity questionnaire, datasets of HR-pQCT, DXA, and medical history from each participant were obtained. The study procedure was approved by the local clinical research ethics committee (CRE.2014.310). All participants had written informed consent.

2.2. HR-pQCT assessment and image analysis

The distal radius and distal tibia (non-dominant and non-fractured) of participants were scanned by HR-pQCT (XtremeCT I, Scanco Medical AG, Switzerland) using the standard scanning protocol (60 kVp, 900 μA, 100 ms integration time, isotropic voxel size 82 μm). The details of the HR-pQCT operational procedures and the image quality grading of HR-pQCT for motion artifact were illustrated in previous studies (18, 19). Daily phantom calibration was performed to monitor the stability of the HR-pQCT system. Precision error was presented by coefficient of variation (CV). In the cohort, CV of vBMD ranged from 0.26% to 0.88%, and for microarchitectural parameters, it ranged from 0.56% to 3.78% (20). The proportion of good-quality image was 79.3% in the NRS cohort (18). Participants with poor-quality images including obvious image artifacts had been excluded in our study (21).

For image analysis, the volume of interest was differentiated into cortical and trabecular components using a fully automated cortical compartment segmentation technique (22). For bone density analysis, volumetric bone mineral density [vBMD, mg hydroxyapatite (HA)/cm³] for the whole bone, trabecular bone, and cortical bone were calculated. For trabecular bone analysis, bone volume fraction (BV/TV, %) was calculated from trabecular vBMD and the density of fully mineralized bone was assumed to be 1,200 mg HA/cm³ (23). Trabecular thickness (Tb.Th, mm) was derived by standard histomorphometry methods (24). For cortical bone analysis, cortical thickness was measured directly by a distance-transform method and intracortical pores in the binary cortex image were eliminated (19). Bone strength was estimated in terms of stiffness (kN/mm) by micro-finite element analyses (μFEA) on three-dimensional images of the distal radius and tibia using a built-in Image Processing Language software of HR-pQCT (IPL-FE v1.15, Scanco Medical). The specific conditions of the software operation were illustrated in previous studies (18, 19). Whole bone stiffness (kN/mm) was evaluated by performing a uniaxial compression test with 1000 N load and an apparent strain of 1%.

2.3. DXA, muscle function, and physical performance assessments

Body composition was measured by dual-energy x-ray absorptiometry (DXA) (Horizon; Hologic, Bedford, MA, USA). Scanning and analysis were performed under standard operating procedure by qualified technicians. Fat and lean mass (kg) in limbs, trunk, and whole body were measured and recorded. The ratio, percentage (divided by weight in kilograms), and index (normalized by height in meters squared) were calculated. Areal BMD (aBMD, g/cm²) at the proximal femur (femoral neck and total hip) and lumbar spine (L1 to L4) was evaluated. T-score and Z-score for total hip, lumbar spine, and femoral neck were evaluated. Daily phantom calibration was performed to monitor the stability of the DXA system. CV for femoral neck, total hip, and lumbar spine aBMD in our center was 1.36%, 1.19%, and 1.01%, respectively (19). Five elderly volunteers were included to assess the reliability of DXA examination. Briefly, they were asked to perform the examination at baseline and after 6 weeks, and not to participate in interventional research during this period. The intraclass correlation coefficient (ICC) with 95% CI of DXA-measured limb lean mass, total fat mass, aBMD of femoral neck, total hip, and lumbar spine was 0.915 (0.507, 0.991), 0.988 (0.888, 0.999), 0.972 (0.790, 0.997), 0.943 (0.546, 0.994), 0.992 (0.726, 0.999) (p≤0.004) .

Handgrip strength (HGS) was measured by the Smedley dynamometer (model EH101, Camry). Participants were in standing position with elbows in full extension and grasped the dynamometer as hard as they could for three times on both hands (30-s rest between each trial) (25). The highest reading among all trials was taken for analysis. Lower limb strength was also measured by the dynamometer (model EH101, Camry). The test aimed to examine the isometric strength of the quadriceps and hip extensors predominately (26). Participants were in sitting position on a highchair with both legs hung loose. The ankle was tied with the straps that connect with the dynamometer with knee angle at 90°C flexion. The participants then raised the leg forward as hard as they could for three times on both legs (30-s rest between each trial). The highest reading was taken for analysis. Physical performance was assessed by gait speed (GS). Participants walked a total distance of 8 m at usual speed and the time for walking the middle 6-m distance was measured by a manual stopwatch. The GS (m/s) was calculated as the distance (6 m) divided by the average time of the two trials for analysis. Five elderly volunteers were also included to assess the reliability of these tests. The reliability of HGS, lower limb strength, and GS was 0.985 (0.883, 0.998), 0.942 (0.575, 0.994), and 0.990 (0.916, 0.999) (p ≤ 0.003).

2.4. Diagnosis of osteoporosis, sarcopenia, and obesity

With reference to the T-score, BMD was defined as osteoporosis when T-score ≤ −2.5 in either femoral neck, total hip, or total lumbar spine (27). To diagnose sarcopenia, the Asian Working Group for Sarcopenia (AWGS) 2019 criteria were used (28). Those with low muscle mass and low HGS/GS were identified as sarcopenia. DXA-measured appendicular skeletal muscle mass index (ASMI) <7.0 kg/m² in male subjects and <5.4 kg/m² in female subjects were regarded as low muscle mass. Low muscle strength was indicated by HGS <28 kg in men and <18 kg in women, and low physical performance was indicated by GS <1 m/s. A recent consensus recommends using BF% to define obesity, which not only regarded the body composition assessment but was shown to be more suitable for elderly people (29). Therefore, two obesity definitions were utilized, BMI ≥ 25 kg/m² and BF% > 60th percentile of the current cohort (male subjects > 32.7%, female subjects > 43.2%) for assessment (29–31). Those with the presence of sarcopenia and obesity were defined as “sarcopenic obesity (SO)”.

2.5. Calcium intake and physical activity assessments

Daily calcium intake was measured by a validated local food frequency questionnaire (32). Participants were asked to report the frequency and usual consumption amount of each food item in the past year by a trained interviewer; memory recall of the food consumption pattern was aided by photo booklet and different portion reference tools. The total calcium content from the reported food was calculated for analysis. The habitual physical activity over the previous 12 months was assessed by the modified Baecke questionnaire (33). Several questions in each of four categories (work, sport, leisure time, and housework) were asked with score ranging from 1 (never) to 5 (always or often). The total score from the four categories was calculated for the analysis of physical activity assessment.

2.6. Statistical analysis

The consistency and stability of DXA, muscle function, and performance were analyzed by ICC. To study the correlation between muscle/fat and bone structure at different sites, we performed Pearson correlation coefficient test and demonstrated a heatmap according to fat and muscle indicators, as well as bone parameters. Muscle indicators included physical activity and performance. Bone parameters were from the hip and spine by DXA, and distal radius and tibia by HR-pQCT. The block color closer to white represented nonsignificance. The correlation coefficient (r) of 0–0.19, 0.2–0.39, 0.4–0.59, 0.6–0.79, and 0.8–1.0 was regarded as very weak, weak, moderate, strong, and very strong correlation, respectively. Nonparametric regression models were executed to display the bone parameter in different bone sites with increasing BMI, ASMI, or fat mass index (FMI). One-way ANOVAs with Bonferroni tests were performed to compare the muscle, fat, and bone parameters of normal, obese, sarcopenic, and SO groups. Continuous variables were exhibited as mean and standard deviation (SD). Odds ratio (OR) of osteoporosis, osteopenia, lower radius and tibial vBMD of participants with obesity, sarcopenia, and SO were performed with logistic regression. Normal people were the reference group (OR = 1.00). Briefly, the logistic regression model was chosen, and the diagnosis of osteoporosis was the outcome variable. In Model 1 (unadjusted), only sarcopenia/obesity/SO and normal status were covariates. In Model 2 (adjusted), age, physical activity score, and calcium intake were also considered as covariates. The odds ratio was read from EXP(B) with 95% CI. Low vBMD was defined as lower than 2 SD of the local young reference (20 to 35 years old) mean value. According to the 240 young subjects from the NRS cohort, the cutoffs of lower radial vBMD were 268.7 mg HA/cm³ in male subjects and 267.4 mg HA/cm³ in female subjects, and for tibia, they were 230.5 mg HA/cm³ in male subjects and 248.9 mg HA/cm³ in female subjects. We used a machine learning method (decision tree classification) to predict osteoporosis. Basic demographics, muscle and fat indicators, calcium intake, and fall and medical history were the selected attributes. The ratio of training sets and test sets was 7:3 (34). The specificity, sensitivity, and accuracy of the classification were detected. Python 3.10.1, R 4.0.2, and SPSS Statistics 20 were utilized for all analyses and visualizations. Significant difference was defined as p < 0.05.

3. Results

A total of 400 subjects (42.3% of which were female) from 50 to 84 (mean ± SD: 65.0 ± 8.2) years old were included. Thirteen subjects had missed at least one of the physical performance or activity data, and 38 subjects missed the calcium intake calculation. These participants were not included in the odds ratio and decision tree analyses. All participants had DXA and HR-pQCT parameters.

3.1. The effects of fat and muscle parameters on different sites of bone

The correlation heatmap is shown in Figure 1. For annotations of all parameters, please refer to Supplementary Table S1. Increased age was negatively correlated with bone homeostasis, displayed by the lower BMD at the sites of hip, lumbar spine, tibia, and radius, as well as BV/TV and cortical bone thickness in both tibia and radius, and trabecular bone thickness in radius. The radial cortical bone vBMD was the most relevant parameter, which showed moderately negative correlation with age (r = −0.473, p < 0.001). BMI was positively correlated to various bone parameters, particularly the weight-bearing hip aBMD (r = 0.400, p < 0.001), but it showed no correlation with trabecular bone thickness in both radius and tibia (p > 0.05). Another obesity-related indicator was body fat percentage (BF%). The majority of the correlation coefficient in BF% was opposite to BMI findings. The stiffness of the tibia and radius was reduced with higher BF% (r ≤ −0.563, p < 0.001). Weak negative correlations were found in BMD and BF% in the tibia, hip, and lumbar spine (r = −0.297 to −0.193, p < 0.001), but very weak in radius (r = −0.176, p < 0.001). The absolute body fat mass was positively associated with hip, lumbar spine, and tibial BMD instead of radius. Considering fat distribution, we observed that limb fat mass was only correlated with femoral neck aBMD, but not BMD in other sites (r = 0.165, p = 0.001). Radius and tibial stiffness and trabecular bone thickness as two fracture predictors were negatively correlated with limb fat mass (35). Conversely, trunk fat mass was positively associated with BMD in all sites, as well as trabecular bone thickness and tibial stiffness (p < 0.05). Trunk-to-limb fat ratio showed positive associations with most bone parameters except for cortical bone vBMD in radius.

FIGURE 1

Figure 1 Correlation heatmap of muscle, fat, and bone indicators. The color intensity of the block represented the strong or weak correlation between two parameters. Blue blocks showed the positive correlation (r = 0–0.8), while the red showed the negative correlation (r = −0.8–0). The color closer to white could be regarded as no statistical significance (p > 0.05). TL, trunk-to-limb; LFR, lean-to-fat ratio; T, tibia; R, radius.

Limb muscle mass and ASMI were weakly associated with radial vBMD (r = 0.335, p < 0.001) and strongly associated with total hip aBMD (r = 0.605–0.632, p < 0.001). Both muscle indicators showed obviously positive correlations with bone parameters, especially the stiffness of peripheral bone (r = 0.685–0.803, p < 0.001). However, they were not associated with tibial trabecular bone thickness (p > 0.05). The lean-to-fat ratio in limb and total body had similar results. Additionally, the functional performance was positively correlated with most bone parameters. Limb muscle mass and index had the strongest positive correlation parameters for radial vBMD, and age had the largest negative correlation (p < 0.001). ASMI was most positively related to tibial vBMD, while age and BF% were negatively related to it (p < 0.001). For both total hip and lumbar spine aBMD, limb lean mass has the strongest positive correlation coefficient and BF% had the strongest negative one (p < 0.001).

3.2. The increase rate of BMD with rising BMI, ASMI, and FMI

It has been reported that BMI and fracture risk have a non-linear relationship (36). We used second-order non-linear regression models to explore the relationships between BMI, ASMI (eliminate influence of organs in trunk), FMI, and BMD in different bone sites, shown in Figure 2. With increasing BMI, a slower increase rate of BMD in male subjects at all bone sites and female subjects at hip and lumbar spine was found when BMI was over approximately 25 kg/m² (quadratic coefficient < 0). However, the rising trend of radial and tibial BMD in female subjects did not change significantly (Figures 2A, D). To observe whether the role of different body compositions in BMD was similar to BMI, we used ASMI as muscle index, and FMI as fat index. An increase in ASMI accelerated BMD enhancement in hip and spine in male subjects, and spine in female subjects (quadratic coefficient > 0) (Figures 2B, E). However, when BMI and ASMI increased to a certain threshold (31.9 kg/m² and 7.3 kg/m²), radius BMD started to decrease in male subjects. Consistent with BMI, as FMI increases, there is also a slower rate of increase for hip and spine BMD in both genders, as well as radius BMD in male subjects (quadratic coefficient < 0) (Figures 2C, F). Increasing FMI benefits the tibial BMD increment in both genders (quadratic coefficient > 0).

FIGURE 2

Figure 2 Non-linear regression models of BMD in total hip, lumbar spine, radius, and tibia with the increased body weight, muscle, and fat index. (A–C) Total hip (blue) and lumbar spine (orange) aBMD changes with BMI, ASMI, and FMI in male subjects (solid line) and female subjects (dashed line). (D–F) Radius (blue) and tibias (orange) vBMD changes with BMI, ASMI, and FMI in male subjects (solid line) and female subjects (dashed line). TH, total hip; LS, lumbar spine; R, radius; T, tibia.

3.3. Osteoporosis risk according to sarcopenic and obese status

We separated the older participants into normal, obese, sarcopenic, and sarcopenic obese groups according to sarcopenia and obesity definitions. Tables 1, 2 show the bone parameters in each group when obesity was defined by BMI and BF%, respectively.

TABLE 1

Table 1 Muscle, fat, and bone parameters in normal, obese, sarcopenic, and SO groups when obesity was defined by BMI.

TABLE 2

Table 2 Muscle, fat, and bone parameters in normal, obese, sarcopenic, and SO groups when obesity was defined by BF%.

When obesity was defined by BMI (Table 1), hip aBMD was lowest in people with sarcopenia and highest in people with obesity in both genders (p < 0.05). Male individuals with SO had higher total hip aBMD than individuals with sarcopenia, but female individuals with SO had a similar value compared to the other three groups. Lumbar spine aBMD significantly increased in female subjects with obesity than in normal individuals with sarcopenia, but male subjects had comparable values among the four groups. Poor muscle status did not contribute to the BMD at this site (p > 0.05). Obese and normal groups had similar radial vBMD, and people with sarcopenia had lower values than individuals with obesity (p < 0.05). Similarly, other radius parameters were comparable between obese and normal groups, and most of the parameters were remarkedly decreased in people with sarcopenia compared to those with obesity. In male tibia, the total distal tibial vBMD and cortical bone thickness were increased in individuals with obesity, and sarcopenic and normal groups had similar values. Nevertheless, female subjects had higher tibial cortical bone vBMD in the obese compared to normal groups, but not the total vBMD. Female subjects with sarcopenia had lowest tibial vBMD, cortical bone thickness, and BV/TV than both normal and obese groups. Radial and tibial stiffness values significantly decreased in people with sarcopenia (p < 0.05). Higher stiffness in the obese group was only found in female tibia compared to the normal group. Older individuals with sarcopenic obesity generally had comparable bone parameters compared to other groups. When obesity was defined by BF% (Table 2), male subjects with obesity had increased ASMI, which was not found in female subjects. Limb lean mass and BMD in all bone sites were similar between obese and normal groups (p > 0.05). Compared to the normal group, the sarcopenic group had lower hip aBMD. Decreased spine, radius, and tibia BMD were only found in female subjects with sarcopenia, while male subjects with SO had reduced radius vBMD (p < 0.05).

To detect the effects of sarcopenia and obesity on impaired bone homeostasis, ORs with 95% confidence interval (CI) of osteoporosis, low radius and tibial vBMD, and fall history according to the above grouping are shown in Tables 3 and 4 based on different obesity definitions. Before the model adjustment, BMI-defined obesity was associated with lower risks of osteoporosis and low tibial vBMD (OR < 1.0, p < 0.05), but not low radial vBMD. Sarcopenia was associated with higher risks of osteoporosis, low radius and tibial vBMD, and falls (OR > 1.0, p < 0.05). After adjusted by age, physical activity, and calcium intake, obesity did not significantly affect bone status, while sarcopenia was still independently related to higher risks of osteoporosis, low tibial vBMD, and falls. Individuals with sarcopenic obesity had similar risks of bone impairments and fall compared to normal people. When obesity was defined by BF%, sarcopenia was associated with lower BMD in all bone sites, and sarcopenic obesity was related to higher fall risk (OR > 1.0, p < 0.05). However, obesity was not significantly associated with reduced osteoporosis risks (p > 0.05). In the adjusted model, the role of sarcopenia as a risk factor in low radial vBMD was no longer present.

TABLE 3

Table 3 The odds ratio of osteoporosis, low radial and tibial vBMD, and fall of obesity, sarcopenia, and SO compared to normal groups when obesity was defined by BMI.

TABLE 4

Table 4 The odds ratio of osteoporosis, low radial and tibial vBMD, and fall of obesity, sarcopenia, and SO compared to normal groups when obesity was defined by BF%.

3.4. Contribution of fat and muscle indicators in osteoporosis identification

We conducted the established “white box” decision tree machine learning model to detect osteoporosis based on the weight, height, age, BMI, muscle and fat indicators, medical history (hypertension, diabetes mellitus, dyslipidemia, cardiovascular diseases, fatty liver, autoimmune diseases, and malignancy), calcium intake, Baecke score, smoking, and drinking status. A four-depth model was chosen since the increased depths did not obviously improve the accuracy and may lead to overfit (Figure 3). According to the test, the most significant node (root) was body weight, followed by gait speed and limb lean mass. Baecke score, calcium intake, lower limb strength, trunk-to-limb fat ratio, hypertension status, age, and height were also selected. The decision tree had a specificity of 0.92, a sensitivity of 0.27, and an accuracy of 0.85. The prevalence osteoporosis in people with body weight lower than 57.1 kg was higher (27.3%) compared to those with higher weight (3.4%). People with both lower weight and gait speed (<0.92 m/s) had higher prevalence of osteoporosis (50.0%) compared to those with lower weight and higher gait speed (20.6%). Old people with higher weight but reduced lower limb strength (<12.5 kg) had higher prevalence of osteoporosis (16.7%) than those with higher lower limb strength (2.4%). Old people with lower body weight and gait speed had higher prevalence of osteoporosis. For those with a higher body weight and increased lower limb strength, the prevalence of osteoporosis was lowest.

FIGURE 3

Figure 3 Machine learning decision tree classification of age-related osteoporosis according to basic information, muscle and fat indicators, physical activity, calcium intake, and medical history. Testing showed that the decision tree with a depth of 4 had a specificity of 0.92, a sensitivity of 0.27, and an accuracy of 0.85 to classify osteoporosis individuals. Normal = 0 (yellow), Osteoporosis = 1 (green).

4. Discussion

Increased weight, limb lean mass, ASMI, physical performance, and lean-to-fat ratio were significantly associated with better bone health, while increased BF%, age, and limb fat mass were harmful. BMI as an indicator for obesity definition was not as relevant to bone as the above parameters. Age-related BMD reduction was more obvious in non-weight-bearing bone radius than other bone sites. The body fat mass was not as stronger related to bone parameter as muscle mass. It is reported that BF% can better predict clinical outcomes as compared to BMI in the elderly population, which is still applicable in bone turnover research (37). High BF% and limb fat mass dramatically impairs bone strength, which indicates that fat infiltration in the bone marrow may occur (38). Our results showed that trunk-to-limb fat ratio mildly benefit bone structure, since these populations usually have higher body weight and bone mechanical stimulation. However, abdominal obesity was related to increased hip fracture risk after normalized by BMI (39). The trunk fat mass measured by DXA included not only abdominal fat, but also fat in the pelvis area; therefore, our data could not show the impact of abdominal obesity on bone health. The skeletal muscle has positive roles in bones, as it regulates metabolism and provides mechanical stimulations (loading and straining) (6, 12, 40). Muscle-related indicators including functions were significantly associated with bone parameters. Lower limb strength was better correlated with bone parameters than HGS and gait speed. Total hip BMD was strongly related to limb lean mass, which shows the cause of hip fracture frequently occurring in patients with sarcopenia (16). Stiffness as a predictor of frailty fracture increased along with muscle mass and functional performance (35). BMD as the diagnostic index for osteoporosis was more positively relevant to muscle than fat mass and BMI, and negatively related with age and BF% in both weight-bearing and non-weight-bearing bones. Herein, muscle mass and function had more benefits in BMD and bone stiffness than fat mass.

Similar to the relationship between BMI and fracture risk, BMI and BMD had a non-linear relationship in hip and spine in both genders, and peripheral bones in male subjects only (36). These findings indicate that with BMI increment, the increase rate of BMD gradually flattened out. The radius BMD in male subjects started to decrease after the BMI reached 31.9 kg/m². Since the radius is a non-weight-bearing bone, excessive weight may impair bone homeostasis, but only plays minor mechanical roles (41). Different from BMI and FMI, an increase in ASMI accelerated BMD enhancement in hip and spine in male subjects and spine in female subjects, indicating that a higher ASMI is beneficial to bone homeostasis. The changing trends of FMI and BMD were more similar to BMI with the comparison of ASMI. In female subjects, the deposit of fat was not as adverse as in male subjects. It may be attributed to the sex differences of fat distribution, as pear-shaped body fat distribution in female subjects bring fewer metabolic problems than the central obesity in male subjects (42). In summary, we found from the non-linear models that with an increase in BMI, the rate of increase for BMD generally slows down, except for tibia and radius in female subjects. Fat and muscle mass indexes play different roles in the increase rate. The balance of bone mechanical stimulation and unhealthy metabolic status tilts to the latter with the BMI increment. In non-weight-bearing bones, the downtrend of BMD was found in male subjects with obesity. Excessive muscle instead of fat increment did not slow down the increase rate of BMD in hip and spine.

There is a disease classification for old people that can reflect the muscle and fat compositions (43). Obesity alone, sarcopenia alone, sarcopenic obesity, and normal groups are used to study the importance of muscle and fat on bone. BMI is a commonly used indicator to define obesity. After removing the patients with sarcopenia from the normal BMI population, individuals with BMI-defined obesity only had limited bone parameters better than the healthy control group, including total hip BMD in both genders and tibial BMD in male subjects. People with sarcopenia had lower total hip BMD, as well as radial and tibial stiffness than the normal individuals. Poor muscle status played significantly negative roles in diverse bone sites except for spine, while obese status only slightly improved bone structure in hip. Hence, people with normal BMI and good muscle status have similar bone microstructure and strength compared to people with obesity. Old people with SO had similar bone parameters compared to other groups, indicating that fat prevented negative effects from the loss of muscle mass, and poor muscle status attenuated the protective effects of fat mass. Considering body composition, those with higher BMI might have higher muscle mass instead of fat mass, but are regarded as obese. Increased body fat mass was more related to clinical outcomes in old individuals than BMI (44, 45). Therefore, BF% is recommended to identify obesity in old people. When obesity was defined by BF%, BMD was comparable between obese and normal groups, which indicated that BF%-defined obesity did not reduce the risk of osteoporosis as expected from the “obesity paradox”. Without considering muscle mass, the increased fat mass percentage in old people did not benefit BMD. Although BMI-defined obesity is related to lower risks of osteoporosis in bearing bones, the protective role did not exist after the adjustments of age, physical activity, and calcium intake. However, sarcopenia is still a dependent risk factor of osteoporosis. The lower muscle mass had stronger effects on the development of osteoporosis. BF%-defined obesity was not significantly related to reduced osteoporosis risks in both models. SO instead of sarcopenia was associated with fall risks when obesity was defined by BF%. Since non-weight-bearing bones were less related to obesity and sarcopenia, the prevention of fall might be the best approach to reduce the fracture risk of those sites. From the odds ratio of fall, the prevention or treatment of sarcopenia and SO should be performed for old people to avoid frailty fracture of non-weight-bearing bones. Patients with sarcopenia hiding in the normal BMI population were the confounding factor that created a false impression of the “obesity paradox”; therefore, BF% might be a better indicator for obesity diagnosis in old people. From the current results, the “obesity paradox” did not exist in old people with obesity defined by BF%, which was consistent with basic science studies and underlying molecular mechanisms (11). From the odds ratio results, sarcopenia is more associated with age-related osteoporosis than obesity, which emphasized the role of muscle in osteoporosis (46).

From the decision tree testing, we observed that the sensitivity of osteoporosis classification according to the involved parameters was low. The underlying reason is that this cohort only recruited those with a low risk of abnormal bone metabolism. Numerous known influencing factors of osteoporosis have already been excluded. Our study compared the role of muscle and fat in age-related bone health. Few confounding factors related to secondary osteoporosis were involved in this study. We could still extract information from the current model, since the accuracy is 85%. Body weight is the most important predictor of osteoporosis. Lower limb performances were also good indicators to predict osteoporosis. Those with greater weight and higher lower limb strength had lower risks of osteoporosis. Higher gait speed can protect people with low weight from osteoporosis compared to those with lower gait speed. Higher limb lean mass also slightly contributed to osteoporosis prevention in the prediction model. Fat mass as well as BMI were not automatically included. Muscle and bone are physically close to each other, and our findings demonstrated that they also have a functional crosstalk. It provides implications for clinicians and nutritionists that body weight, muscle function, and mass for pre-onset osteoporosis patients should be focused on, because skeletal muscle not only provides more mechanical stimulation, but also beneficially affects both bone and systemic metabolism.

There are several strengths to our study. To our knowledge, this is the first research to study the muscle and fat effects on bone health together and compare their contributions (Figure 4). We indicated that obesity did not prevent osteoporosis when compared to people with normal BMI or BF% and without sarcopenia. The bone sites with or without weight bearing had different reflections on BMI, muscle, and fat alterations, as well as gender. Our study also had limitations. Our cohort excluded those with disease-induced osteoporosis, and thus, the results might not be necessarily universal for an old population who had secondary osteoporosis. Participants with missing data have also been excluded in several analyses. This is a cross-sectional study that only reveals the association between muscle/fat and bone, but not the causal relationship. Prospective longitudinal or interventional studies are warranted in the future.

FIGURE 4

Figure 4 Flowchart of the study design. A combined analysis was performed to compare the role of muscle and fat in bone homeostasis. Correlation heatmap showed that the increased muscle mass-related indexes had a stronger positive correlation with bone structure compared to fat mass-related indexes. Non-linear regression showed that the increase rate of BMD in hip and spine almost did not slow down with increased ASMI but FMI. The odds ratio with 95% CI indicated that sarcopenia induced high incidence of osteoporosis, while obesity did not effectively prevent it compared to the normal group. The decision tree also showed that greater weight, muscle performance, and limb lean mass can be negative predictors of osteoporosis. The red and purple up arrows indicated more remarkable roles of increased muscle mass in weight-bearing bones’ (hip, spine, and tibia) health than fat mass, while radius was not significantly affected by them.

In conclusion, skeletal muscle mass was positively related to bone structures, while body fat percentage and age were negatively associated with bone health. BMD was non-linearly associated with increased BMI, of which the increase rate generally slowed down in bearing bones. Increased ASMI did not slow down BMD increase rate in lumbar spine and total hip as BMI and FMI, especially in male subjects. BMI-defined obesity failed to significantly reduce the risk of osteoporosis in the adjusted model when the reference group involved people with normal BMI and without sarcopenia. BF%-defined obesity was not associated with improved bone health. Sarcopenia was significantly related to increased risk of osteoporosis. Old people with higher body weight, better lower extremity performance, and larger limb lean mass had lower prevalence of osteoporosis. The improvement of muscle status should be regarded as a clinical recommendation for elderly people to prevent age-related osteoporosis.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Materials, requests to access these datasets should be directed to Ling Qin, bGluZ3FpbkBjdWhrLmVkdS5oaw==.

Ethics statement

The studies involving human participants were reviewed and approved by The Joint Chinese University of Hong Kong – New Territories East Cluster Clinical Research Ethics Committee (Affiliation: The Chinese University of Hong Kong and New Territories East Cluster). The patients/participants provided their written informed consent to participate in this study.

Author contributions

CL: writing—original draft and editing, conceptualization, and methodology. PW: methodology and writing—review and editing. XT: statistical analysis and data visualization. SC: supervision and writing—review and editing. VH: investigation and data curation. W-HC: supervision and writing—review and editing. QL: investigation, supervision, and writing—review and editing. SL: conceptualization and validation. RW: conceptualization, investigation, supervision, and writing—review and editing. All authors contributed to the article and approved the submitted version.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.990442/full#supplementary-material

References

1. Coughlan T, Dockery F. Osteoporosis and fracture risk in older people. Clin Med (Lond) (2014) 14(2):187–91. doi: 10.7861/clinmedicine.14-2-187

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Wong RMY, Wong PY, Liu C, Wong HW, Chung YL, Chow SKH, et al. The imminent risk of a fracture-existing worldwide data: A systematic review and meta-analysis. Osteoporos Int (2022) 33(12):2453–66. doi: 10.1007/s00198-022-06473-0

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Friedman SM, Mendelson DA. Epidemiology of fragility fractures. Clin Geriatr Med (2014) 30(2):175–81. doi: 10.1016/j.cger.2014.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Borgstrom F, Karlsson L, Ortsater G, Norton N, Halbout P, Cooper C, et al. Fragility fractures in Europe: Burden, management and opportunities. Arch Osteoporos (2020) 15(1):59. doi: 10.1007/s11657-020-0706-y

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Setareh A, Williams P, Chastek MS, Sundquist MS, Barrera-Sierra MD, Leader JPMBA, et al. Economic burden of osteoporotic fractures in us managed care enrollees. Am J Managed Care (2020) 26(05):e142–e9. doi: 10.37765/ajmc.2020.43156

CrossRef Full Text | Google Scholar

6. Fassio A, Idolazzi L, Rossini M, Gatti D, Adami G, Giollo A, et al. The obesity paradox and osteoporosis. Eat Weight Disord (2018) 23(3):293–302. doi: 10.1007/s40519-018-0505-2

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Tang H, He J-H, Gu H-B, Zhu K, Lu C-J, Sun L-L, et al. The different correlations between obesity and osteoporosis after adjustment of static mechanical loading from weight and fat free mass. J Musculoskelet Neuronal Interact (2021) 21(3):351–7.

PubMed Abstract | Google Scholar

8. Nielson CM, Srikanth P, Orwoll ES. Obesity and fracture in men and women: An epidemiologic perspective. J Bone Miner Res (2012) 27(1):1–10. doi: 10.1002/jbmr.1486

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Modig K, Erdefelt A, Mellner C, Cederholm T, Talbäck M, Hedström M. "Obesity paradox" holds true for patients with hip fracture: A registry-based cohort study. J Bone Joint Surg Am (2019) 101(10):888–95. doi: 10.2106/jbjs.18.01249

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Jain RK, Vokes T. Fat mass has negative effects on bone, especially in men: A cross-sectional analysis of nhanes 2011-2018. J Clin Endocrinol Metab (2022) 107(6):e2545–52. doi: 10.1210/clinem/dgac040

CrossRef Full Text | Google Scholar

11. Dai B, Xu J, Li X, Huang L, Hopkins C, Wang H, et al. Macrophages in epididymal adipose tissue secrete osteopontin to regulate bone homeostasis. Nat Commun (2022) 13(1):427. doi: 10.1038/s41467-021-27683-w

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Pagnotti GM, Styner M, Uzer G, Patel VS, Wright LE, Ness KK, et al. Combating osteoporosis and obesity with exercise: Leveraging cell mechanosensitivity. Nat Rev Endocrinol (2019) 15(6):339–55. doi: 10.1038/s41574-019-0170-1

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Kaji H. Effects of myokines on bone. Bonekey Rep (2016) 5:826. doi: 10.1038/bonekey.2016.48

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Clynes MA, Gregson CL, Bruyere O, Cooper C, Dennison EM. Osteosarcopenia: Where osteoporosis and sarcopenia collide. Rheumatol (Oxford) (2021) 60(2):529–37. doi: 10.1093/rheumatology/keaa755

CrossRef Full Text | Google Scholar

15. Wong RMY, Wong H, Zhang N, Chow SKH, Chau WW, Wang J, et al. The relationship between sarcopenia and fragility fracture-a systematic review. Osteoporos Int (2019) 30(3):541–53. doi: 10.1007/s00198-018-04828-0

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Hong W, Cheng Q, Zhu X, Zhu H, Li H, Zhang X, et al. Prevalence of sarcopenia and its relationship with sites of fragility fractures in elderly Chinese men and women. PloS One (2015) 10(9):e0138102. doi: 10.1371/journal.pone.0138102

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Prado CM, Siervo M, Mire E, Heymsfield SB, Stephan BC, Broyles S, et al. A population-based approach to define body-composition phenotypes. Am J Clin Nutr (2014) 99(6):1369–77. doi: 10.3945/ajcn.113.078576

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Zhu TY, Yip BH, Hung VW, Choy CW, Cheng KL, Kwok TC, et al. Normative standards for hrpqct parameters in Chinese men and women. J Bone Miner Res (2018) 33(10):1889–99. doi: 10.1002/jbmr.3481

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Huang M, Hung VW, Li TK, Law SW, Wang Y, Chen S, et al. Performance of hr-pqct, dxa, and frax in the discrimination of asymptomatic vertebral fracture in postmenopausal Chinese women. Arch Osteoporos (2021) 16(1):125. doi: 10.1007/s11657-021-00939-0

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Hung VW, Zhu TY, Cheung WH, Fong TN, Yu FW, Hung LK, et al. Age-related differences in volumetric bone mineral density, microarchitecture, and bone strength of distal radius and tibia in Chinese women: A high-resolution pqct reference database study. Osteoporos Int (2015) 26(6):1691–703. doi: 10.1007/s00198-015-3045-x

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Pauchard Y, Liphardt AM, Macdonald HM, Hanley DA, Boyd SK. Quality control for bone quality parameters affected by subject motion in high-resolution peripheral quantitative computed tomography. Bone (2012) 50(6):1304–10. doi: 10.1016/j.bone.2012.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Burghardt AJ, Buie HR, Laib A, Majumdar S, Boyd SK. Reproducibility of direct quantitative measures of cortical bone microarchitecture of the distal radius and tibia by hr-pqct. Bone (2010) 47(3):519–28. doi: 10.1016/j.bone.2010.05.034

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Nishiyama KK, Shane E. Clinical imaging of bone microarchitecture with hr-pqct. Curr Osteoporos Rep (2013) 11(2):147–55. doi: 10.1007/s11914-013-0142-7

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Liu XS, Zhang XH, Sekhon KK, Adams MF, McMahon DJ, Bilezikian JP, et al. High-resolution peripheral quantitative computed tomography can assess microstructural and mechanical properties of human distal tibial bone. J Bone Miner Res (2010) 25(4):746–56. doi: 10.1359/jbmr.090822

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Roberts HC, Denison HJ, Martin HJ, Patel HP, Syddall H, Cooper C, et al. A review of the measurement of grip strength in clinical and epidemiological studies: Towards a standardised approach. Age Ageing (2011) 40(4):423–9. doi: 10.1093/ageing/afr051

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Scott D, Chandrasekara SD, Laslett LL, Cicuttini F, Ebeling PR, Jones G. Associations of sarcopenic obesity and dynapenic obesity with bone mineral density and incident fractures over 5-10 years in community-dwelling older adults. Calcif Tissue Int (2016) 99(1):30–42. doi: 10.1007/s00223-016-0123-9

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Kanis JA, Melton L, Christiansen C, Johnston CC, Khaltaev N. The diagnosis of osteoporosis. J Bone Miner Res (1994) 9(8):1137–41. doi: 10.1002/jbmr.5650090802

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Chen LK, Woo J, Assantachai P, Auyeung TW, Chou MY, Iijima K, et al. Asian Working group for sarcopenia: 2019 consensus update on sarcopenia diagnosis and treatment. J Am Med Dir Assoc (2020) 21(3):300–7.e2. doi: 10.1016/j.jamda.2019.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Donini LM, Busetto L, Bischoff SC, Cederholm T, Ballesteros-Pomar MD, Batsis JA, et al. Definition and diagnostic criteria for sarcopenic obesity: Espen and easo consensus statement. Obes Facts (2022) 15(3):321–35. doi: 10.1159/000521241

PubMed Abstract | CrossRef Full Text | Google Scholar

30. World Health Organization. Regional office for the Western p. In: The Asia-pacific perspective: Redefining obesity and its treatment. Sydney: Health Communications Australia (2000).

Google Scholar

31. Marques LP, Confortin SC, Ono LM, Barbosa AR, d'Orsi E. Quality of life associated with handgrip strength and sarcopenia: Epifloripa aging study. Arch Gerontol Geriatr (2019) 81:234–9. doi: 10.1016/j.archger.2018.12.015

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Woo J, Leung SSF, Ho SC, Lam TH, Janus ED. A food frequency questionnaire for use in the Chinese population in Hong Kong: Description and examination of validity. Nutr Res (1997) 17(11):1633–41. doi: 10.1016/S0271-5317(97)00170-X

CrossRef Full Text | Google Scholar

33. Ho SC, Yu R, Chan S. Comparison of the modified Chinese baecke questionnaire with a 3-day activity diary in a Hong Kong Chinese population. Asia Pac J Public Health (2015) 27(2):NP2358–71. doi: 10.1177/1010539511416805

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Parikh SA, Gomez R, Thirugnanasambandam M, Chauhan SS, De Oliveira V, Muluk SC, et al. Decision tree based classification of abdominal aortic aneurysms using geometry quantification measures. Ann BioMed Eng (2018) 46(12):2135–47. doi: 10.1007/s10439-018-02116-w

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Cheung WH, Hung VW, Cheuk KY, Chau WW, Tsoi KK, Wong RM, et al. Best performance parameters of hr-pqct to predict fragility fracture: Systematic review and meta-analysis. J Bone Miner Res (2021) 36(12):2381–98. doi: 10.1002/jbmr.4449

PubMed Abstract | CrossRef Full Text | Google Scholar

36. De Laet C, Kanis JA, Oden A, Johanson H, Johnell O, Delmas P, et al. Body mass index as a predictor of fracture risk: A meta-analysis. Osteoporos Int (2005) 16(11):1330–8. doi: 10.1007/s00198-005-1863-y

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Chang SH, Beason TS, Hunleth JM, Colditz GA. A systematic review of body fat distribution and mortality in older people. Maturitas (2012) 72(3):175–91. doi: 10.1016/j.maturitas.2012.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Rosen CJ, Bouxsein ML. Mechanisms of disease: Is osteoporosis the obesity of bone? Nat Clin Pract Rheumatol (2006) 2(1):35–43. doi: 10.1038/ncprheum0070

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Sogaard AJ, Holvik K, Omsland TK, Tell GS, Dahl C, Schei B, et al. Abdominal obesity increases the risk of hip fracture. A population-based study of 43,000 women and men aged 60-79 years followed for 8 years. cohort of Norway. J Intern Med (2015) 277(3):306–17. doi: 10.1111/joim.12230

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Pimenta LD, Massini DA, Santos Dd, Vasconcelos CMT, Simionato AR, Gomes LAT, et al. Bone health, muscle strength and lean mass: Relationships and exercise recommendations. Rev Bras Medicina do Esporte (2019) 25(3):245–51. doi: 10.1590/1517-869220192503210258

CrossRef Full Text | Google Scholar

41. Benova A, Tencerova M. Obesity-induced changes in bone marrow homeostasis. Front Endocrinol (Lausanne) (2020) 11:294. doi: 10.3389/fendo.2020.00294

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Karastergiou K, Smith SR, Greenberg AS, Fried SK. Sex differences in human adipose tissues - the biology of pear shape. Biol Sex Differ (2012) 3(1):13. doi: 10.1186/2042-6410-3-13

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Batsis JA, Villareal DT. Sarcopenic obesity in older adults: Aetiology, epidemiology and treatment strategies. Nat Rev Endocrinol (2018) 14(9):513–37. doi: 10.1038/s41574-018-0062-9

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Heitmann BL, Erikson H, Ellsinger BM, Mikkelsen KL, Larsson B. Mortality associated with body fat, fat-free mass and body mass index among 60-Year-Old Swedish men-a 22-year follow-up. the study of men born in 1913. Int J Obes Relat Metab Disord (2000) 24(1):33–7. doi: 10.1038/sj.ijo.0801082

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Javed AA, Aljied R, Allison DJ, Anderson LN, Ma J, Raina P. Body mass index and all-cause mortality in older adults: A scoping review of observational studies. Obes Rev (2020) 21(8):e13035. doi: 10.1111/obr.13035

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Santos VRD, Christofaro DGD, Gomes IC, Junior IFF, Gobbo LA. Relationship between obesity, sarcopenia, sarcopenic obesity, and bone mineral density in elderly subjects aged 80 years and over. Rev Bras Ortop (2018) 53(3):300–5. doi: 10.1016/j.rboe.2017.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar