Estimating apparent maximum muscle stress of trunk extensor muscles in older adults using subject-specific musculoskeletal models.

Maximum muscle stress (MMS) is a critical parameter in musculoskeletal modeling, defining the maximum force that a muscle of given size can produce. However, a wide range of MMS values have been reported in literature, and few studies have estimated MMS in trunk muscles. Due to widespread use of musculoskeletal models in studies of the spine and trunk, there is a need to determine reasonable magnitude and range of trunk MMS. We measured trunk extension strength in 49 participants over 65 years of age, surveyed participants about low back pain, and acquired quantitative computed tomography (QCT) scans of their lumbar spines. Trunk muscle morphology was assessed from QCT scans and used to create a subject-specific musculoskeletal model for each participant. Model-predicted extension strength was computed using a trunk muscle MMS of 100 N/cm2 . The MMS of each subject-specific model was then adjusted until the measured strength matched the model-predicted strength (±20 N). We found that measured trunk extension strength was significantly higher in men. With the initial constant MMS value, the musculoskeletal model generally over-predicted trunk extension strength. By adjusting MMS on a subject-specific basis, we found apparent MMS values ranging from 40 to 130 N/cm2 , with an average of 75.5 N/cm2 for both men and women. Subjects with low back pain had lower apparent MMS than subjects with no back pain. This work incorporates a unique approach to estimate subject-specific trunk MMS values via musculoskeletal modeling and provides a useful insight into MMS variation. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:498-505, 2018.

Spinal Loading Patterns From Biomechanical Modeling Explain the High Incidence of Vertebral Fractures in the Thoracolumbar Region.

Vertebral fractures occur most frequently in the mid-thoracic and thoracolumbar regions of the spine, yet the reasons for this site-specific occurrence are not known. Our working hypothesis is that the locations of vertebral fracture may be explained by the pattern of spine loading, such that during daily activities the mid-thoracic and thoracolumbar regions experience preferentially higher mechanical loading compared to other spine regions. To test this hypothesis, we used a female musculoskeletal model of the full thoracolumbar spine and rib cage to estimate the variation in vertebral compressive loads and associated factor-of-risk (load-to-strength ratio) throughout the spine for 119 activities of daily living, while also parametrically varying spine curvature (high, average, low, and zero thoracic kyphosis models). We found that nearly all activities produced loading peaks in the thoracolumbar and lower lumbar regions of the spine, but that the highest factor-of-risk values generally occurred in the thoracolumbar region of the spine because these vertebrae had lower compressive strength than vertebrae in the lumbar spine. The peaks in compressive loading and factor-of-risk in the thoracolumbar region were accentuated by increasing thoracic kyphosis. Activation of the multifidus muscle fascicles selectively in the thoracolumbar region appeared to be the main contributor to the relatively high vertebral compressive loading in the thoracolumbar spine. In summary, by using advanced musculoskeletal modeling to estimate vertebral loading throughout the spine, this study provides a biomechanical mechanism for the higher incidence of fractures in thoracolumbar vertebrae compared to other spinal regions. © 2017 American Society for Bone and Mineral Research.

Incorporation of CT-based measurements of trunk anatomy into subject-specific musculoskeletal models of the spine influences vertebral loading predictions.

We created subject-specific musculoskeletal models of the thoracolumbar spine by incorporating spine curvature and muscle morphology measurements from computed tomography (CT) scans to determine the degree to which vertebral compressive and shear loading estimates are sensitive to variations in trunk anatomy. We measured spine curvature and trunk muscle morphology using spine CT scans of 125 men, and then created four different thoracolumbar spine models for each person: (i) height and weight adjusted (Ht/Wt models); (ii) height, weight, and spine curvature adjusted (+C models); (iii) height, weight, and muscle morphology adjusted (+M models); and (iv) height, weight, spine curvature, and muscle morphology adjusted (+CM models). We determined vertebral compressive and shear loading at three regions of the spine (T8, T12, and L3) for four different activities. Vertebral compressive loads predicted by the subject-specific CT-based musculoskeletal models were between 54% lower to 45% higher from those estimated using musculoskeletal models adjusted only for subject height and weight. The impact of subject-specific information on vertebral loading estimates varied with the activity and spinal region. Vertebral loading estimates were more sensitive to incorporation of subject-specific spinal curvature than subject-specific muscle morphology. Our results indicate that individual variations in spine curvature and trunk muscle morphology can have a major impact on estimated vertebral compressive and shear loads, and thus should be accounted for when estimating subject-specific vertebral loading. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:2164-2173, 2017.

