Kinematics of the Spine Under Healthy and Degenerative Conditions: A Systematic Review

Towards Wearable and Portable Spine Motion Analysis Through Dynamic Optimization of Smartphone Videos and IMU Data

BACKGROUND

Monitoring spine kinematics is crucial for applications like disease evaluation and ergonomics analysis. However, the small scale of vertebrae and the number of degrees of freedom present significant challenges for noninvasive and convenient spine kinematics estimation.

METHODS

This study developed a dynamic optimization framework for wearable spine motion tracking at the intervertebral joint level by integrating smartphone videos and Inertia Measurement Units (IMUs) with dynamic constraints from a thoracolumbar spine model. Validation involved motion data from 10 healthy males performing static standing, dynamic upright trunk rotations, and gait. This data included rotations of ten IMUs on vertebrae and virtual landmarks from three smartphone videos preprocessed by OpenCap, an application leveraging computer vision for pose estimation. The kinematic measures derived from the optimized solution were compared against simultaneously collected infrared optical marker-based measurements and in vivo literature data. Solutions only based on IMUs or videos were also compared for accuracy evaluation.

RESULTS

The proposed optimization approach closely matched the reference data in the intervertebral or segmental rotation range, demonstrating minimal angular differences across all motions and the highest correlation in 3D rotations (maximal Pearson and intraclass correlation coefficients of 0.92 and 0.94, respectively). Time-series changes of joint angles also aligned well with the optical-marker reference.

CONCLUSION

Dynamic optimization of the spine simulation that integrates IMUs and computer vision outperforms the single-modality method.

SIGNIFICANCE

This markerless 3D spine motion capture method holds potential for spinal health assessment in large cohorts in real-world settings without dedicated laboratories.

The relationship between paraspinal muscle atrophy and degenerative lumbar spondylolisthesis at the L4/5 level

BACKGROUND/CONTEXT

Degenerative lumbar spondylolisthesis (DLS) is a prevalent spinal condition that can result in significant disability. DLS is thought to result from a combination of disc and facet joint degeneration, as well as various biological, biomechanical, and behavioral factors. One hypothesis is the progressive degeneration of segmental stabilizers, notably the paraspinal muscles, contributes to a vicious cycle of increasing slippage.

PURPOSE

To examine the correlation between paraspinal muscle status on MRI and severity of slippage in patients with symptomatic DLS.

STUDY DESIGN/SETTING

Retrospective cross-sectional study at an academic tertiary care center.

PATIENT SAMPLE

Patients who underwent surgery for DLS at the L4/5 level between 2016-2018 were included. Those with multilevel DLS or insufficient imaging were excluded.

OUTCOME MEASURES

The percentage of relative slippage (RS) at the L4/5 level evaluated on standing lateral radiographs. Muscle morphology measurements including functional cross-sectional area (fCSA), body height normalized functional cross-sectional area (HI) of Psoas, erector spinae (ES) and multifidus muscle (MF) and fatty infiltration (FI) of ES and MF were measured on axial MR. Disc degeneration and facet joint arthritis were classified according to Pfirrmann and Weishaupt, respectively.

METHODS

Descriptive and comparative statistics, univariable and multivariable linear regression models were utilized to examine the associations between RS and muscle parameters, adjusting for confounders sex, age, BMI, segmental degeneration, and back pain severity and symptom duration.

RESULTS

The study analyzed 138 out of 183 patients screened for eligibility. The median age of all patients was 69.5 years (IQR 62 to 73), average BMI was 29.1 (SD±5.1) and average preoperative ODI was 46.4 (SD±16.3). Patients with Meyerding-Grade 2 (M2, N=25) exhibited higher Pfirrmann scores, lower MFfCSA and MFHI, and lower BMI, but significantly more fatty infiltration in the MF and ES muscles compared to those with Meyerding Grade 1 (M1). Univariable linear regression showed that each cm2 decrease in MFfCSA was associated with a 0.9%-point increase in RS (95% CI -1.4 to - 0.4, p<.001), and each cm2/m2 decrease in MFHI was associated with an increase in slippage by 2.2%-points (95% CI -3.7 to -0.7, p=.004). Each 1%-point rise in ESFI and MFFI corresponded to 0.17%- (95% CI 0.05-0.3, p=.01) and 0.20%-point (95% CI 0.1-0.3 p<.001) increases in relative slippage, respectively. Notably, after adjusting for confounders, each cm2 increase in PsoasfCSA and cm2/m2 in PsoasHI was associated with an increase in relative slippage by 0.3% (95% CI 0.1-0.6, p=.004) and 1.1%-points (95% CI 0.4-1.7, p=.001). While MFfCSA tended to be negatively associated with slippage, this did not reach statistical significance (p=.105). However, each 1%-point increase in MFFI and ESFI corresponded to increases of 0.15% points (95% CI 0.05-0.24, p=.002) and 0.14% points (95% CI 0.01-0.27, p=.03) in relative slippage, respectively.

CONCLUSION

This study found a significant association between paraspinal muscle status and severity of slippage in DLS. Whereas higher degeneration of the ES and MF correlate with a higher degree of slippage, the opposite was found for the psoas. These findings suggest that progressive muscular imbalance between posterior and anterior paraspinal muscles could contribute to the progression of slippage in DLS.

Effect of patient-specific factors on regeneration in lumbar spine at healthy disc and total disc replacement. Computer simulation

BACKGROUND AND OBJECTIVE

Degenerative diseases of the spine have a negative impact on the quality of life of patients. This study presents the results of numerical modelling of the mechanical behaviour of the lumbar spine with patient-specific conditions at physiological loads. This paper aims to numerically study the influence of degenerative changes in the spine and the presence of an endoprosthesis on the creation of conditions for tissue regeneration.

METHODS

