Diffusion anisotropy in excised normal rat spinal cord measured by NMR microscopy

A conventional spin-echo NMR imaging pulse sequence was used to obtain high-resolution images of excised normal rat spinal cord at 7 and 14 T. It was observed that the large pulsed-field gradients necessary for high-resolution imaging caused a diffusion weighting that dominated the image contrast and that could be used to infer microscopic structural organization beyond that defined by the resolution of the image matrix (i.e., fiber orientation could be assigned based on diffusion anisotropy). Anisotropic diffusion coefficients were therefore measured using apparent diffusion tensor (ADT) imaging to assess more accurately fiber orientations in the spinal cord; structural anisotropy information is portrayed in the six unique images of the complete ADT. To reduce the dimensionality of the data, a trace image was generated using a separate color scale for each of the three diagonal element images of the ADT. This new image retains much of the invariance of the trace to the relative orientations of laboratory and sample axes (inherent to a greyscale trace image) but provides, by the use of color, contrast reflecting diffusion anisotropy. The colored trace image yields a pseudo-three-dimensional view of the rat spinal cord, from which it is possible to deduce fiber orientations.     

In vivo diffusion tensor imaging of thoracic and cervical rat spinal cord at 7 T

In vivo diffusion tensor imaging (DTI) of rat cervical and thoracic spinal cord was performed using a three-element phased array coil at 7 T. The magnetic field was shimmed over the spinal cord in real time using an in-house developed automatic algorithm. Echo planar imaging (EPI)-based diffusion-weighted images (DWIs) were acquired with 21 gradient encoding directions. The DWIs were tensor encoded, and diffusion tensor metrics, fractional anisotropy (FA), mean diffusivity (MD), longitudinal diffusivity (λ₀) and transverse diffusivity (λ) were determined for both white matter (WM) and gray matter (GM). The results on six normal rats indicated no significant differences in the diffusion tensor metrics between thoracic and cervical regions. However, the DTI-derived metrics in cervical spinal cord from our study are somewhat different from the published results in rats. The possible reasons for these differences are suggested.     

An investigation of motion correction algorithms for pediatric spinal cord DTI in healthy subjects and patients with spinal cord injury

Patient and physiological motion can cause artifacts in DTI of the spinal cord which can impact image quality and diffusion indices. The purpose of this investigation was to determine a reliable motion correction method for pediatric spinal cord DTI and show effects of motion correction on DTI parameters in healthy subjects and patients with spinal cord injury. Ten healthy subjects and ten subjects with spinal cord injury were scanned using a 3 T scanner. Images were acquired with an inner field-of-view DTI sequence covering cervical spine levels C1 to C7. Images were corrected for motion using two types of transformation (rigid and affine) and three cost functions. Corrected images and transformations were examined qualitatively and quantitatively using in-house developed code. Fractional anisotropy (FA) and mean diffusivity (MD) indices were calculated and tested for statistical significance pre- and post- motion correction. Images corrected using rigid methods showed improvements in image quality, while affine methods frequently showed residual distortions in corrected images. Blinded evaluation of pre and post correction images showed significant improvement in cord homogeneity and edge conspicuity in corrected images (p < 0.0001). The average FA changes were statistically significant (p < 0.0001) in the spinal cord injury group, while healthy subjects showed less FA change and were not significant. In both healthy subjects and subjects with spinal cord injury, quantitative and qualitative analysis showed the rigid scaled-least-squares registration technique to be the most reliable and effective in improving image quality.     

Dependence of the fractional anisotropy in cervical spine from the number of diffusion gradients,repeated acquisition and voxel size

The aim of this study is to investigate the consequences of using different gradient schemes, number of repeated measurements and voxel size on the fractional anisotropy (FA) value in a diffusion tensor imaging (DTI) sequence on the cervical tract of the spinal cord. Twenty healthy volunteers underwent a total of 86 DTI axial acquisitions performed by using different voxel size and number of diffusion gradient directions (NDGDs). Three different diffusion gradient schemes were applied, named 6, 15 and 32 according to the NDGD. Furthermore, some acquisitions were repeated to investigate the effects of image averaging on FA value.     

Voxelwise analysis of diffusion MRI of cervical spinal cord using tract-based spatial statistics

