Genetic and environmental correlations between bone geometric parameters and body compositions

What are the differences in paraspinal muscle morphometry among degenerative spondylolisthesis patients, isthmic spondylolisthesis patients, and healthy individuals? A propensity score matching analysis

PURPOSE

To compare the morphometry of paraspinal muscles in patients with degenerative spondylolisthesis (DS), isthmic spondylolisthesis (IS), and healthy individuals.

METHODS

Thirty-seven pairs of DS patients were selected using propensity score matching with IS patients, while 37 healthy individuals matched for age, sex, and BMI were selected as controls. The relative cross-sectional area (rCSA), and relative functional cross-sectional area (rfCSA) of paraspinal muscles were measured, and the degree of fatty infiltration (FI) was calculated. Based on occupational differences, the patients were also divided into worker and farmer groups, and the same measurements were taken on them.

RESULTS

At the L3/L4 level, the multifidus (MF) FI was greater in the DS and IS groups than in the control group, the erector spinae (ES) rfCSA was higher in the IS group than in the DS and control groups. At the L4/L5 level, MF rfCSA was smaller in the DS and IS groups than in the control group; ES rfCSA was higher in the IS group than in the DS and control groups. At the L5/S1 level, MF rfCSA was smaller in the DS and IS groups than in the control group; ES rfCSA was higher in the IS group than in the DS group. At the L3/L4, L4/L5 level, MF rfCSA were higher in the worker group than in the farmer group (p < 0.05).

CONCLUSION

The morphological changes in paraspinal muscles in patients with DS were dominated by selective atrophy of the MF, while in patients with IS, the morphological changes in paraspinal muscle showed selective atrophy of the MF accompanied by compensatory hypertrophy of the ES. The surgeon should consider the morphological differences in paraspinal muscle between different types of lumbar spondylolisthesis when establishing the appropriate surgical program.

Does Unilateral Lumbosacral Radiculopathy Affect the Association between Lumbar Spinal Muscle Morphometry and Bone Mineral Density?

Age-related degenerative changes lead to a gradual decrease in bone mineral density (BMD) and muscle mass. We aimed to assess the effects of decreased BMD and lumbar denervation on lumbar spinal muscle morphometry and the relationship between BMD and lumbar spinal muscular morphometry, respectively. Eighty-one patients, aged 50–85 years, diagnosed with unilateral lumbosacral radiculopathy based on electrodiagnostic studies between January 2016 and April 2021 were enrolled. BMD T scores in the lumbar spine and hip were measured using dual-energy X-ray absorptiometry. The cross-sectional area (CSA) of the psoas, multifidus, and erector spinae located in the middle of the lumbar spine, between the L3 and L4 and between the L4 and L5 levels, respectively, was measured using axial MRI. Functional CSA (FCSA) was defined as the CSA of lean muscle mass. Pearson correlation analyses were performed to evaluate the association between BMD T scores and the CSA, FCSA, and the ratio of the FCSA to the CSA (functional ratio) for each side. The CSA of lumbar spinal muscles showed no significant correlation with lumbar BMD. The FCSA and functional ratio of lumbar spinal muscles were significantly correlated with lumbar BMD. There was no correlation between femur BMD and lumbar spinal muscle morphometry.

Interactive effect of sarcopenia and falls on vertebral osteoporotic fracture in patients with rheumatoid arthritis

Patients with rheumatoid arthritis (RA) had higher incidences of sarcopenia, falls, osteoporosis, and vertebral osteoporotic fractures (VOPF). Sarcopenia was associated with longer disease duration, higher disease activity, and more severe RA. The interactive effect of sarcopenia and falls was associated with a higher risk of VOPF in patients with RA.

PURPOSE

Whether sarcopenia and falls are a risk factor for vertebral fracture in RA patients has not been demonstrated. This study aimed to explore the incidence of vertebral osteoporotic fracture (VOPF) and its relationship with sarcopenia and falls in RA patients.

METHODS

A total of 474 RA patients and 156 controls were enrolled in this study. Anteroposterior and lateral X-ray examinations of the vertebral column (T4-L4) were used for the semiquantitative assessment of VOPF. Bone mineral density was measured by dual-energy X-ray absorptiometry. Skeletal muscle mass was measured by direct segmental multifrequency bioelectrical impedance analysis (DSM-BIA method).

RESULTS

RA patients had an increased risk of sarcopenia (62.4% vs 9.0%, x² = 47.478, P < 0.001), falls (30.2% vs 3.2%), osteoporosis (OP) (33.5% vs 12.8%, x² = 134.276, P < 0.001), and VOPF (20.3% vs 3.8%, x² = 47.478, P < 0.001) than controls. Patients with sarcopenia were more likely to have VOPF than RA without sarcopenia (24.0% vs 14.0%, x² = 6.802, P = 0.009). RA with sarcopenia and prior falls had the highest incidences of VOPF (36.7%). Older age (OR = 1.056, P < 0.001, 95% CI 1.030-1.083), falls (OR = 2.043, P = 0.003, 95% CI 1.238-3.371), OP (OR = 1.819, P = 0.034, 95% CI 1.046-3.163), and usage of glucocorticoids (GCs) (OR = 1.862, P = 0.022, 95% CI 1.093-3.172) were risk factors for VOPF in RA patients, while a higher skeletal muscle index (SMI) was a protective factor (OR = 0.754, P = 0.038, 95% CI 0.578-0.984) for VOPF in RA patients.

CONCLUSIONS

The interactive effect of sarcopenia and falls is associated with a higher risk of VOPF in patients with RA.

Newly developed hip geometry parameters are associated with hip fracture

BACKGROUND

The conventional and newly developed geometry parameters at the femoral neck have formed a large geometry profile and their relationship with hip fracture was largely unknown. We aimed to evaluate the relationship between the geometry profile and hip fracture in Chinese.

METHODS

The hip geometry profile contains seven conventional geometry parameters at femoral neck (FN) and thirty newly developed parameters at three sub-regions (Narrow Neck, NN; Intertrochanter, IT; Femoral shaft, FS) of the total hip. Based on 6294 recruited Chinese (≥65 years), 97 subjects with osteoporotic fracture (OF) history and 388 matched controls were selected. The t test, Chi-square statistic, conditional logistic regression model were used.

