A Method to Track 3D Knee Kinematics by Multi-Channel 3D-Tracked A-Mode Ultrasound.

This paper introduces a method for measuring 3D tibiofemoral kinematics using a multi-channel A-mode ultrasound system under dynamic conditions. The proposed system consists of a multi-channel A-mode ultrasound system integrated with a conventional motion capture system (i.e., optical tracking system). This approach allows for the non-invasive and non-radiative quantification of the tibiofemoral joint's six degrees of freedom (DOF). We demonstrated the feasibility and accuracy of this method in the cadaveric experiment. The knee joint's motions were mimicked by manually manipulating the leg through multiple motion cycles from flexion to extension. To measure it, six custom ultrasound holders, equipped with a total of 30 A-mode ultrasound transducers and 18 optical markers, were mounted on various anatomical regions of the lower extremity of the specimen. During experiments, 3D-tracked intra-cortical bone pins were inserted into the femur and tibia to measure the ground truth of tibiofemoral kinematics. The results were compared with the tibiofemoral kinematics derived from the proposed ultrasound system. The results showed an average rotational error of 1.51 ± 1.13° and a translational error of 3.14 ± 1.72 mm for the ultrasound-derived kinematics, compared to the ground truth. In conclusion, this multi-channel A-mode ultrasound system demonstrated a great potential of effectively measuring tibiofemoral kinematics during dynamic motions. Its improved accuracy, nature of non-invasiveness, and lack of radiation exposure make this method a promising alternative to incorporate into gait analysis and prosthetic kinematic measurements later.

Sensitivity of muscle and intervertebral disc force computations to variations in muscle attachment sites.

The current paper aims at assessing the sensitivity of muscle and intervertebral disc force computations against potential errors in modeling muscle attachment sites. We perturbed each attachment location in a complete and coherent musculoskeletal model of the human spine and quantified the changes in muscle and disc forces during standing upright, flexion, lateral bending, and axial rotation of the trunk. Although the majority of the muscles caused minor changes (less than 5%) in the disc forces, certain muscle groups, for example, quadratus lumborum, altered the shear and compressive forces as high as 353% and 17%, respectively. Furthermore, percent changes were higher in the shear forces than in the compressive forces. Our analyses identified certain muscles in the rib cage (intercostales interni and intercostales externi) and lumbar spine (quadratus lumborum and longissimus thoracis) as being more influential for computing muscle and disc forces. Furthermore, the disc forces at the L4/L5 joint were the most sensitive against muscle attachment sites, followed by T6/T7 and T12/L1 joints. Presented findings suggest that modeling muscle attachment sites based on solely anatomical illustrations might lead to erroneous evaluation of internal forces and promote using anatomical datasets where these locations were accurately measured. When developing a personalized model of the spine, certain care should also be paid especially for the muscles indicated in this work.

Twente Spine Model: A thorough investigation of the spinal loads in a complete and coherent musculoskeletal model of the human spine.

Although in vivospinal loads have been previously measured, existing data are limited to certain lumbar and thoracic levels. A detailed investigation of spinal loads would assist with injury prevention and implant design but is unavailable. In this study, we developed a complete and coherent musculoskeletal model of the entire human spine and studied the intervertebral disc compression forces for physiological movements on three anatomical planes. This model incorporates the individual vertebrae at the cervical, thoracic, and lumbar regions, a flexible ribcage, and complete muscle anatomy. Intradiscal pressures were estimated from predicted compressive forces, and these were generally in close agreement with previously measured data. We found that compressive forces at the trunk discs increased during trunk lateral bending and axial rotation of the trunk. During flexion, compressive forces increased in the thoracolumbar and lumbar regions and slightly decreased at the middle thoracic discs. In extension, the forces generally decreased at the thoracolumbar and lumbar discs whereas they slightly increased at the upper and middle thoracic discs. Furthermore, similar to a previous biomechanical model of the cervical spine, our model predicted increased compression forces in neck flexion, lateral bending, and axial rotation, and decreased forces in neck extension.

