Semiautomated 3D Spine Reconstruction from Biplanar Radiographic Images: Prediction of Intervertebral Loading in Scoliotic Subjects

Dataset of Finite Element Models of Normal and Deformed Thoracolumbar Spine

Adult spine deformity (ASD) is prevalent and leads to a sagittal misalignment in the vertebral column. Computational methods, including Finite Element (FE) Models, have emerged as valuable tools for investigating the causes and treatment of ASD through biomechanical simulations. However, the process of generating personalised FE models is often complex and time-consuming. To address this challenge, we present a dataset of FE models with diverse spine morphologies that statistically represent real geometries from a cohort of patients. These models are generated using EOS images, which are utilized to reconstruct 3D surface spine models. Subsequently, a Statistical Shape Model (SSM) is constructed, enabling the adaptation of a FE hexahedral mesh template for both the bone and soft tissues of the spine through mesh morphing. The SSM deformation fields facilitate the personalization of the mean hexahedral FE model based on sagittal balance measurements. Ultimately, this new hexahedral SSM tool offers a means to generate a virtual cohort of 16807 thoracolumbar FE spine models, which are openly shared in a public repository.

Reliability of semi-automated spinal measurement software

PURPOSE

In the treatment of patients with adult spinal deformity, analysis of spinopelvic balance is essential in clinical assessment and surgical planning. There is currently no gold standard for measurement, whether done by hand or with digital software. New semi-automated software exists that purports to increase efficiency, but its reliability is unknown in the literature.

METHODS

Full spine X-rays were retrospectively reviewed from 25 adult patients seen between 2014 and 2017. Patients were included if they had > 5 cm of sagittal imbalance and radiographs of sufficient quality to perform balance measurements, without prior surgical spinal fusion and/or instrumentation. Spinopelvic parameters were measured in two radiographic programs: one with basic, non-spine-specific measurement tools (eUnity, Client Outlook, Waterloo, Canada); and a second with spine-specific semi-automated measurement tools (Sectra, Sectra AB, Linköping, Sweden). Balance parameters included SVA, PI, PT, and LL. Data were compared by examining inter-rater and inter-program reliability using interclass correlation coefficient (ICC).

RESULTS

The subjects' mean age was 67.9 ± 13.8 years old, and 32% were male. The inter-program reliability was strong, with ICC values greater than 0.91 for each parameter. Similarly, there was strong inter-observer reliability with ICC values greater than 0.88. These results persisted on delayed repeat measurement (p < 0.001 for all measurements).

CONCLUSION

There is excellent inter-observer and inter-program reliability between the basic PACS and semi-automated programs. These data demonstrate that the purported efficiency of semi-automated measurement programs does not come at the cost of measurement reliability.

Prediction of trunk muscle activation and spinal forces in adolescent idiopathic scoliosis during simulated trunk motion: A musculoskeletal modelling study

Due to lack of reference validation data, the common strategy in characterizing adolescent idiopathic scoliosis (AIS) by musculoskeletal modelling approach consists in adapting structure and parameters of validated body models of adult individuals with physiological alignments. Until now, only static postures have been replicated and investigated in AIS subjects. When aiming to simulate trunk motion, two critical factors need consideration: how distributing movement along the vertebral motion levels (lumbar spine rhythm), and if neglecting or accounting for the contribution of the stiffness of the motion segments (disc stiffness). The present study investigates the effect of three different lumbar spine rhythms and absence/presence of disc stiffness on trunk muscle imbalance in the lumbar region and on intervertebral lateral shear at different levels of the thoracolumbar/lumbar scoliotic curve, during simulated trunk motions in the three anatomical planes (flexion/extension, lateral bending, and axial rotation). A spine model with articulated ribcage previously developed in AnyBody software and adapted to replicate the spinal alignment in AIS subjects is employed. An existing dataset of 100 subjects with mild and moderate scoliosis is exploited. The results pointed out the significant impact of lumbar spine rhythm configuration and disc stiffness on changes in the evaluated outputs, as well as a relationship with scoliosis severity. Unfortunately, no optimal settings can be identified due to lack of reference validation data. According to that, extreme caution is recommended when aiming to adapt models of adult individuals with physiological alignments to adolescent subjects with scoliotic deformity.

