Surgical instrumentation for the in vivo determination of lumbar spinal segment stiffness and viscoelasticity

The definition of spinal instability is still controversial. For this reason, it is essential to better understand the difference in biomechanical behaviour between healthy and degenerated human spinal segments in vivo. A novel computer-assisted instrument was developed with the objective to characterize the biomechanical parameters of the spinal segment. Investigation of the viscoelastic properties as well as the dynamic spinal stiffness was performed during a minimally invasive procedure (microdiscectomy) on five patients. Measurements were performed intraoperatively and the protocol consisted of a dynamic part, where spinal stiffness was computed, and a static part, where force relaxation of the segment under constant elongation was studied. The repeatability of the measurement procedure was demonstrated with five replicated tests. The spinal segment tissues were found to have viscoelastic properties. Preliminary tests confirmed a decrease in stiffness after decompression surgery. Patients with non-relaxed muscles showed higher stiffness and relaxation rate compared to patients with relaxed muscles, which can be explained by the contraction and relaxation reflex of muscles under fast and then static elongation. The results show the usefulness of the biomechanical characterization of the human lumbar spinal segment to improve the understanding of the contribution of individual anatomical structures to spinal stability.

Volumetry and diffusion tensor imaging of hippocampal subregions in schizophrenia

Previous MRI-volumetric studies in schizophrenic psychoses have demonstrated more or less pronounced volume reductions of the hippocampus in patients. Correspondingly, neuropathological examinations on the brains of schizophrenics showed diverse structural changes of the hippocampus. Employing a high-resolution 3D-MPRAGE sequence, we found volume reductions in most hippocampal subregions of schizophrenic patients, which, however, did not reach significant levels. An analysis of co-registered diffusion tensor imaging (DTI) data revealed significant alterations of the inter-voxel coherences in single hippocampal subdivisions of these patients, supporting the assumption of characteristic microstructural tissue changes relevant for the pathogenesis of schizophrenic psychoses. Our results argue for the usage of additional MRI modalities like DTI in order to detect subtle regional alterations of hippocampal structure in schizophrenics.

Linking brain connectivity across different time scales with electroencephalogram, functional magnetic resonance imaging, and diffusion tensor imaging

Structural and functional connectivity are intrinsic properties of the human brain and represent the amount of cognitive capacities of individual subjects. These connections are modulated due to development, learning, and disease. Momentary adaptations in functional connectivity alter the structural connections, which in turn affect the functional connectivity. Thus, structural and functional connectivity interact on a broad timescale. In this study, we aimed to explore distinct measures of connectivity assessed by functional magnetic resonance imaging and diffusion tensor imaging and their association to the dominant electroencephalogram oscillatory property at rest: the individual alpha frequency (IAF). We found that in 21 healthy young subjects, small intraindividual temporal IAF fluctuations were correlated to increased blood oxygenation level-dependent signal in brain areas associated to working memory functions and to the modulation of attention. These areas colocalized with functionally connected networks supporting the respective functions. Furthermore, subjects with higher IAF show increased fractional anisotropy values in fascicles connecting the above-mentioned areas and networks. Hence, due to a multimodal approach a consistent functionally and structurally connected network related to IAF was observed.

Rotational integral geometry of tensor valuations

We derive a new rotational Crofton formula for Minkowski tensors. In special cases, this formula gives (1) the rotational average of Minkowski tensors defined on linear subspaces and (2) the functional defined on linear subspaces with rotational average equal to a Minkowski tensor. Earlier results obtained for intrinsic volumes appear now as special cases.

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

Fabric-mechanical property relationships of trabecular bone allografts are altered by supercritical CO2 treatment and gamma sterilization

Tissue grafts are implanted in orthopedic surgery every day. In order to minimize infection risk, bone allografts are often delipidated with supercritical CO2 and sterilized prior to implantation. This treatment may, however, impair the mechanical behavior of the bone graft tissue. The goal of this study was to determine clinically relevant mechanical properties of treated/sterilized human trabecular bone grafts, e.g. the apparent modulus, strength, and the ability to absorb energy during compaction. They were compared with results of identical experiments performed previously on untreated/fresh frozen human trabecular bone from the same anatomical site (Charlebois, 2008). We tested the hypothesis that the morphology–mechanical property relationships of treated cancellous allografts are similar to those of fresh untreated bone. The morphology of the allografts was determined by μCT. Subsequently, cylindrical samples were tested in unconfined and confined compression. To account for various morphologies, the experimental data was fitted to phenomenological mechanical models for elasticity, strength, and dissipated energy density based on bone volume fraction (BV/TV) and the fabric tensor determined by MIL. The treatment/sterilization process does not appear to influence bone graft stiffness. However, strength and energy dissipation of the bone grafts were found to be significantly reduced by 36% to 47% and 66% to 81%, respectively, for a broad range of volume fraction (0.14 < BV/TV < 0.39) and degree of anisotropy (1.24 < DA < 2.18). Since the latter properties are strongly dominated by BV/TV, the clinical consequences of this reduction can be compensated by using grafts with lower porosity. The data of this study suggests that an increase of 5–10% in BV/TV is sufficient to compensate for the reduced post-yield mechanical properties of treated/sterilized bone in monotonic compression. In applications where graft stiffness needs to be matched and strength is not a concern, treated allograft with the same BV/TV as an appropriate fresh bone graft may be used.

