INTRODUCTION

Using magnetic resonance imaging (MRI) techniques, researchers studying major depressive disorder (MDD) have found a number of brain regions to be implicated in mood regulation, including prefrontal areas, limbic structures, and subcortical gray matter regions (Drevets et al, 1997; Blood et al, 2010; Seminowicz et al, 2004). Structural MRI studies have demonstrated that gray matter volume alterations within inferior prefrontal regions (including the orbitofrontal cortex), the anterior cingulate, and the hippocampus, are associated with depression (Ballmaier et al, 2004; Sheline et al, 2003; Bae et al, 2006; Bremner et al, 2002). In addition, deep white matter hyperintensities (WMH) have been linked to depression symptomatology, especially in late life depression (Kumar et al, 2002; Sassi et al, 2003).

Diffusion tensor imaging (DTI) is a MRI technique to study water diffusion in tissue, providing vital information on white matter microstructure, such as fiber connectivity and integrity in the healthy and diseased brain. Fractional anisotropy (FA) quantifies the directionally preferential water molecule diffusivity, and is usually used to assess overall white matter integrity. Using voxel- or region-based analyses, several DTI studies have revealed frontal and temporal white matter abnormalities in late-life depression (Nobuhara et al, 2006; Shimony et al, 2009), as well as in non-geriatric MDD (Wu et al, 2011; Ma et al, 2007; Zhu et al, 2011). For example, lower FA has been reported in frontal, temporal, and parietal regions for first-episode, treatment-naive younger adults with MDD (Ma et al, 2007). A more recent study acquired DTI data on middle-aged patients with MDD and controls and analyzed the images using the tract-based spatial statistics technique (Kieseppä et al, 2010). Suggestive areas of decreased FA were found in the right mid-cingulate cortex and posterior body of the corpus callosum in MDD patients compared with controls. In a somewhat related DTI study, altered white matter integrity in the anterior cingulate was found in patients with acute coronary syndrome who developed persistent depression when compared with patients who were not depressed (Rapp et al, 2010).

Results from treatment studies are mixed and at times contradictory (Taylor et al, 2008; Alexopoulos et al, 2002, 2008). For example, Taylor et al (2008) investigated DTI measures in subjects who underwent a 12-week open-label trial of sertraline, a selective serotonin reuptake inhibitor using a region-based approach. FA was measured in regions of interest placed in the white matter of the dorsolateral prefrontal cortex, anterior cingulate cortex, and corpus callosum. Subjects who did not remit with sertraline showed higher FA values in the superior frontal gyri and anterior cingulate cortices bilaterally. By contrast, in another treatment study, the authors reported lower FA values in multiple frontal limbic brain areas as well as select temporal and parietal regions associated with poor antidepressant response in late life depression (Alexopoulos et al, 2008). Although the precise relationship of FA alterations to treatment response is unclear, FA may represent the neuroanatomical substrates that are related to the pathophysiology of major depression.

On the basis of the evidence from the above volumetric and voxel-based DTI studies, suggesting the role of select regions of the prefrontal cortex, hippocampus, and cingulate gyrus in mood regulation and major depression, we hypothesize that microstructural white matter abnormalities may be present in two association fibers connecting these brain structures: the uncinate fasciculus (UF) and the cingulum bundle (CB). The UF is a ventral associative bundle that connects the anterior temporal lobe, including the amygdala and hippocampus, with the medial and lateral orbitofrontal cortex. The CB is a medial associative bundle that runs within the cingulate gyrus superior to the corpus callosum. It contains fibers of different lengths, the longest of which runs from the anterior temporal gyrus to the orbitofrontal cortex and the shortest (ie, U-fibers) connects the medial frontal, parietal, occipital, temporal lobes, and cingulate cortex. Both the UF and CB tracts, considered part of the limbic system, are thought to be involved in emotion processing, attention, and memory. Thus, a detailed examination of the tracts that connect these regions has implications for understanding the management and treatment of major depression.

