Reach Us +1-504-608-2390
Voxel Based Analysis of Diffusion Indices in Patients with Primary Open-Angle Glaucoma Using Tract-Based Spatial Statistics | OMICS International
ISSN: 2155-9570
Journal of Clinical & Experimental Ophthalmology

Like us on:

Make the best use of Scientific Research and information from our 700+ peer reviewed, Open Access Journals that operates with the help of 50,000+ Editorial Board Members and esteemed reviewers and 1000+ Scientific associations in Medical, Clinical, Pharmaceutical, Engineering, Technology and Management Fields.
Meet Inspiring Speakers and Experts at our 3000+ Global Conferenceseries Events with over 600+ Conferences, 1200+ Symposiums and 1200+ Workshops on Medical, Pharma, Engineering, Science, Technology and Business
All submissions of the EM system will be redirected to Online Manuscript Submission System. Authors are requested to submit articles directly to Online Manuscript Submission System of respective journal.

Voxel Based Analysis of Diffusion Indices in Patients with Primary Open-Angle Glaucoma Using Tract-Based Spatial Statistics

Manuel A. Schmidt1*#, Sultan Haider3,4,6#, Angelika Mennecke1, Richard Johnson4, Joachim Hornegger3,4,5, Georg Michelson2, Arnd Doerfler1 and Tobias Engelhorn1
1Department of Neuroradiology, University of Erlangen-Nuremberg, Germany
2Department of Ophthalmology, University of Erlangen-Nuremberg, Germany
3Pattern Recognition Lab, University of Erlangen-Nuremberg, Germany
4Central Institute of Medical Engineering, University of Erlangen-Nuremberg, Germany
5Erlangen Graduate School in Advanced Optical Technologies (SAOT), University of Erlangen-Nuremberg, Germany
6Siemens Healthcare, Erlangen, Germany
#These authors contributed equally to this work and thus share first authorship
Corresponding Author : Manuel A. Schmidt, MD
Department of Neuroradiology, Friedrich-Alexander-University Erlangen-Nuremberg
Schwabachanlage 6, 91054 Erlangen, Germany
Tel: +49 9131 85 44821
Fax: +49 9131 85 36179
E-mail: [email protected]
Received February 05, 2015; Accepted March 26, 2015; Published March 30, 2015
Citation: Schmidt MA, Haider S, Mennecke A, Johnson R, Hornegger J, et al. (2015) Voxel Based Analysis of Diffusion Indices in Patients with Primary Open-Angle Glaucoma Using Tract-Based Spatial Statistics. J Clin Exp Ophthalmol 6:415. doi: 10.4172/2155-9570.1000415
Copyright: © 2015 Schmidt MA, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Related article at
DownloadPubmed DownloadScholar Google

Visit for more related articles at Journal of Clinical & Experimental Ophthalmology


Background and purpose: Diffusion tensor imaging can provide information about microstructural changes of the visual pathway in glaucoma patients. Tract-based spatial statistics (TBSS) has been established for analysis of fractional anisotropy (FA) and has successfully been used for analysis of DTI data of glaucoma patients. However there are certain DTI indices such as mean diffusivity (MD), radial diffusivity (RD), axial diffusivity (AD) and mode of anisotropy (MO), which can also be analyzed using TBSS. Hence, we tested whether TBSS is useful for analysis of non-FA data in glaucoma patients and for identification of microstructural changes in patients with primary openangle glaucoma characterized by these non-FA DTI indices.
Material and methods: Our study comprised 46 patients with primary open-angle glaucoma (POAG) and a healthy, age-matched control group of 23 subjects. TBSS was used for analysis of FA as well as MO, MD, RD and AD. Results: For an automated analysis of non-FA diffusion indices we could establish a workflow using TBSS in addition to its primary purpose, the analysis of FA. Our results reveal clusters of voxels with overlapping decreased FA and increased RD and MD in the optic radiation as well as in visual association tracts.
Conclusion: TBSS is a useful tool not only for analysis of FA but also of non-FA diffusion indices in glaucoma patients. Overlapping decreased FA and increased MD and RD can be found in POAG patients not only in the optic radiation but also in visual association tracts, suggestive for significant neurodegeneration in POAG.

