Altered helical axis patterns of the lumbar spine indicate increased instability with disc degeneration

Increased spinal kinematic variability in people with chronic low back pain revealed by alterations in helical axis parameters

Changes in spine kinematics are common in people with chronic low back pain (CLBP) and this includes changes in trunk angular displacement and angular velocity. The helical axis (HA) of motion is an approach to investigate three-dimensional variability of joint kinematics. In this study we investigate whether the variability of trunk movement is modified in people with CLBP by measuring the dispersion of HA parameters during repeated trunk movements performed at different movement speed. Nineteen people with CLBP and twenty asymptomatic volunteers performed repetitive continuous trunk movements at three different speeds. Two parameters of the HA were extracted to characterise movement variability at the lumbo-sacral and thoraco-lumbar regions: mean angle (MA) and mean distance (MD). Two-Way mixed ANOVA showed significantly higher MA and MD (p<0.001) especially at the thoraco-lumbar region for those with CLBP. Interestingly, this was not consistent across all directions or speed of movement; higher HA parameters for those with CLBP at the lumbo-sacral region was mainly observed during flexion/extension cycles. In addition, there was a speed and group interaction effect during rotational cycles (p=0.010) which showed higher thoraco-lumbar MA values for those with CLBP during the faster speed (p=0.029, mean dif.(95 % CI) = 2.28, (0.247;4.328)) and slower speed condition (p=0.003, mean dif.(95 % CI) = 2.78, (1.009;4.565)). This study shows that people with CLBP move their spine in a more variable way, a characteristic that could be influenced by speed and direction of trunk movement. This could reflect an adaptive behaviour to long-lasting pain.

In vivo kinematic study of lumbar center of rotation under different loads

BACKGROUND

Dual fluoroscopic imaging system (DFIS) was employed to identify the Center of Rotation(COR) in the lower lumbar spine and determine its relationship with weight bearing.

METHODS

In this study, twenty participants were recruited. A 3D model of each participant's lumbar spine was created using CT images, and their relative positions were determined through DFIS. By integrating CT imaging with DFIS, the kinematic data of the participants' spines during movement were captured. The lower lumbar spine's COR was calculated using the method of perpendicular bisectors.

RESULTS

While flexing and extending, the Center of Rotation (COR) initially moved downward with increasing load, followed by upward movement as the load further increased. During flexion and extension, the COR coordinates of L3-4 at 0 kg, 5 kg and 10 kg are(0.3549 ± 0.2176,0.0177 ± 0.1317),(0.0598 ± 0.2095,-0.1806 ± 0.1719),(0.1427 ± 0.1440,-0.0911 ± 0.2722); The center of rotation coordinates of L4-5 at 0 kg, 5 kg and 10 kg are(0.0566 ± 0.2693,-0.0727 ± 0.2132),(0.0964 ± 0.2671,-0.2037 ± 0.2299),(0.1648 ± 0.1520,-0.0049 ± 0.1641). The anterior-posterior position of the COR shifted posteriorly with increasing weight-bearing. During lateral bending, the center of rotation coordinates of L3-4 at 0 kg, 5 kg and 10 kg are(0.0745 ± 0.1229,0.0966 ± 0.3403) (-0.0438 ± 0.1281,0.1161 ± 0.1584), (-0.0464 ± 0.1517,0.1320 ± 0.2730); The center of rotation coordinates of L4-5 at 0 kg, 5 kg and 10 kg are(-0.0314 ± 0.1411,-0.0355 ± 0.2088), (-0.0764 ± 0.3135,0.0105 ± 0.3230),(-0.0376 ± 0.1701,0.0285 ± 0.2395). Throughout the lateral bending exercises, the upper and lower COR positions increased as the load increased, while the left and right COR positions remained unaffected by the load increment. The COR height differed between flexion and lateral bending. We observed variations in the COR position of the lumbar spine during lateral bending and flexion-extension movements. This enhanced our comprehension of coupled motion patterns within the lumbar spine.

CONCLUSIONS

Position of the lumbar spine COR changes with variations in the load. During different movements, the COR location of the lower lumbar spine varied. This finding suggests the presence of distinct motion patterns in the lower lumbar spine. As the load increases, the lumbar COR position changes significantly. Abnormal movement patterns of the lower lumbar spine under different loads may be one of the factors that accelerate lumbar disc degeneration.

Automated magnetic resonance imaging-based grading of the lumbar intervertebral disc and facet joints

Background

