Evaluation of diffusional anisotropy and microscopic structure in skeletal muscles using magnetic resonance

The pulsed-gradient spin-echo (PGSE) nuclear magnetic resonance (NMR) method is used for detecting the diffusion of water molecules in biological tissues. Because tissues generally have diffusional anisotropy, their diffusion properties are denoted by a tensor. In this study, we evaluated the diffusional anisotropy and microscopic structure in atrophied skeletal muscles using the PGSE NMR method. The left sciatic nerve was severed in twelve 9-week-old rats. Neurotomy caused neurogenic muscular atrophy at the left gastrocnemius. At 2, 4 and 8 weeks after neurotomy, magnetic resonance signals were selectively acquired from a 2 x 2 x 2 mm(3) voxel, which was located on the left gastrocnemius. The diffusion tensor, the mean diffusivity (MD) and the fractional anisotropy (FA) were calculated from the signals. A theoretical model of the diffusion in muscles was derived from Tanner's equation. The muscle fiber diameter was estimated by fitting the model to the measured signals. The measurements were also performed for normal rats as controls. No significant difference was found in the MD and the estimated intracellular diffusion coefficient between the control group and the denervated group. The denervated group had significantly higher FA compared with the control group (P<.05). The estimated muscle fiber diameter of the denervated group was significantly smaller than the estimated value of the control group (P<.05). These differences were found at 8 weeks after neurotomy. The proposed method is effective for evaluating changes in the microscopic structure of skeletal muscles.     

Simultaneous magnetic resonance imaging of diffusion anisotropy and diffusion gradient

The theory of diffusion gradient-weighted MRI (DGWI) is presented in this paper. The Bloch-Torrey equation was modified to include the effect of intravoxel spatial-location variation of water diffusion (diffusion gradient) on MRI signal, in addition to the effect of intravoxel spatial-direction variation of water diffusion (diffusion anisotropy). An analytical solution for a diffusion-encoding spin-echo pulse sequence was derived. Unlike water diffusion which attenuates the image signal intensity, this newly derived solution relates the spatial gradient of the water diffusion with the phase of the image signal. This novel MRI technique directly measures both the water diffusion and its spatial gradient, and thus offers a noninvasive imaging tool to simultaneously investigate the intravoxel inhomogeneity and anisotropy of tissue structures. In addition, as demonstrated with our preliminary data, this new method may be utilized to delineate the interfaces of tissues with different diffusion. This method is an extension of the successful diffusion tensor MRI (DTI), but requires no additional data acquisition. In addition to the measured diffusion tensor, this new method provides measurements of the spatial derivatives of the three principal diffusivities of the tensor, thereby providing additional information for improving white matter fiber tractography.     

Minimal gradient encoding for robust estimation of diffusion anisotropy

This study has investigated the relationship between the noise sensitivity of measurement by magnetic resonance imaging (MRI) of the diffusion tensor (D) of water and the number N of diffusion-weighting (DW) gradient directions, using computer simulations of strongly anisotropic fibers with variable orientation. The DW directions uniformly sampled the diffusion ellipsoid surface. It is shown that the variation of the signal-to-noise ratio (SNR) of three ideally rotationally invariant scalars of D due to variable fiber orientation provides an objective quantitative measure for the diffusion ellipsoid sampling efficiency, which is independent of the SNR value of the baseline signal obtained without DW; the SNR variation decreased asymptotically with increasing N. The minimum number N(0) of DW directions, which minimized the SNR variation of the three scalars of D was determined, thereby achieving the most efficient ellipsoid sampling. The resulting time efficient diffusion tensor imaging (DTI) protocols provide robust estimation of diffusion anisotropy in the presence of noise and can improve the repeatability/reliability of DTI experiments when there is high variability in the orientation of similar anisotropic structures, as for example, in studies which require repeated measurement of one individual, intersubject comparisons or multicenter studies.     

描述人体内水分子扩散各向异性特征的新方法

通过核磁共振扩散张量成像(DTI)得到的特定值域的扩散各向异性指数(DAI) 可用于揭示水分子扩散椭球的形态学特征, 定量反映被成像物体内部水分子扩散的优势方向和强度, 间接得到被成像物体内部的组织结构信息. DAI的可靠性直接影响对DTI数据的分析和理解. 本文基于扩散张量椭球的几何学信息, 提出利用扩散椭球几何比(EGR)定量描述水分子扩散的各向异性程度. 通过蒙特卡罗模拟实验和对人脑DTI数据进行分析, 并与当前广泛应用的水分子扩散各向异性分数(FA)和近期文献提出的扩散椭球面积比(EAR)进行对比. 实验发现EGR在不同级别噪声影响下的对比度效果和抗噪性都优于FA及EAR. 而且EGR 加入了体积修正, 增强了盘形扩散张量情况下的敏感性, 能够更好地鉴别神经纤维束交叉情况, 对于各向异性扩散程度较高的白质深层和相对均质的表层都有较好的量化区分结果. 关键词： 扩散系数 各向异性扩散 扩散张量成像 扩散椭球几何比     

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.     

