Kyphosis and paraspinal muscle composition in older men: a cross-sectional study for the osteoporotic fractures in men (MrOS) research group

Background The prevalence of hyperkyphosis is increased in older men; however, risk factors other than age and vertebral fractures are not well established. We previously reported that poor paraspinal muscle composition contributes to more severe kyphosis in a cohort of both older men and women. Methods To specifically evaluate this association in older men, we conducted a cross-sectional study to evaluate the association of paraspinal muscle composition and degree of thoracic kyphosis in an analytic cohort of 475 randomly selected participants from the Osteoporotic Fractures in Men (MrOS) study with baseline abdominal quantitative computed tomography (QCT) scans and plain thoracic radiographs. Baseline abdominal QCT scans were used to obtain abdominal body composition measurements of paraspinal muscle and adipose tissue distribution. Supine lateral spine radiographs were used to measure Cobb angle of kyphosis. We examined the linear association of muscle volume, fat volume and kyphosis using loess plots. Multivariate linear models were used to investigate the association between muscle and kyphosis using total muscle volume, as well as individual components of the total muscle volume, including adipose and muscle compartments alone, controlling for age, height, vertebral fractures, and total hip bone mineral density (BMD). We examined these associations among those with no prevalent vertebral fracture and those with BMI < 30 kg/m2. Results Among men in the analytic cohort, means (SD) were 74 (SD = 5.9) years for age, and 37.5 (SD = 11.9) degrees for Cobb angle of kyphosis. Men in the lowest tertile of total paraspinal muscle volume had greater mean Cobb angle than men in the highest tertile, although test of linear trend across tertiles did not reach statistical significance. Neither lower paraspinal skeletal muscle volume (p-trend = 0.08), or IMAT (p-trend = 0.96) was associated with greater kyphosis. Results were similar among those with no prevalent vertebral fractures. However, among men with BMI < 30 kg/m2, those in the lowest tertile of paraspinal muscle volume had greater adjusted mean kyphosis (40.0, 95% CI: 37.8 – 42.1) compared to the highest tertile (36.3, 95% CI: 34.2 – 38.4). Conclusions These results suggest that differences in body composition may potentially influence kyphosis.


Background
Hyperkyphosis, an exaggerated curvature in the thoracic spine, is commonly observed in the older population although its causes have not been well investigated. Greater degrees of kyphosis impairs physical function [1][2][3][4][5] and well-being [6], and increases load on the vertebral bodies [7,8] that potentially could increase the risk for vertebral compression fractures. A recent study confirms the need for more comprehensive assessment of health outcomes in older adults with greater degree of kyphosis [9]. Hyperkyphosis is considered a problem that primarily affects older women, even though the prevalence of hyperkyphosis in older men is estimated to range from 15-40%, depending upon how kyphosis is defined [6,10]. While commonly recognized risk factors for hyperkyphosis in women include advanced age, low bone mass and underlying vertebral fractures [11], only 1/3 of older women with the most severe kyphosis have vertebral compression fractures (VCFs) [12], suggesting other factors contribute to kyphosis that have not been well defined. Furthermore, it is not known whether men have similar risk factors for hyperkyphosis.
Recently Katzman, et al. determined that paraspinal trunk muscle composition is associated with hyperkyphosis in a cohort of older community-dwelling men and women [10]. Among individuals with a greater degree of kyphosis, there was more adipose accumulation within the paraspinal muscles. Imaging methods that include computed tomography and MRI are the most accurate methods available for quantifying skeletal muscle composition in vivo, and these methods permit measurement of skeletal muscle and adipose tissues within the muscle fascia [13]. While MRI does not expose an individual to radiation, it overestimates the amount of adipose tissue within the erector spinae muscle fascia [14], suggesting that CT imaging of the erector spinae may be more precise than MRI. Attenuation refers to the energy of the x-ray beam as it passes through the tissue, and whereas all tissue attenuates the x-ray beam, muscle attenuates more than adipose. Attenuation indicates the density of muscle with lower attenuation indicative of greater lipid infiltration [15]. An attenuation value, called a Hounsfield Unit (HU), is assigned to each volume element (or voxel) in the CT image. Previous work demonstrated that muscle and adipose tissues have distinct HU ranges [15][16][17][18]. Thus, CT imaging permits segmentation of the total muscle volume contained within the fascial borders into its skeletal muscle and adipose tissue components.
Skeletal muscle composition changes with aging and differs according to sex. During aging, adipose tissue volume increases, skeletal muscle tissue volume declines, and the remaining muscle tissue is susceptible to lipid infiltration [19]. Younger adults have greater muscle volume and higher muscle attenuation than older adults [20]. Moreover, sex also affects muscle composition of the skeletal muscles, particularly in the paraspinal muscles, and men have higher paraspinal muscle attenuation than women [21,22]. Given these age-related changes in paraspinal muscle composition and the differences in paraspinal attenuation in men compared with women, further investigation of the association of paraspinal muscle composition with kyphosis in older men is warranted.
We conducted a cross-sectional study to investigate differences in degree of thoracic kyphosis according to paraspinal muscle composition in older men using data from the Osteoporotic Fractures in Men (MrOS) Study. We hypothesized that both lower paraspinal muscle volume and higher paraspinal intermuscular adipose volume would be associated with greater degree of Cobb angle of kyphosis, independent of known and potential covariates age, weight, height, abdominal fat, hip bone mineral density, and prevalent vertebral fractures. We further hypothesized that obese men would have more paraspinal muscle volume than men with BMI < 30 kg/m 2 that would have a protective effect on Cobb angle of kyphosis, thus warranting examining these associations according to BMI.

