Schizophrenia Patients Demonstrate Both Inter-Voxel Level and Intra-Voxel Level White Matter Alterations

Abstract

Fractional anisotropy (FA) and mean diffusivity (MD) are the most frequently used metrics to investigate white matter impairments in mental disorders. However, these two metrics are derived from intra-voxel analyses and only reflect the diffusion properties solely within the voxel unit. Local diffusion homogeneity (LDH) is a newly developed inter-voxel metric which quantifies the local coherence of water molecule diffusion in a model-free manner. In this study, 94 schizophrenia patients and 91 sex- and age-matched healthy controls underwent diffusion tensor imaging (DTI) examinations. White matter integrity was assessed by FA, MD and LDH. Group differences in these metrics were compared using tract-based spatial statistics (TBSS). Compared with healthy controls, schizophrenia patients exhibited reduced FA and increased MD in the corpus callosum, cingulum, internal capsule, fornix and widespread superficial white matter in the frontal, parietal, occipital and temporal lobes. We also found decreased LDH in the corpus callosum, cingulum, internal capsule and fornix in schizophrenia. Our findings suggest that both intra-voxel and inter-voxel diffusion metrics are able to detect impairments in the anisotropic white matter regions, and intra-voxel diffusion metrics could detect additional impairments in the widespread isotropic white matter regions in schizophrenia.

About the Authors:

Chuanjun Zhuo

Contributed equally to this work with: Chuanjun Zhuo, Xiaolei Ma, Hongru Qu

Affiliations Tianjin Anding Hospital, Tianjin Mental Health Center, Tianjin, China, Department of Radiology and Tianjin Key Laboratory of Functional Imaging, Tianjin Medical University General Hospital, Tianjin, 300052, China, Tianjin Anning Hospital, Tianjin, 300300, China, Wenzhou seventh people's Hospital, Wenzhou, 325000, China

Xiaolei Ma

Contributed equally to this work with: Chuanjun Zhuo, Xiaolei Ma, Hongru Qu

Affiliation: Tianjin Anning Hospital, Tianjin, 300300, China

Hongru Qu

Contributed equally to this work with: Chuanjun Zhuo, Xiaolei Ma, Hongru Qu

Affiliation: Tianjin Anning Hospital, Tianjin, 300300, China

Lina Wang

Affiliation: Tianjin Anding Hospital, Tianjin Mental Health Center, Tianjin, China

Feng Jia

Affiliation: Tianjin Anding Hospital, Tianjin Mental Health Center, Tianjin, China

Chunli Wang

* E-mail: [email protected]

Affiliation: Tianjin Anding Hospital, Tianjin Mental Health Center, Tianjin, China

Introduction

Diffusion tensor imaging (DTI) is commonly used to investigate white matter (WM) alterations in the human brain [1,2]. Based on the diffusion tensor model, fractional anisotropy (FA) and mean diffusivity (MD) are the most frequently used metrics to characterize the diffusion of water molecules [3,4]. FA represents a normalized ratio of diffusion directionality, while MD quantifies the bulk mobility of water molecules; thus, these indices provide information on white matter integrity from different perspectives. Although FA and MD are the two most commonly used metrics for evaluating alterations in white matter integrity, they are intra-voxel diffusion metrics and only reflect diffusion properties solely within the voxel unit [3,5,6].

Several studies have reported that inter-voxel diffusion metrics are informative and useful [5–7]. These inter-voxel metrics are able to reveal significant region-specific white matter alterations in some neuropsychical diseases and provide useful information for computer-aided diagnosis, although inter-voxel and intra-voxel metrics have shown inconsistent findings [7–9].

In a prior study, a novel inter-voxel metric was proposed to quantify the local coherence of water molecule diffusion in a model-free manner, which is referred to as the local diffusion homogeneity (LDH) [9]. Because the profile of water molecule diffusion is related to the underlying white matter micro-structural properties, the LDH is supposed to reflect the local coherence of fiber orientations, myelination, diameter or density along each white matter tract. Some of these micro-structural factors cannot be interpreted by the traditional diffusion metrics, such as FA and MD. Therefore, LDH can serve as a complementary marker for investigating white matter alterations of the brain and can provide additional information for understanding the pathological mechanisms of neuropsychical diseases [9].

