A Dynamic Radiographic Imaging Study of Lumbar Intervertebral Disc Morphometry and Deformation In Vivo

Cyclic Mechanism Affects Lumbar Spine Creep Response

PURPOSE

This study aims to explore how cyclic loading influences creep response in the lumbar spine under combined flexion-compression loading.

METHODS

Ten porcine functional spinal units (FSUs) were mechanically tested in cyclic or static combined flexion-compression loading. Creep response between loading regimes was compared using strain-time histories and linear regression. High-resolution computed tomography (µCT) visualized damage to FSUs. Statistical methods, ANCOVA and ANOVA, assessed differences in behavior between loading regimes.

RESULTS

Cyclic and static loading regimes exhibited distinct creep response patterns and biphasic response. ANCOVA and ANOVA analyses revealed significant differences in slopes of creep behavior in both linear phases. Cyclic tests consistently showed endplate fractures in µCT imaging.

CONCLUSION

The study reveals statistically significant differences in creep response between cyclic and static loading regimes in porcine lumbar spinal units under combined flexion-compression loading. The observed biphasic behavior suggests distinct phases of tissue response, indicating potential shifts in load transfer mechanisms. Endplate fractures in cyclic tests suggest increased injury risk compared to static loading. These findings underscore the importance of considering loading conditions in computational models and designing preventive measures for occupations involving repetitive spinal loading.

Regional structure-function relationships of lumbar cartilage endplates

Cartilage endplates (CEPs) act as protective mechanical barriers for intervertebral discs (IVDs), yet their heterogeneous structure-function relationships are poorly understood. This study addressed this gap by characterizing and correlating the regional biphasic mechanical properties and biochemical composition of human lumbar CEPs. Samples from central, lateral, anterior, and posterior portions of the disc (n = 8/region) were mechanically tested under confined compression to quantify swelling pressure, equilibrium aggregate modulus, and hydraulic permeability. These properties were correlated with CEP porosity and glycosaminoglycan (s-GAG) content, which were obtained by biochemical assays of the same specimens. Both swelling pressure (142.79 ± 85.89 kPa) and aggregate modulus (1864.10 ± 1240.99 kPa) were found to be regionally dependent (p = 0.0001 and p = 0.0067, respectively) in the CEP and trended lowest in the central location. No significant regional dependence was observed for CEP permeability (1.35 ± 0.97 * 10-16 m4/Ns). Porosity measurements correlated significantly with swelling pressure (r = -0.40, p = 0.0227), aggregate modulus (r = -0.49, p = 0.0046), and permeability (r = 0.36, p = 0.0421), and appeared to be the primary indicator of CEP biphasic mechanical properties. Second harmonic generation microscopy also revealed regional patterns of collagen fiber anchoring, with fibers inserting the CEP perpendicularly in the central region and at off-axial directions in peripheral regions. These results suggest that CEP tissue has regionally dependent mechanical properties which are likely due to the regional variation in porosity and matrix structure. This work advances our understanding of healthy baseline endplate biomechanics and lays a groundwork for further understanding the role of CEPs in IVD degeneration.

Intervertebral Disc Elastography to Relate Shear Modulus and Relaxometry in Compression and Bending

Intervertebral disc degeneration is the most recognized cause of low back pain, characterized by the decline of tissue structure and mechanics. Image-based mechanical parameters (e.g., strain, stiffness) may provide an ideal assessment of disc function that is lost with degeneration but unfortunately remains underdeveloped. Moreover, it is unknown whether strain or stiffness of the disc may be predicted by MRI relaxometry (e.g. T1 or T2), an increasingly accepted quantitative measure of disc structure. In this study, we quantified T1 and T2 relaxation times and in-plane strains using displacement-encoded MRI within the disc under physiological levels of compression and bending. We then estimated shear modulus in orthogonal image planes and compared these values to relaxation times and strains within regions of the disc. Intratissue strain depended on the loading mode, and shear modulus in the nucleus pulposus was typically an order of magnitude lower than the annulus fibrosis, except in bending, where the apparent stiffness depended on the loading. Relative shear moduli estimated from strain data derived under compression generally did not correspond with those from bending experiments, with no correlations in the sagittal plane and only 4 of 15 regions correlated in the coronal plane, suggesting that future inverse models should incorporate multiple loading conditions. Strain imaging and strain-based estimation of material properties may serve as imaging biomarkers to distinguish healthy and diseased discs. Additionally, image-based elastography and relaxometry may be viewed as complementary measures of disc structure and function to assess degeneration in longitudinal studies.

