Spinal cord–fragment interactions following burst fracture: an in vitro model

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.

In vitro assessment of the role of the nucleus pulposus in the mechanism of vertebral body fracture under dynamic compressive loading using high-speed cineradiography

Traumatic Spinal Cord Injuries (TSCI) have a disastrous effect on the physical and mental health of both the patients and their relatives. Around 15 % of these injuries are caused by burst fractures, a sub-type of compressive fractures of the vertebral body. The transient dynamics of these fracture have been studied through in vitro experiments coupled with numerical simulations, but no direct observation have ever been made of their genesis and evolution and the behaviour of the nucleus pulposus under compressive loading has only been hypothesized. The purpose of this study was to evaluate the interactions between the vertebral body and the nucleus pulposus under dynamic compressive loading using high-speed cineradiography. A radiopaque agent was injected into the nuclei pulposi of 4 young porcine thoraco-lumbar and lumbar cadaveric segments, and a dynamic compressive load was applied to them using a servo-hydraulic bench-test. The compression process was filmed with a custom high-speed fluoroscope. The nucleus pulposus loaded the vertebral endplate up to 14,142 ± 486 N, before fracturing it and diffusing into the vertebral body. Then, internal pressure seemingly built up until an outward projection of the nucleus pulposus, at an antero-posterior velocity up to 2.9 m.s^-1, or until retroprojection of bony fragments into the spinal canal. These results directly corroborate the hypotheses previously made by other studies and stress the unprecedented advantages of using high-speed cineradiography for the study of complex fractures genesis and evolution.Clinical Relevance- Methodology and results from this study would provide an unprecedented insight on the genesis and transient evolution of complex spinal fractures.

Dynamics of spinal cord compression with different patterns of thoracolumbar burst fractures: Numerical simulations using finite element modelling

BACKGROUND

In thoracolumbar burst fractures, spinal cord primary injury involves a direct impact and energy transfer from bone fragments to the spinal cord. Unfortunately, imaging studies performed after the injury only depict the residual bone fragments position and pattern of spinal cord compression, with little insight on the dynamics involved during traumas. Knowledge of underlying mechanisms could be helpful in determining the severity of the primary injury, hence the extent of spinal cord damage and associated potential for recovery. Finite element models are often used to study dynamic processes, but have never been used specifically to simulate different severities of thoracolumbar burst fractures.

METHODS

Previously developed thoracolumbar spine and spinal cord finite element models were used and further validated, and representative vertebral fragments were modelled. A full factorial design was used to investigate the effects of comminution of the superior fragment, presence of an inferior fragment, fragments rotation and velocity, on maximum Von Mises stress and strain, maximum major strain, and pressure in the spinal cord.

FINDINGS

Fragment velocity clearly was the most influential factor. Fragments rotation and presence of an inferior fragment increased pressure, but rotation decreased both strains outputs. Although significant for both strains outputs, comminution of the superior fragment isn't estimated to influence outputs.

INTERPRETATION

This study is the first, to the authors' knowledge, to examine a detailed spinal cord model impacted in situ by fragments from burst fractures. This numeric model could be used in the future to comprehensively link traumatic events or imaging study characteristics to known spinal cord injuries severity and potential for recovery.

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.

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.

In Vitro Models as a Platform to Investigate Traumatic Brain Injury

Traumatic Brain Injury (TBI) remains a significant cause of mortality and morbidity, affecting individuals of all age groups. Much remains to be learned about its complex pathophysiology, with a view to designing effective neuroprotective strategies to protect sublethally injured brain tissue that would otherwise die in secondary injury processes. Experimental in vivo models offer the potential to study TBI in the laboratory, however, treatments that were neuroprotective in animals have, thus far, largely failed to translate in human clinical studies. In vitro models of neurotrauma can be used to study specific pathophysiological cascades — individually and without confounding factors — and to test potential neuroprotective strategies. These in vitro models include transection, compression, barotrauma, acceleration, hydrodynamic, chemical injury and cell-stretch methodologies. Various cell culture systems can also be utilised, including brain-on-a-chip, immortalised cell lines, primary cultures, acute preparations and organotypic cultures. Potential positive outcomes of the increased use of in vitro platforms to study TBI would be the refinement of in vivo experiments, as well as enhanced translation of the results into clinically meaningful neuroprotective strategies for the future. In addition, the replacement of in vivo experiments by suitable in vitro studies would lead to a welcome reduction in the numbers of animal procedures in this ethically-challenging field.