Many previous studies have confirmed that alterations of FA and MD can be found in multiple white matter regions in patients with schizophrenia [10–16]. Together, the findings suggest that impairments of white matter integrity are pathological features of schizophrenia. However, as mentioned above, all of these studies focused on the intra-voxel metrics. Therefore, using the inter-voxel metric to explore white matter alterations in patients with schizophrenia may provide additional information on the neural mechanisms of schizophrenia.

In the current study, we reconstructed the white matter skeleton of the brain using tract-based spatial statistics (TBSS) [17] and compared differences in the inter-voxel (LDH) and the intra-voxel (FA and MD) metrics within the skeleton between schizophrenia patients (n = 94) and healthy controls (n = 91). We hypothesized that both intra-voxel and inter-voxel metrics are able to detect white matter alterations in schizophrenia patients and that the distribution of abnormal white matter regions detected by the two types of metrics may be different.

Methods

Subjects

A total of 106 schizophrenia patients and 94 healthy subjects were enrolled in the present study. This study was approved by the Ethics Committee of Tianjin Medical University General Hospital, and written informed consent was obtained from each subject before study enrollment. Diagnosis for each schizophrenia patient was determined by the consensus of two clinical psychiatrists using the Structured Clinical Interview for DSM-IV (SCID). Inclusion criteria were age (16–60 years) and right-handedness. Exclusion criteria for all subjects were MRI contraindications, pregnancy, and histories of systemic medical illness, central nervous system disorder and head trauma, and substance abuse within the last 3 months or lifetime history of substance abuse or dependence. Additional exclusion criteria for healthy controls included a history of psychiatric disease and first-degree relatives with a psychotic disorder. After image quality assessment slice by slice by two professional radiologists, 12 patients and 3 controls were excluded due to poor quality of imaging data. The final sample included 94 schizophrenia patients and 91 healthy controls. Psychotic symptoms were quantified with the Positive and Negative Syndrome Scale (PANSS) [18].

MRI Data Acquisition

MRI was performed using a 3.0-Tesla MR system (Discovery MR750, General Electric, Milwaukee, WI, USA). Tight but comfortable foam padding was used to minimize head motion, and earplugs were used to reduce scanner noise. DTI data were acquired by a spin-echo single-shot echo planar imaging (EPI) sequence with the following parameters: repetition time = 5800ms; echo time = 77ms; matrix = 128 × 128; field of view = 256 mm × 256 mm; in-plane resolution = 2 mm × 2 mm; slice thickness = 3mm without gap; 48 axial slices; b = 1000 s/mm² and 25 encoding diffusion directions; 10 non-diffusion-weighted images (b = 0 s/mm²). The total acquisition time for DKI was 5 min and 54 s.

Calculation of FA, MD and LDH

Eddy current-induced distortion and motion artifacts in the DTI dataset were corrected using affine alignment of each diffusion-weighted image to the b = 0 image. After skull-stripping, we calculated the diffusion tensor matrix for each subject in the native space. Next, three pairs of eigenvalues and eigenvectors were obtained from the diffusion tensor matrix. FA and MD images were calculated based on the three eigenvalues. We also calculated the LDH images (the pre-defined neighborhood was 27 voxels) according to a prior study [9]. All of the procedures described above were implemented by the PANDA toolbox [19].

Tract-Based Spatial Statistics

