The flexion-extension profile of lumbar spine in 100 healthy volunteers

Spinal loads during dynamic full flexion and return to standing posture in different age and sex groups: A musculoskeletal model study

During forward flexion, spine motion varies due to age and sex differences. Previous studies showed that lumbar/pelvis range of flexion (RoF) and lumbo-pelvic ratio (L/P) are age/sex dependent. How variation of these parameters affects lumbar loading in a normal population requires further assessment. We aimed to estimate lumbar loads during dynamic flexion-return cycle and the differences in peak loads (compression) and corresponding trunk inclinations due to variation in lumbar/pelvis RoF and L/P. Based on in vivo L/P (0.11-3.44), temporal phases of flexion (early, middle, and later), the lumbar (45-55°) and hip (60-79°) RoF; full flexion-return cycles of six seconds were reconstructed for three age groups (20-35, 36-50 and 50+ yrs.) in both sexes. Six inverse dynamic analyses were performed with a 50th percentile model, and differences in peak loads and corresponding trunk inclinations were calculated. Peak loads at L4-L5 were 179 N higher in younger males versus females, but 228 N and 210 N lower in middle-aged and older males, respectively, compared to females. Females exhibited higher trunk inclinations (6°-20°) than males across all age groups. Age related differences in L4-L5 peak loads and corresponding trunk inclinations were found up to 415 N and 19° in males and 152 N and 13° in females. With aging, peak loads were reduced in males but were found non-monotonic in females, whereas trunk inclinations at peak loads were reduced in both sexes from young to middle/old age groups. In conclusion, lumbar loading and corresponding trunk inclinations varied notably due to age/sex differences. Such data may help distinguishing normal or pathological condition of the lumbar spine.

Development and evaluation of sex-specific thoracolumbar spine finite element models to study spine biomechanics

Musculoskeletal disorders and low back pain (LBP) are common global afflictions, with a higher prevalence observed in females. However, the cause of many LBP cases continues to elude researchers. Current approaches seldom consider differences in male and female spines. Thus, this study aimed to compare the load distribution between male and female spines through finite element modeling. Two finite element models of the spine, one male and one female, were developed, inclusive of sex-specific geometry and material properties. The models consisted of the vertebrae, intervertebral discs (IVD), tendons, surrounding spinal muscles, and thoracolumbar fascia and were subjected to loading conditions simulating flexion and extension. Following extensive validation against published literature, intersegmental rotation, IVD stress, and vertebral body stress were evaluated. The female model demonstrated increased magnitudes for rotation and stresses when compared to the male model. Results suggest that the augmented stresses in the female model indicate an increased load distribution throughout the spine compared to the male model. These findings may corroborate the higher prevalence of LBP in females. This study highlights the importance of using patient- and sex-specific models for patient analyses and care.

Numerical investigation of intra-abdominal pressure and spinal load-sharing upon the application of an abdominal belt

Chronic low back pain patients may experience spinal instability. Abdominal belts (ABs) have been shown to improve spine stability, trunk stiffness, and resiliency to spinal perturbations. However, research on the contributing mechanisms is inconclusive. ABs may increase intra-abdominal pressure (IAP) and reduce paraspinal soft tissue contribution to spine stability without increasing spinal compressive loads. A finite element model (FEM) of the spine inclusive of the T1-S1 vertebrae, intervertebral discs (IVDs), ribcage, pelvis, soft tissues, and abdominal cavity, without active muscle forces was developed. An identical FEM with an AB was developed. Both FEMs underwent trunk flexion. Following validation, the models' intervertebral rotation (IVR), IAP, IVD pressure, and tensile stress in the multifidus (MF), erector spinae (ES), and thoracolumbar fascia (TLF) were compared. The inclusion of an AB resulted in a 3.8 kPa IAP increase, but a decreased average soft tissue tensile stress of 0.28 kPa. The TLF withstood the majority of tension being transferred across the paraspinal soft tissues (>70 %). The average IVR in the AB model decreased by 10 %, with the lumbar spine experiencing the largest reduction. The lumbar IVDs of the AB model likewise showed a 31 % reduction in average IVD pressure. Using an AB improved trunk bending stiffness, primarily in the lumbar spine. Wearing an AB had minimal effect on reducing tensile stress in theES. The skewed stress distribution towards the TLF suggests its large contribution to spine stability and the potential advantage in unloading the structure when wearing an AB, measured herein at8 %.

Influence of the thoracolumbar junction flexibility on the risk of adding-on after posterior vertebral arthrodesis for thoracic idiopathic adolescent scoliosis

PURPOSE

The objective was to analyze the role of the thoracolumbar sagittal flexibility on the outcome after posterior spinal fusion of Lenke 1 and 2 adolescent idiopathic scoliosis with last touched vertebra as the lowest instrumented vertebra.

METHODS

We included 105 thoracic AIS patients who had a posterior spinal fusion with a 2 years minimum follow-up. Thoracolumbar junction flexibility was assessed on dynamic sagittal X-rays and compared to the standing position. Adding-on was defined according to radiographic Wang criteria. The junction was considered flexible if the variability from the static position to flexion and/or extension was greater than 10°.

RESULTS

Mean age of the patients was 14 ± 2 years. The preoperative mean Cobb angle was 61 ± 12.7° and 27.5 ± 7.7° after surgery. Mean follow-up was 3.1 years. Twenty-nine patients (28%) developed an adding-on. Thoracolumbar junction range of motion was higher (p = 0.017) with higher flexibility in flexion (p < 0.001) in the no adding-on group. In no adding-on group, 53 patients (70%) had a flexible thoracolumbar junction, and 23 patients (30%) had a stiff thoracolumbar junction in flexion and flexible in extension. In adding-on group, 27 patients (93%) had a stiff thoracolumbar junction, and 2 patients (7%) had a flexible junction in flexion and stiff in extension.

CONCLUSION

The flexibility of the thoracolumbar junction is a determining factor in the surgical outcome after posterior spinal fusion for AIS and should be considered in correlation with the frontal and sagittal alignment of the spine.

The Variability of Lumbar Sequential Motion Patterns: Observational Study

BACKGROUND

Physiological motion of the lumbar spine is a topic of interest for musculoskeletal health care professionals since abnormal motion is believed to be related to lumbar complaints. Many researchers have described ranges of motion for the lumbar spine, but only few have mentioned specific motion patterns of each individual segment during flexion and extension, mostly comprising the sequence of segmental initiation in sagittal rotation. However, an adequate definition of physiological motion is still lacking. For the lower cervical spine, a consistent pattern of segmental contributions in a flexion-extension movement in young healthy individuals was described, resulting in a definition of physiological motion of the cervical spine.

OBJECTIVE

This study aimed to define the lumbar spines' physiological motion pattern by determining the sequence of segmental contribution in sagittal rotation of each vertebra during maximum flexion and extension in healthy male participants.

METHODS

Cinematographic recordings were performed twice in 11 healthy male participants, aged 18-25 years, without a history of spine problems, with a 2-week interval (time point T1 and T2). Image recognition software was used to identify specific patterns in the sequence of segmental contributions per individual by plotting segmental rotation of each individual segment against the cumulative rotation of segments L1 to S1. Intraindividual variability was determined by testing T1 against T2. Intraclass correlation coefficients were tested by reevaluation of 30 intervertebral sequences by a second researcher.

RESULTS

No consistent pattern was found when studying the graphs of the cinematographic recordings during flexion. A much more consistent pattern was found during extension, especially in the last phase. It consisted of a peak in rotation in L3L4, followed by a peak in L2L3, and finally, in L1L2. This pattern was present in 71% (15/21) of all recordings; 64% (7/11) of the participants had a consistent pattern at both time points. Sequence of segmental contribution was less consistent in the lumbar spine than the cervical spine, possibly caused by differences in facet orientation, intervertebral discs, overprojection of the pelvis, and muscle recruitment.

CONCLUSIONS

In 64% (7/11) of the recordings, a consistent motion pattern was found in the upper lumbar spine during the last phase of extension in asymptomatic young male participants. Physiological motion of the lumbar spine is a broad concept, influenced by multiple factors, which cannot be captured in a firm definition yet.

