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.

Development of a new ultrasound-based system for tracking motion of the human lumbar spine: Reliability, stability and repeatability during forward bending movement trials

The aim of this study was to develop a new method for quantifying intersegmental motion of the spine in an instrumented motion segment L4–L5 model using ultrasound image post-processing combined with an electromagnetic device. A prospective test–retest design was employed, combined with an evaluation of stability and within- and between-day intra-tester reliability during forward bending by 15 healthy male patients. The accuracy of the measurement system using the model was calculated to be ± 0.9° (standard deviation = 0.43) over a 40° range and ± 0.4 cm (standard deviation = 0.28) over 1.5 cm. The mean composite range of forward bending was 15.5 ± 2.04° during a single trial (standard error of the mean = 0.54, coefficient of variation = 4.18). Reliability (intra-class correlation coefficient = 2.1) was found to be excellent for both within-day measures (0.995–0.999) and between-day measures (0.996–0.999). Further work is necessary to explore the use of this approach in the evaluation of biomechanics, clinical assessments and interventions.

Functional and structural of the erector spinae muscle during isometric lumbar extension

Study Design Cross-sectional study. Objectives To compare erector spinae (ES) muscle fatigue between chronic non-specific lower back pain (CNLBP) sufferers and healthy subjects from a biomechanical perspective during fatiguing isometric lumbar extensions. Background Paraspinal muscle maximal contraction and fatigue are used as a functional predictor for disabilities. The simplest method to determine muscle fatigue is by evaluating the evolution during specific contractions, such as isometric contractions. There are no studies that evaluate the evolution of the ES muscle during fatiguing isometric lumbar extensions and analyse functional and architectural variables. Methods In a pre-calibrated system, participants performed a maximal isometric extension of the lumbar spine for 5 and 30 seconds. Functional variables (torque and muscle activation) and architecture (pennation angle and muscle thickness) were measured using a load cell, surface electromyography and ultrasound, respectively. The results were normalised and a reliability study of the ultrasound measurement was made. Results: The ultrasound measurements were highly reliable, with Cronbach’s alpha values ranging from 0.951 0.981. All measured variables shown significant differences before and after fatiguing isometric lumbar extension. Conclusion During a lumbar isometric extension test, architecture and functional variables of the ES muscle could be analised using ultrasound, surface EMG and load cell. In adition, during an endurance test, ES muscle suffers an acute effect on architectural and functional variables.

Ability to discriminate between healthy and low back pain sufferers using ultrasound during maximum lumbar extension

Objective To analyze the ability to discriminate between healthy individuals and individuals with chronic nonspecific low back pain (CNLBP) by measuring the relation between patient-reported outcomes and objective clinical outcome measures of the erector spinae (ES) muscles using an ultrasound during maximal isometric lumbar extension. Design Cross-sectional study with screening and diagnostic tests with no blinded comparison. Setting University laboratory. Participants Healthy individuals (n=33) and individuals with CNLBP (n=33). Interventions Each subject performed an isometric lumbar extension. With the variables measured, a discriminate analysis was performed using a value ≥6 in the Roland and Morris disability questionnaire (RMDQ) as the grouping variable. Then, a logistic regression with the functional and architectural variables was performed. A new index was obtained from each subject value input in the discriminate multivariate analysis. Main Outcome Measures Morphologic muscle variables of the ES muscle were measured through ultrasound images. The reliability of the measures was calculated through intraclass correlation coefficients (ICCs). The relation between patient-reported outcomes and objective clinical outcome measures was analyzed using a discriminate function from standardized values of the variables and an analysis of the reliability of the ultrasound measurement. Results The reliability tests show an ICC value >.95 for morphologic and functional variables. The independent variables included in the analysis explained 42% (P=.003) of the dependent variable variance. Conclusions The relation between objective variables (electromyography, thickness, pennation angle) and a subjective variable (RMDQ ≥6) and the capacity of this relation to identify CNLBP within a group of healthy subjects is moderate. These results should be considered by clinicians when treating this type of patient in clinical practice.

Comparison of silicate-substituted calcium phosphate (Actifuse®) with rhBMP-2 (Infuse®) in posterolateral instrumented lumbar fusion