Proton diffusion tensor spectroscopy of metabolites in human muscle in vivo

PURPOSE To study the apparent diffusivity and its directionality for metabolites of skeletal muscle in humans in vivo by (1) H magnetic resonance spectroscopy. METHODS The diffusion tensors were determined on a 3 Tesla MR system using optimized acquisition and processing methods including an adapted STEAM sequence with orientation-dependent diffusion weighting, pulse-triggering with individually adapted delays, eddy-current correction schemes, median filtering, and simultaneous prior-knowledge fitting of all related spectra. RESULTS The average apparent diffusivities, as well as the fractional anisotropies of taurine (ADCav  = 0.74 × 10(-3) s/mm(2) , FA = 0.46), creatine (ADCav  = 0.41 × 10(-3)  s/mm(2) , FA = 0.33), trimethylammonium compounds (ADCav  = 0.48 × 10(-3)  s/mm(2) , FA = 0.34), carnosine (ADCav  = 0.46 × 10(-3)  s/mm(2) , FA = 0.47), and water (ADCav  = 1.5 × 10(-3)  s/mm(2) , FA = 0.36) were estimated. The diffusivities of most metabolites and water were significantly different from each other. Diffusion was found to be anisotropic and the diffusion tensors showed tensor correlation coefficients close to 1 and were hence found to be essentially coaligned. The magnitudes of apparent metabolite diffusivities were largely ordered according to molecular weight, with taurine as the smallest molecule diffusing fastest, both along and across the fiber direction. CONCLUSION Diffusivities, directional dependence of diffusion and fractional anisotropies of (1) H MRS-visible muscle metabolites were presented. It was shown that metabolites share diffusion directionality with water and have similar fractional anisotropies, hinting at similar diffusion barriers. Magn Reson Med, 2014. © 2014 Wiley Periodicals, Inc.

Human intervertebral disc stiffness correlates better with the Otsu threshold computed from axial T2 map of its posterior annulus fibrosus than with clinical classifications

Degeneration of the intervertebral disc, sometimes associated with low back pain and abnormal spinal motions, represents a major health issue with high costs. A non-invasive degeneration assessment via qualitative or quantitative MRI (magnetic resonance imaging) is possible, yet, no relation between mechanical properties and T2 maps of the intervertebral disc (IVD) has been considered, albeit T2 relaxation time values quantify the degree of degeneration. Therefore, MRI scans and mechanical tests were performed on 14 human lumbar intervertebral segments freed from posterior elements and all soft tissues excluding the IVD. Degeneration was evaluated in each specimen using morphological criteria, qualitative T2 weighted images and quantitative axial T2 map data and stiffness was calculated from the load-deflection curves of in vitro compression, torsion, lateral bending and flexion/extension tests. In addition to mean T2, the OTSU threshold of T2 (TOTSU), a robust and automatic histogram-based method that computes the optimal threshold maximizing the distinction of two classes of values, was calculated for anterior, posterior, left and right regions of each annulus fibrosus (AF). While mean T2 and degeneration schemes were not related to the IVDs' mechanical properties, TOTSU computed in the posterior AF correlated significantly with those classifications as well as with all stiffness values. TOTSU should therefore be included in future degeneration grading schemes.

Synchrotron x ray induced axonal transections in the brain of rats assessed by high-field diffusion tensor imaging tractography