TRIAL REGISTRATION

ClinicalTrials.gov NCT03737227; https://clinicaltrials.gov/ct2/show/NCT03737227.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

RR2-10.2196/14741.

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

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

Biomechanical evaluation of the thoracolumbar spine comparing healthy and irregular thoracic and lumbar curvatures

BACKGROUND

The geometric alignment of the spine plays an integral role in stability, biomechanical loading, and consequently, pain, and a range of healthy sagittal curvatures has been identified. Spinal biomechanics when sagittal curvature is outside the optimal range remains a debate and may provide insight into the load distribution throughout the spinal column.

METHOD

A thoracolumbar spine model (Healthy) was developed. Thoracic and lumbar curvatures were adjusted by 50% to create models with varying sagittal profiles: hypolordotic (HypoL), hyperlordotic (HyperL), hypokyphotic (HypoK), and hyperkyphotic (HyperK). In addition, lumbar spine models were constructed for the former three profiles. The models were subjected to loading conditions simulating flexion and extension. Following validation, intervertebral disc stresses, vertebral body stresses, disc heights, and intersegmental rotations were compared across all models.

RESULTS

Overall trends demonstrated that HyperL and HyperK models had a noticeable reduction in disc height and greater vertebral body stresses compared to the Healthy model. In comparison, the HypoL and HypoK models displayed opposite trends. Considering the lumbar models, the HypoL model had reduced disc stresses and flexibility, while the contrary was observed in the HyperL model. Results indicate that the models with excessive curvature may be subjected to greater stress magnitudes, while the straighter spine models may reduce these stresses.

CONCLUSIONS

Finite element modeling of spine biomechanics demonstrated that variations in sagittal profiles influence the load distribution and range of motion of the spine. Considering patient-specific sagittal profiles in finite element modeling may provide valuable insight for biomechanical analyses and targeted treatments.

Lumbosacral transitional vertebrae alter the distribution of lumbar mobility–Preliminary results of a radiographic evaluation

Background

Lumbo-sacral transitional vertebrae (LSTV) are one of the most common congenital variances of the spine. They are associated with an increased frequency of degeneration in the cranial adjacent segment. Hypermobility and concomitant increased loads are discussed as a possible reason for segmental degeneration. We therefore examined the lumbar and segmental motion distribution in patients with LSTV with flexion-extension radiographs.

Methods

A retrospective study of 51 patients with osteochondrosis L5/S1 with flexion and extension radiographs was performed. Of these, 17 patients had LSTV and were matched 1:1 for age and sex with patients without LSTV out of the collective of the remaining 34 patients. The lumbar and segmental range of motion (RoM) by segmental lordosis angle and the segmental wedge angle were determined. Normal distribution of parameters was observed by Kolmogorov-Smirnov-test. Parametric data were compared by paired T-test. Non-parametric data were compared by Wilcoxon-rank-sum-test. Correlations were observed using Spearman’s Rank correlation coefficient. A p-value <0.05 was stated as statistically significant.

Results

Patients with LSTV had mean age of 52.2±10.9, control group of 48.9±10.3. Both groups included 7 females and 10 males. Patients with LSTV presented with reduced RoM of the lumbar spine (LSTV 37.3°±19.2°, control 52.1°±20.5°, p = 0.065), however effects were statistically insignificant. LSTV significantly decreased segmental RoM in the transitional segment (LSTV 1.8°±2.7°, control 6.7°±6.0°, p = 0.003). Lumbar motion distribution differed significantly; while RoM was decreased in the transitional segment, (LSTV 5.7%, control 16.2%, p = 0.002), the distribution of lumbar motion to the cranial adjacent segment was increased (LSTV 30.7%, control 21.6%, p = 0.007).

Conclusion

Patients with LSTV show a reduced RoM in the transitional segment and a significantly increased motion distribution to the cranial adjacent segment in flexion-extension radiographs. The increased proportion of mobility in the cranial adjacent segment possibly explain the higher rates of degeneration within the segment.

Biomechanical Properties of Paraspinal Muscles Influence Spinal Loading—A Musculoskeletal Simulation Study

Paraspinal muscles are vital to the functioning of the spine. Changes in muscle physiological cross-sectional area significantly affect spinal loading, but the importance of other muscle biomechanical properties remains unclear. This study explored the changes in spinal loading due to variation in five muscle biomechanical properties: passive stiffness, slack sarcomere length (SSL), in situ sarcomere length, specific tension, and pennation angle. An enhanced version of a musculoskeletal simulation model of the thoracolumbar spine with 210 muscle fascicles was used for this study and its predictions were validated for several tasks and multiple postures. Ranges of physiologically realistic values were selected for all five muscle parameters and their influence on L4-L5 intradiscal pressure (IDP) was investigated in standing and 36° flexion. We observed large changes in IDP due to changes in passive stiffness, SSL, in situ sarcomere length, and specific tension, often with interesting interplays between the parameters. For example, for upright standing, a change in stiffness value from one tenth to 10 times the baseline value increased the IDP only by 91% for the baseline model but by 945% when SSL was 0.4 μm shorter. Shorter SSL values and higher stiffnesses led to the largest increases in IDP. More changes were evident in flexion, as sarcomere lengths were longer in that posture and thus the passive curve is more influential. Our results highlight the importance of the muscle force-length curve and the parameters associated with it and motivate further experimental studies on in vivo measurement of those properties.

How Can Biomechanical Multibody Models of Scoliosis Be Accurate in Simulating Spine Movement Behavior While Neglecting the Changes of Spinal Length?

This paper studies how biomechanical multibody models of scoliosis can neglect the changes of spinal length and yet be accurate in reconstructing spinal columns. As these models with fixed length comprise rigid links interconnected by rotary joints, they resemble polygonal chains that approximate spine curves with a finite number of line segments. In mathematics, using more segments with shorter lengths can result in more accurate curve approximations. This raises the question of whether more accurate spine curve approximations by increasing the number of links/joints can yield more accurate spinal column reconstructions. For this, the accuracy of spine curve approximation was improved consistently by increasing the number of links/joints, and its effects on the accuracy of spinal column reconstruction were assessed. Positive correlation was found between the accuracy of spine reconstruction and curve approximation. It was shown that while increasing the accuracy of curve approximations, the representation of scoliosis concavity and its side-to-side deviations were improved. Moreover, reconstruction errors of the spine regions separated by the inflection vertebrae had minimal impacts on each other. Overall, multibody scoliosis models with fixed spinal lengths can benefit from the extra rotational joints that contribute toward the accuracy of spine curve approximation. The outcome of this study leads to concurrent accuracy improvement and simplification of multibody models; joint-link configurations can be independently defined for the regions separated by the inflection vertebrae, enabling local optimization of the models for higher accuracy without unnecessary added complexity to the whole model.

The examination of stress shielding in a finite element lumbar spine inclusive of the thoracolumbar fascia

Despite the prevalence of low back pain (LBP) in society, the pathomechanism of LBP continues to elude researchers. LBP patients have demonstrated morphological and material property changes to their lumbar soft tissues, potentially leading to irregular load sharing within the lumbar spine. This study aims to analyze potential stress shielding consequential of augmented soft tissue properties via the comparison of a healthy and LBP finite element models. The models developed in this study include the vertebrae, intervertebral discs and soft tissues from L1-S1. Soft tissue morphology and material properties for the LBP model were augmented to reflect documented clinical findings. Model validation preceded testing and was confirmed through comparison to the available literature. Relative to the healthy model, the LBP model demonstrated an increase in stress by 15.6%, with 99.8% of this stress increase being distributed towards the thoracolumbar fascia. The majority of stress skewed towards the fascia may indicate a potential stress allocation bias whereby the lumbar muscles are unable to receive regular loading, leading to stress shielding. This load allocation bias and subsequent stress shielding may potentially contribute to the progression and pathomechanism of LBP but prospective studies would be required to make that link.

Investigation of physiological stress shielding within lumbar spinal tissue as a contributor to unilateral low back pain: A finite element study

