The distribution of tissue damage in the spinal cord is influenced by the contusion velocity

Residual volume of extruded disc material and residual spinal cord compression measured on postoperative MRI do not predict neurological outcomes in dogs following decompressive surgery for thoracolumbar intervertebral disc extrusion

Published studies on the validity of using quantitative MRI measures of pre- and postoperative spinal cord (SC) compression as prognostic indicators for dogs undergoing surgery for intervertebral disc extrusion (IVDE) are currently limited. The aim of this retrospective analytical study was to describe the volume of postoperative residual extradural material (VREM) and the ratio of the cross-sectional area (CSA) of maximum SC compression to the CSA of SC in a compression-free intervertebral space as MRI measures of preoperative and postoperative compression (residual spinal cord compression, RSCC), and to compare these measures between the neurological outcome in a group of dogs. Inclusion criteria were dogs that underwent surgery for thoracolumbar IVDE, were imaged pre- and immediately postoperatively by MRI, and had a neurological follow-up examination 2 to 5 weeks postoperatively. Two blinded observers independently performed measurements in pre- and postoperative MRI studies. Dogs were classified into positive outcome (PO) and negative outcome (NO) groups based on follow-up neurologic examination scores. Seventeen dogs were included (12 PO, 5 NO). Interobserver agreement for MRI measurements was good to excellent (ICCs: 0.76-0.97). The prevalence of residual extradural material in postoperative MRI studies was 100%. No significant differences in mean preoperative SC compression, mean RSCC, mean SC decompression, or VREM were found between outcome groups (P = .25; P = .28; P = .91, P = .98). In conclusion, neither postoperative VREM nor RSCC could predict successful neurological outcomes.

The biomechanical implications of neck position in cervical contusion animal models of SCI

Large animal contusion models of spinal cord injury are an essential precursor to developing and evaluating treatment options for human spinal cord injury. Reducing variability in these experiments has been a recent focus as it increases the sensitivity with which treatment effects can be detected while simultaneously decreasing the number of animals required in a study. Here, we conducted a detailed review to explore if head and neck positioning in a cervical contusion model of spinal cord injury could be a factor impacting the biomechanics of a spinal cord injury, and thus, the resulting outcomes. By reviewing existing literature, we found evidence that animal head/neck positioning affects the exposed level of the spinal cord, morphology of the spinal cord, tissue mechanics and as a result the biomechanics of a cervical spinal cord injury. We posited that neck position could be a hidden factor contributing to variability. Our results indicate that neck positioning is an important factor in studying biomechanics, and that reporting these values can improve inter-study consistency and comparability and that further work needs to be done to standardize positioning for cervical spinal cord contusion injury models.

Dynamic changes in mechanical properties of the adult rat spinal cord after injury

Spinal cord injury (SCI), a debilitating medical condition that can cause irreversible loss of neurons and permanent paralysis, currently has no cure. However, regenerative medicine may offer a promising treatment. Given that numerous regenerative strategies aim to deliver cells and materials in the form of tissue-engineered therapies, understanding and characterising the mechanical properties of the spinal cord tissue is very important. In this study, we have systematically characterised the spatiotemporal changes in elastic stiffness (elastic modulus, Pa) and viscosity (drop in peak force, %) of injured rat thoracic spinal cord tissues at distinct time points after crush injury using the indentation technique. Our results demonstrate that in comparison with uninjured spinal cord tissue, the injured tissues exhibited lower stiffness (median 3281 Pa versus 9632 Pa; P < 0.001) but demonstrated elevated viscosity (median 80% versus 57%; P < 0.001) at 3 days postinjury. Between 4 and 6 weeks after SCI, the overall viscoelastic properties of injured tissues returned to baseline values. At 12 weeks after SCI, in comparison with uninjured tissue, the injured spinal cord tissues displayed a significant increase in both elasticity (median 13698 Pa versus 9920 Pa; P < 0.001) and viscosity (median 64% versus 58%; P < 0.001). This work constitutes the first quantitative mapping of spatiotemporal changes in spinal cord tissue elasticity and viscosity in injured rats, providing a mechanical basis of the tissue for future studies on the development of biomaterials for SCI repair. STATEMENT OF SIGNIFICANCE: Spinal cord injury (SCI) is a devastating disease often leading to permanent paralysis. While enormous progress in understanding the molecular pathomechanisms of SCI has been made, the mechanical properties of injured spinal cord tissue have received considerably less attention. This study provides systematic characterization of the biomechanical evolution of rat spinal cord tissue after SCI using a microindentation test method. We find spinal cord tissue behaves significantly softer but more viscous immediately postinjury. As time passes, the lesion site gradually returns to baseline values and then displays pronounced increased viscoelastic properties. As host tissue mechanical properties are a crucial consideration for any biomaterial implanted into central nervous system, our results may have important implications for further studies of SCI repair.

Quantifying Intraparenchymal Hemorrhage after Traumatic Spinal Cord Injury: A Review of Methodology

