Arin M. Ellingson

Disc degeneration assessed by quantitative T2* (T2 Star) correlated with functional lumbar mechanics

Study Design. Experimental correlation study design to quantify features of disc health, including signal intensity and distinction between the annulus fibrosus and nucleus pulposus, with T2* magnetic resonance imaging (MRI) and correlate with the functional mechanics in corresponding motion segments. Objective. Establish the relationship between disc health assessed by quantitative T2* MRI and functional lumbar mechanics. Summary of Background Data. Degeneration leads to altered biochemistry in the disc, affecting the mechanical competence. Clinical routine MRI sequences are not adequate in detecting early changes in degeneration and fails to correlate with pain or improve patient stratification. Quantitative T2* relaxation time mapping probes biochemical features and may offer more sensitivity in assessing disc degeneration. Methods. Cadaveric lumbar spines were imaged using quantitative T2* mapping, as well as conventional T2-weighted MRI sequences. Discs were graded by the Pfirrmann scale, and features of disc health, including signal intensity (T2* intensity area) and distinction between the annulus fibrosus and nucleus pulposus (transition zone slope), were quantified by T2*. Each motion segment was subjected to pure moment bending to determine range of motion (ROM), neutral zone (NZ), and bending stiffness. Results. T2* intensity area and transition zone slope were significantly correlated with flexion ROM (P = 0.015; P = 0.002), ratio of NZ/ROM (P = 0.010; P = 0.028), and stiffness (P = 0.044; P = 0.026), as well as lateral bending NZ/ROM (P = 0.005; P = 0.010) and stiffness (P = 0.022; P = 0.029). T2* intensity area was also correlated with lateral bending ROM (P = 0.023). Pfirrmann grade was only correlated with lateral bending NZ/ROM (P = 0.001) and stiffness (P = 0.007). Conclusion. T2* mapping is a sensitive quantitative method capable of detecting changes associated with disc degeneration. Features of disc health quantified with T2* predicted altered functional mechanics of the lumbar spine better than traditional Pfirrmann grading. This new methodology and analysis technique may enhance the assessment of degeneration and enable greater patient stratification for therapeutic strategies. Level of Evidence: N/A

Inter-laboratory variability in in vitro spinal segment flexibility testing

In vitro spine flexibility testing has been performed using a variety of laboratory-specific loading apparatuses and conditions, making test results across laboratories difficult to compare. The application of pure moments has been well established for spine flexibility testing, but to our knowledge there have been no attempts to quantify differences in range of motion (ROM) resulting from laboratory-specific loading apparatuses. Seven fresh-frozen lumbar cadaveric motion segments were tested intact at four independent laboratories. Unconstrained pure moments of 7.5 Nm were applied in each anatomic plane without an axial preload. At laboratories A and B, pure moments were applied using hydraulically actuated spinal loading fixtures with either a passive (A) or controlled (B) XY table. At laboratories C and D, pure moments were applied using a sliding (C) or fixed ring (D) cable-pulley system with a servohydraulic test frame. Three sinusoidal load-unload cycles were applied at laboratories A and B while a single quasistatic cycle was applied in 1.5 Nm increments at laboratories C and D. Non-contact motion measurement systems were used to quantify ROM. In all test directions, the ROM variability among donors was greater than single-donor ROM variability among laboratories. The maximum difference in average ROM between any two laboratories was 1.5° in flexion-extension, 1.3° in lateral bending and 1.1° in axial torsion. This was the first study to quantify ROM in a single group of spinal motion segments at four independent laboratories with varying pure moment systems. These data support our hypothesis that given a well-described test method, independent laboratories can produce similar biomechanical outcomes.

Biomechanical analysis of pedicle screw thread differential design in an osteoporotic cadaver model.

BACKGROUND Pedicle screw fixation, the standard surgical care for posterior stabilization in the thoraco-lumbar spine has a high rate of failure in osteoporotic individuals. Screw design factors and insertion techniques have been shown to influence the biomechanical performance of pedicle screws. Our objective was to investigate the biomechanical characteristics of pedicle screw fixation in osteoporotic bone by comparing standard screws with newly designed differential crest thickness dual lead screws. METHODS An in-vitro spinal-level paired factorial study design was used to examine thoraco-lumbar spine biomechanical outcomes for differential pedicle screw thread designs. Six cadaveric human spines (T8-L5) were tested for six groups (n=20) consisting of 2 different crest thickness and 3 different insertion techniques. Bone mineral density was assessed and peak insertion torque measured while placing one screw of new design and control on the contralateral side. Screw pullout properties were measured from classical American Society for Testing and Materials protocols. FINDINGS The screws designed specifically for osteoporotic bone showed significantly larger insertion torque compared with the standard screw design irrespective of insertion technique. Much of the variability in pullout failure and stiffness was explained by bone mineral density. The osteoporotic screws of different crest thickness were statistically similar to each other in all outcome measures. INTERPRETATION Compared with standard pedicle screws, the dual lead osteoporotic-specific pedicle screws demonstrated significantly larger insertion torques and similar pullout properties. Non-significant increased biomechanical strength was observed for thin crest compared to thick crest dual lead pedicle screws indicating their enhanced purchase in osteoporotic bone.

