Analysis of normal-appearing white matter of multiple sclerosis by tensor-based two-compartment model of water diffusion.

OBJECTIVES: To compare the significance of the two-compartment model, considering diffusional anisotropy with conventional diffusion analyzing methods regarding the detection of occult changes in normal-appearing white matter (NAWM) of multiple sclerosis (MS).
METHODS: Diffusion-weighted images (nine b-values with six directions) were acquired from 12 healthy female volunteers (22-52 years old, median 33 years) and 13 female MS patients (24-48 years old, median 37 years). Diffusion parameters based on the two-compartment model of water diffusion considering diffusional anisotropy was calculated by a proposed method. Other parameters including diffusion tensor imaging and conventional apparent diffusion coefficient (ADC) were also obtained. They were compared statistically between the control and MS groups.
RESULTS: Diffusion of the slow diffusion compartment in the radial direction of neuron fibers was elevated in MS patients (0.121 × 10(-3) mm2/s) in comparison to control (0.100 × 10(-3) mm(2)/s), the difference being significant (P = 0.001). The difference between the groups was not significant in other comparisons, including conventional ADC and fractional anisotropy (FA) of diffusion tensor imaging.
CONCLUSION: The proposed method was applicable to clinically acceptable small data. The parameters obtained by this method improved the detectability of occult changes in NAWM compared to the conventional methods. KEY POINTS: • Water diffusion was compared between the controls and multiple sclerosis patients. • A two-compartment model, considering diffusional anisotropy was selected for water diffusion analysis. • Axial and radial diffusion of fast and slow diffusion components were evaluated. • A new method was developed to obtain the metrics stably. • The metrics indicated high detectability of slight differences between the groups.

Introduction

In multiple sclerosis (MS), diffusion analysis of water molecules using magnetic resonance imaging (MRI) is a well-known quantitative technique for revealing specific microstructural changes [1-3]. It enables the apparent diffusion coefficient (ADC) in the tissues to be calculated by the following mono-exponential equation using two different b-values:where Sb and S0 indicate the signals with and without motion-probing gradients (MPG), b indicates the b-value, and D indicates ADC [4]. Its high sensitivity to occult tissue damage in normal-appearing white matter (NAWM) of MS has been well discussed [5-7].

On the other hand, it is also well known that signal changes do not exactly follow the above equation (Eq. 1) in vivo. Thus, another model that considers two separate diffusion components (fast and slow components) has been well discussed. This two-compartment model is given by this bi-exponential equation:where Df and Ds indicate the diffusion coefficients of the fast and slow diffusion components, respectively, and fs indicates the fraction of the slow diffusion compartment [4]. This equation (Eq. 2) is known to fit the experimental signal changes better than Eq. 1, and thus is often adopted for discussing in vivo water diffusion [4, 8–12]. Of note, this equation is also well known as an intra-voxel incoherent motion (IVIM) model, which usually implements the division of water diffusion from water perfusion in low-b-value data. On the other hand, the two-compartment model used in this study divides the fast and slow diffusion components from the data acquired in low to very high b-values. In this respect, the models are different.

Diffusion tensor imaging (DTI) is a widely used technique that provides anisotropy of the diffusion of water molecules which furnishes the indices of ADC in three dimensions. A previous study found that DTI provides additional information to the usual ADC to quantify the extent and pathological severity of the structural changes occurring within MS plaques and in NAWM [13].

The idea and the technique of DTI may also be applied to Df, Ds, and fs of the two-compartment model (Eq. 2) separately to add information of diffusion anisotropy to them. It will enable us to observe the parameters separately in axial and radial directions, which may provide detailed information of the changes in cell structure and character. It is discussed that water permeability and diffusion restriction of cell membrane is reflected in these parameters [4]. In this regard, it should be better observed in the specific radial direction that is perpendicular to the membrane. It may be useful for detailed diagnosis of MS, because those membrane characteristics may change during the accumulation of diffuse chronic inflammation and demyelination of this disease [1]. However, the technique has not been well discussed. One reason might be that the post-processing operation becomes rather complicated for obtaining adequate results without great noise and errors if the technique is applied straightforward to the small data acquired in clinics. In addition, the complicated post-processing is time-consuming, which also makes it a problem for daily use. In this study, we proposed a new data processing procedure to compensate for these problems. The purpose of this study is to compare the significance of the two-compartment model considering diffusional anisotropy obtained by this proposed method with conventional diffusion analysis regarding the detection of occult changes in NAWM of MS patients.

Materials and methods

