Imaging Three Dimensional Temporal Diffusion Spectrum Dispersion Profiles in the Brain
Dan Wu1, Frances J Northington2, and Jiangyang Zhang1,3

1Radiology, Johns Hopkins University School of Medicine, BALTIMORE, MD, United States, 2Pediatrics, Johns Hopkins University School of Medicine, BALTIMORE, MD, United States, 3Radiology, New York University School of Medicine, New Yourk, NY, United States

Synopsis

The dispersion profile of the temporal diffusion spectrum has been linked to key properties of tissue microstructures, however, its directional variance has not been shown. In this study, we extended the conventional one-dimensional dispersion profile to three-dimensional profile, and characterized its directionality with a tensor representation. The temporal diffusion dispersion (TDD) tensor demonstrated unique contrasts that reflected distinct microstructural organization in the mouse brain, and the high anisotropy from TDD tensors correlated with anisotropic structural arrangements, e.g., in the crossing fiber regions. The TDD contrasts are also sensitive to disrupted microstructures in a neonatal mouse model of hypoxic-ischemic injury.

Introduction

Diffusion MRI is a sensitive tool to detect tissue microstructural organization in the brain under normal and pathological conditions. With oscillating gradient diffusion MRI, it is now possible to sample a wide range of the temporal diffusion spectrum to gather a wealth of information on tissue microstructural organization1,2,3. In particular, the dispersion profile of the temporal diffusion spectrum4, which can be approximated with the rate of frequency dependent change for the low to moderate frequency range (e.g., < 300 Hz), has been linked to cell size and surface volume ratio5,6. In this study, we investigated the dispersion profiles along different directions in the mouse brain and explore the potential relationships between the observed directional anisotropy with underlying microstructural organization.

Methods

Neonatal C57BL/6 mice (n=6) at postnatal day 10 were subjected to unilateral HI using the Vannucci model7, and sacrificed at 24hrs after injury and perfusion fixed with 4% PFA. Ex vivo diffusion MRI experiments were performed on an 11.7 T spectrometer with a 10 mm transceiver volume coil. Pulsed gradient spin-echo (PGSE) data were acquired with δ/Δ=2.6/20ms, and oscillating gradient spin-echo (OGSE) data were acquired at frequencies (f) of 50 Hz, 100 Hz, and 200 Hz. All data were acquired with TE/TR = 60/1000 ms, NA=2, 30 diffusion directions, b = 1000 s/mm2, and resolution of 0.125 mm isotropic. Diffusion tensors were estimated from diffusion MRI data at each frequency and the datasets were co-registered with affined registration based on the b0 images. Along each diffusion-encoding direction, the rate of apparent diffusion coefficient (ADC) increase with respect to the square root of frequency was estimated, and the results along 30 directions were fitted to a tensor model, which called the temporal diffusion dispersion (TDD) tensor here. The scalar indices such as apparent diffusion coefficient (ADC), fractional anisotropy (FA), and colormap were calculated from the TDD tensor in the same way as conventional diffusion tensor.

Results

Fig. 1 showed the rates of ADC change with the square root of frequency measured in several gray and white matter structures in the contralateral side of the mouse brain where the tissues were intact (n=6), and the change rates demonstrated directional anisotropy along x/y/z over the 0 (PGSE) to 200 Hz range (Fig. 1D). The trace map of TDD tensor highlighted regions with high density of neurons in the mouse brain, as reported before8. The FA map of TDD tensor highlighted several regions with distinct microstructural organization (Fig. 2). For example, several large white matter structures, e.g., the cerebral peduncle, had high diffusion anisotropy in both conventional PGSE/OGSE as well as in TDD (Table1), with the primary eigen-vector of the TDD tensor mostly parallel to that of PGSE tensor. Interesting, the splenium of the corpus callosum showed a relatively low FA, whereas as the region below the corpus callosum—the dorsal hippocampal commissure (dhc) showed high FA values from the TDD tensor compared to the PGSE FA values (Table 1). Fiber orientation density (FOD) reconstructed using constraint spherical deconvolution9 based on the PGSE data showed that the dhc contained crossing fibers with two fiber group along left-right and anterior-posterior directions (orange arrows in Fig. 2). While the PGSE-colormap showed an AP-oriented eigen-vector (blue) in the dhc, the TDD tensor indicated a LR-oriented eigen-vector (red) in this region. In the hippocampus, a layer (yellow arrows in Fig. 2) with low PGSE FA values showed high FA values from TDD tensor (Table1). This layer is approximately in the location or stratum oriens, in which several groups of axons and dendrites with different directions mixed together. In the ipsilateral side of the hypoxia-ischemia injured mouse brains, the well-organized hippocampal layers were disrupted (Fig. 3). The hippocampal injury was marked by reduction of trace and FA from the TDD tensor.