Multifidus Denervation After Radiofrequency Ablation of the Medial Nerve Alters the Biomechanics of the Spine-A Computational Study

Radiofrequency ablation of the medial branch is commonly used to treat chronic low back pain involving facet joints, which accounts for 12% to 37% of the total cases of chronic low back pain. An adverse effect of this procedure is the denervation of the multifidus muscle, which may lead to its atrophy which can affect the spine and possibly disc degeneration. This study aims to quantify changes in joint angles and loading caused by multifidus denervation after radiofrequency ablation. AnyBody model of the torso was used to evaluate intervertebral joints in flexion, lateral bending, and torsion. Force-dependent kinematics was used to calculate joint angles and forces. These dependent variables were investigated in intact multifidus, unilateral, and bilateral ablations of L3L4, L4L5, and L5S1 joints. The results showed pronounced angular joint changes, especially in bilateral ablations in flexion, when compared with other cases. The same changes' trend from intact to unilaterally then bilaterally ablated multifidus occurred in joint angles of lateral bending. Meanwhile, joint forces were not adversely affected. These results suggest that multifidus denervation after radiofrequency ablation affects spinal mechanics. Such changes may be associated with abnormal tissue deformations and stresses that can potentially alter their mechanobiology and homeostasis, thereby possibly affecting the health of the spine.

Effect of Unilateral Knee Extension Restriction on the Lumbar Region during Gait

Unilateral knee extension restriction might change trunk alignment and increase mechanical load on the lumbar region during walking. We aimed to clarify lumbar region mechanical load during walking with restricted knee extension using a musculoskeletal model simulation. Seventeen healthy adult males were enrolled in this study. Participants walked 10 m at a comfortable velocity with and without restricted right knee extension of 15° and 30° using a knee brace. L4-5 joint moment, joint reaction force, and muscle forces around the lumbar region during walking were calculated for each condition. Peaks of kinetic data were compared among three gait conditions during 0%-30% and 50%-80% of the right gait cycle. Lumbar extension moment at early stance of the bilateral lower limbs was significantly increased in the 30° restricted condition (p ≤ 0.021). Muscle force of the multifidus showed peaks at stance phase of the contralateral side during walking, and the erector spinae showed force peaks at early stance of the bilateral lower limb. Muscle force of the multifidus and erector spinae increased with increasing degree of knee flexion (p ≤ 0.010), with a large effect size (η ² = 0.273-0.486). The joint force acting on L4-5 showed two peaks at early stance of the bilateral lower limbs during the walking cycle. The anterior and vertical joint force on L4-5 increased by 14.2%-36.5% and 10.0%-23.0% in walking with restricted knee extension, respectively (p ≤ 0.010), with a large effect size (η ² = 0.149-0.425). Restricted knee joint extension changed trunk alignment and increased the muscle force and the vertical and anterior joint force on the L4-5 joint during walking; this tendency became more obvious with increased restriction angle. Our results provide important information for therapists engaged in the rehabilitation of patients with knee contracture.

Model-Based Radiation-Free Assessment of Scoliosis: A Principle Validation Study

Purpose

Adolescent scoliosis is one of the common pediatric spinal diseases which has a high risk of progression due to the rapid growth of the skeleton during the growing stage therefore needs regular clinical monitoring including X-rays. Because X-rays could lead to ionizing radiation-related health problems, an ionizing radiation-free, non-invasive method is presented here to estimate the degree of scoliosis and to potentially support the medical assessment.

Methods

The radiation-free body scanner provides a 3D surface scan of the torso. A basic 3D structure of the human ribcage and vertebral column was modeled and simulated with computer-aided design software and finite element method calculation. For comparison with X-rays, courses of vertebral columns derived from 3D torso images and 3D models were analyzed with respect to their apex positions and angles.