Participants
The Osteoporotic Fractures in Men (MrOS) cohort consists of 5,994 men aged 65 or older who were recruited from 6 academic medical centers in Birmingham, AL; Minneapolis, MN; Palo Alto, CA; Pittsburgh, PA; Portland, OR; and San Diego, CA [23,24]. The primary purpose of MrOS has been to describe the risk factors for osteoporotic fractures in older men. Funding for MrOS began in July 1999 and the study is now its 13th year. At baseline, all men were referred for plain thoracic and lumbar radiographs and 3786 men were referred for quantitative computed tomography (QCT) scans of the hips, lumbar spine and abdomen. Abdominal body composition measures were performed for a random sample of 667 participants in the sample of 3786 who had been referred for QCT scans at baseline [25]. The analytic cohort for our study included the 475 randomly selected participants from the MrOS study with both baseline abdominal body composition measures from QCT scans and plain thoracic radiographs.
Approval of the conduct of the MrOS study was obtained from the Institutional Review Board of each of the participating medical centers and written informed consent was obtained from all study participants.

Muscle variables
Baseline abdominal QCT scans were used to obtain abdominal body composition measurements of abdominal wall muscle and adipose tissue distribution. Abdominal scans were obtained using a standardized protocol which specified scanning from the mid-L3 to the mid-L5 vertebra at 5-mm slice thickness with participants in the supine position. We completed abdominal body composition measurements on the scan at the L4-L5 intervertebral disc level. Scans were processed by two specialists at Oregon Health and Science University (OHSU) using a standardized protocol as previously described [26]. Briefly, a spline tool was used to trace the outer and inner fascial borders of the abdominal wall musculature comprising 6 bilateral muscle groups -the rectus abdominis, external oblique, internal oblique, transverse abdominis, psoas and paraspinals. The paraspinal muscle group included the muscles on each side of the vertebral body in the lumbar region located lateral to the lumbar vertebral spinous process and posterior to the lumbar vertebral transverse processes comprising the longissimus, iliocostalis and multifidus muscles. Intermuscular adipose tissue (IMAT) included adipose tissue located around muscle bundles within each muscle group (Figure 1). Attenuation ranges used to estimate volumes (cm 3 ) of muscle were 0 to 100 HU [15,18] and of IMAT were −190 to −30 HU [16,17]. Thus, for all muscle groups studied, we had three volume variables: skeletal muscle only, IMAT, and total volume which is the sum of skeletal muscle plus IMAT volume. Muscle attenuation (HU) was obtained for the total volume and the skeletal muscle only volume. To account for variability in image attenuation across scanners, the HU in the scans were scaled to the calibration standard for water density (0 mg/cm 3 ). Interand intra-rater reliability was monitored throughout image processing with intraclass correlation coefficients (ICC). Final intra-rater ICCs for all measures were ≥0.98 for rater 1 and ≥0.95 for rater 2, and final inter-rater ICCs were all ≥0.96.

