High-field magnetic resonance imaging of structural alterations in first-episode, drug-naive patients with major depressive disorder.

Previous structural imaging studies have found evidence of brain morphometric changes in patients with major depressive disorder (MDD), but these studies rarely excluded compounding effects of certain important factors, such as medications and long duration of illnesses. Furthermore, the neurobiological mechanism of the macroscopic findings of structural alterations in MDD patients remains unclear. In this study, we utilized magnetization transfer imaging, a quantitative measure of the macromolecular structural integrity of brain tissue, to identify biophysical alterations, which are represented by a magnetization transfer ratio (MTR), in MDD patients. To ascertain whether MTR changes occur independent of volume loss, we also conduct voxel-based morphometry (VBM) analysis. The participants included 27 first-episode, drug-naive MDD patients and 28 healthy controls matched for age and gender. Whole-brain voxel-based analysis was used to compare MTR and gray matter volume across groups and to analyse correlations between MTR and age, symptom severity, and illness duration. The patients exhibited significantly lower MTR in the left superior parietal lobule and left middle occipital gyrus compared with healthy controls, which may be related to the attentional and cognitive dysfunction in MDD patients. The VBM analysis revealed significantly increased gray matter volume in right postcentral gyrus in MDD patients. These findings in first-episode, drug-naive MDD patients may reflect microstructural gray matter changes in the parietal and occipital cortices close to illness onset that existed before volume loss, and thus potentially provide important new insight into the early neurobiology of depression.

Introduction

As one of the leading worldwide causes of disability, major depressive disorder (MDD) is characterized by persistent, pervasive feelings of sadness, guilt and worthlessness, and often results in an increased risk of suicide.[1] MDD causes significant individual suffering and impairs social and occupational functioning, thus resulting in a major public health and economic burden. In vivo studies of the brains of MDD patients serve as an important strategy for learning about the neurobiology of depression and are helpful for developing effective therapies.

In recent decades, numerous structural magnetic resonance imaging (MRI) studies have sought to identify the key brain areas involved in the pathogenesis of depressive symptoms. Voxel-based morphometry (VBM) studies have reported that MDD is associated with widespread local abnormalities in many brain regions, such as the frontal gyrus,[2, 3, 4] insula,[4] temporal lobes,[4] anterior cingulate cortex,[3, 4, 5] striatum[6] and hippocampus.[7] Larger volumes of white matter in the inferior parietal lobule[8] and periventricular[9] regions and deep or subcortical[10] white matter hyperintensity have also been reported. However, these reports are inconsistent in detail, probably because of potential confounders associated with medications and/or long duration of illnesses. Thus, studies of first-episode, drug-naive MDD are a starting point for assessing brain structure before it is influenced by potential confounders and thereby provide core information relevant to models of pathogenesis.[11]

The neurobiological mechanism of the macroscopic structural alterations in MDD patients remains unclear. A post-mortem study has demonstrated reduced levels of the metabotropic glutamate receptor 5 monomer protein in the prefrontal cortex in cases of depression.[12] MDD patients have also showed reduced density of neuronal cell bodies with large cell body size in cortical layers 2 through 5 of the orbitofrontal cortex (OFC), and in layers 2, 3 and 6 of the dorsolateral prefrontal cortex (dlPFC), concurrent with an increased density of small body size neurons in layer 3 (OFC) and layers 3 and 6 (dlPFC).[13, 14] Observations of reduced glial cell density in the dlPFC and anterior cingulate cortex in the post-mortem brain of depressed patients provide further evidence of glial changes in patients with MDD.[15, 16] However, it is difficult to characterize the histopathological changes in vivo that underlie the macroscopic abnormalities. Techniques with sufficient neuropathological sensitivity to relatively subtle macromolecules are needed.