Intraparenchymal hemorrhage (IPH) after a traumatic injury has been associated with poor neurological outcomes. Although IPH may result from the initial mechanical trauma, the blood and its breakdown products have potentially deleterious effects. Further, the degree of IPH has been correlated with injury severity and the extent of subsequent recovery. Therefore, accurate evaluation and quantification of IPH following traumatic spinal cord injury (SCI) is important to define treatments' effects on IPH progression and secondary neuronal injury. Imaging modalities, such as magnetic resonance imaging (MRI) and ultrasound (US), have been explored by researchers for the detection and quantification of IPH following SCI. Both quantitative and semiquantitative MRI and US measurements have been applied to objectively assess IPH following SCI, but the optimal methods for doing so are not well established. Studies in animal SCI models (rodent and porcine) have explored US and histological techniques in evaluating SCI and have demonstrated the potential to detect and quantify IPH. Newer techniques using machine learning algorithms (such as convolutional neural networks [CNN]) have also been studied to calculate IPH volume and have yielded promising results. Despite long-standing recognition of the potential pathological significance of IPH within the spinal cord, quantifying IPH with MRI or US is a relatively new area of research. Further studies are warranted to investigate their potential use. Here, we review the different and emerging quantitative MRI, US, and histological approaches used to detect and quantify IPH following SCI.

A Biomimetic Nonwoven-Reinforced Hydrogel for Spinal Cord Injury Repair

In clinical trials, new scaffolds for regeneration after spinal cord injury (SCI) should reflect the importance of a mechanically optimised, hydrated environment. Composite scaffolds of nonwovens, self-assembling peptides (SAPs) and hydrogels offer the ability to mimic native spinal cord tissue, promote aligned tissue regeneration and tailor mechanical properties. This work studies the effects of an aligned electrospun nonwoven of P11-8—enriched poly(ε-caprolactone) (PCL) fibres, integrated with a photo-crosslinked hydrogel of glycidylmethacrylated collagen (collagen-GMA), on neurite extension. Mechanical properties of collagen-GMA hydrogel in compression and shear were recorded, along with cell viability. Collagen-GMA hydrogels showed J-shaped stress–strain curves in compression, mimicking native spinal cord tissue. For hydrogels prepared with a 0.8-1.1 wt.% collagen-GMA concentration, strain at break values were 68 ± 1–81 ± 1% (±SE); maximum stress values were 128 ± 9–311 ± 18 kPa (±SE); and maximum force values were 1.0 ± 0.1–2.5 ± 0.1 N (±SE). These values closely mimicked the compression values for feline and porcine tissue in the literature, especially those for 0.8 wt.%. Complex shear modulus values fell in the range 345–2588 Pa, with the lower modulus hydrogels in the range optimal for neural cell survival and growth. Collagen-GMA hydrogel provided an environment for homogenous and three-dimensional cell encapsulation, and high cell viability of 84 ± 2%. Combination of the aligned PCL/P11-8 electrospun nonwoven and collagen-GMA hydrogel retained fibre alignment and pore structure, respectively, and promoted aligned neurite extension of PC12 cells. Thus, it is possible to conclude that scaffolds with mechanical properties that both closely mimic native spinal cord tissue and are optimal for neural cells can be produced, which also promote aligned tissue regeneration when the benefits of hydrogels and electrospun nonwovens are combined.

Pathophysiology, Classification and Comorbidities after Traumatic Spinal Cord Injury

The spinal cord is a conduit within the central nervous system (CNS) that provides ongoing communication between the brain and the rest of the body, conveying complex sensory and motor information necessary for safety, movement, reflexes, and optimization of autonomic function. After a spinal cord injury (SCI), supraspinal influences on the spinal segmental control system and autonomic nervous system (ANS) are disrupted, leading to spastic paralysis, pain and dysesthesia, sympathetic blunting and parasympathetic dominance resulting in cardiac dysrhythmias, systemic hypotension, bronchoconstriction, copious respiratory secretions and uncontrolled bowel, bladder, and sexual dysfunction. This article outlines the pathophysiology of traumatic SCI, current and emerging methods of classification, and its influence on sensory/motor function, and introduces the probable comorbidities associated with SCI that will be discussed in more detail in the accompanying manuscripts of this special issue.

"White cord syndrome after cervical or thoracic spinal cord decompression. Haemodynamic complication or mechanical damage? An understimated nosographic entity"

The ischemia-reperfusion mechanism is believed to be responsible for parenchymal damage caused by temporary hypoperfusion and worsened by the subsequent attempt of reperfusion. This represents a true challenge for physicians of several fields, including neurosurgeons. A limited number of papers have shed the light on a rare pathological condition that affects patients experiencing an unexplained neurological deficit after spine surgery, the so-called "white cord syndrome". This entity is believed to be caused by an "ischemia-reperfusion" injury on the spinal cord, documented by a post-operative intramedullary hyperintensity on T2 weighted MRI sequences. To date, the cases of white cord syndrome reported in literature mostly refer to cervical spine surgery. However, the analysis of several reviews focusing on spine surgery outcome suggest that post-operative neurological deficits of new onset could be charged to a mechanism of ischemia-reperfusion, even if the physiopathology of this event is seldom explored or at least discussed. The same neuroradiological finding can suggest a mechanical damage due to surgical inappropriate manipulation. On this purpose, we performed a systematic revision of literature with the aim to identify and analyze all the factors potentially contributing to ischemic-reperfusion damage of the spinal cord that may potentially complicate any spinal surgery, without distinction between cervical or thoracic segment. Finally, we believe that post-operative neurological deficit after spinal surgery constituting the "white cord syndrome", could be underreported, while both neurosurgeons and patients should be fully aware of this rare but potentially devasting complication burdening cervical and thoracic spine surgery.

