Aberrant intrinsic connectivity of hippocampus and amygdala overlap in the fronto-insular and dorsomedial-prefrontal cortex in major depressive disorder

Introduction

Major depressive disorder (MDD) is one of the most common psychiatric disorders with a lifetime prevalence of about 16% (Kessler et al., 2003). It is characterized by episodes of prolonged depressed mood, reduced energy, impaired cognition, vegetative symptoms, and suicidal tendencies with suicide rates of nearly 4% (American Psychiatric Association, 2000). Recent work indicates that MDD is associated with dysfunction of emotion/emotion regulation-related circuitry that includes the hippocampus (HP) and amygdale (AY) among other brain regions (Drevets et al., 2008). The strong consistency of structural and functional HP and AY changes observed across a number of studies highlights the prominent role of these two brain regions in the pathophysiology of MDD (Rigucci et al., 2010; Kempton et al., 2011; MacQueen and Frodl, 2011; Murray et al., 2011). The current study assessed HP and AY pathophysiology in MDD with special focus on the HP’s and AY’s intrinsic functional connectivity (iFC), an index of the synchrony between ongoing HP/AY activity and the ongoing activity of other brain areas.

The HP and AY are located adjacent to one another within the medial temporal lobes, and interact extensively with one another as well as within widely distributed cortical-subcortical circuits (Squire et al., 2004; LeDoux, 2007; Aggleton, 2012). Together, they contribute to several cognitive-behavioral functions, particularly to emotional memory (LaBar and Cabeza, 2006), which is impaired in MDD (e.g., increased negative memory bias in MDD) (Hamilton and Gotlib, 2008).

Resting-state functional MRI (rs-fMRI), which is commonly used to analyze intrinsic connectivity (Fox and Raichle, 2007), has revealed specific iFC patterns for both HP and AY in healthy subjects (Kahn et al., 2008; Etkin et al., 2009; Roy et al., 2009). These HP and AY iFC patterns are altered in patients with MDD (Cao et al., 2012; Tang et al., 2012). For example, resting-state activity within the posterior HP is positively correlated with areas of the “task-negative” network (Fox et al., 2005; Fox and Raichle, 2007) and negatively correlated with resting-state activity within a widely distributed set of brain regions involved in a large range of tasks [i.e., “task-positive” network including the so-called salience and central executive network (Fox et al., 2005; Dosenbach et al., 2007; Seeley et al., 2007)]. The default mode, salience, and central executive networks are intrinsic networks of synchronized, ongoing activity that are associated with specific symptoms (e.g., rumination) when disrupted in MDD (Greicius et al., 2007; Sheline et al., 2009, 2010; Hamilton et al., 2012, 2013). Recent research has demonstrated aberrant hippocampal iFC with areas of the task-negative network in MDD (Cao et al., 2012). Prior research has also shown that resting-state activity within the basolateral AY is positively correlated with the activity of primary areas (such as those of the visual or sensorimotor system), the medial prefrontal cortex (PFC), and of the medial temporal lobe, but negatively correlated with areas of the salience network [i.e., the dorsomedial PFC and fronto-insular operculum (FIO)] (Etkin et al., 2009; Roy et al., 2009; Veer et al., 2010). This last finding is also supported by findings in rats (i.e., negative correlations between ongoing amygdala activity and activity in a frontolimbic circuit), suggesting the robustness of the AY intrinsic connectivity pattern across species (Liang et al., 2012). Recent research has also observed aberrant amygdala iFC with the PFC in first-episode MDD (Tang et al., 2012). Due to the intimate interaction between the HP and AY as well as the partly overlapping cortical iFC patterns of the HP and AY, the question arises as to whether aberrant HP and AY iFC to the cortex might be linked in MDD. This question is important for a better understanding of the pathophysiology of MDD, because a correspondence between aberrant HP and AY iFC may link the medial temporal lobe changes of MDD to the aberrant cortical intrinsic networks in MDD (Menon, 2011; Hamilton et al., 2013). Based on the studies described above, we hypothesized that MDD patients would show a spatially overlapping pattern of aberrant HP and AY intrinsic connectivity with the cortex.

To address this hypothesis, we assessed 21 patients with MDD and 20 age, gender, and education matched healthy controls using blood oxygenation level dependent (BOLD) rs-fMRI signal fluctuations. The primary outcome measure for this study was ß-maps reflecting region-of-interest (ROI) based positive and negative BOLD correlations as a surrogate of iFC. ROIs were positioned bilaterally within the basolateral AY and body of HP. These regions have extensive cortical iFC covering intrinsic networks such as default mode, salience, and central executive network (Kahn et al., 2008; Etkin et al., 2009; Roy et al., 2009). Since we were not interested in lateralized effects of left and right AY/HP, iFC, group comparisons for bilateral AY and HP connectivity were framed by an analysis of variance (ANOVA) approach with factors group and brain side, in which the main effect of group was the effect of interest. Prior research has observed structural changes in the HP and AY of MDD patients (Kempton et al., 2011). Therefore, covariates for HP and AY gray matter (GM) volume were included in the analysis to ensure our findings were independent of differences in HP and AY structure between MDD and healthy control groups.

Materials and Methods

Participants