As voxel-based DTI analyses may suffer from inter-subject mis-coregistration and lower spatial resolution because of smoothing in image pre-processing, we used a DTI-tractography technique to probe tract-specific white matter microstructural integrity of UF and CB in middle-aged major depression. By directly extracting and comparing tracts of interest, this DTI-tractography-based approach may provide high detection power in investigating white matter abnormalities. In addition to the standard diffusivity measures of white matter integrity, such as FA and mean diffusivity (MD), we explored tract-specific measures that quantify and summarize the geometric properties of the reconstructed tracts. We also computed two additional measures including the axial diffusivity (the degree of water diffusion along the direction parallel to the fiber bundles) and radial diffusivity (water diffusion perpendicular to the axonal wall). These two measures provide more information on the relative diffusivity of water molecules beyond standard FA, and have been utilized to differentiate among various white matter microstructural changes (Song et al, 2003; Metwalli et al, 2010; Bennett et al, 2010). We hypothesized that white matter abnormalities will be present in the UF and CB tracts of MDD subjects when compared with healthy controls (HC). Moreover, given previous reports of age effects on regional FA values that point toward compromised axonal morphology in aging (Sullivan et al, 2010) and studies of hemispheric asymmetries in brain structure (Gong et al, 2005), we assessed age and laterality effects for UF and CB tracts in a series of exploratory analyses using FA measures.

MATERIALS AND METHODS

Subjects

A total of 21 subjects with MDD and 21 age-matched HC were recruited from the greater Chicago area through flyers and local advertisements. This study was approved by the University of Illinois-Chicago Institutional Review Board, and written informed consent was obtained from each participant. The inclusion criteria for all subjects were right-handedness, age between 30 and 65 years, medication-naive or anti-depressant free for at least 2 weeks, and no history of diabetes, cardiac, or neurological diseases. The exclusion criteria included as folows: schizophrenia, bipolar or any psychotic disorders; current and past alcohol, or substance abuse; history of anxiety disorder outside of major depressive episodes; history of severe eating disorder and somatoform disorders.

All eligible subjects were assessed by a trained research assistant with the Structured Clinical Interview for Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV) (First et al, 1997). The severity of depression was quantified by a board certified/board eligible psychiatrist (AK or OA) using the 17-item Hamilton Depression Rating Scale (HAM-D) (Hamilton, 1967). At the time of enrollment, MDD subjects met DSM-IV criteria for MDD and required a score of 15 or greater on the HAM-D. Subjects were also administered the Center for Epidemiologic Studies of Depression (CES-D) Scale as an independent measure of depression severity (Radloff, 1977). The CES-D was used for correlation analyses as the HAM-D was the measure used in the determination of subject eligibility.

MRI Data Acquisition

All brain MRI data were acquired using Philips 3.0T Achieva scanner (Philips Medical Systems, Best, The Netherlands) with 8-channel SENSE (sensitivity encoding) head-coil. Subjects were fitted with soft ear plugs, positioned comfortably in the coil, and instructed to remain still. Foam pads were used to minimize head motion. For each subject, DTI images were acquired using single-shot spin-echo echo-planar imaging sequence (field of view (FOV)=240 mm; voxel size=0.83 × 0.83 × 2.2 mm; TR/TE=6994/71 ms; flip angle=9°). In all, 67 contiguous axial slices aligned to the anterior commissure–posterior commissure line were collected along 32 gradient directions with b=700 s/mm2 and one minimally diffusion-weighted scan (the b0 image). Parallel imaging technique was utilized with factor at 2.5 to reduce scanning time to 4 min. High-resolution three-dimensional T1-weighted images were acquired with MPRAGE (magnetization prepared rapid acquisition gradient echo) sequence (FOV=240 mm; 134 contiguous axial slices; TR/TE=8.4/3.9 ms; flip angle=8°; voxel size=1.1 × 1.1 × 1.1 mm). In addition, a T2-weighted FLAIR (fluid-attenuated inversion recovery) scan was also acquired with turbo spin echo sequence (FOV=240 mm; 67 contiguous axial slices; TR/TI/TE=11 000/2800/68 ms; voxel size=0.83 × 0.83 × 2.2 mm).

Image Processing