A numerical model of the mechanical behaviour of lumbar spine at healthy and after total disc replacement under low-energy impacts equivalent to physiological loads is presented. The model is based on the movable cellular automaton method (discrete elements), where the mechanical behaviour of bone tissue is described using the Biot poroelasticity accounting for the presence and transfer of interstitial biological fluid. The nutritional pathways of the intervertebral disc in cases of healthy and osteoporotic bone tissues were predicted based on the analysis of the simulation results according to the mechanobiological principles.

RESULTS

Simulation of total disc replacement showed that osseointegration of the artificial disc plates occurs only in healthy bone tissue. With total disc replacement in a patient with osteoporosis, there is an area of increased risk of bone resorption in the near-contact area, approximately 1 mm wide, around the fixators. Dynamic loads may improve the osseointegration of the implant in pathological conditions of the bone tissue.

CONCLUSIONS

The results obtained in the case of healthy spine and osteoporotic bone tissues correspond to the experimental data on biomechanics and possible methods of IVD regeneration from the position of mechanobiological principles. The results obtained with an artificial disc (with keel-type fixation) showed that the use of this type of endoprosthesis in healthy bone tissues allows to reproduce the function of the natural intervertebral disc and does not contribute to the development of neoplastic processes. In the case of an artificial disc with osteoporosis of bone tissues, there is a zone with increased risk of tissue resorption and development of neoplastic processes in the area near the contact of the implant attachment. This circumstance can be compensated by increasing the loading level.

Dynamic segmental kinematics of the lumbar spine during diagnostic movements

Background: In vivo measurements of segmental-level kinematics are a promising avenue for better understanding the relationship between pain and its underlying, multi-factorial basis. To date, the bulk of the reported segmental-level motion has been restricted to single plane motions. Methods: The present work implemented a novel marker set used with an optical motion capture system to non-invasively measure dynamic, 3D in vivo segmental kinematics of the lower spine in a laboratory setting. Lumbar spinal kinematics were measured for 28 subjects during 17 diagnostic movements. Results: Overall regional range of motion data and lumbar angular velocity measurement were consistent with previously published studies. Key findings from the work included measurement of differences in ascending versus descending segmental velocities during functional movements and observations of motion coupling paradigms in the lumbar spinal segments. Conclusion: The work contributes to the task of establishing a baseline of segmental lumbar movement patterns in an asymptomatic cohort, which serves as a necessary pre-requisite for identifying pathological and symptomatic deviations from the baseline.

Biomechanics after spinal decompression and posterior instrumentation

PURPOSE

The aim of this study was to elucidate segmental range of motion (ROM) before and after common decompression and fusion procedures on the lumbar spine.

METHODS

ROM of fourteen fresh-frozen human cadaver lumbar segments (L1/2: 4, L3/4: 5, L5/S1: 5) was evaluated in six loading directions: flexion/extension (FE), lateral bending (LB), lateral shear (LS), anterior shear (AS), axial rotation (AR), and axial compression/distraction (AC). ROM was tested with and without posterior instrumentation under the following conditions: 1) native 2) after unilateral laminotomy, 3) after midline decompression, and 4) after nucleotomy.

RESULTS

Median native ROM was FE 6.8°, LB 5.6°, and AR 1.7°, AS 1.8 mm, LS 1.4 mm, AC 0.3 mm. Unilateral laminotomy significantly increased ROM by 6% (FE), 3% (LB), 12% (AR), 11% (AS), and 8% (LS). Midline decompression significantly increased these numbers to 15%, 5%, 21%, 20%, and 19%, respectively. Nucleotomy further increased ROM in all directions, most substantially in AC of 153%. Pedicle screw fixation led to ROM decreases of 82% in FE, 72% in LB, 42% in AR, 31% in AS, and 17% in LS. In instrumented segments, decompression only irrelevantly affected ROM.

CONCLUSIONS

The amount of posterior decompression significantly impacts ROM of the lumbar spine. The here performed biomechanical study allows creation of a simplified rule of thumb: Increases in segmental ROM of approximately 10%, 20%, and 50% can be expected after unilateral laminotomy, midline decompression, and nucleotomy, respectively. Instrumentation decreases ROM by approximately 80% in bending moments and accompanied decompression procedures only minorly destabilize the instrumentation construct.

Spatio-temporal clustering of lumbar intervertebral flexion interactions in 127 asymptomatic individuals

The purpose of this study was to categorize asymptomatic participants based on the clustering of spatial and temporal intervertebral kinematic variables during lumbar flexion. Lumbar segmental interactions (L2-S1) were evaluated in 127 asymptomatic participants during flexion using fluoroscopy. First, four variables were identified consisting of: 1. Range of motion (ROMC), 2. Peaking time of the first derivative for separate segmentation (PTFDs), 3. Peaking magnitude of the first derivative (PMFD), and 4. Peaking time of the first derivative for stepwise (grouped) segmentation (PTFDss). These variables were used to cluster and order the lumbar levels. The number of participants required to constitute a cluster was chosen as 7. Participants formed eight (ROMC), four (PTFDs), eight (PMFD), and four (PTFDss) clusters, which included 85%, 80%, 77%, and 60% of them, respectively, according to the above features. For all clustering variables, angle time series of some lumbar levels showed significant differences between clusters. However, in general, all clusters could be categorized based on the segmental mobility contexts into three main groups as incidental macro clusters: the upper (L2-L4 > L4-S1), middle (L2-L3 < L3-L5 > L5-S1) and lower (L2-L4 < L4-S1) domains. There are spatial and temporal segmental interactions and between-subject variability in asymptomatic participants. In addition, the differences in angle time series among the clusters have provided evidence of feedback control strategies, while the stepwise segmentation facilitates consideration of the lumbar spine as a system and provides supplementary information about segmental interactions. Clinically, these facts could be taken into account when considering any intervention, but especially fusion surgery.

Facet Joint Orientation/Tropism Could Be Associated with Fatty Infiltration in the Lumbar Paraspinal Muscles

BACKGROUND