INTRODUCTION

The pathomechanism of low back pain (LBP) remains unknown. Unilateral LBP patients have demonstrated ipsilateral morphological and material property changes within the lumbar soft tissues, potentially leading to asymmetric tissue loading. Through the comparison of healthy and unilateral LBP validated finite element models (FEMs), this study investigates potential stress shielding consequential of spinal tissue property augmentation.

METHODS

Two FEMs of the musculoskeletal system - one demonstrating healthy and unilateral LBP conditions - were developed undergoing 30-degree flexion. FEMs included the vertebrae, intervertebral discs, and soft tissues from L1-S1. Material properties selected for the soft tissues were retrieved from published literature. To reflect unilateral LBP, the paraspinal morphology was atrophied, while the tissue moduli were increased. The symptomatic thoracolumbar fascia (TLF) was uniformly increased. Validation of the models preceded testing.

RESULTS

Model validation in spinal flexion was accomplished through comparison to literature. Compared to the healthy model, the unilateral LBP multifidus (MF), longissimus thoracis (LT), and TLF exhibited average tension changes of +7.9, -5.1, and +9.3%, respectively. Likewise, the symptomatic MF, LT, and TLF exhibited tension changes of +19.0, -10.4, and +16.1% respectively, whereas the asymptomatic MF, LT, and TLF exhibited -4.0, -2.0, and +0.4% changes in tension, respectively.

CONCLUSION

Relative to the healthy tissues, the symptomatic LBP soft tissues demonstrated a 19.5 kPa increase in stress, with 99.8% of this increase distributed towards the TLF, suggesting a load allocation bias within the symptomatic unilateral LBP tissues. Consequentially, symptomatic paraspinal muscles may be unable to withstand loading, leading to stress shielding.

ISSLS Prize in Bioengineering Science 2021: in vivo sagittal motion of the lumbar spine in low back pain patients-a radiological big data study

PURPOSE

We investigated the flexion-extension range of motion and centre of rotation of lumbar motion segments in a large population of 602 patients (3612 levels), and the associations between lumbar motion and other parameters such as sex, age and intervertebral disc degeneration.

METHODS

Lumbar radiographs in flexion-extension of 602 patients suffering from low back pain and/or suspect instability were collected; magnetic resonance images were retrieved and used to score the degree of disc degeneration for a subgroup of 354 patients. Range of motion and centre of rotation were calculated for all lumbosacral levels with in-house software allowing for high degree of automation. Associations between motion parameters and age, sex, spinal level and disc degeneration were then assessed.

RESULTS

The median range of motion was 6.6° (range 0.1-28.9°). Associations between range of motion and age as well as spinal level, but not sex, were found. Disc degeneration determined a consistent reduction in the range of motion. The centre of rotation was most commonly located at the centre of the lower endplate or slightly lower. With progressive degeneration, centres of rotation were increasingly dispersed with no preferential directions.

CONCLUSION

This study constitutes the largest analysis of the in vivo lumbar motion currently available and covers a wide range of clinical scenarios in terms of age and degeneration. Findings confirmed that ageing determines a reduction in the mobility independently of degeneration and that in degenerative levels, centres of rotation are dispersed around the centre of the intervertebral space.

Characterization of biomaterials intended for use in the nucleus pulposus of degenerated intervertebral discs

Biomaterials for regeneration of the intervertebral disc must meet complex requirements conforming to biological, mechanical and clinical demands. Currently no consensus on their characterization exists. It is crucial to identify parameters and their method of characterization for accurate assessment of their potential efficacy, keeping in mind the translation towards clinical application. This review systematically analyses the characterization techniques of biomaterial systems that have been used for nucleus pulposus (NP) restoration and regeneration. Substantial differences in the approach towards assessment became evident, hindering comparisons between different materials with respect to their suitability for NP restoration and regeneration. We have analysed the current approaches and identified parameters necessary for adequate biomaterial characterization, with the clinical goal of functional restoration and biological regeneration of the NP in mind. Further, we provide guidelines and goals for their measurement. STATEMENT OF SIGNIFICANCE: Biomaterials intended for restoration of regeneration of the nucleus pulposus within the intervertebral disc must meet biological, biomechanical and clinical demands. Many materials have been investigated, but a lack of consensus on which parameters to evaluate leads to difficulties in comparing materials as well as mostly partial characterization of the materials in question. A gap between current methodology and clinically relevant and meaningful characterization is prevalent. In this article, we identify necessary methods and their implementation for complete biomaterial characterization in the context of clinical applicability. This will allow for a more unified approach to NP-biomaterials research within the field as a whole and enable comparative analysis of novel materials yet to be developed.

A novel method for prediction of postoperative global sagittal alignment based on full-body musculoskeletal modeling and posture optimization

Associations between spinal sagittal balance and pain and disability are well documented. Reciprocal changes after spinal surgery might be critical for the outcomes, but assessing their extent remains a challenge. This paper proposes a method for predicting full-body sagittal alignment including reciprocal changes in response to spinal fusion, based on musculoskeletal modeling and inverse-inverse dynamics approach. An established body model (AnyBody) was used, with fused segments modeled as rigid. Posture was optimized based on muscle expenditure minimization, following the concept of the cone of economy. The data of adult spinal fusion patients were obtained retrospectively from an ongoing clinical study. Patient spino-pelvic alignment, body weight and height, age- and pathology-related muscle deterioration, and underwent treatment details were represented in the model. Predicted postural changes were compared to follow-up radiographs to evaluate method validity. Twenty-one cases were analyzed in this preliminary study (age range = 48-74; number of fused segments 1-14). The model predictions correlated well with the radiographic measures at follow-up: TPA, r = 0.83; ΔPILL, r = 0.90; LL, r = 0.90; TK, r = 0.77. The model demonstrated high accuracy in predicting sagittal imbalance (positive predictive value = 1.00, negative predictive value = 0.75). The presented method for patient- and treatment-specific postoperative posture prediction can be used to guide preoperative planning of spinal fusion, but more extensive validation is needed.

Lumbar Intervertebral Motion Analysis During Flexion and Extension Cinematographic Recordings in Healthy Male Participants: Protocol

Background

Physiological motion of the lumbar spine is a subject of interest for musculoskeletal health care professionals, as abnormal motion is believed to be related to lumbar conditions and complaints. Many researchers have described ranges of motion for the lumbar spine, but only a few have mentioned specific motion patterns of each individual segment during flexion and extension. These motion patterns mostly comprise the sequence of segmental initiation in sagittal rotation. However, an adequate definition of physiological motion of the lumbar spine is still lacking. The reason for this is the reporting of different ranges of motion and sequences of segmental initiation in previous studies. Furthermore, due to insufficient fields of view, none of these papers have reported on maximum flexion and extension motion patterns of L1 to S1. In the lower cervical spine, a consistent pattern of segmental contributions was recently described. In order to understand physiological motion of the lumbar spine, it is necessary to systematically study motion patterns, including the sequence of segmental contribution, of vertebrae L1 to S1 in healthy individuals during maximum flexion and extension.

Objective

This study aims to define the lumbar spines’ physiological motion pattern of vertebrae L1, L2, L3, L4, L5, and S1 by determining the sequence of segmental contribution and the sequence of segmental initiation of motion in sagittal rotation of each vertebra during maximum flexion and extension. The secondary endpoint will be exploring the possibility of analyzing the intervertebral horizontal and vertical translation of each vertebra during maximum flexion and extension.

Methods

Cinematographic recordings will be performed on 11 healthy male participants, aged 18-25 years, without a history of spine problems. Cinematographic flexion and extension recordings will be made at two time points with a minimum 2-week interval in between.

Results

The study has been approved by the local institutional medical ethical committee (Medical Research Ethics Committee of Zuyderland and Zuyd University of Applied Sciences) on September 24, 2018. Inclusion of participants will be completed in 2020.

Conclusions

If successful, these physiological motion patterns can be compared with motion patterns of patients with lumbar conditions before or after surgery. Ultimately, researchers may be able to determine differences in biomechanics that can potentially be linked to physical complaints like low back pain.

Trial Registration