Comparison of numerical methods for cerebrospinal fluid representation and fluid-structure interaction during transverse impact of a finite element spinal cord model

Spinal cord impacts can have devastating consequences. Computational models can investigate such impacts but require biofidelic numerical representations of the neural tissues and fluid-structure interaction with cerebrospinal fluid. Achieving this biofidelity is challenging, particularly for efficient implementation of the cerebrospinal fluid in full computational human body models. The goal of this study was to assess the biofidelity and computational efficiency of fluid-structure interaction methods representing the cerebrospinal fluid interacting with the spinal cord, dura, and pia mater using experimental pellet impact test data from bovine spinal cords. Building on an existing finite element model of the spinal cord and pia mater, an orthotropic hyperelastic constitutive model was proposed for the dura mater and fit to literature data. The dura mater and cerebrospinal fluid were integrated with the existing finite element model to assess four fluid-structure interaction methods under transverse impact: Lagrange, pressurized volume, smoothed particle hydrodynamics, and arbitrary Lagrangian-Eulerian. The Lagrange method resulted in an overly stiff mechanical response, whereas the pressurized volume method over-predicted compression of the neural tissues. Both the smoothed particle hydrodynamics and arbitrary Lagrangian-Eulerian methods were able to effectively model the impact response of the pellet on the dura mater, outflow of the cerebrospinal fluid, and compression of the spinal cord; however, the arbitrary Lagrangian-Eulerian compute time was approximately five times higher than smoothed particle hydrodynamics. Crucial to implementation in human body models, the smoothed particle hydrodynamics method provided a computationally efficient and representative approach to model spinal cord fluid-structure interaction during transverse impact.

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.

Motor, sensitive, and vegetative recovery in rats with compressive spinal-cord injury after combined treatment with erythropoietin and whole-body vibration

BACKGROUND

Physical therapy with whole body vibration (WBV) following compressive spinal cord injury (SCI) in rats restores density of perisomatic synapses, improves body weight support and leads to a better bladder function. The purpose of the study was to determine whether the combined treatment with WBV plus erythropoietin (EPO) would further improve motor, sensory and vegetative functions after SCI in rats.

METHODS

Severe compressive SCI at low thoracic level was followed by a single i.p. injection of 2,5μg (250 IU) human recombinant EPO. Physical therapy with WBV started on 14th day after injury and continued over a 12-week post injury period. Locomotor recovery, sensitivity tests and urinary bladder scores were analysed at 1, 3, 6, 9, and 12 weeks after SCI. The closing morphological measurements included lesion volume and numbers of axons in the preserved perilesional neural tissue bridges (PNTB).

RESULTS

Assessment of motor performance sensitivity and bladder function revealed no significant effects of EPO when compared to the control treatments. EPO treatment neither reduced the lesion volume, nor increased the number of axons in PNTB.

CONCLUSIONS

The combination of WBV + EPO exerts no positive effects on hind limbs motor performance and bladder function after compressive SCI in rats.

Numbers of Axons in Spared Neural Tissue Bridges But Not Their Widths or Areas Correlate With Functional Recovery in Spinal Cord-Injured Rats