A Novel Ultrasound-Based Lower Extremity Motion Tracking System.

Tracking joint motion of the lower extremity is important for human motion analysis. In this study, we present a novel ultrasound-based motion tracking system for measuring three-dimensional (3D) position and orientation of the femur and tibia in 3D space and quantifying tibiofemoral kinematics under dynamic conditions. As ultrasound is capable of detecting underlying bone surface noninvasively through multiple layers of soft tissues, an integration of multiple A-mode ultrasound transducers with a conventional motion tracking system provides a new approach to track the motion of bone segments during dynamic conditions. To demonstrate the technical and clinical feasibilities of this concept, an in vivo experiment was conducted. For this purpose the kinematics of healthy individuals were determined in treadmill walking conditions and stair descending tasks. The results clearly demonstrated the potential of tracking skeletal motion of the lower extremity and measuring six-degrees-of-freedom (6-DOF) tibiofemoral kinematics and related kinematic alterations caused by a variety of gait parameters. It was concluded that this prototyping system has great potential to measure human kinematics in an ambulant, non-radiative, and noninvasive manner.

Feasibility of A-mode ultrasound based intraoperative registration in computer-aided orthopedic surgery: A simulation and experimental study.

PURPOSE: A fast and accurate intraoperative registration method is important for Computer-Aided Orthopedic Surgery (CAOS). A-mode ultrasound (US) is able to acquire bone surface data in a non-invasive manner. To utilize A-mode US in CAOS, a suitable registration algorithm is necessary with a small number of registration points and the presence of measurement errors. Therefore, we investigated the effects of (1) the number of registration points and (2) the Ultrasound Point Localization Error (UPLE) on the overall registration accuracy. METHODS: We proposed a new registration method (ICP-PS), including the Iterative Closest Points (ICP) algorithm and a Perturbation Search algorithm. This method enables to avoid getting stuck in the local minimum of ICP iterations and to find the adjacent global minimum. This registration method was subsequently validated in a numerical simulation and a cadaveric experiment using a 3D-tracked A-mode US system. RESULTS: The results showed that ICP-PS outperformed the standard ICP algorithm. The registration accuracy improved with the addition of ultrasound registration points. In the numerical simulation, for 25 sample points with zero UPLE, the averaged registration error of ICP-PS reached 0.25 mm, while 1.71 mm for ICP, decreasing by 85.38%. In the cadaver experiment, using 25 registration points, ICP-PS achieved an RMSE of 2.81 mm relative to 5.84 mm for the ICP, decreasing by 51.88%. CONCLUSIONS: The simulation approach provided a well-defined framework for estimating the necessary number of ultrasound registration points and acceptable level of UPLE for a given required level of accuracy for intraoperative registration in CAOS. ICP-PS method is suitable for A-mode US based intraoperative registration. This study would facilitate the application of A-mode US probe in registering the point cloud to a known shape model, which also has the potential for accurately estimating bone position and orientation for skeletal motion tracking and surgical navigation.

Measuring relative positions and orientations of the tibia with respect to the femur using one-channel 3D-tracked A-mode ultrasound tracking system: A cadaveric study.

The purpose of this study is to investigate the technical feasibility of measuring relative positions and orientations of the tibia with respect to the femur in an in-vitro experiment by using a 3D-tracked A-mode ultrasound system and to determine its accuracy of angular and translational measurements. As A-mode ultrasound is capable of detecting bone surface through soft tissue in a non-invasive manner, the combination of a single A-mode ultrasound transducer with an optical motion tracking system provides the possibility for digitizing the 3D locations of bony points at different anatomical regions on the thigh and the shank. After measuring bony points over a large area of both the femur and tibia, the bone models of the femur and tibia that were segmented from CT or MRI images were registered to the corresponding bony points. Then the relative position of the tibia with respect to the femur could be obtained and the angular and translational components could also be measured. A cadaveric experiment was conducted to assess its accuracy compared to the reference measurement obtained by optical markers fixed to intra-cortical bone pins placed in the femur and tibia. The results showed that the ultrasound system could achieve 0.49 ±â€¯0.83°, 0.85 ±â€¯1.86° and 1.85 ±â€¯2.78° (mean ±â€¯standard deviation) errors for Flexion-Extension, Adduction-Abduction and External-Internal rotations, respectively, and -2.22 ±â€¯3.62 mm, -2.80 ±â€¯2.35 mm and -1.44 ±â€¯2.90 mm errors for Anterior-Posterior, Proximal-Distal and Lateral-Medial translations, respectively. It was concluded that this technique is feasible and facilitates the integration of arrays of A-mode ultrasound transducers with an optical motion tracking system for non-invasive dynamic tibiofemoral kinematics measurement.