Cobb angle of kyphosis measurements
The Cobb angle of kyphosis was measured from supine lateral spine radiographs taken during baseline visit 1. We used the modified Cobb method [27] with a fixed cut-off of T4 and T12, largely because T1 to T3 are usually not well visualized on lateral spine films due to the overlying shoulders and scapulae interfering with the projection. Technicians placed 6 points corresponding to the 4 corners of the vertebral body and the midpoints of the endplates on each vertebral body from T4 to T12. From the superior surface of T4 and inferior surface of T12, a computerized digitization program erected perpendicular lines whose intersection was the kyphotic angle ( Figure 2). If for any reason T4 or T12 were not visible, the next adjacent visible vertebra (T5 if T4 not visible or T11 if T12 not visible) was used as an alternative. The ICC for Cobb angle of kyphosis has been previously reported 0.99 [27].

Vertebral fracture measurements
Vertebral fractures were adjudicated from baseline lateral lumbar and thoracic spine radiographs based upon a previously developed protocol [28]. The SpineAnalyzer™ (Optasia Medical Ltd., Cheadle, UK) workflow tool was used to automate placement of 6-point morphometric points on each radiograph in order to efficiently identify the vertebral bodies for evaluation. An expert radiologist then graded each vertebra as normal (0), mild (1), moderate (2), severe (3) using the well-established, highly reliable semi-quantitative Genant method [29]. Prevalent vertebral fractures were defined as grade 2 or 3.

Anthropometric measurements
Body weight and height were taken at baseline by an examiner using standard equipment, including a Harpenden stadiometer and balance beam scale while wearing indoor street clothing without shoes [23]. Body mass index (BMI) was calculated as weight divided by the square of height in meters (kg/m 2 ). Bone mineral density (BMD) was also measured at baseline in the proximal femur and lumbar spine using DXA measured by Hologic QDR 4500 densitometers.

Statistical analysis
We built models to control for potential confounding using an empirical approach. We examined the association between an a priori group of potential confounders and Cobb angle and based on the change in the regression coefficient, we used p value <0.1 to select confounders including age, height, vertebral fractures, and total hip bone mineral density (BMD). Loess plots were used initially to determine the shape of the relations of muscle volume, adipose volume and Cobb angle. Because associations were not linear, muscle mass measurements were divided into quintiles, quartiles and tertiles to investigate non-linear trends in kyphosis, but due to the limited sample size, we used tertiles for our analysis. Multivariate linear models were used to investigate the association between the muscle variables and kyphosis using total muscle volume, as well as individual components of the total muscle volume that included the IMAT and skeletal muscle compartments alone. All multivariable analyses controlled for age, height, vertebral fractures, and total hip BMD. These associations were also investigated separately for the total abdominal (rectus abdominis, transversus abdominis, internal obliques and external obliques), psoas and paraspinal muscle groups. We examined these associations among those with no prevalent vertebral fracture (n = 429) and among those with BMI < 30 kg/m 2 (n = 394). We conducted sensitivity analyses with different adjustment variables as well as an alternate definition of hyperkyphosis >40 degrees.

