A COMPREHENSIVE EXAMINATION OF WHITE MATTER TRACTS AND CONNECTOMETRY IN MAJOR DEPRESSIVE DISORDER.

BACKGROUND: Major depressive disorder (MDD) is a debilitating disorder characterized by widespread brain abnormalities. The literature is mixed as to whether or not white matter abnormalities are associated with MDD. This study sought to examine fractional anisotropy (FA) in white matter tracts in individuals with MDD using diffusion tensor imaging (DTI).
METHODS: 139 participants with MDD and 39 healthy controls (HC) in a multisite study were included. DTI scans were acquired in 64 directions and FA was determined in the brain using four methods: region of interest (ROI), tract-based spatial statistics (TBSS), and diffusion tractography. Diffusion connectometry was used to identify white matter pathways associated with MDD.
RESULTS: There were no significant differences when comparing FA in MDD and HC groups using any method. In the MDD group, there was a significant relationship between depression severity and FA in the right medial orbitofrontal cortex, and between age of onset of MDD and FA in the right caudal anterior cingulate cortex using the ROI method. There was a significant relationship between age of onset and connectivity in the thalamocortical radiation, inferior longitudinal fasciculus, and cerebellar tracts using diffusion connectometry.
CONCLUSIONS: The lack of group differences in FA and connectometry analysis may result from the clinically heterogenous nature of MDD. However, the relationship between FA and depression severity may suggest a state biomarker of depression that should be investigated as a potential indicator of response. Age of onset may also be a significant clinical feature to pursue when studying white matter tracts.

INTRODUCTION

Major depressive disorder (MDD) is a debilitating psychiatric disorder that is highly prevalent in the US population. MDD is characterized by a variable combination of affective and cognitive symptoms, and not surprisingly is associated with widespread structural and functional brain differences. In particular, there has been a great deal of attention focused on circuits connecting frontal and subcortical brain regions.1, 2, 3 Evidence suggests that MDD is associated with structural deficits in midline frontal brain regions, particularly in the anterior cingulate cortex (ACC), medial orbitofrontal cortex (MOFC), and dorsomedial prefrontal cortex (DMPFC).2, 4 Abnormal circuits connecting these frontal regions with subcortical regions are hypothesized to underlie the neurobiology of depression,1 leading to a “disconnection syndrome” resulting from white matter deficits.5, 6 In other words, a “disconnection syndrome” occurs when frontal and subcortical brain regions cannot effectively communicate, which is hypothesized to lead to the manifestation of depressive symptoms. One way to examine this phenomenon is to study white matter tracts, which are essential for communication between different brain regions.

Diffusion tensor imaging (DTI) can be used to characterize white matter tracts in the human brain in vivo. Fractional anisotropy (FA) is a DTI‐derived measure reflecting the directionality of water diffusion, which is used as an indication of the location and strength of white matter tracts. FA in certain regions has been reported to be lower in MDD suggesting white matter fiber disruption or disorganization. For example, low FA in MDD has been reported in white matter adjacent to frontal,7, 8, 9 temporal,8 parietal,8, 9 and cingulate cortices.10 Low FA has also been reported in white matter tracts, including the inferior longitudinal fasciculus,11, 12 superior longitudinal fasciculus,9, 13 internal capsule,10, 12, 14, 15 external capsule,14, 15 and corpus callosum.13, 14, 15 We recently reported extensive white matter deficits in frontal medial brain regions in MDD compared to healthy controls (HCs),16 mirroring findings reported by a recent meta‐analysis.17 However, not all studies find group differences in FA between MDD and HCs18, 19 and one study reported low FA only in those with melancholic MDD.20

Using diffusion tractography methods, it is possible to examine FA in individual tracts to determine which tracts may be affected in MDD, rather than discrete brain regions. Song and colleagues21 found lower FA in the right solitary tract of MDD subjects than HC subjects. However, another study reported lower FA in MDD in only one of seven tracts investigated (the parahippocampal cingulum) and only in male subjects.22 Most studies do not find a relationship between FA and clinical characteristics.3, 7, 8, 11, 14, 15, 18, 21 However, two studies showed a negative correlation between depression severity and FA in the left anterior limb of the internal capsule.10, 12