In situ comparison of A-mode ultrasound tracking system and skin-mounted markers for measuring kinematics of the lower extremity.

Skin-mounted marker based motion capture systems are widely used in measuring the movement of human joints. Kinematic measurements associated with skin-mounted markers are subject to soft tissue artifacts (STA), since the markers follow skin movement, thus generating errors when used to represent motions of underlying bone segments. We present a novel ultrasound tracking system that is capable of directly measuring tibial and femoral bone surfaces during dynamic motions, and subsequently measuring six-degree-of-freedom (6-DOF) tibiofemoral kinematics. The aim of this study is to quantitatively compare the accuracy of tibiofemoral kinematics estimated by the ultrasound tracking system and by a conventional skin-mounted marker based motion capture system in a cadaveric experimental scenario. Two typical tibiofemoral joint models (spherical and hinge models) were used to derive relevant kinematic outcomes. Intra-cortical bone pins equipped with optical markers were inserted in the tibial and femoral bones to serve as a reference to provide ground truth kinematics. The ultrasound tracking system resulted in lower kinematic errors than the skin-mounted markers (the ultrasound tracking system: maximum root-mean-square (RMS) error 3.44° for rotations and 4.88 mm for translations, skin-mounted markers with the spherical joint model: 6.32° and 6.26 mm, the hinge model: 6.38° and 6.52 mm). Our proposed ultrasound tracking system has the potential of measuring direct bone kinematics, thereby mitigating the influence and propagation of STA. Consequently, this technique could be considered as an alternative method for measuring 6-DOF tibiofemoral kinematics, which may be adopted in gait analysis and clinical practice.

Inducing targeted failure in cadaveric testing of 3-segment spinal units with and without simulated metastases.

We propose an experimental setup and protocol able to induce targeted failure of the middle vertebra in 3-segment spinal units and to capture the specimens' deformation in their post-failure state. Sixteen 3-segment spinal units with and without artificial metastases were destructively tested in axial compression using one of two failure criteria; either: (A) A clear drop in force (>10-15% of peak force) (n = 4); or (B) A minimum displacement of 5 mm (n = 12). Subsequently, the specimens were fully fixated in polymethylmethacrylate (PMMA), thereby preserving their post-failure state. Pre- and post-experiment computed tomography (CT) scans were acquired to determine the occurrence of failure in one of the vertebral bodies. All specimens were successfully fixated in their post-failure state. When applying failure criterion A, two specimens showed signs of failure. When applying failure criterion B, all specimens showed signs of failure; in 9 out of 12 specimens this occurred in the middle vertebrae only. In conclusion, this research provides an experimental setup and protocol able to induce targeted failure of 3-segment spinal units and to capture the specimens' deformation in their post-failure state. Furthermore, this study illustrates the importance of an adequate failure criterion for successful simulation of vertebral fractures in an experimental setup.

Coupled motions in human and porcine thoracic and lumbar spines.