RESULTS

Three geometric parameters (endocortical diameter, ED; cortical thickness, CT; buckling ratio, BR) have consistent differences at all the sites between the cases and controls (p < 0.01). Conventional geometry parameters (e.g., cross-sectional area, CSA; BR) and the newly developed parameters (e.g., NN_ED, NN_Outer Diameter, IT_ED) were identified as the risk factors of hip fracture independent of BMD. The additional predictive ability of the hip geometric parameters, over BMD alone, (receiver operating characteristic curve (ROC) analysis) was improved. Especially at NN, the area of ROC used single NN_BMD is only 0.596, but increased rapidly at 0.705 when combined with the hip geometric parameters.

CONCLUSIONS

This study found that three newly developed hip geometry parameters are associated with hip fracture. The results will increase our understanding of the determinants of fracture and provide potential clue for future prevention of fracture in Chinese Population.

Effect of race, age, and gender on lumbar muscle volume and fat infiltration in the degenerative spine

BACKGROUND

The quantity and quality of spinal muscles in patients with degenerative spinal diseases and various backgrounds such as age, gender, or race is unclear. We quantitatively evaluated the cross-sectional area (CSA) and fatty degeneration of the muscles around the spine, using magnetic resonance imaging (MRI) in patients with degenerative spinal disease, and studied the effects of age, gender, and race in multicenter retrospective study.

METHODS

The subjects were Caucasian and Asian patients with degenerative lumbar disease who underwent L4-5 single level spinal fusion surgery at centers in the United States and Japan. Using preoperative axial T2 MRI at the L4-5 disc level, the cross-sectional areas of the psoas and paraspinal muscles were measured. Fat infiltration was measured using the threshold method, and percent fat area (%FA) was calculated for each muscle. The muscle/disc area ratio (MDAR) was used to control for size differences per patient. T-test, Pearson's correlation coefficient, partial correlation, and multiple linear regression were used for statistical analysis.

RESULTS

In total, 140 patients (53 men; 87 women; mean age, 69.2 years) were analyzed. Age was similar in Caucasians (n = 64) and Asians (n = 76). MDARs were larger in Caucasians for paraspinal and psoas muscles (p < 0.005). Percent FA of psoas was similar in Caucasians and Asians, but greater in the paraspinal muscles of Asians (p < 0.05). After controlling for race and gender, age was correlated negatively with MDAR (p < 0.001) and positively with %FA (p < 0.001). In the multiple linear regression analysis, age, gender, and race were independently affected by MDAR and %FA.

CONCLUSIONS

Lumbar muscle mass and quality were affected by age, gender, and race, independently, in patients with degenerative lumbar disease.

Multifidus muscle fatty infiltration as an index of dysfunction in patients with single-segment degenerative lumbar spinal stenosis: A case-control study based on propensity score matching

The multifidus muscle morphology and its relation to the function of patients with degenerative lumbar spinal stenosis (DLSS) remains unclear. This study aimed to investigate the multifidus muscle morphology in patients with DLSS and to determine its relations to the patients function. Sixty-two patients with single-segment DLSS at L4-5 and sixty control patients with non-spinal-derived low back pain were retrospectively enrolled and further matched based on propensity scores. The Oswestry Disability Index (ODI) and bodily pain using the Short-Form Health Survey were evaluated. The cross-sectional area (CSA), CSA of fatty free (CSAF), and fatty infiltration rate [FIR; i.e., (1- CSAF/CSA) × 100%] of the multifidus muscle were measured on magnetic resonance images using ImageJ software. Adjustment for confounders was performed using generalized linear models. The FIR at L5-S1 in controls was statistically significant but slightly less than the DLSS group. The between-groups difference was 5% (p < 0.001), and 2.8% (p = 0.036) in the complete and matching cohorts, respectively, after adjustment. Statistically significant differences were not observed in other multifidus muscle parameters between the groups. FIR > 20% at L5-S1 was independently associated with ODI ≥ 41 in patients with DLSS [Retaining demography as control block or not, Odds ratio (OR) = 8.4, p = 0.023; OR = 12.3, p = 0.030]. The multifidus muscle at L5-S1 demonstrated slightly greater fatty infiltration in patients with L4-5 single-segment DLSS than controls. Significant fatty infiltration in the multifidus muscle at L5-S1 may be correlated with poor function in patients with L4-5 single-segment DLSS.

Lumbar lordosis reduction and disc bulge may correlate with multifidus muscle fatty infiltration in patients with single-segment degenerative lumbar spinal stenosis

OBJECTIVES

To investigate the correlation between fatty infiltration in the multifidus muscle related to the involved nerve root and structural parameters associated with stenosis in patients with degenerative lumbar spinal stenosis (DLSS).

PATIENTS AND METHODS

Sixty-four patients with single-segment DLSS at L4-5 were retrospectively enrolled. The fatty infiltration rate (FIR) of the multifidus muscle at L5-S1, lumbar lordosis and the cross-sectional area (CSA) of the structural parameters at L4-5, such as dural sac, disc bulge, ligamentum flava and vertebral body of L5 were measured on magnetic resonance images using ImageJ software. All enrolled patients were divided into an FIR < 25 % group and an FIR ≥ 25 % group according to the FIR of the multifidus muscle at L5-S1. The propensity scores matching and adjustment of potential covariates were performed to reduce the confounding bias between the two groups.

RESULTS

Lumbar lordosis in the FIR ≥ 25 % group was significantly lower than that in the FIR＜25 % group in both cohorts. The mean differences in lumbar lordosis of 14.16 degrees between the two groups in the complete cohort and of 14.23 degrees in the matched cohort remained significant after adjustment. The disc bulge CSA/ vertebral body CSA in the FIR ≥ 25 % group was greater than that in the FIR＜25 % group in both cohorts. The mean differences in the disc bulge CSA/ vertebral body CSA between the two groups of 0.67 in the complete cohort and 0.96 in the matched cohort were statistically significant after adjustment. There was no significant difference in the dural sac CSA/ vertebral body CSA and ligamentum flava CSA/ vertebral body CSA between the two groups in either cohort regardless of adjustment. Logistic regression analysis for FIR ≥ 25 % in the multifidus muscle at L5-S1 exhibited that the disc bulge CSA/ vertebral body CSA were independent risk factors with odds ratio (OR) of 8.52, while lumbar lordosis were independent protective factors (OR = 0.72).

