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, 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. This article is protected by copyright. All rights reserved