Twenty-one patients with recurrent MDD [female/male: 11/10; mean age: 51.0 (SD 15.0)] and 20 healthy controls [female/male: 11/9; mean age: 49.6 (SD 13.9)] participated in this study (Table 1). After receiving the approval from the medical ethical board of Technische Universität München (TUM), all participants provided informed consent in accordance with the Human Research Committee guidelines of the Klinikum Rechts der Isar, TUM. All patients were inpatients recruited from the Department of Psychiatry, TUM. Controls were recruited from the larger community by word-of-mouth advertising and group-matched by age, gender, and education to MDD patients. Participants’ examination included medical history, psychiatric interview, psychometric assessment, and blood tests for patients. Psychiatric diagnoses were based on DSM-IV (American Psychiatric Association, 2000). The Structured Clinical Interview for DSM-IV (CID-I; First et al., 1996) was used to assess the presence of psychiatric diagnoses. The severity of clinical symptoms was measured with the Hamilton Rating Scale for Depression (HAM-D) (Hamilton, 1960) as well as the Beck Depression Inventory (Beck et al., 1961). The global level of social, occupational, and psychological functioning was measured with the Global Assessment of Functioning (GAF) scale (Endicott et al., 1976). Psychiatrists Dirk Schwerthöffer and Martin Scherr performed clinical-psychometric assessment. They have been professionally trained for SCID interviews with inter-rater reliability for diagnosis of more than 95%.

TABLE 1

Table 1. Demographic and clinical data.

Major depressive disorder was the primary diagnosis for all patients. All patients met criteria for a current major depressive episode with an average episode length of 16.2 weeks (range 8–31 weeks). The average age of MDD onset was 35 years, and all patients experienced their first MDD episode before 50 years of age. The average total duration of MDD was 14.7 years and on average, patients had experienced five depressive episodes. Four patients had a positive family history of MDD. Thirteen MDD patients had co-morbid diagnoses: six generalized anxiety disorder, two somatization disorder, and five avoidant or dependent personality disorders. Patients with psychotic symptoms, schizophrenia, bipolar disorder, and substance abuse were excluded from this study. Additional exclusion criteria were age below 18 or above 70 years, pregnancy, neurological or internal systemic diseases, and general contraindications for MRI. One patient was free of any psychotropic medication during MRI assessment. Five patients were treated by antidepressant mono-therapy (two cases: citalopram 30 mg/day (mean dose), two cases: sertraline 200 mg/day, one case: mirtazapine 30 mg/day), 10 patients by dual-therapy (four cases: citalopram 37.5 mg/day; mirtazapine 30 mg/day, one case: citalopram 40 mg/day; venlafaxine 225 mg/day, one case: citalopram 30 mg/day; quetiapine 200 mg/day, one case: sertraline 200 mg/day; mirtazapine 30 mg/day, three cases: venlafaxine 225 mg/day; mirtazapine 30 mg/day), and five patients by tri-therapy (two cases: citalopram 30 mg/day; venlafaxine 187.5 mg/day; amisulprid 200 mg/day, two cases: citalopram 30 mg/day; mirtazapine 30 mg/day; quetiapine 200 mg/day, one case: venlafaxine 22 mg/day; mirtazapine 30 mg/day; quetiapine 200 mg/day). All healthy controls were free of any current or past neurological or psychiatric disorder or psychotropic medication.

All participants underwent a 5-min structural MRI and 10-min rs-fMRI with the instruction to keep their eyes closed and not to fall asleep. Subjects were questioned via intercom immediately after the rs-fMRI scan to verify they remained awake during the scan. Before and after scanning, a medical examination of patients validated their condition was stable and determined whether they experienced any odd sensations during the scanning. No unusual sensations were reported, and no patient dropped out during the scanning session.

MRI Data Acquisition

MRI was performed on a 3-T whole body MR scanner (Achieva, Philips, Netherlands) using an eight-channel phased-array head coil. fMRI data were collected using a gradient echo EPI sequence (TE = 35 ms, TR = 2000 ms, flip angle = 82°, FoV = 220 mm × 220 mm, matrix = 80 × 80, 32 slices, slice thickness = 4 mm, and 0 mm interslice gap; 10 min of scanning resulting in 300 volumes). T1-weighted structural data were obtained using a magnetization-prepared rapid acquisition gradient echo sequence (TE = 4 ms, TR = 9 ms, TI = 100 ms, flip angle = 5°, FoV = 240 mm × 240 mm, matrix = 240 × 240, 170 slices, voxel size = 1 mm × 1 mm × 1 mm).

rs-fMRI Data Analysis

The first three functional images of each subject’s dataset were discarded due to magnetization effects. The remaining rs-fMRI data were preprocessed in SPM8 (Wellcome Department of Cognitive Neurology, London) including head motion correction, spatial normalization into the stereotactic space of the Montreal Neurological Institute (MNI) with isotropic voxels of 3 mm × 3 mm × 3 mm, and spatial smoothing with a 6 mm × 6 mm × 6 mm Gaussian kernel to reduce spatial noise. To ensure data quality, particularly concerning motion-induced artifacts, temporal signal-to-noise ratio (tSNR) and point-to-point head motion were estimated for each subject (Murphy et al., 2007; Van Dijk et al., 2012). Excessive head motion (cumulative motion translation or rotation >3 mm or 3° and mean point-to-point translation or rotation >0.15 mm or 0.1°) was applied as exclusion criterion. Point-to-point motion was defined as the absolute displacement of each brain volume compared to its previous volume. None of the participants had to be excluded. Two-sample t-tests yielded no significant differences between groups in mean point-to-point translation or rotation in any direction (p > 0.10) or in tSNR (p > 0.50). Further control for head motion effects was carried out at subject-level iFC analysis.

Subject-level iFC analysis