Correlation of proton MR spectroscopy and diffusion tensor imaging

Proton magnetic resonance spectroscopy ((1)H-MRS) provides indices of neuronal damage. Diffusion tensor imaging (DTI) relates to water diffusivity and fiber tract orientation. A method to compare (1)H-MRS and DTI findings was developed, tested on phantom and applied on normal brain. Point-resolved spectroscopy (T(R)/T(E)=1500/135) was used for chemical shift imaging of a supraventricular volume of interest of 8 x 8 x 2 cm(3) (64 voxels). In DTI, a segmental spin-echo sequence (T(R)/T(E)=5500/91) was used and slices were stacked to reproduce the slab used in MRS. The spatial distributions of choline and N-acetylaspartate (NAA) correlated to mean fractional anisotropy and apparent diffusion coefficient (ADC) for the inner 6 x 6=36 voxels defined in MRS, most notably NAA and ADC value (r=-.70, P<.00001; correlation across four subjects, 144 data pairs). This is the first association of neuron metabolite contents in volunteers with structure as indicated by DTI.     

Validation of the anisotropy index ellipsoidal area ratio in diffusion tensor imaging

A new diffusion anisotropy index, ellipsoidal area ratio (EAR), was described recently and proved to be less noise-sensitive than fractional anisotropy (FA) by theory and simulation. Here we show that EAR has higher signal-to-noise ratios than FA in average diffusion tensor imaging data from 40 normal subjects. EAR was also more sensitive than FA in detecting white matter abnormalities in a patient with widespread diffuse axonal injury. Monte Carlo simulation showed that EAR's mean values are more biased by noise than FA when anisotropy is small, both for single fiber tracts and when fiber tracts cross. However, the improved signal-to-noise ratio of EAR relative to FA suggests that EAR may be a superior measure of anisotropy both in quantifying both deep white matter with relatively uniform fiber tracts and pericortical white matter structure with relatively low anisotropy and fiber crossings.     

Magnetic field anisotropy based MR tractography

Non-invasive measurements of structural orientation provide unique information regarding the connectivity and functionality of fiber materials. In the present study, we use a capillary model to demonstrate that the direction of fiber structure can be obtained from susceptibility-induced magnetic field anisotropy. The interference pattern between internal and external magnetic field gradients carries the signature of the underlying anisotropic structure and can be measured by MRI-based water diffusion measurements. Through both numerical simulation and experiments, we found that this technique can determine the capillary orientation within 3°. Therefore, susceptibility-induced magnetic field anisotropy may be useful for an alternative tractography method when diffusion anisotropy is small at higher magnetic field strength without the need to rotate the subject inside the scanner.     

The displacement correlation tensor: microstructure, ensemble anisotropy and curving fibers

Experiments with multiple diffusion wave vectors are known to carry more information than what is available from standard diffusion experiments. Here we consider a special case of this class of pulse sequences, the double wave vector diffusion experiment, and use the cumulant expansion of the signal to introduce the displacement correlation tensor. We discuss its physical interpretation and properties, noting in particular that its short time behavior allows determination of the surface to volume ratio of the pore space. We present a general expression for the displacement correlation tensor, and provide explicit expressions for a few model geometries. We then show that the scatter matrix characterizing the orientation distribution of an ensemble of cylinders is simply related to the displacement correlation tensor. This result is generalized to ensembles of pores with arbitrary shapes allowing a precise formulation of the influence of microstructural and ensemble anisotropy on the double wave vector diffusion signal in the Gaussian phase approximation. Finally, as a new application of the double wave vector diffusion signal, we analyze its behavior in a curving fiber, and suggest that the displacement correlation tensor may be used to estimate sub-voxel fiber curvature and deflection angle. The theoretical results are corroborated by computer simulations.     

Effects of random subject rotation on optimised diffusion gradient sampling schemes in diffusion tensor MRI

The choice of the number (N) and orientation of diffusion sampling gradients required to measure accurately the water diffusion tensor remains contentious. Monte Carlo studies have suggested that between 20 and 30 uniformly distributed sampling orientations are required to provide robust estimates of water diffusions parameters. These simulations have not, however, taken into account what effect random subject motion, specifically rotation, might have on optimised gradient schemes, a problem which is especially relevant to clinical diffusion tensor MRI (DT-MRI). Here this question is investigated using Monte Carlo simulations of icosahedral sampling schemes and in vivo data. These polyhedra-based schemes, which have the advantage that large N can be created from optimised subsets of smaller N, appear to be ideal for the study of restless subjects since if scanning needs to be prematurely terminated it should be possible to identify a subset of images that have been acquired with a near optimised sampling scheme. The simulations and in vivo data show that as N increases, the rotational variance of fractional anisotropy (FA) estimates becomes progressively less dependent on the magnitude of subject rotation (), while higher FA values are progressively underestimated as increases. These data indicate that for large subject rotations the B-matrix should be recalculated to provide accurate diffusion anisotropy information.     

