Thresholding to Improve the Specificity of High Spatial and Angular Resolution In Vivo Diffusion-Weighted Tractography to Estimate Brain Stem Connectivity.
Matthew Hey1, Luis Colon-Perez2, William Triplett3, David Fitzgerald4, and Thomas Mareci5

1University of Florida, Gainesville, FL, United States, 2Department of Psychiatry, University of Florida, Gainesville, FL, United States, 3Department of Physical Therapy, University of Florida, Gainesville, FL, United States, 4Department of Neurology, University of Florida College of Medicine, Gainesville, FL, United States, 5Department of Biochemistry and Molecular Biology, University of Florida College of Medicine, Gainesville, FL, United States

Synopsis

The spatial resolution of diffusion-weighted (DWI) images limits the white matter streamline fiber tracks, which can be followed in the brain stem. To address this issue, we introduce a high spatial resolution protocol and the use of a threshold to limit the false positive in streamline track density maps by requiring that a minimum amount of fibers pass through a voxel. This provides increased accuracy in the visualization of streamlines connecting specific regions of the brain stem and may allow the recognition of structural abnormalities due to neurological diseases.

Introduction

The brain stem is an essential component in the central nervous system involved in the regulation of sleep-wake cycle, cardiac functions, and respiratory functions1. DWI and tractography provide an ideal tool to understand changes in the brain stem structure associated with neurological disorders2 since DWI and tractography can provide a comprehensive look into the neural connections of the brain in vivo and with little risk to subject3. One of the biggest challenges facing current methods of tractography is false positives resulting from noise and low spatial resolution of in vivo measurements. Thresholding tractography maps to limit the amount of false positives in the data set may be extremely useful in the diagnosis of neurological diseases. The goal of the protocol was to optimize the acquisition sequence for the brain stem, which can to be performed in less than 20 minutes, without sacrificing the quality of the data. Then limit the number of false positive through the application of a threshold filtering of the calculated streamline tracks.

Methods

Eleven healthy control participants, between the ages of 20 and 40 years, were recruited with no prior history of neurological disorders. Participants were scanned using a 3.0 T Phillips system with a 32-channel head coil. A high-angular-resolution diffusion imaging (HARDI) sequence was acquired. A total 71 directions were obtained; 1 with b = 0, 6 with b = 100 s/mm2 and 64 with b = 1000 s/mm2 using a model electrostatic repulsion to determine the diffusion gradient directions4. The acquired image resolution was 1.7 x 1.7 x 1.7 mm3, TR/TE = 5000/86 ms, and Δ/δ = 42.1/9.1 ms. These images were then interpolated by a factor of 2 using cubic convolution with the CONGRID function in IDL [MH1] for a final image resolution of 0.85 mm isotropic. The fractional anisotropy (FA) and mean diffusivity (MD) were reconstructed from the HARDI data with in-house software written in IDL (Exelis Visual Information Systems, Boulder, CO). The thalamus and the pons were hand drawn using ITKSNAP using the FA and MD maps. The mixture of Wisharts (MOW) model was utilized for reconstruction of the fiber orientations5; this model characterizes the displacement probability function allowing the estimation of multiple fibers per voxel. MOW was calculated four separate times per data set each with a different probability displacement radius of: 13, 16, 18, and 20 microns. Utilizing the maxima of the displacement probability function as fiber orientations, deterministic tractography was performed with 125 seeds per voxel, fiber step size of one-half voxel, and no step-to-step angular deviations greater than 50°. The streamlines connecting the pons and thalamus ROIs was generated using in-house C based software. The streamlines connecting the 2 node-network (i.e. pons and thalamus) was converted into a tract density mask with each voxel containing an intensity value corresponding to the number of streamlines passing through that voxel. Finally a threshold was applied to the density maps for voxels containing less than 1%, 3%, 5%, and 10% of the maximum number of streamline tracts.

Results

Eleven healthy control participants were successfully scanned and the connectivity between the pons and thalamus was estimated with deterministic tractography. The total scan time was 16 minutes and the isotropic resolution of 1.7 x 1.7 x 1.7 mm3 represents a 38.5% reduction in voxel volume compared to traditional in vivo acquisition of 2 x 2 x 2 mm3 with comparable acquisition time. Figure 1 shows the effect of applying a threshold to the density map images. Without a threshold, a large number of voxels exhibit streamlines making to discrimination of any apparent structure difficult. As the threshold increases, small bundles emerge displaying a large number of streamlines. Figure 2 shows a structure that lies within the expected area associated as the locus coeruleus and represents the first time the locus coeruleus has been segmented in vivo with tractography.

Discussion

One of the greatest problems with tractography is the number of false positives that result from the inherent low spatial resolution in DWI data6. These false positives make data interpretation difficult. By optimizing the spatial resolution and thresholding the number of fibers in a voxel, we were able to eliminate extraneous fibers. With this approach, tractography can add in the study anatomical abnormalities due to neurological diseases, such as PTSD, Alzheimer’s, or Parkinson’s.

Acknowledgements

This work was supported in part by Department of Defense USAMRMC/TATRC grant Correlating Sleep Disturbances and Damaged White Matter Tracts in the Brainstem using Diffusion Weighted Imaging (Contract #W81XH-11-1-0454). Data acquisition was performed in the McKnight Brain Institute at the National High Magnetic Field Laboratory’s AMRIS Facility, supported by National Science Foundation Cooperative Agreement No. DMR-1157490, the State of Florida, and the U.S. Department of Energy.

References

1. Nicholls, J. G. & Paton, J. F. Brainstem: neural networks vital for life. Philos Trans R Soc Lond B Biol Sci 364, 2447-2451, doi:10.1098/rstb.2009.0064 (2009).

2. Lazar, M. Mapping brain anatomical connectivity using white matter tractography. NMR in Biomedicine 23, 821-835, doi:10.1002/nbm.1579 (2010).

3. Hagmann, P. et al. MR connectomics: Principles and challenges. Journal of Neuroscience Methods 194, 34-45, doi:10.1016/j.jneumeth.2010.01.014 (2010).

4. Jones, D. K., Horsfield, M. A. & Simmons, A. Optimal strategies for measuring diffusion in anisotropic systems by magnetic resonance imaging. Magn Reson Med 42, 515-525 (1999).

5. Jian, B., Vemuri, B. C., Ozarslan, E., Carney, P. R. & Mareci, T. H. A novel tensor distribution model for the diffusion-weighted MR signal. NeuroImage 37, 164-176 (2007).

6. Thomas, C. et al. in Proc Natl Acad Sci U S A Vol. 111 16574-16579 (2014).

Figures

Figure 1. Density maps of connecting streamlines between the pons and thalamus. Top row transverse and bottom row sagittal views. Images from left to right: Increasing threshold from 0%, 1%, 3%, 5%, and 10% of the maximum number of streamlines where yellow voxels contain a large number of streamlines and red the least.

Figure 2. Effect on threshold in streamline density maps in expanded view of transverse image. Left, streamline density maps without threshold, and right, a streamline density map with a 5% threshold. The blue arrows indication the location of the locus coeruleus.



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