Towards determining soft tissue properties for modelling spine surgery : current progress and challenges

Current complication rates for adolescent scoliosis surgery necessitate the development of better surgical planning tools to improve outcomes. Here we present our approach to developing finite element models of the thoracolumbar spine for deformity surgery simulation, with patient-specific model anatomy based on low-dose pre-operative computed tomography scans. In a first step towards defining patient-specific tissue properties, an initial 'benchmark' set of properties were used to simulate a clinically performed pre-operative spinal flexibility assessment, the fulcrum bending radiograph. Clinical data for ten patients were compared with the simulated results for this assessment and in cases where these data differed by more than 10%, soft tissue properties for the costo-vertebral joint (CVJt) were altered to achieve better agreement. Results from these analyses showed that changing the CVJt stiffness resulted in acceptable agreement between clinical and simulated flexibility in two of the six cases. In light of these results and those of our previous studies in this area, it is suggested that spinal flexibility in the fulcrum bending test is not governed by any single soft tissue structure acting in isolation. More detailed biomechanical characterisation of the fulcrum bending test is required to provide better data for determination of patient-specific soft tissue properties.

Segmental torso masses and gravity-induced coronal plane joint moments in adolescent idiopathic scoliosis

Introduction: Calculating segmental (vertebral level-by-level) torso masses in Adolescent Idiopathic Scoliosis (AIS) patients allows the gravitational loading on the scoliotic spine during relaxed standing to be estimated. This study used supine CT scans of AIS patients to measure segmental torso masses and explored the joint moments in the coronal plane, particularly at the apex of a scoliotic major curve. Methods: Existing low dose CT data from the Paediatric Spine Research Group was used to calculate vertebral level-by-level torso masses and joint moments occurring in the spine for a group of 20 female AIS patients with right sided thoracic curves. The mean age was 15.0 ± 2.7 years and all curves were classified Lenke Type 1 with a mean Cobb angle 52 ± 5.9°. Image processing software, ImageJ (v1.45 NIH USA) was used to create reformatted coronal plane images, reconstruct vertebral level-by-level torso segments and subsequently measure the torso volume corresponding to each vertebral level. Segment mass was then determined by assuming a tissue density of 1.04x103 kg/m3. Body segment masses for the head, neck and arms were taken from published anthropometric data (Winter 2009). Intervertebral joint moments in the coronal plane at each vertebral level were found from the position of the centroid of the segment masses relative to the joint centres with the segmental body mass data. Results and Discussion: The magnitude of the torso masses from T1-L5 increased inferiorly, with a 150% increase in mean segmental torso mass from 0.6kg at T1 to 1.5kg at L5. The magnitudes of the calculated coronal plane joint moments during relaxed standing were typically 5-7 Nm at the apex of the curve, with the highest apex joint torque of 7Nm. The CT scans were performed in the supine position and curve magnitudes are known to be 7-10° smaller than those measured in standing, due to the absence of gravity acting on the spine. Hence, it can be expected that the moments produced by gravity in the standing individual will be greater than those calculated here.

Geometric sensitivity of patient-specific finite element models of the spine to variability in user-selected anatomical landmarks

Software to create individualised finite element (FE) models of the osseoligamentous spine using pre-operative computed tomography (CT) data-sets for spinal surgery patients has recently been developed. This study presents a geometric sensitivity analysis of this software to assess the effect of intra-observer variability in user-selected anatomical landmarks. User-selected landmarks on the osseous anatomy were defined from CT data-sets for three scoliosis patients and these landmarks were used to reconstruct patient-specific anatomy of the spine and ribcage using parametric descriptions. The intra-observer errors in landmark co-ordinates for these anatomical landmarks were calculated. FE models of the spine and ribcage were created using the reconstructed anatomy for each patient and these models were analysed for a loadcase simulating clinical flexibility assessment. The intra-observer error in the anatomical measurements was low in comparison to the initial dimensions, with the exception of the angular measurements for disc wedge and zygapophyseal joint (z-joint) orientation and disc height. This variability suggested that CT resolution may influence such angular measurements, particularly for small anatomical features, such as the z-joints, and may also affect disc height. The results of the FE analysis showed low variation in the model predictions for spinal curvature with the mean intra-observer variability substantially less than the accepted error in clinical measurement. These findings demonstrate that intra-observer variability in landmark point selection has minimal effect on the subsequent FE predictions for a clinical loadcase.

