Anomalous diffusion parameters are sensible to microstructural variations in brain due to aging
Michele Guerreri1,2, Alessandra Caporale1,2, Marco Palombo1,3, Ivan De Berardinis1, Emiliano Macaluso4, Marco Bozzali4, and Silvia Capuani1

1Department of physics, CNR ISC UOS Roma Sapienza, Rome, Italy, 2Department of anatomical, histological, forensic and of the locomotor system science, Morphogenesis & Tissue Homeostasis, Sapienza University, Rome, Italy, 3MIRCen, CEA/DSV/I2BM, Fontenay-Aux Roses, France, 4Neuroimaging Laboratory, Santa Lucia Foundation, Rome, Italy

Synopsis

We investigated the anomalous diffusion (AD) stretched exponential γ-imaging model to overcome the sensitivity limitations of conventional DTI approach based on the assumption of the Gaussian model with regard of displacements of water molecules in tissues. The benefits of this approach are illustrated with an in-vivo diffusion study of the human brain performed on 18 healthy volunteers in the age range (23-70 years). Mean γ (Mγ) and anisotropic γ (Aγ) maps are obtained and compared with DTI maps. The current study suggests that Mγ and Aγ are more sensitive to micro-structural changes caused by normal aging, compared to DTI metrics.

Introduction

DTI investigations of normal older adults demonstrated the great potential of diffusion studies to detect structural changes occurring in the brain that are associated with changes in cognitive ability. In this study we investigated the so-called anomalous diffusion (AD) stretched exponential γ-imaging model1,2 to overcome the sensitivity limitations of conventional DTI approach based on the assumption of the Gaussian model with regard of displacements of water molecules in tissues. The benefits of this approach are illustrated with an in-vivo diffusion study of the human brain performed on 18 healthy volunteers in the age range (23-70 years). Mean γ (Mγ) and anisotropic γ (Aγ) maps are obtained and compared with DTI maps. The results are discussed in the context of Mγ and Aγ potential applications for brain aging monitoring.

Methodological details

The volunteers (11 men and 7 women, with a mean age +/- SD = 37.8 +/- 15.1 years), after providing informed written consent, underwent MRI examination performed on a 3.0 T Siemens Magnetom Allegra (Siemens Medical Solutions, Erlangen, Germany), equipped with a circularly polarized transmit-receive coil, a maximum gradient strength of 400 mT/m and a maximum slew rate of 400 T/m/s. The same MRI protocol was applied to all the subjects, including Diffusion-Weighted Spin Echo-Echo Planar Imaging (DW SE-EPI) with TR/TE = 6400ms/107ms; Δ/δ = 107ms/35ms; bandwidth = 1860 Hz/px; matrix size = 128 x 128, number of axial slices = 32; in-plane resolution 1.8 x 1.8 mm2; slice thickness 3 mm. The diffusion-encoding gradients were applied along 15 non-collinear directions spanning the entire sphere, and 11 different b-values were used (from 200 to 5000 s/mm2) by varying the gradient strength, plus the b0 image, with an anterior-posterior phase encode direction for all the scans. The acquisition time of the DW-experiment was 42 minutes per subject. The images pre-processing was performed with FSL 5.0 (FMRIB Software Library v5.0, FMRIB, Oxford, UK3). Distortions due to eddy currents were corrected by means of EDDY tool. Mean Diffusivity (MD) and Fractional Anisotropy (FA), together with the 3 diffusion tensor eigenvalues (λ1, λ2, λ3) and eigenvectors (V1, V2, V3) maps were obtained by means of FSL DTIFIT routine, considering b-values comprised in between b0 and b = 1500 s/mm2. Mγ, Aγ, $$$\gamma\bot$$$ and $$$\gamma\parallel$$$ were obtained as described in previous work2, by using a custom-made MATLAB script (MATLAB R2012b), which employs a non-linear least square estimation procedure using trust-region reflective algorithm for minimization, and exploits the advantage of parallel computing. The data were spatially smoothed by means of a Gaussian filter with full-width-half-maximum of 3.2 mm. The non linear fit estimates the parameters minimizing the difference between the signal intensity and the theoretical signal, modeled as follows: $$S(b) \propto\prod_1^3\exp\left[-A_ib^\gamma(V_{ix}G_x+V_{iy}G_y+V_{iz}G_z)^{\gamma_{i}}\right]$$ where b, Gx, Gy, and Gz are arrays of dimension 165 x 1 (where 165 derives from the product between the b-values and the diffusion-encoding directions), Ai are the generalized diffusion constants, γi the three values of the anomalous exponent projected along the 3 main axes of the DTI reference frame2. The 3 γi values were used to compute $$$M_{\gamma} = \sum_{i=1}^3\gamma_{i}$$$, $$$A_{\gamma} = \sqrt{\frac{3\sum_{i=1}^3(\gamma_{i}-M_{\gamma})^{2}}{2(\gamma_{1}^{2}+\gamma_{2}^{2}+\gamma_{3}^{2})}}$$$, $$$\gamma_{\parallel} =\gamma_{3}$$$, and $$$\gamma_{\bot} =(\gamma_{1}+\gamma_{2})/2$$$. AD and conventional DTI metrics were computed in selected regions of interest (ROIs)4,5. Specific Atlases as reference spaces were used with TBSS tool6,7. In a second step we re-projected the Atlases onto each single subject-space, and eroded them in order to avoid voxels at the boundaries. The correlation between AD and DTI metrics and subjects' age was quantified with Pearson’s correlation coefficient in each ROI. P-values < 0.05 were considered statistically significant.