Magnetization transfer imaging (MTI) is a MR technique with the potential for providing more neuropathological information in vivo to subtle or early neuropathological changes than volumetric MRI.[17, 18] By utilizing the magnetization exchange between the spins of free water and water bound to macromolecules, MTI creates a contrast between tissues, namely, magnetization transfer ratio (MTR) map. The MTR measures the efficiency of these exchange phenomena, which depends on biophysical integrity of macromolecular protein pools and their local microenvironment.[19] MTI has been used to characterize the status of macromolecular proteins in several clinical brain disorders including depression,[20, 21, 22, 23] schizophrenia[24] and multiple sclerosis, even when a conventional MRI is negative.[25, 26] Previous MTI studies of depression have reported lower MTR in the corpus callosum, caudate nucleus, putamen, and the occipital white matter in late-life depressed patients,[21] and lower MTR in the head of the caudate nucleus in both MDD patients[22] and MDD patients with diabetes.[27] Treatment-resistant depressed patients have been reported to show reduced MTR in the anterior cingulate cortex, insula, caudate tail and amygdala-parahippocampal areas compared with healthy controls.[20]

However, these studies included patients taking antidepressants or patients who were free of antidepressants for at least 2 weeks. Furthermore, most of these studies used only region-of-interest (ROI) methods that focused on specific brain regions. In this study, we use both whole brain and ROI MTI approaches to explore and further characterize the neuropathological abnormalities in vivo in first-episode, drug-naive MDD patients relative to healthy controls without the confounding effects of antidepressant therapy. To ascertain whether MTR changes occurs independent of volume loss, we also conduct VBM analysis.

Materials and methods

Subjects

We recruited patients with 27 MDD from the Department of Psychiatry at West China Hospital of Sichuan University. All patients were assessed by an experienced psychiatrist and determined to meet the Diagnostic and Statistical Manual of Mental Disorders-IV (DSM-IV) criteria for MDD.[28] All patients also met the following inclusion criteria: (1) were between 18 and 60 years; (2) were experiencing first-episode depression; and (3) had a 17-item Hamilton Rating Scale for Depression (HRSD)[29] total score of 18 or higher on the day of their scan. Exclusion criteria included the following: (1) currently taking or had previously taken antidepressant or anti-psychotic drugs; (2) currently taking psychotropic medications, including benzodiazepines, hypnotics or anticonvulsant agents; (3) a history of electroconvulsive therapy; (4) contraindications to the MRI scan; (5) any psychotic disorders other than MDD; and (6) significant medical or neurological illness, including any history of significant head trauma with the loss of consciousness.

We recruited 28 healthy comparison subjects from the local area by poster advertisement who were assessed using the Structured Clinical Interview for the DSM-IV. Exclusion criteria included the same medical and psychiatric factors used to recruit MDD patients, as well as any DSM-IV Axis I disorder or known history of significant psychiatric illness. All participants were of Chinese Han nationality.

The West China Hospital Clinical Trials and Biomedical Ethics Committee of Sichuan University approved the study protocol (registration no. 84). The experiment was performed in accordance with the approved guidelines, and written informed consent was obtained from all participants.

MRI Acquisition

Images were acquired using a 3.0 Tesla General Electric magnetic resonance scanner (EXCITE, General Electric Medical Systems, Milwaukee, WI, USA). Participants were fitted with soft earplugs that were positioned comfortably in the coil, and they were then instructed to relax and remain still. Head motion was minimized using foam pads. Whole-brain MT images were acquired using a three-dimensional fast, low-angle shot sequence. One acquisition was performed with, and another without, the magnetization saturation pulse at 1.5 kHz off-resonance. In this way, MT- and non-MT-weighted images were generated separately. The other sequence parameters were as follows: repetition time/echo time (TR/TE)=37/5 ms; flip angle=15° 50 contiguous axial slices with slice thickness=3 mm; field of view=24 × 24 cm2; and data matrix=320 × 192. High-resolution, three-dimensional, T1-weighted images were acquired using a spoiled gradient recalled sequence (repetition time=8.5 ms, echo time=3.4 ms, fractional anisotropy=12°, 156 axial slices with thickness of 1 mm, axial field of view=240 × 240 mm, data matrix=256 × 256).

Image processing