Since approximately two thirds of epileptic patients are non-eligible for surgery, local axonal fiber transections might be of particular interest for them. Micrometer to millimeter wide synchrotron-generated X-ray beamlets produced by spatial fractionation of the main beam could generate such fiber disruptions non-invasively. The aim of this work was to optimize irradiation parameters for the induction of fiber transections in the rat brain white matter by exposure to such beamlets. For this purpose, we irradiated cortex and external capsule of normal rats in the antero-posterior direction with a 4 mm×4 mm array of 25 to 1000 µm wide beamlets and entrance doses of 150 Gy to 500 Gy. Axonal fiber responses were assessed with diffusion tensor imaging and fiber tractography; myelin fibers were examined histopathologically. Our study suggests that high radiation doses (500 Gy) are required to interrupt axons and myelin sheaths. However, a radiation dose of 500 Gy delivered by wide minibeams (1000 µm) induced macroscopic brain damage, depicted by a massive loss of matter in fiber tractography maps. With the same radiation dose, the damage induced by thinner microbeams (50 to 100 µm) was limited to their paths. No macroscopic necrosis was observed in the irradiated target while overt transections of myelin were detected histopathologically. Diffusivity values were found to be significantly reduced. A radiation dose ≤ 500 Gy associated with a beamlet size of < 50 µm did not cause visible transections, neither on diffusion maps nor on sections stained for myelin. We conclude that a peak dose of 500 Gy combined with a microbeam width of 100 µm optimally induced axonal transections in the white matter of the brain.

Minimally Invasive Intraoperative Stiffness Measurement of Lumbar Spinal Motion Segments

OBJECTIVE: To test a new tiny-tipped intraoperative diagnostic tool that was designed to provide the surgeon with reliable stiffness data on the motion segment during microdiscectomy. A decrease in stiffness after nuclectomy and a measurable influence of muscle tension were assumed. If the influence of muscle tension on the motion segment could at least be ruled out, there should be no difference with regard to stiffness between women and men. If these criteria are met, this new intraoperative diagnostic tool could be used in further studies for objective decision-making regarding additional stabilization systems after microdiscectomy. METHODS: After evaluation of the influence of muscle relaxation during in vivo measurements with a spinal spreader between the spinous processes, 21 motion segments were investigated in 21 patients. Using a standardized protocol, including quantified muscle relaxation, spinal stiffness was measured before laminotomy and after nuclectomy. RESULTS: The decrease in stiffness after microdiscectomy was highly significant. There were no statistically significant differences between men and women. The average stiffness value before discectomy was 33.7 N/mm, and it decreased to 25.6 N/mm after discectomy. The average decrease in stiffness was 8.1 N/mm (24%). CONCLUSION: In the moderately degenerated spine, stiffness decreases significantly after microdiscectomy. Control for muscle relaxation is essential when measuring in vivo spinal stiffness. The new spinal spreader was found to provide reliable data. This spreader could be used in further studies for objective decision-making about additional stabilization systems after microdiscectomy.

Finite Element Based Nonlinear Normalization of Human Lumbar Intervertebral Disc Stiffness to Account for Its Morphology

Disc degeneration, usually associated with low back pain and changes of intervertebral stiffness, represents a major health issue. As the intervertebral disc (IVD) morphology influences its stiffness, the link between mechanical properties and degenerative grade is partially lost without an efficient normalization of the stiffness with respect to the morphology. Moreover, although the behavior of soft tissues is highly nonlinear, only linear normalization protocols have been defined so far for the disc stiffness. Thus, the aim of this work is to propose a nonlinear normalization based on finite elements (FE) simulations and evaluate its impact on the stiffness of human anatomical specimens of lumbar IVD. First, a parameter study involving simulations of biomechanical tests (compression, flexion/extension, bilateral torsion and bending) on 20 FE models of IVDs with various dimensions was carried out to evaluate the effect of the disc's geometry on its compliance and establish stiffness/morphology relations necessary to the nonlinear normalization. The computed stiffness was then normalized by height (H), cross-sectional area (CSA), polar moment of inertia (J) or moments of inertia (Ixx, Iyy) to quantify the effect of both linear and nonlinear normalizations. In the second part of the study, T1-weighted MRI images were acquired to determine H, CSA, J, Ixx and Iyy of 14 human lumbar IVDs. Based on the measured morphology and pre-established relation with stiffness, linear and nonlinear normalization routines were then applied to the compliance of the specimens for each quasi-static biomechanical test. The variability of the stiffness prior to and after normalization was assessed via coefficient of variation (CV). The FE study confirmed that larger and thinner IVDs were stiffer while the normalization strongly attenuated the effect of the disc geometry on its stiffness. Yet, notwithstanding the results of the FE study, the experimental stiffness showed consistently higher CV after normalization. Assuming that geometry and material properties affect the mechanical response, they can also compensate for one another. Therefore, the larger CV after normalization can be interpreted as a strong variability of the material properties, previously hidden by the geometry's own influence. In conclusion, a new normalization protocol for the intervertebral disc stiffness in compression, flexion, extension, bilateral torsion and bending was proposed, with the possible use of MRI and FE to acquire the discs' anatomy and determine the nonlinear relations between stiffness and morphology. Such protocol may be useful to relate the disc's mechanical properties to its degree of degeneration.