This retrospective clinical study was approved by our local ethics committee. It was based on a search of our local clinical database: all the data assessed in this study (diffusion-weighted imaging, DWI) were previously acquired as a part of daily routine for another study (unpublished), and we re-analyzed the data. Written informed consent was obtained from all patients and volunteers at that time.

Subjects

We selected 13 female MS patients (age range 24–48 years, median 37 years) who had relapsing-remitting MS defined according to standard criteria [14]. Those who had either distinguishable brain atrophy or MS plaque at centrum semiovale were not included. The disease duration ranged from 1.8 to 18.4 years, median 6.0 years. In addition, 12 healthy female volunteers (age range 22–52 years, median 33 years) were added to the study as controls.

Imaging procedure and diffusion parameters

Our diffusion-weighted imaging (DWI) data were acquired on a 3 T system (Achieva; Philips Medical Systems, Best, the Netherlands). Twelve different b-values (0, 124, 496, 1116, 1983, 3099, 4463, 6074, 7934, 10,041, 12,397, 15,000 s/mm2), with diffusion encoding in six different directions were applied. Other major imaging parameters are summarized in Table 1. Data from an applied b-value over 8000 s/mm2 were excluded from this study to avoid the effect of image distortion and also because of concern about their low signal-to-noise ratio. We applied two semi-automated masking procedures to the data in order to select solely white matter pixels as much as possible for further assessment. First, to exclude cerebrospinal fluid (CSF) pixels, we excluded the pixels with mean ADC calculated from data of b = 0 and 1116 s/mm2 of over 1.5 × 10−3 mm2/s. Second, to exclude grey matter pixels, we excluded the pixels with fractional anisotropy of diffusion (FA) [13] calculated from the same data of under 0.15. The finally included pixels were regarded as NAWM of MS, because patients who had either distinguishable brain atrophy or MS plaque in FLAIR images were not included. After these operations, smoothing of the remaining image was done by moving average filter (eight connected neighbourhoods).

Table 1

Major parameters of diffusion-weighted imaging sequence

MRI system	3 T Achieva; Philips Medical Systems
Coil	8-channel-array SENSE head coil
b-valuesa	0, 124, 496, 1116, 1983, 3099, 4463, 6074, 7934, 10,041, 12,397, 15,000 s/mm2
Diffusion encoding directions	6
δ /Δ	37.8 /47.3 ms
Acquisition	Single-shot, spin-echo planar imaging
Repetition time /Echo time	4000 /96 ms
Number of signals	one
Section thickness	4 mm
Number of Slices	10 slices
Field of view	256 × 256 mm
Matrixb	64 × 64
Scan time	5 min

aData from an applied b-value over 8000 s/mm2 were excluded from this study

bMatrix was low to avoid low signal-to-noise-ratio in higher b-values

δ: Duration time of the motion probing gradient; Δ, interval time between the motion probing gradients

The diffusion tensor was computed from the data set of b = 0 and 1116 s/mm2. Axial ADC and radial ADC were obtained, which were the first eigenvalue of the diffusion tensor, and the average of the second and third eigenvalues of the tensor. Mean ADC was also obtained by applying the mono-exponential equation to mean DWI (image averaged between the six images of different MPG directions) at the same b-values. Each of the parameters of the two-compartment model in axial or radial direction (axial and radial Ds, Df, and fs) was also estimated following the steps described in the next section (Estimated DWI). In addition, mean Ds, Df, and fs were calculated by applying bi-exponential fitting (Introduction, Eq. 2) to mean DWI (b = 0 to 7934 s/mm2). All the computing operations in this study were performed using our in-house software developed with a commercial analysis package (MATLAB version R2007b, Math Works Inc., Natick, MA, USA).

Estimated DWI (eDWI) and diffusion tensor-based two-compartment model

The estimated DWI (eDWI) method was designed to estimate the virtual b-value-dependent signal change regarding the proper axial or radial diffusion in each pixel. To simulate this virtual imaging, the following steps were carried out with the raw image data set (DWI with nine b-values, 0 to 7994 s/mm2, in six MPG directions).

The diffusion tensor (DTb) and its eigenvalues (λ1b, λ2b, and λ3b) were calculated pixel-by-pixel for every b-value pair that includes b = 0 (e.g., the pair of b = 0 and 124, the pair of b = 0 and 496… and the pair of b = 0 and 7934) (Fig. 1a, b). The first (largest) eigenvalue of each tensor indicates the supposed axial diffusion coefficient of the source pixel at the applied b-value pair, so the first eigenvalue (λ1b) was substituted for D in Eq. 1 (Introduction) to obtain the virtual signal intensity (S ) at that b-value regarding the axial diffusion (Fig. 1c). The whole image series of eDWI regarding axial diffusion was generated by repeating the above steps for all the b-values. Equally, the virtual signal intensity regarding the radial diffusion could be generated by substituting the average of the second and third eigenvalues (λ2b and λ3b) instead of λ1b.