ClinicalTrials.gov NCT03737227; https://clinicaltrials.gov/ct2/show/NCT03737227

International Registered Report Identifier (IRRID)

DERR1-10.2196/14741

Discriminating spatiotemporal movement strategies during spine flexion-extension in healthy individuals

BACKGROUND CONTEXT

The spine is an anatomically complex system with numerous degrees of freedom. Due to this anatomical complexity, it is likely that multiple motor control options exist to complete a given task.

PURPOSE

To identify if distinct spine spatiotemporal movement strategies are utilized in a homogenous sample of young healthy participants.

STUDY DESIGN

Kinematic data were captured from a single cohort of male participants (N=51) during a simple, self-controlled spine flexion-extension task.

METHODS

Thoracic and lumbar flexion-extension data were analyzed to extract the continuous relative phase between each spine subsection. Continuous relative phase data were evaluated using a principal component analysis to identify major sources of variation in spine movement coordination. Unsupervised machine learning (k-means clustering) was used to identify distinct clusters present within the healthy participants sampled. Once distinguished, intersegmental spine kinematics were compared amongst clusters.

RESULTS

The findings of the current work suggest that there are distinct timing strategies that are utilized, within the participants sampled, to control spine flexion-extension movement. These strategies differentiate the sequencing of intersegmental movement and are not discriminable on the basis of simple participant demographic characteristics (ie, age, height, and body mass index), total movement time or range of motion.

CONCLUSIONS

Spatiotemporal spine flexion-extension patterns are not uniform across a population of young healthy individuals.

CLINICAL SIGNIFICANCE

Future work needs to identify whether the motor patterns characterized with this work are driven by distinct neuromuscular activation patterns, and if each given pattern has a varied risk for low back injury.

Five-year follow-up of clinical and radiological outcomes of LP-ESP elastomeric lumbar total disc replacement in active patients

BACKGROUND CONTEXT

The surgical treatment of degenerative disc disease at the lumbar spine may involve fusion. Total disc replacement (TDR) is an alternative treatment to avoid fusion-related adverse events, specifically adjacent segment disease. New generation of elastomeric non-articulating devices has been developed to more effectively replicate the shock absorption and flexural stiffness of native disc.

PURPOSE

To report 5 years clinical and radiographic outcomes, range of motion (ROM), and position of the center of rotation after a viscoelastic lumbar TDR.

STUDY DESIGN

Prospective observational cohort study PATIENT SAMPLE: Sixty-one patients OUTCOME MEASURES: The clinical evaluation was based on visual analog scale (VAS) for pain, Oswestry disability index (ODI) score, short form-36 (SF-36) including the physical component summary (PCS) and the mental component summary (MCS), and general health questionnaire-28 (GHQ28). The radiological outcomes were ROM and position of the center of rotation at the index and the adjacent levels and the adjacent disc height changes.

METHODS

Our study group included 61 consecutive patients with monosegmental disc replacement. We selected patients who could provide a global lumbar spine mobility analysis (intermediate functional activity according to the Baecke score). Hybrid constructs had been excluded. Only the cases with complete clinical and radiological follow-up at 3, 6, 12, 24, and 60 months were included.

RESULTS

There was a significant improvement in VAS (3.3±2.5 vs. 6.6±1.7, p<.001), in ODI (20±17.9 vs. 51.2±14.6, p<.001), GHQ28 (52.6±15.5 vs. 64.2±15.6, p<.001), SF-36 PCS (58.8±4.8 vs. 32.4±3.4, p<.001), and SF-36 MCS (60.7±6 vs. 42.3±3.4, p<.001). The mean location centers of the index level and adjacent discs were comparable to those previously published in asymptomatic patients. According to the definition of Zigler and Delamarter, all of our cases remained grade 0 for adjacent level disc height (within 25% of normal).

CONCLUSIONS

This series reports significant improvement in midterm follow-up after TDR, which is consistent with previously published studies but with a lower rate of revision surgery and no adjacent level disease pathologies. The radiographic assessment of the patients demonstrated the quality of functional reconstruction of the lumbar spine after LP-ESP viscoelastic disc replacement.

Effect of age and sex on lumbar lordosis and the range of motion. A systematic review and meta-analysis

Lumbar lordosis (LL) and the range of motion (RoM) are important physiological measurements when initiating any diagnosis and treatment plan for patients with low back pain. Numerous studies reported differences in LL and the RoM due to age and sex. However, these findings remain contradictory. A systematic review and meta-analysis were performed to synthesize mean values and the differences in LL and the RoM because of age and sex. The quality assessment tool for quantitative studies was applied to assess the methodological quality of the studies included. We identified 2372 papers through electronic (2309) and physical (63) searches. We assessed 218 full-text studies reporting measurements of LL or the RoM. In total, 65 studies were included, and a normative database for LL and the RoM is provided as supplementary material. Among these, 11 were included in the meta-analysis. LL and the RoM displayed non-monotonic variations with significant age and sex differences. Young females showed a significantly greater LL and the range of extension (RoE), whereas young males exhibited a greater range of flexion (RoF). Sex differences in the range of lateral bending (RoLB) were small but were significant for the axial rotation (RoAR). For the RoF, RoE and RoLB, differences because of age were significant among most of the age groups in both sexes, whereas for the RoAR, differences were significant only between the 20s vs the 30s-40s (males) and 40s vs 50s (females). Significant differences because of age/sex were identified. However, the age-dependent reduction in LL and the RoM was non-monotonic and differed in both sexes. These findings will help to better distinguish between functional deficits caused by spinal disorders and natural factors/conditions related to age and sex.

MRI evaluation of the effects of extension exercises on the disc fluid content and location of the centroid of the fluid distribution

BACKGROUND

McKenzie prone press-up exercises have been hypothesised to reduce intradiscal pressure, allowing fluid to be reabsorbed into the disc, which could improve the internal stability and local chemical milieu of the disc, potentially reducing symptoms.

OBJECTIVE

To investigate the immediate effects of prone press-up exercises on lumbar disc fluid content and movement.

DESIGN

Quantification of MRI changes before and after a single exercise session.

METHODS

The mid-sagittal T₂-weighted MR images of 22 volunteers with low back pain were obtained before and immediately after performing press-up exercises. The whole disc and nucleus regions of the L4-5 and L5-S1 discs were then segmented, and their mean signal intensity (MSI) and signal intensity weighted centroid (SIWC) were computed to estimate disc fluid content and displacement.

RESULTS

There were no significant differences between the MSI and the vertical position of the SIWC of the whole disc before and after extension at either disc level (effect size [ES]: -0.23 to 0.09). There was a significant anterior displacement (0.1 ± 5.4 mm) of the location of the SIWC of the disc after extension exercise at L4-5 (ES: 0.22), but not at L5-S1 (ES: 0.00) or at either level for the nucleus region (ES: -0.06; 0.16).

CONCLUSION

Little evidence was found supporting the hypothesis that press-up exercises affect disc fluid content and distribution. Novel parameters reflecting fluid distribution detected similar or larger effects of the extension than MSI. If such exercises are effective in reducing symptoms, it is likely through other mechanisms than by changing fluid content or distribution.

Thoracolumbar spine loading associated with kinematics of the young and the elderly during activities of daily living

Excessive mechanical loading of the spine is a critical factor in vertebral fracture initiation. Most vertebral fractures develop spontaneously or due to mild trauma, as physiological loads during activities of daily living might exceed the failure load of osteoporotic vertebra. Spinal loading patterns are affected by vertebral kinematics, which differ between elderly and young individuals. In this study, the effects of age-related changes in spine kinematics on thoracolumbar spinal segmental loading during dynamic activities of daily living were investigated using combined experimental and modeling approach. Forty-four healthy volunteers were recruited into two age groups: young (N = 23, age = 27.1 ± 3.8) and elderly (N = 21, age = 70.1 ± 3.9). The spinal curvature was assessed with a skin-surface device and the kinematics of the spine and lower extremities were recorded during daily living tasks (flexion-extension and stand-sit-stand) with a motion capture system. The obtained data were used as input for a musculoskeletal model with a detailed thoracolumbar spine representation. To isolate the effect of kinematics on predicted loads, other model properties were kept constant. Inverse dynamics simulations were performed in the AnyBody Modeling System to estimate corresponding spinal loads. The maximum compressive loads predicted for the elderly motion patterns were lower than those of the young for L2/L3 and L3/L4 lumbar levels during flexion and for upper thoracic levels during stand-to-sit (T1/T2-T8/T9) and sit-to-stand (T3/T4-T6/T7). However, the maximum loads predicted for the lower thoracic levels (T9/T10-L1/L2), a common site of vertebral fractures, were similar compared to the young. Nevertheless, these loads acting on the vertebrae of reduced bone quality might contribute to a higher fracture risk for the elderly.