Coupled motions, i.e., motions along axes other than the loaded axis, have been reported to occur in the human spine, and are likely to be influenced by inclined local axes due to the sagittal plane spine curvature. Furthermore, the role of facet joints in such motions is as yet unclear. Therefore, this study aimed at assessing coupled motions in multiple spine sections in vitro, before and after removal of posterior elements. Six elderly human and 6 young porcine spines were sectioned in four segments (high thoracic, mid thoracic, low thoracic and lumbar), each consisting of four vertebrae and three intervertebral discs. Segments were loaded along each of the three axes, and three-dimensional rotations of the middle segment were quantified. Subsequently, posterior elements were removed and the protocol was repeated. To avoid mixed loading between Axial Rotation (AR) and Lateral Bending (LB), in contrast to other studies, local axes at the vertebrae were defined as aligned with the loading device prior to each load application. Expressed as a percentage of motion in the loaded direction, coupled motions were on average larger in human (22.7%, SD = 2.2%) than in porcine (11.9%, SD = 1.2%) spines (p < .001). Largest coupled motions were obtained in AR loading of the lumbar spine segments, with mean magnitudes averaged over coupling axes for human L2-L3 joints of 48.9% (SD = 13.2%), including somewhat more LB (56.4%, SD = 18.6) than FE (41.4%, SD = 14.1%) coupling. For porcine L3-L4 joints average coupling in AR loading was 29.3% (SD = 8.2%). In human segments removal of posterior elements only had substantial effects in the lumbar spine segments, where posterior element removal decreased coupled motion during AR loading, averaged over LB and FE coupling, from 48.9% (SD = 13.2%) to 27.7% (SD = 6.1%), mainly through increased motion in the loaded direction. The present results indicate that coupled motions were largest in the lumbar spine. In human spines, posterior elements only contributed to coupled motions in lumbar axial rotation loading.

Twente spine model: A complete and coherent dataset for musculo-skeletal modeling of the thoracic and cervical regions of the human spine.

Musculo-skeletal modeling could play a key role in advancing our understanding of the healthy and pathological spine, but the credibility of such models are strictly dependent on the accuracy of the anatomical data incorporated. In this study, we present a complete and coherent musculo-skeletal dataset for the thoracic and cervical regions of the human spine, obtained through detailed dissection of an embalmed male cadaver. We divided the muscles into a number of muscle-tendon elements, digitized their attachments at the bones, and measured morphological muscle parameters. In total, 225 muscle elements were measured over 39 muscles. For every muscle element, we provide the coordinates of its attachments, fiber length, tendon length, sarcomere length, optimal fiber length, pennation angle, mass, and physiological cross-sectional area together with the skeletal geometry of the cadaver. Results were consistent with similar anatomical studies. Furthermore, we report new data for several muscles such as rotatores, multifidus, levatores costarum, spinalis, semispinalis, subcostales, transversus thoracis, and intercostales muscles. This dataset complements our previous study where we presented a consistent dataset for the lumbar region of the spine (Bayoglu et al., 2017). Therefore, when used together, these datasets enable a complete and coherent dataset for the entire spine. The complete dataset will be used to develop a musculo-skeletal model for the entire human spine to study clinical and ergonomic applications.

Twente spine model: A complete and coherent dataset for musculo-skeletal modeling of the lumbar region of the human spine.

Musculo-skeletal modeling can greatly help in understanding normal and pathological functioning of the spine. For such models to produce reliable muscle and joint force estimations, an adequate set of musculo-skeletal data is necessary. In this study, we present a complete and coherent dataset for the lumbar spine, based on medical images and dissection measurements from one embalmed human cadaver. We divided muscles into muscle-tendon elements, digitized their attachments at the bones and measured morphological parameters. In total, we measured 11 muscles from one body side, using 96 elements. For every muscle element, we measured three-dimensional coordinates of its attachments, fiber length, tendon length, sarcomere length, optimal fiber length, pennation angle, mass, and physiological cross-sectional area together with the geometry of the lumbar spine. Results were consistent with other anatomical studies and included new data for the serratus posterior inferior muscle. The dataset presented in this paper enables a complete and coherent musculo-skeletal model for the lumbar spine and will improve the current state-of-the art in predicting spinal loading.

