Three-dimensional finite element model of the cervical spinal cord: preliminary results of injury mechanism analysis

Finite element modeling of the human cervical spinal cord and its applications: A systematic review

Background Context

Finite element modeling (FEM) is an established tool to analyze the biomechanics of complex systems. Advances in computational techniques have led to the increasing use of spinal cord FEMs to study cervical spinal cord pathology. There is considerable variability in the creation of cervical spinal cord FEMs and to date there has been no systematic review of the technique. The aim of this study was to review the uses, techniques, limitations, and applications of FEMs of the human cervical spinal cord.

Methods

A literature search was performed through PubMed and Scopus using the words finite element analysis, spinal cord, and biomechanics. Studies were selected based on the following inclusion criteria: (1) use of human spinal cord modeling at the cervical level; (2) model the cervical spinal cord with or without the osteoligamentous spine; and (3) the study should describe an application of the spinal cord FEM.

Results

Our search resulted in 369 total publications, 49 underwent reviews of the abstract and full text, and 23 were included in the study. Spinal cord FEMs are used to study spinal cord injury and trauma, pathologic processes, and spine surgery. Considerable variation exists in the derivation of spinal cord geometries, mathematical models, and material properties. Less than 50% of the FEMs incorporate the dura mater, cerebrospinal fluid, nerve roots, and denticulate ligaments. Von Mises stress, and strain of the spinal cord are the most common outputs studied. FEM offers the opportunity for dynamic simulation, but this has been used in only four studies.

Conclusions

Spinal cord FEM provides unique insight into the stress and strain of the cervical spinal cord in various pathological conditions and allows for the simulation of surgical procedures. Standardization of modeling parameters, anatomical structures and inclusion of patient-specific data are necessary to improve the clinical translation.

Effect of Cervical Stenosis and Rate of Impact on Risk of Spinal Cord Injury During Whiplash Injury

STUDY DESIGN

Finite Element Study.

OBJECTIVE

To determine the risk of spinal cord injury with pre-existing cervical stenosis during a whiplash injury.

SUMMARY OF BACKGROUND DATA

Patients with cervical spinal stenosis are often cautioned on the potential increased risk of spinal cord injury (SCI) from minor trauma such as rear impact whiplash injuries. However, there is no consensus on the degree of canal stenosis or the rate of impact that predisposes cervical SCI from minor trauma.

METHODS

A previously validated three-dimensional finite element model of the human head-neck complex with the spinal cord and activated cervical musculature was used. Rear impact acceleration was applied at 1.8 m/s and 2.6 m/s. Progressive spinal stenosis was simulated at the C5 to C6 segment, from 14 mm to 6 mm, at 2 mm intervals of ventral disk protrusion. Spinal cord von Mises stress and maximum principal strain were extracted and normalized with respect to the 14 mm spine at each cervical spine level from C2 to C7.

RESULTS

The mean segmental range of motion was 7.3 degrees at 1.8 m/s and 9.3 degrees at 2.6 m/s. Spinal cord stress above the threshold for SCI was noted at C5 to C6 for 6 mm stenosis at 1.8 m/s and 2.6 m/s. The segment (C6-C7) inferior to the level of maximum stenosis also showed increasing stress and strain with a higher rate of impact. For 8 mm stenosis, spinal cord stress exceeded SCI thresholds only at 2.6 m/s. Spinal cord strain above SCI thresholds were only noted in the 6 mm stenosis model at 2.6 m/s.

CONCLUSION

Increased spinal stenosis and rate of impact are associated with greater magnitude and spatial distribution of spinal cord stress and strain during a whiplash injury. Spinal canal stenosis of 6 mm was associated with consistent elevation of spinal cord stress and strain above SCI thresholds at 2.6 m/s.

Identification of injury type using somatosensory and motor evoked potentials in a rat spinal cord injury model

The spinal cord is at risk of injury during spinal surgery. If intraoperative spinal cord injury is identified early, irreversible impairment or loss of neurological function can be prevented. Different types of spinal cord injury result in damage to different spinal cord regions, which may cause different somatosensory and motor evoked potential signal responses. In this study, we examined electrophysiological and histopathological changes between contusion, distraction, and dislocation spinal cord injuries in a rat model. We found that contusion led to the most severe dorsal white matter injury and caused considerable attenuation of both somatosensory and motor evoked potentials. Dislocation resulted in loss of myelinated axons in the lateral region of the injured spinal cord along the rostrocaudal axis. The amplitude of attenuation in motor evoked potential responses caused by dislocation was greater than that caused by contusion. After distraction injury, extracellular spaces were slightly but not significantly enlarged; somatosensory evoked potential responses slightly decreased and motor evoked potential responses were lost. Correlation analysis showed that histological and electrophysiological findings were significantly correlated and related to injury type. Intraoperative monitoring of both somatosensory and motor evoked potentials has the potential to identify iatrogenic spinal cord injury type during surgery.

Effects of cervical rotatory manipulation on the cervical spinal cord complex with ossification of the posterior longitudinal ligament in the vertebral canal: A finite element study

Background: There are few studies focusing on biomechanism of spinal cord injury according to the ossification of the posterior longitudinal ligament (OPLL) during cervical rotatory manipulation (CRM). This study aimed to explore the biomechanical effects of CRM on the spinal cord, dura matter and nerve roots with OPLL in the cervical vertebral canal. Methods: Three validated FE models of the craniocervical spine and spinal cord complex were constructed by adding mild, moderate, and severe OPLL to the healthy FE model, respectively. We simulated the static compression of the spinal cord by OPLL and the dynamic compression during CRM in the flexion position. The stress distribution of the spinal cord complex was investigated. Results: The cervical spinal cord experienced higher von Mises stress under static compression by the severe OPLL. A higher von Mises stress was observed on the spinal cord in the moderate and severe OPLL models during CRM. The dura matter and nerve roots had a higher von Mises stress in all three models during CRM. Conclusion: The results show a high risk in performing CRM in the flexion position on patients with OPLL, in that different occupying ratios in the vertebral canal due to OPLL could significantly increase the stress on the spinal cord complex.

Effect of Velocity and Contact Stress Area on the Dynamic Behavior of the Spinal Cord Under Different Testing Conditions

Knowledge of the dynamic behavior of the spinal cord under different testing conditions is critical for our understanding of biomechanical mechanisms of spinal cord injury. Although velocity and contact stress area are known to affect external mechanical stress or energy upon sudden traumatic injury, quantitative investigation of the two clinically relevant biomechanical variables is limited. Here, freshly excised rat spinal-cord–pia-arachnoid constructs were tested through indentation using indenters of different sizes (radii: 0.25, 0.50, and 1.00 mm) at various loading rates ranging from 0.04 to 0.20 mm/s. This analysis found that the ex vivo specimen displayed significant nonlinear viscoelasticity at <10% of specimen thickness depth magnitudes. At higher velocity and larger contact stress area, the cord withstood a higher peak load and exhibited more sensitive mechanical relaxation responses (i.e., increasing amplitude and speed of the drop in peak load). Additionally, the cord became stiffer (i.e., increasing elastic modulus) and softer (i.e., decreasing elastic modulus) at a higher velocity and larger contact stress area, respectively. These findings will improve our understanding of the real-time complex biomechanics involved in traumatic spinal cord injury.

A Hyper-Viscoelastic Continuum-Level Finite Element Model of the Spinal Cord Assessed for Transverse Indentation and Impact Loading

