Investigating the active contractile function of the rat paraspinal muscles reveals unique cross-bridge kinetics in the multifidus.

PURPOSE: Various aspects of paraspinal muscle anatomy, biology, and histology have been studied; however, information on paraspinal muscle contractile function is almost nonexistent, thus hindering functional interpretation of these muscles in healthy individuals and those with low back disorders. The aim of this study was to measure and compare the contractile function and force-sarcomere length properties of muscle fibers from the multifidus (MULT) and erector spinae (ES) as well as a commonly studied lower limb muscle (Extensor digitorum longus (EDL)) in the rat. METHODS: Single muscle fibers (n = 77 total from 6 animals) were isolated from each of the muscles and tested to determine their active contractile function; all fibers used in the analyses were type IIB. RESULTS: There were no significant differences between muscles for specific force (sFo) (p = 0.11), active modulus (p = 0.63), average optimal sarcomere length (p = 0.27) or unloaded shortening velocity (Vo) (p = 0.69). However, there was a significant difference in the rate of force redevelopment (ktr) between muscles (p = < 0.0001), with MULT being significantly faster than both the EDL (p = < 0.0001) and ES (p = 0.0001) and no difference between the EDL and ES (p = 0.41). CONCLUSIONS: This finding suggests that multifidus has faster cross-bridge turnover kinetics when compared to other muscles (ES and EDL) when matched for fiber type. Whether the faster cross-bridge kinetics translate to a functionally significant difference in whole muscle performance needs to be studied further.

Dysfunctional paraspinal muscles in adult spinal deformity patients lead to increased spinal loading.

PURPOSE: Decreased spinal extensor muscle strength in adult spinal deformity (ASD) patients is well-known but poorly understood; thus, this study aimed to investigate the biomechanical and histopathological properties of paraspinal muscles from ASD patients and predict the effect of altered biomechanical properties on spine loading. METHODS: 68 muscle biopsies were collected from nine ASD patients at L4-L5 (bilateral multifidus and longissimus sampled). The biopsies were tested for muscle fiber and fiber bundle biomechanical properties and histopathology. The small sample size (due to COVID-19) precluded formal statistical analysis, but the properties were compared to literature data. Changes in spinal loading due to the measured properties were predicted by a lumbar spine musculoskeletal model. RESULTS: Single fiber passive elastic moduli were similar to literature values, but in contrast, the fiber bundle moduli exhibited a wide range beyond literature values, with 22% of 171 fiber bundles exhibiting very high elastic moduli, up to 20 times greater. Active contractile specific force was consistently less than literature, with notably 24% of samples exhibiting no contractile ability. Histological analysis of 28 biopsies revealed frequent fibro-fatty replacement with a range of muscle fiber abnormalities. Biomechanical modelling predicted that high muscle stiffness could increase the compressive loads in the spine by over 500%, particularly in flexed postures. DISCUSSION: The histopathological observations suggest diverse mechanisms of potential functional impairment. The large variations observed in muscle biomechanical properties can have a dramatic influence on spinal forces. These early findings highlight the potential key role of the paraspinal muscle in ASD.

The Effect of Compression Applied Through Constrained Lateral Eccentricity on the Failure Mechanics and Flexibility of the Human Cervical Spine.

In contrast to sagittal plane spine biomechanics, little is known about the response of the cervical spine to axial compression with lateral eccentricity of the applied force. This study evaluated the effect of lateral eccentricity on the kinetics, kinematics, canal occlusion, injuries, and flexibility of the cervical spine in translationally constrained axial impacts. Eighteen functional spinal units were subjected to flexibility tests before and after an impact. Impact axial compression was applied at one of three lateral eccentricity levels based on percentage of vertebral body width (low = 5%, medium = 50%, high = 150%). Injuries were graded by dissection. Correlations between intrinsic specimen properties and injury scores were examined for each eccentricity group. Low lateral force eccentricity produced predominantly bone injuries, clinically recognized as compression injuries, while medium and high eccentricity produced mostly contralateral ligament and/or disc injuries, an asymmetric pattern typical of lateral loading. Mean compression force at injury decreased with increasing lateral eccentricity (low = 3098 N, medium = 2337 N, and high = 683 N). Mean ipsilateral bending moments at injury were higher at medium (28.3 N·m) and high (22.9 N·m) eccentricity compared to low eccentricity specimens (0.1 N·m), p < 0.05. Ipsilateral bony injury was related to vertebral body area (VBA) (r = -0.974, p = 0.001) and disc degeneration (r = 0.851, p = 0.032) at medium eccentricity. Facet degeneration was correlated with central bony injury at high eccentricity (r = 0.834, p = 0.036). These results deepen cervical spine biomechanics knowledge in circumstances with coronal plane loads.

