The effect of three-dimensional geometrical changes during adolescent growth on the biomechanics of a spinal motion segment

Implications of range of motion requirements for the laxity of ligaments in a lumbar finite element model

Finite element models of the lumbar spine often adopt ligament properties from tensile tests without accounting for possible differences between testing and in situ initial ligament length. Such differences could result in laxities or preloads at the beginning of a simulation that would affect the ligament forces, tangent stiffness, and the posture at which they fail. In vivo and in vitro human experimental data reported laxities or preloads. However, laxities or preloads, which could also result from postural differences, are often neglected in simulation studies. This study proposes a numerical methodology to identify ranges of ligament laxities or preloads compatible with the selected tensile ligament properties, the model, and the range of motion (RoM) the model aims to simulate. The approach assumes that ligaments should remain in a safe elongation range for the complete RoM, and that each ligament should play a significant mechanical role in at least one load case. The methodology was applied to the functional spinal unit (FSU) models using the RoM from healthy subjects and ligament properties from the literature. Without laxity, some ligaments reached their elongation at failure within the RoM. Laxity ranges varied considerably (from -9.2 mm preload to 10.7 mm laxity) and flexion was the most critical load case to determine them. Their effect on the mobility response was also assessed. The effect on the mobility of a FSU was also assessed. While the proposed method cannot determine an exact laxity value, it is simple and it can be applied to any model to identify a plausible range of ligament initial length.

Computational lumbar spine models: A literature review

BACKGROUND

Computational spine models of various types have been employed to understand spine function, assess the risk that different activities pose to the spine, and evaluate techniques to prevent injury. The areas in which these models are applied has expanded greatly, potentially beyond the appropriate scope of each, given their capabilities. A comprehensive understanding of the components of these models provides insight into their current capabilities and limitations.

METHODS

The objective of this review was to provide a critical assessment of the different characteristics of model elements employed across the spectrum of lumbar spine modeling and in newer combined methodologies to help better evaluate existing studies and delineate areas for future research and refinement.

FINDINGS

A total of 155 studies met selection criteria and were included in this review. Most current studies use either highly detailed Finite Element models or simpler Musculoskeletal models driven with in vivo data. Many models feature significant geometric or loading simplifications that limit their realism and validity. Frequently, studies only create a single model and thus can't account for the impact of subject variability. The lack of model representation for certain subject cohorts leaves significant gaps in spine knowledge. Combining features from both types of modeling could result in more accurate and predictive models.

INTERPRETATION

Development of integrated models combining elements from different model types in a framework that enables the evaluation of larger populations of subjects could address existing voids and enable more realistic representation of the biomechanics of the lumbar spine.

Current models to understand the onset and progression of scoliotic deformities in adolescent idiopathic scoliosis: a systematic review

PURPOSE

To create an updated and comprehensive overview of the modeling studies that have been done to understand the mechanics underlying deformities of adolescent idiopathic scoliosis (AIS), to predict the risk of curve progression and thereby substantiate etiopathogenetic theories.

METHODS

In this systematic review, an online search in Scopus and PubMed together with an analysis in secondary references was done, which yielded 86 studies. The modeling types were extracted and the studies were categorized accordingly.

RESULTS

Animal modeling, together with machine learning modeling, forms the category of black box models. This category is perceived as the most clinically relevant. While animal models provide a tangible idea of the biomechanical effects in scoliotic deformities, machine learning modeling was found to be the best curve-progression predictor. The second category, that of artificial models, has, just as animal modeling, a tangible model as a result, but focusses more on the biomechanical process of the scoliotic deformity. The third category is formed by computational models, which are very popular in etiopathogenetic parameter-based studies. They are also the best in calculating stresses and strains on vertebrae, intervertebral discs, and other surrounding tissues.

CONCLUSION

This study presents a comprehensive overview of the current modeling techniques to understand the mechanics of the scoliotic deformities, predict the risk of curve progression in AIS and thereby substantiate etiopathogenetic theories. Although AIS remains to be seen as a complex and multifactorial problem, the progression of its deformity can be predicted with good accuracy. Modeling of AIS develops rapidly and may lead to the identification of risk factors and mitigation strategies in the near future. The overview presented provides a basis to follow this development.

Induction of a representative idiopathic-like scoliosis in a porcine model using a multidirectional dynamic spring-based system