A novel tool to quantify in vivo lumbar spine kinematics and 3D intervertebral disc strains using clinical MRI

Medical imaging modalities that calculate tissue morphology alone cannot provide direct information regarding the mechanical behaviour of load-bearing musculoskeletal organs. Accurate in vivo measurement of spine kinematics and intervertebral disc (IVD) strains can provide important information regarding the mechanical behaviour of the spine, help to investigate the effects of injuries on the mechanics of the spine, and assess the effectiveness of treatments. Additionally, strains can serve as a functional biomechanical marker for detecting normal and pathologic tissues. We hypothesised that combining digital volume correlation (DVC) with 3T clinical MRI can provide direct information regarding the mechanics of the spine. Here, we have developed a novel non-invasive tool for in vivo displacement and strain measurement within the human lumbar spine and we used this tool to calculate lumbar kinematics and IVD strains in six healthy subjects during lumbar extension. The proposed tool enabled spine kinematics and IVD strains to be measured with errors that did not exceed 0.17 mm and 0.5%, respectively. The findings of the kinematics study identified that during extension the lumbar spine of healthy subjects experiences total 3D translations ranging from 1 mm to 4.5 mm for different vertebral levels. The findings of strain analysis identified that the average of the maximum tensile, compressive, and shear strains for different lumbar levels during extension ranged from 3.5% to 7.2%. This tool can provide base-line data that can be used to describe the mechanical environment of healthy lumbar spine, which can help clinicians manage preventative treatments, define patient-specific treatments, and to monitor the effectiveness of surgical and non-surgical interventions.

Non-Invasive Evaluation of Intradiscal Deformation during Axial Loading of the Spine Using Deformation-Field Magnetic Resonance Imaging: A Potential Tool for Micro-Instability Measurements

Degeneration alters the structural components of the disc and its mechanical behavior. Understanding this pathophysiological process is of great importance, as it may lead to back pain. However, non-invasive methods to characterize the disc mechanics in vivo are lacking. Here, a potential method for measurements of the intradiscal deformation under stress is presented. The method utilizes a standard MRI protocol, commercial loading equipment, and registration software. The lumbar spine (L1/L2–L5/S1) of 36 human subjects was imaged with and without axial loading of the spine. The resulting images were registered, and changes in the images during the registration were displayed pixel-by-pixel to visualize the internal deformation of the disc. The degeneration grade, disc height, disc angle and tilt angle were determined and correlated with the deformation using multivariate regression analysis. The largest deformation was found at the lower lumbar spine, and differences in regional behaviors between individual discs were found. Weak to moderate correlations between the deformation and different disc characteristics were found, where the degeneration grade and tilt angle were the main contributing factors. To conclude, the image-based method offers a potential tool to study the pathophysiological process of the disc.

Low-cost luminescent scaffolds-based on thiol chitosans by microwave radiation for vertebral disc repair/theragnostic

This study introduces a new 3D scaffold based on thiolated chitosans with luminescence by microwave radiation using cysteine (Chi_CT_Cys) and 11-mercaptoundecanoic acid (Chi_CT_MUA) for vertebral disc regeneration/theragnostic. These scaffolds were characterized by Raman, PL spectroscopy, swelling, gel-fraction, and morphologies. Cytocompatibility and mechanical behavior were evaluated. Raman showed that disulfide bonds improved the grafting degree (Chi_CT_Cys (1072 ± 136) μmol·g^-1 and Chi_CT_MUA (3245 ± 105) μmol·g^-1). Morphologies showed interesting characteristics. Swelling behavior showed that Chi_CT_MUA presented a slight minor swelling (2101 ± 251) % compared to Chi_CT_Cys (2589 ± 188) %. Differently, gel-fraction showed that the chemical stability of Chi_CT_Cys was worse (29 ± 4) % than Chi_CT_MUA (15 ± 3) %. PL showed a possibility to use theragnostic evaluation of points of greater compression in a vertebral disc. The mechanical behavior of Chi_CT_MUA presented better results ((70 ± 3) MPa) than Chi_CT_Cys ((37 ± 3) MPa). Cytocompatible showed that the scaffolds presented cell viability >90%. Thusly, these 3D scaffolds presented an incredible potential for tissue engineering applications.

