Abnormal Hippocampal Functioning and Impaired Spatial Navigation in Depressed Individuals: Evidence From Whole-Head Magnetoencephalography

Depressed patients were recruited from inpatient and outpatient psychiatric units at the National Institute of Mental Health (NIMH) as part of randomized clinical medication trials for unipolar depression. The patient group consisted of 19 right-handed individuals with a diagnosis of major depressive disorder, currently depressed without psychotic features; diagnosis was confirmed with the Structured Clinical Interview for Axis I DSM-IV Disorders–Patient Version (21). Additional clinical characteristics are presented in Table 1. Patients with a DSM-IV (22) diagnosis of drug or alcohol dependence or abuse within the past 3 months, history of head injury, unstable medical illness, or uncorrected hypo- or hyperthyroidism were excluded. All subjects had been free of psychotropic drugs for at least 2 weeks (5 weeks for fluoxetine) and were medically healthy as determined by medical history, physical and neurological examinations, and laboratory tests. Most patients (89%) had in their lifetime received trials of medications, the most common being selective serotonin reuptake inhibitors (SSRIs) (79%), venlafaxine (42%), duloxetine (37%), and bupropion (37%). During the current major depressive episode the most common trials were SSRIs (63%), duloxetine (26%), and bupropion (26%). No patient had received electroconvulsive therapy. The severity of their current depressive episodes was assessed by using the Montgomery-Åsberg Depression Rating Scale (23). Nine of the 19 patients had a comorbid anxiety disorder, none of which included posttraumatic stress disorder.

Nineteen right-handed healthy comparison subjects were recruited from the local community. Healthy subjects underwent screening at NIMH that included a medical history, physical examination performed by a physician, and a Structured Clinical Interview for DSM-IV (24). Exclusion criteria were current medical illness, axis I DSM-IV psychiatric disorder, or first-degree relative with a history of psychiatric disorders; family history was assessed per self-report during the interview with the clinician. No subject could be currently taking psychotropic medications per self-report. The study was approved by the Combined Neuroscience Institutional Review Board of NIH. After complete description of the study to participants, written informed consent was obtained.

Participants completed 40 trials in which they moved with a fiber-optic joystick around a virtual pool using commercial software (NeuroInvestigations, Inc., Lethbridge, Alb., Canada). Participants began each trial at the pool's edge in one of four locations (pseudorandomly ordered). For half of the trials, participants were instructed to quickly navigate to a hidden platform in a fixed location (hidden condition). If the hidden platform was not found within 20 seconds, it became visible and participants were instructed to navigate to it to finish the trial. This condition requires the use of distal cues on the surrounding walls to learn the fixed location of the hidden platform. For the other half of the trials, participants were instructed to navigate to a visible platform in a variable location (visible condition). As the control condition, navigating to the visible platform does not require use of distal cues. The platform had a variable location across visible trials to limit spatial learning during this condition. Different pool environments were used for the two conditions with the only difference being the distal cues (e.g., door and shelves in one environment versus abstract patterns in the other); each virtual environment was assigned to contain the hidden or visible platform in equal proportions between groups. There were 2-second breaks between trials. The order of conditions was counterbalanced across participants.

Latency to reach the platform from first movement, latency to first movement, path length, and heading error or angular deviation from a direct path to the platform were recorded for each hidden and visible platform trial. Heading error is a point estimate calculated when the path length exceeds one fourth of the pool's diameter. Each of these measures was collapsed into five blocks of four trials. All participants received exposure to the task in a previous session before the MEG session to ensure ample practice. Platform locations were different between sessions.

Neuromagnetic data were recorded with a CTF-OMEGA 275-channel whole-head magnetometer (VSM MedTech, Coquitlam, B.C., Canada) in a magnetically shielded room (Vacuumschmelze GmbH, Hanau, Germany) using synthetic third gradient balancing for active noise cancellation. Magnetic flux density was sampled at 1200 Hz with a bandwidth of 0–300 Hz. Electrical coils attached at three fiducial sites—nasion and right and left preauricular sites—were energized during each MEG recording to track head movements in real time and for offline registration to anatomical MRIs. Radiological markers (IZI Medical Products, Baltimore, Md.) were attached to the same fiducial sites during the MRI scan. High-resolution T1-weighted anatomical MRIs were obtained in a separate session using a 1.5 or 3 Tesla whole-body scanner (GE Signa, Waukesha, Wisc.).

Synthetic aperture magnetometry is a minimum-variance beamformer algorithm that uses the signal covariance from the sensor array for constructing optimum spatial filters at each location in the brain; an optimum spatial filter suppresses all source activity except for that which originates from the location of interest (25). Regional oscillatory changes estimated by this method are highly consistent with fMRI results (26), at cortical and subcortical loci (27), and in healthy individuals and depressed patients (28). Source imaging procedures were identical to those used by Cornwell et al. (18). Signal covariance was calculated between all sensors with data filtered between 4 and 8 Hz for 20 unaveraged 1-second epochs during navigation for each condition. Source space was sampled in 5-mm steps with current sources estimated at each location based on a multisphere model derived from the MRI of each participant. Dual-state images were calculated by dividing integrated power over 1-second windows during navigation by integrated power over 1-second pretrial windows (pseudo-F ratios), yielding estimates of relative increases (or decreases) in theta activity during navigation in each condition. With a constant pretrial window, relative theta power was estimated across 17 1-second windows during navigation with 75% overlap (0.00–1.00 second, 0.25–1.25, …, 4.00–5.00), encompassing the first 5 seconds of each trial before the platform had been reached.

Using Analysis of Functional NeuroImages software (29), individual source image volumes for each condition and time window were manually warped to a Talairach brain template to allow for group analysis in a standardized source space and within-volume normalized using a z-score transformation. The latter procedure was performed to limit the influence of global power that may vary across participants, or systematically between groups, on local power estimates in a particular brain region.