Quantitative T2* (T2 star) relaxation times predict site specific proteoglycan content and residual mechanics of the intervertebral disc throughout degeneration.

Degeneration alters the biochemical composition of the disc, affecting the mechanical integrity leading to spinal instability. Quantitative T2* MRI probes water mobility within the macromolecular network, a potentially more sensitive assessment of disc health. We determined the relationship between T2* relaxation time and proteoglycan content, collagen content, and compressive mechanics throughout the degenerative spectrum. Eighteen human cadaveric lumbar (L4–L5) discs were imaged using T2* MRI. The T2* relaxation time at five locations (nucleous pulposus or NP, anterior annulus fibrosis or AF, posterior AF, inner AF, and outer AF) was correlated with sulfated‐glycosaminoglycan (s‐GAG) content, hydroxyproline content, and residual stress and strain at each location. T2* relaxation times were significantly correlated with s‐GAG contents in all test locations and were particularly strong in the NP (r = 0.944; p < 0.001) and inner AF (r = 0.782; p < 0.001). T2* relaxation times were also significantly correlated with both residual stresses and excised strains in the NP (r = 0.857; p < 0.001: r = 0.816; p < 0.001), inner AF (r = 0.535; p = 0.022: r = 0.516; p = 0.028), and outer AF (r = 0.668; p = 0.002: r = 0.458; p = 0.041). These strong correlations highlight T2* MRIs ability to predict the biochemical and mechanical health of the disc. T2* MRI assessment of disc health is a clinically viable tool showing promise as a biomarker for distinguishing degenerative changes.

Toward mixed method evaluations of scientific visualizations and design process as an evaluation tool

In this position paper we discuss successes and limitations of current evaluation strategies for scientific visualizations and argue for embracing a mixed methods strategy of evaluation. The most novel contribution of the approach that we advocate is a new emphasis on employing design processes as practiced in related fields (e.g., graphic design, illustration, architecture) as a formalized mode of evaluation for data visualizations. To motivate this position we describe a series of recent evaluations of scientific visualization interfaces and computer graphics strategies conducted within our research group. Complementing these more traditional evaluations our visualization research group also regularly employs sketching, critique, and other design methods that have been formalized over years of practice in design fields. Our experience has convinced us that these activities are invaluable, often providing much more detailed evaluative feedback about our visualization systems than that obtained via more traditional user studies and the like. We believe that if design-based evaluation methodologies (e.g., ideation, sketching, critique) can be taught and embraced within the visualization community then these may become one of the most effective future strategies for both formative and summative evaluations.

Visualizing Motion Data in Virtual Reality: Understanding the Roles of Animation, Interaction, and Static Presentation

We present a study of interactive virtual reality visualizations of scientific motions as found in biomechanics experiments. Our approach is threefold. First, we define a taxonomy of motion visualizations organized by the method (animation, interaction, or static presentation) used to depict both the spatial and temporal dimensions of the data. Second, we design and implement a set of eight example visualizations suggested by the taxonomy and evaluate their utility in a quantitative user study. Third, together with biomechanics collaborators, we conduct a qualitative evaluation of the eight example visualizations applied to a current study of human spinal kinematics. Results suggest that visualizations in this style that use interactive control for the time dimension of the data are preferable to others. Within this category, quantitative results support the utility of both animated and interactive depictions for space; however, qualitative feedback suggest that animated depictions for space should be avoided in biomechanics applications.

Trend-Centric Motion Visualization: Designing and Applying a New Strategy for Analyzing Scientific Motion Collections