Spinal shape modulation in a porcine model by a highly flexible and extendable non-fusion implant system.

PURPOSE: In vivo evaluation of scoliosis treatment using a novel approach in which two posterior implants are implanted: XSLAT (eXtendable implant correcting Scoliosis in LAT bending) and XSTOR (eXtendable implant correcting Scoliosis in TORsion). The highly flexible and extendable implants use only small, but continuous lateral forces (XSLAT) and torques (XSTOR), thereby allowing growth and preventing fusion. METHODS: Since (idiopathic) scoliosis does not occur spontaneously in animals, the device was used to induce a spinal deformity rather than correct it. Six of each implants were tested for their ability to induce scoliotic deformations in 12 growing pigs. Each implant spanned six segments and was attached to three vertebrae using sliding anchors. Radiological and histological assessments were done throughout the 8-week study. RESULTS: In all animals, the intended deformation was accomplished. Average Cobb angles were 19° for XSLAT and 6° for XSTOR. Average apical spinal torsion was 0° for XSLAT and 9° for XSTOR. All instrumented segments remained mobile and showed 20 % growth. Moderate degeneration of the facet joints was observed and some debris was found in the surrounding tissue. CONCLUSIONS: The approach accomplished the intended spinal deformation while allowing growth and preventing fusion.

A novel anchoring system for use in a nonfusion scoliosis correction device.

BACKGROUND CONTEXT: Insertion of a pedicle screw in the mid- and high thoracic regions has a serious risk of facet joint damage. Because flexible implant systems require intact facet joints, we developed an enhanced fixation that is less destructive to spinal structures. The XSFIX is a posterior fixation system that uses cables that are attached to the transverse processes of a vertebra. PURPOSE: To determine whether a fixation to the transverse process using the XSFIX is strong enough to withstand the loads applied by the XSLATOR (a novel, highly flexible nonfusion implant system) and thus, whether it is a suitable alternative for pedicle screw fixation. STUDY DESIGN: The strength of a novel fixation system using transverse process cables was determined and compared with the strength of a similar fixation using polyaxial pedicle screws on different vertebral levels. METHODS: Each of the 58 vertebrae, isolated from four adult human cadavers, was instrumented with either a pedicle screw anchor (PSA) system or a prototype of the XSFIX. The PSA consisted of two polyaxial pedicle screws and a 5 mm diameter rod. The XSFIX prototype consisted of two bodies that were fixed to the transverse processes, interconnected with a similar rod. Each fixation system was subjected to a lateral or an axial torque. RESULTS: The PSA demonstrated fixation strength in lateral loading and torsion higher than required for use in the XSLATOR. The XSFIX demonstrated high enough fixation strength (in both lateral loading and torsion), only in the high and midthoracic regions (T10-T12). CONCLUSIONS: This experiment showed that the fixation strength of XSFIX is sufficient for use with the XSLATOR only in mid- and high thoracic regions. For the low thoracic and lumbar region, the PSA is a more rigid fixation. Because the performance of the new fixation system appears to be favorable in the high and midthoracic regions, a clinical study is the next challenge.

Prophylactic vertebroplasty can decrease the fracture risk of adjacent vertebrae: an in vitro cadaveric study.

Adjacent level vertebral fractures are common in patients with osteoporotic wedge fractures, but can theoretically be prevented with prophylactic vertebroplasty. Previous tests on prophylactic vertebroplasties have been performed under axial loading, while in vivo changes in spinal alignment likely cause off-axis loads. In this study we determined whether prophylactic vertebroplasty can also reduce the fracture risk under off-axis loads. In a previous study, we tested vertebral bodies that were loaded axially or 20° off-axis representing vertebrae in an unfractured spine or vertebrae adjacent to a wedge fracture, respectively. In the current study, vertebral failure load and stiffness of our previously tested vertebral bodies were compared to those of a new group of vertebral bodies that were filled with bone cement and then loaded 20° off-axis. These vertebral bodies represented adjacent-level vertebrae with prophylactic bone cement filling. Prophylactic augmentation resulted in failure loads that were comparable to those of the 0° group, and 32% greater than the failure loads of the 20° group. The stiffness of the prophylacticly augmented vertebrae was 21% lower than that of the 0° group, but 27% higher than that of the 20° group. We conclude that prophylactic augmentation can decrease the fracture risk in a malaligned, osteoporotic vertebra. Whether this is enough to actually prevent additional vertebral fractures in vivo remains subject of further study.