BACKGROUND CONTEXT

Scoliosis is a 3D deformity of the spine in which vertebral rotation plays an important role. However, no treatment strategy currently exists that primarily applies a continuous rotational moment over a long period of time to the spine, while preserving its mobility. We developed a dynamic, torsional device that can be inserted with standard posterior instrumentation. The feasibility of this implant to rotate the spine and preserve motion was tested in growing mini-pigs.

PURPOSE

To test the quality and feasibility of the torsional device to induce the typical axial rotation of scoliosis while maintaining growth and mobility of the spine.

STUDY DESIGN

Preclinical animal study with 14 male, 7 month old Gottingen mini-pigs. Comparison of two scoliosis induction methods, with and without the torsional device, with respect to 3D deformity and maintenance of the scoliosis after removal of the implants.

METHODS

Fourteen mini-pigs received either a unilateral tether-only (n=6) or a tether combined with a contralateral torsional device (n=8). X-rays and CT-scans were made post-operative, at 8 weeks and at 12 weeks. Flexibility of the spine was assessed at 12 weeks. In 3 mini-pigs per condition, the implants were removed and the animals were followed until no further correction was expected.

RESULTS

At 12 weeks the tether-only group yielded a coronal Cobb angle of 16.8±3.3°For the tether combined with the torsional device this was 22.0±4.0°. The most prominent difference at 12 weeks was the axial rotation with 3.6±2.8° for the tether-only group compared to 18.1±4.6° for the tether-torsion group. Spinal growth and flexibility remained normal and comparable for both groups. After removal of the devices, the induced scoliosis reduced by 41% in both groups. There were no adverse tissue reactions, implant complications or infections.

CONCLUSION

The present study indicates the ability of the torsional device combined with a tether to induce a flexible idiopathic-like scoliosis in mini-pigs. The torsional device was necessary to induce the typical axial rotation found in human scoliosis.

CLINICAL SIGNIFICANCE

The investigated torsional device could induce apical rotation in a flexible and growing spine. Whether this may be used to reduce a scoliotic deformity remains to be investigated.

Morphology and growth of the pediatric lumbar vertebrae

BACKGROUND CONTEXT

The majority of existing literature describing pediatric lumbar vertebral morphology are limited to characterization of the vertebral bodies, pedicles, and spinal canal and no study has described the rates of growth for any lumbar vertebral structure. While it is known that growth of the lumbar vertebrae results in changes in vertebral shape, the dimension ratios used to quantify these shape changes do not represent the 3D morphology of the vertebral structures. Additionally, many of the previous evaluations of growth and shape are purely descriptive and do not investigate sexual dimorphism or variations across vertebral levels.

PURPOSE

This study aims to establish a database of pediatric lumbar vertebra dimension, growth, and shape data for subjects between and ages of 1 and 19 years.

STUDY DESIGN

A retrospective study of computed tomography (CT) data.

METHODS

Retrospective, abdominal, CT scans of 102 skeletally normal pediatric subjects (54 males, 48 females) between the ages of 1 and 19 years were digitally reconstructed and manually segmented. Thirty surface landmark points (LMPs), 30 vertebral measurements, the centroid size, centroid location, and the local orientation were collected for each lumbar vertebra along with the centroid size of the LMPs comprising each subject's full lumbar spine and their intervertebral disc (IVD) heights. Nonparametric statistics were used to compare dimension values across vertebral levels and between sexes. Linear models with age as the independent variable were used to characterize dimension growth for each sex and vertebral level. Age-dependent quadratic equations were fit to LMP distributions resulting from a generalized Procrustes analysis (GPA) of the vertebrae and fixed effects models were used to investigate differences in model coefficients across levels and between sexes.

RESULTS

Intervertebral level dimension differences were observed across all vertebral structures in both sexes while pedicle widths and IVDs heights were the only measurements found to be sexually dimorphic. Dimension growth rates generally varied across vertebral levels and the growth rates of males were typically larger than those of females. Differences between male and female vertebral shapes were also found for all lumbar vertebral structures.

CONCLUSIONS

To the authors' knowledge, this is the first study to report growth rates for the majority of pediatric lumbar vertebral structures and the first to describe the 3D age-dependent shapes of the pediatric lumbar spine and vertebrae. In addition to providing a quantitative database, the dimension, growth, and shape data reported here would have applications in medical device design, surgical planning, surgical training, and biomechanical modeling.

