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).
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.
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.
Citation: Blood AJ, Iosifescu DV, Makris N, Perlis RH, Kennedy DN, et al. (2010) Microstructural Abnormalities in Subcortical Reward Circuitry of Subjects with Major Depressive Disorder. PLoS ONE 5(11): e13945. doi:10.1371/journal.pone.0013945
Editor: Alessandro Bartolomucci, University of Minnesota, United States of America
Received: May 13, 2010; Accepted: September 16, 2010; Published: November 29, 2010
Copyright: © 2010 Blood et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by a grant to H.C.B. (R01 14118, R21 026104, P20 026002) from NIDA, Bethesda, MD, and a grant (DABK39-03-C-0098; The Phenotype Genotype Project in Addiction and Depression) from the Office of National Drug Control Policy - Counterdrug Technology Assessment Center (ONDCP-CTAC), Washington, D. C. Further support was provided by a grant to A.J.B. (R01 NS052368) from NINDS, by the MGH Department of Radiology to H.C.B., by NCRR (P41RR14075), and by the Mental Illness and Neuroscience Discovery (MIND) Institute. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have included all potential competing interests as follows. This does not alter the authors' adherence to all the PLoS ONE policies on sharing data and materials. Dr. Iosifescu has received research support from Aspect Medical Systems, Forest Laboratories and Janssen Pharmaceutica; he has been a consultant for Forest Laboratories, Gerson Lehrman Group and Pfizer, Inc., and he has been a Speaker for Cephalon, Inc., Eli Lilly and Co., Forest Laboratories and Pfizer, Inc. Dr. Dougherty currently has research support and does consulting with honoraria for Medtronic, Eli Lilly, Brand Ideas, and McNeil, does consulting with honoraria for Reed Elsevier, and has research support from Cyberonics. In the past (more than one year ago), he has had research support and done consulting with honoraria for Northstar Neuroscience, has done consulting with honoraria for Jazz Pharmaceuticals, Wyeth, Bristol Myers Squibb, Trancept Pharmaceuticals, Cyberonics, JK Associates, Inc., American Psychiatric Publishing, Inc., Advocate Health and Hosp. Corp., DHHS/NIH, Leering Swann LLC, Oxford University Press, Professional Practice Group/Psychiatry Syracuse, and Y&R Inc. DBA Sudler & Henessy, and has had research support from Forest and Cephalon. Dr. Rosen is an MRI and PET research sponsor for Siemens. Dr. Fava has received research support from Abbott Laboratories, Alkermes, Aspect Medical Systems, Astra-Zeneca, BioResearch, BrainCells, Inc., Bristol-Myers Squibb Company, Cephalon, Clinical Trial Solutions, LLC, Eli Lilly and Company, Forest Pharmaceuticals Inc., Ganeden, GlaxoSmithKline, J and J Pharmaceuticals, Lichtwer Pharma GmbH, Lorex Pharmaceuticals, NARSAD, NCCAM, NIDA, NIMH, Novartis, Organon Inc., PamLab, LLC, Pfizer Inc, Pharmavite, Roche, Sanofi-Aventis, Shire, Solvay Pharmaceuticals, Inc., Synthelabo, Wyeth-Ayerst Laboratories. He also reported performing advising/consultation with Abbott Laboratories, Amarin, Aspect Medical Systems, Astra-Zeneca, Auspex Pharmaceuticals, Bayer AG, Best Practice Project Management, Inc, BioMarin Pharmaceuticals, Inc.,Biovail Pharmaceuticals, Inc., BrainCells, Inc, Bristol-Myers Squibb Company, Cephalon, Clinical Trials Solutions, LLC, CNS Response, Compellis, Cypress Pharmaceuticals, Dov Pharmaceuticals, Eisai, Inc., Eli Lilly and Company, EPIX Pharmaceuticals, Euthymics Bioscience, Inc., Fabre-Kramer, Pharmaceuticals, Inc., Forest Pharmaceuticals Inc., GlaxoSmithKline, Grunenthal GmBH, Janssen Pharmaceutica, Jazz Pharmaceuticals, J and J Pharmaceuticals, Knoll Pharmaceutical Company, Labopharm, Lorex Pharmaceuticals, Lundbeck, MedAvante Inc., Merck, Methylation Sciences, Neuronetics, Novartis, Nutrition 21, Organon Inc., PamLab, LLC, Pfizer Inc, PharmaStar, Pharmavite, Precision Human Biolaboratory, PsychoGenics, Psylin Neurosciences, Inc., Ridge Diagnostics, Inc., Roche, Sanofi-Aventis, Sepracor, Schering-Plough, Solvay Pharmaceuticals, Inc., Somaxon, Somerset Pharmaceuticals, Synthelabo, Takeda, Tetragenex, TransForm Pharmaceuticals, Inc., Transcept Pharmaceuticals, Vanda Pharmaceuticals Inc, Wyeth-Ayerst Laboratories. He reported performing speaking/publishing for Adamed, Co., Advanced Meeting Partners, American Psychiatric Association, American Society of Clinical Psychopharmacology, Astra-Zeneca, Belvoir, Boehringer-Ingelheim, Bristol-Myers Squibb Company, Cephalon, Eli Lilly and Company, Forest Pharmaceuticals Inc., GlaxoSmithKline, Imedex, Novartis, Organon Inc., Pfizer Inc, PharmaStar, MGH Psychiatry Academy/Primedia, MGH Psychiatry Academy/Reed-Elsevier, UBC, Wyeth-Ayerst Laboratories. He reported having equity holdings in Compellis, receiveing royalty/patent or other income from SPCD and for a combination of azapirones and bupropion in MDD, and copyright royalties for the MGH CPFQ, SFI, ATRQ, DESS, and SAFER. Drs. Blood, Makris, Perlis, Kennedy, Smoller, Gasic, and Breiter have no relevant financial interest in this manuscript. Byoung Woo Kim, Myung Joo Lee, Shirley Wu, Sang Lee, Jesse Calhoun, and Steven Hodge also have no relevant financial interest in this manuscript.
It has been proposed that major depressive disorder (MDD) may result from dysfunction of brain reward/aversion circuitry, , , , , , . Hypothesized in 1999 as a general schema for processing both positive and negative features of potential “goal-objects” or states , 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 , , , , . 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 , , , ,  and synopsized , , . 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 .
The possibility that abnormalities in this reward/aversion circuitry underlie many psychiatric conditions, including MDD, was further schematized by multiple investigators , , , . 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 , ; these subcortical regions are at the core of the animal literature that first identified reward circuitry , , , , . In humans with MDD, there has recently been detection of functional abnormalities in these subcortical regions , , , in addition to well-established abnormalities in target regions receiving subcortical projections , , , , , , , , , , , , , , . 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 , .
We designed the study to take into consideration potential heterogeneity within the MDD cohort , , , , , , , , , , , given that this is a potential confound in group designs . 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.
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 ] 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, ); and (8) the STAI (Spielberger State Trait Anxiety Inventory).
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  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 .
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. . 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 . 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 , . In PD patients, FA was decreased in the SN, which is consistent with neuronal loss in this region , . 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. 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  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.
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 , . 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 ,  and cortical thickness analysis. 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.
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  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 , ; 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.
1. Image Contrast Analysis of MDD Versus Control Subjects
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.
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.
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.
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.
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).
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.
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.
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.
b. Voxel-based contrasts of FA maps for the two MDD subgroups.
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.
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.
(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].
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.
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 , 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. , ], 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 , , 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 ; 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 . 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  along with imaging abnormalities in this and other prefrontal gray matter regions in MDD, , , , 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 , , 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 . 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 , , while the SN processes both reward/aversion and motor information . The SN is thought to participate in modulating reward prediction and expectancy , as well as prediction of aversive and negatively valenced stimuli . In addition, medial portions of the SN project to and receive projections back from the ventral striatum , , which processes reward/aversion information , while more lateral portions project to the dorsal striatum, which processes motor information . 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 . Other efferents of the SN include GABAergic nigro-collicular pathways which have been shown to mediate fear/defense reactions .
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 . 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 , . 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 , our findings are also consistent with reports of altered nucleus accumbens function in MDD , , supporting earlier hypotheses of such an effect .
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 . Although Milak and colleagues  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 , , . 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 , and connects with the idea that abnormal anxiety may reflect altered expectancy , . Elevated trait anxiety in this subgroup is also consistent with findings by Berton & colleagues 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 , and support investigation of interventions based thereon .
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.
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 , , 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 , 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 ,  and more recently human literature in movement disorders  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.
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 .
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).
(0.05 MB DOC)
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.
(0.91 MB TIF)
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.
(0.52 MB TIF)
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.
(0.65 MB TIF)
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.
(0.70 MB TIF)
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.
(0.46 MB TIF)
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.
(0.51 MB TIF)
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).
(0.74 MB TIF)
Permutation tests. Reports the methods and results of the permutation analysis conducted as a complement to the voxel-wise contrasts of FA maps.
(0.12 MB DOC)
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.
(0.03 MB DOC)
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.
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.
- 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.
- 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.
- 3. Forbes EE (2009) Where's the fun in that? Broadening the focus on reward function in depression. Biol Psychiatry 66: 199–200.
- 4. McCabe C, Cowen PJ, Harmer CJ (2009) Neural representation of reward in recovered depressed patients. Psychopharmacology (Berl) 205: 667–677.
- 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.
- 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.
- 7. Nestler EJ, Carlezon WA Jr (2006) The mesolimbic dopamine reward circuit in depression. Biol Psychiatry 59: 1151–1159.
- 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.
- 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.
- 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.
- 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.
- 12. Becerra L, Breiter HC, Wise R, Gonzalez RG, Borsook D (2001) Reward circuitry activation by noxious thermal stimuli. Neuron 32: 927–946.
- 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.
- 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.
- 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.
- 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.
- 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.
- 18. Matsumoto M, Hikosaka O (2009) Two types of dopamine neuron distinctly convey positive and negative motivational signals. Nature 459: 837–841.
- 19. Leknes S, Tracey I (2008) A common neurobiology for pain and pleasure. Nat Rev Neurosci 9: 314–320.
- 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.
- 21. Carlezon WA Jr, Thomas MJ (2009) Biological substrates of reward and aversion: a nucleus accumbens activity hypothesis. Neuropharmacology 56: Suppl 1122–132.
- 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.
- 23. Bogdan R, Pizzagalli DA (2006) Acute stress reduces reward responsiveness: implications for depression. Biol Psychiatry 60: 1147–1154.
- 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.
- 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.
- 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.
- 27. Wise RA, Bozarth MA (1985) Brain mechanisms of drug reward and euphoria. Psychiatr Med 3: 445–460.
- 28. Wise RA (1978) Catecholamine theories of reward: a critical review. Brain Res 152: 215–247.
- 29. Manji HK, Drevets WC, Charney DS (2001) The cellular neurobiology of depression. Nat Med 7: 541–547.
- 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.
- 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.
- 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.
- 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.
- 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.
- 35. Drevets WC, Ongur D, Price JL (1998) Reduced glucose metabolism in the subgenual prefrontal cortex in unipolar depression. Mol Psychiatry 3: 190–191.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 50. Sobin C, Sackeim HA (1997) Psychomotor symptoms of depression. Am J Psychiatry 154: 4–17.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 60. Hamilton M (1967) Development of a rating scale for primary depressive illness. Br J Soc Clin Psychol 6: 278–296.
- 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.
- 62. Deichmann R, Gottfried JA, Hutton C, Turner R (2003) Optimized EPI for fMRI studies of the orbitofrontal cortex. Neuroimage 19: 430–441.
- 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.
- 64. Pierpaoli C, Basser PJ (1996) Toward a quantitative assessment of diffusion anisotropy. Magn Reson Med 36: 893–906.
- 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.
- 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.
- 67. Jenkinson M, Smith S (2001) A global optimisation method for robust affine registration of brain images. Med Image Anal 5: 143–156.
- 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.
- 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.
- 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.
- 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.
- 72. Wainer H (1978) Gapping. Psychometrika 43: 203–212.
- 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.
- 74. Beaulieu C (2002) The basis of anisotropic water diffusion in the nervous system - a technical review. NMR Biomed 15: 435–455.
- 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.
- 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.
- 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.
- 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.
- 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.
- 80. Bains JS, Oliet SH (2007) Glia: 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.
- 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.
- 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.
- 84. Schultz W, Dayan P, Montague PR (1997) A neural substrate of prediction and reward. Science 275: 1593–1599.
- 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.
- 86. Wiedemann G, Pauli P, Dengler W (2001) A priori expectancy bias in patients with panic disorder. J Anxiety Disord 15: 401–412.
- 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.
- 88. Chan CK, Lovibond PF (1996) Expectancy bias in trait anxiety. J Abnorm Psychol 105: 637–647.
- 89. Wise RA (2009) Roles for nigrostriatal—not just mesocorticolimbic—dopamine in reward and addiction. Trends Neurosci 32: 517–524.
- 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.
- 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.
- 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.
- 93. Malberg JE, Eisch AJ, Nestler EJ, Duman RS (2000) Chronic antidepressant treatment increases neurogenesis in adult rat hippocampus. J Neurosci 20: 9104–9110.
- 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.
- 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.