Facet joint orientation (FJO) and facet joint tropism (FJT) are associated with intervertebral disc degeneration and paraspinal muscle atrophy. However, none of the previous studies has evaluated the association of FJO/FJT with fatty infiltration in the multifidus, erector spinae, and psoas muscles at all lumbar levels. In the present study, we aimed to analyze whether FJO and FJT were associated with fatty infiltration in the paraspinal muscles at any lumbar level.

METHODS

Paraspinal muscles and FJO/FJT were evaluated from L1-L2 to L5-S1 intervertebral disc levels on T2-weighted axial lumbar spine magnetic resonance imaging.

RESULTS

Facet joints were more sagittally and coronally oriented at the upper and lower lumbar levels, respectively. FJT was more obvious at lower lumbar levels. The FJT/FJO ratio was higher at upper lumbar levels. Patients with sagittally oriented facet joints at the L3-L4 and L4-L5 levels had fattier erector spinae and psoas muscles at the L4-L5 level. Patients with increased FJT at upper lumbar levels had fattier erector spinae and multifidus at lower lumbar levels. Patients with increased FJT at the L4-L5 level had less fatty infiltration in the erector spinae and psoas at the L2-L3 and L5-S1 levels, respectively.

CONCLUSIONS

Sagittally oriented facet joints at lower lumbar levels could be associated with fattier erector spinae and psoas muscles at lower lumbar levels. The erector spinae at upper lumbar levels and psoas at lower lumbar levels might have become more active to compensate the FJT-induced instability at lower lumbar levels.

Novel force-displacement control passive finite element models of the spine to simulate intact and pathological conditions; comparisons with traditional passive and detailed musculoskeletal models

Passive finite element (FE) models of the spine are commonly used to simulate intact and various pre- and postoperative pathological conditions. Being devoid of muscles, these traditional models are driven by simplistic loading scenarios, e.g., a constant moment and compressive follower load (FL) that do not properly mimic the complex in vivo loading condition under muscle exertions. We aim to develop novel passive FE models that are driven by more realistic yet simple loading scenarios, i.e., in vivo vertebral rotations and pathological-condition dependent FLs (estimated based on detailed musculoskeletal finite element (MS-FE) models). In these novel force-displacement control FE models, unlike the traditional passive FE models, FLs vary not only at different spine segments (T12-S1) but between intact, pre- and postoperative conditions. Intact, preoperative degenerated, and postoperative fused conditions at the L4-L5 segment for five static in vivo activities in upright and flexed postures were simulated by the traditional passive FE, novel force-displacement control FE, and gold-standard detailed MS-FE spine models. Our findings indicate that, when compared to the MS-FE models, the force-displacement control passive FE models could accurately predict the magnitude of disc compression force, intradiscal pressure, annulus maximal von Mises stress, and vector sum of all ligament forces at adjacent segments (L3-L4 and L5-S1) but failed to predict disc shear and facet joint forces. In this regard, the force-displacement control passive FE models were much more accurate than the traditional passive FE models. Clinical recommendations made based on traditional passive FE models should, therefore, be interpreted with caution.

Analysis of the risk factors for early tether breakage following vertebral body tethering in adolescent idiopathic scoliosis

INTRODUCTION

Tether breakage is a common mechanical complication after VBT. When this occurs shortly after surgery, patients may be at higher risk for loss of correction. Aim of this study was to analyze demographic and radiographic parameters that may potentially be risk factors for early tether breakage, as no data are yet available on this topic.

MATERIALS AND METHODS

All skeletally immature patients who underwent VBT and for whom a 1-year follow-up was available were included in the study. Demographic, intraoperative and coronal and sagittal parameters from the preoperative and 1st standing X-rays were collected. Patients were divided in two groups according to the presence or absence of a breakage and the outcomes of interest were compared.

RESULTS

Data from 105 patients were available (age 14.2 ± 1.5, 153 curves). Lumbar curves showed a higher risk of breakage than thoracic ones (71% vs. 29%, P < 0.0001). Overall, preoperative risk factors were a high curve magnitude (MD, mean difference - 4.1°, P = 0.03) and a limited flexibility (MD 8.9%, P = 0.006); postoperative risk factors were a large residual curve (MD - 6.4°, P = 0.0005) and a limited correction (MD 8.4%, P = 0.0005). The same risk factors were identified in thoracic curves, while in lumbar instrumentation only a higher preoperative Cobb angle represented a risk factor for breakage. Age and skeletal maturity did not represent risk factors.

CONCLUSION

The main preoperative risk factors for early tether breakage after VBT are a high curve magnitude and a limited flexibility. A limited curve correction also represents a risk factor for this complication.

Computational model predicts risk of spinal screw loosening in patients

PURPOSE

Pedicle screw loosening is a frequent complication in lumbar spine fixation, most commonly among patients with poor bone quality. Determining patients at high risk for insufficient implant stability would allow clinicians to adapt the treatment accordingly. The aim of this study was to develop a computational model for quantitative and reliable assessment of the risk of screw loosening.

METHODS

A cohort of patient vertebrae with diagnosed screw loosening was juxtaposed to a control group with stable fusion. Imaging data from the two cohorts were used to generate patient-specific biomechanical models of lumbar instrumented vertebral bodies. Single-level finite element models loading the screw in axial or caudo-cranial direction were generated. Further, multi-level models incorporating individualized joint loading were created.

RESULTS

The simulation results indicate that there is no association between screw pull-out strength and the manifestation of implant loosening (p = 0.8). For patient models incorporating multiple instrumented vertebrae, CT-values and stress in the bone were significantly different between loose screws and non-loose screws (p = 0.017 and p = 0.029, for CT-values and stress, respectively). However, very high distinction (p = 0.001) and predictability (R²_Pseudo = 0.358, AUC = 0.85) were achieved when considering the relationship between local bone strength and the predicted stress (loading factor). Screws surrounded by bone with a loading factor higher than 25% were likely to be loose, while the chances of screw loosening were close to 0 with a loading factor below 15%.