Electrospun biodegradable poly(ε-caprolactone) membranes for annulus fibrosus repair: Long-term material stability and mechanical competence

BACKGROUND

Electrospun (ES) poly(ɛ-caprolactone) (PCL) is widely used to provide critical mechanical support in tissue engineering and regenerative medicine applications. Therefore, there is a clear need for understanding the change in the mechanical response of the membranes as the material degrades in physiological conditions.

STUDY DESIGN

ES membranes with fiber diameters from 1.6 to 6.7 μm were exposed to in vitro conditions at 37°C in Dulbecco's modified Eagle's medium (DMEM) or dry for up to 6 months.

METHODS

During this period, the mechanical properties were assessed using cyclic mechanical loading, and material properties such as crystallinity and ester bond degradation were measured.

RESULTS

No significant difference was found for any parameters between samples kept dry and in DMEM. The increase in crystallinity was linear with time, while the ester bond degradation showed an inverse logarithmic correlation with time. All samples showed an increase in modulus with exposure time for the first loading cycle. Modulus changes for the consecutive loading cycles showed a nonlinear relationship to the exposure time that depended on membrane type and maximum strain. In addition, the recovered elastic range showed an expected increase with the maximum strain reached. The mechanical response of ES membranes was compared to experimental tensile properties of the human annulus fibrosus tissue and an in silico model of the intervertebral disk. The modulus of the tested membranes was at the lower range of the values found in literature, while the elastically recoverable strain after preconditioning for all membrane types lies within the desired strain range for this application.

CONCLUSION

The long-term assessment under application-specific conditions allowed to establish the mechanical competence of the electrospun PCL membranes. It can be concluded that with the use of appropriate fixation, the membranes can be used to create a seal on the damaged AF.

Validation of a finite element model of the thoracolumbar spine to study instrumentation level variations in early onset scoliosis correction

Growth-guidance constructs are an alternative to growing rods for the surgical treatment of early onset scoliosis (EOS). Constructs containing ultra-high molecular weight polyethylene (UHMWPE) sublaminar tape have been proposed as an improvement to the traditional Luque trolley. Ideally, a certain minimum number of levels is instrumented, thus offering the best balance between providing adequate spinal fixation and minimizing surgical exposure and spinal mobility reduction. The objective of the current study was to validate a parametric FE model of the thoracolumbar spine including its ability to predict the biomechanical effects of varying the number of levels instrumented with UHMWPE sublaminar tape in a growth-guidance construct for EOS correction. In a first step, the material properties of the L4-L5 segment in the model were calibrated relative to literature data. Next, whole thoracolumbar spine behavior was verified relative to literature data as well. Subsequently, rods, screws, and sublaminar tape were implemented in the model and a simulation of a previously performed in vitro experiment, in which the range of motion (ROM) of porcine spine segments was measured for different tape configurations, was performed. Good agreement between in vitro and FE-results was found for the changes in ROM before and after instrumentation. Good agreement for changes in ROM was obtained when varying the number of instrumented levels as well, indicating that the model can be a useful tool to evaluate the effects of construct composition variations. The present study was limited by the fact that only normal spine curvatures were analyzed and the fact that results of porcine spine experiments were compared to results of human FE models. Nevertheless, the good agreement in results, even at a detailed level, supports the idea that the model can ultimately be used as a pre-operative planning tool to evaluate different construct designs. The FE model of the thoracolumbar spine was successfully validated and was able to capture the biomechanical effect of construct component variations.

Stress distribution changes in growth plates of a trunk with adolescent idiopathic scoliosis following unilateral muscle paralysis: A hybrid musculoskeletal and finite element model