Results

The methods show good results in the estimation of the apex positions of scoliosis. Strong correlations (R = 0.8924) were found between the apex and Cobb angle from X-rays. Similar correlations (R = 0.8087) was obtained between the apex angles extracted from X-rays and the combination of torso scan images with 3D model simulations. Promising agreement was obtained between the spinal trajectories extracted from X-ray and 3D torso images.

Conclusions

Very strong correlations suggest that the apex angle could potentially be used for scoliosis assessment in follow-up examinations in complement to the Cobb angle. However, further improvements of the methods and tests on a larger number of data set are necessary before their introduction into the clinical application.

Computational modelling of the scoliotic spine: A literature review

Scoliosis is a deformity of the spine that in severe cases requires surgical treatment. There is still disagreement among clinicians as to what the aim of such treatment is as well as the optimal surgical technique. Numerical models can aid clinical decision-making by estimating the outcome of a given surgical intervention. This paper provided some background information on the modelling of the healthy spine and a review of the literature on scoliotic spine models, their validation, and their application. An overview of the methods and techniques used to construct scoliotic finite element and multibody models was given as well as the boundary conditions used in the simulations. The current limitations of the models were discussed as well as how such limitations are addressed in non-scoliotic spine models. Finally, future directions for the numerical modelling of scoliosis were addressed.

Accounting for Biomechanical Measures from Musculoskeletal Simulation of Upright Posture Does Not Enhance the Prediction of Curve Progression in Adolescent Idiopathic Scoliosis

A major clinical challenge in adolescent idiopathic scoliosis (AIS) is the difficulty of predicting curve progression at initial presentation. The early detection of progressive curves can offer the opportunity to better target effective non-operative treatments, reducing the need for surgery and the risks of related complications. Predictive models for the detection of scoliosis progression in subjects before growth spurt have been developed. These models accounted for geometrical parameters of the global spine and local descriptors of the scoliotic curve, but neglected contributions from biomechanical measurements such as trunk muscle activation and intervertebral loading, which could provide advantageous information. The present study exploits a musculoskeletal model of the thoracolumbar spine, developed in AnyBody software and adapted and validated for the subject-specific characterization of mild scoliosis. A dataset of 100 AIS subjects with mild scoliosis and in pre-pubertal age at first examination, and recognized as stable (60) or progressive (40) after at least 6-months follow-up period was exploited. Anthropometrical data and geometrical parameters of the spine at first examination, as well as biomechanical parameters from musculoskeletal simulation replicating relaxed upright posture were accounted for as predictors of the scoliosis progression. Predicted height and weight were used for model scaling because not available in the original dataset. Robust procedure for obtaining such parameters from radiographic images was developed by exploiting a comparable dataset with real values. Six predictive modelling approaches based on different algorithms for the binary classification of stable and progressive cases were compared. The best fitting approaches were exploited to evaluate the effect of accounting for the biomechanical parameters on the prediction of scoliosis progression. The performance of two sets of predictors was compared: accounting for anthropometrical and geometrical parameters only; considering in addition the biomechanical ones. Median accuracy of the best fitting algorithms ranged from 0.76 to 0.78. No differences were found in the classification performance by including or neglecting the biomechanical parameters. Median sensitivity was 0.75, and that of specificity ranged from 0.75 to 0.83. In conclusion, accounting for biomechanical measures did not enhance the prediction of curve progression, thus not supporting a potential clinical application at this stage.

Subject-Specific Alignment and Mass Distribution in Musculoskeletal Models of the Lumbar Spine