Does bone cement in percutaneous vertebroplasty act as a stress riser?

STUDY DESIGN: An in vitro cadaveric study. OBJECTIVE: To determine whether percutaneous vertebroplasty (PVP) with a clinically relevant amount of bone cement is capable of causing stress peaks in adjacent-level vertebrae. SUMMARY OF BACKGROUND DATA: It is often suggested that PVP of a primary spinal fracture causes stress peaks in adjacent vertebrae, thereby leading to additional fractures. The in vitro studies that demonstrated this relationship, however, use bigger volumes of bone cement used clinically. METHODS: Ten fresh-frozen vertebrae were loaded until failure, while registering force and displacement as well as the pressure under the lower endplate. After failure, the vertebrae were augmented with clinically relevant amounts of bone cement and then again loaded until failure. The force, displacement, and pressure under the lower endplate were again registered. RESULTS: Stress peaks were not related to the location of the injected bone cement. Both failure load and stiffness were significantly lower after augmentation. CONCLUSION: On the basis of our findings, we conclude that vertebral augmentation with clinically relevant amounts of bone cement does not lead to stress peaks under the endplate. It is therefore unlikely that PVP, in itself, causes detrimental stresses in the adjacent vertebrae, leading to new vertebral fractures. LEVEL OF EVIDENCE: N/A.

Posteriorly directed shear loads and disc degeneration affect the torsional stiffness of spinal motion segments: a biomechanical modeling study.

STUDY DESIGN: Finite element study. OBJECTIVE: To analyze the effects of posterior shear loads, disc degeneration, and the combination of both on spinal torsion stiffness. SUMMARY OF BACKGROUND DATA: Scoliosis is a 3-dimensional deformity of the spine that presents itself mainly in adolescent girls and elderly patients. Our concept of its etiopathogenesis is that an excess of posteriorly directed shear loads, relative to the body's intrinsic stabilizing mechanisms, induces a torsional instability of the spine, making it vulnerable to scoliosis. Our hypothesis for the elderly spine is that disc degeneration compromises the stabilizing mechanisms. METHODS: In an adult lumbar motion segment model, the disc properties were varied to simulate different aspects of disc degeneration. These models were then loaded with a pure torsion moment in combination with either a shear load in posterior direction, no shear, or a shear load in anterior direction. RESULTS: Posteriorly directed shear loads reduced torsion stiffness, anteriorly directed shear loads increased torsion stiffness. These effects were mainly caused by a later (respectively earlier) onset of facet joint contact. Disc degeneration cases with a decreased disc height that leads to slackness of the annular fibers and ligaments caused a significantly decreased torsional stiffness. The combination of this stage with posterior shear loading reduced the torsion stiffness to less than half the stiffness of a healthy disc without shear loads. The end stage of disc degeneration increased torsion stiffness again. CONCLUSION: The combination of a decreased disc height, that leads to slack annular fibers and ligaments, and posterior shear loads very significantly affects torsional stiffness: reduced to less than half the stiffness of a healthy disc without shear loads. Disc degeneration, thus, indeed compromises the stabilizing mechanisms of the elderly spine. A combination with posteriorly directed shear loads could then make it vulnerable to scoliosis. LEVEL OF EVIDENCE: N/A.

Toward a more realistic prediction of peri-prosthetic micromotions.