The stabilising role of the costovertebral joints and spinous processes on thoracic spine stability

Introduction The importance of in vitro biomechanical testing in today’s understanding of spinal pathology and treatment modalities cannot be stressed enough. Different studies have used differing levels of dissection of their spinal segments for their testing protocols[1, 2]. The aim of this study was to assess the impact of removing the costovertebral joints and partial resection of the spinous process sequentially, on the stiffness of the immature thoracic bovine spinal segment. Materials and Methods Thoracic spines from 6-8 week old calves were used. Each spine was dissected and divided into motion segments with 5cm of attached rib on each side and full spinous processes including levels T4-T11 (n=28). They were potted in polymethylemethacrylate. An Instron Biaxial materials testing machine with a custom made jig was used for testing. The segments were tested in flexion/extension, lateral bending and axial rotation at 37⁰C and 100% humidity, using moment control to a maximum 1.75 Nm with a loading rate of 0.3 Nm per second. They were first tested intact for ten load cycles with data collected from the tenth cycle. Progressive dissection was performed by removing first the attached ribs, followed by the spinous process at its base. Biomechanical testing was carried out after each level of dissection using the same protocol. Statistical analysis of the data was performed using repeated measures ANOVA. Results In combined flexion/extension there was a significant reduction in stiffness of 16% (p=0.002). This was mainly after resection of the ribs (14%, p=0.024) and mainly occurred in flexion where stiffness reduced by 22% (p=0.021). In extension, stiffness dropped by 13% (p=0.133). However there was no further significant change in stiffness on resection of the spinous process (<1%) (p=1.00). In lateral bending there was a significant decrease in stiffness of 13% (p<0.001). This comprised a drop of 11% on resection of the ribs (p=0.009) and a further 8% on resection of the spinous process (p=0.014). There was no difference between left and right bending. In axial rotation there was no significant change in stiffness after each stage of dissection (p=0.253). There was no difference between left and right rotation. Conclusion The costovertebral joints play a significant role in providing stability to the bovine thoracic spine in both flexion/extension and lateral bending, whereas the spinous processes play a minor role. Both elements have little effect on axial rotation stability.

Comparison of eight published static finite element models of the intact lumbar spine : predictive power of models improves when combined together

Finite element (FE) model studies have made important contributions to our understanding of functional biomechanics of the lumbar spine. However, if a model is used to answer clinical and biomechanical questions over a certain population, their inherently large inter-subject variability has to be considered. Current FE model studies, however, generally account only for a single distinct spinal geometry with one set of material properties. This raises questions concerning their predictive power, their range of results and on their agreement with in vitro and in vivo values. Eight well-established FE models of the lumbar spine (L1-5) of different research centres around the globe were subjected to pure and combined loading modes and compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges, and their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with published median in vitro values. However, the ranges of predictions were larger and exceeded those reported in vitro, especially for the facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with measured in vivo values. In light of high inter-subject variability, the generalization of results of a single model to a population remains a concern. This study demonstrated that the pooled median of individual model results, similar to a probabilistic approach, can be used as an improved predictive tool in order to estimate the response of the lumbar spine.

Comparison of eight published lumbar spine finite element models