The effect of posture on lumbar muscle morphometry from upright MRI.

PURPOSE: To assess the effect of upright, seated, and supine postures on lumbar muscle morphometry at multiple spinal levels and for multiple muscles. METHODS: Six asymptomatic volunteers were imaged (0.5 T upright open MRI) in 7 postures (standing, standing holding 8 kg, standing 45° flexion, seated 45° flexion, seated upright, seated 45° extension, and supine), with scans at L3/L4, L4/L5, and L5/S1. Muscle cross-sectional area (CSA) and muscle position with respect to the vertebral body centroid (radius and angle) were measured for the multifidus/erector spinae combined and psoas major muscles. RESULTS: Posture significantly affected the multifidus/erector spinae CSA with decreasing CSA from straight postures (standing and supine) to seated and flexed postures (up to 19%). Psoas major CSA significantly varied with vertebral level with opposite trends due to posture at L3/L4 (increasing CSA, up to 36%) and L5/S1 (decreasing CSA, up to 40%) with sitting/flexion. For both muscle groups, radius and angle followed similar trends with decreasing radius (up to 5%) and increasing angle (up to 12%) with seated/flexed postures. CSA and lumbar lordosis had some correlation (multifidus/erector spinae L4/L5 and L5/S1, r = 0.37-0.45; PS L3/L4 left, r = - 0.51). There was generally good repeatability (average ICC(3, 1): posture = 0.81, intra = 0.89, inter = 0.82). CONCLUSION: Changes in multifidus/erector spinae muscle CSA likely represent muscles stretching between upright and seated/flexed postures. For the psoas major, the differential level effect suggests that changing three-dimensional muscle morphometry with flexion is not uniform along the muscle length. The muscle and spinal level-dependent effects of posture and spinal curvature correlation, including muscle CSA and position, highlight considering measured muscle morphometry from different postures in spine models.

Repeatability of a Dislocation Spinal Cord Injury Model in a Rat-A High-Speed Biomechanical Analysis.

Dislocation is the most common, and severe, spinal cord injury (SCI) mechanism in humans, yet there are few preclinical models. While dislocation in the rat model has been shown to produce unique outcomes, like other closed column models it exhibits higher outcome variability. Refinement of the dislocation model will enhance the testing of neuroprotective strategies, further biomechanical understanding, and guide therapeutic decisions. The overall objective of this study is to improve biomechanical repeatability of a dislocation SCI model in the rat, through the following specific aims: (i) design new injury clamps that pivot and self-align to the vertebrae; (ii) measure intervertebral kinematics during injury using the existing and redesigned clamps; and (iii) compare relative motion at the vertebrae-clamp interface to determine which clamps provide the most rigid connection. Novel clamps that pivot and self-align were developed based on the quantitative rat vertebral anatomy. A dislocation injury was produced in 34 rats at C4/C5 using either the existing or redesigned clamps, and a high-speed X-ray device recorded the kinematics. Relative motion between the caudal clamp and C5 was significantly greater in the existing clamps compared to the redesigned clamps in dorsoventral translation and sagittal rotation. This study demonstrates that relative motions can be of magnitudes that likely affect injury outcomes. We recommend such biomechanical analyses be applied to other SCI models when repeatability is an issue. For this dislocation model, the results show the importance of using clamps that pivot and self-align to the vertebrae.

Limitations of current in vitro test protocols for investigation of instrumented adjacent segment biomechanics: critical analysis of the literature.