CONCLUSION

The use of a biomechanics-based score for risk assessment of implant fixation failure might represent a paradigm shift in addressing screw loosening after spondylodesis surgery.

Force Distribution Within Spinal Tissues During Posterior to Anterior Spinal Manipulative Therapy: A Secondary Analysis

Background

Previous studies observed that the intervertebral disc experiences the greatest forces during spinal manipulative therapy (SMT) and that the distribution of forces among spinal tissues changes as a function of the SMT parameters. However, contextualized SMT forces, relative to the ones applied to and experienced by the whole functional spinal unit, is needed to understand SMT’s underlying mechanisms.

Aim

To describe the percentage force distribution between spinal tissues relative to the applied SMT forces and total force experienced by the functional unit.

Methods

This secondary analysis combined data from 35 fresh porcine cadavers exposed to a simulated 300N SMT to the skin overlying the L3/L4 facet joint via servo-controlled linear motor actuator. Vertebral kinematics were tracked optically using indwelling bone pins. The functional spinal unit was then removed and mounted on a parallel robotic platform equipped with a 6-axis load cell. The kinematics of the spine during SMT were replayed by the robotic platform. By using serial dissection, peak and mean forces induced by the simulated SMT experienced by spinal structures in all three axes of motion were recorded. Forces experienced by spinal structures were analyzed descriptively and the resultant force magnitude was calculated.

Results

During SMT, the functional spinal unit experienced a median peak resultant force of 36.4N (IQR: 14.1N) and a mean resultant force of 25.4N (IQR: 11.9N). Peak resultant force experienced by the spinal segment corresponded to 12.1% of the total applied SMT force (300N). When the resultant force experienced by the functional spinal unit was considered to be 100%, the supra and interspinous ligaments experienced 0.3% of the peak forces and 0.5% of the mean forces. Facet joints and ligamentum flavum experienced 0.7% of the peak forces and 3% of the mean forces. Intervertebral disc and longitudinal ligaments experienced 99% of the peak and 96.5% of the mean forces.

Conclusion

In this animal model, a small percentage of the forces applied during a posterior-to-anterior SMT reached spinal structures in the lumbar spine. Most SMT forces (over 96%) are experienced by the intervertebral disc. This study provides a novel perspective on SMT force distribution within spinal tissues.

A novel approach for tetrahedral-element-based finite element simulations of anisotropic hyperelastic intervertebral disc behavior

Intervertebral discs are microstructurally complex spinal tissues that add greatly to the flexibility and mechanical strength of the human spine. Attempting to provide an adjustable basis for capturing a wide range of mechanical characteristics and to better address known challenges of numerical modeling of the disc, we present a robust finite-element-based model formulation for spinal segments in a hyperelastic framework using tetrahedral elements. We evaluate the model stability and accuracy using numerical simulations, with particular attention to the degenerated intervertebral discs and their likely skewed and narrowed geometry. To this end, 1) annulus fibrosus is modeled as a fiber-reinforced Mooney-Rivlin type solid for numerical analysis. 2) An adaptive state-variable dependent explicit time step is proposed and utilized here as a computationally efficient alternative to theoretical estimates. 3) Tetrahedral-element-based FE models for spinal segments under various loading conditions are evaluated for their use in robust numerical simulations. For flexion, extension, lateral bending, and axial rotation load cases, numerical simulations reveal that a suitable framework based on tetrahedral elements can provide greater stability and flexibility concerning geometrical meshing over commonly employed hexahedral-element-based ones for representation and study of spinal segments in various stages of degeneration.

Kinematics of the Cervical Spine Under Healthy and Degenerative Conditions: A Systematic Review

Knowledge of spinal kinematics is essential for the diagnosis and management of spinal diseases. Distinguishing between physiological and pathological motion patterns can help diagnose these diseases, plan surgical interventions and improve relevant tools and software. During the last decades, numerous studies based on diverse methodologies attempted to elucidate spinal mobility in different planes of motion. The authors aimed to summarize and compare the evidence about cervical spine kinematics under healthy and degenerative conditions. This includes an illustrated description of the spectrum of physiological cervical spine kinematics, followed by a comparable presentation of kinematics of the degenerative cervical spine. Data was obtained through a systematic MEDLINE search including studies on angular/translational segmental motion contribution, range of motion, coupling and center of rotation. As far as the degenerative conditions are concerned, kinematic data regarding disc degeneration and spondylolisthesis were available. Although the majority of the studies identified repeating motion patterns for most motion planes, discrepancies associated with limited sample sizes and different imaging techniques and/or spine configurations, were noted. Among healthy/asymptomatic individuals, flexion extension (FE) and lateral bending (LB) are mainly facilitated by the subaxial cervical spine. C4-C5 and C5-C6 were the major FE contributors in the reported studies, exceeding the motion contribution of sub-adjacent segments. Axial rotation (AR) greatly depends on C1-C2. FE range of motion (ROM) is distributed between the atlantoaxial and subaxial segments, while AR ROM stems mainly from the former and LB ROM from the latter. In coupled motion rotation is quantitatively predominant over translation. Motion migrates caudally from C1-C2 and the center of rotation (COR) translocates anteriorly and superiorly for each successive subaxial segment. In degenerative settings, concurrent or subsequent lesions render the association between diseases and mobility alterations challenging. The affected segments seem to maintain translational and angular motion in early and moderate degeneration. However, the progression of degeneration restrains mobility, which seems to be maintained or compensated by adjacent non-affected segments. While the kinematics of the healthy cervical spine have been addressed by multiple studies, the entire nosological and kinematic spectrum of cervical spine degeneration is partially addressed. Large-scale in vivo studies can complement the existing evidence, cover the gaps and pave the way to technological and clinical breakthroughs.

Posterior spinal instrumentation and decompression with or without cross-link?

Background