Fig. 1

Procedure for generating axial eDWI. The schema illustrates the procedures to calculate estimated axial diffusion-weighted images (axial eDWI) and shows the virtual b-value-dependent signal change regarding axial diffusion. a The source images are the series of DWI that were acquired with multiple b-values including zero, and multiple motion probing gradient (MPG). Each b-value except zero requires no less than six MPG encoding directions (six in this study). b Diffusion tensor and its eigenvalues (λ1, λ2, and λ3) are calculated for every b-pair of b = 0 and each b-value. c Estimate of axial DWI (eDWI). Virtual signal intensity (Sk) at a certain b-value (b = k) could be obtained as illustrated by a mono-exponential equation with variables of signal intensity at b = 0 (S0), b-value (k), and first eigenvalue (λ1) of the diffusion tensor calculated by the b-value pair of b = 0 and b = k. This process needs to be repeated for every b-value to estimate the whole b-value-dependent signal change

The two-compartment model metrics considering diffusion anisotropy, which were Ds, Df, and fs in proper axial and radial directions, were calculated by fitting bi-exponential curve (Introduction, Eq. 2) to these axial and radial eDWIs.

Comparison between control and MS patients

Two image slices at the level of the centrum semiovale were selected for comparison. The slices were defined as the adjacent two slices right above the bilateral lateral ventricles.

For comparison, first the axial and radial eDWI signal intensities were averaged in each case separately to demonstrate the b-value-dependent signal changes of those virtual images. Second, all the calculated diffusion parameters were averaged in each case for statistical comparison. The control and MS patient groups were compared using the Mann-Whitney rank sum test. P-value under 0.0125 (=0.05/4) was considered significant to avoid type 1 errors in the multiplicity of statistical analysis.

Map images of the proposed diffusion parameters were also made.

Results

The b-value-dependent signal changes of axial eDWI were consistent between the groups, whereas those of radial eDWI differed between the groups in the high b-values (Fig. 2). This result matched the statistical results in which the median of radial Ds was higher in the MS group (0.121 × 10−3 mm2/s) compared to the control group (0.100 × 10−3 mm2/s) with significant difference (P = 0.001), while there was no significant difference either in axial Ds or in Df and fs in both axial and radial directions (Table 2). The mapped images were inhomogeneous, but the results were also consistent with the statistical results (Table 2), where radial Ds was higher in the MS group compared to control, but the other metric maps were not different between the groups (Fig. 3).

Fig. 2

b-value-dependent signal changes in axial and radial eDWI. The semi log graphs indicate the virtual b-value-dependent signal change (estimated diffusion weighted image: eDWI) regarding axial and radial diffusion, respectively. Data points and error bars indicate the averages and standard deviations of the group. The b-value-dependent signal changes are consistent between the groups in axial diffusion, but in radial diffusion there was a difference in higher b-values, where the signal change of the MS group was greater than that in the control group

Table 2

Results of statistical comparison

	Axial				Radial				Mean
	Ds	Df	fs	ADC	Ds	Df	fs	ADC	Ds	Df	fs	ADC	FA
Control	0.276	2.03	0.280	0.825	0.100	1.09	0.309	0.546	0.150	1.29	0.291	0.733	0.268
(min, MAX)	(0.229, 0.378)	(1.75, 2.43)	(0.267, 0.358)	(0.779, 0.885)	(0.0665, 0.119)	(0.966, 1.28)	(0.248, 0.375)	(0.524, 0.564)	(0.144, 0.186)	(1.18, 1.75)	(0.275, 0.311)	(0.697, 0.762)	(0.235, 0.363)
MS	0.267	1.91	0.288	0.818	0.121	1.16	0.321	0.558	0.158	1.35	0.300	0.739	0.274
(min, MAX)	(0.211, 0.403)	(1.49, 6.16)	(0.221, 0.411)	(0.791, 0.858)	(0.0777, 0.154)	(1.04, 1.40)	(0.293, 0.400)	(0.507, 0.606)	(0.124, 0.184)	(1.15, 1.64)	(0.266, 0.308)	(0.708, 0.801)	(0.241, 0.311)
P-valuea	0.849	0.644	0.892	0.978	0.00100	0.0208	0.183	0.109	0.121	0.242	0.265	0.314	0.126