In biomechanics studies, researchers collect, via experiments or simulations, datasets with hundreds or thousands of trials, each describing the same type of motion (e.g., a neck flexion-extension exercise) but under different conditions (e.g., different patients, different disease states, pre- and post-treatment). Analyzing similarities and differences across all of the trials in these collections is a major challenge. Visualizing a single trial at a time does not work, and the typical alternative of juxtaposing multiple trials in a single visual display leads to complex, difficult-to-interpret visualizations. We address this problem via a new strategy that organizes the analysis around motion trends rather than trials. This new strategy matches the cognitive approach that scientists would like to take when analyzing motion collections. We introduce several technical innovations making trend-centric motion visualization possible. First, an algorithm detects a motion collections trends via time-dependent clustering. Second, a 2D graphical technique visualizes how trials leave and join trends. Third, a 3D graphical technique, using a median 3D motion plus a visual variance indicator, visualizes the biomechanics of the set of trials within each trend. These innovations are combined to create an interactive exploratory visualization tool, which we designed through an iterative process in collaboration with both domain scientists and a traditionally-trained graphic designer. We report on insights generated during this design process and demonstrate the tools effectiveness via a validation study with synthetic data and feedback from expert musculoskeletal biomechanics researchers who used the tool to analyze the effects of disc degeneration on human spinal kinematics.

Instantaneous helical axis methodology to identify aberrant neck motion

BACKGROUND Neck pain afflicts 30-50% of the U.S. population annually; however we currently have poor diagnostic differentiation techniques to inform individualized treatment. Planar neck kinematics has been shown to be correlated with neck pain, but neck motion is much more complex than pure planar activities. Our objective was to define a methodology for determining aberrant neck kinematics and assess it. METHODS We examined a complex neck kinematic activity of neck circumduction and computed the pathway of motion using the instantaneous helical axis approach in 81 patients with non-specific neck pain and in 20 non-matched symptom free subjects. Neck circumduction, or rolling of the head, represents a complex neck kinematic activity, investigating the innate coupled motion of the cervical spine at the end ranges of motion in all directions. Instance of discontinuities in the helical axis patterns, or folds, were identified and labeled as occurrences of aberrant motion. FINDINGS The instances of aberrant motion, or folds, which are nearly non-existent in the healthy sample group, are present in both the pre- and post-treatment neck pain patients. Following a treatment intervention of the symptomatic patients, pain and neck disability index decreased significantly (P<0.001) concomitant with a decrease in the number of folds (P=0.021). INTERPRETATION The present study highlights a new technique using an instantaneous helical axis approach to detect subtle abnormalities in the pathway of motion of the head about the trunk, during a neck circumduction exercise.

Comparative role of disc degeneration and ligament failure on functional mechanics of the lumbar spine.

Abstract Understanding spinal kinematics is essential for distinguishing between pathological conditions of spine disorders, which ultimately lead to low back pain. It is of high importance to understand how changes in mechanical properties affect the response of the lumbar spine, specifically in an effort to differentiate those associated with disc degeneration from ligamentous changes, allowing for more precise treatment strategies. To do this, the goals of this study were twofold: (1) develop and validate a finite element (FE) model of the lumbar spine and (2) systematically alter the properties of the intervertebral disc and ligaments to define respective roles in functional mechanics. A three-dimensional non-linear FE model of the lumbar spine (L3-sacrum) was developed and validated for pure moment bending. Disc degeneration and sequential ligament failure were modelled. Intersegmental range of motion (ROM) and bending stiffness were measured. The prediction of the FE model to moment loading in all three planes of bending showed very good agreement, where global and intersegmental ROM and bending stiffness of the model fell within one standard deviation of the in vitro results. Degeneration decreased ROM for all directions. Stiffness increased for all directions except axial rotation, where it initially increased then decreased for moderate and severe degeneration, respectively. Incremental ligament failure produced increased ROM and decreased stiffness. This effect was much more pronounced for all directions except lateral bending, which is minimally impacted by ligaments. These results indicate that lateral bending may be more apt to detect the subtle changes associated with degeneration, without being masked by associated changes of surrounding stabilizing structures.

Anatomical 2D/3D shape-matching in virtual reality: A user interface for quantifying joint kinematics with radiographic imaging

We introduce a virtual reality 3D user interface (3DUI) for anatomical 2D/3D shape-matching, a challenging task that is part of medical imaging processes required by biomechanics researchers. Manual shape-matching can be thought of as a nuanced version of classic 6 degree-of-freedom docking tasks studied in the 3DUI research community. Our solution combines dynamic gain for precise translation and rotation from 6 degree-of-freedom tracker input, constraints based on both 2D and 3D data, and immersive visualization and visual feedback.