This study aimed to investigate changes occurred in the stress distribution in the growth plates (GPs) of a trunk with adolescent idiopathic scoliosis (AIS) following unilateral muscle paralysis. We hypothesized that weakening the appropriately chosen muscles on the concave side can decelerate AIS deformity progression. Muscle forces and reaction loads were estimated by an optimization-driven musculoskeletal (MS) model of adolescents with a normal- and an AIS trunk, and then applied on the finite element model of GPs of L1 through L4. Different set patterns of 95% reduction in the strength of the concave-side longissimus thoracis pars thoracic (LGPT), multifidus lumborum (MFL), and LGPT + MFL muscles were performed in the MS models. Results of this study showed that weakening of the concave-side MFL and LGPT muscles rendered a 35% correction in the symptomatic axial rotation of the AIS spine, and a reduction of about 25% in the compressive von Mises stress on the concave side of GPs, respectively, which can decelerate the deformity progression. It was observed that unilateral muscle weakening caused a compensatory activation of the rest of muscles to retain the spine stability. The intradiscal pressures and ratio between the rotations toward either side of the scoliotic spine, found here, matched well with some recent in-vivo investigations. One of the applications of the stability-based MS model of AIS spine with unilaterally weakened muscles presented in this study is to optimize the performance of the currently used braces. To fortify the presented therapeutic approach, experiments should be done on scoliotic animals.

Lumbar spinal ligament characteristics extracted from stepwise reduction experiments allow for preciser modeling than literature data

Lumbar ligaments play a key role in stabilizing the spine, particularly assisting muscles at wide-range movements. Hence, valid ligament force–strain data are required to generate physiological model predictions. These data have been obtained by experiments on single ligaments or functional units throughout the literature. However, contrary to detailed spine geometries, gained, for instance, from CT data, ligament characteristics are often inattentively transferred to multi-body system (MBS) or finite element models. In this paper, we use an elaborated MBS model of the lumbar spine to demonstrate how individualized ligament characteristics can be obtained by reversely reenacting stepwise reduction experiments, where the range of motion (ROM) was measured. We additionally validated the extracted characteristics with physiological experiments on intradiscal pressure (IDP). Our results on a total of in each case 160 ROM and 49 IDP simulations indicated superiority of our procedure (seven and eight outliers) toward the incorporation of classical literature data (on average 71 and 31 outliers).

Relative Nucleus Pulposus Area and Position Alters Disc Joint Mechanics

Aging and degeneration of the intervertebral disc are noted by changes in tissue composition and geometry, including a decrease in nucleus pulposus (NP) area. The NP centroid is positioned slightly posterior of the disc's centroid, but the effect of NP size and location on disc joint mechanics is not well understood. We evaluated the effect of NP size and centroid location on disc joint mechanics under dual-loading modalities (i.e., compression in combination with axial rotation or bending). A finite element model was developed to vary the relative NP area (NP:Disc area ratio range = 0.21 - 0.60). We also evaluated the effect of NP position by shifting the NP centroid anteriorly and posteriorly. Our results showed that compressive stiffness and average first principal strains increased with NP size. Under axial compression, stresses are distributed from the NP to the annulus, and stresses were redistributed towards the NP with axial rotation. Moreover, peak stresses were greater for discs with a smaller NP area. NP centroid location had a greater impact on intradiscal pressure during flexion and extension, where peak pressures in the posterior annulus under extension was greater for discs with a more posteriorly situated NP. In conclusion, the findings from this study highlight the importance of closely mimicking NP size and location in computational models that aim to understand stress/strain distribution during complex loading and for developing repair strategies that aim to recapitulate the mechanical behavior of healthy discs.

The rib cage stiffens the thoracic spine in a cadaveric model with body weight load under dynamic moments

The thoracic spine presents a challenge for biomechanical testing. With more segments than the lumbar and cervical regions and the integration with the rib cage, experimental approaches to evaluate the mechanical behavior of cadaveric thoracic spines have varied widely. Some researchers are now including the rib cage intact during testing, and some are incorporating follower load techniques in the thoracic spine. Both of these approaches aim to more closely model physiological conditions. To date, no studies have examined the impact of the rib cage on thoracic spine motion and stiffness in conjunction with follower loads. The purpose of this research was to quantify the mechanical effect of the rib cage on cadaveric thoracic spine motion and stiffness with a follower load under dynamic moments. It was hypothesized that the rib cage would increase stiffness and decrease motion of the thoracic spine with a follower load. Eight fresh-frozen human cadaveric thoracic spines with rib cages (T1-T12) were loaded with a 400 N compressive follower load. Dynamic moments of ± 5 N m were applied in lateral bending, flexion/extension, and axial rotation, and the motion and stiffness of the specimens with the rib cage intact have been previously reported. This study evaluated the motion and stiffness of the specimens after rib cage removal, and compared the data to the rib cage intact condition. Range-of-motion and stiffness were calculated for the upper, middle, and lower segments of the thoracic spine. Range-of-motion significantly increased with the removal of the rib cage in lateral bending, flexion/extension, and axial rotation by 63.5%, 63.0%, and 58.8%, respectively (p < 0.05). Neutral and elastic zones increased in flexion/extension and axial rotation, and neutral zone stiffness decreased in axial rotation with rib cage removal. Overall, the removal of the rib cage increases the range-of-motion and decreases the stiffness of cadaveric thoracic spines under compressive follower loads in vitro. This study suggests that the rib cage should be included when testing a cadaveric thoracic spine with a follower load to optimize clinical relevance.