Seeds for the iFC analysis were positioned in the basolateral AY (Roy et al., 2009) and the body of HP (Kahn et al., 2008) (Figure 1). Both of these seeds are characterized by an extensive neocortical iFC pattern, which covers the above-mentioned intrinsic networks involved in MDD (Greicius et al., 2007; Sheline et al., 2010; Hamilton et al., 2011). Left and right basolateral AY were defined by the Anatomy Toolbox (Wellcome Department of Cognitive Neurology, London) in SPM and converted to ROIs using MarsBaR ¹. The Anatomy Toolbox integrates probabilistic cytoarchitectonic maps derived from human post-mortem studies into the SPM environment. MarsBaR was used to create spherical ROIs (4 mm radius) with center coordinates (−24 −18 −18 and +24 −18 −18) in the body of the posterior HP (Kahn et al., 2008). After Butterworth bandpass-filtering of all voxel time courses for the frequency range from 0.009 to 0.08 Hz, we extracted voxel time courses of seed ROIs and reduced them to ROI-representative time courses by singular value decomposition, respectively. Each time course was put into a first-level fixed-effects general linear model in SPM8, and four separate iFC analyses (i.e., left AY, right AY, left HP, and right HP) were performed for each subject. For each model, additional regressors for global GM, white matter (WM), cerebral spinal fluid (CSF) BOLD-signal, and six movement parameters (three translational and three rotational directions) were included as covariates of no interest.

FIGURE 1

Figure 1. Regional hippocampus (HP) and amygdala (AY) volumes. Seeds for the iFC analysis were positioned in the basolateral AY (yellow) and the body of HP (green). Regional voxel-based morphometry (VBM) volumes from the basolateral AY and body of HP were averaged across hemispheres and compared between patients with major depressive disorder (MDD) and healthy controls (Con) using two-sample t-tests (p < 0.05). HP (p = 0.011) and AY (p = 0.019) volumes were reduced in MDD patients.

Group level iFC analysis

Group analyses were performed using contrast images from the subject-level iFC analysis in two separate flexible factorial models of ANOVA that included covariates of no interest [sex, age, and seed voxel-based morphometry (VBM) volume (see below for VBM analysis)]. The first ANOVA model was performed with factors group (with levels MDD group and healthy control group) and ROI (with levels left and right AY). The second ANOVA model was performed with factors group and ROI (with levels left and right HP). For both ANOVA models, appropriate post hoc t-tests were used to reveal the positively and negatively correlated iFC pattern for HP and AY that was representative of each study group (Roy et al., 2009). For both ANOVA models, the main effect of group (and corresponding post hoc t-tests to reveal direction of change) was the effect of interest. Statistical thresholds were set to p < 0.05 whole brain corrected for false discovery rate (FDR). Reported voxel coordinates correspond to standardized MNI space.

Voxel-Based Morphometry

To determine structural HP and AY changes in MDD and to control for their effects on iFC results, VBM of structural MRI data was performed. VBM analysis followed a protocol described in previous work (Sorg et al., 2012). VBM8 toolbox ² was used for data preprocessing and analysis. Images were corrected for bias-field inhomogeneity, registered using linear (12-parameter affine) and non-linear transformations, and tissue-classified into GM, WM, and CSF within the same generative model (Ashburner and Friston, 2005). GM images were modulated to account for volume changes based on the normalization process. We only considered non-linear volume changes so that subsequent analyses did not have to account for differences in head size. Finally images were smoothed with an 8-mm (FWHM) Gaussian kernel. Regional VBM volumes of HP and AY (mean of left and right HP and AY, respectively) were calculated for the body of HP and basolateral AY. Group differences were assessed by two-sample t-tests (p < 0.05). Regional VBM volumes were included as covariates-of-no interest in the group iFC analysis.

Results

Reduced HP and AY Volume in Patients

Regional volumes of both body of HP and basolateral AY were evaluated by comparing regional VBM volumes between patients with MDD and healthy subjects by the use of two-sample t-test. Both HP (p = 0.011) and AY (p = 0.019) volumes were reduced in patients (Figure 1).

Overlapping Reduced Intrinsic Connectivity of HP and AY in Patients

Figure 2 presents spatial maps of BOLD FC for the body of HP and basolateral AY in MDD patients and healthy controls. Maps are based on post hoc t-tests of different ANOVA models, which included corresponding covariates of regional seed VBM volumes (p < 0.05, FDR for multiple comparisons; Table 2). Results are largely consistent with previous findings (Kahn et al., 2008; Etkin et al., 2009; Roy et al., 2009; Cao et al., 2012). Ongoing AY activity was positively correlated with BOLD activity within the medial temporal lobe and in areas of primary sensorimotor and visual cortex in both patients and controls. HP activity was positively correlated with BOLD activity within regions of the default mode network such as the medial PFC and posterior cingulate cortex. Both AY and HP activity was negatively correlated with BOLD activity within regions of the salience network (such as dorsomedial PFC and FIO) and inferior parietal lobule. Differences between MDD and healthy participants were observed for HP and AY-FC within the dorsomedial PFC and FIO. Specifically, patients showed less negative BOLD-correlation between these brain regions and HP and AY. Further, these changes in HP and AY iFC overlapped (p < 0.05; FDR corrected; Figure 3, Table 2).

FIGURE 2

Figure 2. Intrinsic functional connectivity (iFC) patterns of hippocampus (HP) and amygdala (AY) in patients and healthy controls. Individual spatial ß-maps representing BOLD correlations of ongoing left and right HP and AY activity of patients with major depressive disorder (MDD) and healthy controls (Con), were analyzed using two ANOVA models with factors for group and hemisphere for the HP and AY, respectively. Yellow and green maps, which were superimposed on a single-subject high resolution T1 image, represent results of corresponding post hoc t-tests that reflect positive and negative FC (Pos-FC and Neg-FC, respectively) for each group and seed region (p < 0.05, FDR corrected; bars represent range of t-values).