Posterior lumbar instrumentation requires sufficient primary stiffness to ensure bony fusion and to avoid pseudarthrosis, screw loosening, or implant failure. To enhance primary construct stiffness, transverse cross-link (CL) connectors attached to the vertical rods can be used. Their effect on the stability of a spinal instrumentation with simultaneous decompression is yet not clear. This study aimed to evaluate the impact of CL augmentation on single-level lumbar instrumentation stiffness after gradual decompression procedures.

Methods

Seventeen vertebral segments (6 L1/2, 6 L3/4, 5 L5/S1) of 12 fresh-frozen human cadavers were instrumented with a transpedicular screw–rod construct following the traditional pedicle screw trajectory. Range of motion (ROM) of the segments was sequentially recorded before and after four procedures: (A) instrumented before decompression, (B) instrumented after unilateral laminotomy, (C) instrumented after midline bilateral laminotomy, and (D) instrumented after unilateral facetectomy (with transforaminal lumbar interbody fusion [TLIF]). Each test was performed with and without CL augmentation. The motion between the cranial and caudal vertebrae was evaluated in all six major loading directions: flexion/extension (FE), lateral bending (LB), lateral shear (LS), anterior shear (AS), axial rotation (AR), and axial compression/distraction (AC).

Results

ROM was significantly reduced with CL augmentation in AR by Δ0.03–0.18° (7–12%) with a significantly higher ROM reduction after more extensive decompression. Furthermore, slight reductions in FE and LB were observed; these reached statistical significance for FE after facetectomy and TLIF insertion only (Δ0.15; 3%). The instrumentation levels did not reveal any subgroup differences.

Conclusion

CL augmentation reduces AR-ROM by 7–12% in single-level instrumentation of the lumbar spine, with the effect increasing along with the extensiveness of the decompression technique. In light of the discrete absolute changes, CL augmentation may be warranted for highly unstable vertebral segments rather than for standard single-level posterior spinal fusion and decompression.

A Reference Database of Standardised Continuous Lumbar Intervertebral Motion Analysis for Conducting Patient-Specific Comparisons

Lumbar instability has long been thought of as the failure of lumbar vertebrae to maintain their normal patterns of displacement. However, it is unknown what these patterns consist of. Research using quantitative fluoroscopy (QF) has shown that continuous lumbar intervertebral patterns of rotational displacement can be reliably measured during standing flexion and return motion using standardised protocols and can be used to assess patients with suspected lumbar spine motion disorders. However, normative values are needed to make individualised comparisons. One hundred and thirty-one healthy asymptomatic participants were recruited and performed guided flexion and return motion by following the rotating arm of an upright motion frame. Fluoroscopic image acquisition at 15fps was performed and individual intervertebral levels from L2-3 to L5-S1 were tracked and analysed during separate outward flexion and return phases. Results were presented as proportional intervertebral motion representing these phases using continuous means and 95%CIs, followed by verification of the differences between levels using Statistical Parametric Mapping (SPM). A secondary analysis of 8 control participants matched to 8 patients with chronic, non-specific low back pain (CNSLBP) was performed for comparison. One hundred and twenty-seven asymptomatic participants’ data were analysed. Their ages ranged from 18 to 70 years (mean 38.6) with mean body mass index 23.8 kg/m² 48.8% were female. Both the flexion and return phases for each level evidenced continuous change in mean proportional motion share, with narrow confidence intervals, highly significant differences and discrete motion paths between levels as confirmed by SPM. Patients in the secondary analysis evidenced significantly less L5-S1 motion than controls (p < 0.05). A reference database of spinal displacement patterns during lumbar (L2-S1) intersegmental flexion and return motion using a standardised motion protocol using fluoroscopy is presented. Spinal displacement patterns in asymptomatic individuals were found to be distinctive and consistent for each intervertebral level, and to continuously change during bending and return. This database may be used to allow continuous intervertebral kinematics to drive dynamic models of joint and muscular forces as well as reference values against which to make patient-specific comparisons in suspected cases of lumbar spine motion disorders.

Subject-Specific Alignment and Mass Distribution in Musculoskeletal Models of the Lumbar Spine

Musculoskeletal modeling is a well-established method in spine biomechanics and generally employed for investigations concerning both the healthy and the pathological spine. It commonly involves inverse kinematics and optimization of muscle activity and provides detailed insight into joint loading. The aim of the present work was to develop and validate a procedure for the automatized generation of semi-subject-specific multi-rigid body models with an articulated lumbar spine. Individualization of the models was achieved with a novel approach incorporating information from annotated EOS images. The size and alignment of bony structures, as well as specific body weight distribution along the spine segments, were accurately reproduced in the 3D models. To ensure the pipeline’s robustness, models based on 145 EOS images of subjects with various weight distributions and spinopelvic parameters were generated. For validation, we performed kinematics-dependent and segment-dependent comparisons of the average joint loads obtained for our cohort with the outcome of various published in vivo and in situ studies. Overall, our results agreed well with literature data. The here described method is a promising tool for studying a variety of clinical questions, ranging from the evaluation of the effects of alignment variation on joint loading to the assessment of possible pathomechanisms involved in adjacent segment disease.

Biomechanical effects of lumbar fusion surgery on adjacent segments using musculoskeletal models of the intact, degenerated and fused spine

Adjacent segment disorders are prevalent in patients following a spinal fusion surgery. Postoperative alterations in the adjacent segment biomechanics play a role in the etiology of these conditions. While experimental approaches fail to directly quantify spinal loads, previous modeling studies have numerous shortcomings when simulating the complex structures of the spine and the pre/postoperative mechanobiology of the patient. The biomechanical effects of the L4–L5 fusion surgery on muscle forces and adjacent segment kinetics (compression, shear, and moment) were investigated using a validated musculoskeletal model. The model was driven by in vivo kinematics for both preoperative (intact or severely degenerated L4–L5) and postoperative conditions while accounting for muscle atrophies. Results indicated marked changes in the kinetics of adjacent L3–L4 and L5–S1 segments (e.g., by up to 115% and 73% in shear loads and passive moments, respectively) that depended on the preoperative L4–L5 disc condition, postoperative lumbopelvic kinematics and, to a lesser extent, postoperative changes in the L4–L5 segmental lordosis and muscle injuries. Upper adjacent segment was more affected post-fusion than the lower one. While these findings identify risk factors for adjacent segment disorders, they indicate that surgical and postoperative rehabilitation interventions should focus on the preservation/restoration of patient’s normal segmental kinematics.