Robust voxelwise analysis using tract-based spatial statistics (TBSS) together with permutation statistical method is standardly used in analyzing diffusion tensor imaging (DTI) of brain. A similar analytical method could be useful when studying DTI of cervical spinal cord.Based on anatomical data of sixty-four healthy volunteers, white (WM) and gray matter (GM) masks were created and subsequently registered into DTI space. Using TBSS, two skeleton types were created (single line and dilated for WM as well as GM). From anatomical data, percentage rates of overlap were calculated for all skeletons in relation to WM and GM masks.Voxelwise analysis of fractional anisotropy values depending on age and sex was conducted. Correlation of fraction anisotropy values with age of subjects was also evaluated. The two WM skeleton types showed a high overlap rate with WM masks (~94%); GM skeletons showed lower rates (56% and 42%, respectively, for single line and dilated). WM and GM areas where fraction anisotropy values differ between sexes were identified (p < .05). Furthermore, using voxelwise analysis such WM voxels were identified where fraction anisotropy values differ depending on age (p < .05) and in these voxels linear dependence of fraction anisotropy and age (r = −0.57, p < .001) was confirmed by regression analysis. This dependence was not proven when using WM anatomical masks (r = −0.21, p = .10).The analytical approach presented shown to be useful for group analysis of DTI data for cervical spinal cord.     

Structural changes in glutamate cell swelling followed by multiparametric q-space diffusion MR of excised rat spinal cord

Diffusion in the extracellular and intracellular spaces (ECS and ICS, respectively) was evaluated in excised spinal cords, before and after cell swelling induced by glutamate, by high b-value q-space diffusion MR of specific markers and water. The signal decays of deuterated tetramethylammonium (TMA-d(12)) chloride, an exogenous marker of the ECS, and N-acetyl aspartate (NAA), an endogenous marker of the ICS, were found to be non-mono-exponential at all diffusion times. The signal decays of these markers were found to depend on the diffusion time and the cell swelling induced by the glutamate. It was found, for example, that the mean displacements of the apparent fast and slow diffusion components of TMA-d(12) are 7.21 +/- 0.11 and 1.16 +/- 0.05 microm, respectively at a diffusion time of 496 ms. After exposure of the spinal cords to 10 mM of glutamate, these values decreased to 6.62 +/- 0.13 and 1.01 +/- 0.05 microm, respectively. The mean displacement of NAA, however, showed a less pronounced opposite trend and increased after cell swelling induced by exposure to glutamate. q-Space diffusion MR of water was found to be sensitive to exposure to glutamate, and q-space diffusion MRI showed that a more pronounced decrease in the apparent diffusion coefficient and the mean displacement of water is observed in the gray matter (GM) of the spinal cord. All these changes demonstrate that diffusion MR is indeed sensitive to structural changes caused by cell swelling induced by glutamate. Multiparametric high b-value q-space diffusion MR is useful for obtaining microstructural information in neuronal tissues.     

Phase-aligned multiple spin-echo averaging: a simple way to improve signal-to-noise ratio of in vivo mouse spinal cord diffusion tensor image

Purpose

To improve signal-noise-ratio of in vivo mouse spinal cord diffusion tensor imaging using-phase aligned multiple spin-echo technique.

Material and methods

In vivo mouse spinal cord diffusion tensor imaging maps generated by multiple spin-echo and conventional spin-echo diffusion weighting were examined to demonstrate the efficacy of multiple spin-echo diffusion sequence to improve image quality and throughput. Effects of signal averaging using complex, magnitude and phased images from multiple spin-echo diffusion weighting were also assessed. Bayesian probability theory was used to generate phased images by moving the coherent signals to the real channel to eliminate the effect of phase variation between echoes while preserving the Gaussian noise distribution. Signal averaging of phased multiple spin-echo images potentially solves both the phase incoherence problem and the bias of the elevated Rician noise distribution in magnitude image. The proposed signal averaging with Bayesian phase-aligned multiple spin-echo images approach was compared to the conventional spin-echo data acquired with doubling the scan time. The diffusion tensor imaging parameters were compared in the mouse contusion spinal cord injury. Significance level (p-value) and effect size (Cohen’s d) were reported between the control and contused spinal cord to inspect the sensitivity of each approach in detecting white matter pathology.

Results

Compared to the spin-echo image, the signal-noise-ratio increased to 1.84-fold using the phased image averaging and to 1.30-fold using magnitude image averaging in the spinal cord white matter. Multiple spin-echo phased image averaging showed improved image quality of the mouse spinal cord among the tested methods. Diffusion tensor imaging metrics obtained from multiple spin-echo phased images using three echoes and two averages closely agreed with those derived by spin-echo magnitude data with four averages (two times more in acquisition time). The phased image averaging correctly reflected pathological features in contusion spinal cord injury.