Finite Element (FE) modelling of spinal cord response to impact can provide unique insights into the neural tissue response and injury risk potential. Yet, contemporary human body models (HBMs) used to examine injury risk and prevention across a wide range of impact scenarios often lack detailed integration of the spinal cord and surrounding tissues. The integration of a spinal cord in contemporary HBMs has been limited by the need for a continuum-level model owing to the relatively large element size required to be compatible with HBM, and the requirement for model development based on published material properties and validation using relevant non-linear material data. The goals of this study were to develop and assess non-linear material model parameters for the spinal cord parenchyma and pia mater, and incorporate these models into a continuum-level model of the spinal cord with a mesh size conducive to integration in HBM. First, hyper-viscoelastic material properties based on tissue-level mechanical test data for the spinal cord and hyperelastic material properties for the pia mater were determined. Secondly, the constitutive models were integrated in a spinal cord segment FE model validated against independent experimental data representing transverse compression of the spinal cord-pia mater complex (SCP) under quasi-static indentation and dynamic impact loading. The constitutive model parameters were fit to a quasi-linear viscoelastic model with an Ogden hyperelastic function, and then verified using single element test cases corresponding to the experimental strain rates for the spinal cord (0.32–77.22 s⁻¹) and pia mater (0.05 s⁻¹). Validation of the spinal cord model was then performed by re-creating, in an explicit FE code, two independent ex-vivo experimental setups: 1) transverse indentation of a porcine spinal cord-pia mater complex and 2) dynamic transverse impact of a bovine SCP. The indentation model accurately matched the experimental results up to 60% compression of the SCP, while the impact model predicted the loading phase and the maximum deformation (within 7%) of the SCP experimental data. This study quantified the important biomechanical contribution of the pia mater tissue during spinal cord deformation. The validated material models established in this study can be implemented in computational HBM.

The biomechanical effect of preexisting different types of disc herniation in cervical hyperextension injury

OBJECTIVE

Preexisting severe cervical spinal cord compression is a significant risk factor in cervical hyperextension injury, and the neurological function may deteriorate after a slight force to the forehead. There are few biomechanical studies regarding the influence of pathological factors in hyperextension loading condition. The aim of this study is to analyze the effects of preexisting different types of cervical disc herniation and different degrees of compression on the spinal cord in cervical hyperextension.

METHOD

A 3D finite element (FE) model of cervical spinal cord was modeled. Local type with median herniation, local type with lateral herniation, diffuse type with median herniation, and diffuse type with lateral herniation were simulated in neutral and extention positions. The compressions which were equivalent to 10%, 20%, 30%, and 40% of the sagittal diameter of the spinal cord were modeled.

RESULTS

The results of normal FE model were consistent with those of previous studies. The maximum von Mises stresses appeared in the pia mater for all 32 loading conditions. The maximum von Mises stresses in extension position were much higher than in neutral position. In most cases, the maximum von Mises stresses in diffuse type were higher than in local type.

CONCLUSION

Cervical spinal cord with preexisting disc herniation is more likely to be compressed in hyperextension situation than in neutral position. Diffuse type with median herniation may cause more severe compression with higher von Mises stresses concentrated at the anterior horn and the peripheral white matter, resulting in acute central cord syndrome from biomechanical point of view.

The corticospinal tract structure of collagen/silk fibroin scaffold implants using 3D printing promotes functional recovery after complete spinal cord transection in rats

No effective treatment has been established for nerve dysfunction caused by spinal cord injury (SCI). Orderly axonal growth at the site of spinal cord transection and creation of an appropriate biological microenvironment are important for functional recovery. To axially guiding axonal growth, designing a collagen/silk fibroin scaffold fabricated with 3D printing technology (3D-C/SF) emulated the corticospinal tract. The normal collagen/silk fibroin scaffold with freeze-drying technology (C/SF) or 3D-C/SF scaffold were implanted into rats with completely transected SCI to evaluate its effect on nerve repair during an 8-week observation period. Electrophysiological analysis and locomotor performance showed that the 3D-C/SF implants contributed to significant improvements in the neurogolical function of rats compared to C/SF group. By magnetic resonance imaging, 3D-C/SF implants promoted a striking degree of axonal regeneration and connection between the proximal and distal SCI sites. Compared with C/SF group, rats with 3D-C/SF scaffold exhibited fewer lesions and disordered structures in histological analysis and more GAP43-positive profiles at the lesion site. The above results indicated that the corticospinal tract structure of 3D printing collagen/silk fibroin scaffold improved axonal regeneration and promoted orderly connections within the neural network, which could provided a promising and innovative approach for tissue repair after SCI.

Correlating Tissue Mechanics and Spinal Cord Injury: Patient-Specific Finite Element Models of Unilateral Cervical Contusion Spinal Cord Injury in Non-Human Primates

Non-human primate (NHP) models are the closest approximation of human spinal cord injury (SCI) available for pre-clinical trials. The NHP models, however, include broader morphological variability that can confound experimental outcomes. We developed subject-specific finite element (FE) models to quantify the relationship between impact mechanics and SCI, including the correlations between FE outcomes and tissue damage. Subject-specific models of cervical unilateral contusion SCI were generated from pre-injury MRIs of six NHPs. Stress and strain outcomes were compared with lesion histology using logit analysis. A parallel generic model was constructed to compare the outcomes of subject-specific and generic models. The FE outcomes were correlated more strongly with gray matter damage (0.29 < R2 < 0.76) than white matter (0.18 < R2 < 0.58). Maximum/minimum principal strain, Von-Mises and Tresca stresses showed the strongest correlations (0.31 < R2 < 0.76) with tissue damage in the gray matter while minimum principal strain, Von-Mises stress, and Tresca stress best predicted white matter damage (0.23 < R2 < 0.58). Tissue damage thresholds varied for each subject. The generic FE model captured the impact biomechanics in two of the four models; however, the correlations between FE outcomes and tissue damage were weaker than the subject-specific models (gray matter [0.25 < R2 < 0.69] and white matter [R2 < 0.06] except for one subject [0.26 < R2 < 0.48]). The FE mechanical outputs correlated with tissue damage in spinal cord white and gray matters, and the subject-specific models accurately mimicked the biomechanics of NHP cervical contusion impacts.

Effects of a contusion load on spinal cord with different curvatures

The cervical spine injury is a complicated procedure in the combination of different injury loads and postures. The aim of this study is to investigate the injury mechanism considering different types of cervical curvatures subjected to contusion loads. A finite element model of a cervical spinal cord was constructed. Gray matter, white matter and pia matter were modeled and hyperelastic material properties were assigned. Convergence analysis and validation analysis were carried out. The model was simulated in 3 different spinal curvatures and loaded by 2 directions with 4 compression degrees. The maximum von Mises stress in the whole model was concentrated in the pia matter in all loading cases. When investigating spinal cord injury, the pia matter must be considered. For all three curvatures, the stress in the gray matter and white matter was higher in front-to-back loading condition than that in back-to-front loading condition. The front-to-back impact may cause a larger damage. A back-to-front load damaged the structure around the central canal and a front-to-back contusion load damaged the anterior horn of the spinal cord at most time. From the view of the maximum stress, the lordotic curvature did not show significant buffering effect. However, the pathological curvature had large areas affected and the lordotic curvature showed some benefits to some degree from the view of stress distribution.

Accurate simulation of the herniated cervical intervertebral disc using controllable expansion: a finite element study

Expansions were carried out in finite element (FE) models of disc hernia including symmetric (median, lateral, paramedian) and asymmetric types. In all models, lubricous disk bulging that applied a linear compression to the anterior part of the cord was observed at the posterior surfaces of the expansion zone, respectively. The shape and position of protrusions varyed with the temperature, magnitude, and location of expanding elements. The geometric deformation and stress distribution of the spinal cord increased as the extent of compression grew. This method is in possession of enormous potential in promoting further individualized research of cervical spondylotic myelopathy.

Biomechanics of the effect of subaxial cervical spine degeneration on atlantoaxial complex in idiopathic retro-odontoid pseudotumor development

OBJECTIVES

Retro-odontoid pseudotumor (ROP), with no rheumatoid arthritis, atlantoaxial instability, or other primary diseases, is defined as idiopathic retro-odontoid pseudotumor (IROP). Cervical spine degeneration is associated with IROP development. This study aims to evaluate the effect of cervical spine degeneration on the atlantoaxial complex and find the possible biomechanical mechanism of IROP development.

METHODS

This study was performed using a three-dimensional (3D) finite element (FE) analysis. A degenerated FE model (FEM) and five operation FEMs (C1-C2 fusion, C0-C2 fusion, C0-C3 fusion, C0-C4 fusion, and C1 posterior arch resection) were established based on a normal 3D FEM of the cervical spine including C0-T1 with the main ligaments and muscles. The parameters, including the C1-C2 range of motions (ROMs) and odontoid-related ligaments' stresses in degenerated and operation FEMs, were obtained and compared with those in normal FEM.

RESULTS