PURPOSE: Accelerated degenerative changes at intervertebral levels adjacent to a spinal fusion, the so-called adjacent segment degeneration (ASD), have been reported in many clinical studies. Even though the pathogenesis of ASD is still widely unknown, biomechanical in vitro approaches have often been used to investigate the impact of spinal instrumentation on the adjacent segments. The goal of this review is (1) to summarize the results of these studies with respect to the applied protocol and loads and (2) to discuss if the assumptions made for the different protocols match the patients' postoperative situation. METHODS: A systematic MEDLINE search was performed using the keywords "adjacent", "in vitro" and "spine" in combination. This revealed a total of 247 articles of which 33 met the inclusion criteria. In addition, a mechanical model was developed to evaluate the effects of the current in vitro biomechanical test protocols on the changes in the adjacent segments resulting from different stiffnesses of the "treated" segment. RESULTS: The surgical treatments reported in biomechanical in vitro studies investigating ASD can be categorized into fusion procedures, total disc replacement (TDR), and dynamic implants. Three different test protocols (i.e. flexibility, stiffness, hybrid) with different loading scenarios (e.g. pure moment or eccentric load) are used in current biomechanical in vitro studies investigating ASD. According to the findings with the mechanical model, we found that the results for fusion procedures highly depend on the test protocol and method of load application, whereas for TDR and dynamic implants, most studies did not find significant changes in the adjacent segments, independent of which test protocol was used. CONCLUSIONS: The three test protocols mainly differ in the assumption on the postoperative motion behavior of the patients, which is the main reason for the conflicting findings. However, the protocols have never been validated using in vivo kinematic data. In a parallel review on in vivo kinematics by Malakoutian et al., it was found that the assumption that the patients move exactly the same after fusion implemented with the stiffness- and hybrid protocol does not match the patients' behavior. They showed that the motion of the whole lumbar spine rather tends to decrease in most studies, which could be predicted by the flexibility protocol. However, when the flexibility protocol is used with the "gold standard" pure moment, the difference in the kinematic changes between the cranial and caudal adjacent segment cannot be reproduced, putting the validity of current in vitro protocols into question.

Do in vivo kinematic studies provide insight into adjacent segment degeneration? A qualitative systematic literature review.

PURPOSE: While much evidence suggests that adjacent segment degeneration is merely a manifestation of the natural degenerative process unrelated to any spine fusion, a significant body of literature supports the notion that it is a process due in part to the altered biomechanics adjacent to fused spine segments. The purpose of this study was to review and critically analyze the published literature that investigated the in vivo kinematics of the adjacent segments and entire lumbar spine in patients receiving spinal fusion or motion-preserving devices. METHODS: A systematic review of the PubMed database was conducted, initially identifying 697 studies of which 39 addressed the in vivo kinematics of the segments adjacent to spinal implants or non-instrumented fusion of the lumbar spine. RESULTS: Twenty-nine articles studied fusion, of which three reported a decrease in range of motion of the caudal adjacent segment post-fusion. Examining the rostral adjacent segment, twelve studies observed no change, nine studies found a significant increase, and three studies reported a significant decrease in sagittal plane range of motion. Of the six studies that analyzed motion for the entire lumbar spine as a unit, five studies showed a significant decrease and one study reported no change in global lumbar spine motion. Kinematics of the segment rostral to a total disc replacement was investigated in six studies: four found no change and the results for the other two showed dependence on treatment level. Fifteen studies of non-fusion posterior implants analyzed the motion of the adjacent segment with two studies noting an increase in motion at the rostral level. CONCLUSIONS: There appears to be no overall kinematic changes at the rostral or caudal levels adjacent to a fusion, but some patients (~20-30%) develop excessive kinematic changes (i.e., instability) at the rostral adjacent level. The overall lumbar ROM after fusion appears to decrease after a spinal fusion.

Characterization of the behavior of a novel low-stiffness posterior spinal implant under anterior shear loading on a degenerative spinal model.

