Alan S. Barnett

Water diffusion changes in Wallerian degeneration and their dependence on white matter architecture.

This study investigates water diffusion changes in Wallerian degeneration. We measured indices derived from the diffusion tensor (DT) and T2-weighted signal intensities in the descending motor pathways of patients with small chronic lacunar infarcts of the posterior limb of the internal capsule on one side. We compared these measurements in the healthy and lesioned sides at different levels in the brainstem caudal to the primary lesion. We found that secondary white matter degeneration is revealed by a large reduction in diffusion anisotropy only in regions where fibers are arranged in isolated bundles of parallel fibers, such as in the cerebral peduncle. In regions where the degenerated pathway crosses other tracts, such as in the rostral pons, paradoxically there is almost no change in diffusion anisotropy, but a significant change in the measured orientation of fibers. The trace of the diffusion tensor is moderately increased in all affected regions. This allows one to differentiate secondary and primary fiber loss where the increase in trace is considerably higher. We show that DT-MRI is more sensitive than T2-weighted MRI in detecting Wallerian degeneration. Significant diffusion abnormalities are observed over the entire trajectory of the affected pathway in each patient. This finding suggests that mapping degenerated pathways noninvasively with DT-MRI is feasible. However, the interpretation of water diffusion data is complex and requires a priori information about anatomy and architecture of the pathway under investigation. In particular, our study shows that in regions where fibers cross, existing DT-MRI-based fiber tractography algorithms may lead to erroneous conclusion about brain connectivity.

Comprehensive Approach for Correction of Motion and Distortion in Diffusion-Weighted MRI

Patient motion and image distortion induced by eddy currents cause artifacts in maps of diffusion parameters computed from diffusion‐weighted (DW) images. A novel and comprehensive approach to correct for spatial misalignment of DW imaging (DWI) volumes acquired with different strengths and orientations of the diffusion sensitizing gradients is presented. This approach uses a mutual information‐based registration technique and a spatial transformation model containing parameters that correct for eddy current‐induced image distortion and rigid body motion in three dimensions. All parameters are optimized simultaneously for an accurate and fast solution to the registration problem. The images can also be registered to a normalized template with a single interpolation step without additional computational cost. Following registration, the signal amplitude of each DWI volume is corrected to account for size variations of the object produced by the distortion correction, and the b‐matrices are properly recalculated to account for any rotation applied during registration. Both qualitative and quantitative results show that this approach produces a significant improvement of diffusion tensor imaging (DTI) data acquired in the human brain. Magn Reson Med 51:103–114, 2004. Published 2003 Wiley‐Liss, Inc.

Visualizing and characterizing white matter fiber structure and architecture in the human pyramidal tract using diffusion tensor MRI.

We used diffusion tensor imaging to assess diffusion anisotropy in the pyramidal tract in ten young, and ten elderly subjects (five males and five females in each group). The purpose of this study was to define normative values for anisotropy at different anatomic levels of the brainstem as well as to assess differences due to age, gender, and laterality. In all subjects, anisotropy was highest in the cerebral peduncle, lowest in the caudal pons, and intermediate in the medulla. In the pons and medulla the regional variability was high, with significant differences in anisotropy even between contiguous slices. Multifactorial ANOVA (performed using the average value of anisotropy within each region of interest) revealed that elderly subjects had significantly lower values than young subjects in the cerebral peduncle, with no differences in the pons and medulla. No significant differences in anisotropy due to gender and side were found. The differences in anisotropy at different levels of the brainstem reflect differences in the local architecture of white matter fibers. Anisotropy is high in the cerebral peduncle because fibers have a highly ordered arrangement, while in the pons and medulla, anisotropy is lower because the local fiber architecture is less coherent due to the presence of other fibers and nuclei. The biologic meaning of the intergroup differences in anisotropy is discussed in light of the structure and architecture of the tissue under investigation. We also consider potential sources of artifacts, such as noise and motion, partial volume contamination, anatomic mismatching, and the use of inappropriate statistical tests. We conclude that the age-related decrease in anisotropy in the cerebral peduncle is not artifactual but rather reflects subtle structural changes of the aging white matter. Our study however shows that caution must be exercised in interpreting diffusion anisotropy data.

High Temporal Resolution Diffusion MRI of Global Cerebral Ischemia and Reperfusion