Conclusion

Our in vivo imaging results indicate that averaging the phased multiple spin-echo images yields an 84% signal-noise-ratio increase over the spin-echo images and a 41% gain over the magnitude averaged multiple spin-echo images with equal acquisition time. Current results from the animal model of spinal cord injury suggest that the phased multiple spin-echo images could be used to improve signal-noise-ratio.     

Cervical spondylosis: Evaluation of microstructural changes in spinal cord white matter and gray matter by diffusional kurtosis imaging

Introduction

We investigated microstructural changes in the spinal cord, separately for white matter and gray matter, in patients with cervical spondylosis by using diffusional kurtosis imaging (DKI).

Methods

We studied 13 consecutive patients with cervical myelopathy (15 affected sides and 11 unaffected sides). After conventional magnetic resonance (MR) imaging, DKI data were acquired by using a 3 T MR imaging scanner. Values for fractional anisotropy (FA), apparent diffusion coefficient (ADC), and mean diffusional kurtosis (MK) were calculated and compared between unaffected and affected spinal cords, separately for white matter and gray matter.

Results

Tract-specific analysis of white matter in the lateral funiculus showed no statistical differences between the affected and unaffected sides. In gray matter, only MK was significantly lower in the affected spinal cords than in unaffected spinal cords (0.60 ± 0.18 vs. 0.73 ± 0.13, P = 0.0005, Wilcoxon’s signed rank test).

Conclusions

MK values in the spinal cord may reflect microstructural changes and gray matter damage and can potentially provide more information beyond that obtained with conventional diffusion metrics.     

In vivo diffusion characteristics of rat spinal cord.

Complete apparent diffusion tensor (ADTs) of spinal cord was measured in vivo in nine rats at 2.0 T. Two rotationally invariant parameters, the trace, which is a measure of the mean diffusivity, and the lattice index (LI), which reflects the degree of orientation coherence of tissue, have been estimated from the ADT. The mean white matter (WM) trace value (3.05 +/- 0.26 mm2/sec) was found to be substantially higher than the gray matter (GM) trace (2.36 +/- 0.39 mm2/sec), in contrast with the published results on fixed, excised cord. Statistically significant anisotropic diffusion was observed in both WM and GM, with greater anisotropy in the WM (LI = 0.67 +/- 0.06) than in the GM (LI = 0.51 +/- 0.05).     

Relationships between spinal cord blood flow measured with flow-sensitive alternating inversion recovery (FAIR) and neurobehavioral outcomes in rat spinal cord injury

In the traumatically injured spinal cord, decreased perfusion is believed to contribute to secondary tissue damage beyond the primary mechanical impact, and restoration of perfusion is believed to be a promising therapeutic target. However, methods to monitor spinal cord perfusion non-invasively are limited. Perfusion magnetic resonance imaging (MRI) techniques established for the brain have not been routinely adopted to the spinal cord. The purpose of this study was to examine the relationship between spinal cord blood flow (SCBF) and injury severity in a rat thoracic spinal cord contusion injury (SCI) model using flow-sensitive alternating inversion recovery (FAIR) with two variants of the label position. SCBF as a marker of severity was compared to T₁ mapping and to spinal cord-optimized diffusion weighted imaging (DWI) with filtered parallel apparent diffusion coefficient. Thirty-eight rats underwent a T10 contusion injury with varying severities (8 sham; 10 mild; 10 moderate; 10 severe) with MRI performed at 1 day post injury at the lesion site and follow-up neurological assessments using the Basso, Beattie, Bresnahan (BBB) locomotor scoring up to 28 days post injury. Using whole-cord regions of interest at the lesion epicenter, SCBF was decreased with injury severity and had a significant correlation with BBB scores at 28 days post injury. Importantly, estimates of arterial transit times (ATT) in the injured spinal cord were not altered after injury, which suggests that FAIR protocols optimized to measure SCBF provide more value in the context of acute traumatic injury to the cord. T₁-relaxation time constants were strongly related to injury severity and had a larger extent of changes than either SCBF or DWI measures. These findings suggest that perfusion decreases in the spinal cord can be monitored non-invasively after injury, and multi-parametric MRI assessments of perfusion, diffusion, and relaxation capture unique features of the pathophysiology of preclinical injury.     

High b-value q-space diffusion MRI in myelin-deficient rat spinal cords