FIGURE 3

Figure 3. Overlapping differences in the intrinsic functional connectivity (iFC) of the hippocampus (HP) and amygdala (AY) between patients with major depressive disorder (MDD) and healthy controls (Con). Individual spatial ß-maps, which represent BOLD correlations of ongoing left and right HP and AY BOLD activity of patients with major depressive disorder (MDD) and healthy controls (Con), were analyzed using two ANOVA models with factors for group and hemisphere for AY as well as HP. Yellow and green maps of (A,B) were superimposed on a single-subject’s high resolution T1 image, and represent the main effect of group on the FC of the AY and HP, respectively (p < 0.05, FDR corrected; bars represent range of t-values). Bar graphs on the right side reflect averaged iFC of AY and HP for MDD patients and healthy controls, respectively (two-sample t-tests, p < 0.05). In (C), regions of overlapping reduced FC of AY (A) and HP (B) are shown in red. The HP and AY have reduced FC in the dorsomedial-prefrontal cortex and fronto-insular operculum in patients.

TABLE 2

Table 2. Regions of intrinsic functional connectivity (iFC) with amygdala (AY) and hippocampus (HP).

Discussion

The present study investigated differences in functional connectivity of HP and AY in patients with MDD and healthy controls to provide new insight into the pathophysiology of MDD. We found that ongoing BOLD activity within the body of HP and basolateral AY was negatively correlated with activity within the dorsomedial PFC and FIO. However, the negative FC between these brain regions was weaker in MDD, independent of patients’ reduced HP and AY GM volumes. These findings provide the first evidence that the HP and AY show a corresponding reduction in iFC to other brain regions in MDD. Since dorsomedial PFC and FIO are critical for intrinsic network interactions (i.e., default mode, salience, and central executive network), which are impaired in MDD, these findings suggest that aberrant HP and AY intrinsic connectivity may interact with aberrant intrinsic network activity in MDD.

Overlapping Aberrant HP and AY Intrinsic Connectivity to Prefrontal Cortex and Fronto-Insular Operculum in MDD

The iFC analyses of the present study revealed that during the resting-state, both HP and AY BOLD-functional connectivity with dorsomedial PFC and FIO was reduced in MDD (Figure 3; Table 2). Healthy controls showed a strong negative correlation between basolateral AY and the dorsolateral PFC, dorsomedial PFC, FIO, and inferior parietal lobule (Figure 2). However, AY-FC with the dorsomedial PFC, FIO, and inferior frontal gyrus were substantially reduced in patients. A very similar, and importantly, overlapping pattern was observed in the iFC analysis for the HP (Figure 3; Table 2). While healthy control participants showed a strong negative correlation between HP and the PFC (i.e., dorsomedial PFC and FIO) (Figure 2), patients showed weaker FC between the HP and these brain regions. Reduced FC of both AY and HP to PFC and FIO cannot be explained by decreased AY and HP volumes in patients, as these volumes were included in the iFC analyses as covariates of no interest. Each of these brain regions has been implicated in emotion-related processes (Ochsner and Gross, 2005; Wood et al., 2012). Further, these brain regions appear to play a role in the regulation of emotion (Ochsner and Gross, 2005). For example, both medial and lateral regions of the PFC support processes related to the cognitive control of emotion (Ochsner and Gross, 2005; Kalisch, 2009; Winecoff et al., 2011). Furthermore, AY-FC during emotion regulation was decreased within the FIO and dorsomedial PFC in young MDD patients, perfectly matching our finding (Perlman et al., 2012).

Similarly, it has been demonstrated that there is a robust decrease of resting-state functional connectivity of bilateral amygdala, anterior insula, and anterior cingulate cortex within the salience network and ventromedial prefrontal cortex (vmPFC) and temporal poles in medication-free MDD patients with no comorbidity (Veer et al., 2010). In addition, Cao et al. (2012) found that drug-naïve patients with MDD showed selectively less negative BOLD-correlation between ongoing HP and lateral PFC activity, while Tang et al. (2012) observed reduced FC between AY and ventral PFC. The findings from both of these prior studies are consistent with our current results. In summary, our findings demonstrate there is significant overlap in the aberrant FC of the HP and AY to both the dorsomedial PFC and FIO in MDD.

A key finding in the present study is the reduction in both HP and AY iFC with the dorsomedial PFC in MDD patients. The HP and AY have well established roles in learning, memory, and particularly emotional memory (LaBar and Cabeza, 2006; LeDoux, 2007). Impaired emotional memory may play a critical role in the pathogenesis of MDD (Segal et al., 1996; Kendler et al., 2001; Robinson and Sahakian, 2008). Further, recent research indicates that depression is associated with aberrant iFC in the PFC (referred to as the dorsal nexus) within several intrinsic networks (Sheline et al., 2010). More specifically, depression is associated with an increase in iFC between the dorsomedial PFC and default mode, salience, and executive control networks (Sheline et al., 2010; Scheidegger et al., 2012). This prior work suggests that aberrant functional connectivity between the dorsomedial PFC and these intrinsic brain networks may be an important aspect of the pathophysiology of depression. These changes in functional connectivity may provide a mechanism that explains the co-occurrence of depressive symptoms, such as increased self-focus, that are mediated by these distinct brain networks (Sheline et al., 2010). The present study extends this prior work by demonstrating that both the HP and AY also show aberrant iFC with the dorsomedial PFC in MDD patients. The changes observed in patients’ HP and AY iFC with the dorsomedial PFC (i.e., the dorsal nexus) may provide a mechanism by which disruptions in emotional memory act synergistically with other cognitive and behavioral symptoms of depression.