Degeneration of both intervertebral discs (IVDs) and facet joints in the lumbar spine has been associated with low back pain, but whether and how IVD/joint degeneration contributes to pain remains an open question. Joint degeneration can be identified by pairing T1 and T2 magnetic resonance imaging (MRI) with analysis techniques such as Pfirrmann grades (IVD degeneration) and Fujiwara scores (facet degeneration). However, these grades are subjective, prompting the need to develop an automated technique to enhance inter-rater reliability. This study introduces an automated convolutional neural network (CNN) technique trained on clinical MRI images of IVD and facet joints obtained from public-access Lumbar Spine MRI Dataset. The primary goal of the automated system is to classify health of lumbar discs and facet joints according to Pfirrmann and Fujiwara grading systems and to enhance inter-rater reliability associated with these grading systems.

Methods

Performance of the CNN on both the Pfirrmann and Fujiwara scales was measured by comparing the percent agreement, Pearson's correlation and Fleiss kappa value for results from the classifier to the grades assigned by an expert grader.

Results

The CNN demonstrates comparable performance to human graders for both Pfirrmann and Fujiwara grading systems, but with larger errors in Fujiwara grading. The CNN improves the reliability of the Pfirrmann system, aligning with previous findings for IVD assessment.

Conclusion

The study highlights the potential of using deep learning in classifying the IVD and facet joint health, and due to the high variability in the Fujiwara scoring system, highlights the need for improved imaging and scoring techniques to evaluate facet joint health. All codes required to use the automatic grading routines described herein are available in the Data Repository for University of Minnesota (DRUM).

Dynamic in vivo 3D atlantooccipital kinematics during multiplanar physiologic motions

Normal biomechanics of the upper cervical spine, particularly at the atlantooccipital joint, remain poorly characterized. The purpose of this study was to determine the intervertebral kinematics of the atlantooccipital joint (occiput-C1) during three-dimensional in vivo physiologic movements. Twenty healthy young adults performed dynamic flexion/extension, axial rotation, and lateral bending while biplane radiographs were collected at 30 images per second. Motion at occiput-C1 was tracked using a validated volumetric model-based tracking process that matched subject-specific CT-based bone models to the radiographs. The occiput-C1 total range of motion (ROM) and helical axis of motion (HAM) was calculated for each movement. During flexion/extension, the occiput-C1 moved almost exclusively in-plane (ROM: 17.9 ± 6.9°) with high variability in kinematic waveforms (6.3°) compared to the in-plane variability during axial rotation (1.4°) and lateral bending (0.9°) movements. During axial rotation, there was small in-plane motion (ROM: 4.2 ± 2.5°) compared to out-of-plane flexion/extension (ROM: 12.7 ± 5.4°). During lateral bending, motion occurred in-plane (ROM: 9.0 ± 3.1°) and in the plane of flexion/extension (ROM: 7.3 ± 2.7°). The average occiput-C1 axis of rotation intersected the sagittal and coronal planes 7 mm to 18 mm superior to the occipital condyles. The occiput-C1 axis of rotation pointed 60° from the sagittal plane during axial rotation but only 10° from the sagittal plane during head lateral bending. These novel results are foundational for future work on upper cervical spine kinematics.

Transforaminal endoscopic lumbar discectomy with two-segment foraminoplasty for the treatment of very highly migrated lumbar disc herniation: a retrospective analysis

BACKGROUND

The surgical resection of very highly migrated lumbar disc herniation (VHM-LDH) is technically challenging owing to the absence of technical guidelines. Hence, in the present study, we introduced the transforaminal endoscopic lumbar discectomy (TELD) with two-segment foraminoplasty to manage VHM-LDH and evaluated its radiographic and midterm clinical outcomes.

MATERIALS AND METHODS

The present study is a retrospective analysis of 33 consecutive patients with VHM-LDH who underwent TELD with two-segment foraminoplasty. The foraminoplasty was performed on two adjacent vertebrae on the basis of the migration direction of disc fragments to fully expose the disc fragments and completely decompress the impinged nerve root. The operation duration, blood loss, intra- and postoperative complications, and recurrences were recorded. Additionally, imageological observations were evaluated immediately after the procedure via magnetic resonance image and computerized tomography. Clinical outcomes were evaluated by calculating the visual analog scale (VAS) score and Oswestry Disability Index (ODI). The MacNab criterion was reviewed to assess the patients' opinions on treatment satisfaction. The resection rate of bony structures were quantitatively evaluated on postoperative image. The segmental stability was radiologically evaluated at least a year after the surgery. Additionally, surgery-related and postoperative complications were evaluated.