The relationships between various parameters of tissue damage and subsequent functional recovery after spinal cord injury (SCI) are not well understood. Patients may regain micturition control and walking despite large postinjury medullar cavities. The objective of this study was to establish possible correlations between morphological findings and degree of functional recovery after spinal cord compression at vertebra Th8 in rats. Recovery of motor (Basso, Beattie, Bresnahan, foot-stepping angle, rump-height index, and ladder climbing), sensory (withdrawal latency), and bladder functions was analyzed at 1, 3, 6, 9, and 12 weeks post-SCI. Following perfusion fixation, spinal cord tissue encompassing the injury site was cut in longitudinal frontal sections. Lesion lengths, lesion volumes, and areas of perilesional neural tissue bridges were determined after staining with cresyl violet. The numbers of axons in these bridges were quantified after staining for class III β-tubulin. We found that it was not the area of the spared tissue bridges, which is routinely determined by magnetic resonance imaging (MRI), but the numbers of axons in them that correlated with functional recovery after SCI (Spearman's ρ > 0.8; p < 0.001). We conclude that prognostic statements based only on MRI measurements should be considered with caution.

Mechanical properties of spinal cord grey matter and white matter in confined compression

To better understand the link between spinal cord impact and the resulting tissue damage, computational models are often used. These models typically simulate the spinal cord as a homogeneous and isotropic material. Recent research suggests that grey and white matter tissue differences and directional differences, i.e. anisotropy, are important to predict spinal cord damage. The objective of this research was to characterize the mechanical properties of spinal cord grey and white matter tissue in confined compression. Spinal cords (n = 12) were harvested immediately following euthanasia from Yorkshire and Yucatan pigs. The spinal cords were flash frozen (60 s at -80 °C) and prepared into four types of test samples: grey matter axial, grey matter transverse, white matter axial, and white matter transverse. Each sample type was thawed, and subsequently tested in confined compression within 6 h of euthanasia. Samples were compressed to 10% strain at a quasi-static strain rate (0.001/sec) and allowed to relax for 120 s. A quasi-linear viscoelastic model combining a first-order exponential with a 1-term Prony series characterized the loading and relaxation responses respectively. The effect of tissue type (grey matter vs. white matter), direction (axial vs. transverse), and their interaction were evaluated with a two-way ANOVA (p < 0.05) with peak stress, aggregate modulus, and relaxation time as dependent variables. This study found grey matter to be 1.6-2 times stiffer than white matter and both grey and white matter were isotropic in compression. These findings should be emphasized when studying SCI biomechanics using computational models.

Effect of experimental, morphological and mechanical factors on the murine spinal cord subjected to transverse contusion: A finite element study

Finite element models combined with animal experimental models of spinal cord injury provides the opportunity for investigating the effects of the injury mechanism on the neural tissue deformation and the resulting tissue damage. Thus, we developed a finite element model of the mouse cervical spinal cord in order to investigate the effect of morphological, experimental and mechanical factors on the spinal cord mechanical behavior subjected to transverse contusion. The overall mechanical behavior of the model was validated with experimental data of unilateral cervical contusion in mice. The effects of the spinal cord material properties, diameter and curvature, and of the impactor position and inclination on the strain distribution were investigated in 8 spinal cord anatomical regions of interest for 98 configurations of the model. Pareto analysis revealed that the material properties had a significant effect (p<0.01) for all regions of interest of the spinal cord and was the most influential factor for 7 out of 8 regions. This highlighted the need for comprehensive mechanical characterization of the gray and white matter in order to develop effective models capable of predicting tissue deformation during spinal cord injuries.

Effect of Velocity and Duration of Residual Compression in a Rat Dislocation Spinal Cord Injury Model

Early decompression of the traumatically injured and persistently compressed spinal cord is intuitively beneficial for neurological outcome. Despite considerable pre-clinical evidence of a neurological benefit to early decompression, the effect of early surgical decompression in clinical spinal cord injury (SCI) remains less clear. The discrepancy between pre-clinical and clinical results may be due to differences between the biomechanical variables used in pre-clinical animal models and the biomechanical conditions occurring in clinical injuries. These pre-clinical variables include region of spinal cord, velocity of impact, and injury mechanism. In this study, the effect of velocity and duration of residual compression on injury severity were evaluated using a novel, rodent model of cervical dislocation SCI. Fifty-two male Sprague-Dawley rats were included in five groups: two timings of decompression (24 min, 240 min), two velocities (10 mm/sec, 500 mm/sec), and a sham group. All injuries involved a 1.45-mm dorsal dislocation of the C6 vertebra relative to C5 with subsequent residual compression of 0.8 mm. Animals were evaluated for motor function using the Martinez open field, grip strength, and grooming tests for 6 weeks post-injury. Immunohistochemistry and histology following sacrifice were conducted with counts for NeuN- and choline acetyltransferase (ChAT)-positive neurons, and length of cavitation. Behavioral testing and histological analysis revealed that injuries induced by the high velocity were consistently more severe than those induced by the low velocity, with behavioral correlations ranging between 0.46 and 0.58 (p < 0.05). Longer duration of residual compression did not produce significantly more severe injuries as measured by functional tests and histology. These findings demonstrate that the velocity of the initial traumatic impact may be a more important factor than duration of residual compression in determining SCI severity in a dislocation model of SCI.

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.

Development of a traumatic cervical dislocation spinal cord injury model with residual compression in the rat

BACKGROUND

Preclinical spinal cord injury models do not represent the wide range of biomechanical factors seen in human injuries, such as spinal level, injury mechanism, velocity of spinal cord impact, and residual compression. These factors may be responsible for differences observed between experimental and clinical study results, especially related to the controversial issue of timing of surgical decompression.

NEW METHOD

Somatosensory Evoked Potentials were used to: a) characterize residual compression depths in a dislocation model, and b) evaluate the physiological effect of whether or not the spinal cord was decompressed following the initial injury, prior to the application of residual compression. Modifications to vertebral clamps and the development of a novel surgical frame allowed us to conduct surgical and injury procedures in a controlled manner without the risk of additional damage to the spinal cord. Behavioural outcomes were evaluated following varying dislocation displacements, in addition to the survivability of 4 h of residual compression following a traumatic injury.

RESULTS

Residual compression immediately following the initial dislocation demonstrated significantly different electrophysiological response compared to when the residual compression was delayed.

COMPARISON WITH EXISTING METHOD

There are currently no other residual compression models that utilize a dislocation injury mechanism. Many residual compression studies have demonstrated the effectiveness of early decompression, however the compression of the spinal cord is often not representative of clinical traumatic injuries. Preclinical studies typically model residual compression using a sustained force through quasi-static application, when human injuries often occur at high velocities, followed by a sustained displacement occlusion of the spinal canal.

CONCLUSIONS

This study has validated several novel procedural approaches and injury parameters, and provided critical details to implement in the development of a traumatic cervical dislocation SCI model with residual compression.

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.

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.

Radiography used to measure internal spinal cord deformation in an in vivo rat model

Little is known about the internal mechanics of the in vivo spinal cord during injury. The objective of this study was to develop a method of tracking internal and surface deformation of in vivo rat spinal cord during compression using radiography. Since neural tissue is radio-translucent, radio-opaque markers were injected into the spinal cord. Two tantalum beads (260 µm) were injected into the cord (dorsal and ventral) at C5 of nine anesthetized rats. Four beads were glued to the lateral surface of the cord, caudal and cranial to the injection site. A compression plate was displaced 0.5 mm, 2 mm, and 3 mm into the spinal cord and lateral X-ray images were taken before, during, and after each compression for measuring bead displacements. Potential bead migration was monitored for by comparing displacements of the internal and glued surface beads. Dorsal beads moved significantly more than ventral beads with a range in averages of 0.57-0.71 mm and 0.31-0.35 mm respectively. Bead displacements during 0.5 mm compressions were significantly lower than 2 mm and 3 mm compressions. There was no statistically significant migration of the internal beads. The results indicate the merit of this technique for measuring in vivo spinal cord deformation. The pattern of bead displacements illustrates the complex internal and surface deformations of the spinal cord during transverse compression. This information is needed for validating physical and finite element spinal cord surrogates and to define relationships between loading parameters, internal cord deformation, and biological and functional outcomes.