CONCLUSIONS

The disc bulge at the stenosis segment and lumbar lordosis reduction may be correlated with fatty infiltration in the multifidus muscles at L5-S1 in patients with L4-5 single-segment DLSS.

Association of MRI-defined lumbar paraspinal muscle mass and slip percentage in degenerative and isthmic spondylolisthesis: A multicenter, retrospective, observational study

The objective of this study is to investigate the role of paraspinal muscles in the progression of different types of spondylolisthesis by examining the correlation between cross-sectional area (CSA) of lumbar paraspinal muscle and slip percentage (SP) in degenerative spondylolisthesis and isthmic spondylolisthesis.A multicenter retrospective analysis was carried out including 219 subjects diagnosed with lumbar spondylolisthesis. Using T2-weighted axial magnetic resonance imgaging, CSAs of the psoas major (PM), multifidus (MU), and erector spinae were measured and divided by L5 vertebral body (VB) CSA. SP was measured using sagittal T2-weighted images. Correlations between muscle CSA ratio and SP were calculated in each group. Regression analysis was performed to predict the influence of each muscle CSA/VB CSA ratio on SP.No significant correlation was found in the degenerative spondylolisthesis group between any of the muscle CSA ratios and SP. Both PM/VB ratio (r = -0.24, P = .021) and MU/VB ratio (r = -0.26, P = .012) were negatively correlated with SP in the isthmic spondylolisthesis group. MU had more influence on SP than PM in the isthmic spondylolisthesis group (regression coefficient MU/VB: -8.08, PM/VB: -4.34).Both PM and MU muscle CSA ratios were negatively correlated with SP in the isthmic group. MU had more influence on SP than PM. No muscles had any correlations with SP in the degenerative group. This discrepancy between the two groups suggests that exercise programs or interventions regarding the segmental stability of isthmic spondylolisthesis and degenerative spondylolisthesis should be distinguished in clinical practice.Clinical Research Information Service of Korea Centers for Disease control and Prevention, KCT0002588. Registered on 12 December 2017, https://cris.nih.go.kr/cris/search/search_result_st01.jsp?seq=10702.

Fat Infiltration in the Multifidus Muscle as a Predictor of Prognosis After Decompression and Fusion in Patients with Single-Segment Degenerative Lumbar Spinal Stenosis: An Ambispective Cohort Study Based on Propensity Score Matching

OBJECTIVE

To determine whether fat infiltration in the multifidus muscle would predict surgical prognosis in patients with degenerative lumbar spinal stenosis (DLSS).

METHODS

This ambispective cohort study enrolled 118 consecutive patients undergoing surgery for L4-5 single-segment DLSS. Fat infiltration rate (FIR) on magnetic resonance images of the multifidus muscle at L5-S1 were measured using ImageJ software. The enrolled patients were divided into FIR <25% and FIR ≥25% groups according to their FIR of the multifidus muscle at L5-S1. The 2 groups of patients who finished follow-up were further matched for the baseline covariates based on propensity scores. Patients' reported outcomes including the visual analog scale score for back pain and leg pain, and the Oswestry Disability Index (ODI) score were compared between groups at follow-up and further adjusted using generalized linear models.

RESULTS

Patients in the FIR <25% group showed statistically significantly greater reduction in ODI at 6 and 18 months after surgery than did patients in the FIR ≥25% group in either cohort regardless of adjustment; however, the 2-point between-group difference was smaller than the predefined minimum clinically important difference. In addition, more patients in the FIR <25% group achieved clinically significant improvement in ODI than those in the FIR ≥25% group in either complete cohort or matching cohort (63.8% vs. 21.1%, P < 0.001; 70.3% vs. 24.1%, P < 0.001, respectively) before and after adjustment (63.3% vs. 27.8%, P < 0.001; 69.1% vs. 31.0%, P < 0.001, respectively).

CONCLUSIONS

Fat infiltration in multifidus muscle at L5-S1 could be a potential predictor of functional improvement after surgery in patients with L4-5 single-segment DLSS.

The genetics of bone mass and susceptibility to bone diseases

Osteoporosis is characterized by low bone mass and an increased risk of fracture. Genetic factors, environmental factors and gene-environment interactions all contribute to a person's lifetime risk of developing an osteoporotic fracture. This Review summarizes key advances in understanding of the genetics of bone traits and their role in osteoporosis. Candidate-gene approaches dominated this field 20 years ago, but clinical and preclinical genetic studies published in the past 5 years generally utilize more-sophisticated and better-powered genome-wide association studies (GWAS). High-throughput DNA sequencing, large genomic databases and improved methods of data analysis have greatly accelerated the gene-discovery process. Linkage analyses of single-gene traits that segregate in families with extreme phenotypes have led to the elucidation of critical pathways controlling bone mass. For example, components of the Wnt-β-catenin signalling pathway have been validated (in both GWAS and functional studies) as contributing to various bone phenotypes. These notable advances in gene discovery suggest that the next decade will witness cataloguing of the hundreds of genes that influence bone mass and osteoporosis, which in turn will provide a roadmap for the development of new drugs that target diseases of low bone mass, including osteoporosis.

Lumbar paraspinal muscle morphometry and its correlations with demographic and radiological factors in adult isthmic spondylolisthesis: a retrospective review of 120 surgically managed cases