Study Design This was a randomised controlled trial in patients with degenerative disc disease (DDD) who underwent instrumented posterolateral lumbar fusion (PLF) surgery. Objective The aim of this study was to assess the efficacy of the bone grafting substitute, silicate-substituted calcium phosphate (SiCaP) compared with bone morphogenetic protein (rhBMP-2) and to evaluate clinical outcomes over a period of two years. Methods Patients undergoing PLF surgery for DDD at a single centre were recruited and randomised to one of two groups; SiCaP (n=9) or rhBMP-2 (n=10). One patient withdrew prior to randomisation and another from the rhBMP-2 group after randomisation. The radiological and clinical outcomes were examined and compared. Fusion was assessed at 12 months with computed tomography (CT) and plain radiographs. Clinical outcomes were evaluated by recording measures of pain, quality of life, disability and neurological status from six weeks to two years postoperatively. Results In the SiCaP and rhBMP-2 groups, fusion was observed in 9/9 and 8/9 patients respectively. Pain and disability scores were reduced and quality of life increased in both groups. Leg pain, disability and satisfaction scores were similar between the groups at each postoperative time point, however, back pain was less at six weeks and quality of life was higher at six months in the SiCaP group than the rhBMP-2 group. Conclusions SiCaP and rhBMP-2 were comparable in terms of achieving successful bone growth and fusion. Both groups similarly alleviated pain and improved quality of life, neurological, satisfaction and return to work outcomes following PLF surgery.