T2 FLAIR images were visually inspected to rule out cases with serious WMH, especially in the regions containing the uncinate and cingulum tracts. Visual inspection of all DTI image data was also conducted to ensure quality, and data with serious artifacts from substantial movements were removed. All diffusion weighted images (32 gradient directions) were first coregistered onto the b0 images using the automatic image registration (AIR) algorithm with affine transformation to minimize eddy currents and remove any potential small bulk motions that occurred during the scans (Woods et al, 1998). The diffusion directions in the diffusion gradient table were transformed by the rotation matrix calculated by AIR, to correct changes in section angulation due to coregistration. Diffusion tensor calculation and fiber tracking were then carried out using the DtiStudio software (Jiang et al, 2006; Laboratory of Brain Anatomical MRI, Johns Hopkins Medical Institute, Baltimore, MD). At each voxel, the signals from the 32 diffusion-weighted images were fitted to obtain the six elements of the diffusion tensor. The diffusion tensors were then diagonalized to obtain three eigenvalues (λ1, λ2, and λ3) and three eigenvectors (v1, v2, and v3). The eigenvector (v1) associated with the largest eigenvalue (λL: axial diffusivity) was assumed to represent the local fiber direction. Radial diffusivity (λT) is the mean of two minor eigenvalues (λ2+λ3)/2. To quantify the relative degree of anisotropy in a voxel, FA is calculated by the square root of the sum of squares (SRSS) of the diffusivity differences, divided by the SRSS of the diffusivities:

fiber tracking was performed with the fiber assignment by continuous tractography method (Jiang et al, 2006). For each subject, tractography was first performed with brute-force tracking on the whole brain by initiating tracts at each voxel. Fiber tracking was stopped when FA value falls <0.15 or a turning angle becomes >60°. To obtain a specific tract, a region of interest (ROI) was carefully traced and delineated. The extraction of UF followed the method described by (Wakana et al, 2007). We first identified and selected the most posterior coronal slice in which the temporal lobe is separated from the frontal lobe. A two-ROI approach was used with the first ROI including the entire temporal lobe and the second ROI the entire fiber projections in the frontal lobe. For the CB, we used a one-ROI approach, as recommended by (Catani and Thiebaut de Schotten, 2008). Because the majority of the fibers are short U-shaped fibers connecting adjacent gyri, the use of a two-ROI approach may exclude the majority in these short fibers from the analysis. We first examined the three axial slices above the corpus callosum and identified the slice with the largest continuous cingulum. The seeding ROI was then defined on the corresponding color-coded map with a single cigar-shaped region delineating the contour of the cingulum on this slice. Figure 1 demonstrates sample tract profiles superimposed on the corresponding FA map in a sagittal view. Inter-rater reliability on fiber extraction was performed by authors AZ and AL on a subset of five subjects (randomly selected from the entire sample of 42 subjects) and an intraclass correlation of >99.4% was reached in terms of the FA values for all the tracts. To ensure consistency in extracting fiber tracts, all reconstructed fibers were visually inspected for quality assurance by an experienced psychiatrist and brain imaging researcher (AL). None of the scans were rejected or considered unusable.

Figure 1
figure 1

Tract visualization superimposed on sagittal FA maps: (a) left uncinate fasciculus; (b) right uncinate fasciculus; (c) left cingulum bundle; and (d) right CB. The insert for each of the panel displays the regions of interest drawn on the color map.

PowerPoint slide

Quantitative Tract-Specific Measures

After tracts are obtained, the corresponding fiber data were used to calculate tract-specific metrics with a wide range of measures to quantify white matter characteristics along the tracts. These tract-specific measures can be categorized into two types: the diffusivity measures and the fiber geometric metrics (Table 1).

Table 1 Quantitative Tract-Specific Measures in Two Categories: Diffusivity and Fiber Geometry

Diffusivity measures included mean FA, axial (λL), and radial diffusivity (λT), MD. For each tract, these measures were extracted along the voxels that the fibers travel through and weighted by the number of reconstructed fibers in the particular voxel. All these values were then averaged to obtain a single number for each tract.

Fiber geometric measures included number of fibers (NNF), fiber volume (NFV), average fiber length (NAFL), and FA-weighted NAFL (NAFLFA). These geometric measures were normalized by the total intracranial volume, which was computed from the subject's high-resolution MPRAGE images using an automatic segmentation procedure in the FreeSurfer software (Fischl et al, 2004; Segonne et al, 2004; Jovicich et al, 2009; Martinos Center for Biomedical Imaging, Charlestown, Massachusetts).

Statistical Analyses