Compared to normal FEM, degenerated FEM had reduced C3-C7 ROMs and increased C1-C2 ROMs and odontoid-related ligaments' stresses. After internal fixation, C1-C2 ROMs and most odontoid-related ligaments' stresses were greatly decreased, but with no significant differences among C0-C2, C0-C3, C0-C4, and C1-C2 fusion models. For the C1 posterior arch resection model, C1-C2 ROMs and most odontoid-related ligaments' stresses increased, compared with normal FEM.

CONCLUSIONS

Cervical spine degeneration plays an important part in IROP development in biomechanics. Atlantoaxial complex compensates for cervical spine degeneration, with increased C1-C2 ROMs and odontoid-related ligaments' stresses. Atlantoaxial fusion or short segmental occipitocervical fusion can effectively reduce the stress and should be considered in IROP treatment.

Adjacent-level biomechanics after single-level anterior cervical interbody fusion with anchored zero-profile spacer versus cage-plate construct: a finite element study

BACKGROUND

The development of adjacent segment degeneration (ASD) following ACDF is well established. There is no analytical study related to effects of plate profile on the biomechanics of the adjacent-level after ACDF. This study aimed to test the effects of plate profile on the adjacent-level biomechanics after single-level anterior cervical discectomy and fusion (ACDF).

METHODS

A three-dimensional finite element model (FEM) of an intact C2-T1 segment was built and validated. From this intact model, two instrumentation models were constructed with the anchored zero-profile spacer or the standard plate-interbody spacer after a C5-C6 corpectomy and fusion. Motion patterns, the stresses in the disc, the endplate, and the facet joint at the levels cephalad and caudal to the fusion were assessed.

RESULTS

Compared with the normal condition, the biomechanical responses in the adjacent levels were increased after fusion. Relative to the intact model, the average increase of range of motion (ROM) and stresses in the endplate, the disc, and the facet of the zero-profile spacer fusion model were slightly lower than that of the standard plate-interbody spacer fusion model. The kinematics ROM and stress variations above fusion segment were larger than that below. The biomechanical features of the adjacent segment after fusion were most affected during extension.

CONCLUSIONS

The FE analysis indicated that plate profile may have an impact on the biomechanics of the adjacent-level after a single-level ACDF. The impact may be long-term and cumulative. The current findings may help explain the decreasing incidence of ASD complications in the patients using zero-profile spacer compared with the patients using cage and plate construct.

Diffusion tensor imaging predicting neurological repair of spinal cord injury with transplanting collagen/chitosan scaffold binding bFGF

Prognosis and treatment evaluation of spinal cord injury (SCI) are still in the long-term research stage. Prognostic factors for SCI treatment need effective biomarker to assess therapeutic effect. Quantitative diffusion tensor imaging (DTI) may become a potential indicators for assessing SCI repair. However, its correlation with the results of locomotor function recovery and tissue repair has not been carefully studied. The aim of this study was to use quantitative DTI to predict neurological repair of SCI with transplanting collagen/chitosan scaffold binding basic fibroblast growth factor (bFGF). To achieve our research goals, T10 complete transection SCI model was established. Then collagen/chitosan mixture adsorbed with bFGF (CCS/bFGF) were implanted into rats with SCI. At 8 weeks after modeling, implanting CCS/bFGF demonstrated more significant improvements in locomotor function according to Basso-Beattie-Bresnahan (BBB) score, inclined-grid climbing test, and electrophysiological examinations. DTI was carried out to evaluate the repair of axons by diffusion tensor tractgraphy (DTT), fractional anisotropy (FA) and apparent diffusion coefficient (ADC), a numerical measure of relative white matter from the rostral to the caudal. Parallel to locomotor function recovery, the CCS/bFGF group could significantly promote the regeneration of nerve fibers tracts according to DTT, magnetic resonance imaging (MRI), Bielschowsky's silver staining and immunofluorescence staining. Positive correlations between imaging and locomotor function or histology were found at all locations from the rostral to the caudal (P < 0.0001). These results demonstrated that DTI might be used as an effective predictor for evaluating neurological repair after SCI in experimental trails and clinical cases.

Basic biomechanics of spinal cord injury - How injuries happen in people and how animal models have informed our understanding

The wide variability, or heterogeneity, in human spinal cord injury is due partially to biomechanical factors. This review summarizes our current knowledge surrounding the patterns of human spinal column injury and the biomechanical factors affecting injury. The biomechanics of human spinal injury is studied most frequently with human cadaveric models and the features of the two most common injury patterns, burst fracture and fracture dislocation, are outlined. The biology of spinal cord injury is typically studied with animal models and the effects of the most relevant biomechanical factors - injury mechanism, injury velocity, and residual compression, are described. Tissue damage patterns and behavioural outcomes following dislocation or distraction injury mechanisms differ from the more commonly used contusion mechanism. The velocity of injury affects spinal cord damage, principally in the white matter. Ongoing, or residual compression after the initial impact does affect spinal cord damage, but few models exist that replicate the clinical scenario. Future research should focus on the effects of these biomechanical factors in different preclinical animal models as recent data suggests that treatment outcomes may vary between models.

Engineering approaches to understanding mechanisms of spinal column injury leading to spinal cord injury

BACKGROUND

The mechanical interactions occurring between the spinal column and spinal cord during an injury event are complex and variable, and likely have implications for the clinical presentation and prognosis of the individual.

METHODS

The engineering approaches that have been developed to better understand spinal column and cord interactions during an injury event are discussed. These include injury models utilising human and animal cadaveric specimens, in vivo anaesthetised animals, finite element models, inanimate physical systems and combinations thereof.

FINDINGS

The paper describes the development of these modelling approaches, discusses the advantages and disadvantages of the various models, and the major outcomes that have had implications for spinal cord injury research and clinical practice.

INTERPRETATION

The contribution of these four engineering approaches to understanding the interaction between the biomechanics and biology of spinal cord injury is substantial; they have improved our understanding of the factors contributing to the spinal column disruption, the degree of spinal cord deformation or motion, and the resultant neurological deficit and imaging features. Models of the injury event are challenging to produce, but technological advances are likely to improve these models and, consequently, our understanding of the mechanical context in which the biological injury occurs.

The denticulate ligament - Tensile characterisation and finite element micro-scale model of the structure stabilising spinal cord

BACKGROUND

Damage to the spinal cord is one of the most debilitating pathologies with considerable health, economic and social impact. Improved prevention, treatment and rehabilitation after spinal cord injury (SCI) requires the complex biomechanics of the spinal cord with all its structural elements and the injury mechanism to be understood. This comprehensive understanding will also allow development of models and tools enabling better diagnosis, surgical treatment with increased safety and efficacy and possible development of regenerative therapies. The denticulate ligaments play an important role in stabilising spinal cord within the spinal canal. They participate in spinal cord movements and play a role in determining the stress distribution during physiological but also traumatic loading. We present detailed tensile characterisation of the denticulate ligaments and a Finite Element micro-scale model of the ligament relating its structure with the distribution of stress under physiological loading.

METHOD

Denticulate ligaments were dissected from cervical spinal levels from 6 porcine cervical specimens with fragments of the pia and dura mater and characterised in terms of their geometry and response to uniaxial tensile loading. The stress-strain characteristics were recorded until rupture of the ligament, ultimate parameters and Young's moduli were determined. The parametric micro-structural Finite Element model was constructed based on literature microscope and histological images of a denticulate ligament as a phenomenological representation of the complex microstructure of a soft tissue. The model was validated against the experimental data.

RESULTS

Stress-strain characteristics obtained in tensile test were typical for a soft tissue behaviour. No statistically relevant differences in ultimate strength, strain and Young's moduli were observed between the ligaments harvested from different vertebral levels. Average ultimate tensile stress was 1.26 ± 0.20 MPa at strain 0.51 ± 0.00, rupturing force (1.01 ± 0.21 N) was in agreement with results obtained previously. The Finite Element model accurately predicted the extension-load behaviour of the denticulate ligament in elastic regime. The micro-scale structural representation enabled capturing deformation modes representative of the experimentally observed behaviour.

CONCLUSIONS

The presented stress-strain characteristics of the denticulate ligaments add valuable data to the understanding of the biomechanics of the spinal cord and enable development of more accurate models. The developed micro-scale model was capable of capturing biomechanical response of collagenous tissue under tensile loading, it can be applied for the prediction of other soft tissues behaviours. The denticulate ligament model should be included into future spinal cord models to fully represent the complex system's biomechanics and enable development of surgical aid tools to improve patient outcomes and future regenerative therapies.