Results
The subject characteristics of the overall MrOS cohort, means (standard deviation), were 74 (SD = 5.9) years for age, 174 (SD = 6.8) cms for height, 0.96 (SD = 0.14) g/cm 2 for total hip BMD, 27 (SD = 3.8) kg/m 2 for BMI, and 7.5% had prevalent vertebral fractures. In the analytic cohort of 475 men, means (SD) were 74 (SD = 5.9) years for age and 37.5 (SD = 11.9) degrees for Cobb angle of kyphosis (Table 1). Men with BMI < 30 kg/m 2 were older, had greater Cobb angle, lower hip BMD, and lower total abdominal muscle volume, but no difference in total paraspinal or total psoas muscle volume. When we separated total muscle volume into its skeletal muscle tissue and IMAT components, skeletal muscle volume but not attenuation in all of the abdominal wall muscle groups was greater, and IMAT in the abdominal and psoas muscle groups was lower in men with BMI < 30 kg/m 2 .
Men in the lowest tertile of total paraspinal muscle volume had greater mean Cobb angle than men in the highest tertile (Table 2), although the test of linear trend across tertiles did not reach statistical significance. When we further separated total paraspinal volume into its skeletal muscle tissue and IMAT components, neither lower paraspinal skeletal muscle volume (p-trend = 0.08), or IMAT (p-trend = 0.96) was associated with greater kyphosis. Results were similar among those with no prevalent vertebral fractures. However, among men with BMI < 30 kg/m 2 , those in the lowest tertile of paraspinal skeletal muscle volume had greater adjusted mean kyphosis (40. (Table 3), and no significant associations of the adjusted mean Cobb angle of kyphosis according to abdominal and psoas skeletal muscle and IMAT volumes, p > .05 (data not shown).
The sensitivity analysis where we replaced paraspinal IMAT with paraspinal muscle attenuation in models that also included paraspinal muscle volume showed similar findings. Specifically, neither paraspinal skeletal muscle volume or paraspinal muscle attenuation showed significant associations with Cobb angle, with the exception of paraspinal muscle volume among those with men with BMI < 30 (data not shown). Consistent with the results in Table 2, among men with BMI < 30 kg/m 2 , having less paraspinal muscle volume is associated with greater Cobb angle. In the multivariate adjusted model, when Cobb angle was dichotomized > 40 degrees, there was no association with paraspinal muscle volume or adipose volume, p > 0.05 (data not shown).

Discussion
We found that paraspinal muscle volume influenced kyphosis among those with BMI < 30 kg/m 2 . Lower paraspinal muscle volume may reduce the capacity of paraspinal muscles to extend and stabilize the spine, resulting in increased kyphosis. In contrast, muscle volume in the abdominal and psoas muscles that predominantly flex the spine, did not influence kyphosis. Furthermore, we found that paraspinal muscle volume, rather than paraspinal muscle attenuation, a measure of intramuscular adipose tissue, had greater influence on kyphosis than IMAT in our cohort of older men.
These results are in contrast to our previous finding among older men and women in the Health Aging and Body Composition (Health ABC) Study that paraspinal muscle attenuation, not muscle volume, is associated with kyphosis. It is possible these differences could be because the Health ABC Study included men and women, and our MrOS cohort included men only who have greater muscle volume and attenuation than women. Our results could also be explained by differences in image acquisition and processing of the computed tomography scans that were used to obtain the muscle composition measurements; there are no standard methods, making comparisons of the Health ABC and MrOS findings difficult. Another hindrance to comparing results is that in our previously reported Health ABC study, muscle volume was categorized as the volume of non-bone, nonadipose tissue within the fascial plane of the specified muscle groups, whereas in our current MrOS study, muscle volume was defined as muscle and adipose tissue within the plane of the fascia inclusive of intramuscular adipose tissue. Furthermore,   IMAT was not measured or included in the Health ABC analysis, whereas we controlled for IMAT in our analyses. However, when we did investigate paraspinal muscle attenuation in MrOS, it was not associated with Cobb angle of kyphosis. When we performed a sensitivity analysis and replaced IMAT with paraspinal muscle attenuation, the results did not appreciably change. In addition, while we used a continuous measure of Cobb angle of kyphosis in MrOS, in Health ABC we used a dichotomous variable of hyperkyphosis defined as >40 degrees, which could reflect a threshold effect of paraspinal muscle attenuation on Cobb angle of kyphosis. Nonetheless, we found no evidence of a threshold effect within the MrOS cohort. Finally, Health ABC measurements of Cobb angle were obtained from supine computed tomography scout scans, whereas in MrOS we obtained these measurements from supine lateral spine radiographs, and we do not know how well these measurements compare. Vertebral fractures, commonly thought to be the cause of excessive kyphosis, did not influence the association between paraspinal muscle volume, paraspinal IMAT and Cobb angle of kyphosis. These results provide further evidence that factors other than vertebral fractures are responsible for greater degree of kyphosis. We found that paraspinal muscle volume, a factor that is potentially modifiable with exercise, is associated with Cobb angle of kyphosis in our cohort of older men with BMI <30 kg/m 2 versus high BMI.
Whereas decreased paraspinal muscle attenuation, not muscle volume, has been associated with decline in physical performance among older adults [30,31], in our study, we found reduced paraspinal muscle volume, not muscle attenuation, influenced kyphosis. Resistance exercise trials targeting lower extremity muscle groups in older adults have reported that highintensity resistance exercise produces large increase in muscle strength [32,33]. Improvement in lower extremity muscle strength following high intensity strengthening has also been associated with increased muscle volume and muscle attenuation [32][33][34][35][36]. Several randomized controlled trials of exercise have included back extensor muscle strengthening to reduce kyphosis [37][38][39][40][41], and it is possible that strengthening the back extensor muscles could improve paraspinal muscle volume and attenuation. Unfortunately the effects of back extensor muscle strengthening on paraspinal muscle volume or attenuation have not been investigated, and it is not known whether improving paraspinal muscle volume or attenuation reduces kyphosis.
While the phenotype of a smaller-framed older woman with hyperkyphosis is well-recognized, our findings suggest this phenotype also exists in older men. Body mass index, a proxy for human body fat based on an individual's weight and height, influenced kyphosis in our cohort of men. We found 83% of the men had BMI < 30 kg/m 2 , and among these men, low paraspinal muscle volume was associated with greater Cobb angle of kyphosis. In fact, among those with BMI < 30 kg/m 2 and in the lowest tertile of paraspinal muscle volume, the mean kyphosis was 40.0 degrees (95% CI = 37.8-42.1), the threshold for defining hyperkyphosis in older adults [42,43]. Hyperkyphosis, low paraspinal muscle volume and BMI < 30 kg/m 2 may reflect an underlying geriatric syndrome that impacts both older men and women. This kyphotic phenotype may place older men at risk for vertebral compression fractures due to the increased spinal load and impaired functional outcomes that are associated with greater kyphosis. While obesity did appear to have a protective effect on kyphosis in our cohort, intraorgan fat has been associated with poor health outcomes in older adults [44]. We are not advocating for obesity as a means to control excessive kyphosis, rather the need for more comprehensive assessment of health status among older adults with hyperkyphosis.

Limitations
This study has a number of strengths including that the MrOS cohort includes community-dwelling older men residing in multiple geographic areas of the United States who were not preselected for hyperkyphosis. The kyphosis and the paraspinal muscle measurements were assessed at baseline with high reliability in a subgroup of study subjects. However, this study also has a number of limitations including possible misclassification of muscle measurements depending upon the exact location of the QCT scan slice. The paraspinal muscle group does increase in size closer its attachment sites, therefore our paraspinal muscle volume estimates measured at the L4-L5 intervertebral disc space may be underestimated, as would any association of paraspinal muscle volume and kyphosis found in this study. Another potential limitation is that we could not measure muscle volumes in the thoracic spine due to the CT technology used. However, our results suggest that lumbar paraspinal muscles serve as a proxy for muscle degeneration in the upper spine as well. We also used BMD measurements from the hip even though the muscle and kyphosis measurements were made in the spine. Hip BMD is often used in place of spine BMD because there is less artifact in the hip and the lumbar spine BMD from DXA can be overestimated because of extravertebral calcification (e.g. osteophytes). Hip BMD thus provides a more accurate measurement, particularly in elderly subjects. Moreover, when we replaced hip BMD with volumetric measurements of spine BMD with QCT, there was no appreciable difference in the results. The MrOS population may have degenerative spine conditions that could increase kyphosis and our regression analysis did not include adjustment for these degenerative changes. However, this is beyond the scope of this study. We chose to measure global Cobb angle using T4 and T12 endplates rather than other methods such as the centroid method that is less affected by endplate tilt [45] because the global method using the endplates of T4 and T12 is reliable and one of the most commonly used radiographic methods for measuring thoracic kyphosis.
In future studies, we may consider using other methods that reflect changes in the middle of the curve. Finally, supine measurements of Cobb angle of kyphosis underestimate the degree of kyphosis compared to standing measurements, particularly among those with greater degree of kyphosis [27]. However, this would serve to underestimate any effects we found.

Conclusions
In summary, we found that paraspinal muscle volume influenced kyphosis among older men with BMI < 30 kg/m 2 , and lower paraspinal muscle volume was associated with greater degree of thoracic kyphosis. Interestingly, prevalent vertebral fractures commonly thought to be the cause of excessive kyphosis, did not influence the association between paraspinal muscle volume, paraspinal IMAT and Cobb angle of kyphosis. Additional studies are now warranted that evaluate if low paraspinal muscle volume contributes to worsening of kyphosis in this population.