Due to its ability to represent intricate systems with material nonlinearities as well as irregular loading, boundary, geometrical and material domains, the finite element (FE) method has been recognized as an important computational tool in spinal biomechanics. Current FE models generally account for a single distinct spinal geometry with one set of material properties despite inherently large inter-subject variability. The uncertainty and high variability in tissue material properties, geometry, loading and boundary conditions has cast doubt on the reliability of their predictions and comparability with reported in vitro and in vivo values. A multicenter study was undertaken to compare the results of eight well-established models of the lumbar spine that have been developed, validated and applied for many years. Models were subjected to pure and combined loading modes and their predictions were compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges; their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with previously published median in vitro values. However, the ranges of predictions were larger and exceeded the in vitro ranges, especially for facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with in vivo values. The simulations yielded median facet joint forces of 0 N in flexion, 38 N in extension, 14 N in lateral bending and 60 N in axial rotation that could not be validated due to the paucity of in vivo facet joint forces. In light of high inter-subject variability, one must be cautious when generalizing predictions obtained from one deterministic model. This study demonstrates however that the predictive power increases when FE models are combined together. The median of individual numerical results can hence be used as an improved tool in order to estimate the response of the lumbar spine.

Compressive and shear force generated in the lumbar spine of female rowers

Rowers have and accrue greater lumbar spine bone mineral density (BMD) associated with mechanical loading produced during rowing. The aim of this study was to estimate the mechanical loading generated at the lumbar spine (LS) that is apparently providing an osteogenic benefit. The cohort comprised 14 female rowers (average age: 19.7yrs; height: 170.9 cm, weight: 59.5 kg) and 14 female matched controls (average age: 20.9 m yrs; height: 167.5 cm; weight: 58.1 kg). BMD was assessed using the Hologic QDR 2000+ bone densitometer, indicating higher lumbar spine BMD in the rowers compared to the control subjects (1,069 +/- 0.1 vs. 1,027 +/- 0.1 g/cm2). No significant difference existed for BMD at any other site. All rowers performed a six-minute simulated race on a Concept II rowing ergometer. Mechanical loading generated at the lumbar spine during this task was assessed using a two-dimensional model of the spine, enabling the calculation of the compressive and shear forces at L4/L5. The shear force was the joint reaction force perpendicular to the spine at the L4/L5 joint. Peak compressive and shear force at the lumbar spine of the rowers were 2,694 +/- 609 (N) and 660 +/- 117 (N), respectively. Peak compressive force at the LS relative to body weight was 4.6 times body weight. The literature would suggest that forces of this magnitude, generated at the LS during maximal rowing, may be contributing to the site specific higher LS BMD found in the rowers.

Spinal inflammation in the absence of sacroiliac joint inflammation on magnetic resonance imaging in patients with active nonradiographic axial spondyloarthritis

Objective: To evaluate the presence of spinal inflammation with and without sacroiliac (SI) joint inflammation on magnetic resonance imaging (MRI) in patients with active nonradiographic axial spondyloarthritis (SpA), and to compare the disease characteristics of these subgroups. Methods: ABILITY-1 is a multicenter, randomized, controlled trial of adalimumab versus placebo in patients with nonradiographic axial SpA classified using the Assessment of SpondyloArthritis international Society axial SpA criteria. Baseline MRIs were centrally scored independently by 2 readers using the Spondyloarthritis Research Consortium of Canada (SPARCC) method for the SI joints and the SPARCC 6-discovertebral unit method for the spine. Positive evidence of inflammation on MRI was defined as a SPARCC score of >2 for either the SI joints or the spine. Results: Among patients with baseline SPARCC scores, 40% had an SI joint score of >2 and 52% had a spine score of >2. Forty-nine percent of patients with baseline SI joint scores of <2, and 58% of those with baseline SI joint scores of >2, had a spine score of >2. Comparison of baseline disease characteristics by baseline SI joint and spine scores showed that a greater proportion of patients in the subgroup with a baseline SPARCC score of >2 for both SI joints and spine were male, and patients with spine and SI joint scores of <2 were younger and had shorter symptom duration. SPARCC spine scores correlated with baseline symptom duration, and SI joint scores correlated negatively with the baseline Bath Ankylosing Spondylitis Disease Activity Index, but neither correlated with the baseline Ankylosing Spondylitis Disease Activity Score, total back pain, the patient's global assessment of disease activity, the Bath Ankylosing Spondylitis Functional Index, morning stiffness, nocturnal pain, or C-reactive protein level. Conclusion: Assessment by experienced readers showed that spinal inflammation on MRI might be observed in half of patients with nonradiographic axial SpA without SI joint inflammation.