Vertebral Volumetric Bone Density and Strength Are Impaired in Women With Low-Weight and Atypical Anorexia Nervosa.

Context: Areal bone mineral density (BMD) is lower, particularly at the spine, in low-weight women with anorexia nervosa (AN). However, little is known about vertebral integral volumetric BMD (Int.vBMD) or vertebral strength across the AN weight spectrum, including "atypical" AN [body mass index (BMI) ≥18.5 kg/m2]. Objective: To investigate Int.vBMD and vertebral strength, and their determinants, across the AN weight spectrum. Design: Cross-sectional observational study. Setting: Clinical research center. Participants: 153 women (age 18 to 45): 64 with low-weight AN (BMI <18.5 kg/m2; 58% amenorrheic), 44 with atypical AN (18.5≤BMI<23 kg/m2; 30% amenorrheic), 45 eumenorrheic controls (19.2≤BMI<25 kg/m2). Measures: Int.vBMD and cross-sectional area (CSA) by quantitative computed tomography of L4; estimated vertebral strength (derived from Int.vBMD and CSA). Results: Int.vBMD and estimated vertebral strength were lowest in low-weight AN, intermediate in atypical AN, and highest in controls. CSA did not differ between groups; thus, vertebral strength (calculated using Int.vBMD and CSA) was driven by Int.vBMD. In AN, Int.vBMD and vertebral strength were associated positively with current BMI and nadir lifetime BMI (independent of current BMI). Int.vBMD and vertebral strength were lower in AN with current amenorrhea and longer lifetime amenorrhea duration. Among amenorrheic AN, Int.vBMD and vertebral strength were associated positively with testosterone. Conclusions: Int.vBMD and estimated vertebral strength (driven by Int.vBMD) are impaired across the AN weight spectrum and are associated with low BMI and endocrine dysfunction, both current and previous. Women with atypical AN experience diminished vertebral strength, partially due to prior low-weight and/or amenorrhea. Lack of current low-weight or amenorrhea in atypical AN does not preclude compromise of vertebral strength.

Evaluation of a new approach to compute intervertebral disc height measurements from lateral radiographic views of the spine.

PURPOSE: Current standard methods to quantify disc height, namely distortion compensated Roentgen analysis (DCRA), have been mostly utilized in the lumbar and cervical spine and have strict exclusion criteria. Specifically, discs adjacent to a vertebral fracture are excluded from measurement, thus limiting the use of DCRA in studies that include older populations with a high prevalence of vertebral fractures. Thus, we developed and tested a modified DCRA algorithm that does not depend on vertebral shape. METHODS: Participants included 1186 men and women from the Framingham Heart Study Offspring and Third Generation Multidetector CT Study. Lateral CT scout images were used to place 6 morphometry points around each vertebra at 13 vertebral levels in each participant. Disc heights were calculated utilizing these morphometry points using DCRA methodology and our modified version of DCRA, which requires information from fewer morphometry points than the standard DCRA. RESULTS: Modified DCRA and standard DCRA measures of disc height are highly correlated, with concordance correlation coefficients above 0.999. Both measures demonstrate good inter- and intra-operator reproducibility. 13.9 % of available disc heights were not evaluable or excluded using the standard DCRA algorithm, while only 3.3 % of disc heights were not evaluable using our modified DCRA algorithm. CONCLUSIONS: Using our modified DCRA algorithm, it is not necessary to exclude vertebrae with fracture or other deformity from disc height measurements as in the standard DCRA. Modified DCRA also yields identical measurements to the standard DCRA. Thus, the use of modified DCRA for quantitative assessment of disc height will lead to less missing data without any loss of accuracy, making it a preferred alternative to the current standard methodology.

Vertebral Strength and Estimated Fracture Risk Across the BMI Spectrum in Women.