Is L5-S1 motion segment different from the rest? A radiographic kinematic assessment of 72 patients with chronic low back pain

PURPOSE

The relationship between biomechanical instability and degenerative changes in the lumbar spine in chronic low back pain (CLBP) patients remains controversial. The main objective of this retrospective radiographical study was to evaluate changes in kinematics at different lumbar levels (in particular the L5-S1 level) with progressive grades of disc degeneration and facet joint osteoarthritis in CLBP patients.

METHODS

Using standing neutral and dynamic flexion/extension (Fx/Ex) radiographs of the lumbar spine, in vivo segmental kinematics at L1-L2 through L5-S1 were evaluated in 72 consecutive CLBP patients. Disc degeneration was quantified using changes in signal intensity and central disc height on mid-sagittal T2-weighted magnetic resonance (MR) scans. Additionally, the presence or absence of facet joint osteoarthritis was noted on T2-weighted axial MR scans.

RESULTS

Disc degeneration and facet joint osteoarthritis occurred independent of each other at the L5-S1 level (p = 0.188), but an association was observed between the two at L4-L5 (p < 0.001) and L3-L4 (p < 0.05) levels. In the absence of facet joint osteoarthritis, the L5-S1 segment showed a greater range of motion (ROM) in Ex (3.3° ± 3.6°) and a smaller ROM in Fx (0.6° ± 4.2°) compared with the upper lumbar levels (p < 0.05), but the differences diminished in the presence of it. In the absence of facet joint osteoarthritis, no change in L5-S1 kinematics was observed with progressive disc degeneration, but in its presence, restabilisation of the L5-S1 segment was observed between mild and severe disc degeneration states.

CONCLUSION

The L5-S1 motion segment exhibited unique degenerative and kinematic characteristics compared with the upper lumbar motion segments. Disc degeneration and facet joint osteoarthritis occurred independent of each other at the L5-S1 level, but not at the other lumbar levels. Severe disc degeneration in the presence of facet joint osteoarthritis biomechanically restabilised the L5-S1 motion segment.

Uneven intervertebral motion sharing is related to disc degeneration and is greater in patients with chronic, non-specific low back pain: an in vivo, cross-sectional cohort comparison of intervertebral dynamics using quantitative fluoroscopy

PURPOSE

Evidence of intervertebral mechanical markers in chronic, non-specific low back pain (CNSLBP) is lacking. This research used dynamic fluoroscopic studies to compare intervertebral angular motion sharing inequality and variability (MSI and MSV) during continuous lumbar motion in CNSLBP patients and controls. Passive recumbent and active standing protocols were used and the relationships of these variables to age and disc degeneration were assessed.

METHODS

Twenty patients with CNSLBP and 20 matched controls received quantitative fluoroscopic lumbar spine examinations using a standardised protocol for data collection and image analysis. Composite disc degeneration (CDD) scores comprising the sum of Kellgren and Lawrence grades from L2-S1 were obtained. Indices of intervertebral motion sharing inequality (MSI) and variability (MSV) were derived and expressed in units of proportion of lumbar range of motion from outward and return motion sequences during lying (passive) and standing (active) lumbar bending and compared between patients and controls. Relationships between MSI, MSV, age and CDD were assessed by linear correlation.

RESULTS

MSI was significantly greater in the patients throughout the intervertebral motion sequences of recumbent flexion (0.29 vs. 0.22, p = 0.02) and when flexion, extension, left and right motion were combined to give a composite measure (1.40 vs. 0.92, p = 0.04). MSI correlated substantially with age (R = 0.85, p = 0.004) and CDD (R = 0.70, p = 0.03) in lying passive investigations in patients and not in controls. There were also substantial correlations between MSV and age (R = 0.77, p = 0.01) and CDD (R = 0.85, p = 0.004) in standing flexion in patients and not in controls.

CONCLUSION

Greater inequality and variability of motion sharing was found in patients with CNSLBP than in controls, confirming previous studies and suggesting a biomechanical marker for the disorder at intervertebral level. The relationship between disc degeneration and MSI was augmented in patients, but not in controls during passive motion and similarly for MSV during active motion, suggesting links between in vivo disc mechanics and pain generation.

Multi-segmental thoracic spine kinematics measured dynamically in the young and elderly during flexion

In contrast to the cervical and lumbar region, the normal kinematics of the thoracic spine have not been thoroughly investigated. The aim of this study was to characterize normal multi-segmental continuous motion of the whole thoracolumbar spine, during a flexion maneuver, in young and elderly subjects. Forty-two healthy volunteers were analyzed: 21 young (age=27.00±3.96) and 21 elderly (age=70.1±3.85). Spinal motion was recorded with a motion-capture system and analyzed using a 3rd order polynomial function to approximate spinal curvature throughout the motion sequence. The average motion profiles of the two age groups were characterized. Flexion timing of the thoracic region of the spine, as compared to the lumbar spine and hips, was found to be different in the two age groups (p=0.011): a delayed/sequential motion type was observed in most of the young, whereas mostly a simultaneous motion pattern was observed in the elderly subjects. A similar trend was observed in flexion of the lower thoracic segments (p=0.017). Differences between age groups were also found for regional and segmental displacements and velocities. The reported characterization of the thoracic spine kinematics may in the future support identification of abnormal movement or be used to improve biomechanical models of the spine.

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.

The role of sacral slope in lumbosacral fusion: a biomechanical study

OBJECT

Abnormal sacral slope (SS) has shown to increase progression of spondylolisthesis, yet there exists a paucity in biomechanical studies investigating its role in the correction of adult spinal deformity, its influence on lumbosacral shear, and its impact on the instrumentation selection process. This in vitro study investigates the effect of SS on 3 anterior lumbar interbody fusion constructs in a biomechanics laboratory.

METHODS

Nine healthy, fresh-frozen, intact human lumbosacral vertebral segments were tested by applying a 550-N axial load to specimens with an initial SS of 20° on an MTS Bionix test system. Testing was repeated as SS was increased to 50°, in 10° increments, through an angulated testing fixture. Specimens were instrumented using a standalone integrated spacer with self-contained screws (SA), an interbody spacer with posterior pedicle screws (PPS), and an interbody spacer with anterior tension band plate (ATB) in a randomized order. Stiffness was calculated from the linear portion of the load-deformation curve. Ultimate strength was also recorded on the final construct of all specimens (n = 3 per construct) with SS of 40°.

RESULTS

Axial stiffness (N/mm) of the L5–S1 motion segment was measured at various angles of SS: for SA 292.9 ± 142.8 (20°), 277.2 ± 113.7 (30°), 237.0 ± 108.7 (40°), 170.3 ± 74.1 (50°); for PPS 371.2 ± 237.5 (20°), 319.8 ± 167.2 (30°), 280.4 ± 151.7 (40°), 233.0 ± 117.6 (50°); and for ATB 323.9 ± 210.4 (20°), 307.8 ± 125.4 (30°), 249.4 ± 126.7 (40°), 217.7 ± 99.4 (50°). Axial compression across the disc space decreased with increasing SS, indicating that SS beyond 40° threshold shifted L5–S1 motion into pure shear, instead of compression-shear, defining a threshold. Trends in ultimate load and displacement differed from linear stiffness with SA > PPS > ATB.

CONCLUSIONS