aMann-Whitney rank sum test. Value under 0.0125 (=0.05 /4) was considered significant in this study

Data denote medians and ranges. Ds, Df, and ADC values are written in units of ×10−3 mm2/s

Ds, diffusion coefficient of slow diffusion compartment; Df: diffusion coefficient of fast diffusion compartment; fs: fraction of slow diffusion component; ADC: apparent diffusion coefficient calculated from b-value pair of 0 and 1116 s2/mm; MS, multiple sclerosis

Fig. 3

Sample images of mapped radial Ds, Df, and fs. The sample maps indicate the diffusion parameters of diffusion tensor-based two-compartment model regarding the radial diffusion: Ds (diffusion coefficient of slow diffusion component), Df (diffusion coefficient of fast diffusion component), and fs (fraction of slow diffusion component). The maps (coloured area) are superimposed on a b = 0 DWI image. The maps are restricted to the white matter area, as the cerebral spinal fluid pixels and grey matter pixels were semi-automatically excluded in the masking process (see section of Materials and methods: Diffusion parameters). Colors indicate the values of the scale shown in the colour bar below. Radial Ds was higher in the MS group compared to control, but the other metric maps were not different between the groups. Of note, the parametric images were obtained at 64 × 64 pixels, but the pixels outside the brain parenchyma were trimmed

The differences were not statistically significant between the groups in all conventional DWI and DTI metrics, which were axial, radial, and mean ADC, and FA (Table 2).

Discussion

The proposed method in this study, which was a combination of the “estimated DWI” based on DTI and bi-exponential curve fitting, was designed as a simplified post-processing method to obtain axial and radial Ds, Df, and fs. The previously reported conventional method [10, 12] requires many steps for this process: bi-exponential curve fitting for each MPG direction separately; obtaining the diffusion tensors of Ds, Df, and fs separately from its results; calculating conventional DTI to define axial and radial directions; and finally, generating each metrics by projecting each tensor to the axial and radial directions. On the other hand, the proposed method requires only the diffusion tensor at each b-value paired with b = 0 and two bi-exponential curve fittings (axial and radial). The bi-exponential fitting procedure is usually the most time-consuming and error-generating procedure, so the merit of the proposed method may increase if the number of MPG encoding directions increases.

The estimated b-value-dependent signal change (eDWI) in the proposed method depends on the premise that the direction of the largest eigenvalue of the diffusion tensor is identical to the axial direction in all the b-value pairs used (e.g., pairs of b = 0 and each b-value). Previous studies have reported that the major eigenvectors of the fast and slow diffusions tend to be strongly co-aligned with the known fibre tracts in the regions of high anisotropy [10, 12]. This may support our method, as we selected the ROI by including only pixels with FA of more than 0.15. In the comparison between the control and MS groups, the b-value-dependent signal changes of eDWI were consistent in the axial direction, while a difference was seen in the high b-value areas in the radial direction (Fig. 2). These findings were confirmed by statistical comparison, in which radial Ds was higher in MS with a profoundly significant difference, but the differences in other comparisons were not significant (Table 2); Ds is the diffusion coefficient of the slower diffusion that is reflected most in the signal change at high b-values. Neither of conventional DWI and DTI metrics (axial, radial, and mean ADC by mono-exponential fitting, and FA) showed significant differences between the groups. Therefore, our suggested method may be more sensitive to the occult changes of MS. Furthermore, none of mean Ds, Df, or fs, which were the two-compartment model parameters without diffusion anisotropy information, showed significant differences. Thus, we may consider that the combination of diffusion anisotropy and two-compartment model was important in detecting the changes. The mean value of radial Ds is approximately 50 % less than that of axial Ds, so minimal changes in radial Ds might have been concealed by the larger axial Ds when they were not observed separately.

On the other hand, the map image of Ds, Df, and fs in the control group was still inhomogeneous (Fig. 3). One plausible reason is that there were only six MPG encoding directions in these data, a logical minimum number for post-processing. This might have made the results vulnerable to small noises and outliers even though the model was simplified to fit small data. Of note, the mean Ds, Df, and fs were in good agreement with the values reported earlier [9, 12, 15], indicating that the source image data may have had an acceptable quality (Table 2).

This study has several limitations. First, the sample size was not large. Second, the acquired image resolution was low (64 × 64), which might have caused partial volume effects. Third, the number of MPG encoding directions was the theoretical minimum, six. Vulnerability due to the small data set may have caused some errors throughout the study. Further study with larger sample size and more data may support our results. However, it is important that even the small amount of data did provide useful results with our proposed method, at least with respect to practical application in daily clinics.