There are several limitations in the DTI literature that need to be addressed. First, there is only one group that has reported DTI data with a large MDD sample (N = 95).3, 21 Other studies report on samples ranging from 148 to 45 MDD subjects.12, 13 Second, about half of the studies in the literature used 13 or fewer directions,7, 19, 23 and only one group used 61 directions.13 FA can vary depending on the number of directions,24 and a minimum of 30 unique directions has been shown to provide adequate resolution for estimating FA.25 Third, some studies report data that were not corrected for multiple comparisons7, 8, 19, 26 or report a small cluster size.11 Finally, most studies use only one method to examine FA, either region of interest (ROI),7 voxel based analysis (VBA),8, 27 tract‐based spatial statistics (TBSS),13, 19 or tractography.3, 28

To address shortcomings in previous studies, we examined FA in the largest MDD sample to date using 64‐direction DTI, using multiple approaches: ROI, TBSS, and probabilistic tractography. Additionally, differences in white matter were explored using diffusion connectometry,29 which aimed to map the trajectories of affected tracts.30, 31 Each technique brings a unique contribution to this investigation. Namely, the ROI approach is focused on a small number of a priori brain regions based on the literature. A further advantage of the ROI approach is that the data are obtained in DTI space, so there is no warping of the brain to fit a template. TBSS is a voxel‐based approach to examine white matter tracts. Although the data are derived from a standardized template, this approach allows for a broad investigation of the entire brain. Probabilistic tractography is another focused technique that allows investigators to identify specific tracts by identifying seed ROIs and estimating how many tracts extend from those ROIs. The benefit of this approach is that one can examine white matter indices along a specific tract, rather than discrete parts of the brain that are assessed using ROI and TBSS approaches. Finally, diffusion connectometry differs from probabilistic tractography in that probabilistic tractography defines connectivity by the number of “tracks” or “streamlines,” whereas connectometry uses density of diffusing spins. Conceptually, probabilistic tractography aims to find a difference in tracks, whereas connectometry tracks the difference in voxels that have substantial correlation with the study variable, aiming to identify the entire affected section.

By addressing all the limitations in the literature, we expect the results from this study to provide more definitive results regarding the possibility of white matter deficits in MDD. We hypothesized that MDD controls would have lower FA in white matter tracts in the frontal cortex than HCs. Although there is no supportive literature to develop a specific hypothesis about the connectometry analysis, we expect to identify white matter deficits in the frontal cortex in MDD compared to HCs which is consistent with the FA literature.

MATERIALS AND METHODS

SUBJECTS

