A New Method to Construct and Validate an Optic Radiation Template using Probabilistic Tractography
Chenyu Wang1,2, Alexander Klistorner1,3,4, Linda Ly1,2, Ruth A Oliver1,2, and Michael H Barnett1,2

1Sydney Neuroimaging Analysis Centre, Camperdown, Sydney, Australia, 2Brain and Mind Centre, University of Sydney, Camperdown, Sydney, Australia, 3Ophthalmology, Save Sight Institute, University of Sydney, Sydney, Australia, 4Australian School of Advanced Medicine, Macquarie University, Sydney, Australia

Synopsis

In-vivo delineation the white matter (WM) fibre bundles with diffusion-weighted imaging (DWI) and tractography facilitates the study of tract-specific damage due to neurological disease. However, this approach is limited by suboptimal or absent DWI datasets. Accurate atlas based white matter segmentation may provide an alternative method that can be broadly applied to existing neuroimaging datasets . Using the optic radiation as an example, we propose a new framework for the mapping and validation of white matter tracts in healthy subjects using probabilistic tractography.

Purpose

To construct and validate an atlas based template of the optic radiation.

Methods

34 healthy participants, 24 female and 9 male, mean age 37.18 (SD=14.62), were divided into two groups: Group A (19 subjects) for constructing optic radiation (OR) template and Group B (15 Subjects) for validating the OR template construction. MRI data was acquired for all participants from a 3.0T GE MR750 scanner (General Electric, Milwaukee, WI, USA), using an 8 channel head coil. For each exam, whole brain 64-directions diffusion weighted imaging was acquired with 2 mm isotropic acquisition matrix (TR/TE=8325/86 ms, b = 1000 s/mm2, number of b0s = 2); Additionally, Group A was acquired with Sagittal IR-FSPGR (TR/TE/TI =7.2/2.7/450 ms, 1 mm isotropic acquisition matrix, FOV = 256mm), and Group B was acquired with Axial IR-FSPGR (TR/TE/TI =8.2/3.2/900 ms, 1 mm isotropic acquisition matrix, FOV = 256mm). All diffusion weighted imaging (DWI) was corrected for motion, eddy-current and EPI susceptibility distortion, prior to the tensor reconstruction and T1 structure imaging co-registration in MrDiffusion package (MrVista, Stanford University, http://web.stanford.edu/group/vista/cgi-bin/wiki/index.php/MrVista). Seeding points at the lateral geniculate nucleus (LGN) and calcarine sulcus were used to reconstruct OR using probabilistic tractography with Contrack1 (Figure 1). Bilateral OR was reconstructed for each participant, and constrained by individual white matter masks derived from T1 using SIENAX (FMRIB Software Library; www.fmrib.ox.ac.uk/fsl). ICBM 2009a nonlinear Asymmetric template (http://www.bic.mni.mcgill.ca/ServicesAtlases/ICBM152NLin2009) was used for mapping the final OR template. After non-brain tissue was removed (Brain Extraction Tool, FMRIB software Library), all T1 images were co-registered non-linearly to ICBM 2009a template using ANTS (Advanced Normalization Tools, http://picsl.upenn.edu/software/ants) to obtain the transformation maps. Individual fibre-based ORs were converted to binary ROI, mapped into ICBM space and averaged to construct the OR template. The initial OR template derived from group A was presented as a probability map. Group B were used to validate the OR template and determine the optimal probability threshold for obtaining the binarized OR template. Briefly, the OR probability map was thresholded from 10% to 90% with 1% increments; all thresholded OR templates were then mapped on to Group B individually using the inversed deformation maps, and compared with probabilistic tractography based individual ORs using the Dice Similarity Coefficient (DSC). Finally the threshold that produced the largest mean DSC in group B was used to obtain the final OR template.

Results

Reconstruction of the OR using probabilistic tractography was successful in all subjects, except for a unilateral OR from one subject whose imaging was degraded by imaging artefact. The OR template, constructed from Group A, is presented as a probability map in Figure 2a. Topographical conformity with a histologically determined atlas-based OR2 was confirmed visually (Figure 2c). Meyer loop, in particular, was well defined. When the OR probability map was applied to a validating dataset (Group B), the highest mean DSC (SD), 0.703 (0.047) and 0.716 (0.067), corresponded with a probability threshold level of 40% and 44% for the left and right OR respectively (table 1). The OR template, determined by these thresholds, was compared to individual probabilistic tract-based OR reconstructions (Figure 3). While excellent correspondence between the two methods was found, a relative under-estimation of the Meyer loop component of the OR was noted in tract-based reconstructions (figure 4).

Discussion

Nonlinear registration of a tractography generated OR template facilitates accurate estimation of the tract in individual cases without DWI. In our study, a threshold of ~40% provides the best approximation of template-based OR to individually defined tract-based OR. In addition, template-based tractography provides better representation of Meyer loop in some subject; however, care should be taken in interpreting the size and position of the posterior part of the OR.

Conclusion

We propose a new pipeline to generate and validate an OR template using probabilistic tractography. We have demonstrated topographic concordance with a histological atlas and robust estimation when compared to individual tractography-based reconstruction. The variance between template based and tractography based OR reconstructions is most apparent in the vicinity of the Meyer loops and posterior part of the OR. Future work should include validation of our work in a larger dataset derived from multiple sites.

Acknowledgements

No acknowledgement found.

References

1. Sherbondy, A.J., Dougherty RF, Ben-Shachar M, Napel S & Wandell BA. (2008). ConTrack: Finding the most likely pathways between brain regions using diffusion tractography. Journal of Vision, 1–3.

2. Bürgel, U., Schormann, T., Schleicher, A., & Zilles, K. (1999). Mapping of histologically identified long fiber tracts in human cerebral hemispheres to the MRI volume of a reference brain: position and spatial variability of the optic radiation. NeuroImage, 10(5), 489–99.

Figures

Figure 1. Probabilistic tractography based optic radiation reconstruction pipeline

Figure 2. (a) OR tempate constructed from 19 healthy controls, presented as a probability map; (b) OR template thresholded at L:40%, R:44% from (a); (c) OR probability map from Jülich histological atlas2.

Figure 3. Comparison of template based OR reconstruction (Red) and Tract based OR reconstruction (Blue) for all participants in Group B. Left and right OR templates thresholded at 40% and 44% from OR probabilistic template.

Figure 4. Comparison of template based to tract based OR. Red-yellow colour indicates areas of likely overestimation and blue-lightblue areas of likely underestimation by the template based method, when compared to individual tractography methods. Greater intensity indicates higher degree overestimation/underestimation.

Table 1. Individual DSC in group B when the thresholds are applied, achieving the maximum group mean DSC.



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