The HP and AY also showed aberrant iFC with the FIO in patients compared to healthy participants. More specifically, both the HP and AY showed reduced iFC with the FIO. The FIO is considered the hub of a neural circuit that detects salient events (i.e., the salience network), facilitates access to attention and working memory resources, and interacts with other intrinsic networks (i.e., the default mode and central executive network) to generate appropriate responses to salient stimuli (Seeley et al., 2007; Menon, 2011). Psychopathological processes may be mediated by dysfunction of the salience network. Given that a minimal amount of salience network activity is necessary for healthy cognitive-emotional functioning, a decrease in network activity may result in psychopathology (e.g., depression). For example, prior work suggests the FIO initiates adaptive emotion regulation and network interaction processes that are triggered by the rumination of depressed individuals (Hamilton et al., 2011). In the present study, the HP and AY showed impaired iFC to the FIO in MDD patients. This reduction in iFC to the salience network may provide a mechanism for the impaired emotional memory associated with the HP and AY dysfunction that contributes to impaired salience detection. Thus, disruption of salience detection processes via decreased HP and AY iFC to the FIO may play a critical role in the pathophysiology of MDD.

Voxel-based morphometry findings from the current study replicate prior research investigating HP and AY volume changes associated with MDD, indicating the representative character of our patient sample and thereby of our iFC results. We observed AY atrophy in MDD patients compared to controls (Figure 1). Findings from prior research investigating changes in AY volume associated with MDD have been mixed (Frodl et al., 2003, 2008a,b; Kempton et al., 2011). In addition, we observed significant HP atrophy in MDD patients compared to healthy controls (Figure 1). Similar results have been demonstrated in prior MDD neuroimaging studies (Sheline et al., 1996; Kronmuller et al., 2008). These findings provide further evidence that structural changes in the HP are associated with MDD.

Limitations and Methodological Issues

Patients in our study were treated with antidepressant medications. It is difficult to statistically control for antidepressant effects because there is no reliable tool available to make drug effects comparable across different antidepressants [in contrast to antipsychotic drugs, which can be mapped on chlorpromazine equivalents CPZ and included as additional covariates in analyses (e.g., Sorg et al., 2012)]. Due to the complexity of this issue, the impact of antidepressants on iFC is not completely understood (Bruhl et al., 2010; Delaveau et al., 2011). Recent studies suggest antidepressant medications normalize brain function (Anand et al., 2005; Fu et al., 2007; Heller et al., 2013). Consistent with this suggested normalization of function, the reduction in HP and AY iFC in medicated patients in the present study is largely consistent with previous findings in treatment-naïve patients with MDD (Cao et al., 2012; Tang et al., 2012). For example Cao et al. (2012) reported weaker negative BOLD-correlations between the HP and the middle frontal gyrus in non-medicated patients that are similar to our present findings in medicated patients. However, although our findings fit previous results, we cannot exclude medication effects for our results. Therefore, the present results should be interpreted with caution. Future studies in non-medicated patients are necessary to replicate and extend our findings.

We used standard rs-fMRI analysis procedures to calculate our main outcome measure of ROI-based FC ß-maps. These procedures included global signal regression to remove physiological noise (such as respiratory and cardiac based signals) from the resting-state signal. However, there is an ongoing debate regarding whether global signal regression induces artificial BOLD correlations, particularly negative correlations (Murphy et al., 2007; Chang and Glover, 2009; Fox et al., 2009; Chai et al., 2012). We decided to use global signal regression for the following reasons: (i) recent studies have demonstrated the biological origin of negatively correlated BOLD FC during rest that is independent of global signal regression application (Chai et al., 2012; Wong et al., 2012). (ii) Particularly for the AY, negative BOLD correlation with frontolimbic areas have recently been found in awake rats, independent of global signal regression (Liang et al., 2012). This finding demonstrates the robustness of the negative BOLD-correlation network seeded by AY across species, as well as its independence from global signal regression. (iii) Previous studies in healthy controls and MDD used global signal regression when investigating the iFC of the medial temporal lobe (Roy et al., 2009; Cao et al., 2012). Therefore, to better compare our findings with these prior studies, we used the same approach.

Conclusion

In patients with MDD, reduced HY and AY intrinsic connectivity overlaps within the dorsomedial PFC and FIO. These regions not only support the expression and regulation of emotion, which is impaired in MDD, but also regulate intrinsic brain network interactions. Since these networks and their interactions are disrupted in MDD, our finding suggests that there is a link between the medial temporal lobe and intrinsic network pathophysiology in MDD that is mediated aberrant HP and AY intrinsic connectivity.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by German Federal Ministry of Education and Research (BMBF 01EV0710 to Afra M. Wohlschläger), German research foundation (DFG) grants No. DR 445/4-1 and DR 445/5-2 to Alexander Drzezga), the Kommission für Klinische Forschung of the Klinikum Rechts der Isar der Technischen Universität München (KKF 8765162 to Christian Sorg) and the German Academic Exchange service (DAAD) for financial support of this project (Masoud Tahmasian). We are grateful to the participants of the study and the staff of the Department of Psychiatry and Neuroradiology for their help in recruitment and data collection.

Footnotes

1. ^http://marsbar.sourceforge.net
2. ^http://dbm.neuro.uni-jena.de/vbm.html

References