At larger SSs, bilateral pedicle screw constructs with spacers were the most stable; however, none of the constructs were significantly stiffer than intact segments. For load to failure, the integrated spacer performed the best; this may be due to angulations of integrated plate screws. Increasing SS significantly reduced stiffness, which indicates that surgeons need to consider using more aggressive fixation techniques.

Sagittal plane lumbar intervertebral motion during seated flexion-extension radiographs of 658 asymptomatic nondegenerated levels

OBJECT

Evaluation of lumbar stability is fundamentally dependent on a clear understanding of normal lumbar motion. There are inconsistencies in reported lumbar motion across previously published studies, and it is unclear which provide the most reliable reference data. New technology now allows valid and reliable determination of normal lumbar intervertebral motion (IVM). The object of this study was to provide normative reference data for lumbar IVM and center of rotation (COR) using validated computer-assisted measurement tools.

METHODS

Sitting flexion-extension radiographs were obtained in 162 asymptomatic volunteers and then analyzed using a previously validated and widely used computerized image analysis method. Each lumbar level was subsequently classified as “degenerated” or “nondegenerated” using the Kellgren-Lawrence classification. Of the 803 levels analyzed, 658 were nondegenerated (Kellgren-Lawrence grade < 2). At each level of the lumbar spine, the magnitude of intervertebral rotation and translation, the ratio of translation per degree of rotation (TPDR), and the position of the COR were calculated in the nondegenerative cohort. Translations were calculated in millimeters and percentage endplate width.

RESULTS

All parameters were significantly dependent on the intervertebral level. The upper limit of the 95% CIs for anteroposterior intervertebral translation in this asymptomatic cohort ranged from 2.1 mm (6.2% endplate width) to 4.6 mm (13.3% endplate width). Intervertebral rotation upper limits ranged from 16.3° to 23.5°. The upper limits for TPDR ranged from 0.49% to 0.82% endplate width/degree. The COR coordinates were clustered in level-dependent patterns.

CONCLUSIONS

New normal values for IVM, COR, and the ratio of TPDR in asymptomatic nondegenerative lumbar levels are proposed, providing a reference for future interpretation of sagittal plane motion in the lumbar spine.

Attainment rate as a surrogate indicator of the intervertebral neutral zone length in lateral bending: an in vitro proof of concept study

Background

Lumbar segmental instability is often considered to be a cause of chronic low back pain. However, defining its measurement has been largely limited to laboratory studies. These have characterised segmental stability as the intrinsic resistance of spine specimens to initial bending moments by quantifying the dynamic neutral zone. However these measurements have been impossible to obtain in vivo without invasive procedures, preventing the assessment of intervertebral stability in patients. Quantitative fluoroscopy (QF), measures the initial velocity of the attainment of intervertebral rotational motion in patients, which may to some extent be representative of the dynamic neutral zone. This study sought to explore the possible relationship between the dynamic neutral zone and intervertebral rotational attainment rate as measured with (QF) in an in vitro preparation. The purpose was to find out if further work into this concept is worth pursuing.

Method

This study used passive recumbent QF in a multi-segmental porcine model. This assessed the intrinsic intervertebral responses to a minimal coronal plane bending moment as measured with a digital force guage. Bending moments about each intervertebral joint were calculated and correlated with the rate at which global motion was attained at each intervertebral segment in the first 10° of global motion where the intervertebral joint was rotating.

Results

Unlike previous studies of single segment specimens, a neutral zone was found to exist during lateral bending. The initial attainment rates for left and right lateral flexion were comparable to previously published in vivo values for healthy controls. Substantial and highly significant levels of correlation between initial attainment rate and neutral zone were found for left (Rho = 0.75, P = 0.0002) and combined left-right bending (Rho = 0.72, P = 0.0001) and moderate ones for right alone (Rho = 0.55, P = 0.0012).

Conclusions

This study found good correlation between the initial intervertebral attainment rate and the dynamic neutral zone, thereby opening the possibility to detect segmental instability from clinical studies. However the results must be treated with caution. Further studies with multiple specimens and adding sagittal plane motion are warranted.

Development of a novel radiographic measure of lumbar instability and validation using the facet fluid sign

BACKGROUND

Lumbar spinal instability is frequently referenced in clinical practice and the scientific literature despite the lack of a standard definition or validated radiographic test. The Quantitative Stability Index (QSI) is being developed as a novel objective test for sagittal plane lumbar instability. The QSI is calculated using lumbar flexion-extension radiographs. The goal of the current study was to use the facet fluid sign on MRI as the "gold standard" and determine if the QSI is significantly different in the presence of the fluid sign.

METHODS

Sixty-two paired preoperative MRI and flexion-extension exams were obtained from a large FDA IDE study. The MRI exams were assessed for the presence of a facet fluid sign, and the QSI was calculated from sagittal plane intervertebral rotation and translation measurements. The QSI is based on the translation per degree of rotation (TPDR) and is calculated as a Z-score. A QSI > 2 indicates that the TPDR is > 2 std dev above the mean for an asymptomatic and radiographically normal population. The reproducibility of the QSI was also tested.

RESULTS

The mean difference between trained observers in the measured QSI was between -0.28 and 0.36. The average QSI was significantly (P = 0.047, one-way analysis of variance) higher at levels with a definite fluid sign (2.3±3.2 versus 0.60±2.4).

CONCLUSIONS

Although imperfect, the facet fluid sign observed may be the best currently available test for lumbar spine instability. Using the facet fluid sign as the "gold standard" the current study documents that the QSI can be expected to be significantly higher in the presence of the facet fluid sign. This supports that QSI might be used to test for sagittal plane lumbar instability.

CLINICAL RELEVANCE

A validated, objective and practical test for spinal instability would facilitate research to understand the importance of instability in diagnosis and treatment of low-back related disorders.

Distribution of Schmorl nodes in the lumbar spine and their relationship with lumbar disk degeneration and range of motion

STUDY DESIGN

A kinematic magnetic resonance imaging study.

OBJECTIVE

To investigate the distribution of Schmorl nodes (SNs) in the lumbar spine in healthy adults, and determine the association with lumbar disk degeneration and lumbar spine motion.

SUMMARY OF BACKGROUND DATA

SNs have been associated with several pathologies of the lumbar spine, although it has been demonstrated that they also occur in the healthy adult population without a clearly identified cause. A thorough understanding of SN distribution may help reveal reasons for their formation. How disk degeneration and lumbar spine motion relate to SNs is poorly understood.

METHODS

Kinematic magnetic resonance images (0.6 T) were available for 1179 healthy individuals from 15 to 85 years of age. Spine specialists performed computer-based measurements. All parameters were measured and calculated automatically using the eRAD PACS Viewer (eRAD Inc., version 6.2.1.1). Lumbar disk degeneration was documented according to the Pfirrmann classification system. Lumbar spine lordosis was quantified as the angle between the inferior endplate of L1 and superior endplate of S1. The level of significance was defined as P ≤ 0.05. The distribution of SNs along the lumbar spine and their relationship with age and sex was investigated using the single factor analysis of variance χ test. The relationship between SNs, age group, disk location, and overall grades of lumbar disk degeneration were investigated by multiple logistic regression analysis. Lumbar spine motion was compared between patients with and without SNs via independent t test among 585 individuals with qualified kinematic images. Multiple logistic regression analysis was performed on associations of lumbar motion range among the SN population.

RESULTS

The prevalence of SNs in our study population was 28.4%, and SNs were observed to be present more frequently in males (34.6%) than in females (20.2%) (ρ< 0.01). There was no significant difference in the incidence of SNs between age groups (ρ= 0.18). SNs were more common at the L2 and L3 vertebral bodies (14.3% and 14.4%), whereas SNs were least common at S1 vertebral bodies (1.5%). The highest incidences of SNs presentation was on disks with degeneration grade III (41.9%) and grade IV (45.3%). SN occurrence, aging, and disk location were positively correlated with lumbar disk degeneration grade. The lumbar spine range of motion was significantly different between individuals with and without SNs (31.4° vs. 37.9°, ρ< 0.01). The frequency of SNs was associated with decreased lumbar range of motion in all age groups except 51 to 60 years and 61 to 70 years.