Somewhat paradoxically, fracture risk, which depends on applied loads and bone strength, is elevated in both anorexia nervosa and obesity at certain skeletal sites. Factor-of-risk (Φ), the ratio of applied load to bone strength, is a biomechanically based method to estimate fracture risk; theoretically, higher Φ reflects increased fracture risk. We estimated vertebral strength (linear combination of integral volumetric bone mineral density [Int.vBMD] and cross-sectional area from quantitative computed tomography [QCT]), vertebral compressive loads, and Φ at L4 in 176 women (65 anorexia nervosa, 45 lean controls, and 66 obese). Using biomechanical models, applied loads were estimated for: 1) standing; 2) arms flexed 90°, holding 5 kg in each hand (holding); 3) 45° trunk flexion, 5 kg in each hand (lifting); 4) 20° trunk right lateral bend, 10 kg in right hand (bending). We also investigated associations of Int.vBMD and vertebral strength with lean mass (from dual-energy X-ray absorptiometry [DXA]) and visceral adipose tissue (VAT, from QCT). Women with anorexia nervosa had lower, whereas obese women had similar, Int.vBMD and estimated vertebral strength compared with controls. Vertebral loads were highest in obesity and lowest in anorexia nervosa for standing, holding, and lifting (p < 0.0001) but were highest in anorexia nervosa for bending (p < 0.02). Obese women had highest Φ for standing and lifting, whereas women with anorexia nervosa had highest Φ for bending (p < 0.0001). Obese and anorexia nervosa subjects had higher Φ for holding than controls (p < 0.03). Int.vBMD and estimated vertebral strength were associated positively with lean mass (R = 0.28 to 0.45, p ≤ 0.0001) in all groups combined and negatively with VAT (R = -[0.36 to 0.38], p < 0.003) within the obese group. Therefore, women with anorexia nervosa had higher estimated vertebral fracture risk (Φ) for holding and bending because of inferior vertebral strength. Despite similar vertebral strength as controls, obese women had higher vertebral fracture risk for standing, holding, and lifting because of higher applied loads from higher body weight. Examining the load-to-strength ratio helps explain increased fracture risk in both low-weight and obese women.

Incorporating Six Degree-of-Freedom Intervertebral Joint Stiffness in a Lumbar Spine Musculoskeletal Model-Method and Performance in Flexed Postures.

Intervertebral translations and rotations are likely dependent on intervertebral stiffness properties. The objective of this study was to incorporate realistic intervertebral stiffnesses in a musculoskeletal model of the lumbar spine using a novel force-dependent kinematics approach, and examine the effects on vertebral compressive loading and intervertebral motions. Predicted vertebral loading and intervertebral motions were compared to previously reported in vivo measurements. Intervertebral joint reaction forces and motions were strongly affected by flexion stiffness, as well as force-motion coupling of the intervertebral stiffness. Better understanding of intervertebral stiffness and force-motion coupling could improve musculoskeletal modeling, implant design, and surgical planning.

Development and Validation of a Musculoskeletal Model of the Fully Articulated Thoracolumbar Spine and Rib Cage.

We developed and validated a fully articulated model of the thoracolumbar spine in opensim that includes the individual vertebrae, ribs, and sternum. To ensure trunk muscles in the model accurately represent muscles in vivo, we used a novel approach to adjust muscle cross-sectional area (CSA) and position using computed tomography (CT) scans of the trunk sampled from a community-based cohort. Model predictions of vertebral compressive loading and trunk muscle tension were highly correlated to previous in vivo measures of intradiscal pressure (IDP), vertebral loading from telemeterized implants and trunk muscle myoelectric activity recorded by electromyography (EMG).

Vertebral size, bone density, and strength in men and women matched for age and areal spine BMD.

To explore the possible mechanisms underlying sex-specific differences in skeletal fragility that may be obscured by two-dimensional areal bone mineral density (aBMD) measures, we compared quantitative computed tomography (QCT)-based vertebral bone measures among pairs of men and women from the Framingham Heart Study Multidetector Computed Tomography Study who were matched for age and spine aBMD. Measurements included vertebral body cross-sectional area (CSA, cm(2) ), trabecular volumetric BMD (Tb.vBMD, g/cm(3) ), integral volumetric BMD (Int.vBMD, g/cm(3) ), estimated vertebral compressive loading and strength (Newtons) at L3 , the factor-of-risk (load-to-strength ratio), and vertebral fracture prevalence. We identified 981 male-female pairs (1:1 matching) matched on age (± 1 year) and QCT-derived aBMD of L3 (± 1%), with an average age of 51 years (range 34 to 81 years). Matched for aBMD and age, men had 20% larger vertebral CSA, lower Int.vBMD (-8%) and Tb.vBMD (-9%), 10% greater vertebral compressive strength, 24% greater vertebral compressive loading, and 12% greater factor-of-risk than women (p < 0.0001 for all), as well as higher prevalence of vertebral fracture. After adjusting for height and weight, the differences in CSA and volumetric bone mineral density (vBMD) between men and women were attenuated but remained significant, whereas compressive strength was no longer different. In conclusion, vertebral size, morphology, and density differ significantly between men and women matched for age and spine aBMD, suggesting that men and women attain the same aBMD by different mechanisms. These results provide novel information regarding sex-specific differences in mechanisms that underlie vertebral fragility.