Although brain ischemia has been extensively studied using diffusion-weighted magnetic resonance imaging, most studies performed so far have not had adequate time resolution to follow the temporal changes in the water apparent diffusion coefficient (ADC) in hyperacute ischemia. Using diffusion echo planar imaging, we obtained ADC maps (calculated from measurements made with 8 b-values) with a time resolution of 43 s in a feline model of global brain ischemia and reperfusion. Different protocols were performed: 10-min hypoperfusion, 10- and 22-min ischemia followed by reperfusion, and cardiac arrest. ADC values were obtained from white matter of the internal capsule and from the thalamus. Cortical gray matter measurements were not deemed reliable due to the close proximity of CSF in the cortical sulci. Following occlusion, the ADC declined in the thalamus to <2 SD of its normal baseline value within 1.5–2.5 min. This decay was exponential with a time constant (τ ± SD) of 6.0 ± 2.6 min; no further decrease in the ADC was observed 10 min following ischemia. Following reperfusion, in animals that showed ADC recovery, the ADC began increasing immediately, returning to its preischemic value in ∼15 min. No significant ADC changes were observed during hypoperfusion. Following cardiac arrest, the decay of ADC was more rapid in the thalamus (τ = 2.6 ± 0.6 min) than in white matter (τ = 6.6 ± 1.8 min). We observed that the ADC at 40 min after cardiac arrest was similar to the ADC at 10 min after ischemia. Given that all animals subjected to 10-min ischemic episodes showed ADC recovery with reperfusion, doubt is cast on whether it is possible to define a threshold value of the ADC below which brain tissue is irreversibly damaged. Finally, despite variability in the time constants of the ADC decay induced by ischemia, the ADC values at 10 min were very similar in all the animals. This suggests that when blood flow is diminished sufficiently to induce an ADC reduction, differences in perfusion affect the rapidity of the decrease but not the final asymptotic value reached.

Generalized reconstruction of phase contrast MRI: Analysis and correction of the effect of gradient field distortions†

To characterize gradient field nonuniformity and its effect on velocity encoding in phase contrast (PC) MRI, a generalized model that describes this phenomenon and enables the accurate reconstruction of velocities is presented. In addition to considerable geometric distortions, inhomogeneous gradient fields can introduce deviations from the nominal gradient strength and orientation, and therefore spatially‐dependent first gradient moments. Resulting errors in the measured phase shifts used for velocity encoding can therefore cause significant deviations in velocity quantification. The true magnitude and direction of the underlying velocities can be recovered from the phase difference images by a generalized PC velocity reconstruction, which requires the acquisition of full three‐directional velocity information. The generalized reconstruction of velocities is applied using a matrix formalism that includes relative gradient field deviations derived from a theoretical model of local gradient field nonuniformity. In addition, an approximate solution for the correction of one‐directional velocity encoding is given. Depending on the spatial location of the velocity measurements, errors in velocity magnitude can be as high as 60%, while errors in the velocity encoding direction can be up to 45°. Results of phantom measurements demonstrate that effects of gradient field nonuniformity on PC‐MRI can be corrected with the proposed method. Magn Reson Med 50:791–801, 2003. Published 2003 Wiley‐Liss, Inc.

Investigation of Anatomical Thalamo-Cortical Connectivity and fMRI Activation in Schizophrenia

The purpose of this study was to examine measures of anatomical connectivity between the thalamus and lateral prefrontal cortex (LPFC) in schizophrenia and to assess their functional implications. We measured thalamocortical connectivity with diffusion tensor imaging (DTI) and probabilistic tractography in 15 patients with schizophrenia and 22 age- and sex-matched controls. The relationship between thalamocortical connectivity and prefrontal cortical blood-oxygenation-level-dependent (BOLD) functional activity as well as behavioral performance during working memory was examined in a subsample of 9 patients and 18 controls. Compared with controls, schizophrenia patients showed reduced total connectivity of the thalamus to only one of six cortical regions, the LPFC. The size of the thalamic region with at least 25% of model fibers reaching the LPFC was also reduced in patients compared with controls. The total thalamocortical connectivity to the LPFC predicted working memory task performance and also correlated with LPFC BOLD activation. Notably, the correlation with BOLD activation was accentuated in patients as compared with controls in the ventral LPFC. These results suggest that thalamocortical connectivity to the LPFC is altered in schizophrenia with functional consequences on working memory processing in LPFC.

Comparison of EPI Distortion Correction Methods in Diffusion Tensor MRI Using a Novel Framework