Aggleton, J. P. (2012). Multiple anatomical systems embedded within the primate medial temporal lobe: implications for hippocampal function. Neurosci. Biobehav. Rev. 36, 1579–1596. doi:10.1016/j.neubiorev.2011.09.005

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

American Psychiatric Association. (2000). Diagnostic and Statistical Manual of Mental Disorders (DSM-IV-TR), 4th Edn. Washington, DC: American Psychiatric Association.

Anand, A., Li, Y., Wang, Y., Wu, J., Gao, S., Bukhari, L., et al. (2005). Activity and connectivity of brain mood regulating circuit in depression: a functional magnetic resonance study. Biol. Psychiatry 57, 1079–1088. doi:10.1016/j.biopsych.2005.02.021

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ashburner, J., and Friston, K. J. (2005). Unified segmentation. Neuroimage 26, 839–851. doi:10.1016/j.neuroimage.2005.02.018

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Beck, A. T., Ward, C. H., Mendelson, M., Mock, J., and Erbaugh, J. (1961). An inventory for measuring depression. Arch. Gen. Psychiatry 4, 561–571. doi:10.1001/archpsyc.1961.01710120031004

CrossRef Full Text

Bruhl, A. B., Kaffenberger, T., and Herwig, U. (2010). Serotonergic and noradrenergic modulation of emotion processing by single dose antidepressants. Neuropsychopharmacology 35, 521–533. doi:10.1038/npp.2009.159

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cao, X., Liu, Z., Xu, C., Li, J., Gao, Q., Sun, N., et al. (2012). Disrupted resting-state functional connectivity of the hippocampus in medication-naive patients with major depressive disorder. J. Affect. Disord. 141, 194–203. doi:10.1016/j.jad.2012.03.002

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Chai, X. J., Castanon, A. N., Ongur, D., and Whitfield-Gabrieli, S. (2012). Anticorrelations in resting state networks without global signal regression. Neuroimage 59, 1420–1428. doi:10.1016/j.neuroimage.2011.08.048

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Chang, C., and Glover, G. H. (2009). Effects of model-based physiological noise correction on default mode network anti-correlations and correlations. Neuroimage 47, 1448–1459. doi:10.1016/j.neuroimage.2009.05.012

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Delaveau, P., Jabourian, M., Lemogne, C., Guionnet, S., Bergouignan, L., and Fossati, P. (2011). Brain effects of antidepressants in major depression: a meta-analysis of emotional processing studies. J. Affect. Disord. 130, 66–74. doi:10.1016/j.jad.2010.09.032

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Dosenbach, N. U., Fair, D. A., Miezin, F. M., Cohen, A. L., Wenger, K. K., Dosenbach, R. A., et al. (2007). Distinct brain networks for adaptive and stable task control in humans. Proc. Natl. Acad. Sci. U.S.A. 104, 11073–11078. doi:10.1073/pnas.0704320104

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Drevets, W. C., Price, J. L., and Furey, M. L. (2008). Brain structural and functional abnormalities in mood disorders: implications for neurocircuitry models of depression. Brain Struct. Funct. 213, 93–118. doi:10.1007/s00429-008-0189-x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Endicott, J., Spitzer, R. L., Fleiss, J. L., and Cohen, J. (1976). The global assessment scale. A procedure for measuring overall severity of psychiatric disturbance. Arch. Gen. Psychiatry 33, 766–771. doi:10.1001/archpsyc.1976.01770060086012

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Etkin, A., Prater, K. E., Schatzberg, A. F., Menon, V., and Greicius, M. D. (2009). Disrupted amygdalar subregion functional connectivity and evidence of a compensatory network in generalized anxiety disorder. Arch. Gen. Psychiatry 66, 1361–1372. doi:10.1001/archgenpsychiatry.2009.104

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

First, M. B., Spitzer, R. L., Gibbon, M., and Williams, J. B. W. (1996). Structured Clinical Interview for DSM-IV Axis I Disorders: Clinician Version (SCID-CV). Washington, DC: American Psychiatric Press, Inc.

Fox, M. D., and Raichle, M. E. (2007). Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat. Rev. Neurosci. 8, 700–711. doi:10.1038/nrn2201

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fox, M. D., Snyder, A. Z., Vincent, J. L., Corbetta, M., Van Essen, D. C., and Raichle, M. E. (2005). The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc. Natl. Acad. Sci. U.S.A. 102, 9673–9678. doi:10.1073/pnas.0504136102

CrossRef Full Text

Fox, M. D., Zhang, D., Snyder, A. Z., and Raichle, M. E. (2009). The global signal and observed anticorrelated resting state brain networks. J. Neurophysiol. 101, 3270–3283. doi:10.1152/jn.90777.2008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Frodl, T., Jager, M., Smajstrlova, I., Born, C., Bottlender, R., Palladino, T., et al. (2008a). Effect of hippocampal and amygdala volumes on clinical outcomes in major depression: a 3-year prospective magnetic resonance imaging study. J. Psychiatry Neurosci. 33, 423–430.

Pubmed Abstract | Pubmed Full Text