RESULTS

The average age of the patients was 56.87 ± 7.77 years, with a mean follow-up of 20.95 ± 2.09 months. The pain was relieved in all patients immediately after the surgery. The VAS score and ODI decreased significantly at each postoperative follow-up compared with those observed before the surgery (P < 0.05). The mean operation duration, blood loss, and hospital stay were 56.17 ± 16.21 min, 10.57 ± 6.92 mL, and 3.12 ± 1.23 days, respectively. No residual disc fragments, iatrogenic pedicle fractures, and segmental instability were observed in the postoperative images. For both up- and down- migrated herniation in the upper lumbar region, the upper limit value of resection percentage for the cranial SAP, caudal SAP, and pedicle was 33%, 30%, and 34%, respectively; while those in the lower lumbar region was 42%, 36%, and 46%, respectively. At the last follow-up, the satisfaction rate of the patients regarding the surgery was 97%. Surgery-related complications including dural tear, nerve root injury, epidural hematoma, iatrogenic pedicle fractures, and segmental instability were not observed. One patient (3%) suffered from the recurrence of LDH 10 months after the initial surgery and underwent revision surgery.

CONCLUSIONS

The TELD with two-segment foraminoplasty is safe and effective for VHM-LDH management. Proper patient selection and efficient endoscopic skills are required for applying this technique to obtain satisfactory outcomes.

Cervical facet capsular ligament mechanics: Estimations based on subject-specific anatomy and kinematics

Background

To understand the facet capsular ligament's (FCL) role in cervical spine mechanics, the interactions between the FCL and other spinal components must be examined. One approach is to develop a subject-specific finite element (FE) model of the lower cervical spine, simulating the motion segments and their components' behaviors under physiological loading conditions. This approach can be particularly attractive when a patient's anatomical and kinematic data are available.

Methods

We developed and demonstrated methodology to create 3D subject-specific models of the lower cervical spine, with a focus on facet capsular ligament biomechanics. Displacement-controlled boundary conditions were applied to the vertebrae using kinematics extracted from biplane videoradiography during planar head motions, including axial rotation, lateral bending, and flexion-extension. The FCL geometries were generated by fitting a surface over the estimated ligament-bone attachment regions. The fiber structure and material characteristics of the ligament tissue were extracted from available human cervical FCL data. The method was demonstrated by application to the cervical geometry and kinematics of a healthy 23-year-old female subject.

Results

FCL strain within the resulting subject-specific model were subsequently compared to models with generic: (1) geometry, (2) kinematics, and (3) material properties to assess the effect of model specificity. Asymmetry in both the kinematics and the anatomy led to asymmetry in strain fields, highlighting the importance of patient-specific models. We also found that the calculated strain field was largely independent of constitutive model and driven by vertebrae morphology and motion, but the stress field showed more constitutive-equation-dependence, as would be expected given the highly constrained motion of cervical FCLs.

Conclusions

The current study provides a methodology to create a subject-specific model of the cervical spine that can be used to investigate various clinical questions by coupling experimental kinematics with multiscale computational models.

The Lumbar Facet Capsular Ligament Becomes More Anisotropic and the Fibers Become Stiffer With Intervertebral Disc and Facet Joint Degeneration

Degeneration of the lumbar spine, and especially how that degeneration may lead to pain, remains poorly understood. In particular, the mechanics of the facet capsular ligament may contribute to low back pain, but the mechanical changes that occur in this ligament with spinal degeneration are unknown. Additionally, the highly nonlinear, heterogeneous, and anisotropic nature of the facet capsular ligament makes understanding mechanical changes more difficult. Clinically, magnetic resonance imaging (MRI)-based signs of degeneration in the facet joint and the intervertebral disc (IVD) correlate. Therefore, this study examined how the nonlinear, heterogeneous mechanics of the facet capsular ligament change with degeneration of the lumbar spine as characterized using MRI. Cadaveric human spines were imaged via MRI, and the L2-L5 facet joints and IVDs were scored using the Fujiwara and Pfirrmann grading systems. Then, the facet capsular ligament was isolated and biaxially loaded. The nonlinear mechanical properties of the ligament were obtained using a nonlinear generalized anisotropic inverse mechanics analysis (nGAIM). Then a Holzapfel-Gasser-Ogden (HGO) model was fit to the stress-strain data obtained from nGAIM. The facet capsular ligament is stiffer and more anisotropic at larger Pfirrmann grades and higher Fujiwara scores than at lower grades and scores. Analysis of ligament heterogeneity showed all tissues are highly heterogeneous, but no distinct spatial patterns of heterogeneity were found. These results show that degeneration of the lumbar spine including the facet capsular ligament appears to be occurring as a whole joint phenomenon and advance our understanding of lumbar spine degeneration.