The Transverse Isotropy of Spinal Cord White Matter Under Dynamic Load

The rostral-caudally aligned fiber-reinforced structure of spinal cord white matter (WM) gives rise to transverse isotropy in the material. Stress and strain patterns generated in the spinal cord parenchyma following spinal cord injury (SCI) are multidirectional and dependent on the mechanism of the injury. Our objective was to develop a WM constitutive model that captures the material transverse isotropy under dynamic loading. The WM mechanical behavior was extracted from the published tensile and compressive experiments. Combinations of isotropic and fiber-reinforcing models were examined in a conditional quasi-linear viscoelastic (QLV) formulation to capture the WM mechanical behavior. The effect of WM transverse isotropy on SCI model outcomes was evaluated by simulating a nonhuman primate (NHP) contusion injury experiment. A second-order reduced polynomial hyperelastic energy potential conditionally combined with a quadratic reinforcing function in a four-term Prony series QLV model best captured the WM mechanical behavior (0.89 < R2 < 0.99). WM isotropic and transversely isotropic material models combined with discrete modeling of the pia mater resulted in peak impact forces that matched the experimental outcomes. The transversely isotropic WM with discrete pia mater resulted in maximum principal strain (MPS) distributions which effectively captured the combination of ipsilateral peripheral WM sparing, ipsilateral injury and contralateral sparing, and the rostral/caudal spread of damage observed in in vivo injuries. The results suggest that the WM transverse isotropy could have an important role in correlating tissue damage with mechanical measures and explaining the directional sensitivity of the spinal cord to injury.

Mechanical Design and Analysis of a Unilateral Cervical Spinal Cord Contusion Injury Model in Non-Human Primates

Non-human primate (NHP) models of spinal cord injury better reflect human injury and provide a better foundation to evaluate potential treatments and functional outcomes. We combined finite element (FE) and surrogate models with impact data derived from in vivo experiments to define the impact mechanics needed to generate a moderate severity unilateral cervical contusion injury in NHPs (Macaca mulatta). Three independent variables (impactor displacement, alignment, and pre-load) were examined to determine their effects on tissue level stresses and strains. Mechanical measures of peak force, peak displacement, peak energy, and tissue stiffness were analyzed as potential determinants of injury severity. Data generated from FE simulations predicted a lateral shift of the spinal cord at high levels of compression (>64%) during impact. Submillimeter changes in mediolateral impactor position over the midline increased peak impact forces (>50%). Surrogate cords established a 0.5 N pre-load protocol for positioning the impactor tip onto the dural surface to define a consistent dorsoventral baseline position before impact, which corresponded with cerebrospinal fluid displacement and entrapment of the spinal cord against the vertebral canal. Based on our simulations, impactor alignment and pre-load were strong contributors to the variable mechanical and functional outcomes observed in in vivo experiments. Peak displacement of 4 mm after a 0.5N pre-load aligned 0.5-1.0 mm over the midline should result in a moderate severity injury; however, the observed peak force and calculated peak energy and tissue stiffness are required to properly characterize the severity and variability of in vivo NHP contusion injuries.

A Unilateral Cervical Spinal Cord Contusion Injury Model in Non-Human Primates (Macaca mulatta)

The development of a non-human primate (NHP) model of spinal cord injury (SCI) based on mechanical and computational modeling is described. We scaled up from a rodent model to a larger primate model using a highly controllable, friction-free, electronically-driven actuator to generate unilateral C6-C7 spinal cord injuries. Graded contusion lesions with varying degrees of functional recovery, depending upon pre-set impact parameters, were produced in nine NHPs. Protocols and pre-operative magnetic resonance imaging (MRI) were used to optimize the predictability of outcomes by matching impact protocols to the size of each animal's spinal canal, cord, and cerebrospinal fluid space. Post-operative MRI confirmed lesion placement and provided information on lesion volume and spread for comparison with histological measures. We evaluated the relationships between impact parameters, lesion measures, and behavioral outcomes, and confirmed that these relationships were consistent with our previous studies in the rat. In addition to providing multiple univariate outcome measures, we also developed an integrated outcome metric describing the multivariate cervical SCI syndrome. Impacts at the higher ranges of peak force produced highly lateralized and enduring deficits in multiple measures of forelimb and hand function, while lower energy impacts produced early weakness followed by substantial recovery but enduring deficits in fine digital control (e.g., pincer grasp). This model provides a clinically relevant system in which to evaluate the safety and, potentially, the efficacy of candidate translational therapies.

Influence of sagittal and axial types of ossification of posterior longitudinal ligament on mechanical stress in cervical spinal cord: A finite element analysis

BACKGROUND