Biomechanical Contributions of Spinal Structures with Different Degrees of Disc Degeneration

STUDY DESIGN

Biomechanical cadaveric study.

OBJECTIVE

The aim of this study was to evaluate the effect of degeneration on biomechanical properties of the passive structures of the lumbar spine.

SUMMARY OF BACKGROUND DATA

Although the load apportionment among the passive structures in healthy spines follows well-defined contribution patterns, it remains unknown how this load distribution and sagittal preload changes by degenerative processes of the intervertebral disc (IVD).

METHODS

Fifty lumbar spinal segments were tested in a displacement-controlled stepwise reduction study in flexion, extension, axial rotation, lateral bending, anterior, posterior and lateral shear. The intertransverse ligaments (ITLs), supraspinous and interspinous ligaments (ISL&SSL), facet joint capsules (FJC), facet joints (FJ), ligamentum flavum (LF), posterior longitudinal ligament (PLL), anterior longitudinal ligament (ALL), and spondylophytes were subsequently reduced. The results were set in relation to IVD-degeneration, quantified with Pfirrmann classification.

RESULTS

In flexion, a load redistribution from LF (-28% n.s.) and PLL (-13% n.s.) towards the IVD (+9%, n.s.) is observed comparing grade 2 to 5 IVD degeneration, whereas in all other loading directions, a reduction of IVD-contribution from -12% to -53% is recorded. In axial rotation, anterior and lateral shear, more load is shared by the FJ (+4% n.s., +23% ∗, +13% n.s.). The preload of the ALL, LF, PLL, and IVD is reduced ranging from -0.06 Nm to -0.37 Nm.

CONCLUSION

IVD degeneration is related to notable load-redistributions between the passive spinal structures. With further degeneration, reduced contribution of the LF and PLL and higher loads on the IVD are observed in flexion. In the other tested loading directions, the relative load on the IVD is reduced, whereas higher FJ-exposure in axial rotation, anterior and lateral shear is observed. Furthermore, the preload of the spinal structures is reduced. These observations can further the understanding of the degenerative cascade in the spine.Level of Evidence: N/A.

Intervertebral disc degeneration relates to biomechanical changes of spinal ligaments

BACKGROUND CONTEXT

The ligamentum flavum (LF), the inter- and supraspinous ligament (ISL&SSL) and the intertransverse ligament (ITL) are relevant spinal structures for segmental stability. The biomechanical effect of degeneration and aging on their biomechanical properties remains largely unknown.

PURPOSE

The aim of this study was to assess the material properties of the ITL, ISL&SSL and LF and to correlate parameters of biomechanical function with LF-thickness, intervertebral disc (IVD) degeneration and age.

STUDY DESIGN

Biomechanical cadaveric study.

METHODS

MRI- and CT-scans of 50 human lumbar segments (Th12-L5) were used to assess the ISL (acc. to Keorochana), the grade of IVD degeneration (acc. to Pfirrmann) and to quantify LF-thickness. The ITL, ISL&SSL and LF were resected in the neutral position of the spinal segment with a specifically developed method to conserve initial strain. Ramp to failure testing was performed (0.5 mm/s) to record initial tension, slack length, stiffness and ultimate strength. The relationship between the biomechanical characteristics and age and radiological parameters were analyzed. There are no study-specific conflicts of interest and no external funding was received for this study.

RESULTS

With aging, a significant reduction in initial tension (r=-0.5, p<.01) and ultimate strength (r=-0.41, p<.01) of the LF was observed, while the effect on LF-stiffness and the characteristics of the other ligaments was non-significant. IVD-degeneration was correlated with a significant reduction in stiffness (r=-0.47, p=.001; r=-0.36, p=.01) and ultimate strength (r=-0.3, p=.04; r=-0.36, p=.01) of the LF and ISL&SSL respectively and a significant reduction in initial tension (r=-0.4, p<.01) of the LF. For the ITL, no significant correlation was observed. Comparing Pfirrman 2 to 5, this reduction was 40% to 80% for stiffness 60% to 70% for ultimate strength and 88% for initial tension of the LF. ISL&SSL-stiffness between Kerorochana grade A and D differed significantly (p=.03), while all other comparisons were non-significant (p>.05). LF-thickness did not correlate with the biomechanical properties of the LF (p>.05).

CONCLUSIONS

Aging is primarily related to biomechanical changes to the LF. IVD-degeneration is related to a relevant reduction in stiffness and ultimate strength of the LF and ISL&SSL, with a similar trend for the ITL. The ISL-specific Keorochana grading system provides only minimal biomechanical information and LF-thickness does not provide biomechanical information.

CLINICAL SIGNIFICANCE

Patient age and the degenerative state of the IVD can be used to evaluate the biomechanical characteristics of the dorsal spinal ligaments, which can be helpful in selecting the optimal surgical procedure (e.g. in decompression surgery) for a specific situation.

Is a cross-connector beneficial for single level traditional or cortical bone trajectory pedicle screw instrumentation?