PURPOSE: Dynamic implants have been developed to address potential adjacent level effects due to rigid instrumentation. Rates of revision surgeries may be reduced by using improved implants in the primary surgery. Prior to clinical use, implants should be rigorously tested ex vivo. The objective of our study was to characterize the load-sharing and kinematic behavior of a novel low-stiffness spinal implant. METHODS: A human cadaveric model of degenerative spondylolisthesis was tested in shear. Lumbar functional spinal units (N = 15) were tested under a static 300 N axial compression force and a cyclic anterior shear force (5-250 N). Translation was tracked with a motion capture system. A novel implant was compared to three standard implants with shear stiffness ranging from low to high. All implants were instrumented with strain gauges to measure the supported shear force. Each implant was affixed to each specimen, and the specimens were tested intact and in two progressively destabilized states. RESULTS: Specimen condition and implant type affected implant load-sharing and specimen translation (p < 0.0001). Implant load-sharing increased across all degeneration-simulating specimen conditions and decreased across the three standard implants (high- to low-stiffness). Translation increased with the three standard implants (trend). The novel implant behaved similarly to the medium-stiffness implant (p > 0.2). CONCLUSIONS: The novel implant behaved similarly to the medium-stiffness implant in both load-sharing and translation despite having a different design and stiffness. Complex implant design and specimen-implant interaction necessitate pre-clinical testing of novel implants. Further in vitro testing in axial rotation and flexion-extension is recommended as they are highly relevant loading directions for non-rigid implants.

Moment measurements in dynamic and quasi-static spine segment testing using eccentric compression are susceptible to artifacts based on loading configuration.

The tolerance of the spine to bending moments, used for evaluation of injury prevention devices, is often determined through eccentric axial compression experiments using segments of the cadaver spine. Preliminary experiments in our laboratory demonstrated that eccentric axial compression resulted in "unexpected" (artifact) moments. The aim of this study was to evaluate the static and dynamic effects of test configuration on bending moments during eccentric axial compression typical in cadaver spine segment testing. Specific objectives were to create dynamic equilibrium equations for the loads measured inferior to the specimen, experimentally verify these equations, and compare moment responses from various test configurations using synthetic (rubber) and human cadaver specimens. The equilibrium equations were verified by performing quasi-static (5 mm/s) and dynamic experiments (0.4 m/s) on a rubber specimen and comparing calculated shear forces and bending moments to those measured using a six-axis load cell. Moment responses were compared for hinge joint, linear slider and hinge joint, and roller joint configurations tested at quasi-static and dynamic rates. Calculated shear force and bending moment curves had similar shapes to those measured. Calculated values in the first local minima differed from those measured by 3% and 15%, respectively, in the dynamic test, and these occurred within 1.5 ms of those measured. In the rubber specimen experiments, for the hinge joint (translation constrained), quasi-static and dynamic posterior eccentric compression resulted in flexion (unexpected) moments. For the slider and hinge joints and the roller joints (translation unconstrained), extension ("expected") moments were measured quasi-statically and initial flexion (unexpected) moments were measured dynamically. In the cadaver experiments with roller joints, anterior and posterior eccentricities resulted in extension moments, which were unexpected and expected, for those configurations, respectively. The unexpected moments were due to the inertia of the superior mounting structures. This study has shown that eccentric axial compression produces unexpected moments due to translation constraints at all loading rates and due to the inertia of the superior mounting structures in dynamic experiments. It may be incorrect to assume that bending moments are equal to the product of compression force and eccentricity, particularly where the test configuration involves translational constraints and where the experiments are dynamic. In order to reduce inertial moment artifacts, the mass, and moment of inertia of any loading jig structures that rotate with the specimen should be minimized. Also, the distance between these structures and the load cell should be reduced.

Muscle short-range stiffness behaves like a maxwell element, not a spring: Implications for joint stability.