Effects of Different Angles of the Traction Table on Lumbar Spine Ligaments: A Finite Element Study

Background

The traction bed is a noninvasive device for treating lower back pain caused by herniated intervertebral discs. In this study, we investigated the impact of the traction bed on the lower back as a means of increasing the disc height and creating a gap between facet joints.

Methods

Computed tomography (CT) images were obtained from a female volunteer and a three-dimensional (3D) model was created using software package MIMICs 17.0. Afterwards, the 3D model was analyzed in an analytical software (Abaqus 6.14). The study was conducted under the following traction loads: 25%, 45%, 55%, and 85% of the whole body weight in different angles.

Results

Results indicated that the loading angle in the L3–4 area had 36.8%, 57.4%, 55.32%, 49.8%, and 52.15% effect on the anterior longitudinal ligament, posterior longitudinal ligament, intertransverse ligament, interspinous ligament, and supraspinous ligament, respectively. The respective values for the L4–5 area were 32.3%, 10.6%, 53.4%, 56.58%, and 57.35%. Also, the body weight had 63.2%, 42.6%, 44.68%, 50.2%, and 47.85% effect on the anterior longitudinal ligament, posterior longitudinal ligament, intertransverse ligament, interspinous ligament, and supraspinous ligament, respectively. The respective values for the L4–5 area were 67.7%, 89.4%, 46.6%, 43.42% and 42.65%. The authenticity of results was checked by comparing with the experimental data.

Conclusions

The results show that traction beds are highly effective for disc movement and lower back pain relief. Also, an optimal angle for traction can be obtained in a 3D model analysis using CT or magnetic resonance imaging images. The optimal angle would be different for different patients and thus should be determined based on the decreased height of the intervertebral disc, weight and height of patients.

Osseointegration improves bone-implant interface of pedicle screws in the growing spine: a biomechanical and histological study using an in vivo immature porcine model

PURPOSE

Implant failure is a frequent complication in corrective surgery for early onset scoliosis, since considerable forces are acting on small and fragile vertebrae. Osseointegration showing biomechanical and histological improvement in bone-implant interface (BII) after dental implant placement has been well investigated. However, there are no studies regarding osseointegration in immature vertebral bone. The purpose was to evaluate the timecourse of biomechanical and histological changes at BII after pedicle screw placement using in vivo immature porcine model.

METHODS

Ten immature porcine were instrumented with titanium pedicle screws in the thoracic spine. After a 0-, 2-, 4-, and 6-month survival periods, the spines were harvested at the age of 12 months. Histological evaluation of BII was conducted by bone volume/tissue volume (BV/TV) and bone surface/implant surface (BS/IS) measurements. Bone mineral density (BMD) measurement and biomechanical testing of BII were done.

RESULTS

Contact surface and bone volume around the screw threads were significantly increased over the time. BV/TV and BS/IS were improved with statistically significant differences between 0- and ≥4-month (p ≤ 0.001) periods. BMD in all subjects was determined to be the same (p ≥ 0.350). Pullout strength was also increased over time with significant differences between 0- and ≥2-month (p ≤ 0.011) periods.

CONCLUSION

Improved stability at BII caused by osseointegration was confirmed by in vivo immature porcine model. A two-stage operation is proposed based on the osseointegration theory, in which an implant is installed in advance in the vertebrae at the first stage and deformity correction surgery is performed after sufficient stability is obtained by osseointegration at a later stage.

The Effectiveness of Percutaneous Vertebroplasty Is Determined by the Patient-Specific Bone Condition and the Treatment Strategy

Purpose