In this study, we explore the effect of the lack of myelin on the diffusion characteristics and diffusion anisotropy obtained from high b-value q-space diffusion-weighted MRI (q-space DWI) in excised rat spinal cords. Twenty-one-day-old myelin-deficient (md) mutant (N=6) and control rats (N=6) were used in this study. The MRI protocol included multi-slice T(1), T(2), proton density (PD) MR images and high b-value q-space diffusion MRI measured perpendicular and parallel to the fibers of the spine. q-Space displacement and probability maps, in both directions, as well as displacement anisotropy maps, were computed from the diffusion data. At the end of the MRI protocol, representative spinal cords from both groups were subjected to electron microscopy (EM). The md spinal cords show different gray/white matter contrast in the T(1), T(2) and PD MR images as compared with controls. In addition, the mean displacement extracted from the high b-value q-space diffusion data was found to be dramatically higher in the white matter (WM) of the md spinal cords than the controls when diffusion was measured perpendicular and parallel to the fibers of the spine. However, interestingly, at the diffusion time used in the present study, the difference in the WM displacement anisotropies of the two groups was not found to be statistically significant. Myelin was found to have a pronounced effect on the diffusion characteristics of water in WM but less so on the diffusion anisotropy observed at the diffusion time used in the present study.     

q-Space high b value diffusion MRI of hemi-crush in rat spinal cord: evidence for spontaneous regeneration

The development of the damage following hemi-crush trauma in rat spinal cord was studied ex vivo using high b value (bmax = 1 x 10(7) s cm(-2)) q-space diffusion weighted MRI (DWI) at five days, ten days and six weeks post-trauma. Rat spinal cord trauma, produced by hemi-crush of 15s and 60s duration, was studied. The water signal decay in these diffusion experiments was found to be non mono-exponential and was analyzed using the q-space approach. The q-space MRI parameters were compared with T1 and T2 MR images, behavioral tests and histopathological osmium staining. A very good anatomical correlation was found between the q-space MRI parameters and the osmium staining. Interestingly, we found that in the 15s hemi-crush model significant recovery was observed in both the q-space MR images and the osmium staining six weeks post-trauma. However, in the 60s hemi-crush trauma model very little recovery was observed. These results paralleled those obtained from behavioral tests demonstrating that partial spontaneous recovery seems to occur in the 15s hemi-crush spinal cord model, which should be taken in consideration when using it to evaluate new therapies.     

Cervical spinal cord MTR histogram analysis in multiple sclerosis using a 3D acquisition and a B-spline active surface segmentation technique

The application of a three-dimensional magnetization transfer (MT) sequence and B-spline active surface segmentation method to produce MT histograms of the cervical spinal cord in a pilot study of controls and multiple sclerosis (MS) patients is presented. Subjects' cervical spinal cords were imaged with (a) a volume-acquired inversion-prepared fast spoiled gradient echo sequence and (b) a volume-acquired noninversion-prepared fast spoiled gradient echo MT sequence. The images were segmented using the B spline active surface technique and MT histograms were produced from the MT images. The method was sensitive enough to detect differences between seven MS patients and 10 controls in mean MT ratio (42.4 pu versus 44.0 pu, p = 0.03) and peak location (45.2 versus 46.8, p = 0.03). The spinal cord volumes obtained from the two sequences were associated with each other (parameter estimate 0.972, 95% confidence intervals 0.742, 1.202, p < 0.001).     

In vivo echo-planar imaging of rat spinal cord

An integrated approach to echo-planar imaging of rat spinal cord in vivo with a small field of view (FOV) is presented. This protocol is based on a multishot interleaved echo-planar imaging (EPI) sequence and includes: 1) use of an inductively coupled implantable coil for improved signal-to-noise ratio (SNR); 2) three-dimensional (3D) automatic shimming of the magnetic field over the spinal cord; and 3) post-acquisition data processing using a multireference scan for minimizing image artifacts. Some of the practical issues in implementing this protocol are discussed. This imaging protocol will be useful in characterizing the spinal cord pathology using techniques that are otherwise time-consuming, such as diffusion tensor imaging.     

In vivo visualization of displacement-distribution-derived parameters in q-space imaging