Calibration and validation of a novel hybrid model of the lumbosacral spine in ArtiSynth-The passive structures

In computational biomechanics, two separate types of models have been used predominantly to enhance the understanding of the mechanisms of action of the lumbosacral spine (LSS): Finite element (FE) and musculoskeletal multibody (MB) models. To combine advantages of both models, hybrid FE-MB models are an increasingly used alternative. The aim of this paper is to develop, calibrate, and validate a novel passive hybrid FE-MB open-access simulation model of a ligamentous LSS using ArtiSynth. Based on anatomical data from the Male Visible Human Project, the LSS model is constructed from the L1-S1 rigid vertebrae interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, and facet joints. A mesh convergence study, sensitivity analyses, and systematic calibration were conducted with the hybrid functional spinal unit (FSU) L4/5. The predicted mechanical responses of the FSU L4/5, the lumbar spine (L1-L5), and the LSS were validated against literature data from in vivo and in vitro measurements and in silico models. Spinal mechanical responses considered when loaded with pure moments and combined loading modes were total and intervertebral range of motions, instantaneous axes and centers of rotation, facet joint contact forces, intradiscal pressures, disc bulges, and stiffnesses. Undesirable correlations with the FE mesh were minimized, the number of crisscrossed collagen fiber rings was reduced to five, and the individual influences of specific anatomical structures were adjusted to in vitro range of motions. Including intervertebral motion couplings for axial rotation and nonlinear stiffening under increasing axial compression, the predicted kinematic and structural mechanics responses were consistent with the comparative data. The results demonstrate that the hybrid simulation model is robust and efficient in reproducing valid mechanical responses to provide a starting point for upcoming optimizations and extensions, such as with active skeletal muscles.

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.

Three dimensional finite element analysis used to study the influence of the stress and strain of the operative and adjacent segments through different foraminnoplasty technique in the PELD: Study protocol clinical trial (SPIRIT Compliant)

INTRODUCTION

Percutaneous endoscopic lumbar disectomy (PELD) is one of the most popular minimally invasive techniques of spinal surgery in recent years. At present, there are 2 main surgical approaches in PELD: foraminal approach and interlaminar approach. What's more, foraminoplasty is a necessary step for both approaches. However, there are few biomechanical studies on the formation of different parts of the intervertebral foramen. The aim of this study is to explore the effects of different foraminoplasty methods on the biomechanics of the corresponding and adjacent segments of the lumbar through a 3-dimensional finite element model analysis.

METHODS

We established a normal 3-dimensional finite element mode of L3 to L5, simulated lumbar percutaneous endoscopy by doing cylindrical excision of bone whose diameter was 7.5 mm on the L5 superior articular process and the L4 inferior articular process, respectively, so that we obtained 3 models: the first one was normal lumbar model, the second one was the L4 inferior articular process shaped model, and the third one was the L5 superior articular process shaped model. We compared the biomechanics of the intervertebral disc of L3/4 and L4/5 when they were in the states of forward flexion, backward extension, left and right flexion, and left and right rotation on specific loading condition.

DISCUSSION

If the outcomes indicate the trial is feasible and there is evidence that one of the foraminoplasty technique may make few differences in biomechanics of corresponding lumbar intervertebral disc, we will proceed to a definitive trial to test the best way to foraminplasty, which could make biomechanical influence as little as possible.

TRIAL REGISTRATION

Chinese Clinical Trial Registry, ChiCTR1900026973. Registered on September 27, 2019.

Analysis of clinical effect and radiographic outcomes of Isobar TTL system for two-segment lumbar degenerative disease: a retrospective study

Background

Nonfusion fixation is an effective way to treat lumbar degeneration. In the present study, we analyzed the clinical effects and radiographic outcomes of the Isobar TTL system used to treat two-segment lumbar degenerative disease.

Method

Forty-one patients diagnosed with two-segment lumbar degenerative disease underwent surgical implantation of the Isobar TTL dynamic stabilization system (n = 20) or a rigid system (n = 21) from January 2013 to June 2017. The mean follow-up time was 23.6 (range 15–37) months. Clinical results were evaluated with the Oswestry Disability Index (ODI), modified Macnab criteria, and the visual analog score (VAS). Radiographic evaluations included the height of the intervertebral space and the range of motion (ROM) of the treated and adjacent segments. The intervertebral disc signal was classified using the modified Pfirrmann grading system and the University of California at Los Angeles (UCLA) system.