Vertebral fragility fractures are often treated by injecting bone cement into the collapsed vertebral bodies (vertebroplasty). The mechanisms by which vertebroplasty induces pain relief are not completely understood yet and recent debates cast doubt over the outcome of the procedure. The controversy is intensified by inconsistent results of randomized clinical trials and biomechanical studies that have investigated the effectiveness or the change in biomechanical response due to the reinforcement. The purpose of this study was to evaluate the effectiveness of vertebroplasty, by varying the relevant treatment parameters and (a) computationally predicting the improvement of the fracture risk depending on the chosen treatment strategy, and (b) identifying the determinants of a successful treatment.

Methods

A Finite Element model with a patient-specific failure criterion and direct simulation of PMMA infiltration in four lumbar vertebrae was used to assess the condition of the bone under compressive load before and after the virtual treatment, simulating in a total of 12000 virtual treatments.

Results

The results showed that vertebroplasty is capable of reducing the fracture risk by magnitudes, but can also have a detrimental effect. Effectiveness was strongly influenced by interactions between local bone quality, cement volume and injection location. However, only a moderate number of the investigated treatment strategies were able to achieve the necessary improvement for preventing a fracture.

Conclusions

We conclude that the effectiveness of vertebroplasty is sensitive to the patient’s condition and the treatment strategy.

A database of lumbar spinal mechanical behavior for validation of spinal analytical models

Data from two experimental studies with eight specimens each of spinal motion segments and/or intervertebral discs are presented in a form that can be used for comparison with finite element model predictions. The data include the effect of compressive preload (0, 250 and 500N) with quasistatic cyclic loading (0.0115Hz) and the effect of loading frequency (1, 0.1, 0.01 and 0.001Hz) with a physiological compressive preload (mean 642N). Specimens were tested with displacements in each of six degrees of freedom (three translations and three rotations) about defined anatomical axes. The three forces and three moments in the corresponding axis system were recorded during each test. Linearized stiffness matrices were calculated that could be used in multi-segmental biomechanical models of the spine and these matrices were analyzed to determine whether off-diagonal terms and symmetry assumptions should be included. These databases of lumbar spinal mechanical behavior under physiological conditions quantify behaviors that should be present in finite element model simulations. The addition of more specimens to identify sources of variability associated with physical dimensions, degeneration, and other variables would be beneficial. Supplementary data provide the recorded data and Matlab® codes for reading files. Linearized stiffness matrices derived from the tests at different preloads revealed few significant unexpected off-diagonal terms and little evidence of significant matrix asymmetry.

Hypothesis about an existent biomechanical cause-effect relationship between Schëuermann's kyphosis and scoliosis

Schëuermanńs kyphosis is usually observed with a mild idiopathic scoliosis, and there is parity between these two diseases. The aim of this work is to establish a hypothesis about the existence of a biomechanical causal relationship between Schëuermann's kyphosis and scoliosis. To achieve this, a literature review was conducted. A simple mechanical model of the passive thoracolumbar subsystem was created to support part of the discussion. This mechanical model describes the passive thoracolumbar subsystem under ideal conditions of equilibrium. After giving consideration to the system under these conditions, some of the geometrical changes that may be found in Schëuermanńs kyphosis are considered. Next, this work discusses the evolution of the spine, taking into account its relationship with stable equilibrium, which the passive subsystem tends toward. We hypothesized about the postural response of the body to compensate for possible situations of imbalance. In conclusion, it can be found that a change in the alignment of the spine may occur due to the postural adaptation of the body to an inadequate mechanical situation that may lead to scoliotic deformity of the spine.

Mechanism of right thoracic adolescent idiopathic scoliosis at risk for progression; a unifying pathway of development by normal growth and imbalance

Adolescent idiopathic scoliosis is regarded as a multifactorial disease and none of the many suggested causal etiologies have yet prevailed. I will suggest that adolescent idiopathic scoliosis has one common denominator, namely that initial curve development is mediated through one common normal physiological pathway of thoracic rotational instability. This is a consequence of gender specific natural growth of the passive structural components of thoracic spinal tissues for the adolescent female. This causes an unbalanced mechanical situation, which progresses if the paravertebral muscles cannot maintain spinal alignment. The alteration in the coronal plane with the lateral curve deformity is an uncoupling effect due to a culmination of a secondary, temporary sagittal plane thoracic flattening and of a primary, temporary transverse plane rotational instability for the adolescent female. Treatment of adolescent idiopathic scoliosis should address this physiological pathway and the overall treatment strategy is early intervention with strengthening of thoracic rotational stability for small curve adolescent idiopathic scoliosis.