A study of rotationally invariant and symmetric indices of diffusion anisotropy.

This study investigated the properties of a class of rotationally invariant and symmetric (relative to the principal diffusivities) indices of the anisotropy of water self-diffusion, namely fractional anisotropy (FA), relative anisotropy (RA), and volume ratio (VR), with particular emphasis to their measurement in brain tissues. A simplified theoretical analysis predicted significant differences in the sensitivities of the anisotropy indices (AI) over the distribution of the principal diffusivities. Computer simulations were used to investigate the effects on AI image quality of three magnetic resonance (MR) diffusion tensor imaging (DTI) acquisition schemes, one being novel: the schemes were simulated on cerebral model fibres varying in shape and spatial orientation. The theoretical predictions and the results of the simulations were corroborated by experimentally determined spatial maps of the AI in a normal feline brain in vivo. We found that FA mapped diffusion anisotropy with the greatest detail and SNR whereas VR provided the strongest contrast between low- and high-anisotropy areas at the expense of increased noise contamination and decreased resolution in anisotropic regions. RA proved intermediate in quality. By sampling the space of the effective diffusion ellipsoid more densely and uniformly and requiring the same total imaging time as the published schemes, the novel DTI scheme achieved greater rotational invariance than the published schemes, with improved noise characteristics, resulting in improved image quality of the AI examined. Our findings suggest that significant improvements in diffusion anisotropy mapping are possible and provide criteria for the selection of the most appropriate AI for a particular application.     

Changes in diffusion tensor imaging (DTI) eigenvalues of skeletal muscle due to hybrid exercise training

Several studies have proposed the cell membrane as the main water diffusion restricting factor in the skeletal muscle cell. We sought to establish whether a particular form of exercise training (which is likely to affect only intracellular components) could affect water diffusion. The purpose of this study is to characterise prospectively the changes in diffusion tensor imaging (DTI) eigenvalues of thigh muscle resulting from hybrid training (HYBT) in patients with non-alcoholic fatty liver disease (NAFLD). Twenty-one NAFLD patients underwent HYBT for 30 minutes per day, twice a week for 6 months. Patients were scanned using DTI of the thigh pre- and post-HYBT. Fractional anisotropy (FA), apparent diffusion coefficient (ADC), the three eigenvalues lambda 1 (λ1), λ2, λ3, and the maximal cross sectional area (CSA) were measured in bilateral thigh muscles: knee flexors (biceps femoris (BF), semitendinosus (ST), semimembranous (SM)) and knee extensors (medial vastus (MV), intermediate vastus (IV), lateral vastus (LV), and rectus femoris (RF)), and compared pre- and post-HYBT by paired t-test. Muscle strength of extensors (P < 0.01), but not flexors, increased significantly post-HYBT. For FA, ADC and eigenvalues, the overall picture was of increase. Some (P < 0.05 in λ2 and P < 0.01 in λ1) eigenvalues of flexors and all (λ1-λ3) eigenvalues of extensors increased significantly (P < 0.01) post-HYBT. HYBT increased all 3 eigenvalues. We suggest this might be caused by enlargement of muscle intracellular space.     

A framework for quality control and parameter optimization in diffusion tensor imaging: theoretical analysis and validation

In this communication, a theoretical framework for quality control and parameter optimization in diffusion tensor imaging (DTI) is presented and validated. The approach is based on the analytical error propagation of the mean diffusivity (D(av)) obtained directly from the diffusion-weighted data acquired using rotationally invariant and uniformly distributed icosahedral encoding schemes. The error propagation of a recently described and validated cylindrical tensor model is further extrapolated to the spherical tensor case (diffusion anisotropy approximately 0) to relate analytically the precision error in fractional tensor anisotropy (FA) with the mean diffusion-to-noise ratio (DNR). The approach provided simple analytical and empirical quality control measures for optimization of diffusion parameter space in an isotropic medium that can be tested using widely available water phantoms.     

磁共振扩散张量成像中扩散敏感梯度磁场方向分布方案的研究进展