High-Speed Fluoroscopy to Measure Dynamic Spinal Cord Deformation in an In Vivo Rat Model

Although spinal cord deformation is thought to be a predictor of injury severity, few researchers have investigated dynamic cord deformation, in vivo, during impact. This is needed to establish correlations among impact parameters, internal cord deformation, and histological and functional outcomes. Relying on surface deformations alone may not sufficiently represent spinal cord deformation. The objective of this study was to develop a high-speed fluoroscopic method of tracking the surface and internal cord deformations of rat spinal cord during experimental cord injury. Two radio-opaque beads were injected into the cord at C5/6 in the dorsal and ventral white matter. Four additional beads were glued to the surface of the cord. Dynamic bead displacement was tracked during a dorsal impact (130 mm/sec, 1 mm depth) by high-speed radiographic imaging at 3000 FPS, laterally. The internal spinal cord beads displaced significantly more than the surface beads in the ventral direction (1.1-1.9 times) and more than most surface beads in the cranial direction (1.2-1.5 times). The dorsal beads (internal and surface) displaced more than the ventral beads during all impacts. The bead displacement pattern implies that the spinal cord undergoes complex internal and surface deformations during impact. Residual displacement of the internal beads was significantly greater than that of the surface beads in the cranial-caudal direction but not the dorsoventral direction. Finite element simulation confirmed that the additional bead mass likely had little effect on the internal cord deformations. These results support the merit of this technique for measuring in vivo spinal cord deformation.

Viscoelasticity of spinal cord and meningeal tissues

Compared to the outer dura mater, the mechanical behavior of spinal pia and arachnoid meningeal layers has received very little attention in the literature. This is despite experimental evidence of their importance with respect to the overall spinal cord stiffness and recovery following compression. Accordingly, inclusion of the mechanical contribution of the pia and arachnoid maters would improve the predictive accuracy of finite element models of the spine, especially in the distribution of stresses and strain through the cord's cross-section. However, to-date, only linearly elastic moduli for what has been previously identified as spinal pia mater is available in the literature. This study is the first to quantitatively compare the viscoelastic behavior of isolated spinal pia-arachnoid-complex, neural tissue of the spinal cord parenchyma, and intact construct of the two. The results show that while it only makes up 5.5% of the overall cross-sectional area, the thin membranes of the innermost meninges significantly affect both the elastic and viscous response of the intact construct. Without the contribution of the pia and arachnoid maters, the spinal cord has very little inherent stiffness and experiences significant relaxation when strained. The ability of the fitted non-linear viscoelastic material models of each condition to predict independent data within experimental variability supports their implementation into future finite element computational studies of the spine.

STATEMENT OF SIGNIFICANCE

The neural tissue of the spinal cord is surrounded by three fibrous layers called meninges which are important in the behavior of the overall spinal-cord-meningeal construct. While the mechanical properties of the outermost layer have been reported, the pia mater and arachnoid mater have received considerably less attention. This study is the first to directly compare the behavior of the isolated neural tissue of the cord, the isolated pia-arachnoid complex, and the construct of these individual components. The results show that, despite being very thin, the inner meninges significantly affect the elastic and time-dependent response of the spinal cord, which may have important implications for studies of spinal cord injury.

Comparison of in vivo and ex vivo viscoelastic behavior of the spinal cord

Despite efforts to simulate the in vivo environment, post-mortem degradation and lack of blood perfusion complicate the use of ex vivo derived material models in computational studies of spinal cord injury. In order to quantify the mechanical changes that manifest ex vivo, the viscoelastic behavior of in vivo and ex vivo porcine spinal cord samples were compared. Stress-relaxation data from each condition were fit to a non-linear viscoelastic model using a novel characterization technique called the direct fit method. To validate the presented material models, the parameters obtained for each condition were used to predict the respective dynamic cyclic response. Both ex vivo and in vivo samples displayed non-linear viscoelastic behavior with a significant increase in relaxation with applied strain. However, at all three strain magnitudes compared, ex vivo samples experienced a higher stress and greater relaxation than in vivo samples. Significant differences between model parameters also showed distinct relaxation behaviors, especially in non-linear relaxation modulus components associated with the short-term response (0.1-1 s). The results of this study underscore the necessity of utilizing material models developed from in vivo experimental data for studies of spinal cord injury, where the time-dependent properties are critical. The ability of each material model to accurately predict the dynamic cyclic response validates the presented methodology and supports the use of the in vivo model in future high-resolution finite element modeling efforts.

STATEMENT OF SIGNIFICANCE

Neural tissues (such as the brain and spinal cord) display time-dependent, or viscoelastic, mechanical behavior making it difficult to model how they respond to various loading conditions, including injury. Methods that aim to characterize the behavior of the spinal cord almost exclusively use ex vivo cadaveric or animal samples, despite evidence that time after death affects the behavior compared to that in a living animal (in vivo response). Therefore, this study directly compared the mechanical response of ex vivo and in vivo samples to quantify these differences for the first time. This will allow researchers to draw more accurate conclusions about spinal cord injuries based on ex vivo data (which are easier to obtain) and emphasizes the importance of future in vivo experimental animal work.

Analysis of stress application at the thoracolumbar junction and influence of vertebral body collapse on the spinal cord and cauda equina

The thoracolumbar junction comprises the spinal cord, nerve roots and the cauda equina, exhibiting unique anatomical features that may give rise to a diverse array of symptoms under conditions of injury, thus complicating the diagnosis of compressive disorders. The present study aimed to examine varying degrees and forms of compression at this level of the spinal cord using a two-dimensional model to calculate the relationship of these variables to injury. The degree of compression was expressed as a percentage of the spinal canal that was occupied. Results were compared with findings from clinical observations to assess the validity of the model. Analysis revealed that higher levels of compression/spinal canal occupation are associated with the presence of neurological symptoms. This finding was consistent with clinical data. Results of the present analysis warrant further research involving evaluation of compression with respect to other parameters, such as blood flow, as well as more anatomically accurate three-dimensional analysis.

Intubation biomechanics: validation of a finite element model of cervical spine motion during endotracheal intubation in intact and injured conditions

OBJECTIVE

Because of limitations inherent to cadaver models of endotracheal intubation, the authors’ group developed a finite element (FE) model of the human cervical spine and spinal cord. Their aims were to 1) compare FE model predictions of intervertebral motion during intubation with intervertebral motion measured in patients with intact cervical spines and in cadavers with spine injuries at C-2 and C3–4 and 2) estimate spinal cord strains during intubation under these conditions.

METHODS

The FE model was designed to replicate the properties of an intact (stable) spine in patients, C-2 injury (Type II odontoid fracture), and a severe C3–4 distractive-flexion injury from prior cadaver studies. The authors recorded the laryngoscope force values from 2 different laryngoscopes (Macintosh, high intubation force; Airtraq, low intubation force) used during the patient and cadaver intubation studies. FE-modeled motion was compared with experimentally measured motion, and corresponding cord strain values were calculated.

RESULTS

FE model predictions of intact intervertebral motions were comparable to motions measured in patients and in cadavers at occiput–C2. In intact subaxial segments, the FE model more closely predicted patient intervertebral motions than did cadavers. With C-2 injury, FE-predicted motions did not differ from cadaver measurements. With C3–4 injury, however, the FE model predicted greater motions than were measured in cadavers. FE model cord strains during intubation were greater for the Macintosh laryngoscope than the Airtraq laryngoscope but were comparable among the 3 conditions (intact, C-2 injury, and C3–4 injury).

CONCLUSIONS

The FE model is comparable to patients and cadaver models in estimating occiput–C2 motion during intubation in both intact and injured conditions. The FE model may be superior to cadavers in predicting motions of subaxial segments in intact and injured conditions.

Biomechanical analysis of brachial plexus injury: Availability of three-dimensional finite element model of the brachial plexus