Musculoskeletal modeling is a well-established method in spine biomechanics and generally employed for investigations concerning both the healthy and the pathological spine. It commonly involves inverse kinematics and optimization of muscle activity and provides detailed insight into joint loading. The aim of the present work was to develop and validate a procedure for the automatized generation of semi-subject-specific multi-rigid body models with an articulated lumbar spine. Individualization of the models was achieved with a novel approach incorporating information from annotated EOS images. The size and alignment of bony structures, as well as specific body weight distribution along the spine segments, were accurately reproduced in the 3D models. To ensure the pipeline’s robustness, models based on 145 EOS images of subjects with various weight distributions and spinopelvic parameters were generated. For validation, we performed kinematics-dependent and segment-dependent comparisons of the average joint loads obtained for our cohort with the outcome of various published in vivo and in situ studies. Overall, our results agreed well with literature data. The here described method is a promising tool for studying a variety of clinical questions, ranging from the evaluation of the effects of alignment variation on joint loading to the assessment of possible pathomechanisms involved in adjacent segment disease.

Spinal Compressive Forces in Adolescent Idiopathic Scoliosis With and Without Carrying Loads: A Musculoskeletal Modeling Study

The pathomechanisms of curve progression in adolescent idiopathic scoliosis (AIS) remain poorly understood and biomechanical data are limited. A deeper insight into spinal loading could provide valuable information toward the improvement of current treatment strategies. This work therefore aimed at using subject-specific musculoskeletal full-body models of patients with AIS to predict segmental compressive forces around the curve apex and to investigate how these forces are affected by simulated load carrying. Models were created based on spatially calibrated biplanar radiographic images from 24 patients with mild to moderate AIS and validated by comparing predictions of paravertebral muscle activity with reported values from in vivo studies. Spinal compressive forces were predicted during unloaded upright standing as well as standing with external loads of 10, 15, and 20% of body weight (BW) applied to the scapulae to simulate carrying a backpack in the regular way on the back as well as in front of the body and over the shoulder on the concave and convex sides of the scoliotic curve. The predicted muscle activities around the curve apex were higher on the convex side for the erector spinae (ES) and multifidi (MF) muscles, which was comparable to the EMG-based in vivo measurements from the literature. In terms of spinal loading, the implementation of spinal deformity resulted in a 10% increase of compressive force at the curve apex during unloaded upright standing. Apical compressive forces further increased by 50–62% for a simulated 10% BW load and by 77–94% and 103–128% for 15% and 20% BW loads, respectively. Moreover, load-dependent compressive force increases were the lowest in the regular backpack and the highest in the frontpack and convex conditions, with concave side-carrying forces in between. The predictions indicated increased segmental compressive forces during unloaded upright standing, which could be ascribed to the scoliotic deformation. When carrying loads, compressive forces further increased depending on the carrying mode and the weight of the load. These results can be used as a basis for further studies investigating segmental loading in AIS patients during functional activities. Models can thereby be created using the same approach as proposed in this study.

Dependence of lumbar loads on spinopelvic sagittal alignment: An evaluation based on musculoskeletal modeling

Still little is known about how spinopelvic alignment affects spinal load distribution. Musculoskeletal modeling can potentially help to discover associations between spine alignment and risk factors of spinal disorders (e.g. disc herniation, vertebral fracture, spondylolisthesis, low back pain). The present study exploited the AnyBody full-body musculoskeletal model to assess the relation between lumbar loads and spinopelvic alignment in the sagittal plane. The model was evaluated in the standing position. The simulated postures were set using spinopelvic parameters gleaned from the literature and characterizing the healthy adult population. The parameters were: sagittal vertical axis, Roussouly lumbar type, sacral slope, and pelvic incidence. A total of 2772 configurations were simulated based on the following measurements: compression force and anterior shear at levels L4L5 and L5S1; multifidus, longissimus spinae, and rectus abdominis muscle forces. Changes in global sagittal alignment, lumbar typology, and sacral inclination, but not in pelvic incidence, were found to affect intervertebral loads in the lumbar spine and spinal muscle activation. Considering these changes would be advantageous for clinical evaluation, due to the recognized relation between altered loads and risk of disc herniation, vertebral fracture, spondylolisthesis, and low back pain. Musculoskeletal modeling proved to be a valuable biomechanical tool to non-invasively investigate the relation between internal loads and anatomical parameters.