磁共振扩散张量成像可以定量无创研究人体内水分子在三维空间中的各向异性扩散规律,进而获取重要的病理及生理信息.为了得到水分子各向异性扩散信息,需要按照一定的方案依次施加不同方向的扩散敏感梯度磁场,测量水分子在这些方向上的扩散系数用以估算扩散张量.扩散张量成像测量结果的准确程度受梯度磁场方向分布方案的影响,本文对扩散敏感梯度磁场方向分布方案进行综述,包括完全随机方案、启发式方案、规则多面体式方案和数值优化方案等,分析这些方案的优势与局限性,并提出需进一步研究的问题.     

An evaluation of the time dependence of the anisotropy of the water diffusion tensor in acute human ischemia

We have performed MRI examinations to determine the water diffusion tensor in the brain of six patients who were admitted to the hospital within 12 h after the onset of cerebral ischemic symptoms. The examinations have been carried out immediately after admission, and thereafter at varying intervals up to 90 days post admission. Maps of the trace of the diffusion tensor, the fractional anisotropy and the lattice index, as well as maps of cerebral blood perfusion parameters, were generated to quantitatively assess the character of the water diffusion tensor in the infarcted area. In patients with significant perfusion deficits and substantial lesion volume changes, four of six cases, our measurements show a monotonic and significant decrease in the diffusion anisotropy within the ischemic lesion as a function of time. We propose that retrospective analysis of this quantity, in combination with brain tissue segmentation and cerebral perfusion maps, may be used in future studies to assess the severity of the ischemic event.     

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.     

Evidence of red cell alignment in the magnetic field of an NMR spectrometer based on the diffusion tensor of water

The alignment of human erythrocytes in aqueous suspensions in the magnetic field B(0) (called the z-direction) of an NMR spectrometer was shown by calculating the diffusion tensor for water in the sample. The diffusion was measured using a pulsed-field-gradient spin-echo NMR method. The extent of diffusion anisotropy for water was exemplified by the values of the apparent diffusion coefficients with erythrocytes of normal shape and volume: for a typical experiment the values for the x-, y-, and z-directions were (6.88 +/- 0.17) x 10(-10), (7.07 +/- 0.17) x 10(-10), and (10.20 +/- 0.17) x 10(-10) m(2) s(-1), respectively. Cells in hypo- and hyperosmotic media were also studied and they too showed the anisotropy of the apparent diffusion coefficients but the extents were different. A new method of data analysis was developed using the Standard Add-On Packages in a Mathematica program. The experimental findings support evidence of erythrocyte alignment that was previously obtained with a high-field-gradient q-space method.     

Diffusion tensor encoding schemes optimized for white matter fibers with selected orientations

A method to produce gradient encoding schemes that minimize the noise of diffusion tensor imaging (DTI) indices for selected fiber orientations has been developed. The accuracy of DTI measurements depends on the gradient encoding scheme used. Most current acquisition schemes contain diffusion directions uniformly distributed in 3D space in order to provide equal noise levels for fibers in any orientation. However, when considering specific fiber bundles such as the corticospinal tract (CST) or parts of fiber bundles, the range of fiber orientations of interest may be limited. We hypothesized that, when studying fiber tracts with a limited range of orientations, measuring diffusion in directions that are uniformly distributed in 3D space may be suboptimal for the noise levels of various DTI indices. Therefore, we first used simulations to determine six diffusion directions that minimize the noise of DTI measurements for selected fiber orientations. The resulting optimized set of directions was then tested on the right CST of a healthy human subject, and its performance was compared with that of conventional acquisition strategies. Both the simulations and the experiments on the human subject demonstrated that the new scheme significantly reduced the standard deviation of DTI indices for tensors with primary eigenvectors within a selected range of orientations.     

Directional sensitivity of anomalous diffusion in human brain assessed by tensorial fractional motion model

Anisotropic diffusion in the nervous system is most commonly modeled by apparent diffusion tensor, which is based on regular diffusion theory. However, the departure of diffusion-induced signal attenuation from a mono-exponential form implies that there is anomalous diffusion. Recently, a novel diffusion NMR theory based on the fractional motion (FM) model, which is an anomalous diffusion model, has been proposed. While the FM model has been applied to both healthy subjects and tumor patients, its anisotropy in the nervous system remains elusive. In this study, this issue was addressed by measuring the FM-related parameters in 12 non-collinear directions. A metric to quantify the directional deviation was derived. Furthermore, the FM-related parameters were modeled as tensors and analyzed in analogy with the conventional diffusion tensor imaging (DTI). Experimental results, which were obtained for 15 healthy subjects at 3T, exhibited pronounced anisotropy of the FM-related parameters, although the effects were smaller than the apparent diffusion coefficient (ADC). The tensorial nature for α, which is the Noah exponent in the FM model, showed behavior similar to the ADC, especially the principal eigenvector for α aligned with the dominant white matter fiber directions. The Hurst exponent H in the FM model, however, showed no correlation with the major fiber directions. The anisotropy of the FM model may provide complementary information to DTI and may have potential for tractography and detecting brain abnormalities.