Discussion and conclusion

It is known that temporal diffusion spectrum provides information on the spatial scale of tissue microstructural organization, such as cell size and surface-to-volume ratio, etc5,6. It is also well known that, in several brain regions, the spatial organizations of microstructures are different along different directions, which may lead to the high anisotropy in the TDD tensor. When these unique microstructural organizations were disrupted by injuries, e.g., hypoxic ischemic injury in this study, both the conventional tensor and TTD tensor showed changes in anisotropy. In summary, the new method reported here revealed an interesting new contrast, which potentially offers additional information about the microstructural organization.

Acknowledgements

No acknowledgement found.

References

1. Stepisnik J. Time-dependent self-diffusion by NMR spin-echo. Physica B. 1993;183:343–350.

2. Parsons EC Jr, Does MD, Gore JC. Temporal diffusion spectroscopy: theory and implementation in restricted systems using oscillating gradients. Magn Reson Med. 2006;55:75–84.

3. Does MD, Parsons EC, Gore JC. Oscillating gradient measurements of water diffusion in normal and globally ischemic rat brain. Magn Reson Med. 2003;49:206–215.

4. Gore J, Xu J, Colvin D, Yankeelov T, et al. Characterization of tissue structure at varying length scales using temporal diffusion spectroscopy. NMR Biomed. 2010;23(7):745-56.

5. Xu J, Does MD, Gore JC. Dependence of temporal diffusion spectra on microstructural properties of biological tissues. Magnetic Reson Imaging. 2011;29:380–390.

6. Novikov DS, Kiselev VG. Surface-to-volume ratio with oscillating gradients. J Magn Reson. 2011;210:141–145.

7. Rice JE III, Vannucci RC, Brierley JB. The influence of immaturity on hypoxic-ischemic brain damage in the rat. Ann Neurol. 1981;9:131–141.

8. Aggarwal M, Jones MV, Calabresi PA, Mori S, Zhang J. Probing mouse brain microstructure using oscillating gradient diffusion MRI. Magn Reson Med. 2012;67:98–109.

9. Tournier D, Calamante F, Connelly A. Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution. NeuroImage. 2007;35(4):1489-1472.

Figures

Fig. 1: (A-C) Directional anisotropy of apparent diffusion coefficients (ADCs) measured along the x, y, and z axes with PGSE and OGSE (50Hz, 100Hz, and 200Hz) in several mouse brain structures. ΔfDi denotes the rate of ADC change with frequency along the ith direction. (D) The estimated ΔfDx, ΔfDy, ΔfDz in these structures indicated directional anisotropy.

Fig.2: (A) Trace, FA, and colormap images obtained from TDD (left) and conventional PGSE tensors (right) at the same location. (B) Nissl stained section of the mouse brain at a similar slice. (C) PGSE-based FOD map, reconstructed at a region as indicated by a box in (A). The color-coded arrows point to the corresponding structures in (A-C).

Table 1: FA values obtained from the TDD, PGSE and OGSE tensors.

Fig. 3: Trace (top row) and FA (bottom row) images obtained from PGSE (left), OGSE (middle), and the TDD (right) tensors in a representative hypoxia-ischemia injured neonatal mouse brain. The arrows point to the injured ipsilateral hippocampus.



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