Results

The clinical results, including the ODI and VAS, showed that there was significant improvement in the two groups after implantation and that the difference between the two groups was not significant. In addition, the clinical efficacy indicated by the modified Macnab criteria for the two groups was similar. Radiological outcomes included the height of the intervertebral space, lumbar mobility, and intervertebral disc signal. The height of the intervertebral space of the upper adjacent segment L2/3 in the rigid group was significantly lower than that in the Isobar TTL group at the last follow-up. Furthermore, the number of ROMs of the fixed-segment L3/4 in the Isobar TTL group was significantly less than that before implantation, suggesting that the fixed-segment ROMs in the Isobar TTL group were limited. In addition, the ROM of the upper adjacent segment L2/3 in the last follow-up of the rigid group increased significantly, while that of the Isobar TTL group did not change after implantation. Finally, the incidence of adjacent-segment degeneration (ASD) was significantly greater in the rigid group than in the Isobar TTL group according to the UCLA system.

Conclusion

The Isobar TTL system can be clinically effective for treating two-segment lumbar degenerative disease. Compared with rigid fixation, the Isobar TTL system yielded better radiographic outcomes and maintained the mobility of the treated segments with less impact on the proximal adjacent segment.

Influence of Deviated Centers of Rotation on Kinematics and Kinetics of a Lumbar Functional Spinal Unit: An In Vitro Study

Background

Center of rotation (COR) has been used for assessing spinal motion quality. However, the biomechanical influence of COR deviation towards different directions during flexion-extension (FE) remains largely unknown. This study aimed to investigate the alteration in the range of motion (ROM), compressive force, shear force, and neutral zone size (NZ) in a lumbar functional spinal unit (FSU), caused by the deviated COR in different directions during FE.

Material/Methods

Twelve human cadaveric lumbar FSUs (6 for L2–L3, 6 for L4–L5) were tested in a 6-degree-of-freedom servo-hydraulic load frame. These FSUs were firstly applied a 7.5 Nm pure moment to perform FE to obtain their natural COR during FE. Subsequently, they were subjected to FE around 9 established deviated CORs with 6 Nm cyclical loading.

Results

It was found that the ROM and NZ increased significantly when the COR moved from the superior plane to the inferior plane for the L2–L3 unit and when the COR located in the superior plane compared with the inferior plane for the L4–L5 unit. The compressive forces for both FSUs demonstrated significant changes caused by COR shift in the same horizontal plane, while the shear forces demonstrated significant changes caused by COR shift in the same vertical plane.

Conclusions

The ROM, NZ, and shear force of FSU are sensitive to the vertical COR shift, while the compressive force of FSU is highly sensitive to the horizontal COR shift. Additionally, the kinematics and kinetics of the L2–L3 unit are more sensitive to COR location than those of the L4–L5 unit.

The role of the facet capsular ligament in providing spinal stability

Low back pain (LBP) is the most common type of pain in America, and spinal instability is a primary cause. The facet capsular ligament (FCL) encloses the articulating joints of the spine and is of particular interest due to its high innervation - as instability ensues, high stretch values likely are a cause of this pain. Therefore, this work investigated the FCL's role in providing stability to the lumbar spine. A previously validated finite element model of the L4-L5 spinal motion segment was used to simulate pure moment bending in multiple planes. FCL failure was simulated and the following outcome measures were calculated: helical axes of motion, range of motion (ROM), bending stiffness, facet joint space, and FCL stretch. ROM increased, bending stiffness decreased, and altered helical axis patterns were observed with the removal of the FCL. Additionally, a large increase in FCL stretch was measured with diminished FCL mechanical competency, providing support that the FCL plays an important role in spinal stability.

Relationships between lumbar inter-vertebral motion and lordosis in healthy adult males: a cross sectional cohort study

Background

Intervertebral motion impairment is widely thought to be related to chronic back disability, however, the movements of inter-vertebral pairs are not independent of each other and motion may also be related to morphology. Furthermore, maximum intervertebral range of motion (IV-RoMmax) is difficult to measure accurately in living subjects. The purpose of this study was to explore possible relationships between (IV-RoMmax) and lordosis, initial attainment rate and IV-RoMmax at other levels during weight-bearing flexion using quantitative fluoroscopy (QF).