OBJECTIVE The objective of this study was to assess the cross-sectional areas (CSAs) of lumbar paraspinal muscles in adults with isthmic spondylolisthesis (IS), to compare them with those in the normative population, and to evaluate their correlations with demographic factors and MRI changes in various spinal elements. METHODS The authors conducted a retrospective study of patients who had undergone posterior lumbar interbody fusion for IS, and 2 of the authors acting as independent observers calculated the CSAs of various lumbar paraspinal muscles (psoas, erector spinae [ES], multifidus [MF]) on preoperative axial T2-weighted MR images from the L-3 to L-5 vertebral levels and computed the CSAs as ratios with respect to the corresponding vertebral body areas. These values were then compared with those in an age- and sex-matched normative population and were analyzed with respect to age, sex, duration of symptoms, grade of listhesis, and various MRI changes at the level of the listhesis (pedicle signal change, disc degeneration, and facetal arthropathy). RESULTS Compared with values in normative controls, the mean CSA value for the ES muscle was significantly higher in the study cohort of 120 patients (p = 0.002), whereas that for the MF muscle was significantly lower (p = 0.009), and more so in the patients with PSC (p = 0.002). Magnetic resonance imaging signal change in the pedicle was seen in half of the patients, all of whom demonstrated a Type 2 change. Of the variables tested in a multivariate analysis, age independently predicted lower area values for all 3 muscles (p ≤ 0.001), whereas female sex predicted a lower mean psoas area value (p < 0.001). None of the other variables significantly predicted any of the muscle area values. A decrease in the mean MF muscle area value alone was associated with a significantly increased likelihood of a PSC (p = 0.039). CONCLUSIONS Compared with normative controls, patients with IS suffer selective atrophy of their MF muscle, whereas their ES muscle undergoes a compensatory hypertrophy. Advancing age has a detrimental effect on the areas of the lumbar PSMs, whereas female sex predisposes to a decreased psoas muscle area. Multifidus muscle atrophy correlates with PSC, indicating the role of this deep stabilizer in the biomechanical stability of spondylolisthetic spines. This may be of clinical significance in targeted physiotherapy programs during the conservative management of IS.

The relationship between lower limb bone and muscle in military recruits, response to physical training, and influence of smoking status

The relationship between bone and skeletal muscle mass may be affected by physical training. No studies have prospectively examined the bone and skeletal muscle responses to a short controlled exercise-training programme. We hypothesised that a short exercise-training period would affect muscle and bone mass together. Methods: Femoral bone and Rectus femoris Volumes (RF_VOL) were determined by magnetic resonance imaging in 215 healthy army recruits, and bone mineral density (BMD) by Dual X-Ray Absorptiometry (DXA) and repeated after 12 weeks of regulated physical training. Results: Pre-training, RF_VOL was smaller in smokers than non-smokers (100.9 ± 20.2 vs. 108.7 ± 24.5, p = 0.018; 96.2 ± 16.9 vs. 104.8 ± 21.3, p = 0.002 for dominant/non-dominant limbs), although increases in RF_VOL with training (of 14.2 ± 14.5% and 13.2 ± 15.6%] respectively, p < 0.001) were independent of prior smoking status. Pre-training RF_VOL was related to bone cortical volume (r² = 0.21 and 0.30, p < 0.001 for dominant and non-dominant legs), and specifically to periosteal (r² = 0.21 and 0.23, p < 0.001) volume. Pre-training dominant RF_VOL was independently associated with Total Hip BMD (p < 0.001). Training-related increases in RF_VOL and bone volumes were related. Whilst smokers demonstrated lower muscle mass than non-smokers, differences were abolished with training. Training-related increases in muscle mass were related to increases in periosteal bone volume in both dominant and non-dominant legs.

Genetic and environmental variances of bone microarchitecture and bone remodeling markers: a twin study

All genetic and environmental factors contributing to differences in bone structure between individuals mediate their effects through the final common cellular pathway of bone modeling and remodeling. We hypothesized that genetic factors account for most of the population variance of cortical and trabecular microstructure, in particular intracortical porosity and medullary size - void volumes (porosity), which establish the internal bone surface areas or interfaces upon which modeling and remodeling deposit or remove bone to configure bone microarchitecture. Microarchitecture of the distal tibia and distal radius and remodeling markers were measured for 95 monozygotic (MZ) and 66 dizygotic (DZ) white female twin pairs aged 40 to 61 years. Images obtained using high-resolution peripheral quantitative computed tomography were analyzed using StrAx1.0, a nonthreshold-based software that quantifies cortical matrix and porosity. Genetic and environmental components of variance were estimated under the assumptions of the classic twin model. The data were consistent with the proportion of variance accounted for by genetic factors being: 72% to 81% (standard errors ∼18%) for the distal tibial total, cortical, and medullary cross-sectional area (CSA); 67% and 61% for total cortical porosity, before and after adjusting for total CSA, respectively; 51% for trabecular volumetric bone mineral density (vBMD; all p < 0.001). For the corresponding distal radius traits, genetic factors accounted for 47% to 68% of the variance (all p ≤ 0.001). Cross-twin cross-trait correlations between tibial cortical porosity and medullary CSA were higher for MZ (rMZ  = 0.49) than DZ (rDZ  = 0.27) pairs before (p = 0.024), but not after (p = 0.258), adjusting for total CSA. For the remodeling markers, the data were consistent with genetic factors accounting for 55% to 62% of the variance. We infer that middle-aged women differ in their bone microarchitecture and remodeling markers more because of differences in their genetic factors than differences in their environment.

Longitudinal relationships of metabolic syndrome and obesity with kidney function: Healthy Twin Study

BACKGROUND

To examine cross-sectional and longitudinal relationships, including genetic and environmental correlations, of metabolic syndrome (MetS) and obesity with kidney function.

METHODS

Subjects were 3,437 Korean adults of the Healthy Twin Study for cross-sectional relationships and 1,881 participants for longitudinal relationships (follow-up interval 3.7 ± 1.4 years). Obesity (body mass index ≥ 25 vs. <25 kg/m(2)), MetS, and chronic kidney disease [CKD; estimated glomerular filtration rate (eGFR) (<60 mL/min/1.73 m(2) using the modification of diet in renal disease study equation)] were categorized at baseline and follow-up.

RESULTS