There are few studies focusing on the prediction of stress distribution according to the types of ossification of the posterior longitudinal ligament, which can be fundamental information associated with clinical aspects such as the relationship between stress level and neurological symptom severity. In this study, the influence of sagittal and axial types of ossification of the posterior longitudinal ligament on mechanical stress in the cervical spinal cord was investigated.

METHODS

A three-dimensional finite element model of the cervical spine with spinal cord was developed and validated. The von Mises stresses in the cord and the reduction in cross-sectional areas and volume of the cord were investigated for various axial and sagittal types according to the occupying ratio of ossification of the posterior longitudinal ligament in the spinal canal.

FINDINGS

The influence of axial type was less than that of the sagittal type, even though the central type showed higher maximum stresses in the cord, especially for the continuous type. With a 60% occupying ratio of ossification of the posterior longitudinal ligament, the maximum stress was significantly high and the cross-sectional area of the spinal cord was reduced by more than 30% of the intact area regardless of sagittal or axial types. Finally, a higher level of sagittal extension would increase the peak cord tissue stress, which would be related to the neurological dysfunction and tissue damage.

INTERPRETATION

Quantitative investigation of biomechanical characteristics such as mechanical stress may provide fundamental information for pre-operative planning of treatment for ossification of the posterior longitudinal ligament.

Impact depth and the interaction with impact speed affect the severity of contusion spinal cord injury in rats

Spinal cord injury (SCI) biomechanics suggest that the mechanical factors of impact depth and speed affect the severity of contusion injury, but their interaction is not well understood. The primary aim of this work was to examine both the individual and combined effects of impact depth and speed in contusion SCI on the cervical spinal cord. Spinal cord contusions between C5 and C6 were produced in anesthetized rats at impact speeds of 8, 80, or 800 mm/s with displacements of 0.9 or 1.5 mm (n=8/group). After 7 days postinjury, rats were assessed for open-field behavior, euthanized, and spinal cords were harvested. Spinal cord tissue sections were stained for demyelination (myelin-based protein) and tissue sparing (Luxol fast blue). In parallel, a finite element model of rat spinal cord was used to examine the resulting maximum principal strain in the spinal cord during impact. Increasing impact depth from 0.9 to 1.5 mm reduced open-field scores (p<0.01) above 80 mm/s, reduced gray (GM) and white matter (WM) sparing (p<0.01), and increased the amount of demyelination (p<0.01). Increasing impact speed showed similar results at the 1.5-mm impact depth, but not the 0.9-mm impact depth. Linear correlation analysis with finite element analysis strain showed correlations (p<0.001) with nerve fiber damage in the ventral (R(2)=0.86) and lateral (R(2)=0.74) regions of the spinal cord and with WM (R(2)=0.90) and GM (R(2)=0.76) sparing. The results demonstrate that impact depth is more important in determining the severity of SCI and that threshold interactions exist between impact depth and speed.

Novel nanometer scaffolds regulate the biological behaviors of neural stem cells

Ideal tissue-engineered scaffold materials regulate proliferation, apoptosis and differentiation of cells seeded on them by regulating gene expression. In this study, aligned and randomly oriented collagen nanofiber scaffolds were prepared using electronic spinning technology. Their diameters and appearance reached the standards of tissue-engineered nanometer scaffolds. The nanofiber scaffolds were characterized by a high swelling ratio, high porosity and good mechanical properties. The proliferation of spinal cord-derived neural stem cells on novel nanofiber scaffolds was obviously enhanced. The proportions of cells in the S and G2/M phases noticeably increased. Moreover, the proliferation rate of neural stem cells on the aligned collagen nanofiber scaffolds was high. The expression levels of cyclin D1 and cyclin-dependent kinase 2 were increased. Bcl-2 expression was significantly increased, but Bax and caspase-3 gene expressions were obviously decreased. There was no significant difference in the differentiation of neural stem cells into neurons on aligned and randomly oriented collagen nanofiber scaffolds. These results indicate that novel nanofiber scaffolds could promote the proliferation of spinal cord-derived neural stem cells and inhibit apoptosis without inducing differentiation. Nanofiber scaffolds regulate apoptosis and proliferation in neural stem cells by altering gene expression.

Characterization of a Novel, Magnetic Resonance Imaging-Compatible Rodent Model Spinal Cord Injury Device

Rodent models of acute spinal cord injury (SCI) are often used to investigate the effects of injury mechanism, injury speed, and cord displacement magnitude, on the ensuing cascade of biological damage in the cord. However, due to its small size, experimental observations have largely been limited to the gross response of the cord. To properly understand the relationship between mechanical stimulus and biological damage, more information is needed about how the constituent tissues of the cord (i.e., gray and white matter) respond to injurious stimuli. To address this limitation, we developed a novel magnetic resonance imaging (MRI)-compatible test apparatus that can impose either a contusion-type or dislocation-type acute cervical SCI in a rodent model and facilitate MR-imaging of the cervical spinal cord in a 7 T MR scanner. In this study, we present the experimental performance parameters of the MR rig. Utilizing cadaveric specimens and static radiographs, we report contusion magnitude accuracy that for a desired 1.8 mm injury, a nominal 1.78 mm injury (SD = 0.12 mm) was achieved. High-speed video analysis was employed to determine the injury speeds for both mechanisms and were found to be 1147 mm/s (SD = 240 mm/s) and 184 mm/s (SD = 101 mm/s) for contusion and dislocation injuries, respectively. Furthermore, we present qualitative pilot data from a cadaveric trial, employing the MR rig, to show the expected results from future studies.