INTRODUCTION: Muscles play a critical role in supporting joints during activities of daily living, owing, in part, to the phenomenon of short-range stiffness. Briefly, when an active muscle is lengthened, bound cross-bridges are stretched, yielding forces greater than what is predicted from the force length relationship. For this reason, short-range stiffness has been proposed as an attractive mechanism for providing joint stability. However, there has yet to be a forward dynamic simulation employing a cross-bridge model, that demonstrates this stabilizing role. Therefore, the purpose of this investigation was to test whether Huxley-type muscle elements, which exhibit short-range stiffness, can stabilize a joint while at constant activation. METHODS: We analyzed the stability of an inverted pendulum (moment of inertia: 2.7 kg m2) supported by Huxley-type muscle models that reproduce the short-range stiffness phenomenon. We calculated the muscle forces that would provide sufficient short-range stiffness to stabilize the system based in minimizing the potential energy. Simulations consisted of a 50 ms long, 5 Nm square-wave perturbation, with numerical simulations carried out in ArtiSynth. RESULTS: Despite the initial analysis predicting shared activity of antagonist and agonist muscles to maintain stable equilibrium, the inverted pendulum model was not stable, and did not maintain an upright posture even with fully activated muscles. DISCUSSION & CONCLUSION: Our simulations suggested that short-range stiffness cannot be solely responsible for joint stability, even for modest perturbations. We argue that short-range stiffness cannot achieve stability because its dynamics do not behave like a typical spring. Instead, an alternative conceptual model for short-range stiffness is that of a Maxwell element (spring and damper in series), which can be obtained as a first-order approximation to the Huxley model. We postulate that the damping that results from short-range stiffness slows down the mechanical response and allows the central nervous system time to react and stabilize the joint. We speculate that other mechanisms, like reflexes or residual force enhancement/depression, may also play a role in joint stability. Joint stability is due to a combination of factors, and further research is needed to fully understand this complex system.

Compressive follower load influences cervical spine kinematics and kinetics during simulated head-first impact in an in vitro model.

Current understanding of the biomechanics of cervical spine injuries in head-first impact is based on decades of epidemiology, mathematical models, and in vitro experimental studies. Recent mathematical modeling suggests that muscle activation and muscle forces influence injury risk and mechanics in head-first impact. It is also known that muscle forces are central to the overall physiologic stability of the cervical spine. Despite this knowledge, the vast majority of in vitro head-first impact models do not incorporate musculature. We hypothesize that the simulation of the stabilizing mechanisms of musculature during head-first osteoligamentous cervical spine experiments will influence the resulting kinematics and injury mechanisms. Therefore, the objective of this study was to document differences in the kinematics, kinetics, and injuries of ex vivo osteoligamentous human cervical spine and surrogate head complexes that were instrumented with simulated musculature relative to specimens that were not instrumented with musculature. We simulated a head-first impact (3 m/s impact speed) using cervical spines and surrogate head specimens (n = 12). Six spines were instrumented with a follower load to simulate in vivo compressive muscle forces, while six were not. The principal finding was that the axial coupling of the cervical column between the head and the base of the cervical spine (T1) was increased in specimens with follower load. Increased axial coupling was indicated by a significantly reduced time between head impact and peak neck reaction force (p = 0.004) (and time to injury (p = 0.009)) in complexes with follower load relative to complexes without follower load. Kinematic reconstruction of vertebral motions indicated that all specimens experienced hyperextension and the spectrum of injuries in all specimens were consistent with a primary hyperextension injury mechanism. These preliminary results suggest that simulating follower load that may be similar to in vivo muscle forces results in significantly different impact kinetics than in similar biomechanical tests where musculature is not simulated.

Biomechanical Properties of Paraspinal Muscles Influence Spinal Loading-A Musculoskeletal Simulation Study.

Paraspinal muscles are vital to the functioning of the spine. Changes in muscle physiological cross-sectional area significantly affect spinal loading, but the importance of other muscle biomechanical properties remains unclear. This study explored the changes in spinal loading due to variation in five muscle biomechanical properties: passive stiffness, slack sarcomere length (SSL), in situ sarcomere length, specific tension, and pennation angle. An enhanced version of a musculoskeletal simulation model of the thoracolumbar spine with 210 muscle fascicles was used for this study and its predictions were validated for several tasks and multiple postures. Ranges of physiologically realistic values were selected for all five muscle parameters and their influence on L4-L5 intradiscal pressure (IDP) was investigated in standing and 36° flexion. We observed large changes in IDP due to changes in passive stiffness, SSL, in situ sarcomere length, and specific tension, often with interesting interplays between the parameters. For example, for upright standing, a change in stiffness value from one tenth to 10 times the baseline value increased the IDP only by 91% for the baseline model but by 945% when SSL was 0.4 µm shorter. Shorter SSL values and higher stiffnesses led to the largest increases in IDP. More changes were evident in flexion, as sarcomere lengths were longer in that posture and thus the passive curve is more influential. Our results highlight the importance of the muscle force-length curve and the parameters associated with it and motivate further experimental studies on in vivo measurement of those properties.