CONCLUSION

SNs have a high incidence in individuals without persistent lumbar disorders and were found in disks at all degrees of degeneration. SNs occurrence were positively associated with lumbar disk degeneration In addition, the presence of SNs was correlated with decreased overall lumbar motion across all age groups.

In vitro biomechanical evaluation of single impulse and repetitive mechanical shockwave devices utilized for spinal manipulative therapy

Mechanical shockwave therapy devices have been in clinical use for almost 40 years. While most often used to treat back pain, our understanding of their biomechanical performance is very limited. From biomechanical studies we know that biological tissue is viscoelastic and preferably excited around its resonance frequency. Targeting these frequencies has been the focus in extracorporeal shock wave lithotripsy, but these concepts are relatively new in orthopedic and rehabilitation therapies. The exact mechanism by which shockwave therapy acts is not known. Knowledge of the performance characteristics of these devices, correlated with clinical outcome studies, may lead to better patient selection, improvement of device functionality, and knowledge of the underlying working principals of therapy. The objectives of this study were to determine the ability of several commercial shockwave devices to achieve a desired thrust profile in a benchtop setting, determine the thrust profile in a clinical analog, and determine the influence of operator experience level on device performance. We conducted two different types of testing: (1) bench testing to evaluate the devices themselves, and (2) clinical equivalent testing to determine the influence of the operator. The results indicated a significant dependence of thrust output on the compliance of the test media. The Activator V-E device matched the ideal half-sine thrust profile to 94%, followed by the Impulse device (84%), the Activator IV/FS (74%), and the Activator II (48%). While most devices deviated from the ideal profile on the return path, the Impulse device exhibited a secondary peak. Moreover, the Activator V-E device provided evidence that the device performs consistently despite operator experience level. This has been a major concern in manual spinal manipulation. Based on our results, a hyper-flexible spine would receive a lower peak thrust force than a hypo-flexible spine at the same power setting. Furthermore, a hand-held operation further reduced the peak thrust force as it increased the system compliance. However, that influence was dissimilar for the different devices. Although controlled clinical trials are needed to determine the correlation between thrust profile and clinical outcome, already ongoing clinical studies indicate an improved patient satisfaction due to reduced treatment pain when devices are used with a thrust characteristic closer to an ideal sine wave.

Electromyographic assessment of trunk and shoulder muscles during a Pilates pull-up exercise

This study compares surface electromyographic activity of the internal oblique, rectus abdominis, multifidus, iliocostalis, anterior deltoids during the pull-up on a lower and on a higher difficulty level. We assessed nine adults with previous experience in Pilates. The root mean square (RMS) values were normalized by maximum isometric contraction for each participant. During the ascent phase, the low spring position showed a significantly higher RMS than the high spring position of 8.9% for deltoid, 17.2% for internal oblique, 22.3% for rectus abdominis, 4.1% for iliocostalis, and 5.6% for multifidus, and in the descent phase, the RMS in the lower spring exceeded significantly the high spring position in 1.6% for the deltoid, 10% for internal oblique, 31.4% for rectus abdominis and 11.4% for iliocostalis. There was no predominance of abdominal muscles over the shoulder muscle in any spring position. The pull-up exercise can be a useful choice for the core and anterior deltoid muscles strengthening.

Proportional lumbar spine inter-vertebral motion patterns: a comparison of patients with chronic, non-specific low back pain and healthy controls

INTRODUCTION

Identifying biomechanical subgroups in chronic, non-specific low back pain (CNSLBP) populations from inter-vertebral displacements has proven elusive. Quantitative fluoroscopy (QF) has excellent repeatability and provides continuous standardised inter-vertebral kinematic data from fluoroscopic sequences allowing assessment of mid-range motion. The aim of this study was to determine whether proportional continuous IV rotational patterns were different in patients and controls. A secondary aim was to update the repeatability of QF measurement of range of motion (RoM) for inter-vertebral (IV) rotation.

METHODS AND MATERIALS

Fluoroscopic sequences were recorded of passive, recumbent coronal and sagittal motion, which was controlled for range and velocity. Segments L2-5 in 40 primary care CNSLBP patients and 40 matched controls were compared. Patients also completed the von Korff Chronic Pain Grade and Roland and Morris Disability Questionnaire. Sequences were processed using automated image tracking algorithms to extract continuous inter-vertebral rotation data. These were converted to continuous proportional ranges of rotation (PR). The continuous proportional range variances were calculated for each direction and combined to produce a single variable representing their fluctuation (CPRV). Inter- and intra-rater repeatability were also calculated for the maximum IV-RoM measurements obtained during controlled trunk motion to provide an updated indication of the reliability and agreement of QF for measuring spine kinematics.

RESULTS

CPRV was significantly higher in patients (0.011 vs. 0.008, Mann-Whitney two-sided p = 0.008), implying a mechanical subgroup. Receiver operating characteristic curve analysis found its sensitivity and specificity to be 0.78 % (60-90) and 0.55 % (37-73), respectively (area under the curve 0.672). CPRV was not correlated with pain severity or disability. The repeatability of maximum inter-vertebral range was excellent, but range was only significantly greater in patients at L4-5 in right side bending (p = 0.03).

CONCLUSION

The variation in proportional motion between lumbar vertebrae during passive recumbent trunk motion was greater in patients with CNSLBP than in matched healthy controls, indicating that biomechanical factors in passive structures play a part.

Capturing Three-Dimensional In Vivo Lumbar Intervertebral Joint Kinematics Using Dynamic Stereo-X-Ray Imaging

Availability of accurate three-dimensional (3D) kinematics of lumbar vertebrae is necessary to understand normal and pathological biomechanics of the lumbar spine. Due to the technical challenges of imaging the lumbar spine motion in vivo, it has been difficult to obtain comprehensive, 3D lumbar kinematics during dynamic functional tasks. The present study demonstrates a recently developed technique to acquire true 3D lumbar vertebral kinematics, in vivo, during a functional load-lifting task. The technique uses a high-speed dynamic stereo-radiography (DSX) system coupled with a volumetric model-based bone tracking procedure. Eight asymptomatic male participants performed weight-lifting tasks, while dynamic X-ray images of their lumbar spines were acquired at 30 fps. A custom-designed radiation attenuator reduced the radiation white-out effect and enhanced the image quality. High resolution CT scans of participants' lumbar spines were obtained to create 3D bone models, which were used to track the X-ray images via a volumetric bone tracking procedure. Continuous 3D intervertebral kinematics from the second lumbar vertebra (L2) to the sacrum (S1) were derived. Results revealed motions occurring simultaneously in all the segments. Differences in contributions to overall lumbar motion from individual segments, particularly L2–L3, L3–L4, and L4–L5, were not statistically significant. However, a reduced contribution from the L5–S1 segment was observed. Segmental extension was nominally linear in the middle range (20%–80%) of motion during the lifting task, but exhibited nonlinear behavior at the beginning and end of the motion. L5–S1 extension exhibited the greatest nonlinearity and variability across participants. Substantial AP translations occurred in all segments (5.0 ± 0.3 mm) and exhibited more scatter and deviation from a nominally linear path compared to segmental extension. Maximum out-of-plane rotations (<1.91 deg) and translations (<0.94 mm) were small compared to the dominant motion in the sagittal plane. The demonstrated success in capturing continuous 3D in vivo lumbar intervertebral kinematics during functional tasks affords the possibility to create a baseline data set for evaluating the lumbar spinal function. The technique can be used to address the gaps in knowledge of lumbar kinematics, to improve the accuracy of the kinematic input into biomechanical models, and to support development of new disk replacement designs more closely replicating the natural lumbar biomechanics.

The LP-ESP(®) lumbar disc prosthesis with 6 degrees of freedom: development and 7 years of clinical experience