Severity of spinal cord injury in adult and infant rats after vertebral dislocation depends upon displacement but not speed

Spinal cord injury (SCI) is less common in children than in adults, but in children it is generally more severe. Spinal loading conditions (speed and displacement) are also thought to affect SCI severity, but the relationship between these parameters is not well understood. This study aimed to investigate the effects of vertebral speed and displacement on the severity of SCI in infants and adults using a rodent model of vertebral dislocation. Thoracolumbar vertebral dislocation was induced in anaesthetized infant rats (∼30 g, 13-15 days postnatal, n=40) and adult rats (∼250 g, n=57). The 12th thoracic vertebra was secured, whereas the first lumbar vertebra was dislocated laterally. Dislocation speed and magnitude were varied independently and scaled between adults and infants (Adults: 100-250mm/s, 4-10mm; Infants: 40-100mm/s, 1.6-4mm). At 5 h post-injury, rats were euthanized and spinal cords harvested. Spinal cord sections were stained to detect hemorrhage (hematoxylin and eosin) and axonal injury (β-amyloid precursor protein). For each millimeter increase in vertebral displacement, normalized hemorrhage volume increased by 1.9×10(-3) mm(3) (p=0.028) and normalized area of axonal injury increased by 2.2×10(-1)mm(2) (p<0.001). Normalized hemorrhage volume was 3.3×10(-3) mm(3) greater for infants than for adults (p<0.001). Magnitude of dislocation was found to have a different effect on the normalized area of axonal injury in adults than in infants (p=0.003). Speed of dislocation was not found to have a significant effect on normalized hemorrhage volume (p=0.427) or normalized area of axonal injury (p=0.726) independent of displacement for the range of speeds tested. The findings of this study suggest that both age and amount of spinal motion are key factors in the severity of acute SCI.

Development of surrogate spinal cords for the evaluation of electrode arrays used in intraspinal implants

We report the development of a surrogate spinal cord for evaluating the mechanical suitability of electrode arrays for intraspinal implants. The mechanical and interfacial properties of candidate materials (including silicone elastomers and gelatin hydrogels) for the surrogate cord were tested. The elastic modulus was characterized using dynamic mechanical analysis, and compared with values of actual human spinal cords from the literature. Forces required to indent the surrogate cords to specified depths were measured to obtain values under static conditions. Importantly, to quantify surface properties in addition to mechanical properties normally considered, interfacial frictional forces were measured by pulling a needle out of each cord at a controlled rate. The measured forces were then compared to those obtained from rat spinal cords. Formaldehyde-crosslinked gelatin, 12 wt% in water, was identified as the most suitable material for the construction of surrogate spinal cords. To demonstrate the utility of surrogate spinal cords in evaluating the behavior of various electrode arrays, cords were implanted with two types of intraspinal electrode arrays (one made of individual microwires and another of microwires anchored with a solid base), and cord deformation under elongation was evaluated. The results demonstrate that the surrogate model simulates the mechanical and interfacial properties of the spinal cord, and enables in vitro screening of intraspinal implants.

Spinal cord injury with unilateral versus bilateral primary hemorrhage--effects of glibenclamide

In spinal cord injury (SCI), block of Sur1-regulated NC(Ca-ATP) channels by glibenclamide protects penumbral capillaries from delayed fragmentation, resulting in reduced secondary hemorrhage, smaller lesions and better neurological function. All published experiments demonstrating a beneficial effect of glibenclamide in rat models of SCI have used a cervical hemicord impact calibrated to produce primary hemorrhage located exclusively ipsilateral to the site of impact. Here, we tested the hypothesis that glibenclamide also would be protective in a model with more extensive, bilateral primary hemorrhage. We studied the effect of glibenclamide in 2 rat cervical hemicord contusion models with identical impact force (10 g, 25 mm), one with the impactor positioned laterally to yield unilateral primary hemorrhage (UPH), and the other with the impactor positioned more medially, yielding larger, bilateral primary hemorrhages (BPH) and 6-week lesion volumes that were 45% larger. Functional outcome measures included: modified (unilateral) Basso, Beattie, and Bresnahan scores, angled plane performance, and rearing times. In the UPH model, the effects of glibenclamide were similar to previous observations, including a functional benefit as early as 24h after injury and 6-week lesion volumes that were 57% smaller than controls. In the BPH model, glibenclamide exerted a significant benefit over controls, but the functional benefit was smaller than in the UPH model and 6-week lesion volumes were 33% smaller than controls. We conclude that glibenclamide is beneficial in different models of cervical SCI, with the magnitude of the benefit depending on the magnitude and extent of primary hemorrhage.

Independent evaluation of the effects of glibenclamide on reducing progressive hemorrhagic necrosis after cervical spinal cord injury