All statistical analyses were performed in the SPSS software, version 18.0 (SPSS, Somers, NY). Individual one-way analyses of covariance (ANCOVA) were used to assess group differences in all tract measures. Given the reported impact of age and sex on neuroanatomical measurements across various neuroimaging modalities (Barnes et al, 2010), these variables were used as covariates. We planned for post-hoc Pearson's product moment correlations (two-tailed) to further explore the relationship between depression severity scale CES-D and results of the aforementioned ANCOVAs in the MDD group. These analyses were also controlled for age and sex. We did not perform correlation with HAM-D as it was used in our inclusion criteria.

Exploratory analyses for age and laterality effect were performed for tract FA only. Pearson's product moment correlations (two-tailed) were conducted between tract FA and age. Analysis of laterality was performed on FA values of the left and right tract by paired-sample t-test. And then an asymmetry index (A) was computed which quantifies the differences between FA of the left and right tract, for both UF and CB using the following formula:

where FAL is the FA derived from the left UF or left CB and FAR is the FA of the right UF or right CB.

RESULTS

Demographic and Clinical Data

Table 2 summarizes the demographical and clinical characteristics of the MDD subjects and HC, including CES-D total score, Hamilton score, and duration of current symptoms in months. There were no significant differences between the MDD group and HC group in either age or sex. As expected, there were significant differences between the MDD group and HC across measures of depressive symptomatology (p<0.01).

Table 2 Demographic and Clinical Characteristics of Depressed (MDD) and Healthy Control Subjects

Between-Group Comparison

A significant FA reduction was observed for the right UF (F(1, 40)=6.86, p=0.01) in MDD subjects compared with HC. Within the right UF, radial diffusivity was significantly increased in the depressed group (F(1, 40)=4.46, p=0.04) whereas axial diffusivity remained relatively unchanged compared with HC. Figures 2a and b display FA and radial diffusivity respectively for all four tracts in both the MDD and HC groups. Post-hoc analyses revealed that FA and radial diffusivity in the right UF are negatively correlated (r=−0.735, p<0.001). Furthermore, the use of radial diffusivity as a covariate eliminated the group difference in FA for the right UF (F(1, 39)=2.27, p=0.140), suggesting a single physiological process responsible for FA and radial diffusivity changes. For MD, no significant between-group differences were detected for any of the four tracts. Similarly, the MDD and control groups did not differ on any tract geometric measure, including NNF, average length, and NFV for any of the four tracts.

Figure 2
figure 2

Mean values of: (a) FA; (b) radial diffusivity (λT) for MDD (N=21; in orange), and healthy control (HC) (N=21; in blue) groups. Error bars represent ±SD. aSignificant reduction in FA for the right UF (F(1, 40)=6.86, p=0.01). bSignificant increase in λT for the right UF (F(1, 40)=4.46, p=0.04).

PowerPoint slide

Correlations of Tract Measures to Depression Severity

To further investigate the relationship between depression severity and tract-specific measures of the right UF in the MDD group, two-tailed partial Pearson's correlation controlled for age and sex were performed between CES-D total score and tract-specific diffusivity as well as geometric measures. Results revealed that in the right UF, the normalized NNF (r=−0.53, p=0.02) significantly correlated with CES-D total score; diffusivity measures (eg, FA), average length measures, and NFV did not. Figure 3 shows the scatter plot of CES-D vs NNF in the right UF. In addition, in our exploratory analysis, we did not find any significant correlation between HAM-D and any of the tract measures. It should be noted that any potential correlation with HAM-D could be attenuated by the fact that the HAM-D scores were artificially thresholded at 15 and higher (as 15 was the cutoff point for a depressed subject to be eligible for the study).

Figure 3
figure 3

Scatter plot of CES-D score vs NNF in the right UF.

PowerPoint slide

Age Effect

A general trend of decreased FA was observed with aging, in all four tracts. Bilateral FA measures for the CB showed a significant negative correlation with age across the entire sample (left: r=−0.46, p=0.00; right: r=−0.49, p=0.00) and within the HC group only (left: r=−0.62, p=0.00; right: r=−0.58, p=0.00). Neither left nor right side of UF tract showed an association with age.

Laterality Effect

Paired-sample t-test on FA values showed that CB demonstrated significant laterality effect (L>R) in all three subject groups (MDD, HC, and the combined group; p<0.001). UF showed no significant laterality effect (L=R). Table 3 shows the asymmetry index values for the UF and the CB.