Understanding how axial loads on the spine influence segmental biomechanics for idiopathic scoliosis patients: A magnetic resonance imaging study

Background Segmental biomechanics of the scoliotic spine are important since the overall spinal deformity is comprised of the cumulative coronal and axial rotations of individual joints. This study investigates the coronal plane segmental biomechanics for adolescent idiopathic scoliosis patients in response to physiologically relevant axial compression. Methods Individual spinal joint compliance in the coronal plane was measured for a series of 15 idiopathic scoliosis patients using axially loaded magnetic resonance imaging. Each patient was first imaged in the supine position with no axial load, and then again following application of an axial compressive load. Coronal plane disc wedge angles in the unloaded and loaded configurations were measured. Joint moments exerted by the axial compressive load were used to derive estimates of individual joint compliance. Findings The mean standing major Cobb angle for this patient series was 46°. Mean intra-observer measurement error for endplate inclination was 1.6°. Following loading, initially highly wedged discs demonstrated a smaller change in wedge angle, than less wedged discs for certain spinal levels (+ 2,+1,− 2 relative to the apex, (p < 0.05)). Highly wedged discs were observed near the apex of the curve, which corresponded to lower joint compliance in the apical region. Interpretation While individual patients exhibit substantial variability in disc wedge angles and joint compliance, overall there is a pattern of increased disc wedging near the curve apex, and reduced joint compliance in this region. Approaches such as this can provide valuable biomechanical data on in vivo spinal biomechanics of the scoliotic spine, for analysis of deformity progression and surgical planning.

Gravity-induced coronal plane joint moments in adolescent idiopathic scoliosis

Background Adolescent Idiopathic Scoliosis is the most common type of spinal deformity, and whilst the risk of progression appears to be biomechanically mediated (larger deformities are more likely to progress), the detailed biomechanical mechanisms driving progression are not well understood. Gravitational forces in the upright position are the primary sustained loads experienced by the spine. In scoliosis they are asymmetrical, generating moments about the spinal joints which may promote asymmetrical growth and deformity progression. Using 3D imaging modalities to estimate segmental torso masses allows the gravitational loading on the scoliotic spine to be determined. The resulting distribution of joint moments aids understanding of the mechanics of scoliosis progression. Methods Existing low-dose CT scans were used to estimate torso segment masses and joint moments for 20 female scoliosis patients. Intervertebral joint moments at each vertebral level were found by summing the moments of each of the torso segment masses above the required joint. Results The patients’ mean age was 15.3 years (SD 2.3; range 11.9 – 22.3 years); mean thoracic major Cobb angle 52° (SD 5.9°; range 42°-63°) and mean weight 57.5 kg (SD 11.5 kg; range 41 – 84.7 kg). Joint moments of up to 7 Nm were estimated at the apical level. No significant correlation was found between the patients’ major Cobb angles and apical joint moments. Conclusions Patients with larger Cobb angles do not necessarily have higher joint moments, and curve shape is an important determinant of joint moment distribution. These findings may help to explain the variations in progression between individual patients. This study suggests that substantial corrective forces are required of either internal instrumentation or orthoses to effectively counter the gravity-induced moments acting to deform the spinal joints of idiopathic scoliosis patients.

Does flexibility of the cervical and thoracic spine influence forward head posture in healthy individuals?

Background: Forward head posture (FHP) is a commonly reported deviation from a neutral neck posture, usually implying a protracted head position in the sagittal plane. Habitual FHP has been associated with a higher incidence of painful neck disorders, changes in joint mobility and muscle behaviour within the cervicothoracic regions. One factor that has received attention in the literature with regards to FHP is flexibility of the neck. A number of previous studies have previously investigated the relationship between neck flexibility and neck posture under different conditions, but at present this relationship is unclear.