Frodl, T. S., Koutsouleris, N., Bottlender, R., Born, C., Jager, M., Scupin, I., et al. (2008b). Depression-related variation in brain morphology over 3 years: effects of stress? Arch. Gen. Psychiatry 65, 1156–1165. doi:10.1001/archpsyc.65.10.1156

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Frodl, T., Meisenzahl, E. M., Zetzsche, T., Born, C., Jager, M., Groll, C., et al. (2003). Larger amygdala volumes in first depressive episode as compared to recurrent major depression and healthy control subjects. Biol. Psychiatry 53, 338–344. doi:10.1016/S0006-3223(02)01474-9

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fu, C. H., Williams, S. C., Brammer, M. J., Suckling, J., Kim, J., Cleare, A. J., et al. (2007). Neural responses to happy facial expressions in major depression following antidepressant treatment. Am. J. Psychiatry 164, 599–607. doi:10.1176/appi.ajp.164.4.599

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Greicius, M. D., Flores, B. H., Menon, V., Glover, G. H., Solvason, H. B., Kenna, H., et al. (2007). Resting-state functional connectivity in major depression: abnormally increased contributions from subgenual cingulate cortex and thalamus. Biol. Psychiatry 62, 429–437. doi:10.1016/j.biopsych.2006.09.020

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hamilton, J. P., Chen, M. C., and Gotlib, I. H. (2013). Neural systems approaches to understanding major depressive disorder: an intrinsic functional organization perspective. Neurobiol. Dis. 52, 4–11. doi:10.1016/j.nbd.2012.01.015

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hamilton, J. P., Etkin, A., Furman, D. J., Lemus, M. G., Johnson, R. F., and Gotlib, I. H. (2012). Functional neuroimaging of major depressive disorder: a meta-analysis and new integration of base line activation and neural response data. Am. J. Psychiatry 169, 693–703. doi:10.1176/appi.ajp.2012.11071105

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hamilton, J. P., Furman, D. J., Chang, C., Thomason, M. E., Dennis, E., and Gotlib, I. H. (2011). Default-mode and task-positive network activity in major depressive disorder: implications for adaptive and maladaptive rumination. Biol. Psychiatry 70, 327–333. doi:10.1016/j.biopsych.2011.02.003

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hamilton, J. P., and Gotlib, I. H. (2008). Neural substrates of increased memory sensitivity for negative stimuli in major depression. Biol. Psychiatry 63, 1155–1162. doi:10.1016/j.biopsych.2007.12.015

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hamilton, M. (1960). A rating scale for depression. J. Neurol. Neurosurg. Psychiatr. 23, 56–62. doi:10.1136/jnnp.23.1.56

CrossRef Full Text

Heller, A. S., Johnstone, T., Light, S. N., Peterson, M. J., Kolden, G. G., Kalin, N. H., et al. (2013). Relationships between changes in sustained fronto-striatal connectivity and positive affect in major depression resulting from antidepressant treatment. Am. J. Psychiatry 170, 197–206. doi:10.1176/appi.ajp.2012.12010014

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kahn, I., Andrews-Hanna, J. R., Vincent, J. L., Snyder, A. Z., and Buckner, R. L. (2008). Distinct cortical anatomy linked to subregions of the medial temporal lobe revealed by intrinsic functional connectivity. J. Neurophysiol. 100, 129–139. doi:10.1152/jn.00077.2008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kalisch, R. (2009). The functional neuroanatomy of reappraisal: time matters. Neurosci. Biobehav. Rev. 33, 1215–1226. doi:10.1016/j.neubiorev.2009.06.003

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kempton, M. J., Salvador, Z., Munafo, M. R., Geddes, J. R., Simmons, A., Frangou, S., et al. (2011). Structural neuroimaging studies in major depressive disorder. Meta-analysis and comparison with bipolar disorder. Arch. Gen. Psychiatry 68, 675–690. doi:10.1001/archgenpsychiatry.2011.60

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kendler, K. S., Thornton, L. M., and Gardner, C. O. (2001). Genetic risk, number of previous depressive episodes, and stressful life events in predicting onset of major depression. Am. J. Psychiatry 158, 582–586. doi:10.1176/appi.ajp.158.4.582

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kessler, R. C., Berglund, P., Demler, O., Jin, R., Koretz, D., Merikangas, K. R., et al. (2003). The epidemiology of major depressive disorder: results from the National Comorbidity Survey Replication (NCS-R). JAMA 289, 3095–3105. doi:10.1001/jama.289.23.3095

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kronmuller, K. T., Pantel, J., Gotz, B., Kohler, S., Victor, D., Mundt, C., et al. (2008). Life events and hippocampal volume in first-episode major depression. J. Affect. Disord. 110, 241–247. doi:10.1016/j.jad.2008.01.022

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

LaBar, K. S., and Cabeza, R. (2006). Cognitive neuroscience of emotional memory. Nat. Rev. Neurosci. 7, 54–64. doi:10.1038/nrn1825

CrossRef Full Text

LeDoux, J. (2007). The amygdala. Curr. Biol. 17, R868–R874. doi:10.1016/j.cub.2007.08.005

CrossRef Full Text

Liang, Z., King, J., and Zhang, N. (2012). Anticorrelated resting-state functional connectivity in awake rat brain. Neuroimage 59, 1190–1199. doi:10.1016/j.neuroimage.2011.08.009

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

MacQueen, G., and Frodl, T. (2011). The hippocampus in major depression: evidence for the convergence of the bench and bedside in psychiatric research? Mol. Psychiatry 16, 252–264. doi:10.1038/mp.2010.80

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Menon, V. (2011). Large-scale brain networks and psychopathology: a unifying triple network model. Trends Cogn. Sci. (Regul. Ed.) 15, 483–506. doi:10.1016/j.tics.2011.08.003

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Murphy, K., Bodurka, J., and Bandettini, P. A. (2007). How long to scan? The relationship between fMRI temporal signal to noise ratio and necessary scan duration. Neuroimage 34, 565–574. doi:10.1016/j.neuroimage.2006.09.032

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Murray, E. A., Wise, S. P., and Drevets, W. C. (2011). Localization of dysfunction in major depressive disorder: prefrontal cortex and amygdala. Biol. Psychiatry 69, e43–e54. doi:10.1016/j.biopsych.2010.09.041

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ochsner, K. N., and Gross, J. J. (2005). The cognitive control of emotion. Trends Cogn. Sci. (Regul. Ed.) 9, 242–249. doi:10.1016/j.tics.2005.03.010