Investigation of geometric deformations of the lumbar disc during axial body rotations

Background

Quantitative data on in vivo vertebral disc deformations are critical for enhancing our understanding of spinal pathology and improving the design of surgical materials. This study investigated in vivo lumbar intervertebral disc deformations during axial rotations under different load-bearing conditions.

Methods

Twelve healthy subjects (7 males and 5 females) between the ages of 25 and 39 were recruited. Using a combination of a dual fluoroscopic imaging system (DFIS) and CT, the images of L3–5 segments scanned by CT were transformed into three-dimensional models, which matched the instantaneous images of the lumbar spine taken by a double fluorescent X-ray system during axial rotations to reproduce motions. Then, the kinematic data of the compression and shear deformations of the lumbar disc and the coupled bending of the vertebral body were obtained.

Results

Relative to the supine position, the average compression deformation caused by rotation is between + 10% and − 40%, and the shear deformation is between 17 and 50%. Under physiological weightbearing loads, different levels of lumbar discs exhibit similar deformation patterns, and the deformation patterns of left and right rotations are approximately symmetrical. The deformation patterns change significantly under a 10 kg load, with the exception of the L3–4 disc during the right rotation.

Conclusion

The deformation of the lumbar disc was direction-specific and level-specific during axial rotations and was affected by extra weight. These data can provide new insights into the biomechanics of the lumbar spine and optimize the parameters of artificial lumbar spine devices.

ISSLS Prize in Bioengineering Science 2021: in vivo sagittal motion of the lumbar spine in low back pain patients-a radiological big data study

PURPOSE

We investigated the flexion-extension range of motion and centre of rotation of lumbar motion segments in a large population of 602 patients (3612 levels), and the associations between lumbar motion and other parameters such as sex, age and intervertebral disc degeneration.

METHODS

Lumbar radiographs in flexion-extension of 602 patients suffering from low back pain and/or suspect instability were collected; magnetic resonance images were retrieved and used to score the degree of disc degeneration for a subgroup of 354 patients. Range of motion and centre of rotation were calculated for all lumbosacral levels with in-house software allowing for high degree of automation. Associations between motion parameters and age, sex, spinal level and disc degeneration were then assessed.

RESULTS

The median range of motion was 6.6° (range 0.1-28.9°). Associations between range of motion and age as well as spinal level, but not sex, were found. Disc degeneration determined a consistent reduction in the range of motion. The centre of rotation was most commonly located at the centre of the lower endplate or slightly lower. With progressive degeneration, centres of rotation were increasingly dispersed with no preferential directions.

CONCLUSION

This study constitutes the largest analysis of the in vivo lumbar motion currently available and covers a wide range of clinical scenarios in terms of age and degeneration. Findings confirmed that ageing determines a reduction in the mobility independently of degeneration and that in degenerative levels, centres of rotation are dispersed around the centre of the intervertebral space.

Lumbar Intervertebral Kinematics During an Unstable Sitting Task and Its Association With Standing-Induced Low Back Pain

People developing transient low back pain during standing have altered control of their spine and hips during standing tasks, but the transfer of these responses to other tasks has not been assessed. This study used video fluoroscopy to assess lumbar spine intervertebral kinematics of people who do and do not develop standing-induced low back pain during a seated chair-tilting task. A total of 9 females and 8 males were categorized as pain developers (5 females and 3 males) or nonpain developers (4 females and 5 males) using a 2-hour standing exposure; pain developers reported transient low back pain and nonpain developers did not. Participants were imaged with sagittal plane fluoroscopy at 25 Hz while cyclically tilting their pelvises anteriorly and posteriorly on an unstable chair. Intervertebral angles, relative contributions, and anterior-posterior translations were measured for the L3/L4, L4/L5, and L5/S1 joints and compared between sexes, pain groups, joints, and tilting directions. Female pain developers experienced more extension in their L5/S1 joints in both tilting directions compared with female nonpain developers, a finding not present in males. The specificity in intervertebral kinematics to sex-pain group combinations suggests that these subgroups of pain developers and nonpain developers may implement different control strategies.