Cerebrospinal Fluid Pressures Resulting From Experimental Traumatic Spinal Cord Injuries in a Pig Model

Despite considerable effort over the last four decades, research has failed to translate into consistently effective treatment options for spinal cord injury (SCI). This is partly attributed to differences between the injury response of humans and rodent models. Some of this difference could be because the cerebrospinal fluid (CSF) layer of the human spine is relatively large, while that of the rodents is extremely thin. We sought to characterize the fluid impulse induced in the CSF by experimental SCIs of moderate and high human-like severity, and to compare this with previous studies in which fluid impulse has been associated with neural tissue injury. We used a new in vivo pig model (n = 6 per injury group, mean age 124.5 days, 20.9 kg) incorporating four miniature pressure transducers that were implanted in pairs in the subarachnoid space, cranial, and caudal to the injury at 30 mm and 100 mm. Tissue sparing was assessed with Eriochrome Cyanine and Neutral Red staining. The median peak pressures near the injury were 522.5 and 868.8 mmHg (range 96.7–1430.0) and far from the injury were 7.6 and 36.3 mmHg (range 3.8–83.7), for the moderate and high injury severities, respectively. Pressure impulse (mmHg.ms), apparent wave speed, and apparent attenuation factor were also evaluated. The data indicates that the fluid pressure wave may be sufficient to affect the severity and extent of primary tissue damage close to the injury site. However, the CSF pressure was close to normal physiologic values at 100 mm from the injury. The high injury severity animals had less tissue sparing than the moderate injury severity animals; this difference was statistically significant only within 1.6 mm of the epicenter. These results indicate that future research seeking to elucidate the mechanical origins of primary tissue damage in SCI should consider the effects of CSF. This pig model provides advantages for basic and preclinical SCI research due to its similarities to human scale, including the existence of a human-like CSF fluid layer.

The pressure distribution of cerebrospinal fluid responds to residual compression and decompression in an animal model of acute spinal cord injury

STUDY DESIGN

In vivo large animal (pig) model study of cerebrospinal fluid (CSF) pressures after acute experimental spinal cord injury (SCI).

OBJECTIVE

To determine how the CSF pressure (CSFP) and CSF pulse pressure amplitude (CSFPPA) cranial and caudal to the injury site change after an acute SCI with subsequent thecal occlusion and decompression.

SUMMARY OF BACKGROUND DATA

Lowering intrathecal pressure via CSF drainage is currently instituted to prevent ischemia-induced SCI during thoracoabdominal aortic aneurysm surgery and was recently investigated as a potential intervention for acute traumatic SCI. However, in SCI patients, persistent extradural compression commonly occludes the subarachnoid space. This may generate a CSFP differential across the injury site, which cannot be appreciated with lumbar catheter pressure measurements.

METHODS

Anesthetized pigs were subjected to an acute contusive SCI at T11 and 8 hours of sustained compression (n = 12), or sham surgery (n = 2). CSFP was measured cranial and caudal to the injury site, using miniature pressure transducers, during compression and for 6 hours after decompression.

RESULTS

The cranial-caudal CSFP differential increased (mean, 0.39 mm Hg/h), predominantly due to increased cranial pressure. On decompression, cranial CSFP decreased (mean, -1.16 mm Hg) and caudal CSFP increased (mean, 0.65 mm Hg). The CSFP differential did not change significantly after decompression. Cranial CSFPPA was greater than caudal CSFPPA, but this differential did not change during compression. On decompression, the caudal CSFPPA increased in some but not all animals.

CONCLUSION

Although extradural compression exists at the site of injury, lumbar CSFP may not accurately indicate CSFP cranial to the injury. Decompression may provide immediate, though perhaps partial, resolution of the pressure differential. CSFPPA was not a consistent indicator of decompression in this animal model. These findings may have implications for the design of future clinical protocols in which CSFP is monitored after acute SCI.

Gross morphological changes of the spinal cord immediately after surgical decompression in a large animal model of traumatic spinal cord injury

STUDY DESIGN

Quantitative in vivo ultrasound imaging study of spinal cord and dura morphology after acute experimental spinal cord injury (SCI) and decompression in a pig model.

OBJECTIVE