The prevalence and incidence of chronic kidney disease were 2.5 and 3.3 %, respectively. Compared to individuals without obesity and MetS, prevalent CKD was associated with MetS regardless of weight status [adjusted odds ratio (AOR) 4.19 for those with MetS but without obesity; AOR 4.63 for those with MetS and obesity]. Incident CKD was associated with obesity regardless of baseline metabolic status (AOR 2.03 for those with obesity but without MetS; AOR 2.85 for those with obesity and MetS). MetS at follow-up was associated with incident CKD regardless of baseline MetS (AOR 2.42-2.52). Sex-adjusted bivariate analyses show inverse environmental correlations of the number of MetS components and BMI at baseline, with eGFR at baseline and follow-up (ρ E, -0.26 to -0.42, P < 0.001).

CONCLUSIONS

MetS predicts prevalent CKD regardless of obesity, and obesity predicts incident CKD regardless of baseline MetS. Incident CKD is also associated with MetS at follow-up regardless of baseline MetS. These associations appear to be explained by shared environmental factors.

Paraspinal muscle morphometry in cervical spondylotic myelopathy and its implications in clinicoradiological outcomes following central corpectomy

Object

The objective of this study was to assess the cross-sectional areas (CSAs) of the superficial, deep flexor (DF), and deep extensor (DE) paraspinal muscles in patients with cervical spondylotic myelopathy (CSM), and to evaluate their correlations with functional status and sagittal spinal alignment changes following central corpectomy with fusion and plating.

Methods

In this retrospective study of 67 patients who underwent central corpectomy with fusion and plating for CSM, the CSAs of the paraspinal muscles were calculated on the preoperative T2-weighted axial MR images and computed as ratios with respect to the corresponding vertebral body areas (VBAs) and as flexor/extensor CSA ratios. These ratios were then compared with those in the normative population and analyzed with respect to various clinicoradiological factors, including pain status, Nurick grade, and segmental angle change at follow-up (SACF).

Results

The mean CSA values for all muscle groups and the DF/DE ratio were significantly lower in the study cohort compared with an age- and sex-matched normative study group (p < 0.001). Among various independent variables tested in a multivariate regression analysis, increasing age and female sex significantly predicted a lower total extensor CSA/VBA ratio (p < 0.001), while a longer duration of symptoms significantly predicted a greater total flexor/total extensor CSA ratio (p = 0.02). In patients undergoing single-level corpectomy, graft subsidence had a positive correlation with SACF in all patients (p < 0.05), irrespective of the preoperative segmental angle and curvature, while in patients undergoing 2-level corpectomy, graft subsidence demonstrated such a correlation only in the subgroup with lordotic curvatures (p = 0.02). Among the muscle area ratios, the DF/DE ratio demonstrated a negative correlation with SACF in the subgroup with preoperative straight or kyphotic segmental angles (p = 0.04 in the single corpectomy group, p = 0.01 in the 2-level corpectomy group). There was no correlation of any of the muscle ratios with change in Nurick grade.

Conclusions

Patients with CSM demonstrate significant atrophy in all the flexor and extensor paraspinal muscles, and also suffer a reduction in the protective effect of a strong DF/DE CSA ratio. Worsening of this ratio significantly correlates with greater segmental kyphotic change in some patients. A physiological mechanism based on DF dysfunction is discussed to elucidate these findings that have implications in preventive physiotherapy and rehabilitation of patients with CSM. Considering that the influence of a muscle ratio was significant only in patients with hypolordosis, a subgroup that is known to have facetal ligament laxity, it may also be postulated that ligamentous support supersedes the influence of paraspinal muscles on postoperative sagittal alignment in CSM.

Variation in the MC4R gene is associated with bone phenotypes in elderly Swedish women