Adult brachial plexus injuries frequently lead to significant and permanent physical disabilities. Investigating the mechanism of the injury using biomechanical approaches may lead to further knowledge with regard to preventing brachial plexus injuries. However, there are no reports of biomechanical studies of brachial plexus injuries till date. Therefore, the present study used a complex three-dimensional finite element model (3D-FEM) of the brachial plexus to analyze the mechanism of brachial plexus injury and to assess the validity of the model. A complex 3D-FEM of the spinal column, dura mater, spinal nerve root, brachial plexus, rib bone and cartilage, clavicle, scapula, and humerus were conducted. Stress was applied to the model based on the mechanisms of clinically reported brachial plexus injuries: Retroflexion of the cervical, lateroflexion of the cervical, rotation of the cervical, and abduction of the upper limb. The present study analyzed the distribution and strength of strain applied to the brachial plexus during each motion. When the cervical was retroflexed or lateroflexed, the strain was focused on the C5 nerve root and the upper trunk of the brachial plexus. When the upper limb was abducted, strain was focused on the C7 and C8 nerve roots and the lower trunk of the brachial plexus. The results of brachial plexus injury mechanism corresponded with clinical findings that demonstrated the validity of this model. The results of the present study hypothesized that the model has a future potential for analyzing pathological conditions of brachial plexus injuries and other injuries or diseases, including that of spine and spinal nerve root.

Damage Accumulation Modeling and Rate Dependency of Spinal Dura Mater

As the strongest of the meningeal tissues, the spinal dura mater plays an important role in the overall behavior of the spinal cord-meningeal complex (SCM). It follows that the accumulation of damage affects the dura mater's ability to protect the cord from excessive mechanical loads. Unfortunately, current computational investigations of spinal cord injury (SCI) etiology typically do not include postyield behavior. Therefore, a more detailed description of the material behavior of the spinal dura mater, including characterization of damage accumulation, is required to comprehensively study SCI. Continuum mechanics-based viscoelastic damage theories have been previously applied to other biological tissues; however, the current work is the first to report damage accumulation modeling in a tissue of the SCM complex. Longitudinal (i.e., cranial-to-caudal long-axis) samples of ovine cervical dura mater were tensioned-to-failure at one of three strain rates (quasi-static, 0.05/s, and 0.3/s). The resulting stress–strain data were fit to a hyperelastic continuum damage model to characterize the strain-rate-dependent subfailure and failure behavior. The results show that the damage behavior of the fibrous and matrix components of the dura mater are strain-rate dependent, with distinct behaviors when exposed to strain rates above that experienced during normal voluntary neck motion suggesting the possible existence of a protective mechanism.

Effects of Low Bone Mineral Status on Biomechanical Characteristics in Idiopathic Scoliotic Spinal Deformity

BACKGROUND

Low bone mass in patients with adolescent idiopathic scoliosis has been well reported. Poor bone quality was regarded as a new and unique prognostic factor in aggravating curve progression. However, the potential biomechanical correlation between them remains unclear.

METHODS

Three-dimensional finite element models of idiopathic scoliotic spine with different bone mineral status were created for axial loading simulation. An axial load of 3 different body weights was applied on different bone mineral mass models. The mechanical responses of the vertebral cortical and cancellous bone, facet joints, end plate, and intervertebral disc were analyzed.

RESULTS

Accompanied with the low bone mineral status, thoracic scoliosis produced asymmetric and higher stress in the cortical bone, lumbar facet joints, and end plate at the concave side of the thoracic structure curve. Stress increased in the disc at the apex of the scoliosis, whereas it mildly decreased in the L4-5 and L5-S1 disc. Body weight gain increased the stress in scoliotic spine structures in all bone mineral statues.

CONCLUSIONS

Biomechanical simulations indicated that low bone mineral mass might aggravate curve progression and induce more serious lumbar compensatory scoliosis in patients with adolescent idiopathic scoliosis. Weight gain was also a risk factor for curve progression.

Increased stress and strain on the spinal cord due to ossification of the posterior longitudinal ligament in the cervical spine under flexion after laminectomy

Myelopathy in the cervical spine due to cervical ossification of the posterior longitudinal ligament could be induced by static compression and/or dynamic factors. It has been suggested that dynamic factors need to be considered when planning and performing the decompression surgery on patients with the ossification of the posterior longitudinal ligament. A finite element model of the C2-C7 cervical spine in the neutral position was developed and used to generate flexion and extension of the cervical spine. The segmental ossification of the posterior longitudinal ligament on the C5 was assumed, and laminectomy was performed on C4-C6 according to a conventional surgical technique. For various occupying ratios of the ossified ligament between 20% and 60%, von-Mises stresses, maximum principal strains in the spinal cord, and cross-sectional area of the cord were investigated in the pre-operative and laminectomy models under flexion, neutral position, and extension. The results were consistent with previous experimental and computational studies in terms of stress, strain, and cross-sectional area. Flexion leads to higher stresses and strains in the cord than the neutral position and extension, even after decompression surgery. These higher stresses and strains might be generated by residual compression occurring at the segment with the ossification of the posterior longitudinal ligament. This study provides fundamental information under different neck positions regarding biomechanical characteristics of the spinal cord in cervical ossification of the posterior longitudinal ligament.

Rodent, large animal and non-human primate models of spinal cord injury

In this narrative review we aimed to assess the usefulness of the different animal models in identifying injury mechanisms and developing therapies for humans suffering from spinal cord injury (SCI). Results obtained from rodent studies are useful but, due to the anatomical, molecular and functional differences, confirmation of these findings in large animals or non-human primates may lead to basic discoveries that cannot be made in rodent models and that are more useful for developing treatment strategies in humans. SCI in dogs can be considered as intermediate between rodent models and human clinical trials, but the primate models could help to develop appropriate methods that might be more relevant to humans. Ideally, an animal model should meet the requirements of availability and repeatability as well as reproduce the anatomical features and the clinical pathological changing process of SCI. An animal model that completely simulates SCI in humans does not exist. The different experimental models of SCI have advantages and disadvantages for investigating the different aspects of lesion development, recovery mechanisms and potential therapeutic interventions. The potential advantages of non-human primate models include genetic similarities, similar caliber/length of the spinal cord as well as biological and physiological responses to injury which are more similar to humans. Among the potential disadvantages, high operating costs, infrastructural requirements and ethical concerns should be considered. The translation from experimental repair strategies to clinical applications needs to be investigated in future carefully designed studies.

Cervical ossification of the posterior longitudinal ligament: factors affecting the effect of posterior decompression

OBJECTIVE

Decompression procedures for cervical myelopathy of ossification of the posterior longitudinal ligament (OPLL) are anterior decompression with fusion, laminoplasty, and posterior decompression with fusion. Preoperative and postoperative stress analyses were performed for compression from hill-shaped cervical OPLL using 3-dimensional finite element method (FEM) spinal cord models.

METHODS

Three FEM models of vertebral arch, OPLL, and spinal cord were used to develop preoperative compression models of the spinal cord to which 10%, 20%, and 30% compression was applied; a posterior compression with fusion model of the posteriorly shifted vertebral arch; an advanced kyphosis model following posterior decompression with the spinal cord stretched in the kyphotic direction; and a combined model of advanced kyphosis following posterior decompression and intervertebral mobility. The combined model had discontinuity in the middle of OPLL, assuming the presence of residual intervertebral mobility at the level of maximum cord compression, and the spinal cord was mobile according to flexion of vertebral bodies by 5°, 10°, and 15°.

RESULTS

In the preoperative compression model, intraspinal stress increased as compression increased. In the posterior decompression with fusion model, intraspinal stress decreased, but partially persisted under 30% compression. In the advanced kyphosis model, intraspinal stress increased again. As anterior compression was higher, the stress increased more. In the advanced kyphosis  + intervertebral mobility model, intraspinal stress increased more than in the only advanced kyphosis model following decompression. Intraspinal stress increased more as intervertebral mobility increased.

CONCLUSION

In high residual compression or instability after posterior decompression, anterior decompression with fusion or posterior decompression with instrumented fusion should be considered.

Strain rate dependent behavior of the porcine spinal cord under transverse dynamic compression