Navigation performance was assessed as follows: 1) latency to reach the platform from the first recorded movement, 2) path length to the platform, and 3) heading error to the platform. Latency to first movement was also analyzed as a control measure. Gender and age were added as nuisance factors in our analyses given evidence of their associations with spatial navigation performance (30, 31). We performed 2×2×2×5 mixed factorial ANCOVAs, with age as a continuous covariate, group and gender as between-group factors, and condition (hidden versus visible platform) and block (1 to 5) as repeated-measures variables. Significant interactions between group and condition were followed up with simple effect analyses for each condition. For cases in which differences emerged between patients and healthy subjects, post hoc analyses were run to determine whether there was evidence of differences within the patient group between those with and those without comorbid anxiety disorders. Two patients' path length and heading error data were lost to computer malfunction.

Group SAM analyses of 4–8 Hz theta power consisted of 2 (patients versus comparison subjects) × 2 (hidden versus visible) mixed factorial ANOVAs. We confined our subsequent simple effect analyses to functional regions of interest showing significant 2-way interactions at a nominal p value of <0.05. Within these regions, we compared theta activity between groups for the hidden condition and the visible condition using Student's t tests and estimated false discovery rates (32). A single correction for multiple comparisons was performed that incorporated all voxels within these regions and all time windows, such that we kept the overall false discovery rate of the source imaging results below 5%. As with the behavioral data, post hoc analyses were run that compared depressed patients with and those without comorbid anxiety disorders specifically within clusters in medial temporal cortices showing differential activation between groups.

We determined what regional activity correlated at the individual level with navigation performance on hidden trials. Spearman correlations were calculated at each voxel for each time window for each group separately. A nonparametric approach was deemed most appropriate given the small samples. We identified only those regions with at least two contiguous voxels showing negative correlations, commonly in both groups but irrespective of the time window, surviving a statistical threshold of p<0.005. Given previous results in left posterior hippocampal and parahippocampal cortices (18), we considered negative correlations only insofar as they reflect higher regional theta activity being associated with better navigation performance (e.g., shorter mean latencies).

There were significant group-by-condition interactions on mean latency to reach the platform (F=6.89, df=1, 33, p<0.05) and path length (F=6.06, df=1, 31, p<0.05) (Figure 1) but not heading error nor latency to first movement. For the hidden platform condition, patients were slower than comparison subjects to reach the platform (F=10.77, df=1, 33, p<0.005) and they also took longer paths relative to comparison subjects (F=13.28, df=1, 31, p<0.005). For the visible platform condition, there was no difference between patients and healthy subjects in mean latency, but patients did take longer paths to the platform relative to comparison subjects (F=15.76, df=1, 31, p<0.001). For mean heading error, there was a significant group difference (F=8.86, df=1, 31, p<0.01), with patients performing worse than comparison subjects in both conditions. There was also a significant group difference in latency to first movement (F=4.70, df=1, 33, p<0.05), with patients initiating movements faster than comparison subjects in both conditions. Post hoc comparisons between patients with and without comorbid anxiety disorders in the hidden platform condition revealed no differences in latency or path length. Regarding gender differences, men showed shorter path lengths and smaller heading errors than women (F=16.52, df=1, 31, p<0.001 and F=7.69, df=1, 31, p<0.01, respectively). There was a difference between men and women in latency to reach the platform that approached significance (F=3.44, df=1, 33, p<0.10) but not in latency to first movement.

Table 2 presents regions that showed significant group differences in theta activity for the hidden and visible platform conditions. Healthy subjects exhibited greater theta activity than depressed patients in medial temporal cortices in five time windows during the hidden condition, all in the right anterior region (Figure 2). Peak activity differences varied between the hippocampus at 1.25–2.25 seconds and 3.25–4.25 seconds relative to trial onset, and anterior parahippocampal cortices at 0.75–1.75, 2.75–3.75, and 3.00–4.00 seconds. In all cases, group differences reflected a combination of greater theta activity in healthy subjects and lesser activity in patients relative to pretrial baseline activity. There were no instances in which patients exhibited greater theta activity than comparison subjects in medial temporal cortices during the hidden condition.

Other regions showing greater theta activity in healthy subjects relative to patients included multiple loci in lateral prefrontal cortices. Conversely, greater theta activity in the cerebellum and posterior cortical regions was observed in patients relative to comparison subjects. Post hoc analyses found no significant differences in theta activity between patients with and without comorbid anxiety disorders within regions showing differences between depressed and healthy subjects when corrected for multiple comparisons. Finally, few regions exhibited group differences in activity for the visible condition. Healthy subjects exhibited greater activity in dorsal prefrontal cortical regions and, in one instance, left parahippocampal cortices relative to patients. Patients exhibited greater theta activity than comparison subjects in the cerebellum and posterior parietal cortices.

Correlation analyses demonstrated that left posterior hippocampal/parahippocampal theta activity was negatively correlated with mean latency to reach the hidden platform (i.e., greater theta was associated with better performance) in both patients and comparison subjects (Figure 3). Time windows in which these correlations were observed differed between groups, with patients showing this correlation at 3.00–4.00 seconds and healthy subjects at 1.25–2.25 seconds relative to trial onset. Other regions with at least two contiguous voxels showing negative correlations of a similar magnitude commonly across groups in any time window included the following: left precuneus (Brodmann's area [BA] 7; local maxima [x,y,z]= –4, –49, 42), right precentral gyrus (BA 4; x,y,z=21, –23, 59), right precuneus (BA 7; x,y,z=3, –37, 45) and right middle frontal gyrus (BA 10; x,y,z=36, 43, 22).