The cortical bone trajectory (CBT) has been introduced with the aim of better screw hold, however, screw-rod constructs with this trajectory might provide less rigidity in lateral bending (LB) and axial rotation (AR) compared to the constructs with the traditional trajectory (TT). Therefore, the addition of a horizontal cross-connector could be beneficial in counteracting this possible inferiority. The aim of this study was to compare the primary rigidity of TT with CBT screw-rod constructs and to quantify the effect of cross-connector-augmentation in both. Spines of four human cadavers (T9 –L5) were cropped into 15 functional spine units (FSU). Eight FSUs were instrumented with TT and seven FSUs with CBT pedicle screws. The segments were tested in six loading directions in three configurations: uninstrumented, instrumented with and without cross-connector. The motion between the cranial and caudal vertebra was recorded. The range of motion (ROM) between the CBT and the TT group did not differ significantly in either configuration. Cross-connector -augmentation did reduce the ROM in AR (16.3%, 0.27°, p = 0.02), LB (2.9%, 0.07°, p = 0.03) and flexion-extension FE (2.3%, 0.04°, p = 0.02) for the TT group and in AR (20.6%, 0.31°, p = 0.01) for the CBT-group. The primary rigidity of TT and CBT single level screw-rod constructs did not show significant difference. The minimal reduction of ROM due to cross-connector-augmentation seems clinically not relevant. Based on the findings of these study there is no increased necessity to use a cross-connector in a CBT-construct.

Tropism of Sub-Axial Cervical Facet Joints Is Not Related to Segmental Movement during Active Movement or Therapist-Perceived Symptomatic Locations

Tropism, or asymmetry, of facet joints in the cervical spine has been found to be related to degenerative changes of the joints and discs. Clinicians often assume that differences in segmental mobility are related to tropism. The aims of this study were to determine the relationship between asymmetry of facet joints in the sub-axial cervical spine and (1) segmental mobility and (2) spinal levels perceived by therapists to have limited mobility. Eighteen participants with idiopathic neck pain had MRIs of their cervical spine in neutral and at the end of active rotation. Angular movement and translational movement of each motion segment was calculated from 3D segmentations of the vertebrae. A plane was fitted to the facet on each side. Tropism was considered to be the difference in the orientation of the facet planes and ranged from 1 to 30° with a median of 7.7°. No relationships were found between the extent of tropism and either segmental movement or locations deemed to be symptomatic. Tropism in the sub-axial cervical spine does not appear to be related to segmental mobility in rotation or to levels deemed to be symptomatic.

A Dynamic Optimization Approach for Solving Spine Kinematics While Calibrating Subject-Specific Mechanical Properties

This study aims to propose a new optimization framework for solving spine kinematics based on skin-mounted markers and estimate subject-specific mechanical properties of the intervertebral joints. The approach enforces dynamic consistency in the entire skeletal system over the entire time-trajectory while personalizing spinal stiffness. 3D reflective markers mounted on ten vertebrae during spine motions were measured in ten healthy volunteers. Biplanar X-rays were taken during neutral stance of the subjects wearing the markers. Calculated spine kinematics were compared to those calculated using inverse kinematics (IK) and IK with imposed generic kinematic constraints. Calculated spine kinematics compared well with standing X-rays, with average root mean square differences of the vertebral body center positions below 10.1 mm and below [Formula: see text] for joint orientation angles. For flexion/extension and lateral bending, the lumbar rotation distribution patterns, as well as the ranges of rotations matched in vivo literature data. The approach outperforms state-of-art IK and IK with constraints methods. Calculated ratios reflect reduced spinal stiffness in low-resistance zone and increased stiffness in high-resistance zone. The patterns of calibrated stiffness were consistent with previously reported experimentally determined patterns. This approach will further our insight into spinal mechanics by increasing the physiological representativeness of spinal motion simulations.

Region- and degeneration dependent stiffness distribution in intervertebral discs derived by shear wave elastography

Information on the local stiffness characteristics of the intervertebral disc (IVD) is crucial for the understanding of its structure-function properties in health and disease and may improve numerical modeling. Previous studies have attempted to map local tissue stiffness by sectioning the disc and performing mechanical testing on these discrete tissue units, which is technically challenging and may bias the results. Shear wave elastography (SWE) represents a nondestructive alternative that can provide spatially continuous elasticity estimates. We investigated the feasibility of SWE for human intervertebral disc elasticity mapping in a laboratory setting. To this end, global spinal segment mechanical behavior was determined in 6 loading directions and served as ground truth data for the validation of the approach. Subsequently, the cranial spinal vertebra was removed and shear wave elastographic scans of the IVD were acquired. SWE-measurements were reconstructed into three-dimensional elastographic maps, discretized into distinct IVD regions and correlated with global segment mechanical parameters. SWE-derived Young's modulus estimates were compared among different regions and as a function of their state of degeneration. We found annulus shear wave speed to be moderately correlated with segment mechanical behavior irrespective of the loading direction whereas shear wave speed in the nucleus pulposus showed a very weak association (mean (SD) absolute Pearson correlation coefficients: 0.51 (0.14) and 0.17 (0.12), respectively). Young's modulus mapping of the intervertebral disc revealed stiffness to be highest in the ventral annulus with a stiffness decrease both circumferentially towards the dorsal aspect as well as towards the center of the disc. SWE hence provides a valid alternative to disc sectioning and piecewise mechanical testing.

Biomechanical contribution of spinal structures to stability of the lumbar spine-novel biomechanical insights

BACKGROUND CONTEXT

The contribution of anatomical structures to the stability of the spine is of great relevance for diagnostic, prognostic and therapeutic evaluation of spinal pathologies. Although a plethora of literature is available, the contribution of anatomical structures is still not well understood.

PURPOSE

We aimed to quantify the biomechanical relevance of each of the passive spinal structure trough deliberate biomechanical test series using a stepwise reduction approach on cadavers.

STUDY DESIGN

Biomechanical cadaveric study.

METHODS