Methods

Continuous QF motion sequences were recorded during controlled active sagittal flexion of 60° in 18 males (mean age 27.6 SD 4.4) with no history of low back pain in the previous year. IV-RoMmax, lordotic angle, and initial attainment rate at all inter-vertebral levels from L2-S1 were extracted. Relationships between IV-RoMmax and the other variables were explored using correlation coefficients, and simple linear regression was used to determine the effects of any significant relationships. Within and between observer repeatability of IV-RoMmax and initial attainment rate measurements were assessed in a sub-set of ten participants, using the intra-class correlation coefficient (ICC) and standard error of measurement (SEM).

Results

QF measurements were highly repeatable, the lowest ICC for IV-RoMmax, being 0.94 (0.80–0.99) and highest SEM (0.76°). For initial attainment rate the lowest ICC was 0.84 (0.49–0.96) and the highest SEM (0.036). The results also demonstrated significant positive and negative correlations between IV-RoMmax and IV-RoMmax at other lumbar levels (r = −0.64–0.65), lordosis (r = −0.52–0.54), and initial attainment rate (r = −0.64–0.73). Simple linear regression analysis of all significant relationships showed that these predict between 28 and 42 % of the variance in IV-RoMmax.

Conclusions

This study found weak to moderate effects of individual kinematic variables and lumbar lordosis on IV-RoMmax at other intervertebral levels. These effects, when combined, may be important when such levels are being considered by healthcare professionals as potential sources of pain generation. Multivariate investigations in larger samples are warranted.

Electronic supplementary material

The online version of this article (doi:10.1186/s12891-016-0975-1) contains supplementary material, which is available to authorized users.

Comparative role of disc degeneration and ligament failure on functional mechanics of the lumbar spine

Understanding spinal kinematics is essential for distinguishing between pathological conditions of spine disorders, which ultimately lead to low back pain. It is of high importance to understand how changes in mechanical properties affect the response of the lumbar spine, specifically in an effort to differentiate those associated with disc degeneration from ligamentous changes, allowing for more precise treatment strategies. To do this, the goals of this study were twofold: (1) develop and validate a finite element (FE) model of the lumbar spine and (2) systematically alter the properties of the intervertebral disc and ligaments to define respective roles in functional mechanics. A three-dimensional non-linear FE model of the lumbar spine (L3-sacrum) was developed and validated for pure moment bending. Disc degeneration and sequential ligament failure were modelled. Intersegmental range of motion (ROM) and bending stiffness were measured. The prediction of the FE model to moment loading in all three planes of bending showed very good agreement, where global and intersegmental ROM and bending stiffness of the model fell within one standard deviation of the in vitro results. Degeneration decreased ROM for all directions. Stiffness increased for all directions except axial rotation, where it initially increased then decreased for moderate and severe degeneration, respectively. Incremental ligament failure produced increased ROM and decreased stiffness. This effect was much more pronounced for all directions except lateral bending, which is minimally impacted by ligaments. These results indicate that lateral bending may be more apt to detect the subtle changes associated with degeneration, without being masked by associated changes of surrounding stabilizing structures.

Use of the alpha shape to quantify finite helical axis dispersion during simulated spine movements

In biomechanical studies examining joint kinematics the most common measurement is range of motion (ROM), yet other techniques, such as the finite helical axis (FHA), may elucidate the changes in the 3D motion pathology more effectively. One of the deficiencies with the FHA technique is in quantifying the axes generated throughout a motion sequence. This study attempted to solve this issue via a computational geometric technique known as the alpha shape, which bounds a set of point data within a closed boundary similar to a convex hull. The purpose of this study was to use the alpha shape as an additional tool to visualize and quantify FHA dispersion between intact and injured cadaveric spine movements and compare these changes to the gold-standard ROM measurements. Flexion-extension, axial rotation, and lateral bending were simulated with five C5-C6 motion segments using a spinal loading simulator and Optotrak motion tracking system. Specimens were first tested intact followed by a simulated injury model. ROM and the FHAs were calculated post-hoc, with alpha shapes and convex hulls generated from the anatomic planar intercept points of the FHAs. While both ROM and the boundary shape areas increased with injury (p<0.05), no consistent geometric trends in the alpha shape growth were identified. The alpha shape area was sensitive to the alpha value chosen and values examined below 2.5 created more than one closed boundary. Ultimately, the alpha shape presents as a useful technique to quantify sequences of joint kinematics described by scatter plots such as FHA intercept data.