Microstructural Abnormalities in Subcortical

Full text

Turn on search term navigation

About the Authors:

Anne J. Blood

Affiliations Depression Clinic and Research Program, Mood and Motor Control Laboratory, Addiction Research Program, Laboratory of Neuroimaging and Genetics, Department of Psychiatry, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America, Motivation and Emotion Neuroscience Collaboration (MENC) and Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America, Department of Neurology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

Dan V. Iosifescu

Nikos Makris

Roy H. Perlis

David N. Kennedy

Affiliations Center for Morphometric Analysis and Center for Integrative Informatics, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America, Department of Neurology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

Darin D. Dougherty

Affiliation: Depression Clinic and Research Program, Mood and Motor Control Laboratory, Addiction Research Program, Laboratory of Neuroimaging and Genetics, Department of Psychiatry, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

Byoung Woo Kim

Myung Joo Lee

Shirley Wu

Sang Lee

Jesse Calhoun

Steven M. Hodge

Affiliations Depression Clinic and Research Program, Mood and Motor Control Laboratory, Addiction Research Program, Laboratory of Neuroimaging and Genetics, Department of Psychiatry, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America, Motivation and Emotion Neuroscience Collaboration (MENC) and Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America, Center for Morphometric Analysis and Center for Integrative Informatics, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

Maurizio Fava

Bruce R. Rosen

Affiliation: Motivation and Emotion Neuroscience Collaboration (MENC) and Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

Jordan W. Smoller

Gregory P. Gasic

Hans C. Breiter

* E-mail: [email protected]

for the Phenotype Genotype Project on Addiction and Mood Disorders

Introduction

It has been proposed that major depressive disorder (MDD) may result from dysfunction of brain reward/aversion circuitry[1], [2], [3], [4], [5], [6], [7]. Hypothesized in 1999 as a general schema for processing both positive and negative features of potential “goal-objects” or states [8], a generalized reward/aversion system that processes the salience of stimuli across a continuum of aversion and reward was described in multiple publications between 1996 and 2001 [9], [10], [11], [12], [13]. This work identified an extended set of brain regions as variably processing a continuum between positive and negative valence and intensity information, along with category and incidence information from goal-objects. These observations have been extensively replicated [14], [15], [16], [17], [18] and synopsized [19], [20], [21]. The recent report of a law-like relationship between patterns of approach and avoidance behavior to rewarding and aversive stimuli further argues that the systems processing this information do not function independently of one another [22].

The possibility that abnormalities in this reward/aversion circuitry underlie many psychiatric conditions, including MDD, was further schematized by multiple investigators [1], [6], [23], [24]. Animal models of MDD strongly support this thesis and have hypothesized dopaminergic midbrain nuclei and the medial forebrain bundle (MFB) to be involved with the illness [7], [25]; these subcortical regions are at the core of the animal literature that first identified reward circuitry [26], [27], [28], [29], [30]. In humans with MDD, there has recently been detection of functional abnormalities in these subcortical regions [2], [31], [32], in addition to well-established abnormalities in target regions receiving subcortical projections [33], [34], [35], [36], [37], [38], [39], [40], [41], [42], [43], [44], [45], [46], [47]. However, there has not yet been an evaluation of microstructural integrity in the midbrain nuclei or MFB.

The current study aimed to evaluate microstructural integrity of subcortical brain reward/aversion circuitry in cohorts of subjects with and without MDD, using approaches that maximized detection sensitivity in subcortical regions. Our primary hypotheses focused on evaluating brain microstructure in the midbrain ventral tegmental area/substantia nigra (VTA/SN), the medial forebrain bundle (MFB), and the amygdalofugal pathway (AFP). Our secondary a priori hypotheses were based on a general evaluation of the broad set of reward/aversion circuitry implicated in MDD [8], [48].

We designed the study to take into consideration potential heterogeneity within the MDD cohort [33], [37], [49], [50], [51], [52], [53], [54], [55], [56], [57], given that this is a potential confound in group designs [58]. Such an approach may not only improve methods, but may also lead to identification of subgroups of MDD. We hypothesized there would be microstructural heterogeneity in our cohort that could not be predicted a priori by functional imaging abnormalities or symptom profiles since MDD has been hypothesized to be a systems-level disorder, and it is possible that microstructural abnormalities at different points in a distributed circuit could all lead to similar abnormalities in behavior/mood. To address this, we allowed the DTI data itself to drive potential segregation across patients.

Methods

Subjects

All subjects signed written informed consent prior to participation in this study, and the study was approved by the Institutional Review Board of Massachusetts General Hospital (Partners Human Research Committee). All experiments were conducted in accordance with the principles of the Declaration of Helsinki.

22 of 44 subjects met DSM-IV criteria for Major Depressive Disorder (MDD) diagnosed by physician-administered Structured Clinical Interview for DSM-IV Axis I Disorders - Patient Edition [SCID-P [59]] and were between the ages of 18 and 65 (mean age = 36.3±12.1 years; 12/22 females, mean educational history of 15.6±2.6 years, 19 Caucasian and 3 African American, 20/22 right-handed). These patients were matched one-to-one with 22 healthy volunteers on the following criteria: age (within 5 years), years of education (within 5 years), gender, self-reported race, and handedness (control mean age = 35.3±11.6 years; mean educational history of 15.7±2.1 years; 12/22 females, 19 Caucasian and 3 African American, 20/22 right-handed. Age and years of education did not differ significantly across groups (age: F = 0.001, p = 0.979;.education: F = 0.085, p = 0.772). This tight matching was done because each of these factors may potentially influence neural structure and function. All subjects were drawn from a larger study evaluating cocaine addicted, depressed and control subjects [The Massachusetts General Hospital (MGH) Phenotype Genotype Project in Addiction and Depression (PGP; http://pgp.mgh.harvard.edu)]. The 22 MDD subjects were drawn from a larger cohort of individuals with MDD (n = 47) because they met quality assurance criteria, including (a) minimal residual motion artifacts after motion correction of DTI images, (b) absence of MR susceptibility artifacts, and (c) availability of control subjects which met the strict matching criteria described above, and who also had DTI data without motion or susceptibility artifacts.

MDD subjects were excluded if they met DSM-IV criteria for primary psychotic disorders, bipolar disorder, eating disorders, substance abuse disorders, generalized anxiety disorder, panic disorder, PTSD or OCD by SCID interview (current or lifetime); healthy volunteers were excluded if they met DSM-IV criteria for any Axis I psychiatric disorder by SCID interview.

Additional exclusion criteria for both MDD subjects and healthy volunteers were: 1) currently suicidal or at risk for suicide in the judgment of the investigator; 2) pregnant women; 3) carrying a medical device incompatible with MRI (e.g., metal implants such as surgical clips or pacemakers) or significant claustrophobia or weight that would make MRI unfeasible; 4) serious medical illness including a known history of HIV-1+ status; 5) Subjects with insulin dependent diabetes mellitus (IDDM) or subjects with noninsulin dependent diabetes mellitus (NIDDM) and abnormal Hemoglobin A1C; 6) severe respiratory compromise; 7) history of head trauma with neurological sequelae, including multiple concussions and/or history of stroke; 8) history of seizure disorder, delirium, dementia, or mental disorders due to general medical conditions; 9) clinical or laboratory evidence of uncontrolled hypothyroidism or hyperthyroidism; and 10) subjects which, in the opinion of the Management Group running the PGP, were not able to participate safely in this study. In addition, subjects in the larger PGP study were screened for Hepatitis C (by Hepatitis C+ titer); no subjects included in this study tested positive for Hepatitis C. Three MDD subjects and one control subject had a history of tobacco use (one current smoker in MDD cohort, one current smoker in the control cohort, and two previous smokers in the MDD cohort).

In addition to screening for exclusion criteria, we also performed for all subjects (MDD and healthy controls): (1) the Structured Clinical Interview for DSM-IV for Axis I disorders (SCID-I); (2) a medical history and concurrent medication status (see Table 1 for medication status in relation to MDD subgrouping and VTA/SN FA values); (3) sociodemographic information; (4) Edinburgh Handedness assay; (5) blood and urinary analysis; (6) Inventory of Depressive Symptomatology – Self Report (IDS-SR); (7) 31-item Hamilton Rating Scale-Depression (HAM-D, [60]); and (8) the STAI (Spielberger State Trait Anxiety Inventory).

[Figure omitted. See PDF.]

Table 1. Antidepressant History and Status for MDD Cohort.

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

DTI Scanning Protocol

During each scanning session, a high-resolution (2 mm isotropic) whole head DTI scan was acquired on a Siemens 3.0 Tesla Sonata Magnet System (Siemens AG, Medical Solutions, Erlangen, Germany); TR = 24 s; TE = 81 ms; slice thickness = 2 mm, 60 slices total, 128×128 matrix, 256×256 mm FOV, 6 averages, 6 noncolinear directions with b value = 700 s/mm2, and 1 image with b-value = 0 s/mm2. DTI scans were acquired using auto-align software [61] to normalize brain image slice orientation across subjects. Slices were situated parallel to the AC-PC line, and parallel to the inside curve of the FOC to minimize signal distortion in this region [62].

DTI Image Preprocessing and Registration

1. DTI preprocessing.