Results and Conclusions

Color-maps showing the R-values of Pearson analysis are displayed in Fig.1-2, respectively. Only zones in which statistical correlation was significant are reported. Hot and cold colors represent positive and negative correlation, respectively. A statistically significant decrease of Mγ and increase of Aγ as a function of age has been observed (see Fig. 3), suggesting that water diffusion becomes more anomalous with aging. This behavior may depend on a global increase of barriers and compartments among white matter fibers and bundles, due to degradation of axons and related macromolecular structures with aging. The current study suggests that Mγ and Aγ are more sensitive to micro-structural changes caused by normal aging, compared to DTI metrics.

Acknowledgements

No acknowledgement found.

References

[1] Palombo M, et al. J Chem Phys 2011;135:034504 [2] De Santis S, et al. Magn Reson Med 2011;65:1043-1052 [3] Jenkinson M, et al. NeuroImage 2012; 62:782-90 [4] Michielse S et al. Neuroimage 2010; 52(4):1190-1201 [5] Bartzokis G et al. Archival Report 2013; 72(12):1026-1034 [6] Smith SM et al. Neuroimage 2006; 31(4):1487-1505 [7] Cherubini A et al. Magn Reson Med 2009; 61(5):1066-1072

Figures

Axial slices of FA maps projected on a common space, with superimposed ROIs where the Pearson's correlation between AD metrics (mean anomalous exponent, Mγ, its anisotropy, Aγ, radial and longitudinal anomalous exponent, $$$\gamma_{\bot}$$$, $$$\gamma_{\parallel}$$$ ) and age was significant. Color-bars indicate the local correlation coefficient, R

Axial slices of FA maps projected on a common space, with superimposed ROIs where the Pearson's correlation between DTI metrics (mean diffusivity, MD, radial diffusivity,$$$D_{\bot}$$$, longitudinal diffusivity, $$$D_{\parallel}$$$) and age was significant. Color-bars indicate the local correlation coefficient, R

Scatter plots of DTI and AD parameters averaged over ROIs in the corpus callosum of each patient, towards the respective age. Parameters ranging from 0 to 1 (FA, $$$M_{\gamma}$$$, $$$A_{\gamma}$$$) share the same units-scale, while MD scale is reported on the right side. The line indicates significant linear correlation



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
3090