To study the morphological changes of the spinal cord and dura immediately after surgical decompression for acute SCI.

SUMMARY OF BACKGROUND DATA

Surgical decompression for traumatic SCI is currently a topic of debate. After decompression, relief of bony impingement on the thecal sac and spinal cord can be confirmed intraoperatively. However, postoperative imaging often reveals that the cord has swollen to fill the subarachnoid space. Little is known about the extent and timing of this morphological response.

METHODS

Yucatan miniature pigs received sham surgery (N = 1) or a moderate (N = 6, 20 g, 2.3 m/s) or high (N = 6, 20 g, 4.7 m/s) severity weight-drop SCI followed by 8 hours of sustained compression (100 g) and 6 hours of postdecompression monitoring. Sagittal-plane ultrasound images were used to quantify spinal cord, dura, and subarachnoid space dimensions preinjury and once per hour after decompression.

RESULTS

Animals with a moderate SCI exhibited a residual cord deformation of up to 0.64 mm within 10 minutes of decompression, which tended to resolve during 6 hours because of tissue relaxation and swelling. For animals with high-severity SCIs, cord swelling was immediate and resulted in occlusion of the subarachnoid space within 10 minutes to 5 hours, whereas this occurred for only half of the moderate injury group.

CONCLUSION

Decompression of an acute SCI may result in residual cord deformation followed by gradual swelling or immediate swelling leading to subarachnoid occlusion. The response is dependent on initial injury severity. These observations may partly explain the lack of benefit of decompression in some patients and suggest a need to reduce cord swelling to optimize the clinical outcome after acute SCI.

Development of a large-animal model to measure dynamic cerebrospinal fluid pressure during spinal cord injury

Object

Spinal cord injury (SCI) often results in considerable permanent neurological impairment, and unfortunately, the successful translation of effective treatments from laboratory models to human patients is lacking. This may be partially attributed to differences in anatomy, physiology, and scale between humans and rodent models. One potentially important difference between the rodent and human spinal cord is the presence of a significant CSF volume within the intrathecal space around the human cord. While the CSF may “cushion” the spinal cord, pressure waves within the CSF at the time of injury may contribute to the extent and severity of the primary injury. The objective of this study was to develop a model of contusion SCI in a miniature pig and establish the feasibility of measuring spinal CSF pressure during injury.

Methods

A custom weight-drop device was used to apply thoracic contusion SCI to 17 Yucatan miniature pigs. Impact load and velocity were measured. Using fiber optic pressure transducers implanted in the thecal sac, CSF pressures resulting from 2 injury severities (caused by 50-g and 100-g weights released from a 50-cm height) were measured.

Results

The median peak impact loads were 54 N and 132 N for the 50-g and 100-g injuries, respectively. At a nominal 100 mm from the injury epicenter, the authors observed a small negative pressure peak (median −4.6 mm Hg [cranial] and −5.8 mm Hg [caudal] for 50 g; −27.6 mm Hg [cranial] and −27.2 mm Hg [caudal] for 100 g) followed by a larger positive pressure peak (median 110.5 mm Hg [cranial] and 77.1 mm Hg [caudal] for 50 g; 88.4 mm Hg [cranial] and 67.2 mm Hg [caudal] for 100 g) relative to the preinjury pressure. There were no significant differences in peak pressure between the 2 injury severities or the caudal and cranial transducer locations.

Conclusions

A new model of contusion SCI was developed to measure spinal CSF pressures during the SCI event. The results suggest that the Yucatan miniature pig is an appropriate model for studying CSF, spinal cord, and dura interactions during injury. With further development and characterization it may be an appropriate in vivo largeanimal model of SCI to answer questions regarding pathological changes, therapeutic safety, or treatment efficacy, particularly where humanlike dimensions and physiology are important.

The importance of fluid-structure interaction in spinal trauma models

While recent studies have demonstrated the importance of the initial mechanical insult in the severity of spinal cord injury, there is a lack of information on the detailed cord-column interaction during such events. In vitro models have demonstrated the protective properties of the cerebrospinal fluid, but visualization of the impact is difficult. In this study a computational model was developed in order to clarify the role of the cerebrospinal fluid and provide a more detailed picture of the cord-column interaction. The study was validated against a parallel in vitro study on bovine tissue. Previous assumptions about complete subdural collapse before any cord deformation were found to be incorrect. Both the presence of the dura mater and the cerebrospinal fluid led to a reduction in the longitudinal strains within the cord. The division of the spinal cord into white and grey matter perturbed the bone fragment trajectory only marginally. In conclusion, the cerebrospinal fluid had a significant effect on the deformation pattern of the cord during impact and should be included in future models. The type of material models used for the spinal cord and the dura mater were found to be important to the stress and strain values within the components, but less important to the fragment trajectory.