Osteoporosis is characterized by reduced bone mineral density (BMD) and increased fracture risk. Fat mass is a determinant of bone strength and both phenotypes have a strong genetic component. In this study, we examined the association between obesity associated polymorphisms (SNPs) with body composition, BMD, Ultrasound (QUS), fracture and biomarkers (Homocysteine (Hcy), folate, Vitamin D and Vitamin B12) for obesity and osteoporosis. Five common variants: rs17782313 and rs1770633 (melanocortin 4 receptor (MC4R); rs7566605 (insulin induced gene 2 (INSIG2); rs9939609 and rs1121980 (fat mass and obesity associated (FTO) were genotyped in 2 cohorts of Swedish women: PEAK-25 (age 25, n = 1061) and OPRA (age 75, n = 1044). Body mass index (BMI), total body fat and lean mass were strongly positively correlated with QUS and BMD in both cohorts (r² = 0.2–0.6). MC4R rs17782313 was associated with QUS in the OPRA cohort and individuals with the minor C-allele had higher values compared to T-allele homozygotes (TT vs. CT vs. CC: BUA: 100 vs. 103 vs. 103; p = 0.002); (SOS: 1521 vs. 1526 vs. 1524; p = 0.008); (Stiffness index: 69 vs. 73 vs. 74; p = 0.0006) after adjustment for confounders. They also had low folate (18 vs. 17 vs. 16; p = 0.03) and vitamin D (93 vs. 91 vs. 90; p = 0.03) and high Hcy levels (13.7 vs 14.4 vs. 14.5; p = 0.06). Fracture incidence was lower among women with the C-allele, (52% vs. 58%; p = 0.067). Variation in MC4R was not associated with BMD or body composition in either OPRA or PEAK-25. SNPs close to FTO and INSIG2 were not associated with any bone phenotypes in either cohort and FTO SNPs were only associated with body composition in PEAK-25 (p≤0.001). Our results suggest that genetic variation close to MC4R is associated with quantitative ultrasound and risk of fracture.

The genetic pleiotropy of musculoskeletal aging

Musculoskeletal aging is detrimental to multiple bodily functions and starts early, probably in the fourth decade of an individual's life. Sarcopenia is a health problem that is expected to only increase as a greater portion of the population lives longer; prevalence of the related musculoskeletal diseases is similarly expected to increase. Unraveling the biological and biomechanical associations and molecular mechanisms underlying these diseases represents a formidable challenge. There are two major problems making disentangling the biological complexity of musculoskeletal aging difficult: (a) it is a systemic, rather than "compartmental," problem, which should be approached accordingly, and (b) the aging per se is neither well defined nor reliably measurable. A unique challenge of studying any age-related condition is a need of distinguishing between the "norm" and "pathology," which are interwoven throughout the aging organism. We argue that detecting genes with pleiotropic functions in musculoskeletal aging is needed to provide insights into the potential biological mechanisms underlying inter-individual differences insusceptibility to the musculoskeletal diseases. However, exploring pleiotropic relationships among the system's components is challenging both methodologically and conceptually. We aimed to focus on genetic aspects of the cross-talk between muscle and its "neighboring" tissues and organs (tendon, bone, and cartilage), and to explore the role of genetics to find the new molecular links between skeletal muscle and other parts of the "musculoskeleton." Identification of significant genetic variants underlying the musculoskeletal system's aging is now possible more than ever due to the currently available advanced genomic technologies. In summary, a "holistic" genetic approach is needed to study the systems's normal functioning and the disease predisposition in order to improve musculoskeletal health.

Bivariate genome-wide association analyses of femoral neck bone geometry and appendicular lean mass

OBJECTIVE

Femoral neck geometric parameters (FNGPs), such as periosteal diameter (W), cross-sectional area (CSA), cortical thickness (CT), buckling ratio (BR), and section modulus (Z), are highly genetically correlated with body lean mass. However, the specific SNPs/genes shared by these phenotypes are largely unknown.

METHODS

To identify the specific SNPs/genes shared between FNGPs and appendicular lean mass (ALM), we performed an initial bivariate genome-wide association study (GWAS) by scanning ∼690,000 SNPs in 1,627 unrelated Han Chinese adults (802 males and 825 females) and a follow-up replicate study in 2,286 unrelated US Caucasians.

RESULTS

We identified 13 interesting SNPs that may be important for both FNGPs and ALM. Two SNPs, rs681900 located in the HK2 (hexokinase 2) gene and rs11859916 in the UMOD (uromodulin) gene, were bivariately associated with FNGPs and ALM (p = 7.58×10(-6) for ALM-BR and p = 2.93×10(-6) for ALM-W, respectively). The associations were then replicated in Caucasians, with corresponding p values of 0.024 for rs681900 and 0.047 for rs11859916. Meta-analyses yielded combined p values of 3.05×10(-6) and 2.31×10(-6) for rs681900 and rs11859916, respectively. Our findings are consistent with previous biological studies that implicated HK2 and UMOD in both FNGPs and ALM. Our study also identified a group of 11 contiguous SNPs, which spanned a region of ∼130 kb, were bivariately associated with FNGPs and ALM, with p values ranging from 3.06×10(-7) to 4.60×10(-6) for ALM-BR. The region contained two neighboring miRNA coding genes, MIR873 (MicroRNA873) and MIR876 (MicroRNA876).

CONCLUSION

Our study implicated HK2, UMOD, MIR873 and MIR876, as pleiotropic genes underlying variation of both FNGPs and ALM, thus suggesting their important functional roles in co-regulating both FNGPs and ALM.

Evidence for pleiotropic factors in genetics of the musculoskeletal system

There are both theoretical and empirical underpinnings that provide evidence that the musculoskeletal system develops, functions, and ages as a whole. Thus, the risk of osteoporotic fracture can be viewed as a function of loading conditions and the ability of the bone to withstand the load. Both bone loss (osteoporosis) and muscle wasting (sarcopenia) are the two sides of the same coin, an involution of the musculoskeletal system. Skeletal loads are dominated by muscle action; both bone and muscle share environmental, endocrine and paracrine influences. Muscle also has an endocrine function by producing bioactive molecules, which can contribute to homeostatic regulation of both bone and muscle. It also becomes clear that bone and muscle share genetic determinants; therefore the consideration of pleiotropy is an important aspect in the study of the genetics of osteoporosis and sarcopenia. The aim of this review is to provide an additional evidence for existence of the tight genetic co-regulation of muscles and bones, starting early in development and still evident in aging. Recently, important papers were published, including those dealing with the cellular mechanisms and anatomic substrate of bone mechanosensitivity. Further evidence has emerged suggesting that the relationship between skeletal muscle and bone parameters extends beyond the general paradigm of bone responses to mechanical loading. We provide insights into several pathways and single genes, which apparently have a biologically plausible pleiotropic effect on both bones and muscles; the list is continuing to grow. Understanding the crosstalk between muscles and bones will translate into a conceptual framework aimed at studying the pleiotropic genetic relationships in the etiology of complex musculoskeletal disease. We believe that further progress in understanding the common genetic etiology of osteoporosis and sarcopenia will provide valuable insight into important biological underpinnings for both musculoskeletal conditions. This may translate into new approaches to reduce the burden of both conditions, which are prevalent in the elderly population.

Genome-wide copy number variation association study suggested VPS13B gene for osteoporosis in Caucasians

Osteoporotic fracture (OF) is a serious outcome of osteoporosis. Important risk factors for OF include reduced bone mineral density and unstable bone structure. This genome-wide copy number variation association study suggested VPS13B gene for osteoporosis in Caucasians.

INTRODUCTION

Bone mineral density (BMD) and femoral neck cross-sectional geometric parameters (FNCSGPs) are under strong genetic control. DNA copy number variation (CNV) is an important source of genetic diversity for human diseases. This study aims to identify CNVs associated with BMD and FNCSGPs.

METHODS

Genome-wide CNV association analyses were conducted in 1,000 unrelated Caucasian subjects for BMD at the spine, hip, femoral neck, and for three FNCSGPs -cortical thickness (CT), cross-section area (CSA), and buckling ratio (BR). BMD was measured by dual energy X-ray absorptiometry (DEXA). CT, CSA, and BR were estimated using DEXA measurements. Affymetrix 500K arrays and copy number analysis tool was used to identify CNVs.

RESULTS

A CNV in VPS13B gene was significantly associated with spine, hip and FN BMDs, and CT, CSA, and BR (p < 0.05). Compared to subjects with two copies of the CNV, carriers of one copy had an average of 14.6%, 12.4%, and 13.6% higher spine, hip, and FN BMD, 20.0% thicker CT, 10.6% larger CSA, and 12.4% lower BR. Thus, a decrease of the CNV consistently produced stronger bone, thereby reducing osteoporotic fracture risk.

CONCLUSIONS

VPS13B gene, via affecting BMD and FNCSGPs, is a novel osteoporosis risk gene.

Bivariate whole-genome linkage scan for bone geometry and total body fat mass

To quantify the genetic correlations between total body fat mass (TBFM) and femoral neck geometric parameters (FNGPs) and, if possible, to detect the specific genomic regions shared by them, bivariate genetic analysis and bivariate whole-genome linkage scan were carried out in a large Caucasian population. All the phenotypes studied were significantly controlled by genetic factors (P < 0.001) with the heritabilities ranging from 0.45 to 0.68. Significantly genetic correlations were found between TBFM and CSA (cross-section area), W (sub-periosteal diameter), Z (section modulus) and CT (cortical thickness) except between TBFM and BR (buckling ratio). The peak bivariate LOD scores were 3.23 (20q12), 2.47 (20p11), 3.19 (6q27), 1.68 (20p12), and 2.47 (7q11) for the five pairs of TBFM and BR, CSA, CT, W, and Z in the entire sample, respectively. Gender-specific bivariate linkage evidences were also found for the five pairs. 6p25 had complete pleiotropic effects on the variations of TBFM & Z in the female sub-population, and 6q27 and 17q11 had coincident linkages for TBFM & CSA and TBFM & Z in the entire population. We identified moderate genetic correlations and several shared genomic regions between TBFM and FNGPs in a large Caucasian population.

Bivariate genome-wide linkage analysis of femoral bone traits and leg lean mass: Framingham study

The risk of osteoporotic fracture is a function of both applied muscle mass and bone tissue distribution. Leg lean mass (LLM) and femoral bone geometry are both known to have substantial genetic components. Therefore, we estimated shared heritability (h(2)) and performed linkage analysis to identify chromosomal regions governing both LLM and bone geometry. A genome-wide scan (using 636 microsatellite markers) for linkage analyses was performed on 1346 adults from 327 extended families of the Framingham study. DXA measures were LLM, femoral neck length, neck-shaft angle (NSA), subperiosteal width, cross-sectional area (CSA), and section modulus (Z) at the femoral narrow neck and shaft (S) regions. Variance component linkage analysis was performed on normalized residuals (adjusted for age, height, BMI, and estrogen status in women). The results indicated substantial h(2) for LLM (0.42 +/- 0.07) that was comparable to bone geometry traits. Phenotypic correlations between LLM and bone geometry phenotypes ranged from 0.033 with NSA (p > 0.05) to 0.251 with S_Z (p < 0.001); genetic correlations ranged from 0.087 (NSA, p > 0.05) to 0.454 (S_Z, p < 0.001). Univariate linkage analysis of covariate-adjusted LLM identified no chromosomal regions with LOD scores >or=2.0; however, bivariate analysis identified two loci with LOD scores >3.0, shared by LLM with S_CSA on chromosome 12p12.3-12p13.2, and with NSA, on 14q21.3-22.1. In conclusion, we identified chromosomal regions potentially linked to both LLM and femoral bone geometry. Identification and subsequent characterization of these shared loci may further elucidate the genetic contributions to both osteoporosis and sarcopenia.

Genetics of the musculoskeletal system: a pleiotropic approach

The risk of osteoporotic fracture can be viewed as a function of loading conditions and the ability of the bone to withstand the load. Skeletal loads are dominated by muscle action. Recently, it has become clear that bone and muscle share genetic determinants. Involution of the musculoskeletal system manifests as bone loss (osteoporosis) and muscle wasting (sarcopenia). Therefore, the consideration of pleiotropy is an important aspect in the study of the genetics of osteoporosis and sarcopenia. This Perspective will provide the evidence for a shared genetic influence on bone and muscle. We will start with an overview of accumulating evidence that physical exercise produces effects on the adult skeleton, seeking to unravel some of the contradictory findings published thus far. We will provide indications that there are pleiotropic relationships between bone structure/mass and muscle mass/function. Finally, we will offer some insights and practical recommendations as to the value of studying shared genetic factors and will explore possible directions for future research. We consider several related questions that together comprise the general paradigm of bone responses to mechanical loading and the relationship between muscle strength and bone parameters, including the genetic factors that modulate these responses. We believe that further progress in understanding the common genetic etiology of osteoporosis and sarcopenia will provide valuable insight into important biological underpinnings for both conditions and may translate into new approaches to reduce the burdens of both conditions through improved diagnosis, prevention, and early targeted treatment.

Genetic and environmental influence on structural strength of weight-bearing and non-weight-bearing bone: a twin study

A bivariate genetic analysis among 217 older female twin pairs showed that, although the structural strength of tibia and radius are mainly regulated by same genetic and environmental factors, the tibia is more affected by environment.

INTRODUCTION

The habitual loading environment of the bone may modulate the relative contribution of genetic and environmental factors to bone structure. The purpose of this study was to estimate the contribution of the common and site-specific genetic and environmental factors to interindividual variation in compressive structural strength of the weight-bearing tibia and non-weight-bearing radius.

MATERIALS AND METHODS

pQCT scans were obtained from both members of 103 monozygotic (MZ) and 114 dizygotic (DZ) 63- to 76-yr-old female twin pairs to estimate the compressive strength of the distal tibia and distal radius. Quantitative genetic models were used to decompose the phenotypic variance into additive genetic, shared environmental, and individual environmental effects at each bone site and to study whether these bone sites share genetic or environmental effects.

RESULTS

The MZ and DZ twins did not differ in mean age, height, weight, or bone structural strength. The age-adjusted Cholesky model showed that additive genetic factors accounted for 83% (95% CI, 77-88%) of the variance in radial strength and 61% (95% CI, 52-69%) of the variance in tibial strength, and these were fully correlated. A shared environmental factor accounted for 15% (95% CI, 10-20%) of tibial strength. An individual environmental factor accounted for 17% (95% CI, 12-23%) of the variance in radial strength and 10% (95% CI, 5-17%) of the variance in tibial strength. The relative contribution of an individual environmental factor specific to tibial strength was 14% (95% CI, 11-18%).

CONCLUSIONS

The results suggest that, in older women, the majority of the individual differences in the compressive structural strength of the forearm and leg are regulated by genetic and environmental factors that are common to both bone sites. However, the relative importance of environmental factors was greater for the weight-bearing tibia than for the non-weight-bearing radius. Thus, the heritability of bone strength seems to vary between skeletal sites according to differences in the typical loading environment.

Pleiotropy and heterogeneity in the expression of bone strength-related phenotypes in extended pedigrees

Genetic analysis in 3,535 relative pairs from extended multigenerational families of African heritage showed that volumetric BMD is a highly heritable polygenic trait that is under compartment-specific genetic regulation. The majority of the phenotypic variation in bone size and volumetric BMD also seems to be strongly influenced by distinct genes for each trait.

INTRODUCTION

BMD and bone size contribute to bone strength and the risk of fracture. Little is known about the genetic architecture of QCT measures of volumetric BMD and bone size. We studied the contribution of genes, shared genes (pleiotropy), and shared environment to cortical and trabecular volumetric BMD and bone size using variance components analysis.

MATERIALS AND METHODS

A total of 471 individuals >or=18 yr of age (mean, 43 yr) from eight multigenerational Afro-Caribbean families (mean family size > 50; 3535 relative pairs) underwent a peripheral QCT scan of the radius and tibia and anthropometry.

RESULTS

Strong positive genetic correlations were observed for trabecular or cortical BMD measured at the tibia and radius (rho(G) > 0.82, p < 0.01), but not between trabecular and cortical BMD measured within the same anatomical site. Genetic correlations between volumetric BMD and bone length and circumference were also not statistically significant.

CONCLUSIONS

BMD is a highly heritable polygenic trait that is under compartment-specific genetic regulation. The majority of the phenotypic variation in skeletal size and density seems to be strongly influenced by distinct sets of genes for each trait.

Linkage of genes to total lean body mass in normal women

BACKGROUND

Total lean body mass (LEAN-tot) is one of the three major components of body weight. Its deterioration is a risk factor for frailty. Despite this, there are few studies examining the contribution of genetic factors.

OBJECTIVE

Our objective was to examine the contribution of genetic factors for LEAN-tot variation, including a genome-wide search for the genes.

RESEARCH METHODS

Dual-energy x-ray absorptiometry measurements of LEAN-tot were obtained from each of the 3180 United Kingdom females (509 monozygotic and 1081 dizygotic twin pairs). Contribution of genetic factors was assessed using variance component analysis. A genome-wide linkage analysis was performed on the dizygotic twins using a modified version of the Haseman-Elston method.

RESULTS

Age, body height, total fat, and bone mass were correlated with LEAN-tot, and commonly explained 52% of the LEAN-tot variation. The crude heritability estimate was 74.0 +/- 4.0%, after adjustment for the aforementioned factors; 65.2 +/- 4.6% was attributable to independent genetic effects. Significant (P < 0.001) genetic correlations were found between LEAN-tot and bone mass, and LEAN-tot and total fat. Adjusted only for age, LEAN-tot showed no significant linkage. After adjustment for all covariates, significant linkage (LOD = 4.49 and 3.62) was observed at chromosome 12q24.3 and 14q22.3, respectively. Additional peaks of interest were on 7p15.3-15.1 (LOD = 2.86) and 8p22 (LOD = 2.83).

CONCLUSIONS

LEAN-tot measured by dual-energy x-ray absorptiometry is highly heritable, independent of other body measures. This first genomic search for genes associated with the lean component of body mass suggests significant linkage to quantitative trait loci on chromosomes 12 and 14.

Bivariate whole genome linkage analysis for femoral neck geometric parameters and total body lean mass

A genome-wide bivariate analysis was conducted for femoral neck GPs and TBLM in a large white sample. We found QTLs shared by GPs and TBLM in the total sample and the sex-specific samples. QTLs with potential pleiotropy were also disclosed.

INTRODUCTION

Previous studies have suggested that femoral neck cross-section geometric parameters (FNCS-GPs), including periosteal diameter (W), cross-sectional area (CSA), cortical thickness (CT), buckling ratio (BR), and section modulus (Z), are genetically correlated with total body lean mass (TBLM). However, the shared genetic factors between them are unknown.

MATERIALS AND METHODS

To identify the specific QTLs shared by FNCS-GPs and TBLM, we performed bivariate whole genome linkage analysis (WGLA) in a large sample of 451 white families made up of 4498 subjects.

RESULTS

Multipoint bivariate linkage analyses for 22 autosomes showed evidence of suggestive or significant linkages (thresholds of LOD = 2.3 and 3.7, respectively) to chromosomes 3q12 and 20q13 in the entire sample, 6p25 and 10q24 in women, and 4p15, 5q34-35 and 7q21 in men. Two-point linkage analyses for chromosome X showed strong linkage to Xp22.13, Xp11.4, Xq22.3, Xq23-24, and Xq25. Complete pleiotropy was identified on 10q24 and 5q35 for TBLM and BR in women and for TBLM and CT in men, respectively. Furthermore, chromosomes 5q34-35, 7q21, 10q24, 20q13, Xp22.13, Xp11.4, and Xq25 are also of importance because of their linkage to multiple trait pairs. For example, linkage to chromosome 10q24 was found for TBLM x W (LOD = 2.31), TBLM x CT (LOD = 2.51), TBLM x CSA (LOD = 2.51), TBLM x BR (LOD = 2.64), and TBLM x Z (LOD = 2.55) in women.

CONCLUSIONS

In this study, we identified several genomic regions (e.g., 3q12 and 20q13) that seem to be linked to both FNCS-GPs and TBLM. These regions are of interesting because they may harbor genes that may contribute to variation in both FNCS-GPs and TBLM.