Variations of CT-based trunk muscle attenuation by age, sex, and specific muscle.

BACKGROUND: Fat accumulation in muscle may contribute to age-related declines in muscle function and is indicated by reduced attenuation of x-rays by muscle tissue in computed tomography scans. Reduced trunk muscle attenuation is associated with poor physical function, low back pain, and increased hyperkyphosis in older adults. However, variations in trunk muscle attenuation with age, sex and between specific muscles have not been investigated. METHODS: A cross-sectional examination of trunk muscle attenuation in computed tomography scans was performed in 60 younger (35-50 years) and 60 older (75-87 years) adults randomly selected from participants in the Framingham Heart Study Offspring and Third Generation Multidetector Computed Tomography Study. Computed tomography attenuation of 11 trunk muscles was measured at vertebral levels T8 and L3, and the effects of age, sex, and specific muscle on computed tomography attenuation of trunk muscles were determined. RESULTS: Muscle attenuation varied by specific muscle (p < .001), was lower in older adults (p < .001), and was generally lower in women than in men (p < .001), although not in all muscles. Age-related differences in muscle attenuation varied with specific muscle (p < .001), with the largest age differences occurring in the paraspinal and abdominal muscles. CONCLUSIONS: Trunk muscle attenuation is lower in older adults than in younger adults in both women and men, but such age-related differences vary widely between muscle groups. The reasons that some muscles exhibit larger age-related differences in fat content than others should be further explored to better understand age-related changes in functional capacity and postural stability.

The effect of thoracic kyphosis and sagittal plane alignment on vertebral compressive loading.

To better understand the biomechanical mechanisms underlying the association between hyperkyphosis of the thoracic spine and risk of vertebral fracture and other degenerative spinal pathology, we used a previously validated musculoskeletal model of the spine to determine how thoracic kyphosis angle and spinal posture affect vertebral compressive loading. We simulated an age-related increase in thoracic kyphosis (T(1) -T(12) Cobb angle 50-75 degrees) during two different activities (relaxed standing and standing with 5-kg weights in the hands) and three different posture conditions: (1) an increase in thoracic kyphosis with no postural adjustment (uncompensated posture); (2) an increase in thoracic kyphosis with a concomitant increase in pelvic tilt that maintains a stable center of mass and horizontal eye gaze (compensated posture); and (3) an increase in thoracic kyphosis with a concomitant increase in lumbar lordosis that also maintains a stable center of mass and horizontal eye gaze (congruent posture). For all posture conditions, compressive loading increased with increasing thoracic kyphosis, with loading increasing more in the thoracolumbar and lumbar regions than in the mid-thoracic region. Loading increased the most for the uncompensated posture, followed by the compensated posture, with the congruent posture almost completely mitigating any increases in loading with increased thoracic kyphosis. These findings indicate that both thoracic kyphosis and spinal posture influence vertebral loading during daily activities, implying that thoracic kyphosis measurements alone are not sufficient to characterize the impact of spinal curvature on vertebral loading.

Regressions for estimating muscle parameters in the thoracic and lumbar trunk for use in musculoskeletal modeling.

Musculoskeletal modeling requires information on muscle parameters such as cross-sectional area (CSA) and moment arms. A variety of previous studies have reported muscle parameters in the trunk based on in vivo imaging, but there remain gaps in the available data as well as limitations in the generalizability of such data. Specifically, available trunk muscle CSA data is very limited for older adults, lacking entirely in the thoracic region. In addition, previous studies have made measurements in groups of healthy volunteers or hospital patients who may not be representative of the population in general. Finally, such studies have not reported data for the major muscles connecting the upper limb to the thoracic trunk. In this study, muscle morphology measurements were made for major muscles present in the trunk between vertebral levels T6 and L5 using quantitative computed tomography scans from a community-based sample of 100 men and women aged 36-87. We present regression equations to predict trunk muscle CSA and position relative to the vertebral body in the transverse plane from sex, age, height and weight at vertebral levels T6 to L5. Regressions were also developed for predicting anatomical CSA and muscle moment arms, which were estimated using literature data on muscle line of action. This work thus provides a resource for estimating muscle parameters in the general population for musculoskeletal modeling of the thoraco-lumbar trunk.