OBJECTIVE: This study aimed to explore the potential of in vivo q-space imaging in the differentiation between different cerebral water components. MATERIALS AND METHODS: Diffusion-weighted imaging was performed in six directions with 32 equally spaced q values and a maximum b value of 6600 s/mm(2). The shape of the signal-attenuation curve and the displacement propagator were examined and compared with a normal distribution using the kurtosis parameter. Maps displaying kurtosis, fast and slow components of the apparent diffusion coefficients, fractional anisotropy and directional diffusion were calculated. The displacement propagator was further described by the full width at half and at tenth maximum and by the probability density of zero displacement P(0). Three healthy volunteers and three patients with previously diagnosed multiple sclerosis (MS) were examined. RESULTS: Simulations indicated that the kurtosis of a signal-attenuation curve can determine if more than one water component is present and that care must be taken to select an appropriate threshold. It was possible to distinguish MS plaques in both signal and diffusional kurtosis maps, and in one patient, plaques of different degree of demyelinization showed different behavior. DISCUSSION: Our results indicate that in vivo q-space analysis is a potential tool for the assessment of different cerebral water components, and it might extend the diagnostic interpretation of data from diffusion magnetic resonance imaging.     

Diffusion of particles over strongly anisotropic surface with traps

We investigate surface diffusion in a system of particles adsorbed on a two-dimensional strongly anisotropic lattice. There are two kinds of the lattice sites - ordinary sites and deep traps. Particles adsorbed in the ordinary sites can migrate over the surface, but particles adsorbed in traps are immobile. These particles do not move over the surface and they obstacle also the mobile particles migration (surface defects). Using kinetic Monte Carlo simulations we obtained coverage dependencies of the tracer, jump, and chemical diffusion coefficients. The coefficients are rather sensitive to the defect concentration. Even small admixture of the defects decreases drastically the fast diffusion. The effect is rather specific: strong dependence of the pre-exponential factor on the defect concentration and almost independent activation energy. The defect influence on the slow diffusion is weak. It results in strong decreasing of the surface diffusion anisotropy with the defect concentration. Such unusual behavior of the diffusion coefficients was observed in many experimental investigations of the surface diffusion of lithium, cesium, potassium, and strontium over strongly anisotropic W(1 1 2) and Mo(1 1 2) planes. It was shown that this specific behavior arises exclusively due to the surface anisotropy, and does not depend on the lateral interaction between the particles.     

Correlation between spinal cord diffusion tensor imaging and postural response latencies in persons with multiple sclerosis: A pilot study

PurposeLonger latency of postural response in multiple sclerosis (MS) may be linked to imbalance and increased likelihood of falls. It may be caused by the compromised microstructural integrity in the spinal cord, as evidenced by slowed somatosensory conduction in the spinal cord. Thus, the purpose of this study is to investigate the correlation between latency of postural responses and microstructural integrity of the cervical spinal cord, the region particularly related to the disease severity in MS, using diffusion tensor imaging (DTI) metrics.MethodsSeventeen persons with MS with mild-to-moderate disease severity were enrolled in this study. Postural response latencies of each patient were measured using electromyography of the tibialis anterior muscle (TA) and gastrocnemius muscle (GN) in response to surface perturbations. Cervical spinal cord DTI images were obtained from each patient. DTI mean, radial, axial diffusivity, and fractional anisotropy (FA) were measured between segments C4 and C6. Correlations of DTI metrics with postural response latencies, expanded disability status scale (EDSS) scores, and 25-foot walk (T25FW) were assessed using the Spearman's rank correlation coefficient at α = 0.05.ResultsLower FA was significantly correlated with longer latencies measured on right TA in response to forward postural perturbations (r = −0.51, p = .04). DTI metrics showed no significant correlations with EDSS scores (r = −0.06–0.09, p = .73–0.95) or T25FW (r = −0.1–0.14, p = .6–0.94). DTI metrics showed no significant differences between subjects with and without spinal cord lesions (p = .2–0.7).ConclusionsOur results showed a significant correlation between lower FA in the cervical spinal cord and longer latencies measured on right TA in response to forward postural perturbations in persons with MS, suggesting that impaired cervical spinal cord microstructure assessed by DTI may be associated with the delayed postural responses.     

Mechanism of proton anisotropic velocity distribution in the solar wind

Although it has been long that spacecraft observed the anisotropy of velocity protons in the solar wind, there is still not a reasonable explanation. In this paper we try to give an explanation from the diffusion plateau of protoncyclotron resonance predicted by the quasi-linear theory for the resonance between the protons and the parallel propagating waves. We consider the effect of dispersion relation on diffusion plateau and notice that the diffusion plateau we have got by using cold plasma dispersion relation accords with the density contours in the velocity phase space detected at 0.3 AU in fast solar wind. For explaining proton distributions obtained in the fast solar wind from 0.7 AU to 1 AU hot plasma dispersion relation should be considered. We also give a theoretical relation of proton thermal anisotropy A and plasma parameter β.