Tract-based spatial statistics; Primary open-angle glaucoma; DTI; Fractional anisotropy; Diffusion indices; Mode of anisotropy; Radial diffusivity; Mean diffusivity; Axial diffusivity
Magnetic resonance diffusion imaging has been used to analyze the visual pathway in glaucoma [1-3]. A tensor model fit to diffusion MRI data provides scalar measurements to perform voxel wise analysis of white matter properties [4]. As glaucoma is more and more comprehended as a neurodegenerative disease, advanced imaging post processing algorithms have been developed for sophisticated analysis of diffusion indices, primarily fractional anisotropy (FA) [5-7].
Tract-based spatial statistics (TBSS) can help to solve crucial issues of spatial alignment, allowing non-biased whole brain cross-subject statistical analysis [8]. In addition to FA, several diffusion indices have been described to be associated with distinct microstructural changes of white-matter tracts: mean diffusivity (MD) is thought to represent cell membrane degradation [9] and demyelination has been shown to be represented by radial diffusivity (RD) [10]. Non-FA diffusion indices have been analyzed with statistical parametric mapping [11]. However there is only one study [12] that conducted a TBSS based analysis; this study comprised only a small amount of patients and performed no analysis of MD and another important diffusion measure, mode of anisotropy.
Thus, the aim of this study was to establish an automated workflow for whole brain analysis of multiple diffusion indices in glaucoma patients and to identify and assess microstructural changes that are represented by changes of diffusion measures compared to healthy, age-matched control individuals.
Materials and Methods
This retrospective study included 69 subjects comprising 46 severely affected patients with primary open-angle glaucoma (POAG ,intraocular pressure of both eyes prior to treatment ≥ 22 mmHg; mean age=63.75 years ± 10.79) and a healthy, age-matched control group of 23 subjects (CONT, mean age=59.54 years ± 14.06). The individuals of the CONT group underwent full ophthalmological examination to exclude an intraocular pressure ≥ 22 mmHg, optic nerve head atrophy and visual disturbances. Patients’ characteristics are summarized in Table 1. The Clinical Investigation Ethics Committee of the University of Erlangen-Nuremberg approved the study protocol and the research was conducted in accordance with the Declaration of Helsinki and informed consent was obtained from all subjects.
MRI Sequence
A 3T high-field scanner (Magnetom Tim Trio, Siemens Healthcare AG, Erlangen, Germany) with gradient field strength up to 45 mT/m (72 mT/m effective) was used. DTI was performed in the axial plane with 4 mm slice thickness using a single-shot, spin echo, echo planar imaging (EPI) diffusion tensor sequence (TR=3400 ms, TE=93 ms, FoV=230 × 230 mm2, acquisition matrix size=256 × 256 reconstructed to 512 × 512, number of signal averages=7, partial Fourier acquisition=60%). Diffusion weighting was carried out with a maximal b-factor of 1000 s/mm2 along 15 icosahedral directions complemented by one scan with b=0.
Image processing and TBSS analysis
The entire image processing and analysis workflow is pictured in Figure 1 and consisted of four steps: (1) Image preprocessing to set up DTI data for tract-based spatial statistics (TBSS), (2) run TBSS, (3) voxel wise statistical analysis and (4) examination and review of results.
Voxel wise statistical analysis of FA was carried out using TBSS [5], part of FSL, a comprehensive library of analysis tools for FMRI, MRI and DTI brain imaging data [13,14]. For Preprocessing DICOM images were converted to NIfTI files (neuroimaging informatics technology initiative) using dcm2nii from the MRIcron-package (
The resulting images were then preprocessed for TBSS analysis, i.e. corrected for eddy currents and brain extracted with bet [15]. We chose a bet2-threshold of 0.2 to remove all non-brain tissue. Next a diffusion tensor model was fit at each voxel to extract the FA, MO, MD and AD maps using DTIFIT. All of these preprocessing steps were carried out using FMRIB diffusion toolbox (FDT) which is also part of FSL ( [14].
During registration, mean FA skeletons were created using the registration to a group wise template (GW) as proposed in [7]. GW was used with 5 affine and 15 non-linear, gradually converging registration-steps and the mean FA skeleton image was thresholded to 0.2. Within the next step, the final transforms found in the previous stage were applied to all subjects to bring them into standard space.
After that, the mean FA image was skeletonized and regarding every individual FA image, the maximum FA values perpendicular to the mean skeleton were projected on the mean skeleton.
This step resulted in a 4D skeletonized FA image, comprising the individual skeletonized FA images. It was used to feed into voxel wise statistics and to observe which FA skeleton voxels are significantly different between the two groups of subjects.
Voxelwise statistics was performed using a two sample t-test with the permutation function randomise included in FSL, number of permutations set to 5000 and threshold-free cluster enhancement, fully corrected for multiple comparisons across space, with p<0.05 resulting and 2 contrasts (CONT>POAG and CONT<POAG) [13,16].
Additionally to FA, TBSS evaluation was applied by this means to non-FA diffusion indices, obtained in the preprocessing step by DTIFIT. This included mean diffusivity (MD = λ 123/3 ), axial diffusivity (AD=λ1), radial diffusivity (RD = λ 23 2 , calculated with fslmaths, a FSL command line utility) and mode of anisotropy (MO). MO describes the shape of the diffusion tensor. It ranges from −1 to +1 and gives complementary information to FA, as the shape ranges from more planar (λ1 ~ λ2 ~ λ3, e.g. in areas of crossing fibers) to linear (λ12 ~ λ3, e.g. areas with one dominating fiber orientation) [16,17].
The data was processed referring to the calculated FA skeleton from group-wise registration. Voxelwise statistics with randomise (FSL's tool for nonparametric permutation inference) were carried out in same manner as for FA (see above).
JHU white-matter tractography atlas was used for annotation [18,19].
Two experienced neuroradiologists performed evaluation of the results.
We found clusters of voxels with decreased FA and MO and increased RD and MD in and outside the visual pathway in POAG compared to controls.
POAG patients exhibited decreased FA in the optic radiation bilaterally as well as in the right corticospinal tract, left anterior thalamic radiation, the corpus callosum and in association tracts, specifically the inferior fronto-occipital fasciculus (bilateral), the right superior and the inferior longitudinal fasciculus (bilateral) and the uncinate fasciculus (bilateral). MO was significantly decreased in the right corticospinal tract and superior corona radiata, right internal capsule and anterior thalamic radiation and in the right uncinate fasciculus and inferior fronto-occipital fasciculus. There were no clusters of voxels with increased FA or MO in POAG patients compared to healthy controls.
RD and MD were increased in nearly all association tracts as well as in projection and commissural fibers. Cluster size and spatial distribution matched to a great extend for RD and MD (Figure 2). For contrast number 2, there were no voxels with decreased RD or MD in POAG patients.
We could not observe significant changes in axial diffusivity. Results are summarized in Table 1.
For the optic radiation we found decreased FA and overlapping increased MD and RD respectively (Figure 3). This was less pronounced for other projection fibers but could also be found in important association tracts related to the visual system.
The aim of this study was to evaluate an automatized approach for analyzing FA and non-FA diffusion indices in POAG patients using TBSS. Therefore we optimized the TBSS pipeline for analysis of patients with a potentially neurodegenerative disease, as it has previously been shown, that brain atrophy has to be taken into account concerning the registration process [7,16,20].
Our results reveal decreased FA in POAG patients in the bilateral optic radiation as well as in association tracts linked to the visual pathway. Furthermore, we show that significant and widespread alteration of diffusion indices, especially MD and RD, in POAG patients compared to healthy control individuals is not limited to the visual pathway but can also be found in projection, association and commissural fibers.
Interestingly, overlapping of decreased FA and increased MD and RD could be found in the bilateral optic radiation as well as in important association fibers for processing of visual information indicating advanced neurodegeneration in POAG patients.
Decreased FA could be found in the bilateral optic radiation in POAG patients, which has previously been described [6,7]. Moreover, we found decreased FA in medial prefrontal association fibers, predominantly the superior and inferior longitudinal fasciculus, the inferior fronto-occipital fasciculus and uncinate fasciculus. The inferior fronto-occipital fasciculus connects the frontal region with inferior and dorsal parietal and occipital areas as an intrahemispheric association fiber running in sagittal orientation. The superior longitudinal fasciculus connects frontotemporal and frontoparietal regions. Both have been described to play a role in visual spatial processing [21] and whilst being affected by degeneration through decreased FA may contribute to visual impairment of patients.
The inferior longitudinal fasciculus runs from the frontal to the ipsilateral occipital lobe and interruption may lead to visual neglect and has been reported to be particularly involved in visual memory [22].
Resection of the uncinate fasciculus in glioma patients has led to a loss of performance in famous face naming [23] and it has been reported to be associated with visual associative learning suggesting that degeneration of the uncinate fasciculus may also be somewhat involved in the pathogenesis of glaucoma.
Our results are in line with previous reports that described reduced gray matter in POAG patients in the calcarine fissure, the postcentral gyrus, the superior frontal gyrus and the inferior frontal gyrus [24].
We found decreased mode of anisotropy values in the right corticospinal tract as well as in the superior corona radiata and the internal capsule and the right anterior thalamic radiation. Additionally, MO was significantly lower partly in the right inferior fronto-occipital and unicinate fasciculus and in the corpus callosum. The mode of anisotropy reflects the shape of the diffusion tensor, ranging from linear to planar. Lower mode of anisotropy indicates more planar diffusion, e.g. in areas of crossing fibers [17]. Reduced FA has been reported in such areas of crossing fibers and may in those cases not be due to axonal degradation but reflect disorganization of white-matter tracts [16]. The absence of decreased MO in the optic radiation and at least partly, in association tracts of the visual pathway, where we found markedly decreased FA values, may thus represent “true” FA reduction due to loss of axonal integrity.
While FA describes overall white matter integrity [9], increased MD (also known as ADC, apparent diffusion coefficient) has been linked to tissue breakdown due to cell membrane degeneration and is associated with increased brain water content [9,25]. We found increased MD in nearly all association tracts as well as in projection and commissural fibers.
Interestingly, overlapping of increased MD and decreased FA was markedly pronounced in the optic radiation and in those association tracts, which are strongly connected to the visual system, like the inferior fronto-occipital fasciculus and superior longitudinal fasciculus (Figure 3).
In a cohort of patients with Alzheimer’s disease, overlapping changes in FA and MD have been reported to be due to gross tissue loss. This might also apply for glaucoma patients and may be a better indicator of white-matter degeneration as decreased FA alone.
Additionally, decrease in RD followed the decrease of MD to a great extent. Decreased RD has been associated with loss of myelin in animal studies [10] and has also been described in patients with Alzheimer’s disease [26,27] and in glaucoma [28]. Thus the overlapping of increased MD and RD with decreased FA is suggestive for advanced neurodegeneration in the optic radiation as well as in association tracts necessary for adequate processing of visual information in glaucoma patients.
For AD, there were also clusters of voxels with increased values in POAG patients, however, those clusters scarcely failed to reach statistical significance. This may be due to dynamic changes, which have been observed regarding axial diffusivity [29] (Table 2).
Glaucoma has commonly been understood as a disease affecting primarily the optic nerve, i.e. the third neuron of the visual pathway. However, as diffusion tensor imaging is more and more used not just for evaluation of the extracranial but also of the intracranial part of the visual pathway [1,30], the idea of a solely affected optic nerve in glaucoma has to be questioned. It has already been proposed that at least in some glaucoma patients a primary degeneration of the fourth neuron of the visual pathway, i.e. the optic radiation, might strongly contribute to disease pathogenesis rather than a secondary affection of the fourth neuron due to transsynaptic degeneration [31].
Our results suggest that neurodegeneration in POAG patients goes far beyond degeneration of the third and fourth neuron of the visual pathway. Decreased FA with simultaneously increased MD and RD is a strong indicator for advanced neurodegeneration in primary openangle glaucoma. The described findings lead to the hypothesis that degeneration of visual association tracts is secondary to the initial degeneration of visual projection fibers, first to fourth neuron, rather than degeneration of the fourth neuron being secondary due to optic nerve atrophy. This notion would be in line with previous reports of an initial injury of the fourth neuron in primary open-angle glaucoma [31].
Disease severity in glaucoma is estimated by a profound ophthalmological examination, e.g. optic nerve head morphology or retinal nerve fiber layer thickness, in clinical practice. DTI of the central visual pathway may contribute to a better overall estimation of disease severity in POAG patients. At least for the optic radiation changes in FA, MD and RD have been shown to correlate with glaucoma severity [28].
Further studies are required to assess changes of these diffusion parameters in visual association tracts with regard to disease severity. Given that an impaired axonal integrity in tracts crucial for processing of visual information leads to an impaired visual field over time, it even seems possible to use DTI of the central visual pathway for an early detection of disease progress, before degeneration has reached the optic nerve or the retina, respectively.
It seems that, at least in certain glaucoma subtypes, i.e. primary open-angle glaucoma, the classical understanding of the pathogenesis has to be rethought, as of glaucoma being a complex neurodegenerative disease.
TBSS is a useful tool not only for analysis of FA but also non-FA diffusion indices in glaucoma patients. Concurrent analysis of different DTI derived indices may provide a better pathophysiological understanding of primary open-angle glaucoma, as overlapping alterations of DTI indices suggest POAG to be a complex neurodegenerative disorder.


Tables and Figures at a glance

Table icon Table icon
Table 1 Table 2


Figures at a glance

Figure Figure Figure
Figure 1 Figure 2 Figure 3
Select your language of interest to view the total content in your interested language
Post your comment

Share This Article

Article Usage

  • Total views: 12300
  • [From(publication date):
    specialissue-2015 - Sep 17, 2019]
  • Breakdown by view type
  • HTML page views : 8505
  • PDF downloads : 3795