CrossRef Full Text

Perlman, G., Simmons, A. N., Wu, J., Hahn, K. S., Tapert, S. F., Max, J. E., et al. (2012). Amygdala response and functional connectivity during emotion regulation: a study of 14 depressed adolescents. J. Affect. Disord. 139, 75–84. doi:10.1016/j.jad.2012.01.044

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Rigucci, S., Serafini, G., Pompili, M., Kotzalidis, G. D., and Tatarelli, R. (2010). Anatomical and functional correlates in major depressive disorder: the contribution of neuroimaging studies. World J. Biol. Psychiatry 11, 165–180. doi:10.1080/15622970903131571

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Robinson, O. J., and Sahakian, B. J. (2008). Recurrence in major depressive disorder: a neurocognitive perspective. Psychol. Med. 38, 315–318. doi:10.1017/S0033291707001249

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Roy, A. K., Shehzad, Z., Margulies, D. S., Kelly, A. M., Uddin, L. Q., Gotimer, K., et al. (2009). Functional connectivity of the human amygdala using resting state fMRI. Neuroimage 45, 614–626. doi:10.1016/j.neuroimage.2008.11.030

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Scheidegger, M., Walter, M., Lehmann, M., Metzger, C., Grimm, S., Boeker, H., et al. (2012). Ketamine decreases resting state functional network connectivity in healthy subjects: implications for antidepressant drug action. PLoS ONE 7:e44799. doi:10.1371/journal.pone.0044799

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Seeley, W. W., Menon, V., Schatzberg, A. F., Keller, J., Glover, G. H., Kenna, H., et al. (2007). Dissociable intrinsic connectivity networks for salience processing and executive control. J. Neurosci. 27, 2349–2356. doi:10.1523/JNEUROSCI.5587-06.2007

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Segal, Z. V., Williams, J. M., Teasdale, J. D., and Gemar, M. (1996). A cognitive science perspective on kindling and episode sensitization in recurrent affective disorder. Psychol. Med. 26, 371–380. doi:10.1017/S0033291700034760

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sheline, Y. I., Barch, D. M., Price, J. L., Rundle, M. M., Vaishnavi, S. N., Snyder, A. Z., et al. (2009). The default mode network and self-referential processes in depression. Proc. Natl. Acad. Sci. U.S.A. 106, 1942–1947. doi:10.1073/pnas.0812686106

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sheline, Y. I., Price, J. L., Yan, Z., and Mintun, M. A. (2010). Resting-state functional MRI in depression unmasks increased connectivity between networks via the dorsal nexus. Proc. Natl. Acad. Sci. U.S.A. 107, 11020–11025. doi:10.1073/pnas.1000446107

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sheline, Y. I., Wang, P. W., Gado, M. H., Csernansky, J. G., and Vannier, M. W. (1996). Hippocampal atrophy in recurrent major depression. Proc. Natl. Acad. Sci. U.S.A. 93, 3908–3913. doi:10.1073/pnas.93.9.3908

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sorg, C., Manoliu, A., Neufang, S., Myers, N., Peters, H., Schwerthoffer, D., et al. (2012). Increased intrinsic brain activity in the striatum reflects symptom dimensions in schizophrenia. Schizophr. Bull. 39, 387–395. doi:10.1093/schbul/sbr184

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Squire, L. R., Stark, C. E., and Clark, R. E. (2004). The medial temporal lobe. Annu. Rev. Neurosci. 27, 279–306. doi:10.1146/annurev.neuro.27.070203.144130

CrossRef Full Text

Tang, Y., Kong, L., Wu, F., Womer, F., Jiang, W., Cao, Y., et al. (2012). Decreased functional connectivity between the amygdala and the left ventral prefrontal cortex in treatment-naive patients with major depressive disorder: a resting-state functional magnetic resonance imaging study. Psychol. Med. 43, 1921–1927. doi:10.1017/S0033291712002759

CrossRef Full Text

Van Dijk, K. R., Sabuncu, M. R., and Buckner, R. L. (2012). The influence of head motion on intrinsic functional connectivity MRI. Neuroimage 59, 431–438. doi:10.1016/j.neuroimage.2011.07.044

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Veer, I. M., Beckmann, C. F., van Tol, M. J., Ferrarini, L., Milles, J., Veltman, D. J., et al. (2010). Whole brain resting-state analysis reveals decreased functional connectivity in major depression. Front. Syst. Neurosci. 4:41. doi:10.3389/fnsys.2010.00041

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Winecoff, A., Labar, K. S., Madden, D. J., Cabeza, R., and Huettel, S. A. (2011). Cognitive and neural contributors to emotion regulation in aging. Soc. Cogn. Affect. Neurosci. 6, 165–176. doi:10.1093/scan/nsq030

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wong, C. W., Olafsson, V., Tal, O., and Liu, T. T. (2012). Anti-correlated networks, global signal regression, and the effects of caffeine in resting-state functional MRI. Neuroimage 63, 356–364. doi:10.1016/j.neuroimage.2012.06.035

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wood, K. H., Ver Hoef, L. W., and Knight, D. C. (2012). Neural mechanisms underlying the conditioned diminution of the unconditioned fear response. Neuroimage 60, 787–799. doi:10.1016/j.neuroimage.2011.12.048

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text