Diffusion weighted images (DWIs) are commonly acquired with Echo-planar imaging (EPI). B0 inhomogeneities affect EPI by producing spatially nonlinear image distortions. Several strategies have been proposed to correct EPI distortions including B0 field mapping (B0M) and image registration. In this study, an experimental framework is proposed to evaluation the performance of different EPI distortion correction methods in improving DT-derived quantities. A deformable registration based method with mutual information metric and cubic B-spline modeled constrained deformation field (BSP) is proposed as an alternative when B0 mapping data are not available. BSP method is qualitatively and quantitatively compared to B0M method using the framework. Both methods can successful reduce EPI distortions and significantly improve the quality of DT-derived quantities. Overall, B0M was clearly superior in infratentorial regions including brainstem and cerebellum, as well as in the ventral areas of the temporal lobes while BSP was better in all rostral brain regions.

Analysis and generalized correction of the effect of spatial gradient field distortions in diffusion-weighted imaging.

Nonuniformities of magnetic field gradients can cause serious artifacts in diffusion imaging. While it is well known that nonlinearities of the imaging gradients lead to image warping, those imperfections can also cause spatially dependent errors in the direction and magnitude of the diffusion encoding. This study shows that the potential errors in diffusion imaging are considerable. Further, we show that retrospective corrections can be applied to reduce these errors. A general mathematical framework was formulated to characterize the contribution of gradient nonuniformities to diffusion experiments. The gradient field was approximated using spherical harmonic expansion, and this approximation was employed (after geometric distortions were eliminated) to predict and correct the errors in diffusion encoding. Before the corrections were made, the experiments clearly revealed marked deviations of the calculated diffusivity for fields of view (FOVs) generally used in diffusion experiments. These deviations were most significant farther away from the magnets isocenter. For an FOV of 25 cm, the resultant errors in absolute diffusivity ranged from approximately –10% to +20%. Within the same FOV, the diffusion‐encoding direction and the orientation of the calculated eigenvectors can be significantly altered if the perturbations by the gradient nonuniformities are not considered. With the proposed correction scheme, most of the errors introduced by gradient nonuniformities can be removed. Magn Reson Med 50:560–569, 2003.

Dependence on diffusion time of apparent diffusion tensor of ex vivo calf tongue and heart.

The time dependence of the apparent diffusion tensor of ex vivo calf heart and tongue was measured for diffusion times (τd) between 32 and 810 ms. The results showed evidence of restricted diffusion in the muscle tissues of both organs. In regions where the myofibers are parallel, the largest eigenvalue (λ1) of the diffusion tensor remained the same for all diffusion times measured, while the other eigenvalues (λ2, λ3) decreased by 29–36% between τd = 32 ms and τd = 400 ms. In regions where the fibers cross, the λ1 also changed, decreasing by 17% between τd = 32 ms and τd = 400 ms. The restricting compartment size and volume fraction were effectively estimated by fitting the time courses of the eigenvalues to a model consisting of a nonrestricted compartment and a cylindrically restricted compartment. To our knowledge, this study is the first demonstrating diffusion time dependence of measured water diffusion tensor in muscular tissue. With improvement in scanning technology, future studies may permit noninvasive, in vivo detection of changes in muscle myoarchitecture due to disease, treatment, and exercise. Magn Reson Med, 2005. Published 2005 Wiley‐Liss, Inc.

Regional distribution of measurement error in diffusion tensor imaging

The characterization of measurement error is critical in assessing the significance of diffusion tensor imaging (DTI) findings in longitudinal and cohort studies of psychiatric disorders. We studied 20 healthy volunteers, each one scanned twice (average interval between scans of 51 +/- 46.8 days) with a single shot echo planar DTI technique. Intersession variability for fractional anisotropy (FA) and Trace (D) was represented as absolute variation (standard deviation within subjects: SDw), percent coefficient of variation (CV) and intra-class correlation coefficient (ICC). The values from the two sessions were compared for statistical significance with repeated measures analysis of variance or a non-parametric equivalent of a paired t-test. The results showed good reproducibility for both FA and Trace (CVs below 10% and ICCs at or above 0.70 in most regions of interest) and evidence of systematic global changes in Trace between scans. The regional distribution of reproducibility described here has implications for the interpretation of regional findings and for rigorous pre-processing. The regional distribution of reproducibility measures was different for SDw, CV and ICC. Each one of these measures reveals complementary information that needs to be taken into consideration when performing statistical operations on groups of DT images.