The accurate description of the mechanical properties of spinal cord tissue benefits to clinical evaluation of spinal cord injuries and is a required input for analysis tools such as finite element models. Unfortunately, available data in the literature generally relate mechanical properties of the spinal cord under quasi-static loading conditions, which is not adapted to the study of traumatic behavior, as neurological tissue adopts a viscoelastic behavior. Thus, the objective of this study is to describe mechanical properties of the spinal cord up to mechanical damage, under dynamic loading conditions. A total of 192 porcine cervical to lumbar spinal cord samples were compressed in a transverse direction. Loading conditions included ramp tests at 0.5, 5 or 50 s−1 and cyclic loading at 1, 10 or 20 Hz. Results showed that spinal cord behavior was significantly influenced by strain rate. Mechanical damage occurred at 0.64, 0.68 and 0.73 strains for 0.5, 5 or 50 s−1 loadings, respectively. Variations of behavior between the tested strain rates were explained by cyclic loading results, which revealed behavior more or less viscous depending on strain rate. Also, a parameter (stress multiplication factor) was introduced to allow transcription of a stress–strain behavior curve to different strain rates. This factor was described and was significantly different for cervical, thoracic and lumbar vertebral heights, and for the strain rates evaluated in this study.

Biomechanical Behaviors in Three Types of Spinal Cord Injury Mechanisms

Clinically, spinal cord injuries (SCIs) are radiographically evaluated and diagnosed from plain radiographs, computed tomography (CT), and magnetic resonance imaging. However, it is difficult to conclude that radiographic evaluation of SCI can directly explain the fundamental mechanism of spinal cord damage. The von-Mises stress and maximum principal strain are directly associated with neurological damage in the spinal cord from a biomechanical viewpoint. In this study, the von-Mises stress and maximum principal strain in the spinal cord as well as the cord cross-sectional area (CSA) were analyzed under various magnitudes for contusion, dislocation, and distraction SCI mechanisms, using a finite-element (FE) model of the cervical spine with spinal cord including white matter, gray matter, dura mater with nerve roots, and cerebrospinal fluid (CSF). A regression analysis was performed to find correlation between peak von-Mises stress/peak maximum principal strain at the cross section of the highest reduction in CSA and corresponding reduction in CSA of the cord. Dislocation and contusion showed greater peak stress and strain values in the cord than distraction. The substantial increases in von-Mises stress as well as CSA reduction similar to or more than 30% were produced at a 60% contusion and a 60% dislocation, while the maximum principal strain was gradually increased as injury severity elevated. In addition, the CSA reduction had a strong correlation with peak von-Mises stress/peak maximum principal strain for the three injury mechanisms, which might be fundamental information in elucidating the relationship between radiographic and mechanical parameters related to SCI.

Stress analysis of the cervical spinal cord: Impact of the morphology of spinal cord segments on stress

OBJECTIVE

Although there are several classifications for cervical myelopathy, these do not take differences between spinal cord segments into account. Moreover, there has been no report of stress analyses for individual segments to date.

METHODS

By using the finite element method, we constructed 3-dimensional spinal cord models comprised of gray matter, white matter, and pia mater of the second to eighth cervical vertebrae (C2-C8). We placed compression components (disc and yellow ligament) at the front and back of these models, and applied compression to the posterior section covering 10%, 20%, 30%, or 40% of the anteroposterior diameter of each cervical spinal cord segment.

RESULTS

Our results revealed that, under compression applied to an area covering 10%, 20%, or 30% of the anteroposterior diameter of the cervical spinal cord segment, sites of increased stress varied depending on the morphology of each cervical spinal cord segment. Under 40% compression, stress was increased in the gray matter, lateral funiculus, and posterior funiculus of all spinal cord segments, and stress differences between the segments were smaller.

CONCLUSION

These results indicate that, under moderate compression, sites of increased stress vary depending on the morphology of each spinal cord segment or the shape of compression components, and also that the variability of symptoms may depend on the direction of compression. However, under severe compression, the differences among the cervical spinal segments are smaller, which may facilitate diagnosis.

Effect of posterior decompression extent on biomechanical parameters of the spinal cord in cervical ossification of the posterior longitudinal ligament

Ossification of the posterior longitudinal ligament is a common cause of the cervical myelopathy due to compression of the spinal cord. Patients with ossification of the posterior longitudinal ligament usually require the decompression surgery, and there is a need to better understand the optimal surgical extent with which sufficient decompression without excessive posterior shifting can be achieved. However, few quantitative studies have clarified this optimal extent for decompression of cervical ossification of the posterior longitudinal ligament. We used finite element modeling of the cervical spine and spinal cord to investigate the effect of posterior decompression extent for continuous-type cervical ossification of the posterior longitudinal ligament on changes in stress, strain, and posterior shifting that occur with three different surgical methods (laminectomy, laminoplasty, and hemilaminectomy). As posterior decompression extended, stress and strain in the spinal cord decreased and posterior shifting of the cord increased. The location of the decompression extent also influenced shifting. Laminectomy and laminoplasty were very similar in terms of decompression results, and both were superior to hemilaminectomy in all parameters tested. Decompression to the extents of C3-C6 and C3-C7 of laminectomy and laminoplasty could be considered sufficient with respect to decompression itself. Our findings provide fundamental information regarding the treatment of cervical ossification of the posterior longitudinal ligament and can be applied to patient-specific surgical planning.

Cervical ossification of the posterior longitudinal ligament: Biomechanical analysis of the influence of static and dynamic factors

OBJECTIVE

Cervical myelopathy due to ossification of the posterior longitudinal ligament (OPLL) is induced by static factors, dynamic factors, or a combination of both. We used a three-dimensional finite element method (3D-FEM) to analyze the stress distributions in the cervical spinal cord under static compression, dynamic compression, or a combination of both in the context of OPLL.

METHODS

Experimental conditions were established for the 3D-FEM spinal cord, lamina, and hill-shaped OPLL. To simulate static compression of the spinal cord, anterior compression at 10, 20, and 30% of the anterior-posterior diameter of the spinal cord was applied by the OPLL. To simulate dynamic compression, the OPLL was rotated 5°, 10°, and 15° in the flexion direction. To simulate combined static and dynamic compression under 10 and 20% anterior static compression, the OPLL was rotated 5°, 10°, and 15° in the flexion direction.

RESULTS

The stress distribution in the spinal cord increased following static and dynamic compression by cervical OPLL. However, the stress distribution did not increase throughout the entire spinal cord. For combined static and dynamic compression, the stress distribution increased as the static compression increased, even for a mild range of motion (ROM).

CONCLUSION

Symptoms may appear under static or dynamic compression only. However, under static compression, the stress distribution increases with the ROM of the responsible level and this makes it very likely that symptoms will worsen. We conclude that cervical OPLL myelopathy is induced by static factors, dynamic factors, and a combination of both.

Determination of the ideal rat model for spinal cord injury by diffusion tensor imaging

Four different spinal cord injury (SCI) models (hemisection, contusion, transection, and segment resection) were produced in male Sprague-Dawley rats to determine the most suitable animal model of SCI by analyzing the changes in diffusion tensor imaging (DTI) parameters both qualitatively and quantitatively in vivo. Radiological examinations were performed before surgery and weekly within 4 weeks after surgery to obtain DTI tractography, MRI routine images, and DTI data of fractional anisotropy (FA) and apparent diffusion coefficient (ADC). The Basso, Beattie, and Bresnahan scale was used to evaluate the locomotor outcomes. We found that DTI tractography tracked nerve fibers and showed conspicuous changes in the injured spinal cord in all the model groups, which confirmed that our modeling was successful. A decrease in FA values and an increase in ADC were observed in all the model groups after surgery. There were significant differences in FA and ADC between weeks 1 and 4 in both hemisection and contusion groups (P<0.05), whereas the differences in the transection and segment resection groups were not as remarkable (P>0.05). Basso, Beattie, and Bresnahan scores further proved the results because of a significant, positive correlation of the scores with FA (R=0.899, P<0.01) and a significant, negative correlation of the scores with ADC (R=-0.829, P<0.01). Therefore, the transection model, which is more quantified and stable within 4 weeks after injury according to the DTI and behavioral evaluation, should be used as the standard model for SCI animal testing.

A novel closed-body model of spinal cord injury caused by high-pressure air blasts produces extensive axonal injury and motor impairments