MR images of all the subjects were first reviewed by a neuro-radiologist to ensure that there were no structural abnormalities or data quality flaws. Data processing and analyses were conducted using the statistical parametric mapping software SPM8 (Wellcome Trust Centre for Neuroimaging, London, UK; http://www.fil.ion.ucl.ac.uk/spm). For each subject, we first co-registered the MT- and non-MT-weighted images using a mutual information registration algorithm. MTR was calculated on a voxel-by-voxel basis as follows: MTR=(M0−Ms)/M0 × 100, where M0 and Ms are the signal intensities without and with, respectively, the saturation pulse applied. As the non-MT images are partially T1-weighted, we directly normalized them to the MNI T1 W template and then used the transformation parameters to normalize the co-registered MTR map. The normalized non-MT images were skull-stripped using the brain extraction tool (http://www.fmrib.ox.ac.uk/fsl/bet/)[30] and were then used as masks to remove non-brain tissues on the normalized MTR maps. Finally, MTR maps were smoothed with a Gaussian kernel of 6-mm full width at half maximum.

Voxel-based magnetization transfer analysis

MTR maps of MDD patients were compared with those of healthy controls using two-tailed two-sample t-tests controlling for age and gender in SPM8 (Wellcome Trust Centre for Neuroimaging). Significance in the resulting statistical maps was set at 0.05 and corrected for multiple comparisons. The correction standard was determined by Monte Carlo simulations, with the following parameters: individual voxel P-value=0.01, 1000 simulations, full width at half maximum=6 mm, applied using the Resting-State fMRI Data Analysis Toolkit of the AlphaSim program (http://www.restfmri.net).[31] In this way, a corrected significance level of P<0.05 was obtained for a minimum cluster size of 132 voxels. To quantify changes in the affected regions, MTR values were extracted using a volume-of-interest approach in SPM. We conducted correlation analyses between the average regional values in these regions and the patients' age, HRSD score and illness duration.

ROI analysis

We were also interested in determining the MTR changes of specific brain regions, such as the bilateral medial OFC, hippocampus, caudate nucleus, globus pallidus and thalamus, as they have been frequently reported in studies of depression.[32, 33, 34] These brain regions were selected as ROIs, and the segmentation was performed with the Freesurfer software package (http://surfer.nmr.mgh.harvard.edu/; Figure 1). In brief, we first conducted the atlas-based analysis using the Desikan–Killiany atlas[35] to segment the T1 image of each subject into 68 cortical regions (34 per hemisphere) and 40 subcortical structures. The T1 image of each subject was then co-registered to his/her non-MT-weighted image. The resultant transformation was used to map the segmented map to the non-MT image space. Finally, we overlaid the segmented map onto the MTR images and extracted the MTR values of bilateral ROIs using ITK-SNAP (www.itksnap.org),[36] and the MTR values of MDD patients were compared with those of healthy controls using two-tailed univariate analysis of covariance controlling for age and gender.

Figure 1

Regions of interest of (a) caudate nucleus, thalamus; (b) medial orbitofrontal cortex (OFC); (c) globus pallidus; (d) hippocampus for magnetization transfer ratio analysis.

VBM analysis

To identify whether morphological abnormalities could be detected in the first-episode, drug-naive MDD patients, T1W SPGR images were used to conduct VBM analysis using Exponentiated Lie algebra (DARTEL)[37] as implemented in SPM8. Two-tailed two-sample t-tests controlling for age and gender were performed to compare group differences of gray matter volume between first-episode, drug-naive MDD patients and healthy controls. Significance in the resulting statistical maps was set at 0.05 and corrected for multiple comparisons using the Resting-State fMRI Data Analysis Toolkit AlphaSim utility same as in voxel-based magnetization transfer analysis.

Results

Demographic and clinical comparisons

Table 1 presents the demographic and clinical characteristics of the study participants. There were no significant differences between the first-episode, drug-naive MDD patients and the healthy controls in gender, age or handedness (P>0.05).

Table 1

Demographics and clinical characteristics of subjects

Characteristics	MDD (n=27)	HC (n=28)	P-value
Age (years), mean (s.d.)	33 (10.8)	33 (11.7)	0.951*
Gender (male/female)	13/14	14/14	0.891**
Handness (right/left)	27/0	28/0	>0.99**
Illness duration (months), mean (s.d.)	79 (86)	NA	—
Total HRSD score, mean (s.d.)	22 (3.4)	NA	—

Abbreviations: HC, healthy controls; HRSD, Hamilton Rating Scale for Depression; MDD, major depressive disorder; NA, not applicable.

* and ** indicate the P-values for the two-sample t-test and χ2-test, respectively.

As presented in Table 2, the first-episode, drug-naive MDD group exhibited lower MTR in the left superior parietal lobule (SPL; Figure 2a) and left middle occipital gyrus (MOG; Figures 2b and c) compared with the healthy controls (P<0.05, corrected for multiple comparisons). No marked regional MTR increase in the MDD group compared with the controls was found. Figure 2d illustrates the changes in the MTR values in the affected regions (t =4.244, P<0.001 for left SPL; t =3.681, P=0.001 and t =3.014, P=0.004 for the two left MOG regions, respectively). No association was found between regional MTR values and patient age, HRSD scores or illness durations (P=0.391 for age and left SPL; P=0.05 and P=0.985 for age and two left MOG regions; P=0.073 for HRSD scores and left SPL; P=0.424 and P=0.507 for HRSD scores and two left MOG regions, respectively; P=0.143 for illness duration and left SPL; P=0.319 and P=0.8 for illness duration and two left MOG regions, respectively).

Table 2

Regions showing differences of MTR and gray matter volume between MDD patients and healthy controls

Regions	Coordinates of peak voxel in MNI space (x, y,z)	Cluster size (voxels)	t scores of peak voxela
MTR analysis
HC>MDD
L superior parietal lobule	−32, −74, 56	401	3.87
L middle occipital gyrus	−30, −100, 2	214	3.52
L middle occipital gyrus	−52, −76, 2	420	3.47
MDD>HC
None
VBM analysis
MDD>HC
R postcentral gyrus	45, −33, 56	825	3.45
HC>MDD
None

Abbreviations: HC, healthy controls; L, left; MDD, major depressive disorder; MNI, Montreal Neurological Institute Coordinate System or Template; MTR, magnetization transfer ratio; R, right; T, statistical value of peak voxel showing differences of MTR and gray matter volume between the two groups; VBM, voxel-based morphometry; x, y, z, coordinates of primary peak locations in the MNI space.

All effects survived a voxel-wise statistical threshold (P<0.05), as corrected for multiple comparisons using AlphaSim program.

Figure 2

Magnetization transfer ratio differences in voxel-based analysis comparisons between patients with major depressive disorder and healthy controls. Images are presented in neurological orientations. Patients with major depressive disorder exhibited a reduced magnetization transfer ratio in (a) the left superior parietal lobule and (b and c) the left middle occipital gyrus relative to healthy controls. Statistical inferences were made with a voxel-level statistical threshold of P<0.05 (corrected). (d) Quantification of MTR values in the affected regions. HC, healthy control; L, left; MDD, major depressive disorder; MOG, middle occipital gyrus; MTR, magnetization transfer ratio; SPL, superior parietal lobule.

A ROI analysis revealed no significant group differences between the MDD patients and healthy controls in the MTR in bilateral medial OFC (left t=1.801, P=0.186; right t=0.238, P=0.628), hippocampus (left t=0.582, P=0.434; right t=0.032, P=0.437), caudate nucleus (left t=0.485, P=0.489; right t=0.178, P=0.675), globus pallidus (left t=1.198, P=0.465; right t=1.227, P=0.24) or thalamus (left t=0.114, P=0.737; right t=0.916, P=0.343) when controlling for age and gender (Supplementary Table 1).

The VBM analysis revealed increased gray matter volume in the right postcentral gyrus in first-episode, drug-naive MDD group compared with the healthy controls (P<0.05, corrected for multiple comparisons; Table 2; Figure 3). No marked gray matter decrease in the MDD group compared with the controls was found.

Figure 3

Gray matter volume differences in voxel-based analysis comparisons between patients with major depressive disorder and healthy controls. Images are presented in neurological orientations. Patients with major depressive disorder exhibited increased gray matter volume in the right postcentral gyrus compared with healthy controls. Statistical inferences were made with a voxel-level statistical threshold of P<0.05 (corrected). R, right.

Discussion

In the present study, we used MTI to explore the microstructural changes in first-episode, drug-naive MDD patients relative to healthy controls. The main findings are that the MDD patients exhibited lower MTR in the left superior parietal lobule and left middle occipital gyrus relative to healthy controls. To our knowledge, this is the first MTI study of first-episode, drug-naive MDD patients using both voxel-based and ROI analyses.

MTI techniques provide semi-quantitative (MTR) metrics to assess the integrity of macromolecules in the brain tissue such as myelin or cellular membranes components (proteins, cholesterol, phosphatidylcholine and galactocerebrosides).[38] The neuropathological correlates of the MTR changes involving gray and white matter are under contention and may be related to multiple neurobiological alterations.[18, 39, 40, 41] Reductions of MTR in gray matter are indicative of decreased size and number of neurons, or dendritic density, and may reflect abnormal cell membrane structure.[39] In the white matter, MTR has been shown to reflect mainly demyelination, some neuro-axonal loss and, to a lesser extent, microglial activation.[42]

In this study, we found the first-episode, drug-naive MDD patients had a significantly lower MTR in the left SPL than did the healthy controls. The SPL is an important brain region of the dorsal attention network, which is hypothesized to have a key role in the top-down, or goal-driven, allocation of attention.[43] This brain region has been consistently associated with the volitional allocation of attention to visual aspects of, and coordinated motor actions to, specific regions in extrapersonal space.[44, 45] Previous functional MRI studies have revealed that shifting the location of attention produces a rapid, transient increase in SPL activation.[46, 47] Impairments in attention may be due to the neuropsychological deficits that correspond to the cognitive criterion ‘impaired ability to think or concentrate' as stipulated for MDD in the DSM-IV.[28] A structural neuroimaging study has reported a reduction in the gray matter volume in the SPL in MDD patients compared with normal subjects.[48] Functional MRI studies have reported that first-episode, drug-naive MDD patients exhibited reduced resting-state activity in the parietal lobe,[49] increased activation in bilateral SPL after escitalopram oxalate treatment during sad facial expression recognition[50] and decreased functional connectivity between the dlPFC and parietal lobe in first-episode MDD patients.[51] Furthermore, decreased glucose metabolism was also observed in the parietal cortex of depressed patients compared with that of the healthy controls.[52] Therefore, our findings are generally consistent with these studies and further reveal that the microstructural gray matter changes in the parietal lobe may have already existed in the early stage of MDD because our group only included first-episode, drug-naive MDD patients. Together with previous findings on the parietal lobe in MDD patients, our result of lower MTR in the superior parietal lobe of MDD patients may improve our understanding of the cerebral mechanisms underlying attentional deficits in depression.

We also identified a lower MTR in the left middle occipital gyrus in first-episode, drug-naive MDD patients compared with that of healthy controls. The middle occipital gyrus is within the visual recognition network[53] and is involved in the perception of facial emotion.[54] Depressed patients are thought to have a negative cognitive bias in their appraisal of people or life events and are believed to be unable to carry out normal interpersonal interactions.[55] Thus, abnormal facial emotion processing may give rise to some of the affective and social symptoms in major depression. Neuroimaging studies have suggested that the occipital cortex may be relevant in MDD. A study comparing patients with psychotic major depression, patients with non-psychotic major depression and a healthy group of individuals revealed greater right occipital activation only in the non-psychotic major depression group compared with healthy controls during a verbal working memory task, thus suggesting that increased effort is required for basic visual processing.[56] Decreased resting-state activity in the left middle occipital gyrus has been reported in MDD patients compared with healthy controls.[57, 58] It has also been reported that regional cerebral blood flow (rCBF) in the parieto-occipital regions was decreased in MDD patients and was significantly increased following pharmacotherapy.[59] A magnetic resonance spectroscopy study found that gamma-aminobutyric acid (GABA) neurotransmission was altered in the occipital cortex of patients with MDD.[60] The molecular mechanism behind the neuroimaging results may be related to the reduced number of GABAergic neurons in the occipital cortex.[61] Consistent with our results, these findings of microstructural gray matter changes alterations in the middle occipital gyrus indicate that the visual recognition network interference may contribute to cognitive dysregulation in MDD.

Whole-brain voxel-based analysis investigates region-specific changes by spatially normalizing brain images and conducting statistical tests at each voxel, which is both reproducible and automatic. Whereas ROI analysis, usually applied in the original image space, allows one to carry out group analyses of functionally corresponding regions with high statistical power and substantially reduces the correction needed for multiple comparisons as is needed for the many thousands of voxels of a whole-brain analysis.[62] In this study, bilateral medial OFC, hippocampus, caudate nucleus, globus pallidus and thalamus were chosen as ROIs based on previous studies and no significant difference of MTR was found between MDD patients and healthy control group that was consistent with the our whole-brain voxel-based analysis. Nonetheless, lower MTR has been demonstrated in gray matter regions in previous MTI studies of MDD.[20, 22, 63] Two MTI studies compared the MTR value changes between MDD patients and healthy controls by manually placing ROIs in the key components of cortical–subcortical circuits—the caudate, thalamic, striatal, orbitofrontal, anterior cingulate and dorsolateral prefrontal regions, and the results showed significant group differences of MTR values in the bilateral caudate[22] and thalamus.[63] However, it should be noted that the MDD patients in these two studies were free of antidepressant medication at least 2 weeks, whereas patients in our study were first-episode and drug-naive. In another whole-brain MTR study, treatment-resistant depressed patients with long mean illness duration were recruited, and the results showed a reduced MTR in the anterior cingulate cortex, insula, caudate tail and amygdala-parahippocampal areas relative to healthy controls.[20] The discrepancy between our results and those of previous studies may be mainly attributed to inter-study differences in sample size, the sociodemographic and clinical characteristics of the participants, such as the various medication statuses and the incompatible illness duration.

The VBM analysis revealed increased gray matter volume in the right postcentral gyrus in the first-episode, drug-naive MDD group compared with the healthy controls. The postcentral gyrus is the primary sensory area of the brain. The somatosensory-related cortices are considered to not only encode bodily sensations but also have a key role in using social cues to understand emotional states of others.[64] A study has reported that the MDD patients showed reduced cerebral activation in the right postcentral gyrus than healthy controls during empathy for others pain.[65] Similarly, another functional MRI study demonstrated that following escitalopram oxalate treatment, first-episode, treatment-naive MDD patients showed decreased brain activation in bilateral postcentral gyrus during emotion recognition.[50] Although the underlying mechanism for increases in gray matter volume in our study is still unknown, one possible account for this effect may be the inflammatory response, as a compensatory effect in the early phase of depression.[66, 67] It is worth mentioning that different regions were implicated when MTI and VBM methods were separately used. MTI provides information on the integrity of macromolecules in the brain tissue,[19] whereas VBM allows for the study of the gray matter morphology. These isolated findings of MTI and VBM analyses suggested that MTR changes in occipital and parietal cortices existed beyond volume loss in the early stage of MDD which needs further longitudinal studies.

Our study has certain limitations. First, the sample size is modest. It is possible that a larger sample size will provide further insight. However, it should be noted that our patient group was a homogeneous group of first-episode MDD patients without medication. Second, MTR is very sensitive to pulse sequences and relaxation properties. This makes it difficult to accurately address the specific neurobiological causes of the identified MTR alterations. Third, as our study uses a cross-sectional design, the question arises whether these differences are a consequence or precondition of MDD. This issue can only be resolved by longitudinal design studies.

In conclusion, this study demonstrates that first-episode, drug-naive MDD patients exhibit structural abnormalities, namely, lower MTR in the left superior parietal lobule and left middle occipital gyrus and increased gray matter volume in right postcentral gyrus that largely avoided the confounding effects from medications and/or long duration of illness. The MTR changes in the parietal and occipital cortices existed beyond volume loss in the early stage of MDD. The results indicate that the attentional and cognitive dysfunction in MDD patients may be associated with microstructural gray matter changes in the parietal and occipital cortices close to the illness onset that existed before volume loss, and thus potentially provide important new insight into the early neurobiology of the depression.