Human L3L4 intervertebral disc mean 3D shape, modes of variation, and their relationship to degeneration

Intervertebral disc mechanics are affected by both disc shape and disc degeneration, which in turn each affect the other; disc mechanics additionally have a role in the etiology of disc degeneration. Finite element analysis (FEA) is a favored tool to investigate these relationships, but limited data for intervertebral disc 3D shape has forced the use of simplified or single-subject geometries, with the effect of inter-individual shape variation investigated only in specialized studies. Similarly, most data on disc shape variation with degeneration is based on 2D mid-sagittal images, which incompletely define 3D shape changes. Therefore, the objective of this study was to quantify inter-individual disc shape variation in 3D, classify this variation into independently-occurring modes using a statistical shape model, and identify correlations between disc shape and degeneration. Three-dimensional disc shapes were obtained from MRI of 13 human male cadaver L3L4 discs. An average disc shape and four major modes of shape variation (representing 90% of the variance) were identified. The first mode represented disc axial area and was significantly correlated to degeneration (R(2)=0.44), indicating larger axial area in degenerate discs. Disc height variation occurred in three distinct modes, each also involving non-height variation. The statistical shape model provides an average L3L4 disc shape for FEA that is fully defined in 3D, and makes it convenient to generate a set of shapes with which to represent aggregate inter-individual variation. Degeneration grade-specific shapes can also be generated. To facilitate application, the model is included in this paper׳s supplemental content.

The Use of Finite Element Models to Assist Understanding and Treatment For Scoliosis: A Review Paper

INTRODUCTION

Scoliosis is a complex spinal deformity whose etiology is still unknown, and its treatment presents many challenges. Finite element modeling (FEM) is one of the analytical techniques that has been used to elucidate the mechanism of scoliosis and the effects of various treatments.

METHODS

A literature review on the application of FEM in scoliosis evaluation and treatment has been undertaken. A literature search was performed in each of three major electronic databases (Google Scholar, Web of Science, and Ovid) using the key words "scoliosis" and "finite element methods/model". Articles using FEM and having a potential impact on clinical practice were included.

RESULTS

A total of 132 abstracts were retrieved. The query returned 105 articles in which the abstracts appeared to correspond to this review's focus, and 85 papers were retained. The current state of the art of FEM related to the biomechanical analysis of scoliosis is discussed in 4 sections: the etiology of adolescent idiopathic scoliosis, brace treatment, instrumentation treatment, and sensitivity studies of FEM. The limitations of FEM and suggested future work are also discussed.

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

STUDY DESIGN

Finite element study.

OBJECTIVE

To analyze the effects of posterior shear loads, disc degeneration, and the combination of both on spinal torsion stiffness.

SUMMARY OF BACKGROUND DATA

Scoliosis is a 3-dimensional deformity of the spine that presents itself mainly in adolescent girls and elderly patients. Our concept of its etiopathogenesis is that an excess of posteriorly directed shear loads, relative to the body's intrinsic stabilizing mechanisms, induces a torsional instability of the spine, making it vulnerable to scoliosis. Our hypothesis for the elderly spine is that disc degeneration compromises the stabilizing mechanisms.

METHODS

In an adult lumbar motion segment model, the disc properties were varied to simulate different aspects of disc degeneration. These models were then loaded with a pure torsion moment in combination with either a shear load in posterior direction, no shear, or a shear load in anterior direction.

RESULTS

Posteriorly directed shear loads reduced torsion stiffness, anteriorly directed shear loads increased torsion stiffness. These effects were mainly caused by a later (respectively earlier) onset of facet joint contact. Disc degeneration cases with a decreased disc height that leads to slackness of the annular fibers and ligaments caused a significantly decreased torsional stiffness. The combination of this stage with posterior shear loading reduced the torsion stiffness to less than half the stiffness of a healthy disc without shear loads. The end stage of disc degeneration increased torsion stiffness again.

CONCLUSION

The combination of a decreased disc height, that leads to slack annular fibers and ligaments, and posterior shear loads very significantly affects torsional stiffness: reduced to less than half the stiffness of a healthy disc without shear loads. Disc degeneration, thus, indeed compromises the stabilizing mechanisms of the elderly spine. A combination with posteriorly directed shear loads could then make it vulnerable to scoliosis.