Participants were enrolled in the EMBARC (Establishing Moderators and Biosignatures of Antidepressant Response for Clinical Care for Depression, NIMH 1U01 MH092250, http://embarc.utsouthwestern.edu/) project at four sites. These sites are the University of Texas Southwestern Medical Center (TX), University of Michigan (UM), Massachusetts General Hospital (MG), and Columbia University Medical Center (CU). The EMBARC study was designed to measure treatment outcomes. However, the current study is a preliminary analysis aimed at comparing baseline neuroimaging data in depressed and HCs. The Institutional Review Board for all four sites approved the protocol, and subjects gave written informed consent.

One hundred thirty‐nine subjects who met Diagnostic and Statistical Manual of Mental Disorders (DSM‐IV‐TR)32 criteria for a current major depressive episode (MDE) in context of MDD and 39 HCs were included. Overall, there were 166 MDD subjects enrolled in the EMBARC study at the time of this analysis. Twenty‐seven MDD subjects were not included in the analysis for the following reasons: they did not complete baseline imaging (N = 10), had poor quality images due to significant head motion (N = 4) and poor segmentation (N = 11). Therefore, there were 139 MDD subjects remaining for analysis. As for the HCs, there were 40 HC subjects enrolled in the EMBARC study at that time of this analysis. One HC was not included due to significant head motion.

MDD and HCs were recruited by flyers, advertisements, and posting on websites. Participants who were potentially eligible were scheduled for an evaluation visit where study eligibility was assessed through psychiatric and medical history, chart review, clinical interview, physical examination, routine blood tests, pregnancy test, and urine toxicology. All participants were between the ages of 18 and 65. MDD participants had to meet DSM‐IV‐TR criteria for a current MDE (early onset, before 30 years) and have a minimum score of 14 on the Quick Inventory of Depressive Symptoms (QIDS‐SR16).33 MDD participants were excluded if they had a lifetime history of bipolar disorder or psychosis, current drug or alcohol abuse/dependence, or a current primary diagnosis of obsessive–compulsive disorder. Patients on psychoactive medication underwent a 3‐week medication washout prior to neuroimaging. The following current comorbid diagnoses were reported in the MDD group: panic disorder (N = 15), obsessive–compulsive disorder (N = 4), social phobia (N = 23), PTSD (N = 11), generalized anxiety disorder (N = 19), anxiety disorder, not otherwise specified (N = 5), bulimia nervosa (N = 2), or other diagnosis (N = 2). HCs were included if they did not have a lifetime history of a mood, anxiety, or psychotic disorder.

CLINICAL MEASURES

Diagnoses were based on the Structured Clinical Interview for DSM‐IV Patient Version (SCID‐IV‐I/P).34 The Hamilton Depression Rating Scale (HAM‐D)35 and QIDS‐SR1633 assessed clinician‐ and self‐rated depression severity, respectively. Scores from clinical measures administered at the evaluation visit are reported.

IMAGE ACQUISITION

All participants underwent a magnetic resonance imaging (MRI) scan at one of the participating sites. Each site had a different 3 T MRI scanner. Two sites operated Philips scanners (TX and UM), one site operated a GE scanner (CU), and one site operated a Siemens scanner (MG). Structural MRIs were acquired with following parameters: TR (repetition time): 5.9–8.2 ms, TE (echo time): 2.4–3.7 ms, Flip angle: 9°–12°, slice thickness: 1 mm, FOV (field of view): 256 × 256 mm, voxel dimensions: 1 × 1 × 1 mm, acquisition matrix: 240 × 240, 256 × 243, or 256 × 256, acceleration factor: 2, and 174–178 sagittal slices. Diffusion images were acquired using a single‐shot echo planar imaging (EPI) sequence. Scan parameters were as follows: TR = 8,310–9,500 ms, TE = 95–96.3 ms, Flip angle = 90 degrees, slice thickness = 2.5 mm, FOV = 240 × 240 mm, voxel dimensions 2.5 × 2.5 × 2.5 mm or 1.9 × 1.9 × 2.5 mm, acquisition matrix = 96 × 96, b value = 1,000 s/mm2, and 64 collinear directions with 1 or 5 nonweighted images. DTI scan time was approximately 10 min.

IMAGE PROCESSING

Diffusion Tensor Imaging (DTI)

All T1 and DTI images were transferred to and processed at Stony Brook University. Each DTI image was run through a series of quality assurance tests for common artifacts, including ghost, ring, slice‐wise intensity, venetian blind, and gradient‐wise motion artifacts.36 Diffusion images were then corrected for distortion induced by gradient coils and simple head motion using the eddy current correction routine within FSL (FMRIB Software Library, http://www.fmrib.ox.ac.uk/fsl/). FSL's Brain Extraction Tool (BET) removed nonbrain tissue from the image. Following this, Camino [http://web4.cs.ucl.ac.uk/research/medic/camino/pmwiki/pmwiki.php] was used to estimate FA computing the least squares fit diffusion tensor with nonlinear optimization using a Levenburg–Marquardt algorithm, constrained to be positive by fitting its Cholesky decomposition.37, 38

ROI Determinations

Anatomical T1 images were cropped and corrected for image nonuniformity using the Bias Field Corrector in Advanced Normalization Tools (ANTS; http://www.picsl.upenn.edu/ANTS/). ROIs were created using the cortical parcellation tool based on the Desikan–Killiany atlas.39 Each white matter voxel was labeled on the basis of the nearest cortical voxel.40 Then, the DTI image was coregistered to the cropped T1 image using ANTS, and the inverse transformation was applied to the cortical map in order to place ROIs into DTI space for analysis. Finally, mean FA values in white matter were computed for each ROI. In order to reduce the number of comparisons, the following bilateral ROIs in the frontal cortex were chosen a priori based on reported white matter deficits in frontal cortical regions in the MDD literature: MOFC, superior frontal cortex (SFC), rostral anterior cingulate cortex (rACC), and caudal anterior cingulate cortex (cACC).2, 6 Figure 1 shows the location of the right cACC and the right MOFC ROIs. We renamed the SFC label “dorsomedial prefrontal cortex” or “DMPFC” to be consistent with the literature and provide a more descriptive name for the region.

Figure 1

Sagittal, coronal, and axial view of the right caudal anterior cingulate cortex region of interest (ACC ROI) (top) and the right medial orbitofrontal cortex (MOFC) ROI (bottom) in a representative subject.

Tract‐based Spatial Statistics (TBSS) Analysis

TBSS analysis41 of the FA diffusion maps included several steps. Individual FA maps were nonlinearly aligned to a common FMRIB58 FA template. Then the aligned FA maps were averaged to obtain a study mean FA map, and a white matter skeleton was computed on the mean FA map. The white matter skeleton was then thresholded (FA ≥ .2) to maintain only white matter (WM) voxels. Each individual FA map was projected onto the common FA skeleton to obtain the individual FA skeleton, on which a voxel‐wise analysis was performed to examine differences between the two populations.

Tractography and Weighted FA

Probabilistic tractography was performed using FMRIB's Diffusion Toolbox (FDT).42 This algorithm computes probabilistic streamlines through each voxel by repetitive sampling from the principal diffusion directions. The outcome of this method is the probability of connections from the seed (MOFC) to all targets in the brain. The MOFC was chosen based on the finding that FA in this region was significantly correlated with depression severity. The algorithm was run with a tract curvature threshold of 0.2 mm, maximum number of steps per sample equal to 2,000, length of each step equal to 0.5 mm, and 5,000 samples. In order to calculate the (weighted) average FA within the defined tracts, the voxel‐based FA measures were multiplied by the probability of connection at each voxel and divided by the sum of probabilities. The two outcomes of interest were weighted average FA and number of tracts. The tractography algorithm failed in 16 MDD subjects, therefore data on 123 MDD subjects are reported. Additionally, two MDD subjects were removed for the “total number of tracts” analysis due to extremely high values (>7 million tracts).

Diffusion MRI Connectometry

We conducted diffusion MRI connectometry29 to further explore group difference. The diffusion data were reconstructed using q‐space diffeomorphic reconstruction,43 with a diffusion sampling length ratio of 1.25. The output resolution was 2 mm. A multiple regression model including age, site of data collection, and sex was used to calculate the T‐score of the difference between control subjects and patients. A deterministic fiber tracking algorithm44 was used to connect local fiber directions with T‐score greater than 2.0, and tracts with connected length greater than 40 mm were collected. This cutoff was determined by the lower bound of the uncinated fasciculus length,45 which is used to estimate the minimum length of a major tract. To estimate the false discovery rate, a total of 1,000 randomized permutations were applied to the group label in order to obtain the null distribution of the tract length. In the MDD group, clinical measures (i.e., depression severity and age of onset of depression) were used in a multiple regression model to determine their relationship to connectometry results. The regression model included site of data collection, age, and sex. The analysis was conducted using publicly available software DSI Studio (http://dsi‐studio.labsolver.org).

STATISTICAL ANALYSIS

Data were statistically evaluated using IBM SPSS Statistics (SPSS Inc., Chicago, Illinois, USA, Version 22.0). Comparisons of demographic and clinical variables were performed with a univariate analysis of covariance (ANCOVA) for continuous variables and chi‐square analyses for categorical variables with group (HC vs. MDD) as the between‐subjects factor with site as a covariate. For the ROI analyses, ANCOVAs were performed with site, age, and sex included as covariates. For the TBSS analyses, we used the general linear model (GLM) method with the permutation‐based nonparametric testing included in the FSL tools. Significance threshold was set to P < .05 corrected for family‐wise errors (FWE) using the threshold‐free cluster enhancement (TFCE) method in FSL's randomized tool. Using TFCE, FWEs are corrected using the null distribution of the maximum test statistic across the image.46 The group comparison was adjusted for (mean centered) age, site and sex. Only clusters with a minimum of 50 voxels were considered significant. For the tractography analyses, an ANCOVA was used to compare MDD and HC groups with site, age, and sex as covariates. For number of tracts, MOFC seed size was included as a covariate. Partial correlations between clinical and diffusion measures were performed in the MDD group using Spearman's Rho for ordinal variables (HAM‐D 17 and QIDS), variables that were not normally distributed (duration of MDE), and variables that were converted into categorical categories (number of MDEs; 1 = 1–2 episodes, 2 = 3–4 episodes, 3 = 5–8 episodes, 4 = 9–20 episodes, and 5 = 24+ episodes. Each category included approximately 15–20% of the sample). Pearson correlation coefficient was used for age of onset correlations. All correlations included age, site, and sex as covariates. Bonferroni correction was used to account for multiple correlations (.05/50 = .001).

RESULTS

DEMOGRAPHICS AND CLINICAL MEASURES

Demographic and clinical data are presented in Table 1. There were no significant differences in age (F(1,175) = 0.12, P = .73) or years of education (F(1,170) = 0.46, P = .50) between HC and MDD groups There were no significant differences in sex (χ2(1) = 0.006, P = .94), race (χ2(4) = 1.25, P = .87), ethnicity (χ2(1) = 1.71, P = .19), or handedness (χ2(2) = 0.24, P = .88) between the HC and MDD groups. Additionally, there were no group differences when comparing sex, race, ethnicity, or handedness by site (all P‐values >.05), with the exception that TX has significantly younger HCs compared to the MDD group (F(1,64) = 7.75, P = .007). As expected, the MDD group scored higher than the HC group on the HAM‐D (F(1,169) = 551.05, P < .001).

Table 1

Demographic and clinical variables

	HC (N = 39)	MDD (N = 139)
Demographic variables
Sex, male; N (%)	14 (36%)	49 (35%)
Age (years), mean	37.0 ± 14.6	37.2 ± 13.1
Range	18–65	18–65
Race, N (%); Caucasian	13 (33%)	43 (31%)
Ethnicity, N (%); Hispanic	3 (8%)	22 (16%)
Handedness, N (%); Right	32 (82%)	118 (85%)
Education (years), mean	15.2 ± 2.3	14.8 ± 2.5
Range	10–20	8–21
Clinical variables
Age of onset of MDE, mean	—	15.8 ± 6.2
Range		2–30
Number of MDEs, median	—	6
Range		1–120
Length of current episode, mean	—	37.6 ± 63.0
Range (months)		1–400
HAM‐D 17‐item*, mean	0.7 ± 0.9	18.7 ± 4.5
Range	0–3	10–32
QIDS‐SR16 *, mean	1.4 ± 1.3	18.3 ± 2.8
Range	0–4	14–27

Data are presented as mean ± standard deviation unless otherwise specified.

HC, healthy control; MDD, major depressive disorder; MDE, major depressive episode; HAM‐D, Hamilton Rating Scale for Depression; QIDS‐SR16, quick inventory of depressive symptoms.

*P < .001 versus HC.

ROI AND TBSS ANALYSES

Table 2 presents mean FA in each of the bilateral ROIs investigated (cACC, rACC, MOFC, and DMPFC). There were no significant differences in mean FA in ROIs when comparing the HC and MDD groups. There were no clusters identified by the TBSS analysis that differentiated the two groups.

Table 2

Fractional anisotropy (FA) in white matter parcellated regions of interest (ROIs) using the Desikan‐Killiany atlas in FreeSurfer

ROI	R/L	HC (N = 39)	MDD (N = 139)	F‐test
Caudal ACC	R	.62 ± .04	.62 ± .05	F(1,173) = 0.20, P =.66
	L	.61 ± .04	.61 ± .05	F(1,173) = 0.02, P =.88
Rostral ACC	R	.58 ± .05	.59 ± .05	F(1,173) = 0.21, P =.65
	L	.55 ± .04	.54 ±. 04	F(1,173) = 0.50, P =.48
Medial orbital frontal	R	.43 ± .03	.43 ± .03	F(1,173) = 0.79, P =.38
	L	.41 ± .04	.40 ± .03	F(1,173) = 1.52, P =.22
Dorsomedial prefrontal	R	.48 ± .03	.49 ± .03	F(1,173) = 0.67, P = .41
	L	.49 ± .03	.49 ± .03	F(1,173) = 0.39, P = .54

F‐test includes age, site and sex as covariates. Data are presented as mean ± standard deviation.

ACC, anterior cingulate cortex; HC, healthy control; L, left; MDD, major depressive disorder; R, right.

TRACTOGRAPHY AND WEIGHTED FA

Probabalistic tractography was used to examine white matter tracts extending from the right MOFC in light of the correlations reported below. Table 3 presents the means for the two outcome measures from the tractography: weighted FA and total number of tracts. Seed size was used as a covariate for total number of tracts. There were no significant differences in these outcome measures when comparing the HC and MDD groups.

Table 3

Outcome measures from probabilistic tractography using the right MOFC as the seed and all voxels in the brain as the target

	HC (N = 39)	MDD (N = 123)*	F‐test
Weighted FA	.32 ± .03	.32 ± .03	F(1,157) = .32, P = .57
Total number of tracts	2,041,410 ± 575,536	2,030,248 ± 523,381	F(1,155) = .12, P = .73

F‐test includes age, site, and sex as covariates. For total number of tracts, seed size is also included as a covariate. Data are presented as mean ± standard deviation.

FA, fractional anisotropy; HC, healthy control; MDD, major depressive disorder.

*Two subjects were removed for the “total number of tracts” analysis due to extremely high values (>7 million tracts).

DIFFUSION MRI CONNECTOMETRY

We used connectometry to identify white matter pathways associated with MDD.29 The density of diffusing spins at various orientations were calculated and regressed against MDD and HC groups, with age, sex, and site included in the regression model. The results did not show any white matter pathways that were significantly associated with MDD.

CLINICAL CORRELATIONS

In order to test the relationship between FA and clinical severity, correlations were performed between DTI and clinical measures in the MDD group (Table 4). There was a significant correlation between FA in the right MOFC and HAM‐D scores (P < .001) and FA in the right cACC and age of onset (P = .003; Fig. 2). There were no significant correlations between tractography outcome measures and clinical measures at the P ≤ .001 level. Diffusion MRI connectometry was also used to identify white matter pathways associated with depression severity and age of onset in the MDD group. The analysis showed a significant negative correlation (FDR = 0.02) between age of onset of depression and connectivity in four fiber pathways: (1) left sensory pathway through the posterior limb of the internal capsule, (2) right thalamocortical radiation, (3) left inferior longitudinal fasciculus, and (4) bilateral cerebellar tracts (Fig. 3). The connectometry analysis showed no significant correlation between the MOFC and HAM‐D scores (FDR > 0.05).

Table 4

Partial correlation coefficients for the association between FA in white matter parcellated ROIs, tractography outcome measures, and clinical measures in the MDD group

	R/L	HAM‐D	QIDS‐SR16	Age of onset	Length of MDE (months)	Number of MDEs
ROI analysis
Caudal ACC	R	−.10	.004	.25*	−.003	−.07
	L	.04	.00	.10	.01	.03
Rostral ACC	R	−.17	−.04	.07	.03	−.09
	L	.04	.00	.05	.01	.03
Medial orbital frontal	R	−.29**	−.17	.23	−.06	.02
	L	−.20	−.10	.15	−.12	.03
Dorsomedial prefrontal	R	−.24	.02	.21	−.08	.09
	L	−.18	.01	.13	−.05	.14
Tractography
Weighted FA in tract extending from the medial orbital frontal seed	R	−.18	−.09	.02	.06	−.03
Number of tracts extending from the medial orbital frontal seed	R	−.09	.23	.03	.03	.01

Age, site, and sex are included as covariates.

ACC, anterior cingulate cortex; FA, fractional anisotropy; HAM‐D, Hamilton Rating Scale for Depression; HP, healthy participant; L, left; MDD, major depressive disorder; MDE, major depressive episode; QIDS‐SR16, quick inventory of depressive symptoms; ROI, regions of interest; R, right.

*P = .003.

**P < .001 indicates a significant correlation.

Figure 2

Top: Scatterplot depicting the relationship between fractional anisotropy (FA) in the right medial orbitofrontal cortex (MOFC) and Hamilton Depression Rating Scale scores (HAM‐D) in the major depressive disorder (MDD) group. Bottom: Scatterplot depicting the relationship between FA in the right caudal anterior cingulate cortex (cACC) and age of onset of depression in the MDD group.

Figure 3

Connectometry results showing significant negative correlations between age of onset of depression in the major depressive disorder (MDD) group and white matter density in the following tracts: (1) left sensory pathway through the posterior limb of the internal capsule (red), (2) right thalamocortical radiation (orange), (3) left inferior longitudinal fasciculus (yellow), and (4) bilateral cerebellar tracts (blue and green). Diffusion tensor imaging (DTI) image slices are presented in radiologic convention, and the tract images are presented in neurological convention.

DISCUSSION

Prevailing theories of MDD include impaired connections in the brain linking cortical and subcortical regions.1, 2, 3 The current DTI study examined white matter tracts in MDD in one of the largest samples published to date, using a 64‐direction protocol. Using multiple methods to examine FA we found no difference in FA comparing individuals with MDD and HCs. This is at odds with a previous paper where we reported decreased FA in MDD in several frontal cortex regions.16 One main difference between the two studies is that the HC group was significantly younger than the MDD group in our previous paper. FA decreases in adulthood47 making this a potentially confounding variable. Although we used age as a covariate in all of our analyses, it is possible that variability due to age could not be completely regressed out.48 The literature on DTI in MDD has been mixed, both the findings and the methods used to examine white matter, making it difficult to draw any concrete conclusions. The current study had groups well matched on all demographic variables, had a larger sample, and had a higher resolution DTI scan than our previous study, strengthening our confidence in our findings.

Of significant interest, we found a relationship between FA in the right MOFC and depression severity, such that low FA was related to high depression scores. The MOFC is thought to be central to reward processing in the brain.49, 50 It is also uniquely situated to integrate sensory and emotional information through its connections with visual, auditory,51 and limbic structures.1 There have been only two other DTI papers that reported a relationship between FA and depression severity in adult‐onset MDD.10, 12 Both groups found that FA in the anterior limb of the internal capsule was negatively correlated with depression severity. One study of late‐life depression found that FA in the MOFC was negatively correlated with depression severity,52 and another found a positive correlation between MOFC white matter lesion density and depression severity,53 supporting our finding. DTI studies that did not find such a relationship between clinical measures and FA had a small sample size, whereas the studies reporting a relationship had 25 MDD10 and 45 MDD subjects.12 When analyzing data from each of the four EMBARC sites separately, only two of the sites showed a significant correlation between the right MOFC and depression severity. These two sites had the largest samples. It is possible that studies with a small sample size limit the range of depression severity scores, making it unlikely to find such a relationship. Although two DTI studies with large samples did not find a significant correlation between FA and depression severity,3, 21 these studies were focused on FA in specific locations21 or broad network based measures of FA.3

Age of depression onset was also significantly related to FA in the right cACC using the ROI method, and to white matter density in the thalamocortical radiation (connecting the thalamus and caudal ACC, see Fig. 2) using diffusion connectometry. This is the first study to report a significant relationship between age of onset of depression and DTI measures; other groups do not find a significant relationship.11, 18, 21 A recent study showed that individuals with pediatric‐onset MDD had a smaller genu of the corpus callosum compared to HCs, whereas corpus callosum volume was not different when comparing adult‐onset MDD and HCs.54 It is possible that suffering from depression while white matter growth is still in progress (i.e., throughout childhood and adolescence) may have a deleterious effect on continued development. There have been a number of studies that support the notion that early life stress adversely affects brain development,55, 56 which is consistent with this finding. It is important to note that age of depression onset was reported as early as 2 years of age in some MDD subjects. Limitations associated with retrospective recall is a significant factor to consider, although we found comparable results when limiting age of onset to a minimum of 10 years (r = .24, P = .017). Additionally, the significance of the relationship between FA in the right cACC using the ROI method (P = .003) was just below our strict significance cut‐off of P < .001. This finding should be interpreted with caution.

The other three tracts identified using the diffusion connectometry method are important for sensory processing. The internal capsule facilitates the communication of sensory information to the cortex, the inferior longitudinal fasciculus connects the orbitofrontal and occipital cortices, and cerebellar tracts provide proprioceptive information about the limbs to the cerebellum. There have been some studies that report low FA in these tracts in MDD,10, 11, 12, 14, 15 no previous study reported a significant correlation between these tracts and age of onset of MDD. Future studies are needed to confirm this relationship and to expand on the clinical and behavioral correlates underlying this relationship.

There is increasing support for evaluating dimensional variables, rather than categorically grouping individuals based on heterogeneous symptoms, in order to identify biomarkers associated with psychiatric symptoms.57 Therefore, the findings from the current study lend support for a relationship between depression severity, age of onset of depression, and white matter integrity. The main goal of the current study is to find moderators of treatment outcome, which doesn't require a biomarker to distinguish controls and MDD subjects. In fact, it is likely that it will be difficult to find a strong biomarker linked to depression as an overarching construct that describes a clinically heterogeneous disorder. Since an important goal of this research is to find endophenotypes or biomarkers that are useful for diagnostic and prognostic purposes, as well as potential biomarkers for personalized treatment, we are hopeful that the current findings point in an interesting direction for further examination.

CONCLUSIONS

Overall, the findings from the current study should be considered a benchmark for future work, taking great care in the recruitment of a large number of participants who are well matched, and the acquisition of high‐resolution data. Future studies should investigate dimensional characteristics within depression samples that may be linked to DTI measures in specific white matter tracts that could provide a connection between behavioral and biological measures. Parsing the clinically heterogeneous groups of depressed individuals, such as isolating those with early‐onset depression, will be of significant value moving forward in order for us to understand how depression may be impacting the development of white matter connections in the brain.

Conflict of interest

There are no conflicts of interest or financial disclosures related to this work.

Disclaimer: Supplementary materials have been peer‐reviewed but not copyedited.

Table S1. Demographic and clinical variables

Table S2. FA in white matter parcellated ROIs

Table S3. Outcome measures from probabilistic tractography

Table S4. Correlation coefficient for relationship between FA in ROIs (top) and tractography (bottom) and clinical measures by site.

Click here for additional data file.