All data processing was performed using Freesurfer software (http://surfer.nmr.mgh.harvard.edu) and FSL (http://www.fmrib.ox.ac.uk/fsl) processing streams. Detailed methods for preprocessing DTI data are described in Salat et al. [63]. Briefly, each tensor volume from the DTI dataset was resampled to the T2 image to correct for remaining eddy current distortion, and to correct for head motion. The fractional anisotropy (FA) metric was derived from the diffusion tensor as previously described [64]. All resulting maps were resampled to 1 mm3 resolution. Images were inspected for residual motion and susceptibility artifacts; subjects with significant artifacts were excluded.

2. Image registration.

FA maps were registered to a Montreal Neurological Institute (MNI) ICBM152 T2 template for the voxel-based contrast and subsequent analyses; all contrast analysis results are reported in MNI coordinates. Registration procedures (both automated and manual) were done blinded to subject and group identification. Given the small size of our areas of interest, our methods included registration techniques designed to maximize intersubject alignment in our areas of interest. Registration and analysis of the VTA/SN, in particular, has been previously shown to be valid and informative, given abnormalities have been detected in the SN in patients with Parkinson's Disease, a population known to exhibit structural pathology in this region [65], [66]. In PD patients, FA was decreased in the SN, which is consistent with neuronal loss in this region [65], [66]. While the Parkinson's studies did not use our directed registration methods, the abnormalities in these patients (i.e. the “signal”) was presumed to be much greater in Parkinson's patients, making it less susceptible to being masked by imprecision in registration of small nuclei by automated methods (i.e. “noise”). Thus, in a disorder where we expect more subtle abnormalities in this region (e.g., MDD), we believe the directed registration methods were required to resolve the signal above the noise.

a. Automated registration. The initial image registration was done using FSL software to perform an automated, 12-degrees of freedom, global affine transformation[67]. The linear affine transformation (in combination with a subsequent landmark-guided manual registration) was selected over non-linear transform methods to minimize loss of signal due to warping, which is relatively greater in the small subcortical regions [68] relevant to our a priori hypotheses.

b. Manual registration. After the automated registration, all registered images were put through an additional manual registration step using Martinos Center Freesurfer software (http://surfer.nmr.mgh.harvard.edu). This step maximized the accuracy of registration of each subject's FA map to ICBM152 T2 space, based on three registration landmarks selected to maximize registration in our a priori areas of evaluation (AOEs) (see Text S1). Figures S1, S2, S3, S4 illustrate anatomical co-localization across MDD and control cohorts overall, and within our primary AOEs.

Data Analysis

1. A priori hypotheses and segmentation of a priori areas of evaluation (AOEs).

Primary and secondary a priori hypotheses corresponded to a priori AOEs for constraining our voxel-based search of patient/control differences. These AOEs constrained which voxels were evaluated in FA map group contrasts, to identify clusters meeting volume and significance thresholds. AOEs were segmented by an anatomist (N.M.) using landmark-based, atlas-guided definitions of the regions (see Text S1 and Figures S1, S2, S3, S4, S5, S6, S7).

a. Primary a priori hypothesis. We predicted MDD subjects would exhibit abnormal brain microstructure in (1) the medial forebrain bundle (MFB) and contiguous lateral nucleus of the hypothalamus (LNH), along with (2) regions feeding the MFB/LNH, including the ventral tegmental area/substantia nigra (VTA/SN), and (3) the amygdalofugal pathway. Follow-up analyses (described below), considered clusters falling within these regions.

b. Secondary a priori hypotheses. Based on imaging findings in the MDD literature, we further hypothesized that MDD subjects would exhibit abnormal brain microstructure in white matter adjacent to orbitofrontal cortex (FOC), anterior cingulate cortex (ACC), paracingulate cortex (PAC), subgenual prefrontal cortex (SGC), and dorsolateral prefrontal cortex (DLPFC).

c. Brain regions outside a priori hypothesized AOEs. We evaluated all other brain regions outside primary and secondary a priori AOEs using a whole-brain Bonferroni correction.

2. Voxel-based contrast of FA maps for MDD subjects versus control subjects.

a. Contrast analysis. A voxel-based contrast analysis was performed between MDD and control subjects, using a two-tailed t-test with Freediffusion software [63], [69]. To minimize the chance of false positives we required (1) that group differences meet a cluster threshold and (2) that the p value of the peak voxel within each cluster meet a correction for the number of voxels in a search volume.

b. Cluster thresholds for group contrast of FA maps. For primary a priori AOEs, the cluster threshold was at least 9 contiguous voxels, with p<0.05 for each voxel. This was increased to 27 contiguous voxels for secondary a priori AOEs, and 81 contiguous voxels for regions outside a priori AOEs. The least stringent threshold exceeded cluster thresholds used with fMRI [30], [70] and cortical thickness analysis[71]. We defined contiguous voxels as voxels sharing an edge (i.e., not just a corner). To be considered within that AOE, greater than 50% of voxels in a cluster had to fall within a segmented a priori AOE.

c. Multiple comparisons correction for group contrast of FA maps. To be considered statistically significant, clusters were required to have a peak voxel meeting a Bonferroni correction for multiple comparisons. Corrections were based on the total number of voxels in the area being considered (418 for primary a priori AOEs; 9,308 for secondary a priori AOEs; 200,000 for the whole brain [i.e., not a priori]), divided by the required cluster size. Thus, the uncorrected p value (reported in Tables 2, 3, 4) was required to be p<0.05/(418/9), or p<0.00108, for primary a priori AOEs, p<0.05/(9,308/27), or p<1.45×10−4, for secondary a priori AOEs, and p<0.05/(200,000/81), or 2.03×10−5, for the rest of the brain. Peak voxels within an order of magnitude of the corrected p value were considered trends toward significance. All regions that met cluster threshold criteria were tabulated; however, the results and discussion sections focus primarily on regions which met full significance criteria.

[Figure omitted. See PDF.]

Table 2. Voxel-Based Contrast for the Whole Cohort (22 MDD Versus 22 Matched Controls).

https://doi.org/10.1371/journal.pone.0013945.t002

[Figure omitted. See PDF.]

Table 3. Post-Hoc Contrast for “Abnormal VTA/SN” MDD Subgroup (10 MDD Versus 10 Matched Controls).

https://doi.org/10.1371/journal.pone.0013945.t003

[Figure omitted. See PDF.]

Table 4. Post-Hoc Contrast for “Normal VTA/SN” MDD Subgroup (12 MDD Versus 12 Matched Controls).

https://doi.org/10.1371/journal.pone.0013945.t004

d. Permutation test as validation of our findings. As a supplemental test to further validate the robustness of findings, we computed a permutation test for regions meeting the cluster threshold. Methods and results for this analysis are reported in Tables S1, S2, S3, S4, S5, S6 in Dataset S1.

e. Laterality test. Given initial findings in the VTA/SN region were lateralized to the right hemisphere in MDD subjects, we assessed laterality of FA in this region for both groups to distinguish between potential presence of lateralized effects in patients, versus potential loss of laterality relative to controls (see Dataset S2).

3. Follow-up analyses.

Follow-up analyses were conducted to assess MDD group heterogeneity, and the relationship between FA and symptom profiles for MDD subjects.

a. Individual subject FA values. Once clusters representing group differences were identified in primary a priori AOEs, we measured mean FA across each cluster for each subject (see Figures S5, S6 for illustrations of these clusters). Mean FA values for individuals were then used in follow-up analyses. To evaluate whether MDD cohort heterogeneity led to false negatives, we included findings which met the cluster, but not the Bonferroni threshold in the initial voxel-based contrast.

b. Gapping analysis for individual FA values in MDD subjects. We used a gapping analysis [72] to test the likelihood that there were gaps in the distribution of FA values in the MDD cohort, suggesting the presence of two or more population distributions of these values. This analysis involved rank ordering individual VTA/SN FA values and calculating the mean gap distance between each of the middle 50% of values (eliminating the top 25% and bottom 25% of values to exclude potential outliers). The ratio of each individual gap distance between adjacent data points was calculated relative to mean gap distance. Finally, we identified the greatest gap ratio within the middle 50% of values and calculated the likelihood that this gap ratio would be observed by chance if data points reflected the distribution of a single population. This likelihood was calculated both for a Gaussian distribution and for a t distribution with df = 4. Each statistic was calculated using 10,000 simulations.

c. Voxel-based contrasts of FA maps for the two MDD subgroups. A voxel-based contrast was calculated for each of the two subgroups of MDD subjects identified in the gapping analysis versus their individually matched controls. Clusters were identified in these voxel-based contrasts and peak voxels corrected for multiple comparisons using the same criteria as in the initial group contrast.

d. VTA/SN FA and symptom measures. We evaluated whether symptom profiles differed across the two MDD subgroups. Specifically, we looked at (1) overall depression severity (total score on IDS-SR), and scores on individual IDS-SR items relating to (2) anhedonia (question 21), (3) sadness (question 5), and (4) psychomotor symptoms (question 30). Although MDD subjects in our cohort had not been clinically subtyped, anxious depression is one of the most common clinical subtypes of MDD [57], [73]; therefore we also evaluated measures of (5) trait anxiety (via the STAI-T). We used two-tailed t-tests to evaluate whether these measures differed between MDD subgroups. Significance was determined using a Bonferroni correction of p<0.05/5 = 0.01 for these analyses.

Results

1. Image Contrast Analysis of MDD Versus Control Subjects

a. Primary a priori areas of evaluation (AOEs) (Table 2; Figure 1).

Significantly elevated FA was detected in MDD subjects within the VTA/SN. Specifically, FA was elevated at the ventral/lateral edge of the right SN adjacent to the cerebral peduncle and overlying the nigral fiber system (striatonigral, nigrostriatal, and corticostriatal fibers; Figure 1). A trend toward reduced FA in MDD subjects was also noted in the left MFB where it passes through the LNH. There were no group differences for the amygdalofugal AOE. Mean FA values for each cluster detected in the voxel-based contrast are reported in Table S2.

[Figure omitted. See PDF.]

Figure 1. FA difference in MDD subjects in primary a priori areas of evaluation.

This figure shows the significant FA difference in MDD subjects versus matched healthy controls for the voxel-based contrast (all subjects: 22 MDD, 22 controls) within areas of evaluation (AOEs) included in our primary a priori hypotheses. MDD subjects exhibited elevated FA in the ventral tegmental area/substantia nigra (VTA/SN), which localized to the ventral/lateral edge of the substantia nigra (SN) adjacent to the cerebral peduncle, and the nigral fiber system (striatonigral, nigrostriatal, and corticostriatal fibers). Translucent green in image on left indicates the a priori AOE for the VTA/SN in the slice shown, which was used to constrain the initial contrast analysis. The color bar indicates the range of p values in this figure, from the threshold (p<0.05) to the order of magnitude of the voxel of peak significance (i.e. smallest p value) in the VTA/SN; color in images is viewed using trilinear interpolation. Warm tones (red, orange, yellow) indicate regions in which MDD subjects exhibited elevated FA relative to control subjects. RH: right hemisphere; LH: left hemisphere.

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

b. Secondary a priori AOEs (Table 2; Figure 2).

There were trends toward significantly reduced FA in two white matter regions underlying DLPFC (Figure 2), and in white matter adjacent to the right ACC.

[Figure omitted. See PDF.]

Figure 2. FA differences in MDD subjects in secondary a priori areas of evaluation.

This figure shows the significant FA differences in MDD subjects versus matched healthy controls for the voxel-based contrast (all subjects: 22 MDD, 22 controls) in a priori regions of evaluation (AOEs) included in our secondary a priori hypotheses. MDD subjects exhibited reduced FA in white matter regions underlying dorsolateral prefrontal cortex (DLPFCwm), bilaterally. Translucent green in images on left indicates the a priori AOE for DLPFCwm in the slices shown, which were used to constrain the initial contrast analysis. The color bar indicates the range of p values in this figure, from the threshold (p<0.05) to the order of magnitude of the voxel of peak significance (i.e. smallest p value) in DLPFCwm; color in images is viewed using trilinear interpolation. Cool tones (blues) indicate regions in which MDD subjects exhibited reduced FA relative to control subjects. RH: right hemisphere; LH: left hemisphere.

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

c. Other regions meeting the cluster threshold (Table 2; Figure 3).

Outside a priori AOEs, MDD subjects exhibited significantly elevated FA in white matter adjacent to right calcarine cortex (Figure 3B), and a trend toward significantly reduced FA in white matter within the left precentral gyrus, below premotor cortex (BA6) (Figure 3A).

[Figure omitted. See PDF.]

Figure 3. FA differences in MDD subjects in regions not in our a priori areas of evaluation.

This figure shows the significant FA differences in MDD subjects versus matched healthy controls for the voxel-based contrast (all subjects: 22 MDD, 22 controls) in regions not included in our a priori hypotheses. MDD subjects exhibited (A) reduced FA in white matter adjacent to right premotor cortex (PMCwm), and (B) elevated FA in white matter adjacent to right calcarine cortex (CCtx). The color bar indicates the range of p values in this figure, from the threshold (p<0.05) to the order of magnitude of the voxel of peak significance (i.e. smallest p value) in PMCwm for reduced FA and from the threshold (p<0.05) to the order of magnitude of the voxel of peak significance in CCtx white matter for elevated FA, color in images is viewed using trilinear interpolation. Warm tones (red, orange, yellow) indicate regions in which MDD subjects exhibited elevated FA relative to control subjects; cool tones (blues) indicate regions in which MDD subjects exhibited reduced FA relative to control subjects. RH: right hemisphere; LH: left hemisphere.

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

2. Follow-up Analyses: Gapping, Subgroup Subtractions, Symptom Correlations

a. Gapping analysis.

A gap in the distribution of mean FA values across MDD subjects was observed for the VTA/SN cluster (Figure 4A). This gap was located at the upper limit of control values; all MDD subjects with values above the gap were outside the range of control values. Gapping analysis indicated it was unlikely that this gap in the MDD VTA/SN FA values would have been observed by chance if this were a homogeneous population (p = 0.001 with a Gaussian distribution; p = 0.005 for a t distribution, df = 4). Ten MDD subjects fell above the gap (“abnormal VTA/SN” subgroup), and 12 fell below it (“normal VTA/SN” subgroup). There were no statistically significant gaps across MDD subjects for the MFB/LNH cluster.

[Figure omitted. See PDF.]

Figure 4. Subgrouping of MDD cohort based on individual FA values.

This figure shows the subgrouping of the MDD cohort based on individual FA values, and double dissociation of VTA/SN and MFB abnormalities across these two subgroups. (A) Scatterplots of FA values from the VTA/SN cluster from MDD subjects (blue) and control subjects (pink). There was a statistically significant gap in the middle of the VTA/SN values for MDD subjects; the values above this gap were all outside the range of control values. When contrasts were calculated separately for MDD subjects in the abnormal VTA/SN (10 MDD and 10 controls) versus the normal VTA/SN (12 MDD and 12 controls) MDD subgroups, (B) abnormal VTA/SN MDD subjects exhibited significantly elevated FA bilaterally in the ventral tegmental area/substantia nigra (VTA/SN), localized to both the SN and the VTA and a trend toward significance in this region in the left hemisphere, while (C) MDD subjects in the normal VTA/SN subgroup did not exhibit any significant FA differences in this region. In contrast, (B) abnormal VTA/SN MDD subjects did not exhibit any significant FA differences overlying the medial forebrain bundle/lateral nucleus of the hypothalamus (MFB/LNH), while (C) normal VTA/SN MDD subjects exhibited a trend toward significant reduction in FA values in this region. Translucent green in images on left indicates a priori AOEs for the (B) VTA/SN and (C) MFB/LNH in the slices shown, which were used to constrain the initial contrast analysis. The color bar indicates the range of p values in this figure, from the threshold (p<0.05) to the order of magnitude of the voxel of peak significance (i.e. smallest p value) in the VTA/SN for elevated FA and from the threshold (p<0.05) to the order of magnitude of the voxel of peak significance in the MFB for reduced FA, color in images is viewed using trilinear interpolation. Warm tones (red, orange, yellow) indicate regions in which MDD subjects exhibited elevated FA relative to control subjects; cool tones (blues) indicate regions in which MDD subjects exhibited reduced FA relative to control subjects. RH: right hemisphere; LH: left hemisphere.

https://doi.org/10.1371/journal.pone.0013945.g004

b. Voxel-based contrasts of FA maps for the two MDD subgroups.

MDD subgroups exhibited two non-overlapping sets of brain microstructural abnormalities when contrasts were conducted separately (Figures 4, 5, 6).

[Figure omitted. See PDF.]

Figure 5. FA difference in the abnormal VTA/SN MDD subgroup in premotor cortex.

This figure shows the significant FA difference in MDD subjects versus matched healthy controls for the voxel-based contrast in the abnormal VTA/SN MDD subgroup for AOEs included in our secondary a priori hypotheses (10 MDD and 10 controls). In addition to elevated FA in the VTA/SN (Figure 4B), MDD subjects in the abnormal VTA/SN subgroup exhibited significantly reduced FA in white matter adjacent to right premotor cortex (PMCwm). The color bar indicates the range of p values in this figure, from the threshold (p<0.05) to the order of magnitude of the voxel of peak significance (i.e. smallest p value) in PMCwm; color in images is viewed using trilinear interpolation. Cool tones (blues) indicate regions in which MDD subjects exhibited reduced FA relative to control subjects. RH: right hemisphere; LH: left hemisphere.

https://doi.org/10.1371/journal.pone.0013945.g005

[Figure omitted. See PDF.]

Figure 6. FA difference in the normal VTA/SN MDD subgroup in secondary a priori areas of evaluation.

This figure shows the FA difference in MDD subjects versus matched healthy controls for the voxel-based contrast in the normal VTA/SN MDD subgroup for AOEs included in our secondary a priori hypotheses (12 MMD and 12 controls). In addition to the trend toward reduced FA in the MFB (Figure 4C), MDD subjects in the normal VTA/SN subgroup also exhibited a trend toward reduced FA in white matter adjacent to the right anterior cingulate (ACCwm). Translucent green in image on left indicates the a priori AOE for ACCwm in the slices shown, which were used to constrain the initial contrast analysis. The color bar indicates the range of p values in this figure, from the threshold (p<0.05) to the order of magnitude of the voxel of peak significance (i.e. smallest p value) in ACCwm; color in images is viewed using trilinear interpolation. Cool tones (blues) indicate regions in which MDD subjects exhibited reduced FA relative to control subjects. RH: right hemisphere; LH: left hemisphere.

https://doi.org/10.1371/journal.pone.0013945.g006

(1) Abnormal VTA/SN MDD subgroup versus matched control subjects (Table 3, Figures 4, 5). For MDD subjects with FA values above the gap, a VTA/SN cluster in the right hemisphere covered a significant proportion of the SN and its peak was centered within the SN (Figure 4B). Voxels in this cluster extended into the lateral VTA. There was also a cluster of significantly elevated FA in the left hemisphere VTA/SN (Figure 4B). This subgroup exhibited significantly reduced FA in white matter underlying right premotor cortex (Figure 5). Mean FA values for clusters are reported in Table S4. This subgroup did not exhibit FA abnormalities within the MFB/LNH a priori AOE (Figure 4C), or within ACC white matter.

(2) Normal VTA/SN MDD subgroup versus matched control subjects (Table 4; Figures 4, 6). The MDD subgroup with VTA/SN FA values below the gap, exhibited a trend toward reduced FA in the left MFB/LNH (Figure 4C), and a trend toward reduced FA in white matter adjacent to right ACC (Figure 6). Mean FA values within clusters are reported in Table S6. This subgroup did not exhibit FA abnormalities within the VTA/SN a priori AOE (Figure 4B) or premotor white matter.

c. Symptom measures in abnormal versus normal VTA/SN MDD subgroups.

Total IDS-SR scores, and scores for items relating to anhedonia, sadness, and psychomotor symptoms did not differ between the normal and abnormal VTA/SN subgroups (Table 5). In contrast, trait anxiety scores (STAI-T) were significantly different between abnormal and normal VTA/SN MDD subgroups [t(2,19) = 2.96, p = 0.0084, normal VTA/SN mean = 51.27±8.96, abnormal VTA/SN mean = 61.80±7.32; Figure 7].

[Figure omitted. See PDF.]

Figure 7. STAI-T (trait anxiety) scores differed between abnormal and normal VTA/SN MDD subgroups.

This figure shows the STAI-T scores for MDD subgroups and for control subjects (n = 10 abnormal VTA/SN MDD; n = 11 normal VTA/SN MDD; n = 19 controls). These scores exhibited a similar distribution to VTA/SN FA values (Figure 4A) across MDD subgroups and controls: Scores for MDD subjects in the abnormal VTA/SN subgroup (left column of blue data points) did not overlap with control subject scores (pink data points), whereas the range of scores for subjects in the normal VTA/SN subgroup (right column of blue data points) overlapped with control scores; scores showed a statistically significant difference between MDD subgroups. The control subject score range also extended below the range for MDD subjects in the normal VTA/SN subgroup, as it did for control VTA/SN FA values.

https://doi.org/10.1371/journal.pone.0013945.g007

[Figure omitted. See PDF.]

Table 5. MDD Symptom Measures in Abnormal Versus Normal VTA/SN MDD Subgroups.

https://doi.org/10.1371/journal.pone.0013945.t005

Discussion

In this study, MDD subjects exhibited brain microstructural differences in subcortical reward/aversion regions, specifically the VTA/SN, compared to control subjects. Individual FA values in the VTA/SN divided the MDD cohort into subgroups with distinct profiles of microstructural abnormalities and different levels of trait anxiety, but no difference in other clinical symptoms of MDD. This subtyping supports a hypothesis that etiology and symptoms of MDD may not match one-to-one, although distinguishing clinical factors may be associated with a given etiology.

Evidence That MDD Is Characterized By Reward/Aversion Circuitry Abnormalities

FA abnormalities in the VTA/SN region.

In the initial voxel-based contrast, FA was elevated in MDD subjects at the border of the SN and cerebral peduncle in a region that, according to anatomic atlases, contains nigrostriatal projection fibers. Because the initial contrast included a number of patients who did not exhibit the abnormality in the VTA/SN, we viewed the voxel-based contrast for the MDD subgroup with abnormal VTA/SN FA as a better reflection of the location and extent of this abnormality. This follow-up analysis revealed a patient/control difference peak centered within the SN, which covered a significant proportion of the SN and extended into the lateral VTA in both hemispheres, suggesting the VTA/SN finding was not a localized effect specific to the nigrostriatal fibers, and that the effect was seen in the nucleus, rather than the adjacent white matter.

Etiological considerations for increased FA in the VTA/SN.

The goal of the current study was to identify regions in which brain microstructure was different in MDD, independent of specific etiology, which can be associated with other known structural abnormalities, clinical or behavioral features in a population. Although DTI evaluates water diffusion properties and therefore can be influenced by many potential etiologies, there are certain etiologies commonly associated with abnormalities in white matter, and emerging evidence for abnormalities underlying FA changes in mixed gray/white matter regions; these can be considered here in relation to generating new hypotheses to test in MDD. Specifically, white matter FA is influenced by the integrity of axons and myelin, as well as directional coherence of axons [74], with reduced FA suggesting either a loss of axonal integrity or coherence. Although the basis for altered FA in gray matter regions has not been as well studied, a number of potential alterations in cell or axonal composition within these regions may alter water diffusion properties; this includes reduced cell density and/or cell loss, altered cell structure, or altered integrity or orientation of the cells or axons projecting out of the region. Altered FA has been previously shown in brainstem nuclei in the presence of known neuronal loss, including altered FA in the SN in Parkinson's Disease [e.g. [65], [66]], confirming the potential for FA to detect pathology-relevant changes in the VTA/SN region.

While neuronal loss, as occurs in Parkinson's Disease, would be expected to lead to reduced FA [65], [66], a reduction in glial density could potentially lead to elevated FA due to an increase in the ratio of axons to cell bodies. Evidence for such an effect is seen in a recent study in a mouse model of Pelizaeus-Merzbacher disease (PMD), in which there is transient astrocytic hypertrophy in females. Hypertrophy in these animals was associated with reduced FA, while FA subsequently increased upon the reversal of hypertrophy [75]; i.e. when glial density decreased, FA increased. Elevated FA in a mixed gray/white matter region has also been directly demonstrated in association with induction of a disease process. A recent study showed that an animal model of febrile seizures exhibited elevated FA in the hippocampus following seizure induction [76]. In humans, such seizures promote hyperexcitability of the limbic system and are accompanied by structural and metabolic abnormalities of the limbic system.

Given that reduced glial density has been demonstrated in MDD in the SGC [77] along with imaging abnormalities in this and other prefrontal gray matter regions in MDD, [34], [78], [79], we suggest that future studies can test the hypothesis that there may also be reduced glial density in the VTA/SN in conjunction with the FA abnormality observed here. Since glia are thought to be critical in synapse formation and synaptic plasticity [80], [81], reduced density and/or loss of glia would have a critical impact on information processing in the reward/aversion circuitry. An alternate possibility is that cell morphology is altered in MDD subjects. Russo and colleagues have shown that chronic drug use leads to a reduction in neuronal size in the VTA in an animal model of chronic opiate addiction [82]. This size reduction was accompanied by a reduction in the rewarding effects of morphine. In MDD, analogous factors, such as chronic stress, might have a similar effect on neurons in this region.

Functional implications of the VTA/SN abnormality for MDD and reward/aversion processing.

Within the VTA/SN region, the VTA predominantly processes reward/aversion information [18], [27], while the SN processes both reward/aversion and motor information [83]. The SN is thought to participate in modulating reward prediction and expectancy [84], as well as prediction of aversive and negatively valenced stimuli [18]. In addition, medial portions of the SN project to and receive projections back from the ventral striatum [30], [70], which processes reward/aversion information [21], while more lateral portions project to the dorsal striatum, which processes motor information [83]. These regions appear to be interconnected in an ascending spiral so that information relating to reward/aversion and motor function is likely mixed in this circuitry [83]. Other efferents of the SN include GABAergic nigro-collicular pathways which have been shown to mediate fear/defense reactions [85].

Taken together, the literature above suggests a hypothesis that altered systems for reward/aversion prediction might lead to an alteration in the capacity of some MDD subjects to assess realistic likelihoods of aversive events [86]. Such individuals, as a consequence, would be likely to show impaired expectancies around negative events and compensate by maintaining higher levels of vigilance/arousal for bad outcomes, which would be clinically observed in the form of higher anxiety [87], [88]. The observation of significantly higher levels of trait anxiety in the abnormal VTA/SN MDD subgroups supports such a hypothesis. Given that the SN and VTA each project to the ventral striatum [89], our findings are also consistent with reports of altered nucleus accumbens function in MDD [2], [32], supporting earlier hypotheses of such an effect [1].

FA abnormalities in premotor cortex.

Independent of a priori regions, MDD subjects exhibited a trend toward reduced FA in white matter underlying premotor cortex in the cohort as a whole, which was significant in the abnormal VTA/SN group. Since a number of tracts run through this region, we cannot be certain that this finding reflects an abnormality in motor fibers. Nevertheless, this finding is of interest in conjunction with the FA abnormality in the SN, given these regions are each key components of motor circuitry, and some MDD subjects in our cohort exhibited psychomotor symptoms.

Evidence For Subtyping of the MDD Cohort Based on DTI Measures and Relationship to Depressive Symptoms

There are three findings arguing that subjects with normal versus abnormal VTA/SN FA fell into two biologically relevant subtypes. The first is the observed gap in the middle of the distribution of FA values for MDD subjects, which coincided with the cut-off for control values. The second is that there was no overlap in regional localization of FA abnormalities for the two subgroups when evaluated separately. The third is the difference in trait anxiety between groups, with higher mean trait anxiety in the abnormal VTA/SN subgroup. The evidence here for MDD subtyping based on features of the VTA/SN resonates with a recent finding of differences in midbrain resting metabolism between groups of patients whose depressive symptoms did versus did not remit in response to antidepressants [31]. Although Milak and colleagues [31] did not test whether there was a division on an individual basis, their findings suggest the possibility that the mictrostructural subgroups in the current study may also predict treatment responsiveness and could be tested as a potential diagnostic/prognostic biomarker. Future studies will be necessary to evaluate the relationship between altered microstructure and treatment responsiveness.

In contrast, no significant differences in two core symptoms of depression (anhedonia and sadness) or overall depression severity were observed between the two MDD subgroups. While this could potentially reflect a power issue, it is also possible that the neural subgroups observed here are an example of how the same illness can arise when related components of the reward/aversion circuitry are “hit” in different places. These data support the idea that it may be the functional system (i.e. reward/aversion) hit, and not the specific etiology of that hit, that determines whether an individual develops MDD [1], [48], [90]. There is, however, potential for vast differences in a subset of clinical symptoms, such as anxiety, and in treatment responsiveness across different biological etiologies.

The higher trait anxiety levels in the abnormal VTA/SN subgroup are consistent with research connecting the VTA/SN to expectancy processing [84], and connects with the idea that abnormal anxiety may reflect altered expectancy [86], [88]. Elevated trait anxiety in this subgroup is also consistent with findings by Berton & colleagues[25] of VTA abnormalities in an animal model of MDD induced by chronic stress. Clinically, our findings in the VTA/SN support hypotheses regarding involvement of the midbrain dopamine reward circuitry in MDD [7], and support investigation of interventions based thereon [91].

Limitations of the Interpretation of Our Findings

Several factors inherent to patient imaging studies must be considered for interpretation of findings such as ours. These include the cohort size in relation to subgrouping, the balance between strengths and weaknesses of our selection of registration methods, the issue of whether microstructural abnormalities are primary or secondary to MDD, and medication status of patients. These factors are discussed below.

Cohort size, subgrouping, and thresholding.

Our search for heterogeneity within the MDD cohort was primarily for the purpose of showing potential reasons previous DTI studies may have exhibited false negatives in subcortical regions. The potential clinical relevance of the subgroups identified is of great interest; however, we emphasize that these findings are preliminary and need to be replicated in larger cohorts. It is also possible that either the cohort size of the subgroups or our conservative thresholds led to false negatives in the current study, particularly in regions that reached the level of a trend (e.g. the MFB). The current study provides justification for future large-scale studies to verify the reproducibility of our findings. Such studies can be conducted using, not only DTI, but also complementary techniques such as post-mortem histological evaluation, which require directed hypotheses such as those generated here.

Registration Methods.

Because our hypotheses focused on specific components of subcortical reward circuitry, we used registration methods which would maximize precise registration of these regions. It is therefore important to emphasize that alternative registration methods might have more accurately or reliably detected abnormalities in cortical white matter regions in MDD subjects. We did, however, focus some of our manual registration points on ventral and medial prefrontal cortical regions to maximize the likelihood that the SGC, FOC, and ACC regions included in our secondary a priori hypotheses would be well registered.

Spatial resolution of our a priori regions.

Small subcortical nuclei and small white matter fiber tracts inherently have lower effective spatial resolution than cortical regions, and this is a limitation particularly when using the spatial resolution of standard DTI sequences (2 mm). In the current study we aimed to optimize the signal at this spatial resolution by using directed registration methods to maximize alignment in our a priori regions. We also note in our Methods that previous studies have successfully detected abnormalities in the VTA/SN region when gross pathology was present [65], [66], suggesting that detection of biological abnormalities in this region should be feasible.

Potential confounds of cardiorespiratory movement.

Because the brainstem is particularly susceptible to motion artifact from arterial or respiratory pulsation [92], it cannot be ruled out that such movements influenced our data. However, in the case of DTI, image acquisition is integrated over a significant period of time (∼10 minutes), and is therefore less susceptible than functional MRI (fMRI) to the effects of such artifacts leading to individual differences. Unless patients exhibited different cardiorespiratory features than controls (or patients differed across subgroups), differences in cardiorespiratory function would be expected to average out, or in the worst case lead to increased variance of the signal in one cohort, reducing the likelihood of a statistically significant finding. Furthermore, all subjects in this study were thoroughly assessed by review of systems and physical exam by a physician, from which no such physical differences were discerned across groups or subgroups.

Are microstructural abnormalities primary or secondary to MDD?

We cannot be certain whether the observed microstructural abnormalities in our MDD subjects were the cause or the result (or both) of the illness or its symptoms. Future studies in larger cohorts and with repeated measures will be necessary to further assess the primacy of abnormalities observed in our study. Such studies will also be needed to evaluate the effects of other factors secondary to MDD, such as tobacco use. Since only three of the 22 MDD subjects in this cohort were current or previous smokers, it is unlikely that our results reflect FA abnormalities secondary to tobacco use. Specifically, this supports the idea that the VTA/SN abnormality could be observed in the absence of tobacco use. Conversely, the control subject who smoked (and had done so for 37 years) exhibited the second lowest VTA/SN FA value (0.24) in the control cohort, indicating that smoking is not sufficient to produce the VTA/SN FA abnormality.

Medication status of MDD subjects and DTI findings.

Because animal literature in depression [93], [94] and more recently human literature in movement disorders [95] suggests that treatment of symptoms of a disorder may influence brain microstructure, it is important to consider whether some of the observations in our study might have been brought about by medication.

Three factors argue that our main findings were not affected by medication status, although future studies will be required to prospectively evaluate whether antidepressants affect brain microstructure as detected by DTI. First, because our patient population was drawn partially from recently diagnosed (i.e. untreated) patients or patients who did not respond to antidepressant medications in the past; the cohort here included a number of patients who had either never taken or were not currently taking antidepressants (Table 1). Second, similar proportions of patients in the normal and abnormal VTA/SN subgroups were on medications (3/10 in the abnormal VTA/SN subgroup; 4/12 in the normal VTA/SN subgroup) at the time of scanning. Third, individual data points indicate there were subjects in each of the two MDD subgroups who had never taken antidepressants so it is unlikely that the MDD DTI subgrouping, which was based on individual rather than group measures, was an effect of medication or medication differences between subgroups.

Conclusion

In this study, MDD subjects exhibited brain microstructural abnormalities in the ventral tegmentum, a primary component of the subcortical reward/aversion circuitry. These abnormalities subdivided the cohort into two subgroups which exhibited similar core depressive symptoms, but differences in trait anxiety. These findings add direct support to the hypothesis that alterations in brain reward/aversion circuitry play a role in the etiology of MDD [1].

Supporting Information

[Figure omitted. See PDF.]

Text S1.

Registration and AOE segmentation methods. Describes the anatomical landmarks used in the manual registration step for FA maps, and the procedures used to segment our a priori areas of evaluation (AOEs).

https://doi.org/10.1371/journal.pone.0013945.s001

(0.05 MB DOC)

Figure S1.

Registration comparison across groups. Examples of (A) coronal and (B) axial slices through average FA maps for control subjects (left column) and MDD subjects (right column) show that registration accuracy of images used in group contrasts was similar between groups and, to illustrate the visibility of anatomy on these images. Coronal images are at y = -18.3 and axial images at z = 3.5 in MNI Talairach coordinates. RH: right hemisphere; LH: left hemisphere.

https://doi.org/10.1371/journal.pone.0013945.s002

(0.91 MB TIF)

Figure S2.

Ventral tegmentum area of evaluation (AOE). This figure shows the primary a priori AOE in the VTA/SN. A priori AOEs were segmented on the average control FA map and then used to constrain the voxels considered in the initial contrast analysis for MDD versus control cohorts. The image on the left in (A) illustrates the visibility of regions and landmarks used for segmentation of the VTA/SN and the three images to the right depict the AOE in translucent green from (left to right) an axial view (the primary orientation in which segmentation was done), a coronal view, and a sagittal view. In (B) the AOE is superimposed on the average MDD FA map at the same locations. RH: right hemisphere; LH: left hemisphere.

https://doi.org/10.1371/journal.pone.0013945.s003

(0.52 MB TIF)

Figure S3.

Medial forebrain bundle area of evaluation (AOE). This figure shows the primary a priori AOE in the medial forebrain bundle (MFB), which also coincided with and included the lateral nucleus of the hypothalamus (LNH). A priori AOEs were segmented on the average control FA map and then used to constrain the voxels considered in the initial contrast analysis for MDD versus control cohorts. The image on the left in (A) illustrates the visibility of regions and landmarks used for segmentation of the MFB/LNH and the three images to the right depict the AOE in translucent green from (left to right) a coronal view (the primary orientation in which segmentation was done), an axial view, and a sagittal view. In (B) the AOE is superimposed on the average MDD FA map at the same locations. RH: right hemisphere; LH: left hemisphere.

https://doi.org/10.1371/journal.pone.0013945.s004

(0.65 MB TIF)

Figure S4.

Amygdalofugal pathway area of evaluation (AOE). This figure shows the primary a priori AOE in the amygdalofugal pathway, which also coincided with and included the substantia inominota. A priori AOEs were segmented on the average control FA map and then used to constrain the voxels considered in the initial contrast analysis for MDD versus control cohorts. The image on the left in (A) illustrates the visibility of regions and landmarks used for segmentation of the amygdalofugal pathway/substantia inominota and the three images to the right depict the AOE in translucent green from (left to right) a coronal view (the primary orientation in which segmentation was done), an axial view, and a sagittal view. In (B) the AOE is superimposed on the average MDD FA map at the same locations. RH: right hemisphere; LH: left hemisphere.

https://doi.org/10.1371/journal.pone.0013945.s005

(0.70 MB TIF)

Figure S5.

Ventral tegmentum group difference cluster. This figure shows the cluster detected within the VTA/SN a priori AOE in the voxel-based contrast analysis of MDD versus control cohorts. This cluster was used to extract values from individual FA maps for use in the follow-up analyses (main text) and permutation test (Dataset SI). (A) Images illustrating the average control FA map (left image), and a priori AOE segmented on that map (right image) that was used to constrain the voxel-based analysis. (B) The cluster of voxels (at right, shown on average control and MDD FA maps, respectively) that met the uncorrected p<0.05 cluster threshold on the p map of the initial MDD versus control voxel-based contrast analysis (left image) and fell within the VTA/SN AOE. Clusters were identified on p maps that had not been smoothed, and thus, p maps are illustrated here without smoothing or interpolation. RH: right hemisphere; LH: left hemisphere.

https://doi.org/10.1371/journal.pone.0013945.s006

(0.46 MB TIF)

Figure S6.

Medial forebrain bundle group difference cluster. This figure shows the cluster detected within the MFB/LNH a priori AOE in the voxel-based contrast analysis of MDD versus control cohorts. This cluster was used to extract values from individual FA maps for use in the follow-up analyses (main text) and permutation test (Dataset SI). (A) Images illustrating the average control FA map (left image), and a priori AOE segmented on that map (right image) that was used to constrain the voxel-based analysis. (B) The cluster of voxels (at right, shown on average control and MDD FA maps, respectively) that met the uncorrected p<0.05 cluster threshold on the p map of the initial MDD versus control voxel-based contrast analysis (left image) and fell within the MFB/LNH AOE. Clusters were identified on p maps that had not been smoothed, and thus, p maps are illustrated here without smoothing or interpolation. RH: right hemisphere; LH: left hemisphere.

https://doi.org/10.1371/journal.pone.0013945.s007

(0.51 MB TIF)

Figure S7.

Examples of secondary a priori areas of evaluation (AOEs). This figure shows slices through each of the five segmentations included in our secondary a priori AOEs, including white matter adjacent to (A) dorsolateral prefrontal cortex (DLPFCwm), (B) anterior cingulate cortex (ACCwm), (C) paracingulate cortex (PACwm), (D) orbitofrontal cortex (FOCwm), and (E) subgenual prefrontal cortex (SGCwm).

https://doi.org/10.1371/journal.pone.0013945.s008

(0.74 MB TIF)

Dataset S1.

Permutation tests. Reports the methods and results of the permutation analysis conducted as a complement to the voxel-wise contrasts of FA maps.

https://doi.org/10.1371/journal.pone.0013945.s009

(0.12 MB DOC)

Dataset S2.

VT abnormalities and laterality. Analysis evaluating whether right hemisphere VT abnormalities in the MDD cohort reflected a loss of normal hemispheric asymmetry in this region or a bilateral change that was simply more significant in one hemisphere than the other.

https://doi.org/10.1371/journal.pone.0013945.s010

(0.03 MB DOC)

Acknowledgments

Within the authorship list, a subset of individuals served on the PGP Publications Committee, which entailed significant management and oversight responsibilities for the project as a whole, and involvement with the development of each manuscript. These individuals were: R. Perlis (Chair of committee), J. Smoller, M. Fava, A. Blood, G. Gasic, H. Breiter.

We would like to thank Mark Vangel, Ph.D. (MGH, HMS) for critical statistical consultation pertaining to the gapping analysis.

Author Contributions

Conceived and designed the experiments: AJB DI RHP BWK MF BR JWS HCB. Performed the experiments: AJB DI RHP DDD BWK MJL SW SL JC JWS GPG HCB. Analyzed the data: AJB DI NM DNK DDD BWK MJL SW SL JC SMH MF JWS GPG HCB. Contributed reagents/materials/analysis tools: AJB NM RHP DNK BWK MJL SMH BR HCB. Wrote the paper: AJB HCB. Contributed revisions to the manuscript drafts produced by AJB and HCB: DI NM RHP DNK DDD BWK MJL SW SL JC SMH MF BR JWS GPG.

Citation: Blood AJ, Iosifescu DV, Makris N, Perlis RH, Kennedy DN, Dougherty DD, et al. (2010) Microstructural Abnormalities in Subcortical Reward Circuitry of Subjects with Major Depressive Disorder. PLoS ONE5(11): e13945. https://doi.org/10.1371/journal.pone.0013945

References

1. Breiter HC, Gasic GP (2004) A General Circuitry of Reward/Aversion and Its Implications for Neuropsychiatric Illness. In: Gazzaniga M, editor. The Cognitive Neurosciences. MIT Press. pp. 1043–1065.HC BreiterGP Gasic2004A General Circuitry of Reward/Aversion and Its Implications for Neuropsychiatric Illness.M. GazzanigaThe Cognitive NeurosciencesMIT Press10431065

2. Pizzagalli DA, Holmes AJ, Dillon DG, Goetz EL, Birk JL, et al. (2009) Reduced caudate and nucleus accumbens response to rewards in unmedicated individuals with major depressive disorder. Am J Psychiatry 166: 702–710.DA PizzagalliAJ HolmesDG DillonEL GoetzJL Birk2009Reduced caudate and nucleus accumbens response to rewards in unmedicated individuals with major depressive disorder.Am J Psychiatry166702710

3. Forbes EE (2009) Where's the fun in that? Broadening the focus on reward function in depression. Biol Psychiatry 66: 199–200.EE Forbes2009Where's the fun in that? Broadening the focus on reward function in depression.Biol Psychiatry66199200

4. McCabe C, Cowen PJ, Harmer CJ (2009) Neural representation of reward in recovered depressed patients. Psychopharmacology (Berl) 205: 667–677.C. McCabePJ CowenCJ Harmer2009Neural representation of reward in recovered depressed patients.Psychopharmacology (Berl)205667677

5. Schlaepfer TE, Cohen MX, Frick C, Kosel M, Brodesser D, et al. (2007) Deep Brain Stimulation to Reward Circuitry Alleviates Anhedonia in Refractory Major Depression. Neuropsychopharmacology. TE SchlaepferMX CohenC. FrickM. KoselD. Brodesser2007Deep Brain Stimulation to Reward Circuitry Alleviates Anhedonia in Refractory Major Depression.Neuropsychopharmacology

6. Tremblay LK, Naranjo CA, Graham SJ, Herrmann N, Mayberg HS, et al. (2005) Functional neuroanatomical substrates of altered reward processing in major depressive disorder revealed by a dopaminergic probe. Arch Gen Psychiatry 62: 1228–1236.LK TremblayCA NaranjoSJ GrahamN. HerrmannHS Mayberg2005Functional neuroanatomical substrates of altered reward processing in major depressive disorder revealed by a dopaminergic probe.Arch Gen Psychiatry6212281236

7. Nestler EJ, Carlezon WA Jr (2006) The mesolimbic dopamine reward circuit in depression. Biol Psychiatry 59: 1151–1159.EJ NestlerWA Carlezon Jr2006The mesolimbic dopamine reward circuit in depression.Biol Psychiatry5911511159

8. Breiter HC, Rosen BR (1999) Functional magnetic resonance imaging of brain reward circuitry in the human. Ann N Y Acad Sci 877: 523–547.HC BreiterBR Rosen1999Functional magnetic resonance imaging of brain reward circuitry in the human.Ann N Y Acad Sci877523547

9. Breiter HC, Etcoff NL, Whalen PJ, Kennedy WA, Rauch SL, et al. (1996) Response and habituation of the human amygdala during visual processing of facial expression. Neuron 17: 875–887.HC BreiterNL EtcoffPJ WhalenWA KennedySL Rauch1996Response and habituation of the human amygdala during visual processing of facial expression.Neuron17875887

10. Blood AJ, Zatorre RJ, Bermudez P, Evans AC (1999) Emotional responses to pleasant and unpleasant music correlate with activity in paralimbic brain regions. Nat Neurosci 2: 382–387.AJ BloodRJ ZatorreP. BermudezAC Evans1999Emotional responses to pleasant and unpleasant music correlate with activity in paralimbic brain regions.Nat Neurosci2382387

11. Breiter HC, Aharon I, Kahneman D, Dale A, Shizgal P (2001) Functional imaging of neural responses to expectancy and experience of monetary gains and losses. Neuron 30: 619–639.HC BreiterI. AharonD. KahnemanA. DaleP. Shizgal2001Functional imaging of neural responses to expectancy and experience of monetary gains and losses.Neuron30619639

12. Becerra L, Breiter HC, Wise R, Gonzalez RG, Borsook D (2001) Reward circuitry activation by noxious thermal stimuli. Neuron 32: 927–946.L. BecerraHC BreiterR. WiseRG GonzalezD. Borsook2001Reward circuitry activation by noxious thermal stimuli.Neuron32927946

13. Delgado MR, Nystrom LE, Fissell C, Noll DC, Fiez JA (2000) Tracking the hemodynamic responses to reward and punishment in the striatum. J Neurophysiol 84: 3072–3077.MR DelgadoLE NystromC. FissellDC NollJA Fiez2000Tracking the hemodynamic responses to reward and punishment in the striatum.J Neurophysiol8430723077

14. Jensen J, McIntosh AR, Crawley AP, Mikulis DJ, Remington G, et al. (2003) Direct activation of the ventral striatum in anticipation of aversive stimuli. Neuron 40: 1251–1257.J. JensenAR McIntoshAP CrawleyDJ MikulisG. Remington2003Direct activation of the ventral striatum in anticipation of aversive stimuli.Neuron4012511257

15. Roitman MF, Wheeler RA, Carelli RM (2005) Nucleus accumbens neurons are innately tuned for rewarding and aversive taste stimuli, encode their predictors, and are linked to motor output. Neuron 45: 587–597.MF RoitmanRA WheelerRM Carelli2005Nucleus accumbens neurons are innately tuned for rewarding and aversive taste stimuli, encode their predictors, and are linked to motor output.Neuron45587597

16. Menon M, Jensen J, Vitcu I, Graff-Guerrero A, Crawley A, et al. (2007) Temporal difference modeling of the blood-oxygen level dependent response during aversive conditioning in humans: effects of dopaminergic modulation. Biol Psychiatry 62: 765–772.M. MenonJ. JensenI. VitcuA. Graff-GuerreroA. Crawley2007Temporal difference modeling of the blood-oxygen level dependent response during aversive conditioning in humans: effects of dopaminergic modulation.Biol Psychiatry62765772

17. Diaconescu AO, Menon M, Jensen J, Kapur S, McIntosh ARDopamine-induced changes in neural network patterns supporting aversive conditioning. Brain Res 1313: 143–161.AO DiaconescuM. MenonJ. JensenS. KapurAR McIntoshDopamine-induced changes in neural network patterns supporting aversive conditioning.Brain Res1313143161

18. Matsumoto M, Hikosaka O (2009) Two types of dopamine neuron distinctly convey positive and negative motivational signals. Nature 459: 837–841.M. MatsumotoO. Hikosaka2009Two types of dopamine neuron distinctly convey positive and negative motivational signals.Nature459837841

19. Leknes S, Tracey I (2008) A common neurobiology for pain and pleasure. Nat Rev Neurosci 9: 314–320.S. LeknesI. Tracey2008A common neurobiology for pain and pleasure.Nat Rev Neurosci9314320

20. Levita L, Hare TA, Voss HU, Glover G, Ballon DJ, et al. (2009) The bivalent side of the nucleus accumbens. Neuroimage 44: 1178–1187.L. LevitaTA HareHU VossG. GloverDJ Ballon2009The bivalent side of the nucleus accumbens.Neuroimage4411781187

21. Carlezon WA Jr, Thomas MJ (2009) Biological substrates of reward and aversion: a nucleus accumbens activity hypothesis. Neuropharmacology 56: Suppl 1122–132.WA Carlezon JrMJ Thomas2009Biological substrates of reward and aversion: a nucleus accumbens activity hypothesis.Neuropharmacology56Suppl 1122132

22. Kim BW, Kennedy DN, Lehar J, Lee MJ, Blood AJ, et al. Recurrent, robust and scalable patterns underlie human approach and avoidance. PLoS One 5: e10613.BW KimDN KennedyJ. LeharMJ LeeAJ BloodRecurrent, robust and scalable patterns underlie human approach and avoidance.PLoS One5e10613

23. Bogdan R, Pizzagalli DA (2006) Acute stress reduces reward responsiveness: implications for depression. Biol Psychiatry 60: 1147–1154.R. BogdanDA Pizzagalli2006Acute stress reduces reward responsiveness: implications for depression.Biol Psychiatry6011471154

24. Breiter HC, Gasic GP, Makris , N (2006) Imaging the neural systems for motivated behavior and their dysfunction in neuropsychiatric illness. In: Deisboeck T, Kresh J, editors. Complex systems science in biomedicine. New York: Springer Verlag. pp. 763–810.HC BreiterGP GasicMakrisN2006Imaging the neural systems for motivated behavior and their dysfunction in neuropsychiatric illness.T. DeisboeckJ. KreshComplex systems science in biomedicineNew YorkSpringer Verlag763810

25. Berton O, McClung CA, Dileone RJ, Krishnan V, Renthal W, et al. (2006) Essential role of BDNF in the mesolimbic dopamine pathway in social defeat stress. Science 311: 864–868.O. BertonCA McClungRJ DileoneV. KrishnanW. Renthal2006Essential role of BDNF in the mesolimbic dopamine pathway in social defeat stress.Science311864868

26. Olds J, Milner P (1954) Positive reinforcement produced by electrical stimulation of septal area and other regions of rat brain. J Comp Physiol Psychol 47: 419–427.J. OldsP. Milner1954Positive reinforcement produced by electrical stimulation of septal area and other regions of rat brain.J Comp Physiol Psychol47419427

27. Wise RA, Bozarth MA (1985) Brain mechanisms of drug reward and euphoria. Psychiatr Med 3: 445–460.RA WiseMA Bozarth1985Brain mechanisms of drug reward and euphoria.Psychiatr Med3445460

28. Wise RA (1978) Catecholamine theories of reward: a critical review. Brain Res 152: 215–247.RA Wise1978Catecholamine theories of reward: a critical review.Brain Res152215247

29. Manji HK, Drevets WC, Charney DS (2001) The cellular neurobiology of depression. Nat Med 7: 541–547.HK ManjiWC DrevetsDS Charney2001The cellular neurobiology of depression.Nat Med7541547

30. Breiter HC, Gollub RL, Weisskoff RM, Kennedy DN, Makris N, et al. (1997) Acute effects of cocaine on human brain activity and emotion. Neuron 19: 591–611.HC BreiterRL GollubRM WeisskoffDN KennedyN. Makris1997Acute effects of cocaine on human brain activity and emotion.Neuron19591611

31. Milak MS, Parsey RV, Lee L, Oquendo MA, Olvet DM, et al. (2009) Pretreatment regional brain glucose uptake in the midbrain on PET may predict remission from a major depressive episode after three months of treatment. Psychiatry Res 173: 63–70.MS MilakRV ParseyL. LeeMA OquendoDM Olvet2009Pretreatment regional brain glucose uptake in the midbrain on PET may predict remission from a major depressive episode after three months of treatment.Psychiatry Res1736370

32. Wacker J, Dillon DG, Pizzagalli DA (2009) The role of the nucleus accumbens and rostral anterior cingulate cortex in anhedonia: integration of resting EEG, fMRI, and volumetric techniques. Neuroimage 46: 327–337.J. WackerDG DillonDA Pizzagalli2009The role of the nucleus accumbens and rostral anterior cingulate cortex in anhedonia: integration of resting EEG, fMRI, and volumetric techniques.Neuroimage46327337

33. Nobuhara K, Okugawa G, Sugimoto T, Minami T, Tamagaki C, et al. (2006) Frontal white matter anisotropy and symptom severity of late-life depression: a magnetic resonance diffusion tensor imaging study. J Neurol Neurosurg Psychiatry 77: 120–122.K. NobuharaG. OkugawaT. SugimotoT. MinamiC. Tamagaki2006Frontal white matter anisotropy and symptom severity of late-life depression: a magnetic resonance diffusion tensor imaging study.J Neurol Neurosurg Psychiatry77120122

34. Drevets WC, Price JL, Simpson JR Jr, Todd RD, Reich T, et al. (1997) Subgenual prefrontal cortex abnormalities in mood disorders. Nature 386: 824–827.WC DrevetsJL PriceJR Simpson JrRD ToddT. Reich1997Subgenual prefrontal cortex abnormalities in mood disorders.Nature386824827

35. Drevets WC, Ongur D, Price JL (1998) Reduced glucose metabolism in the subgenual prefrontal cortex in unipolar depression. Mol Psychiatry 3: 190–191.WC DrevetsD. OngurJL Price1998Reduced glucose metabolism in the subgenual prefrontal cortex in unipolar depression.Mol Psychiatry3190191

36. Taylor WD, MacFall JR, Payne ME, McQuoid DR, Provenzale JM, et al. (2004) Late-life depression and microstructural abnormalities in dorsolateral prefrontal cortex white matter. Am J Psychiatry 161: 1293–1296.WD TaylorJR MacFallME PayneDR McQuoidJM Provenzale2004Late-life depression and microstructural abnormalities in dorsolateral prefrontal cortex white matter.Am J Psychiatry16112931296

37. Keedwell PA, Andrew C, Williams SC, Brammer MJ, Phillips ML (2005) The neural correlates of anhedonia in major depressive disorder. Biol Psychiatry 58: 843–853.PA KeedwellC. AndrewSC WilliamsMJ BrammerML Phillips2005The neural correlates of anhedonia in major depressive disorder.Biol Psychiatry58843853

38. Alexopoulos GS, Kiosses DN, Choi SJ, Murphy CF, Lim KO (2002) Frontal white matter microstructure and treatment response of late-life depression: a preliminary study. Am J Psychiatry 159: 1929–1932.GS AlexopoulosDN KiossesSJ ChoiCF MurphyKO Lim2002Frontal white matter microstructure and treatment response of late-life depression: a preliminary study.Am J Psychiatry15919291932

39. Bae JN, MacFall JR, Krishnan KR, Payne ME, Steffens DC, et al. (2006) Dorsolateral prefrontal cortex and anterior cingulate cortex white matter alterations in late-life depression. Biol Psychiatry 60: 1356–1363.JN BaeJR MacFallKR KrishnanME PayneDC Steffens2006Dorsolateral prefrontal cortex and anterior cingulate cortex white matter alterations in late-life depression.Biol Psychiatry6013561363

40. Frodl T, Meisenzahl EM, Zetzsche T, Born C, Groll C, et al. (2002) Hippocampal changes in patients with a first episode of major depression. Am J Psychiatry 159: 1112–1118.T. FrodlEM MeisenzahlT. ZetzscheC. BornC. Groll2002Hippocampal changes in patients with a first episode of major depression.Am J Psychiatry15911121118

41. Mintun MA, Sheline YI, Moerlein SM, Vlassenko AG, Huang Y, et al. (2004) Decreased hippocampal 5-HT2A receptor binding in major depressive disorder: in vivo measurement with [18F]altanserin positron emission tomography. Biol Psychiatry 55: 217–224.MA MintunYI ShelineSM MoerleinAG VlassenkoY. Huang2004Decreased hippocampal 5-HT2A receptor binding in major depressive disorder: in vivo measurement with [18F]altanserin positron emission tomography.Biol Psychiatry55217224

42. Sheline YI, Wang PW, Gado MH, Csernansky JG, Vannier MW (1996) Hippocampal atrophy in recurrent major depression. Proc Natl Acad Sci U S A 93: 3908–3913.YI ShelinePW WangMH GadoJG CsernanskyMW Vannier1996Hippocampal atrophy in recurrent major depression.Proc Natl Acad Sci U S A9339083913

43. Bremner JD, Narayan M, Anderson ER, Staib LH, Miller HL, et al. (2000) Hippocampal volume reduction in major depression. Am J Psychiatry 157: 115–118.JD BremnerM. NarayanER AndersonLH StaibHL Miller2000Hippocampal volume reduction in major depression.Am J Psychiatry157115118

44. Saxena S, Brody AL, Ho ML, Alborzian S, Ho MK, et al. (2001) Cerebral metabolism in major depression and obsessive-compulsive disorder occurring separately and concurrently. Biol Psychiatry 50: 159–170.S. SaxenaAL BrodyML HoS. AlborzianMK Ho2001Cerebral metabolism in major depression and obsessive-compulsive disorder occurring separately and concurrently.Biol Psychiatry50159170

45. Edwards E, Kornrich W, van Houtten P, Henn FA (1992) In vitro neurotransmitter release in an animal model of depression. Neurochem Int 21: 29–35.E. EdwardsW. KornrichP. van HouttenFA Henn1992In vitro neurotransmitter release in an animal model of depression.Neurochem Int212935

46. Drevets WC, Price JL, Bardgett ME, Reich T, Todd RD, et al. (2002) Glucose metabolism in the amygdala in depression: relationship to diagnostic subtype and plasma cortisol levels. Pharmacol Biochem Behav 71: 431–447.WC DrevetsJL PriceME BardgettT. ReichRD Todd2002Glucose metabolism in the amygdala in depression: relationship to diagnostic subtype and plasma cortisol levels.Pharmacol Biochem Behav71431447

47. Frodl T, Meisenzahl E, Zetzsche T, Bottlender R, Born C, et al. (2002) Enlargement of the amygdala in patients with a first episode of major depression. Biol Psychiatry 51: 708–714.T. FrodlE. MeisenzahlT. ZetzscheR. BottlenderC. Born2002Enlargement of the amygdala in patients with a first episode of major depression.Biol Psychiatry51708714

48. Mayberg HS (2003) Modulating dysfunctional limbic-cortical circuits in depression: towards development of brain-based algorithms for diagnosis and optimised treatment. Br Med Bull 65: 193–207.HS Mayberg2003Modulating dysfunctional limbic-cortical circuits in depression: towards development of brain-based algorithms for diagnosis and optimised treatment.Br Med Bull65193207

49. Fava M, Martinez JM, Greist J, Marangell LB, Brown E, et al. (2007) The efficacy and tolerability of duloxetine in the treatment of anxious versus non-anxious depression: a post-hoc analysis of an open-label outpatient study. Ann Clin Psychiatry 19: 187–195.M. FavaJM MartinezJ. GreistLB MarangellE. Brown2007The efficacy and tolerability of duloxetine in the treatment of anxious versus non-anxious depression: a post-hoc analysis of an open-label outpatient study.Ann Clin Psychiatry19187195

50. Sobin C, Sackeim HA (1997) Psychomotor symptoms of depression. Am J Psychiatry 154: 4–17.C. SobinHA Sackeim1997Psychomotor symptoms of depression.Am J Psychiatry154417

51. Dougherty DD, Weiss AP, Cosgrove GR, Alpert NM, Cassem EH, et al. (2003) Cerebral metabolic correlates as potential predictors of response to anterior cingulotomy for treatment of major depression. J Neurosurg 99: 1010–1017.DD DoughertyAP WeissGR CosgroveNM AlpertEH Cassem2003Cerebral metabolic correlates as potential predictors of response to anterior cingulotomy for treatment of major depression.J Neurosurg9910101017

52. Langguth B, Wiegand R, Kharraz A, Landgrebe M, Marienhagen J, et al. (2007) Pre-treatment anterior cingulate activity as a predictor of antidepressant response to repetitive transcranial magnetic stimulation (rTMS). Neuro Endocrinol Lett 28: 633–638.B. LangguthR. WiegandA. KharrazM. LandgrebeJ. Marienhagen2007Pre-treatment anterior cingulate activity as a predictor of antidepressant response to repetitive transcranial magnetic stimulation (rTMS).Neuro Endocrinol Lett28633638

53. Grimm S, Beck J, Schuepbach D, Hell D, Boesiger P, et al. (2008) Imbalance between left and right dorsolateral prefrontal cortex in major depression is linked to negative emotional judgment: an fMRI study in severe major depressive disorder. Biol Psychiatry 63: 369–376.S. GrimmJ. BeckD. SchuepbachD. HellP. Boesiger2008Imbalance between left and right dorsolateral prefrontal cortex in major depression is linked to negative emotional judgment: an fMRI study in severe major depressive disorder.Biol Psychiatry63369376

54. Surguladze S, Brammer MJ, Keedwell P, Giampietro V, Young AW, et al. (2005) A differential pattern of neural response toward sad versus happy facial expressions in major depressive disorder. Biol Psychiatry 57: 201–209.S. SurguladzeMJ BrammerP. KeedwellV. GiampietroAW Young2005A differential pattern of neural response toward sad versus happy facial expressions in major depressive disorder.Biol Psychiatry57201209

55. Phillips ML, Drevets WC, Rauch SL, Lane R (2003) Neurobiology of emotion perception II: Implications for major psychiatric disorders. Biol Psychiatry 54: 515–528.ML PhillipsWC DrevetsSL RauchR. Lane2003Neurobiology of emotion perception II: Implications for major psychiatric disorders.Biol Psychiatry54515528

56. Abercrombie HC, Schaefer SM, Larson CL, Oakes TR, Lindgren KA, et al. (1998) Metabolic rate in the right amygdala predicts negative affect in depressed patients. Neuroreport 9: 3301–3307.HC AbercrombieSM SchaeferCL LarsonTR OakesKA Lindgren1998Metabolic rate in the right amygdala predicts negative affect in depressed patients.Neuroreport933013307

57. Fava M, Uebelacker LA, Alpert JE, Nierenberg AA, Pava JA, et al. (1997) Major depressive subtypes and treatment response. Biol Psychiatry 42: 568–576.M. FavaLA UebelackerJE AlpertAA NierenbergJA Pava1997Major depressive subtypes and treatment response.Biol Psychiatry42568576

58. Treadway MT, Buckholtz JW, Schwartzman AN, Lambert WE, Zald DH (2009) Worth the ‘EEfRT’? The effort expenditure for rewards task as an objective measure of motivation and anhedonia. PLoS One 4: e6598.MT TreadwayJW BuckholtzAN SchwartzmanWE LambertDH Zald2009Worth the ‘EEfRT’? The effort expenditure for rewards task as an objective measure of motivation and anhedonia.PLoS One4e6598

59. First MB, Spitzer RL, Gibbon M, Williams JBW (1996) Structured Clinical Interview for DSM-IV Axis I Disorders. Patient edition (SCID-I/P) ed. New York: Biometrics Research Department, New York State Psychiatric Institute. MB FirstRL SpitzerM. GibbonJBW Williams1996Structured Clinical Interview for DSM-IV Axis I Disorders. Patient edition (SCID-I/P) ed.New YorkBiometrics Research Department, New York State Psychiatric Institute

60. Hamilton M (1967) Development of a rating scale for primary depressive illness. Br J Soc Clin Psychol 6: 278–296.M. Hamilton1967Development of a rating scale for primary depressive illness.Br J Soc Clin Psychol6278296

61. van der Kouwe AJ, Benner T, Fischl B, Schmitt F, Salat DH, et al. (2005) On-line automatic slice positioning for brain MR imaging. Neuroimage 27: 222–230.AJ van der KouweT. BennerB. FischlF. SchmittDH Salat2005On-line automatic slice positioning for brain MR imaging.Neuroimage27222230

62. Deichmann R, Gottfried JA, Hutton C, Turner R (2003) Optimized EPI for fMRI studies of the orbitofrontal cortex. Neuroimage 19: 430–441.R. DeichmannJA GottfriedC. HuttonR. Turner2003Optimized EPI for fMRI studies of the orbitofrontal cortex.Neuroimage19430441

63. Salat DH, Tuch DS, Greve DN, van der Kouwe AJ, Hevelone ND, et al. (2005) Age-related alterations in white matter microstructure measured by diffusion tensor imaging. Neurobiol Aging 26: 1215–1227.DH SalatDS TuchDN GreveAJ van der KouweND Hevelone2005Age-related alterations in white matter microstructure measured by diffusion tensor imaging.Neurobiol Aging2612151227

64. Pierpaoli C, Basser PJ (1996) Toward a quantitative assessment of diffusion anisotropy. Magn Reson Med 36: 893–906.C. PierpaoliPJ Basser1996Toward a quantitative assessment of diffusion anisotropy.Magn Reson Med36893906

65. Boska MD, Hasan KM, Kibuule D, Banerjee R, McIntyre E, et al. (2007) Quantitative diffusion tensor imaging detects dopaminergic neuronal degeneration in a murine model of Parkinson's disease. Neurobiol Dis 26: 590–596.MD BoskaKM HasanD. KibuuleR. BanerjeeE. McIntyre2007Quantitative diffusion tensor imaging detects dopaminergic neuronal degeneration in a murine model of Parkinson's disease.Neurobiol Dis26590596

66. Chan LL, Rumpel H, Yap K, Lee E, Loo HV, et al. (2007) Case control study of diffusion tensor imaging in Parkinson's disease. J Neurol Neurosurg Psychiatry 78: 1383–1386.LL ChanH. RumpelK. YapE. LeeHV Loo2007Case control study of diffusion tensor imaging in Parkinson's disease.J Neurol Neurosurg Psychiatry7813831386

67. Jenkinson M, Smith S (2001) A global optimisation method for robust affine registration of brain images. Med Image Anal 5: 143–156.M. JenkinsonS. Smith2001A global optimisation method for robust affine registration of brain images.Med Image Anal5143156

68. Smith SM, Jenkinson M, Johansen-Berg H, Rueckert D, Nichols TE, et al. (2006) Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. Neuroimage 31: 1487–1505.SM SmithM. JenkinsonH. Johansen-BergD. RueckertTE Nichols2006Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data.Neuroimage3114871505

69. Tuch DS, Salat DH, Wisco JJ, Zaleta AK, Hevelone ND, et al. (2005) Choice reaction time performance correlates with diffusion anisotropy in white matter pathways supporting visuospatial attention. Proc Natl Acad Sci U S A 102: 12212–12217.DS TuchDH SalatJJ WiscoAK ZaletaND Hevelone2005Choice reaction time performance correlates with diffusion anisotropy in white matter pathways supporting visuospatial attention.Proc Natl Acad Sci U S A1021221212217

70. Aharon I, Etcoff N, Ariely D, Chabris CF, O'Connor E, et al. (2001) Beautiful faces have variable reward value: fMRI and behavioral evidence. Neuron 32: 537–551.I. AharonN. EtcoffD. ArielyCF ChabrisE. O'Connor2001Beautiful faces have variable reward value: fMRI and behavioral evidence.Neuron32537551

71. Makris N, Gasic GP, Kennedy DN, Hodge SM, Kaiser JR, et al. (2008) Cortical thickness abnormalities in cocaine addiction—a reflection of both drug use and a pre-existing disposition to drug abuse? Neuron 60: 174–188.N. MakrisGP GasicDN KennedySM HodgeJR Kaiser2008Cortical thickness abnormalities in cocaine addiction—a reflection of both drug use and a pre-existing disposition to drug abuse?Neuron60174188

72. Wainer H (1978) Gapping. Psychometrika 43: 203–212.H. Wainer1978Gapping.Psychometrika43203212

73. Fava M, Rush AJ, Alpert JE, Balasubramani GK, Wisniewski SR, et al. (2008) Difference in treatment outcome in outpatients with anxious versus nonanxious depression: a STAR*D report. Am J Psychiatry 165: 342–351.M. FavaAJ RushJE AlpertGK BalasubramaniSR Wisniewski2008Difference in treatment outcome in outpatients with anxious versus nonanxious depression: a STAR*D report.Am J Psychiatry165342351

74. Beaulieu C (2002) The basis of anisotropic water diffusion in the nervous system - a technical review. NMR Biomed 15: 435–455.C. Beaulieu2002The basis of anisotropic water diffusion in the nervous system - a technical review.NMR Biomed15435455

75. Harsan LA, Poulet P, Guignard B, Parizel N, Skoff RP, et al. (2007) Astrocytic hypertrophy in dysmyelination influences the diffusion anisotropy of white matter. J Neurosci Res 85: 935–944.LA HarsanP. PouletB. GuignardN. ParizelRP Skoff2007Astrocytic hypertrophy in dysmyelination influences the diffusion anisotropy of white matter.J Neurosci Res85935944

76. Jansen JF, Lemmens EM, Strijkers GJ, Prompers JJ, Schijns OE, et al. (2008) Short- and long-term limbic abnormalities after experimental febrile seizures. Neurobiol Dis. JF JansenEM LemmensGJ StrijkersJJ PrompersOE Schijns2008Short- and long-term limbic abnormalities after experimental febrile seizures.Neurobiol Dis

77. Ongur D, Drevets WC, Price JL (1998) Glial reduction in the subgenual prefrontal cortex in mood disorders. Proc Natl Acad Sci U S A 95: 13290–13295.D. OngurWC DrevetsJL Price1998Glial reduction in the subgenual prefrontal cortex in mood disorders.Proc Natl Acad Sci U S A951329013295

78. Ballmaier M, Toga AW, Blanton RE, Sowell ER, Lavretsky H, et al. (2004) Anterior cingulate, gyrus rectus, and orbitofrontal abnormalities in elderly depressed patients: an MRI-based parcellation of the prefrontal cortex. Am J Psychiatry 161: 99–108.M. BallmaierAW TogaRE BlantonER SowellH. Lavretsky2004Anterior cingulate, gyrus rectus, and orbitofrontal abnormalities in elderly depressed patients: an MRI-based parcellation of the prefrontal cortex.Am J Psychiatry16199108

79. Botteron KN, Raichle ME, Drevets WC, Heath AC, Todd RD (2002) Volumetric reduction in left subgenual prefrontal cortex in early onset depression. Biol Psychiatry 51: 342–344.KN BotteronME RaichleWC DrevetsAC HeathRD Todd2002Volumetric reduction in left subgenual prefrontal cortex in early onset depression.Biol Psychiatry51342344

80. Bains JS, Oliet SH (2007) Glia: they make your memories stick! Trends Neurosci.JS BainsSH Oliet2007Glia: they make your memories stick! Trends Neurosci.

81. Ullian EM, Sapperstein SK, Christopherson KS, Barres BA (2001) Control of synapse number by glia. Science 291: 657–661.EM UllianSK SappersteinKS ChristophersonBA Barres2001Control of synapse number by glia.Science291657661

82. Russo SJ, Bolanos CA, Theobald DE, DeCarolis NA, Renthal W, et al. (2007) IRS2-Akt pathway in midbrain dopamine neurons regulates behavioral and cellular responses to opiates. Nat Neurosci 10: 93–99.SJ RussoCA BolanosDE TheobaldNA DeCarolisW. Renthal2007IRS2-Akt pathway in midbrain dopamine neurons regulates behavioral and cellular responses to opiates.Nat Neurosci109399

83. Haber SN, Fudge JL, McFarland NR (2000) Striatonigrostriatal pathways in primates form an ascending spiral from the shell to the dorsolateral striatum. J Neurosci 20: 2369–2382.SN HaberJL FudgeNR McFarland2000Striatonigrostriatal pathways in primates form an ascending spiral from the shell to the dorsolateral striatum.J Neurosci2023692382

84. Schultz W, Dayan P, Montague PR (1997) A neural substrate of prediction and reward. Science 275: 1593–1599.W. SchultzP. DayanPR Montague1997A neural substrate of prediction and reward.Science27515931599

85. Nobre MJ, Lopes MG, Brandao ML (2004) Defense reaction mediated by NMDA mechanisms in the inferior colliculus is modulated by GABAergic nigro-collicular pathways. Brain Res 999: 124–131.MJ NobreMG LopesML Brandao2004Defense reaction mediated by NMDA mechanisms in the inferior colliculus is modulated by GABAergic nigro-collicular pathways.Brain Res999124131

86. Wiedemann G, Pauli P, Dengler W (2001) A priori expectancy bias in patients with panic disorder. J Anxiety Disord 15: 401–412.G. WiedemannP. PauliW. Dengler2001A priori expectancy bias in patients with panic disorder.J Anxiety Disord15401412

87. Compton RJ, Carp J, Chaddock L, Fineman SL, Quandt LC, et al. (2007) Anxiety and error monitoring: Increased error sensitivity or altered expectations? Brain Cogn. RJ ComptonJ. CarpL. ChaddockSL FinemanLC Quandt2007Anxiety and error monitoring: Increased error sensitivity or altered expectations?Brain Cogn

88. Chan CK, Lovibond PF (1996) Expectancy bias in trait anxiety. J Abnorm Psychol 105: 637–647.CK ChanPF Lovibond1996Expectancy bias in trait anxiety.J Abnorm Psychol105637647

89. Wise RA (2009) Roles for nigrostriatal—not just mesocorticolimbic—dopamine in reward and addiction. Trends Neurosci 32: 517–524.RA Wise2009Roles for nigrostriatal—not just mesocorticolimbic—dopamine in reward and addiction.Trends Neurosci32517524

90. Blood AJ (2008) New hypotheses about postural control support the notion that all dystonias are manifestations of excessive brain postural function. Biosci Hypotheses 1: 14–25.AJ Blood2008New hypotheses about postural control support the notion that all dystonias are manifestations of excessive brain postural function.Biosci Hypotheses11425

91. Guiard BP, El Mansari M, Blier P (2009) Prospect of a Dopamine Contribution in the Next Generation of Antidepressant Drugs: The Triple Reuptake Inhibitors. Curr Drug Targets. BP GuiardM. El MansariP. Blier2009Prospect of a Dopamine Contribution in the Next Generation of Antidepressant Drugs: The Triple Reuptake Inhibitors.Curr Drug Targets

92. Guimaraes AR, Melcher JR, Talavage TM, Baker JR, Ledden P, et al. (1998) Imaging subcortical auditory activity in humans. Hum Brain Mapp 6: 33–41.AR GuimaraesJR MelcherTM TalavageJR BakerP. Ledden1998Imaging subcortical auditory activity in humans.Hum Brain Mapp63341

93. Malberg JE, Eisch AJ, Nestler EJ, Duman RS (2000) Chronic antidepressant treatment increases neurogenesis in adult rat hippocampus. J Neurosci 20: 9104–9110.JE MalbergAJ EischEJ NestlerRS Duman2000Chronic antidepressant treatment increases neurogenesis in adult rat hippocampus.J Neurosci2091049110

94. Santarelli L, Saxe M, Gross C, Surget A, Battaglia F, et al. (2003) Requirement of hippocampal neurogenesis for the behavioral effects of antidepressants. Science 301: 805–809.L. SantarelliM. SaxeC. GrossA. SurgetF. Battaglia2003Requirement of hippocampal neurogenesis for the behavioral effects of antidepressants.Science301805809

95. Blood AJ, Tuch DS, Makris N, Makhlouf ML, Sudarsky LR, et al. (2006) White matter abnormalities in dystonia normalize after botulinum toxin treatment. Neuroreport 17: 1251–1255.AJ BloodDS TuchN. MakrisML MakhloufLR Sudarsky2006White matter abnormalities in dystonia normalize after botulinum toxin treatment.Neuroreport1712511255

Word count: 13733

Show less

© 2010 Blood et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License: https://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.

Abstract

Translate

Background

Previous studies of major depressive disorder (MDD) have focused on abnormalities in the prefrontal cortex and medial temporal regions. There has been little investigation in MDD of midbrain and subcortical regions central to reward/aversion function, such as the ventral tegmental area/substantia nigra (VTA/SN), and medial forebrain bundle (MFB).

Methodology/Principal Findings

We investigated the microstructural integrity of this circuitry using diffusion tensor imaging (DTI) in 22 MDD subjects and compared them with 22 matched healthy control subjects. Fractional anisotropy (FA) values were increased in the right VT and reduced in dorsolateral prefrontal white matter in MDD subjects. Follow-up analysis suggested two distinct subgroups of MDD patients, which exhibited non-overlapping abnormalities in reward/aversion circuitry. The MDD subgroup with abnormal FA values in VT exhibited significantly greater trait anxiety than the subgroup with normal FA values in VT, but the subgroups did not differ in levels of anhedonia, sadness, or overall depression severity.

Conclusions/Significance

These findings suggest that MDD may be associated with abnormal microstructure in brain reward/aversion regions, and that there may be at least two subtypes of microstructural abnormalities which each impact core symptoms of depression.

Details

Title

Microstructural Abnormalities in Subcortical Reward Circuitry of Subjects with Major Depressive Disorder

Author

Blood, Anne J; Iosifescu, Dan V; Makris, Nikos; Perlis, Roy H; Kennedy, David N; Dougherty, Darin D; Byoung Woo Kim; Lee, Myung Joo; Wu, Shirley; Lee, Sang; Calhoun, Jesse; Hodge, Steven M; Fava, Maurizio; Rosen, Bruce R; Smoller, Jordan W; Gasic, Gregory P; Breiter, Hans C; for the Phenotype Genotype Project on Addiction; Mood Disorders

First page

e13945

Section

Research Article

Publication year

2010

Publication date

Nov 2010

Publisher

Public Library of Science

e-ISSN

19326203

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1371/journal.pone.0013945

ProQuest document ID

1292214895

Microstructural Abnormalities in Subcortical Reward Circuitry of Subjects with Major Depressive Disorder

Jump to:

Full text

Abstract

Details

Suggested sources