Diffuse axonal injury is thought to be the basis of the functional impairments stemming from mild traumatic brain injury. To examine how axons are damaged by traumatic events, such as motor vehicle accidents, falls, sports activities, or explosive blasts, we have taken advantage of the spinal cord with its extensive white matter tracts. We developed a closed-body model of spinal cord injury in mice whereby high-pressure air blasts targeted to lower thoracic vertebral levels produce tensile, compressive, and shear forces within the parenchyma of the spinal cord and thereby cause extensive axonal injury. Markers of cytoskeletal integrity showed that spinal cord axons exhibited three distinct pathologies: microtubule breakage, neurofilament compaction, and calpain-mediated spectrin breakdown. The dorsally situated axons of the corticospinal tract primarily exhibited microtubule breakage, whereas all three pathologies were common in the lateral and ventral white matter. Individual axons typically demonstrated only one of the three pathologies during the first 24h after blast injury, suggesting that the different perturbations are initiated independently of one another. For the first few days after blast, neurofilament compaction was frequently accompanied by autophagy, and subsequent to that, by the fragmentation of degenerating axons. TuJ1 immunolabeling and mice with YFP-reporter labeling each revealed more extensive microtubule breakage than did βAPP immunolabeling, raising doubts about the sensitivity of this standard approach for assessing axonal injury. Although motor deficits were mild and largely transient, some aspects of motor function gradually worsened over several weeks, suggesting that a low level of axonal degeneration continued past the initial wave. Our model can help provide further insight into how to intervene in the processes by which initial axonal damage culminates in axonal degeneration, to improve outcomes after traumatic injury. Importantly, our findings of extensive axonal injury also caution that repeated trauma is likely to have cumulative adverse consequences for both brain and spinal cord.

Meta analysis of olfactory ensheathing cell transplantation promoting functional recovery of motor nerves in rats with complete spinal cord transection

OBJECTIVE:

To evaluate the effects of olfactory ensheathing cell transplantation on functional recovery of rats with complete spinal cord transection.

DATA SOURCES:

A computer-based online search of Medline (1989–2013), Embase (1989–2013), Cochrane library (1989–2013), Chinese Biomedical Literature Database (1989–2013), China National Knowledge Infrastructure (1989–2013), VIP (1989–2013), Wanfang databases (1989–2013) and Chinese Clinical Trial Register was conducted to collect randomized controlled trial data regarding olfactory ensheathing cell transplantation for the treatment of complete spinal cord transection in rats.

SELECTION CRITERIA:

Randomized controlled trials investigating olfactory ensheathing cell transplantation and other transplantation methods for promoting neurological functional recovery of rats with complete spinal cord transection were included in the analysis. Meta analysis was conducted using RevMan 4.2.2 software.

MAIN OUTCOME MEASURES:

Basso, Beattie and Bresnahan scores of rats with complete spinal cord transection were evaluated in this study.

RESULTS:

Six randomized controlled trials with high quality methodology were included. Meta analysis showed that Basso, Beattie and Bresnahan scores were significantly higher in the olfactory ensheathing cell transplantation group compared with the control group (WMD = 3.16, 95% CI (1.68, 4.65); P < 0.00001).

CONCLUSION:

Experimental studies have shown that olfactory ensheathing cell transplantation can promote the functional recovery of motor nerves in rats with complete spinal cord transection.

Biomechanical characteristics of the porcine denticulate ligament in different vertebral levels of the cervical spine-preliminary results of an experimental study

BACKGROUND

Few studies exist on the mechanical properties of denticulate ligaments and none report the variation in these properties at different levels of the spine. The aim of this study was to perform an experimental determination of load-extension and stress-strain characteristics of the denticulate ligament and to establish if their properties change at different vertebral levels of the cervical spine.

METHOD

The study was carried out on a total of 98 porcine denticulate ligament samples dissected from seven fresh porcine cervical spinal cord specimens. All of the samples were subjected to an uniaxial tensile test at a speed of 2mm/min, during which the load-extension characteristics were registered.

RESULTS

The analysis revealed a decrease of the failure force in the caudal orientation indicated by significant differences between the C1 (1.04±0.41N) and C7 (0.55±0.12N) vertebral levels (P=0.037). The average ultimate force that broke the denticulate ligaments was 0.88N. The mean value of Young׳s modulus was 2.06MPa with a minimum of 1.31MPa for C7 and maximum of 2.46MPa for C5.

CONCLUSIONS

The values of the denticulate ligament failure force in samples from different cervical vertebrae levels differ significantly. The presented data should be taken into consideration during numerical modelling of the human cervical spinal cord.

Biomechanical analysis of cervical myelopathy due to ossification of the posterior longitudinal ligament: Effects of posterior decompression and kyphosis following decompression

Cervical ossification of the posterior longitudinal ligament (OPLL) results in myelopathy. Conservative treatment is usually ineffective, thus, surgical treatment is required. One of the reasons for the poor surgical outcome following laminoplasty for cervical OPLL is kyphosis. In the present study, a 3-dimensional finite element method (3D-FEM) was used to analyze the stress distribution in preoperative, posterior decompression and kyphosis models of OPLL. The 3D-FEM spinal cord model established in this study consisted of gray and white matter, as well as pia mater. For the preoperative model, 30% anterior static compression was applied to OPLL. For the posterior decompression model, the lamina was shifted backwards and for the kyphosis model, the spinal cord was studied at 10, 20, 30, 40 and 50° kyphosis. In the preoperative model, high stress distributions were observed in the spinal cord. In the posterior decompression model, stresses were lower than those observed in the preoperative model. In the kyphosis model, an increase in the angle of kyphosis resulted in augmented stress on the spinal cord. Therefore, the results of the present study indicated that posterior decompression was effective, but stress distribution increased with the progression of kyphosis. In cases where kyphosis progresses following surgery, detailed follow-ups are required in case the symptoms worsen.

Morphometrics of the entire human spinal cord and spinal canal measured from in vivo high-resolution anatomical magnetic resonance imaging

STUDY DESIGN

Measurements of cervical and thoracolumbar human spinal cord (SC) geometry based on in vivo magnetic resonance imaging and investigation of morphological "invariants."

OBJECTIVE

The current work aims at providing morphological features of the complete in vivo human normal SC and at investigating possible "invariant" parameters that may serve as normative data for individualized study of SC injuries.

SUMMARY OF BACKGROUND DATA

Few in vivo magnetic resonance image-based studies have described human SC morphology at the cervical level, and similar description of the entire SC only relies on postmortem studies, which may be prone to atrophy biases. Moreover, large interindividual variations currently limit the use of morphological metrics as reference for clinical applications or as modeling inputs.

METHODS

Absolute metrics of SC (transverse and anteroposterior diameters, width of anterior and posterior horns, cross-sectional SC area, and white matter percentage) were measured using semiautomatic segmentation of high resolution in vivo T2*-weighted transverse images acquired at 3 T, at each SC level, on healthy young (N = 15) and older (N = 8) volunteers. Robustness of measurements, effects of subject, age, or sex, as well as comparison with previously published postmortem data were investigated using statistical analyses (separate analysis of variance, Tukey-HSD, Bland-Altman). Normalized-to-C3 parameters were evaluated as invariants using a leave-one-out analysis. Spinal canal parameters were measured and occupation ratio border values were determined.

RESULTS

Metrics of SC morphology showed large intra- and interindividual variations, up to 30% and 13%, respectively, on average. Sex had no influence except on posterior horn width (P < 0.01). Age-related differences were observed for anteroposterior diameter and white matter percentage (P < 0.05) and all postmortem metrics were significantly lower than in vivo values (P < 0.001). In vivo normalized SC area and diameters seemed to be invariants (R > 0.74, root-mean-square deviation < 10%). Finally, minimal and maximal occupation ratio were 0.2 and 0.6, respectively.

CONCLUSION

This study presented morphological characteristics of the complete in vivo human SC. Significant differences linked to age and postmortem state have been identified. Morphological "invariants" that could be used to calculate the normally expected morphology accurately, were also identified. These observations should benefit to biomechanical and SC pathology studies.

LEVEL OF EVIDENCE

N/A.

Nonlinear viscoelastic characterization of the porcine spinal cord

Although quasi-static and quasi-linear viscoelastic properties of the spinal cord have been reported previously, there are no published studies that have investigated the fully (strain-dependent) nonlinear viscoelastic properties of the spinal cord. In this study, stress relaxation experiments and dynamic cycling were performed on six fresh porcine lumbar cord specimens to examine their viscoelastic mechanical properties. The stress relaxation data were fitted to a modified superposition formulation and a novel finite ramp time correction technique was applied. The parameters obtained from this fitting methodology were used to predict the average dynamic cyclic viscoelastic behavior of the porcine cord. The data indicate that the porcine spinal cord exhibited fully nonlinear viscoelastic behavior. The average weighted root mean squared error for a Heaviside ramp fit was 2.8 kPa, which was significantly greater (p<0.001) than that of the nonlinear (comprehensive viscoelastic characterization method) fit (0.365 kPa). Further, the nonlinear mechanical parameters obtained were able to accurately predict the dynamic behavior, thus exemplifying the reliability of the obtained nonlinear parameters. These parameters will be important for future studies investigating various damage mechanisms of the spinal cord and studies developing high-resolution finite elements models of the spine.