Anteroposterior shear stiffness of the upper thoracic spine at quasi-static and dynamic loading rates-An in vitro biomechanical study.

To evaluate the biomechanical properties of the upper thoracic spine in anterior-posterior shear loading at various displacement rates. These data broaden our understanding of thoracic spine biomechanics and inform efforts to model the spine and spinal cord injuries. Seven T1-T2 thoracic functional spinal units were loaded non-destructively by a pure shear force up to 200 N, starting from a neutral posture. Tests were run in both posterior and anterior directions, at displacement rates of 1, 10, and 100 mm/s. The three-dimensional motion of the specimen was recorded at 1000 Hz. Individual and averaged load-displacement curves were generated and specimen stiffnesses were calculated. Due to a nonlinear response of the specimens, stiffness was defined separately for both the lower half and the upper half of the specimen range of motion. Specimens were significantly stiffer in the anterior direction than in the posterior direction, across all rates. At low displacements, the anterior stiffness averaged 230 N/mm, 76% higher than the low displacement posterior stiffness of 131 N/mm. At high displacements, anterior stiffness averaged 258 N/mm, 51% stiffer than the high displacement posterior stiffness of 171 N/mm. Shear displacement rate had a small effect on the load response, with the 100 mm/s rate causing a mildly stiffer response at low displacements in the anterior direction. Overall, the load-displacement response exhibited pseudo-quadratic behavior at 1 and 10 mm/s but became more linear at 100 mm/s. The shear stiffness in the upper thoracic spine is greatest in the anterior loading direction, being 51%-76% greater than posterior, most likely due to facet interactions. The effect of the shear displacement rate is low.

The diagnostic precision of computed tomography for traumatic cervical spine injury: An in vitro biomechanical investigation.

BACKGROUND: CT is considered the best method for vertebral fracture detection clinically, but its efficacy in laboratory studies is unknown. Therefore, our objective was to determine the sensitivity, precision, and specificity of high-resolution CT imaging compared to detailed anatomic dissection in an axial compression and lateral bending cervical spine biomechanical injury model. METHODS: 35 three-vertebra human cadaver cervical spine specimens were impacted in dynamic axial compression (0.5 m/s) at one of three lateral eccentricities (low 5% of the spine transverse diameter, middle 50%, high 150%) and two end conditions (19 constrained lateral translation and 16 unconstrained). All specimens were imaged using high resolution CT imaging (246 µm). Two clinicians (spine surgeon and neuroradiologist) diagnosed the vertebral fractures based on 34 discrete anatomical structures using both the CT images and anatomical dissection. FINDINGS: The sensitivity of CT was highest for fractures of the facet joint (59%) and vertebral endplate (57%), and was lowest for pedicle (13%) and lateral mass fractures (23%). The precision of CT was highest for spinous process fractures (83%) and lowest for pedicle (21%), uncinate process and lateral mass (both 23%) fractures. The specificity of CT exceeded 90% for all fractures. The Kappa value between the two reviewers was 0.52, indicating moderate agreement. INTERPRETATION: In this in vitro cervical spine injury model, high resolution CT scanning missed many fractures, notably those of the lateral mass and pedicle. This finding is potentially important clinically, as the integrity of these structures is important to clinical stability and surgical fixation planning.

Automated Segmentation of Spinal Muscles From Upright Open MRI Using a Multiscale Pyramid 2D Convolutional Neural Network.

STUDY DESIGN: Randomized trial. OBJECTIVE: To implement an algorithm enabling the automated segmentation of spinal muscles from open magnetic resonance images in healthy volunteers and patients with adult spinal deformity (ASD). SUMMARY OF BACKGROUND DATA: Understanding spinal muscle anatomy is critical to diagnosing and treating spinal deformity.Muscle boundaries can be extrapolated from medical images using segmentation, which is usually done manually by clinical experts and remains complicated and time-consuming. METHODS: Three groups were examined: two healthy volunteer groups (N = 6 for each group) and one ASD group (N = 8 patients) were imaged at the lumbar and thoracic regions of the spine in an upright open magnetic resonance imaging scanner while maintaining different postures (various seated, standing, and supine). For each group and region, a selection of regions of interest (ROIs) was manually segmented. A multiscale pyramid two-dimensional convolutional neural network was implemented to automatically segment all defined ROIs. A five-fold crossvalidation method was applied and distinct models were trained for each resulting set and group and evaluated using Dice coefficients calculated between the model output and the manually segmented target. RESULTS: Good to excellent results were found across all ROIs for the ASD (Dice coefficient >0.76) and healthy (dice coefficient > 0.86) groups. CONCLUSION: This study represents a fundamental step toward the development of an automated spinal muscle properties extraction pipeline, which will ultimately allow clinicians to have easier access to patient-specific simulations, diagnosis, and treatment.