Poisson's ratio and strain rate dependency of the constitutive behavior of spinal dura mater

Knowledge of the mechanical behavior of spinal dura mater is important for a number of applications including the experimental and computational modeling of physiological phenomena and spinal cord trauma. However, mechanical characterization of dura mater is relatively sparse and is further compounded by the use of the tangent modulus as the sole measure of stiffness. This study aims to provide a more complete description of the mechanical properties of spinal dura mater, including the effect of strain rate. Bovine dura mater was tested under uniaxial tension in both the longitudinal and the circumferential directions at three different strain rates; 0.01, 0.1, and 1.0 s(-1). An Ogden model was fitted to the resulting stress-stretch data. The morphology of the dura mater was assessed using Sirius red and H&E staining. No significant effect of the strain rate was found for the Ogden model parameters. Longitudinal specimens were significantly stronger and more deformable than circumferential samples, probably due to the structural arrangement of the collagen fibers. At low strains, however, the circumferential specimens were stiffer than the longitudinal ones. The findings of this study will allow more complete representations of the spinal dura mater to be developed.

The effect of bone fragment size and cerebrospinal fluid on spinal cord deformation during trauma: an ex vivo study

Object

The purpose of the study was to assess the effect of CSF and the size of the impacting bone fragment area on spinal cord deformation during trauma.

Methods

A transverse impact rig was used to produce repeated impacts on bovine and surrogate cord models. Tests were recorded with high-speed video and performed on specimens with and without CSF and/or dura mater and with 3 different impactor areas.

Results

The CSF layer was found to reduce the maximum cord deformation significantly. A 50% reduction in impact area significantly increased the maximum cord deformation by 20–30%. The surrogate model showed similar trends to the bovine model but with lower absolute deformation values.

Conclusions

Cerebrospinal fluid protects the cord during impact by reducing its deformation. A smaller bone fragment impact area increases the deformation of the cord, in agreement with clinical results, where a higher impact energy—possibly giving rise to smaller fragments—results in a worse neurological deficit.

The effect of cerebrospinal fluid on the biomechanics of spinal cord: an ex vivo bovine model using bovine and physical surrogate spinal cord

STUDY DESIGN

A biomechanical study using ex vivo bovine spinal cord and dura, and a synthetic surrogate spinal cord with bovine dura.

OBJECTIVE

To investigate the effect of cerebrospinal fluid (CSF) on spinal cord deformation characteristics and to evaluate the biofidelity of a new surrogate spinal cord using an ex vivo bovine model of the burst fracture process.

SUMMARY OF BACKGROUND DATA

Spinal cord injury is associated with significant personal, economic and social costs. The role of CSF during the injury event and its effect on the spinal cord deformation and neurologic injury is not well understood. Such knowledge could inform preventative strategies and clinical interventions and aid the development and validation of experimental and computational models.

METHODS

The transverse impact of a propelled bone fragment analogue with bovine and surrogate cord models was recorded with high speed video and the images analyzed to determine deformation trajectories. Each cord specimen was tested in 3 states: with dura and CSF, with dura only, and without dura. The effect of these states on deformation magnitude, duration, and energy loss parameters was assessed. RESULTS.: The estimated spinal cord deformation was significantly reduced, although not eliminated, in the presence of CSF when compared to the bare state. The duration of deformation was generally increased in the presence of CSF, though this difference was not statistically significant. This may indicate a reduction in the cord-fragment interaction force for a given impulse. The dura was found to have no significant effect on deformation parameters for the bovine spinal cord. The deformation of the surrogate cord gave similar trends for the different states in comparison to the bovine cord, but was significantly less than the bovine spinal cord for all conditions.

CONCLUSION

The results indicate that the protective mechanism of CSF may not eliminate cord deformationunder the high energy transverse impact characteristic of a burst fracture. However, CSF may contribute to a lessening of cord deformation and applied force.