A semiautomatic method to identify vertebral end plate lesions (Schmorl's nodes)

Background Context There are differences in definitions of end plate lesions (EPLs), often referred to as Schmorl’s nodes, that may, to some extent, account for the large range of reported prevalence (3.8 to 76%). Purpose To develop a technique to measure the size, prevalence and location of EPLs in a consistent manner. Study Design/Setting This study proposed a method using a detection algorithm which was applied to five adolescent females (average age 15.1 years, range 13.0 to 19.2 years) with idiopathic scoliosis (average major Cobb angle 60°, range 55 to 67°). Methods Existing low-dose, computed tomography scans were segmented semi-automatically to extract 3D morphology of each vertebral endplate. Any remaining attachments to the posterior elements of adjacent vertebrae or endplates were then manually sectioned. An automatic algorithm was used to determine the presence and position of EPLs. Results EPLs were identified in 15 of the 170 (8.8%) endplates analysed with an average depth of 3.1mm. 11/15 of the EPLs were seen in the lumbar spine. The algorithm was found to be most sensitive to changes in the minimum EPL gradient at the edges of the EPL. Conclusions This study describes an imaging analysis technique for consistent measurement of the prevalence, location and size of EPLs. The technique can be used to analyse large populations without observer errors in EPL definitions.

Identifying vertebral end plate lesions (schmorl's nodes) using a semi-automatic method

INTRODUCTION There is a large range in the reported prevalence of end plate lesions (EPLs), sometimes referred to as Schmorl's nodes in the general population (3.8-76%). One possible reason for this large range is the differences in definitions used by authors. Previous research has suggested that EPLs may potentially be a primary disturbance of growth plates that leads to the onset of scoliosis. The aim of this study was to develop a technique to measure the size, prevalence and location of EPLs on Computed Tomography (CT) images of scoliosis patients in a consistent manner. METHODS A detection algorithm was developed and applied to measure EPLs for five adolescent females with idiopathic scoliosis (average age 15.1 years, average major Cobb 60°). In this algorithm, the EPL definition was based on the lesion depth, the distance from the edge of the vertebral body and the gradient of the lesion edge. Existing low-dose, CT scans of the patients' spines were segmented semi-automatically to extract 3D vertebral endplate morphology. Manual sectioning of any attachments between posterior elements of adjacent vertebrae and, if necessary, endplates was carried out before the automatic algorithm was used to determine the presence and position of EPLs. RESULTS EPLs were identified in 15 of the 170 (8.8%) endplates analysed with an average depth of 3.1mm. 73% of the EPLs were seen in the lumbar spines (11/15). A sensitivity study demonstrated that the algorithm was most sensitive to changes in the minimum gradient required at the lesion edge. CONCLUSION An imaging analysis technique for consistent measurement of the prevalence, location and size of EPLs on CT images has been developed. Although the technique was tested on scoliosis patients, it can be used to analyse other populations without observer errors in EPL definitions.

Quantifying progressive anterior overgrowth in the thoracic vertebrae of adolescent idiopathic scoliosis patients

Study design Anterior and posterior vertebral body heights were measured from sequential MRI scans of adolescent idiopathic scoliosis (AIS) patients and healthy controls. Objective To measure changes in vertebral body height over time during scoliosis progression to assess how vertebral body height discrepancies change during growth. Summary of background data Relative anterior overgrowth has been proposed as a potential driver for AIS initiation and progression. This theory proposes that the anterior column grows faster, and the posterior column slower, in AIS patients when compared to healthy controls. There is disagreement in the literature as to whether the anterior vertebral body heights are proportionally greater than posterior vertebral body heights in AIS patients when compared to healthy controls. To some extent, these discrepancies may be attributed to methodological differences. Methods MRI scans of the major curve of 21 AIS patients (mean age 12.5 ± 1.4 years, mean Cobb 32.2 ± 12.8º) and between T4 and T12 of 21 healthy adolescents (mean age 12.1 ± 0.5 years) were captured for this study. Of the 21 AIS patients, 14 had a second scan on average 10.8 ± 4.7 months after the first. Anterior and posterior vertebral body heights were measured from the true sagittal plane of each vertebra such that anterior overgrowth could be quantified. Results The difference between anterior and posterior vertebral body height in healthy, non-scoliotic children was significantly greater than in AIS patients with mild to moderate scoliosis. However there was no significant relationship between the overall anterior-posterior vertebral body height difference in AIS and either severity of the curve or its progression over time. Conclusions Whilst AIS patients have a proportionally longer anterior column than non-scoliotic controls, the degree of anterior overgrowth was not related to the rate of progression or the severity of the scoliotic curve.

Does the addition of Clonidine to lumbar nerve root blocks improve outcome? A randomised, prospective, single-blinded controlled pilot study

Purpose To evaluate if adding clonidine to a standard nerve root block containing local anaesthetic and steroid improved the outcome of patients with severe lumbar nerve root pain secondary to MRI proven lumbar disc prolapse. Methods We undertook a single blind, prospective, randomised controlled trial evaluating 100 consecutive patients with nerve root pain secondary to lumbar disc prolapse undergoing trans-foraminal epidural steroid injection either with or without the addition of clonidine. 50 patients were allocated to each arm of the study. The primary outcome measure was the avoidance of a second procedure- repeat injection or micro-discectomy surgery. Secondary outcome measures were also studied: pain scores for leg and back pain using a visual analogue scale (VAS), the Roland Morris Disability Questionnaire (RMDQ) and the Measure Your Own Medical Outcome Profile (MYMOP). Follow up was carried out at 6 weeks, 6 months and 1 year. Results No serious complications occurred. Of the 50 patients who received the addition of clonidine, 56% were classified as successful injections, with no further intervention required, as opposed to 40% who received the standard injection. This difference did not reach statistical significance (p=0.109, chi-squared test). All secondary measures showed no statistically significant differences between the groups except curiously, the standard group who had been classified as successful had better leg pain relief than the clonidine group (p=0.026) at 1 year. Conclusions This pilot study has shown a 16% treatment effect with adding clonidine to lumbar nerve root blocks and that it is a safe injectate for this purpose.

Lateral bone density variations in the scoliotic spine

Adolescent Idiopathic Scoliosis (AIS) is the most common deformity of the spine, affecting 2-4% of the population. Previous studies have shown that the vertebrae in scoliotic spines undergo abnormal shape changes, however there has been little exploration of how scoliosis affects bone density distribution within the vertebrae. In this study, existing CT scans of 53 female idiopathic scoliosis patients with right-sided main thoracic curves were used to measure the lateral (right to left) bone density profile at mid-height through each vertebral body. Five key bone density profile measures were identified from each normalised bone density distribution, and multiple regression analysis was performed to explore the relationship between bone density distribution and patient demographics (age, height, weight, body mass index (BMI), skeletal maturity, time since Menarche, vertebral level, and scoliosis curve severity). Results showed a marked convex/concave asymmetry in bone density for vertebral levels at or near the apex of the scoliotic curve. At the apical vertebra, mean bone density at the left side (concave) cortical shell was 23.5% higher than for the right (convex) cortical shell, and cancellous bone density along the central 60% of the lateral path from convex to concave increased by 13.8%. The centre of mass of the bone density profile at the thoracic curve apex was located 53.8% of the distance along the lateral path, indicating a shift of nearly 4% toward the concavity of the deformity. These lateral bone density gradients tapered off when moving away from the apical vertebra. Multi-linear regressions showed that the right cortical shell peak bone density is significantly correlated with skeletal maturity, with each Risser increment corresponding to an increase in mineral equivalent bone density of 4-5%. There were also statistically significant relationships between patient height, weight and BMI, and the gradient of cancellous bone density along the central 60% of the lateral path. Bone density gradient is positively correlated with weight, and negatively correlated with height and BMI, such that at the apical vertebra, a unit decrease in BMI corresponds to an almost 100% increase in bone density gradient.

Generation of a 3D proximal femur shape from a single projection 2D radiographic image.

Summary Generalized Procrustes analysis and thin plate splines were employed to create an average 3D shape template of the proximal femur that was warped to the size and shape of a single 2D radiographic image of a subject. Mean absolute depth errors are comparable with previous approaches utilising multiple 2D input projections. Introduction Several approaches have been adopted to derive volumetric density (g cm-3) from a conventional 2D representation of areal bone mineral density (BMD, g cm-2). Such approaches have generally aimed at deriving an average depth across the areal projection rather than creating a formal 3D shape of the bone. Methods Generalized Procrustes analysis and thin plate splines were employed to create an average 3D shape template of the proximal femur that was subsequently warped to suit the size and shape of a single 2D radiographic image of a subject. CT scans of excised human femora, 18 and 24 scanned at pixel resolutions of 1.08 mm and 0.674 mm, respectively, were equally split into training (created 3D shape template) and test cohorts. Results The mean absolute depth errors of 3.4 mm and 1.73 mm, respectively, for the two CT pixel sizes are comparable with previous approaches based upon multiple 2D input projections. Conclusions This technique has the potential to derive volumetric density from BMD and to facilitate 3D finite element analysis for prediction of the mechanical integrity of the proximal femur. It may further be applied to other anatomical bone sites such as the distal radius and lumbar spine.

Biological performance of a polycaprolactone-based scaffold used as fusion cage device in a large animal model of spinal reconstructive surgery

A bioactive and bioresorbable scaffold fabricated from medical grade poly (epsilon-caprolactone) and incorporating 20% beta-tricalcium phosphate (mPCL–TCP) was recently developed for bone regeneration at load bearing sites. In the present study, we aimed to evaluate bone ingrowth into mPCL–TCP in a large animal model of lumbar interbody fusion. Six pigs underwent a 2-level (L3/4; L5/6) anterior lumbar interbody fusion (ALIF) implanted with mPCL–TCP þ 0.6 mg rhBMP-2 as treatment group while four other pigs implanted with autogenous bone graft served as control. Computed tomographic scanning and histology revealed complete defect bridging in all (100%) specimen from the treatment group as early as 3 months. Histological evidence of continuing bone remodeling and maturation was observed at 6 months. In the control group, only partial bridging was observed at 3 months and only 50% of segments in this group showed complete defect bridging at 6 months. Furthermore, 25% of segments in the control group showed evidence of graft fracture, resorption and pseudoarthrosis. In contrast, no evidence of graft fractures, pseudoarthrosis or foreign body reaction was observed in the treatment group. These results reveal that mPCL–TCP scaffolds could act as bone graft substitutes by providing a suitable environment for bone regeneration in a dynamic load bearing setting such as in a porcine model of interbody spine fusion.

Prediction of adjacent vertebral fracture risk associated with vertebroplasty : a computer-based simulation study

Vertebrplasty involved injecting cement into a fractured vertebra to provide stabilisation. There is clinical evidence to suggest however that vertebroplasty may be assocated with a higher risk of adjacent vertebral fracture; which may be due to the change in material properties of the post-procedure vertebra modifying the transmission of mechanical stresses to adjacent vertebrae.

The neuro-muscular and musculo-skeletal characterization of children with joint hypermobility

In children, joint hypermobility (typified by structural instability of joints) manifests clinically as neuro-muscular and musculo-skeletal conditions and conditions associated with development and organization of control of posture and gait (Finkelstein, 1916; Jahss, 1919; Sobel, 1926; Larsson, Mudholkar, Baum and Srivastava, 1995; Murray and Woo, 2001; Hakim and Grahame, 2003; Adib, Davies, Grahame, Woo and Murray, 2005:). The process of control of the relative proportions of joint mobility and stability, whilst maintaining equilibrium in standing posture and gait, is dependent upon the complex interrelationship between skeletal, muscular and neurological function (Massion, 1998; Gurfinkel, Ivanenko, Levik and Babakova, 1995; Shumway-Cook and Woollacott, 1995). The efficiency of this relies upon the integrity of neuro-muscular and musculo-skeletal components (ligaments, muscles, nerves), and the Central Nervous System’s capacity to interpret, process and integrate sensory information from visual, vestibular and proprioceptive sources (Crotts, Thompson, Nahom, Ryan and Newton, 1996; Riemann, Guskiewicz and Shields, 1999; Schmitz and Arnold, 1998) and development and incorporation of this into a representational scheme (postural reference frame) of body orientation with respect to internal and external environments (Gurfinkel et al., 1995; Roll and Roll, 1988). Sensory information from the base of support (feet) makes significant contribution to the development of reference frameworks (Kavounoudias, Roll and Roll, 1998). Problems with the structure and/ or function of any one, or combination of these components or systems, may result in partial loss of equilibrium and, therefore ineffectiveness or significant reduction in the capacity to interact with the environment, which may result in disability and/ or injury (Crotts et al., 1996; Rozzi, Lephart, Sterner and Kuligowski, 1999b). Whilst literature focusing upon clinical associations between joint hypermobility and conditions requiring therapeutic intervention has been abundant (Crego and Ford, 1952; Powell and Cantab, 1983; Dockery, in Jay, 1999; Grahame, 1971; Childs, 1986; Barton, Bird, Lindsay, Newton and Wright, 1995a; Rozzi, et al., 1999b; Kerr, Macmillan, Uttley and Luqmani, 2000; Grahame, 2001), there has been a deficit in controlled studies in which the neuro-muscular and musculo-skeletal characteristics of children with joint hypermobility have been quantified and considered within the context of organization of postural control in standing balance and gait. This was the aim of this project, undertaken as three studies. The major study (Study One) compared the fundamental neuro-muscular and musculo-skeletal characteristics of 15 children with joint hypermobility, and 15 age (8 and 9 years), gender, height and weight matched non-hypermobile controls. Significant differences were identified between previously undiagnosed hypermobile (n=15) and non-hypermobile children (n=15) in passive joint ranges of motion of the lower limbs and lumbar spine, muscle tone of the lower leg and foot, barefoot CoP displacement and in parameters of barefoot gait. Clinically relevant differences were also noted in barefoot single leg balance time. There were no differences between groups in isometric muscle strength in ankle dorsiflexion, knee flexion or extension. The second comparative study investigated foot morphology in non-weight bearing and weight bearing load conditions of the same children with and without joint hypermobility using three dimensional images (plaster casts) of their feet. The preliminary phase of this study evaluated the casting technique against direct measures of foot length, forefoot width, RCSP and forefoot to rearfoot angle. Results indicated accurate representation of elementary foot morphology within the plaster images. The comparative study examined the between and within group differences in measures of foot length and width, and in measures above the support surface (heel inclination angle, forefoot to rearfoot angle, normalized arch height, height of the widest point of the heel) in the two load conditions. Results of measures from plaster images identified that hypermobile children have different barefoot weight bearing foot morphology above the support surface than non-hypermobile children, despite no differences in measures of foot length or width. Based upon the differences in components of control of posture and gait in the hypermobile group, identified in Study One and Study Two, the final study (Study Three), using the same subjects, tested the immediate effect of specifically designed custom-made foot orthoses upon balance and gait of hypermobile children. The design of the orthoses was evaluated against the direct measures and the measures from plaster images of the feet. This ascertained the differences in morphology of the modified casts used to mould the orthoses and the original image of the foot. The orthoses were fitted into standardized running shoes. The effect of the shoe alone was tested upon the non-hypermobile children as the non-therapeutic equivalent condition. Immediate improvement in balance was noted in single leg stance and CoP displacement in the hypermobile group together with significant immediate improvement in the percentage of gait phases and in the percentage of the gait cycle at which maximum plantar flexion of the ankle occurred in gait. The neuro-muscular and musculo-skeletal characteristics of children with joint hypermobility are different from those of non-hypermobile children. The Beighton, Solomon and Soskolne (1973) screening criteria successfully classified joint hypermobility in children. As a result of this study joint hypermobility has been identified as a variable which must be controlled in studies of foot morphology and function in children. The outcomes of this study provide a basis upon which to further explore the association between joint hypermobility and neuro-muscular and musculo-skeletal conditions, and, have relevance for the physical education of children with joint hypermobility, for footwear and orthotic design processes, and, in particular, for clinical identification and treatment of children with joint hypermobility.