Table 3 Asymmetry Index Values for FA of Uncinate and Cingulum

DISCUSSION

For the first time, we have investigated tract-specific white matter alterations of UF and CB in non-geriatric MDD subjects utilizing DTI-tractography-based quantitative measures. The main finding of our study showed that white matter integrity of the right UF is significantly impaired in MDD subjects compared with HC as measured by FA. The decrease of FA in the right UF may be explained by the increase in the radial diffusivity, an index of myelin integrity, without significant alteration in the axial diffusivity. The geometric characteristics of the right UF did not contribute to the group differences in this tract. As concluded in the pioneering work of Song et al (2002, 2003), a decrease of axial diffusivity may reflect axonal injury whereas an increase of radial diffusivity may indicate a demyelination process. Thus, our results (unchanged axial but increased radial diffusivity) suggested that the white matter impairment of the right UF may be primarily demyelination in nature. This demyelination process could be due to neuronal loss, inflammation, or vascular damage as discussed in our previous work (Kumar and Cook, 2002). Lastly, newer reconstruction techniques based on high-angular-resolution-diffusion-weighted imaging are capable of disentangling the confound of variable degrees of sub-voxel fiber crossing, a variable previously known to alter axial and/or radial diffusivity (Tuch, 2004; Leow et al, 2009). Future studies may thus benefit from these novel techniques that allow us to differentiate among various neuropathological processes.

The UF connects the medial and lateral orbitofrontal cortex, the anterior cingulate cortex, and the anterior temporal lobe including the amygdala and hippocampus, all of which are structures known to be important in mood regulation. The FA reduction in the right UF implicates that the regions connected via UF could contribute to the impairment of the white matter integrity. As previously mentioned in the Introduction section, these regions have been previously investigated and shown to be altered in several imaging studies on mood disorders. For example, decreased radial diffusivity in the orbitomedial prefrontal cortex was reported in bipolar subjects relative to HC (Versace et al, 2008) and decreased FA in the anterior cingulate cortex was reported in subjects with acute coronary syndrome who developed persistent depression (Rapp et al, 2010). However, further studies are needed to directly examine which specific regions contribute to the FA reduction in the right UF.

Although the exact function of the UF in MDD is not well understood, a recent study of late-onset major depression found that depression severity significantly correlated with lower FA in white matter lesions of the right UF (Dalby et al, 2010). In contrast, Taylor et al (2007) found early-onset geriatric depressed subjects exhibited lower FA of the left UF compared with mid and late onset or HC subjects. Furthermore, converging lines of evidence implicated the possible involvement of UF in other neuropsychiatric illnesses that often overlap with MDD including anxiety disorders (Phan et al, 2009) and bipolar affective disorder (McIntosh et al, 2008; Versace et al, 2008). Compared with HC subjects, patients with anxiety disorder had significantly lower FA localized to the right UF white matter near the orbitofrontal cortex (Phan et al, 2009). Likewise, patients with bipolar disorder showed similar abnormalities in the UF (McIntosh et al, 2008; Versace et al, 2008; Versace et al, 2010). Although we did not find abnormalities in the left UF in our depressed subjects, a previous study demonstrated that bipolar patients exhibited divergent FA changes in the left vs the right UF (greater FA in the left UF and lower FA in the right UF) (Versace et al, 2008). Further evidence of the involvement of the right amygdala-orbitofrontal cortical functional connectivity was found using fMRI (Versace et al, 2010). Thus, alterations in the structural and functional connectivity of this circuit may represent a marker of mood disturbances in general. Overall, our findings pointing to the right UF involvement in mid-life MDD adds to the growing literature on the importance of this white matter tract to psychiatric illnesses involving mood and affective disturbances.