The following steps were adopted for the TBSS analysis [17] using FMRIB’s diffusion toolbox (FSL 4.0, http://www.fmrib.ox.ac.uk/fsl). All subjects’ FA images were aligned to a FMRIB-58 FA template in Montreal Neurological Institute (MNI) space using a non-linear registration algorithm. After transformation into MNI space, a mean FA image was created and thinned to generate a mean FA skeleton of the white matter tracts. Each subject’s FA image was then projected onto the skeleton via filling the mean FA skeleton with FA values from the nearest relevant tract center by searching perpendicular to the local skeleton structure for maximum FA value. The registration and projection information derived from the FA analysis were then applied to MD and LDH images of each subject to ensure an exact spatial correspondence of the different parameters.

Statistical Analysis

A Chi-square test of Pearson and a t-test of Student were used to test the group differences in sex and age, respectively. Voxel-wise statistical analysis across subjects on the skeleton space was carried out using a permutation-based inference tool for nonparametric statistic (“randomize”, part of FSL). Group comparisons between schizophrenia patients and healthy controls were performed using a general linear model with age and gender as covariates of no interest. The mean FA skeleton was used as a mask, and the number of permutations was set to 5000. The significance threshold was determined with a P < 0.01 (two-tailed) after correcting for family-wise error (FWE) using the threshold-free cluster enhancement (TFCE) [20] option in FSL. Fiber tracts corresponding to the clusters were identified with reference to Johns Hopkins University ICBM-DTI-81 White-Matter Labels provided in FSL toolbox.

Results

Demographic and Clinical Characteristics

The demographic and clinical data of the subjects are summarized in Table 1. There were no significant group differences in sex (Chi-square test, χ² = 1.912, P = 0.167) and age (two-sample t test, t = 0.067, P = 0.947). 85 patients were receiving atypical antipsychotic medications during the MRI examinations and the other 9 patients had never received any medications. For schizophrenia patients, the mean antipsychotic dosage of chlorpromazine equivalents was 462.5 ± 346.8 mg/d; the mean duration of illness was 123.1 ± 98.6 months; the mean scores of PANSS positive sub-scale and negative sub-scale were 17.1 ± 7.8 and 19.5 ± 8.2.

[Figure omitted. See PDF.]

Table 1. Demographic and Clinical Characteristics of the Sample.

https://doi.org/10.1371/journal.pone.0162656.t001

Intergroup Differences in FA and MD

The white matter regions with significant intergroup differences (P < 0.01, two-tailed, FWE corrected) in FA and MD are shown in Figs 1 and 2. Schizophrenia patients demonstrated reduced FA and increased MD compared to healthy controls. The distribution patterns of white matter regions with reduced FA and increased MD were similar. The decreased FA and increased MD were mainly located in the corpus callosum, cingulum, internal capsule, fornix and widespread superficial white matter in the frontal, parietal, occipital and temporal lobes.

[Figure omitted. See PDF.]

Fig 1. White matter regions with significant differences in FA between patients and controls (P < 0.01, two-tailed, FWE corrected).

FA, fractional anisotropy.

https://doi.org/10.1371/journal.pone.0162656.g001

[Figure omitted. See PDF.]

Fig 2. White matter regions with significant differences in MD between patients and controls (P < 0.01, two-tailed, FWE corrected).

MD, mean diffusivity.

https://doi.org/10.1371/journal.pone.0162656.g002

Intergroup Differences in LDH

The white matter regions with significant intergroup differences (P < 0.01, two-tailed, FWE corrected) in LDH are shown in Fig 3. Compared with healthy controls, schizophrenia patients showed decreased LDH mainly in the corpus callosum, cingulum, internal capsule and fornix.

[Figure omitted. See PDF.]

Fig 3. White matter regions with significant differences in LDH between patients and controls (P < 0.01, two-tailed, FWE corrected).

LDH, local diffusion homogeneity.

https://doi.org/10.1371/journal.pone.0162656.g003

Discussion

In this study, white matter impairments in schizophrenia patients were detected by both intra-voxel and inter-voxel metrics. We found that schizophrenia patients exhibited reduced FA and increased MD in the corpus callosum, cingulum, internal capsule, fornix and widespread superficial white matter in the frontal, parietal, occipital and temporal lobes. We also found decreased LDH in the corpus callosum, cingulum, internal capsule and fornix in schizophrenia. Our findings suggest that both intra-voxel and inter-voxel diffusion metrics are able to detect impairments in the anisotropic white matter regions, and intra-voxel diffusion metrics could detect additional impairments in the widespread isotropic white matter regions in schizophrenia.

According to a previous study, intra-voxel and inter-voxel metrics may yield discrepant findings even under the same pathological conditions [8]. FA represents a normalized ratio of diffusion directionality, whereas MD quantifies the bulk mobility of water molecules. Accordingly, FA reflects the degree of alignment of cellular structures within fiber tracts and their structural integrity, while MD reflects the overall situation of the whole molecule diffusion level and diffusion resistance. Limited by some conditions, both FA and MD reflect diffusion properties solely within the voxel. However, LDH, which measures the overall coherence of water molecule diffusion within a neighborhood, can reflect the micro-structural coherence of the underlying white matter fibers. Its alterations may represent disruptions of the local coherence of fiber orientations, myelination, diameter or density along the affected white matter tracts. These differences in the type of diffusion metrics may account for the discrepant findings among previous studies. However, more intriguingly, we did not observe any discrepant findings among the metrics of FA, MD and LDH in this study, which is inconsistent with previous findings in other neurological disorders such as Alzheimer’s disease [21] and stroke [8].

The lack of discrepant findings among these three indices in our study may be explained by several potential factors. First, the fact that schizophrenia patients show widespread white matter impairments was previously confirmed by numerous studies, and these white matter impairments included both inter-voxel and intra-voxel micro-structural aberrant white matter. Hence, in the anisotropic white matter regions such as the corpus callosum, cingulum, internal capsule and fornix, schizophrenia patients exhibited both altered inter-voxel and intra-voxel metrics. However, in the isotropic white matter regions such as white matter in the frontal, parietal, occipital and temporal lobes, the local diffusion homogeneity of the fiber tracts may not be affected or may be lightly affected; therefore, FA and MD, rather than LDH, can detect these alterations. Second, the sensitivity of these three indices is different for specific white matter micro-structural properties, which supports the inability of LDH to detect abnormality in fiber tracts in the frontal, parietal, occipital and temporal lobe regions. This study has several limitations that should be addressed in future work. First, we only focused on the capacity of LDH and FA/MD to detect abnormality in the white matter skeleton reconstructed using the TBSS method. These findings should be validated either in white matter regions using an appropriate atlas or in white matter fiber tracts reconstructed by tractography in future studies. Second, our study lacked correlation analyses between imaging parameters and clinical data. The primary purpose of this study was to test the feasibility of LDH in schizophrenia research. Thus, future studies that evaluate correlations of DTI parameters with clinical characteristics are needed.

In conclusion, our findings suggest that both intra-voxel and inter-voxel diffusion metrics are able to detect impairments in the anisotropic white matter regions, and intra-voxel diffusion metrics could detect additional impairments in the widespread isotropic white matter regions in schizophrenia. Thus, the combination of two types of diffusion metrics may provide complementary information for understanding the underlying pathological changes in schizophrenia.

Author Contributions

1. Conceptualization: CZ.

2. Formal analysis: XM HQ LW CZ FJ CW.

3. Methodology: XM HQ LW CZ.

4. Project administration: CZ.

5. Supervision: CZ.

6. Validation: CZ XM HQ LW.

7. Visualization: XM HQ.

8. Writing – original draft: CZ.

9. Writing – review & editing: CZ.

Citation: Zhuo C, Ma X, Qu H, Wang L, Jia F, Wang C (2016) Schizophrenia Patients Demonstrate Both Inter-Voxel Level and Intra-Voxel Level White Matter Alterations. PLoS ONE 11(9): e0162656. https://doi.org/10.1371/journal.pone.0162656

References

1. Catani M, Howard RJ, Pajevic S, Jones DK (2002) Virtual in vivo interactive dissection of white matter fasciculi in the human brain. Neuroimage 17: 77–94. pmid:12482069

2. Mori S, Kaufmann WE, Davatzikos C, Stieltjes B, Amodei L, Fredericksen K,et al. (2002) Imaging cortical association tracts in the human brain using diffusion-tensor-based axonal tracking. Magn Reson Med 47: 215–223. pmid:11810663

3. Pierpaoli C, Basser PJ (1996) Toward a quantitative assessment of diffusion anisotropy. Magn Reson Med 36: 893–906. pmid:8946355

4. Le Bihan D, Mangin JF, Poupon C, Clark CA, Pappata S, Molko N,et al. (2001) Diffusion tensor imaging: concepts and applications. J Magn Reson Imaging 13: 534–546. pmid:11276097

5. Basser PJ, Mattiello J, LeBihan D (1994) Estimation of the effective self-diffusion tensor from the NMR spin echo. J Magn Reson B 103: 247–254. pmid:8019776

6. Pierpaoli C, Barnett A, Pajevic S, Chen R, Penix LR, Virta A,et al. (2001) Water diffusion changes in Wallerian degeneration and their dependence on white matter architecture. Neuroimage 13: 1174–1185. pmid:11352623

7. Wedeen VJ, Rosene DL, Wang R, Dai G, Mortazavi F, Hagmann P, et al. (2012) The geometric structure of the brain fiber pathways. Science 335: 1628–1634. pmid:22461612

8. Bihel E, Roussel S, Toutain J, Bernaudin M, Touzani O (2011) Diffusion tensor MRI reveals chronic alterations in white matter despite the absence of a visible ischemic lesion on conventional MRI: a nonhuman primate study. Stroke 42: 1412–1419. pmid:21441158

9. Gong G (2013) Local diffusion homogeneity (LDH): an inter-voxel diffusion MRI metric for assessing inter-subject white matter variability. PLoS One 8: e66366. pmid:23776665

10. Abdul-Rahman MF, Qiu A, Sim K (2011) Regionally specific white matter disruptions of fornix and cingulum in schizophrenia. PLoS One 6: e18652. pmid:21533181

11. Ellison-Wright I, Bullmore E (2009) Meta-analysis of diffusion tensor imaging studies in schizophrenia. Schizophr Res 108: 3–10. pmid:19128945

12. Ellison-Wright I, Nathan PJ, Bullmore ET, Zaman R, Dudas RB, Agius M, et al. (2014) Distribution of tract deficits in schizophrenia. BMC Psychiatry 14: 99. pmid:24693962

13. Kunimatsu N, Aoki S, Kunimatsu A, Abe O, Yamada H, Masutani Y, et al. (2012) Tract-specific analysis of white matter integrity disruption in schizophrenia. Psychiatry Res 201: 136–143. pmid:22398298

14. Kyriakopoulos M, Bargiotas T, Barker GJ, Frangou S (2008) Diffusion tensor imaging in schizophrenia. Eur Psychiatry 23: 255–273. pmid:18524546

15. Nazeri A, Chakravarty MM, Felsky D, Lobaugh NJ, Rajji TK, Mulsant BH, et al. (2013) Alterations of superficial white matter in schizophrenia and relationship to cognitive performance. Neuropsychopharmacology 38: 1954–1962. pmid:23591167

16. Rosenberger G, Nestor PG, Oh JS, Levitt JJ, Kindleman G, Bouix S, et al. (2012) Anterior limb of the internal capsule in schizophrenia: a diffusion tensor tractography study. Brain Imaging Behav 6: 417–425. pmid:22415192

17. Smith SM, Jenkinson M, Johansen-Berg H, Rueckert D, Nichols TE, Mackay CE, et al. (2006) Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. Neuroimage 31: 1487–1505. pmid:16624579

18. Kay SR, Fiszbein A, Opler LA (1987) The positive and negative syndrome scale (PANSS) for schizophrenia. Schizophr Bull 13: 261–276. pmid:3616518

19. Cui Z, Zhong S, Xu P, He Y, Gong G (2013) PANDA: a pipeline toolbox for analyzing brain diffusion images. Front Hum Neurosci 7: 42. pmid:23439846

20. Smith SM, Nichols TE (2009) Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. Neuroimage 44: 83–98. pmid:18501637

21. Rose SE, Chen F, Chalk JB, Zelaya FO, Strugnell WE, Benson M, et al. (2000) Loss of connectivity in Alzheimer's disease: an evaluation of white matter tract integrity with colour coded MR diffusion tensor imaging. J Neurol Neurosurg Psychiatry 69: 528–530. pmid:10990518