These experiments were completed as part of an NIH-NINDS contract entitled "Facilities of Research Excellence - Spinal Cord Injury (FORE-SCI) - Replication". Our goal was to replicate pre-clinical data from Simard et al. (2007) showing that glibenclamide, an FDA approved anti-diabetic drug that targets sulfonylurea receptor 1 (SUR1)-regulated Ca(2+) activated, [ATP](i)-sensitive nonspecific cation channels, attenuates secondary intraspinal hemorrhage and secondary neurodegeneration caused by hemicontusion injury in rat cervical spinal cord. In an initial replication attempt, the Infinite Horizons impactor was used to deliver a standard unilateral contusion injury near the spinal cord midline. Glibenclamide was administered continuously via osmotic pump beginning immediately post-SCI. The ability of glibenclamide to limit intraspinal hemorrhage was analyzed at 6, 12 and 24 h post-injury using a colorimetric assay. Acute recovery (24 h) of forelimb function was also assessed. Analysis of data from these initial studies revealed no difference between glibenclamide and vehicle-treated SCI rats. Later, it was determined that differences in primary trauma affect the efficacy of glibenclamide. Indeed, the magnitude and distribution of primary intraspinal hemorrhage was greater when the impact was directed to the dorsomedial region of the cervical hemicord (as in our initial replication experiment), as compared to the dorsolateral spinal cord (as in the Simard et al. experiment). In three subsequent experiments, injury was directed to the dorsolateral spinal cord. In each case, glibenclamide reduced post-traumatic hemorrhage 24-48 h post-injury. In the third experiment, we also assessed function and found that acute reduction of hemorrhage led to improved functional recovery. Thus, independent replication of the Simard et al. data was achieved. These data illustrate that the injury model and type of trauma can determine the efficacy of pre-clinical pharmacological treatments after SCI.

Effect of dynamic stiffness of the substrates on neurite outgrowth by using a DNA-crosslinked hydrogel

Central nervous system tissues, like other tissue types, undergo constant remodeling, which potentially leads to changes in their mechanical stiffness. Moreover, mechanical compliance of central nervous system tissues can also be modified under external load such as that experienced in traumatic brain or spinal cord injury, and during pathological processes. Thus, the neuronal responses to the dynamic stiffness of the microenvironment are of significance. In this study, we induced decrease in stiffness by using a DNA-crosslinked hydrogel, and subjected rat spinal cord neurons to such dynamic stiffness. The neurons respond to the dynamic cues as evidenced by the primary neurite structure, and the response from each neurite property (e.g., axonal length and primary dendrite number) is consistent with the behavior on static gels of same substrate rigidity, with one exception of mean primary dendrite length. The results on cell population distribution confirm the neuronal responses to the dynamic stiffness. Quantification on the focal adhesion kinase expression in the neuronal cell body on dynamic gels suggests that neurons also modify adhesion in coping with the dynamic stiffnesses. The results reported here extend the neuronal mechanosensing capability to dynamic stiffness of extracellular matrix, and give rise to a novel way of engineering neurite outgrowth in time dimension.

Modeling spinal cord contusion, dislocation, and distraction: characterization of vertebral clamps, injury severities, and node of Ranvier deformations

Spinal cord contusion and transection models are widely used for studying spinal cord injury (SCI). Clinically, however, other biomechanical injury mechanisms such as vertebral dislocation and distraction frequently occur, but these injuries are difficult to produce in animals. We mechanically characterize a vertebral clamping strategy that enables the modeling of vertebral dislocation and distraction injuries--in addition to the standard contusion paradigm--in the rat cervical spine. These vertebral clamps have a stiffness of 83.6+/-18.9 N/mm and clamping strength 64.7+/-10.2N which allows injuries to be modeled at high-speed (approximately 100 cm/s). Logistic regression indicated that a moderate-to-severe injury, with an acute mortality rate of 10%, occurs at 2.6 mm of C4/5 dorso-ventral dislocation and 4.1 mm of rostro-caudal distraction between C4 and C5. Injuries produced by dislocation and distraction exhibited features of axonal damage that were absent in contusion injuries. We conducted morphometric analysis at the nodes of Ranvier using immunohistochemistry for potassium channels (Kv1.2) in the juxtaparanodal region. Following distraction injuries, elongated nodes of Ranvier were observed up to 4mm rostral to the lesion. In contrast, contusion injuries produced distortions in nodal geometry which were restricted to the vicinity of the lesion. The greatest deformations in node of Ranvier geometry occurred at the dislocation epicenter. Given the importance of white matter damage in SCI pathology, the distinctiveness of these injury patterns demonstrate that the dislocation and distraction injury models complement existing contusion models. Together, these three animal models span a broader clinical spectrum for more reliably gauging the potential human efficacy of therapeutic strategies.

The development of an improved physical surrogate model of the human spinal cord--tension and transverse compression

To prevent spinal cord injury, optimize treatments for it, and better understand spinal cord pathologies such as spondylotic myelopathy, the interaction between the spinal column and the spinal cord during injury and pathology must be understood. The spinal cord is a complex and very soft tissue that changes properties rapidly after death and is difficult to model. Our objective was to develop a physical surrogate spinal cord with material properties closely corresponding to the in vivo human spinal cord that would be suitable for studying spinal cord injury under a variety of injurious conditions. Appropriate target material properties were identified from published studies and several candidate surrogate materials were screened, under uniaxial tension, in a materials testing machine. QM Skin 30, a silicone elastomer, was identified as the most appropriate material. Spinal cords manufactured from QM Skin 30 were tested under uniaxial tension and transverse compression. Rectangular specimens of QM Skin 30 were also tested under uniform compression. QM Skin 30 produced surrogate cords with a Young's modulus in tension and compression approximately matching values reported for in vivo animal spinal cords (0.25 and 0.20 MPa, respectively). The tensile and compressive Young's modulus and the behavior of the surrogate cord simulated the nonlinear behavior of the in vivo spinal cord.