Although the severity of depression as measured by CES-D in MDD subjects did not correlate with FA, it did negatively correlate with the geometric measure of normalized NNF. To our knowledge, this is the first report of such findings in the literature. Here, we note that the NNF as a measure may be dependent on the protocols or procedures used, and thus can be problematic in multiple- or cross-site studies. Also, NNF summarizes one aspect of the mathematical representation of white matter tracts, rather than the actual count of fibers, and thus its real physiological meaning is not straightforward. Nevertheless, it remains a valid measure as in this paper we used the same procedure/protocol (eg, whole brain tractography by seeding at the center of each voxel with same stopping criteria of FA <0.15 and/or bending angle >60°), which was executed by two of the authors AZ and AL. As FA did not correlate with CES-D in the depressed group (but NNF did), our results support that some of the novel measures (such as NNF) introduced in this paper may be useful imaging metrics that provide complementary information on white matter microstructure when combined with other conventional DTI measures. Although further studies are needed to better understand the neuroanatomical nature of this correlation and its possible clinical implications, evidence from postmortem studies may provide some initial insight. Neuropathological analyses in depressed subjects have shown abnormal reduction in glial cells in the subgenual prefrontal cortex, orbital cortex, dorsal anterolateral prefrontal, and amygdala (Manji et al, 2001). Such reduction may negatively impact glial functions integral to structural plasticity in the brain resulting in a reduction of fiber number in MDD (Coyle and Schwarcz, 2000).

White matter provides the primary structural substrate for brain circuitry by connecting critical cortical and subcortical regions. The UF, in particular, connects the amgydala and hippocampus with the lateral and medial orbitofrontal cortex. Impaired connectivity between these important regions, reflected by decreased FA or reduced fiber number, may have a key role in the pathophysiology of MDD. Furthermore, functional imaging studies demonstrate the importance of intact white matter connectivity to successful mood regulation (Drevets, 2001). Abnormal functional connectivity in the subgenual cingulate and thalamus has been associated with major depression using resting-state functional MRI (Greicius et al, 2007). On the basis of similar functional imaging data, suggesting the over-activation of the cingulate region (Brodmann area 25) in MDD, Seminowicz et al (2004) found that the anatomical connectivity of the subgenual cingulate region targeted with deep brain stimulation for depression supports the hypothesis that treatment efficacy is mediated via effects on a distributed network of frontal, limbic, and visceromotor brain regions (Johansen-Berg et al, 2008). Given that the UF is a key white matter tract coursing through the anterior cingulate cortex and amygdala, our findings may have implications for treatment studies as well.

Exploratory analyses of age and laterality effects revealed that bilateral FA values in the CB correlated with age in HC, and across each subject group the CB exhibited significant asymmetry (L>R). The fact that the age and FA correlation in the combined sample was driven by HC subjects suggests that the presence of depression may attenuate the effect of age on FA of the cingulum. In contrast to the attenuated age effect for the CB, the presence of depression did not alter the laterality properties of either tract. We theorize that greater FA found in the left CB may reflect a higher degree of axonal myelination or orientational/directional coherence of fibers in the left CB in comparison to the right side. Thus, the asymmetry of the CB may reflect the structural and functional differences between the two hemispheres, a finding that is consistent with the literature (Gong et al, 2005; Kubicki et al, 2003). More directed studies of age and laterality effects using DTI-derived measures of white matter integrity and geometry should be conducted.

Lastly, our study is limited by a relatively small sample size. Although previous diffusion imaging studies revealed similar results, the small sample size may limit the generalizability of these findings. In addition, with regards to the DTI acquisition, we applied the well-established protocol optimized for Philips 3T scanner, with a SENSE factor of 2.5 chosen to well balance the scan time and the signal-to-noise ratio (Reich et al, 2006). Note that the sequence did not use the dual spin echo technique (Reese et al, 2003) to reduce eddy current effects which in our study are less prominent with a low b-value of 700 s/mm2. Moreover, any residual effects were further minimized by the coregistation step in image post-processing before tensor fitting.

In summary, based on previous findings suggesting the role of select regions of the frontal cortex, hippocampus, and cingulate gyrus in mood regulation and major depression, we investigated the white matter microstructural integrity of two important association tracts connecting these regions: the UF and the CB in MDD using DTI-tractography. Our study makes the following novel contributions in clarifying the tract-specific abnormalities in major depression: (i) supports the role of the right UF in the pathology of MDD; (ii) demonstrates that radial (but not axial) diffusivity is altered in the right UF despite intact fiber geometry; and (iii) associates the severity of depression in MDD subjects with NNF in the right UF. We believe that this study points to the role of demyelination and reduced fiber number, highlighting the importance of abnormalities in major white matter association tracts in MDD.