LEVEL OF EVIDENCE

N/A.

Vertebral morphology influences the development of Schmorl's nodes in the lower thoracic vertebrae

Schmorl's nodes are the result of herniations of the nucleus pulposus into the adjacent vertebral body and are commonly identified in both clinical and archaeological contexts. The current study aims to identify aspects of vertebral shape that correlate with Schmorl's nodes. Two-dimensional statistical shape analysis was performed on digital images of the lower thoracic spine (T10-T12) of adult skeletons from the late medieval skeletal assemblages from Fishergate House, York, St. Mary Graces and East Smithfield Black Death cemeteries, London, and postmedieval Chelsea Old Church, London. Schmorl's nodes were scored on the basis of their location, depth, and size. Results indicate that there is a correlation between the shape of the posterior margin of the vertebral body and pedicles and the presence of Schmorl's nodes in the lower thoracic spine. The size of the vertebral body in males was also found to correlate with the lesions. Vertebral shape differences associated with the macroscopic characteristics of Schmorl's nodes, indicating severity of the lesion, were also analyzed. The shape of the pedicles and the posterior margin of the vertebral body, along with a larger vertebral body size in males, have a strong association with both the presence and severity of Schmorl's nodes. This suggests that shape and/or size of these vertebral components are predisposing to, or resulting in, vertically directed disc herniation.

Geometry strongly influences the response of numerical models of the lumbar spine--a probabilistic finite element analysis

Typical FE models of the human lumbar spine consider a single, fixed geometry. Such models cannot account for potential effects of the natural variability of the spine's geometry. In this study, we performed a probabilistic uncertainty and sensitivity analysis of a fully parameterized, geometrically simplified model of the L3-L4 segment. We examined the impact of the uncertainty in all 40 geometry parameters, estimated lower and upper bounds for the required sample size and determined the most important geometry parameters. The natural variability of the spine's geometry indeed strongly affects intradiscal pressure, range of motion and facet joint contact forces. Deriving generalized statements from fixed-geometry models as well as transferring those results to different cases thus can easily lead to wrong conclusions and should only be performed with extreme caution. We recommend a sample size of ≈ 100 to obtain reasonable accurate point estimates and a sufficient overview of the remaining uncertainties. Yet, only few parameters, especially those determining the disc geometry (disc height, end-plate width and depth) and the facets' position (intra-articular space, pedicle length, facet angles), proved to be truly important. Accurate measurement and modeling of those structures should therefore be prioritized.

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

STUDY DESIGN

A validated finite element model of an L3-L4 motion segment is used to analyze the effects of interpersonal differences in geometry on spinal stiffness.

OBJECTIVE

The objective of this study is to determine which of the interpersonal variations of the geometry of the spine have a large effect on spinal stiffness. This will improve patient-specific modeling.

SUMMARY OF BACKGROUND DATA

The parameters that define the geometry of a motion segment are vertebral height, disc height, endplate width, endplate depth, spinous process length, transverse process width, nucleus size, lordosis angle, facet area, facet orientation, and the cross-sectional areas of the ligaments. All these parameters differ between patients. The influence of each parameter on spinal stiffness is largely unknown and such knowledge would greatly help in patient-specific modeling of the spine.

METHODS

The range of interpersonal variation of each of the geometric parameters was set at mean±2SD (covering 95% of the population). Subsequently, we determined the effect of each of these ranges on the bending stiffness in flexion, extension, axial rotation, and lateral bending.

RESULTS

Disc height had the largest influence; a maximal disc height reduced the spinal stiffness to 75-86% of the mean motion segment stiffness, and a minimal disc height increased the spinal stiffness to 154-226% of the mean motion segment stiffness. Lordosis angle, transversal and longitudinal facet angle, endplate depth, and area of the capsular ligament also had a substantial influence (>5%) on the stiffness, but considerable less than the influence of the disc height. Ligament areas, nucleus size, spinous process length, and length of processes are of negligible effect (<2%) on the stiffness.

CONCLUSION

The disc height should be accurately determined in patients to estimate the spinal stiffness. Ligament areas, nucleus size, spinous process length, and transverse process width do not need patient-specific modeling.