Mechanical properties of the lamprey spinal cord: uniaxial loading and physiological strain

During spinal cord injury, nerves suffer a strain beyond their physiological limits which damages and disrupts their structure. Research has been done to measure the modulus of the spinal cord and surrounding tissue; however the relationship between strain and spinal cord fibers is still unclear. In this work, our objective is to measure the stress-strain response of the spinal cord in vivo and in vitro and model this response as a function of the number of fibers. We used the larvae lamprey (Petromyzon Marinus), a model for spinal cord regeneration and animal locomotion. We found that physiologically the spinal cord is pre-stressed to a longitudinal strain of 10% and this strain increases to 15% during swimming. Tensile measurements show that uniaxial, longitudinal loading is independent of the meninges. Stress values for uniaxial strains below 18%, are homogeneous through the length of the body. However, for higher uniaxial strains the Head section shows more resistance to longitudinal loading than the Tail. These data, together with the number of fibers obtained from histological sections were used in a composite-material model to obtain the properties of the spinal cord fibers (2.4 MPa) and matrix (0.017 MPa) to uniaxial longitudinal loading. This model allowed us to approximate the percentage of fibers in the spinal cord, establishing a relationship between uniaxial longitudinal strains and spinal cord composition. We showed that there is a proportional relationship between the number of fibers and the properties of the spinal cord at large uniaxial strains.

Finite element study of the mechanical response in spinal cord during the thoracolumbar burst fracture

Background

The mechanical response of the spinal cord during burst fracture was seldom quantitatively addressed and only few studies look into the internal strain of the white and grey matters within the spinal cord during thoracolumbar burst fracture (TLBF). The aim of the study is to investigate the mechanical response of the spinal cord during TLBF and correlate the percent canal compromise (PCC) with the strain in the spinal cord.

Methodology/Principal Findings

A three-dimensional (3D) finite element (FE) model of human T12-L1 spinal cord with visco-elastic property was generated based on the transverse sections images of spinal cord, and the model was validated against published literatures under static uniaxial tension and compression. With the validated model, a TLBF simulation was performed to compute the mechanical strain in the spinal cord with the PCC. Linear regressions between PCC and strain in the spinal cord show that at the initial stage, with the PCC at 20%, and 45%, the corresponding mechanical strains in ventral grey, dorsal grey, ventral white, dorsal white matters were 0.06, 0.04, 0.12, 0.06, and increased to 0.14, 0.12, 0.23, and 0.13, respectively. At the recoiled stage, when the PCC was decreased from 45% to 20%, the corresponding strains were reduced to 0.03, 0.02, 0.04 and 0.03. The strain was correlated well with PCC.

Conclusions/Significance

The simulation shows that the strain in the spinal cord correlated well with the PCC, and the mechanical strains in the ventral regions are higher than those in the dorsal regions of spinal cord tissue during burst fracture, suggesting that the ventral regions of the spinal cord may susceptible to injury than the dorsal regions.

Biomechanical analysis of cervical spondylotic myelopathy: the influence of dynamic factors and morphometry of the spinal cord

OBJECTIVE

Patients with cervical spondylotic myelopathy (CSM) have the same clinical symptoms that vary according to the degree of spinal cord compression and the cross-sectional cord shape. We used a three-dimensional finite element method (3D-FEM) to analyze the stress distributions of the spinal cord with neck extension under three cross-sectional cord shapes.

METHODS

Experimental condition for the 3D-FEM spinal cord, ligamentum flavum, and anterior compression shape (central, lateral, and diffuse types) was established. To simulate neck extension, the spinal cord was extended by 20° and the ligamentum flavum was shifted distally according to movement of the cephalad lamina.

RESULTS

The stress distribution in the spinal cord increased due to invagination of the ligamentum flavum into the neck extension. The range of stress distribution observed for the diffuse type was wider than for the central and lateral types. In addition, the stress distribution in the spinal cord was increased by the pincer movement of the ligamentum flavum and by the anterior compression of the spinal cord. The range of stress distribution observed for the diffuse type under antero-posterior compression was also wider than for the central and lateral types.

CONCLUSION

This simulation model showed that the clinical symptoms of CSM due to compression of the diffuse type may be stronger than for the central and lateral types. Therefore, careful follow-up is recommended for anterior compression of the spinal cord of diffuse type.

Maximum principal strain correlates with spinal cord tissue damage in contusion and dislocation injuries in the rat cervical spine

The heterogeneity of the primary mechanical mechanism of spinal cord injury (SCI) is not currently used to tailor treatment strategies because the effects of these distinct patterns of acute mechanical damage on long-term neuropathology have not been fully investigated. A computational model of SCI enables the dynamic analysis of mechanical forces and deformations within the spinal cord tissue that would otherwise not be visible from histological tissue sections. We created a dynamic, three-dimensional finite element (FE) model of the rat cervical spine and simulated contusion and dislocation SCI mechanisms. We investigated the relationship between maximum principal strain and tissue damage, and compared primary injury patterns between mechanisms. The model incorporated the spinal cord white and gray matter, the dura mater, cerebrospinal fluid, spinal ligaments, intervertebral discs, a rigid indenter and vertebrae, and failure criteria for ligaments and vertebral endplates. High-speed (∼ 1 m/sec) contusion and dislocation injuries were simulated between vertebral levels C3 and C6 to match previous animal experiments, and average peak maximum principal strains were calculated for several regions at the injury epicenter and at 1-mm intervals from +5 mm rostral to -5 mm caudal to the lesion. Average peak principal strains were compared to tissue damage measured previously in the same regions via axonal permeability to 10-kD fluorescein-dextran. Linear regression of tissue damage against peak maximum principal strain for pooled data within all white matter regions yielded similar and significant (p<0.0001) correlations for both contusion (R(2)=0.86) and dislocation (R(2)=0.52). The model enhances our understanding of the differences in injury patterns between SCI mechanisms, and provides further evidence for the link between principal strain and tissue damage.

Presence of anterior compression of the spinal cord after laminoplasty inhibits upper extremity motor recovery in patients with cervical spondylotic myelopathy

STUDY DESIGN

A retrospective single-center study.

OBJECTIVE

To investigate how functional recovery is influenced by anterior compression of the spinal cord (ACS) and instability at the level of ACS after laminoplasty in patients with cervical spondylotic myelopathy.

SUMMARY OF BACKGROUND DATA

There have been many reports that patients whose spinal cord cannot be decompressed sufficiently after laminoplasty are likely to show unsatisfactory neurologic outcomes. Notably, postoperative ACS is well known to cause problems. Clinically, however, it remains unknown how functional recovery is inhibited by postoperative ACS.

METHODS

Sixty-four consecutive patients who underwent expansive laminoplasty for the treatment of myelopathy at our hospital between 1998 and 2005 were reviewed. All 64 patients were available for follow-up. The average follow-up period was 97 months (60-156 months). Patients were divided into 2 groups: the ACS(+) group comprised 16 patients who had ACS 3 years postoperatively, and the ACS(-) group comprised 48 patients with no ACS. Clinical outcome was compared in terms of the Japanese Orthopaedic Association score (mean total score, mean score of each item, and recovery rates).

RESULTS

Demographics were similar between the 2 groups. Mean Japanese Orthopaedic Association score at final follow-up was 12.1 points (recovery rate 34.0%) in the ACS(+) group and 13.8 points (recovery rate 56.6%) in the ACS(-) group, and there was a significant difference in recovery rate between the groups (P < 0.05). Notably, a significant difference was found between the 2 groups in improvement of upper extremity motor function (P < 0.05). In addition, we found that not only the presence of ACS but also postoperative hypermobility of the intervertebral segment with ACS influenced clinical outcome negatively.

CONCLUSION

These results demonstrate that ACS after laminoplasty could be a risk factor for clinical outcome and might prevent improvement in upper extremity motor function in patients with myelopathy.