The finite element (FE) method has become a common tool to evaluate peri-prosthetic micromotions in cementless total hip arthroplasty. Often, only the peak joint load and a selected number of muscle loads are applied to determine micromotions. Furthermore, the applied external constraints are simplified (diaphyseal fixation), resulting in a non-physiological situation. In this study, a scaled musculoskeletal model was used to extract a full set of muscle and hip joint loads occurring during a walking cycle. These loads were applied incrementally to an FE model to analyze micromotions. The relation between micromotions and external loads was investigated, and how micromotions during a full loading cycle compared to those calculated when applying a peak load only. Finally, the effect of external constraints was analyzed (full model vs. diaphyseal fixation and reduced number of muscle loads). Relatively large micromotions were found during the swing phase when the hip joint forces were relatively low. Maximal micromotions, however, did concur with the peak hip joint force. Applying only a peak joint force resulted in peak micromotions similar to those found when full walking cycle loads were applied. The magnitude and direction of the micromotions depended on the applied muscle loads, but not on external constraints.

Can vertebral density changes be explained by intervertebral disc degeneration?

One of the major problems facing the elderly spine is the occurrence of vertebral fractures due to low bone mass. Although typically attributed to osteoporosis, disc degeneration has also been suggested to play a role in vertebral fractures. Existing bone adaptation theories and simulations may explain the biomechanical pathway from a degenerated disc to an increased fracture risk. A finite element model of a lumbar segment was created and calibrated. Subsequently the disc properties were varied to represent either a healthy or degenerated disc and the resulting bone adaptation was simulated. Disc degeneration resulted in a shift of load from the nucleus to the annulus. The resulting bone adaptation led to a dramatically reduced density of the trabecular core and to an increased density in the vertebral walls. Degeneration of just the nucleus, and in particular the dehydration of the nucleus, resulted in most of this bone density change. Additional annulus degeneration had much less of an effect on the density values. The density decrease in the trabecular core as seen in this study matches clinical observations. Therefore, bone remodeling theories can assists in explaining the potential synergistic effects of disc degeneration and osteoporotis in the occurrence of vertebral fractures.

Medial branch blocks and fresh vertebral compression fractures.

Article: Bogduk Nikolai, MacVicar J, Borowczyk J. The pain of vertebral compression fractures can arise in the posterior elements. Pain Med 2010;11:1666-73.

Influence of interpersonal geometrical variation on spinal motion segment stiffness: implications for patient-specific modeling.

STUDY DESIGN: A validated finite element model of an L3-L4 motion segment is used to analyze the effects of interpersonal differences in geometry on spinal stiffness. OBJECTIVE: The objective of this study is to determine which of the interpersonal variations of the geometry of the spine have a large effect on spinal stiffness. This will improve patient-specific modeling. SUMMARY OF BACKGROUND DATA: The parameters that define the geometry of a motion segment are vertebral height, disc height, endplate width, endplate depth, spinous process length, transverse process width, nucleus size, lordosis angle, facet area, facet orientation, and the cross-sectional areas of the ligaments. All these parameters differ between patients. The influence of each parameter on spinal stiffness is largely unknown and such knowledge would greatly help in patient-specific modeling of the spine. METHODS: The range of interpersonal variation of each of the geometric parameters was set at mean±2SD (covering 95% of the population). Subsequently, we determined the effect of each of these ranges on the bending stiffness in flexion, extension, axial rotation, and lateral bending. RESULTS: Disc height had the largest influence; a maximal disc height reduced the spinal stiffness to 75-86% of the mean motion segment stiffness, and a minimal disc height increased the spinal stiffness to 154-226% of the mean motion segment stiffness. Lordosis angle, transversal and longitudinal facet angle, endplate depth, and area of the capsular ligament also had a substantial influence (>5%) on the stiffness, but considerable less than the influence of the disc height. Ligament areas, nucleus size, spinous process length, and length of processes are of negligible effect (<2%) on the stiffness. CONCLUSION: The disc height should be accurately determined in patients to estimate the spinal stiffness. Ligament areas, nucleus size, spinous process length, and transverse process width do not need patient-specific modeling.