Fifty lumbar spinal segments originating from 22 human lumbar cadavers were biomechanically tested in a displacement-controlled stepwise reduction study: the intertransverse ligaments, the supraspinous and interspinous ligaments, the facet joint capsules (FJC), the facet joints (FJ), the ligamentum flavum (LF), the posterior longitudinal ligament (PLL), and the anterior longitudinal ligament were subsequently reduced. In the intact state and after each transection step, the segments were physiologically loaded in flexion, extension, axial rotation (AR), lateral bending (LB) and with anterior (AS), posterior (PS) and lateral shear (LS). Thirty-two specimens with only minor degeneration, representing a reasonably healthy subpopulation, were selected for the here presented evaluation. Quantitative values for load and spinal level dependent contribution patterns for the anatomical structures were derived.

RESULTS

Small variability between of the contribution patterns are observed. The intervertebral disc (IVD) is exposed to about 67% of the applied load in LB and during shear loading, but less by load in flexion, extension and AR (less than 35%). The FJ&FJC are the main stabilizers in AR with 49%, but provide only 10% of the stability in extension. Beside the IVD, the LF and the PLL contribute mainly in flexion (22% and 16%, respectively), while the ALL plays a major role during extension (40%) and also contributes during LB (15%). The contribution of the intertransverse ligaments and the supraspinous and interspinous ligaments are very small in all loading directions (<2% and <6%, respectively).

CONCLUSION

The IVD takes the main load in LB and absorbs shear loading, while the FJ&FJC stabilize AR. The ALL resists extension while LF and PLL stabilize flexion. With the small variability of contribution patterns, suggesting distinct adaptation of the structures to one another, the biomechanical characteristics of one structure have to be put in context of the whole spinal segment.

CLINICAL SIGNIFICANCE

The novel information on load distribution helps predict the biomechanical consequences of surgical procedures in more detail.

Muscle-driven and torque-driven centrodes during modeled flexion of individual lumbar spines are disparate

Lumbar spine biomechanics during the forward-bending of the upper body (flexion) are well investigated by both in vivo and in vitro experiments. In both cases, the experimentally observed relative motion of vertebral bodies can be used to calculate the instantaneous center of rotation (ICR). The timely evolution of the ICR, the centrode, is widely utilized for validating computer models and is thought to serve as a criterion for distinguishing healthy and degenerative motion patterns. While in vivo motion can be induced by physiological active structures (muscles), in vitro spinal segments have to be driven by external torque-applying equipment such as spine testers. It is implicitly assumed that muscle-driven and torque-driven centrodes are similar. Here, however, we show that centrodes qualitatively depend on the impetus. Distinction is achieved by introducing confidence regions (ellipses) that comprise centrodes of seven individual multi-body simulation models, performing flexion with and without preload. Muscle-driven centrodes were generally directed superior–anterior and tail-shaped, while torque-driven centrodes were located in a comparably narrow region close to the center of mass of the caudal vertebrae. We thus argue that centrodes resulting from different experimental conditions ought to be compared with caution. Finally, the applicability of our method regarding the analysis of clinical syndromes and the assessment of surgical methods is discussed.

Center of rotation locations during lumbar spine movements: a scoping review protocol

OBJECTIVE

The objective of this review is to identify and map current literature describing the center of rotation locations and migration paths during lumbar spine movements.

INTRODUCTION

Altered lumber spine kinematics has been associated with pain and injury. Intervertebral segments' center of rotations, the point around which spinal segments rotate, are important for determining the features of lumbar spine kinematics and the potential for increased injury risk during movements. Although many studies have investigated the center of rotations of humans' lumbar spine, no review has summarized and organized the state of the science related to center of rotation locations and migration paths of the lumbar spine during lumbar spine movements.

INCLUSION CRITERIA

This review will consider studies that include human lumbar spines of any age and condition (e.g. heathy, pathological) during lumbar spine movements. Quantitative study designs, including clinical, observational, laboratory biomechanical experimental studies, mathematical and computer modeling studies will be considered. Only studies published in English will be included, and there will be no limit on dates of publication.

METHODS

PubMed, MEDLINE, Embase, the Cochrane Library Controlled Register of Trials, CINAHL, ACM Digital Library, Compendex, Inspec, Web of Science, Scopus, Google Scholar, and dissertation and theses repositories will be searched. After title and abstract screening of identified references, two independent reviewers will screen the full-text of identified studies and extract data. Data will be summarized and categorized, and a comprehensive narrative summary will be presented with the respective results.

The effect of head and gaze orientation on spine kinematics during forward flexion

Compound, or awkward, spine postures have been suggested as a biomechanical risk factor for low back injury. This experiment investigates the influence of head (i.e. head-on-torso) and gaze (i.e. eye-in-head) orientation on three-dimensional (3D) neck and spine range of motion (ROM) during forward flexion movements. To emulate previous experimental protocols and replicate real-world scenarios, a sample of ten young, healthy males (mean ± standard deviation: age: 20.8 ± 1.03 years, height: 180.2 ± 7.36 cm, and mass: 81.9 ± 6.47 kg) completed forward flexion movements with a constrained and unconstrained pelvis, respectively. Surface kinematics were gathered from the head and spine (C₇-S₁). Movements were completed under a baseline condition as well as upward, downward, leftward, and rightward head and gaze orientations. For each condition, mean neck angle and inter-segmental spine (C₇T₁ through L₅S₁) ROM were evaluated. The results demonstrate that directed head and gaze orientations can influence the ROM of specific spine regions during a forward flexion task. With leftward and rightward directed head and gaze orientations, the neck became increasingly twisted and superior thoracic segments (i.e. C₇T₁-T₂T₃) were significantly more twisted during the leftward head orientation condition than the baseline condition. With upward and downward directed head and gaze orientations, a similar effect was observed for neck and superior thoracic (i.e. C₇T₁-T₄T₅) flexion-extension. Interestingly, it was also demonstrated that changes in upward/downward head orientation can also change flexion-extension kinematics of the thoracolumbar region as well (i.e. T₇T₈-L₁L₂), suggesting that head postures requiring neck extension may also promote extension throughout these spine regions. These findings provide evidence for a functional link between changes in neck flexion-extension posture and flexion-extension movement of the thoracolumbar region of the spine.