Temporal Progression of Acute Spinal Cord Injury Mechanisms in a Rat Model: Contusion, Dislocation, and Distraction.

Traumatic spinal cord injuries (SCIs) occur due to different spinal column injury patterns, including burst fracture, dislocation, and flexion-distraction. Pre-clinical studies modeling different SCI mechanisms have shown distinct histological differences between these injuries both acutely (3 h and less) and chronically (8 weeks), but there remains a temporal gap. Different rates of injury progression at specific regions of the spinal cord may provide insight into the pathologies that are initiated by specific SCI mechanisms. Therefore, the objective of this study was to evaluate the temporal progression of injury at specific tracts within the white matter, for time-points of 3 h, 24 h, and 7 days, for three distinct SCI mechanisms. In this study, 96 male Sprague Dawley rats underwent one of three SCI mechanisms: contusion, dislocation, or distraction. Animals were sacrificed at one of three times post-injury: 3 h, 24 h, or 7 days. Histological analysis using eriochrome cyanide and immunostaining for MBP, SMI-312, neurofilament-H (NF-H), and ß-III tubulin were used to characterize white matter sparing and axon and myelinated axon counts. The regions analyzed were the gracile fasciculus, cuneate fasciculus, dorsal corticospinal tract, and ventrolateral white matter. Contusion, dislocation, and distraction SCIs demonstrated distinct damage patterns that progressed differently over time. Myelinated axon counts were significantly reduced after dislocation and contusion injuries in most locations and time-points analyzed (compared with sham). This indicates early myelin damage often within 3 h. Myelinated axon counts after distraction dropped early and did not demonstrate any significant progression over the next 7 days. Important differences in white matter degeneration were identified between injury types, with distraction injuries showing the least variability across time-points These findings and the observation that white matter injury occurs early, and in many cases, without much dynamic change, highlight the importance of injury type in SCI research-both clinically and pre-clinically.

Preliminary investigation of spinal level and postural effects on thoracic muscle morphology with upright open MRI.

OBJECTIVE: Spinal-muscle morphological differences between weight-bearing and supine postures have potential diagnostic, prognostic, and therapeutic applications. While the focus to date has been on cervical and lumbar regions, recent findings have associated spinal deformity with smaller paraspinal musculature in the thoracic region. We aim to quantitatively investigate the morphology of trapezius (TZ), erector spinae (ES) and transversospinalis (TS) muscles in upright postures with open upright MRI and also determine the effect of level and posture on the morphological measures. METHODS: Six healthy volunteers (age 26 ± 6 years) were imaged (0.5 T MROpen, Paramed, Genoa, Italy) in four postures (supine, standing, standing with 30° flexion, and sitting). Two regions of the thorax, middle (T4-T5), and lower (T8-T9), were scanned separately for each posture. 2D muscle parameters such as cross-sectional area (CSA) and position (radius and angle) with respect to the vertebral body centroid were measured for the three muscles. Effect of spinal level and posture on muscle parameters was examined using 2-way repeated measures ANOVA separately for T4-T5 and T8-T9 regions. RESULTS: The TZ CSA was smaller (40%, P = .0027) at T9 than at T8. The ES CSA was larger at T5 than at T4 (12%, P = .0048) and at T9 than at T8 (10%, P = .0018). TS CSA showed opposite trends at the two spinal regions with it being smaller (16%, P = .0047) at T5 than at T4 and larger (11%, P = .0009) at T9 than at T8. At T4-T5, the TZ CSA increased (up to 23%), and the ES and TS CSA decreased (up to 10%) in upright postures compared to supine. CONCLUSION: Geometrical parameters that describe muscle morphology in the thorax change with level and posture. The increase in TZ CSA in upright postures could result from greater activation while upright. The decrease in ES CSA in flexed positions likely represents passive stretching compared to neutral posture.