Spinal circuits can accommodate interaction torques during multijoint limb movements

The dynamic interaction of limb segments during movements that involve multiple joints creates torques in one joint due to motion about another. Evidence shows that such interaction torques are taken into account during the planning or control of movement in humans. Two alternative hypotheses could explain the compensation of these dynamic torques. One involves the use of internal models to centrally compute predicted interaction torques and their explicit compensation through anticipatory adjustment of descending motor commands. The alternative, based on the equilibrium-point hypothesis, claims that descending signals can be simple and related to the desired movement kinematics only, while spinal feedback mechanisms are responsible for the appropriate creation and coordination of dynamic muscle forces. Partial supporting evidence exists in each case. However, until now no model has explicitly shown, in the case of the second hypothesis, whether peripheral feedback is really sufficient on its own for coordinating the motion of several joints while at the same time accommodating intersegmental interaction torques. Here we propose a minimal computational model to examine this question. Using a biomechanics simulation of a two-joint arm controlled by spinal neural circuitry, we show for the first time that it is indeed possible for the neuromusculoskeletal system to transform simple descending control signals into muscle activation patterns that accommodate interaction forces depending on their direction and magnitude. This is achieved without the aid of any central predictive signal. Even though the model makes various simplifications and abstractions compared to the complexities involved in the control of human arm movements, the finding lends plausibility to the hypothesis that some multijoint movements can in principle be controlled even in the absence of internal models of intersegmental dynamics or learned compensatory motor signals.

Coupling between muscle activities and muscle torques during horizontal-planar arm movements with direction reversal

In this study we investigated the hypothesis that the simple set of rules used to explain the modulation of muscle activities during single-joint movements could also be applied for reversal movements of the shoulder and elbow joints. The muscle torques of both joints were characterized by a triphasic impulse. The first impulse of each joint accelerated the limb to the target and was generated by an initial burst of the muscles activated first (primary mover). The second impulse decelerated the limb to the target, reversed movement direction and accelerated the limb back to the initial position, and was generated by an initial burst of the muscles activated second (secondary movers). A third impulse, in each joint, decelerated the limb to the initial position due to the generation of a second burst of the primary movers. The first burst of the primary mover decreased abruptly, and the latency between the activation of the primary and secondary movers varied in proportion with target distances for the elbow, but not for the shoulder muscles. All impulses and bursts increased with target distances and were well coupled. Therefore, as predicted, the bursts of muscle activities were modulated to generate the appropriate level of muscle torque. (C) 2005 Elsevier Ltd. All rights reserved.

Motor readiness and joint torque production in lower limbs of older women fallers and non-fallers

This study aimed to evaluate the motor response time and ability to develop joint torque at the knee and ankle in older women with and without a history of falls, in addition to investigating the effect of aging on these capacities. We assessed 18 young females, 21 older female fallers and 22 older female non-fallers. The peak torque, rate of torque development, rate of electromyography (EMG) rise, reaction time, premotor time and motor time were obtained through a dynamometric assessment and simultaneous electromyography. Surface EMGs of the rectus femoris (RF), vastus lateralis (VL), biceps femoris (BF), gastrocnemius lateralis (GL) and tibialis anterior (TA) muscles were recorded. Knee extension and flexion peak torques were lower in older fallers than in non-fallers. Knee extension and flexion and ankle plantarflexion and dorsiflexion peak torques were lower in both older groups than in the younger group. The rate of EMG rise of the BF and the motor time of the TA were lower and higher, respectively, in older fallers than in the younger adults. The time to reach peak torque in knee extension/flexion and ankle plantarflexion/dorsiflexion and the motor times of the RF, VL, BF and GL were higher in both older groups than in the younger groups. The motor time of the TA during ankle dorsiflexion and the knee extension peak torque were the major predictors of falls in older women, accounting for approximately 28% of the number of falls. Thus, these results further reveal the biomechanical parameters that affect the risk of falls and provide initial findings to support the prescription of exercises in fall prevention programs. © 2013 Elsevier Ltd.