The viscoelastic lumbar disk prosthesis-elastic spine pad (LP-ESP(®)) is an innovative one-piece deformable but cohesive interbody spacer providing 6 full degrees of freedom about the 3 axes, including shock absorption. A 20-year research program has demonstrated that this concept provides mechanical properties very close to those of a natural disk. Improvements in technology have made it possible to solve the problem of the bond between the elastic component and the titanium endplates and to obtain an excellent biostability. The prosthesis geometry allows limited rotation and translation with resistance to motion (elastic return property) aimed at avoiding overload of the posterior facets. The rotation center can vary freely during motion. It thus differs substantially from current prostheses, which are 2- or 3-piece devices involving 1 or 2 bearing surfaces and providing 3 or 5 degrees of freedom. This design and the adhesion-molding technology differentiate the LP-ESP prosthesis from other mono-elastomeric prostheses, for which the constraints of shearing during rotations or movement are absorbed at the endplate interface. Seven years after the first implantation, we can document in a solid and detailed fashion the course of clinical outcomes and the radiological postural and kinematic behavior of this prosthesis.

Development and kinematic verification of a finite element model for the lumbar spine: application to disc degeneration

The knowledge of the lumbar spine biomechanics is essential for clinical applications. Due to the difficulties to experiment on living people and the irregular results published, simulation based on finite elements (FE) has been developed, making it possible to adequately reproduce the biomechanics of the lumbar spine. A 3D FE model of the complete lumbar spine (vertebrae, discs, and ligaments) has been developed. To verify the model, radiological images (X-rays) were taken over a group of 25 healthy, male individuals with average age of 27.4 and average weight of 78.6 kg with the corresponding informed consent. A maximum angle of 34.40° is achieved in flexion and of 35.58° in extension with a flexion-extension angle of 69.98°. The radiological measurements were 33.94 ± 4.91°, 38.73 ± 4.29°, and 72.67°, respectively. In lateral bending, the maximum angles were 19.33° and 23.40 ± 2.39, respectively. In rotation a maximum angle of 9.96° was obtained. The model incorporates a precise geometrical characterization of several elements (vertebrae, discs, and ligaments), respecting anatomical features and being capable of reproducing a wide range of physiological movements. Application to disc degeneration (L5-S1) allows predicting the affection in the mobility of the different lumbar segments, by means of parametric studies for different ranges of degeneration.

Validity and reliability of skin markers for measurement of intersegmental mobility at L2-3 and L3-4 during lateral bending in healthy individuals: a fluoroscopy study

It is clinically important to assess kinematic parameters of lumbar spine movement to increase our understanding of lumbar mobility impairments in patients with low back pain. This is the first step for restoration of motor function. The use of non-invasive surface markers has currently attracted the interests of many researchers but scientific utilization of this technique for clinical research requires validity and reliability studies. The aim of the present study was to examine whether skin markers can be used to measure lumbar motions during lateral bending. Twelve healthy individuals were lying in prone position on the video fluoroscopy table and skin markers were attached over their spinous processes. Fluoroscopy images were taken in two positions of neutral and right lateral bending (RLB). The correlation of the L2-3 and L3-4 angles and lumbar curvature between markers and vertebrae measurements in the neutral and RLB positions was determined by Pearson Correlation Coefficient. The Intraclass correlation coefficient (ICC) was used to measure inter-examiner reliability of measurement in five participants. The results showed high reliability (ranging from 0.94 to 0.99) for angular measurements at L2-3 and L3-4 and lumbar curvature and also significant correlation between angular measurement derived from markers and vertebrae at L2-3 (r = 0.7, p = 0.015), L3-4 and lumbar curvature (r = 0.91 p = 0.001). The results showed that motions of skin markers follow the motions of the assigned underlying lumbar vertebrae. Therefore, skin markers can be confidently used for estimation of lumbar movements during lateral bending.

Measurement of intervertebral motion using quantitative fluoroscopy: report of an international forum and proposal for use in the assessment of degenerative disc disease in the lumbar spine

Quantitative fluoroscopy (QF) is an emerging technology for measuring intervertebral motion patterns to investigate problem back pain and degenerative disc disease. This International Forum was a networking event of three research groups (UK, US, Hong Kong), over three days in San Francisco in August 2009. Its aim was to reach a consensus on how best to record, analyse, and communicate QF information for research and clinical purposes. The Forum recommended that images should be acquired during regular trunk motion that is controlled for velocity and range, in order to minimise externally imposed variability as well as to correlate intervertebral motion with trunk motion. This should be done in both the recumbent passive and weight bearing active patient configurations. The main recommended outputs from QF were the true ranges of intervertebral rotation and translation, neutral zone laxity and the consistency of shape of the motion patterns. The main clinical research priority should initially be to investigate the possibility of mechanical subgroups of patients with chronic, nonspecific low back pain by comparing their intervertebral motion patterns with those of matched healthy controls.

Effects of lumbar artificial disc design on intervertebral mobility: in vivo comparison between mobile-core and fixed-core

Although in theory, the differences in design between fixed-core and mobile-core prostheses should influence motion restoration, in vivo kinematic differences linked with prosthesis design remained unclear. The aim of this study was to investigate the rationale that the mobile-core design seems more likely to restore physiological motion since the translation of the core could help to mimic the kinematic effects of the natural nucleus. In vivo intervertebral motion characteristics of levels implanted with the mobile-core prosthesis were compared with untreated levels of the same population, levels treated by a fixed-core prosthesis, and normal levels (data from literature). Patients had a single-level implantation at L4L5 or L5S1 including 72 levels with a mobile-core prosthesis and 33 levels with a fixed-core prosthesis. Intervertebral mobility characteristics included the range of motion (ROM), the motion distribution between flexion and extension, the prosthesis core translation (CT), and the intervertebral translation (VT). A method adapted to the implanted segments was developed to measure the VT: metal landmarks were used instead of the bony landmarks. The reliability assessment of the VT measurement method showed no difference between three observers (p < 0.001), a high level of agreement (ICC = 0.908) and an interobserver precision of 0.2 mm. Based on this accurate method, this in vivo study demonstrated that the mobile-core prosthesis replicated physiological VT at L4L5 levels but not at L5S1 levels, and that the fixed-core prosthesis did not replicate physiological VT at any level by increasing VT. As the VT decreased when the CT increased (p < 0.001) it was proven that the core mobility minimized the VT. Furthermore, some physiologic mechanical behaviors seemed to be maintained: the VT was higher at implanted the L4L5 level than at the implanted L5S1 level, and the CT appeared lower at the L4L5 level than at the L5S1 level. ROM and motion distribution were not different between the mobile-core prosthesis and the fixed-core prosthesis implanted levels. This study validated in vivo the concept that a mobile-core helps to restore some physiological mechanical characteristics of the VT at the implanted L4L5 level, but also showed that the minimizing effect of core mobility on the VT was not sufficient at the L5S1 level.

Kinematic analysis of dynamic lumbar motion in patients with lumbar segmental instability using digital videofluoroscopy

The study design is a prospective, case-control. The aim of this study was to develop a reliable measurement technique for the assessment of lumbar spine kinematics using digital video fluoroscopy in a group of patients with low back pain (LBP) and a control group. Lumbar segmental instability (LSI) is one subgroup of nonspecific LBP the diagnosis of which has not been clarified. The diagnosis of LSI has traditionally relied on the use of lateral functional (flexion-extension) radiographs but use of this method has proven unsatisfactory.Fifteen patients with chronic low back pain suspected to have LSI and 15 matched healthy subjects were recruited. Pulsed digital videofluoroscopy was used to investigate kinematics of lumbar motion segments during flexion and extension movements in vivo. Intersegmental linear translation and angular displacement, and pathway of instantaneous center of rotation (PICR) were calculated for each lumbar motion segment. Movement pattern of lumbar spine between two groups and during the full sagittal plane range of motion were analyzed using ANOVA with repeated measures design. Intersegmental linear translation was significantly higher in patients during both flexion and extension movements at L5-S1 segment (p < 0.05). Arc length of PICR was significantly higher in patients for L1-L2 and L5-S1 motion segments during extension movement (p < 0.05). This study determined some kinematic differences between two groups during the full range of lumbar spine. Devices, such as digital videofluoroscopy can assist in identifying better criteria for diagnosis of LSI in otherwise nonspecific low back pain patients in hope of providing more specific treatment.