The effect of vertebral level on biomechanical properties of the lumbar paraspinal muscles in a rat model.

INTRODUCTION: Passive mechanical properties of the paraspinal muscles are important to the biomechanical functioning of the spine. In most computational models, the same biomechanical properties are assumed for each paraspinal muscle group, while cross-sectional area or fatty infiltration in these muscles have been reported to differ between the vertebral levels. Two important properties for musculoskeletal modeling are the slack sarcomere length and the tangent modulus. This study aimed to investigate the effect of vertebral level on these biomechanical properties of paraspinal muscles in a rat model. METHODS: The left paraspinal muscles of 13 Sprague-Dawley rats were exposed under anesthesia. Six muscle biopsies were collected from each rat: three from multifidus (one per each of the L1, L3, and L5 levels) and similarly three from longissimus. Each biopsy was cut into two halves. From one half, two to three single muscle fibers and two to six muscle fiber bundles (14 ± 7 fibers surrounded in their connective tissue) were extracted and mechanically tested in a passive state. From the resulting stress-strain data, tangent modulus was calculated as the slope of the tangent at 30% strain and slack sarcomere length (beyond which passive force starts to develop) was recorded. The other half of each biopsy, which represented the muscle at the fascicle level, was snap frozen, sectioned, stained for Collagen I and its area fraction was measured. To evaluate the effect of spinal level on these biomechanical properties of multifidus and longissimus, one-way repeated measures ANOVA (p < 0.05) was performed for tangent modulus and slack sarcomere length, while for collagen I content linear mixed-models analysis was adopted. RESULTS: In total, 192 fibers and 262 fiber bundles were mechanically tested. For both muscle groups, no significant difference in tangent modulus of the single fibers was detected between the three spinal levels (p = 0.9 for multifidus and p = 0.08 for longissimus). Similarly, the tangent modulus values for the fiber bundles were not significantly different between the three spinal levels (p = 0.13 for multifidus and p = 0.49 for longissimus). In both muscle groups, the slack sarcomere lengths were not different among the spinal levels except for multifidus fibers (p = 0.02). Collagen I area fraction in muscle fascicles averaged 6.8% for multifidus and 5.3% for longissimus and was not different between the spinal levels. DISCUSSION: The results of this study highlighted that the tangent modulus, slack sarcomere length, and collagen I content of the lumbar paraspinal muscles are independent of spinal level. This finding provides the basis for the assumption of similar mechanical properties along a paraspinal muscle group.

Estimation and assessment of sagittal spinal curvature and thoracic muscle morphometry in different postures.

Spine models are typically developed from supine clinical imaging data, and hence clearly do not fully reflect postures that replicate subjects' clinical symptoms. Our objectives were to develop a method to: (i) estimate the subject-specific sagittal curvature of the whole spine in different postures from limited imaging data, (ii) obtain muscle lines-of-action in different postures and analyze the effect of posture on muscle fascicle length, and (iii) correct for cosine between the magnetic resonance imaging (MRI) scan plane and dominant fiber line-of-action for muscle parameters (cross-sectional area (CSA) and position). The thoracic spines of six healthy volunteers were scanned in four postures (supine, standing, flexion, and sitting) in an upright MRI. Geometry of the sagittal spine was approximated with a circular spline. A pipeline was developed to estimate spine geometry in different postures and was validated. The lines-of-action for two muscles, erector spinae (ES) and transversospinalis (TS) were obtained for every posture and hence muscle fascicle lengths were computed. A correction factor based on published literature was then computed and applied to the muscle parameters. The maximum registration error between the estimated spine geometry and MRI data was small (average RMSE∼1.2%). The muscle fascicle length increased (up to 20%) in flexion when compared to erect postures. The correction factor reduced muscle parameters (∼5% for ES and ∼25% for TS) when compared to raw MRI data. The proposed pipeline is a preliminary step in subject-specific modeling. Direction cosines of muscles could be used while improving the inputs of spine models.