Abstract
Here, we used diffusion tensor imaging (DTI) and showed that the strength of an axonal pathway identified between the amygdala and prefrontal cortex predicted individual differences in trait anxiety. A functional magnetic resonance imaging (fMRI) functional localizer that has been shown to produce reliable amygdala activation was collected in 20 psychiatrically healthy subjects. Voxelwise regression analyses using this fMRI amygdala reactivity as a regressor were performed on fractional anisotropy images derived from DTI. This analysis identified a white matter pathway between the amygdala and ventromedial prefrontal cortex. Individual differences in the structural integrity of this putative amygdala–prefrontal pathway were inversely correlated with trait anxiety levels (i.e., higher pathway strength predicted lower anxiety). More generally, this study illustrates a strategy for combining fMRI and DTI to identify individual differences in structural pathways that predict behavioral outcomes.
Introduction
Evidence from the nonhuman animal literature indicates that the ventromedial prefrontal cortex (vmPFC), including the medial orbitofrontal cortex and ventral anterior cingulate cortex, is reciprocally connected with the amygdala, both functionally and anatomically (Amaral et al., 1992; Milad and Quirk, 2002; Ghashghaei et al., 2007). In humans, functional magnetic resonance imaging (fMRI) studies have used functional connectivity methods to evaluate the relationship between the amygdala and vmPFC (Pezawas et al., 2005; Hare et al., 2008). These data will be bolstered by evidence in the human of actual structural connections between the amygdala and vmPFC, which can only be inferred from nonhuman animal data.
One method that holds promise for delineating such connections in the human brain is diffusion tensor imaging (DTI). This neuroimaging method measures the strength and direction of water diffusivity (i.e., anisotropy) in brain tissue. A popular scale that measures the degree of anisotropy is fractional anisotropy (FA), which is thought to be modulated by the degree of myelination, axonal membrane thickness and diameter, and/or amount of parallel organization of axons (Basser and Pierpaoli, 1996; Beaulieu, 2002). Thus, FA is interpreted as an indicator of white matter pathway strength, or integrity. Of critical relevance to the present study is a recent DTI study that demonstrated an axonal pathway between the amygdala and vmPFC (Johansen-Berg et al., 2008).
Here, we sought to combine fMRI and DTI by using individual differences in amygdala reactivity as the basis for identifying individual differences in the structural integrity of connections between the amygdala and prefrontal cortex. To this end, we had subjects view fearful versus neutral facial expressions, a task that has been shown to reliably activate the human amygdala. Individual differences in these activations across subjects were regressed against the strength of FA signal values focusing on potential connections with the prefrontal region. Given previous studies demonstrating lower anxiety levels in subjects showing a greater functional coupling between the amygdala and prefrontal cortex (Pezawas et al., 2005), we planned to regress individual differences in the strength of this pathway against the range of trait anxiety scores in our subject sample to determine whether the strength of amygdala–prefrontal structural connectivity inversely predicted individual differences in trait anxiety levels.
Materials and Methods
Subjects.
Twenty healthy volunteers (12 women; 21.1 ± 4.15 years of age; 20 right-handed) were screened for current or past psychiatric illness (axis I or II) using the Structured Clinical Interview for DSM-IV (SCID). No subjects had ever taken psychotropic medications. The study protocol was approved by the Committee for the Protection of Human Subjects at Dartmouth College. Written, informed consent was obtained from the participants before the experiment.
Stimuli.
Faces with fearful and neutral expressions from six different individuals (three males and three females) were used (Ekman and Friesen, 1976). The faces were normalized for size and luminance. All of the stimuli were backprojected onto a screen, on which the subjects viewed using a mirror that was mounted on the head coil.
Procedure.
Subjects were asked to passively view fearful and neutral faces during scanning presented separately in alternating 18 s blocks, interleaved with 18 s blocks showing a single crosshair at the middle of the screen. Within each 18 s face block, a total of 36 fearful or neutral faces were presented on a black background. Each face was shown on the screen for 200 ms with a fixed interstimulus interval of 300 ms, consistent with previous studies (Kim et al., 2003). The order of the faces within a block was pseudorandomized to ensure that the same face was not presented more than twice in a row. The order of face blocks was counterbalanced across subjects. Each scan was 4 min and 14 s long. Immediately after the scanning session, all subjects were asked to rate the valence (scale from 1 to 9: 1, very negative; 3, negative; 5, neither negative nor positive; 7, positive; 9, very positive) and arousal (scale from 1 to 9: 1, very low; 3, low; 5, medium; 7, high; 9, very high) levels of each face that they had seen in the scanner. Valence and arousal ratings from one subject were unavailable because of technical problems. After the scanning sessions, each subject's anxiety level was assessed with the State Trait Anxiety Inventory (STAI-S; STAI-T) (Spielberger et al., 1988) self-report questionnaires. In addition, the subjects completed the Beck Depression Inventory (BDI) (Beck et al., 1961). Handedness was determined with the Edinburgh Handedness Inventory (Oldfield, 1971).
Image acquisition.
All subjects were scanned on a 3.0 tesla Philips Intera Achieva Scanner (Philips Medical Systems) equipped with a SENSE birdcage head coil. Anatomical T1-weighted images were collected using a high-resolution three-dimensional magnetization-prepared rapid gradient echo sequence, with 160 contiguous 1 mm thick sagittal slices [echo time (TE), 4.6 ms; repetition time (TR), 9.8 ms; field of view (FOV), 240 mm; flip angle, 8°; voxel size, 1 × 0.94 × 0.94 mm]. Functional images were acquired using echo-planar T2*-weighted imaging sequence. Each volume consisted of 36 interleaved 3 mm thick slices with 0.5 mm interslice gap (TE, 35 ms; TR, 2000 ms; FOV, 240 mm; flip angle, 90°; voxel size, 3 × 3 × 3.5 mm). Diffusion weighted images were collected using echoplanar imaging with 70 contiguous 2-mm-thick axial slices and 32 noncollinear diffusion gradients (TE, 91 ms; TR, 9013 ms; b value, 1000 s/mm2; FOV, 240 mm; flip angle, 90°; voxel size, 1.875 × 1.875 × 2 mm).
Data analysis.
Anatomical and functional images were processed using Statistical Parametric Mapping software (SPM2) (Wellcome Department of Imaging Neuroscience, London, UK). Amygdala blood oxygen level-dependent (BOLD) signal (β weights) from the significantly activated voxels for the fearful versus neutral contrast were extracted for each subject and used as a predictor variable in subsequent voxelwise regression analysis with FA images (for details, see supplemental Methods, available at www.jneurosci.org as supplemental material). Diffusion-weighted images were preprocessed using the Diffusion Toolbox and FSL software of the Oxford Centre for Functional MRI of the Brain (Smith et al., 2004). Diffusion data were corrected for eddy current and head movement, and then the tensor model was fitted to calculate FA values for each voxel, producing one FA image per subject. To generate spatially normalized individual FA images, each individual image was aligned into Montreal Neurological Institute (MNI)-152 standard space using nonlinear registration (Rueckert et al., 1999). After spatial normalization, FA images were smoothed using a Gaussian kernel of 6 mm full width at half-maximum. Voxels with FA values <0.2 were removed before statistical analyses, and an additional white matter mask derived from the Wake Forest University Pick Atlas (Maldjian et al., 2003) was applied to the images to ensure the results were attributable to white matter.
To generate a statistical map in which local FA values were correlated with amygdala BOLD signals, preprocessed FA images were entered into a general linear model with amygdala BOLD signals from fearful versus neutral contrast as a covariate of interest using SPM2. Given our focus on fibers connecting the amygdala with the vmPFC [Johansen-Berg et al. (2008), their Fig. 4] (we will refer to this as the “amygdala–vmPFC pathway”) and the lack of well defined anatomical region of interest mask for this region, a significance threshold of p < 0.05 corrected for multiple comparisons over the frontal and temporal white matter tissue (∼82,000 mm3, which encompasses the spatial location of the amygdala–vmPFC pathway) was used, as determined by Monte Carlo simulations implemented in AlphaSim within AFNI software (Cox, 1996). For regions about which we had no a priori prediction, a more conservative significance threshold of p < 0.00001 (uncorrected, 34 mm3) was used (Shin et al., 2005). FA values of the significant voxels were extracted from FA images that had been normalized and not smoothed to eliminate the effects of smoothing (Tuch et al., 2005). Extracted FA values from these voxels were used for additional statistical investigations, including correlation analyses with trait anxiety measures (supplemental Methods, available at www.jneurosci.org as supplemental material).
Results
Behavioral data
Mean valence ratings of fearful faces (mean ± SD, 2.91 ± 0.85) were significantly more negative than neutral faces (4.60 ± 0.41; t(18) = −7.14; p < 0.001). Fearful faces (5.19 ± 1.42) were also rated as significantly more arousing than neutral faces (3.38 ± 0.98; t(18) = 6.41; p < 0.001). Descriptive statistics for self-report measures are as follows: STAI-S, 32.85 ± 7.03; STAI-T, 37.15 ± 8.25; BDI, 3.40 ± 3.97. These results show that overall subjects perceived fearful faces as negative and arousing, and all scores for anxiety and depression were within the normal range.
Amygdala reactivity to fearful versus neutral faces
As expected, there were significant BOLD signal increases in the amygdala to fearful faces versus neutral faces. The most prominent signal increases were observed at the lateral border of the left ventral amygdala (MNI −33, 3, −24; t(19) = 3.70; p < 0.05 corrected).
DTI voxelwise regression analysis
Left ventral amygdala responses to fearful versus neutral faces were positively correlated with local FA values signaling the presence of white matter fibers located between the amygdala and the medial prefrontal cortex. Critically, these significantly correlated voxels were located within an amygdala–vmPFC pathway independently identified by Johansen-Berg et al. (2008) using tractography. Significantly correlated voxel clusters extended from the left dorsal amygdala/substantia innominata region (MNI −27, 3, −11; t(18) = 2.14), through the left ventral striatum (MNI −18, 3, −8; t(18) = 2.70), and terminated within the left medial orbitofrontal cortex (MNI −17, 25, −20; t(18) = 4.51) (Fig. 1, Table 1).
The identification of this pathway was specific to its correlation with amygdala responses to fearful faces since a similar correlation was observed to fearful versus the fixation baseline (r = 0.559; p < 0.05) but not to neutral versus fixation (r = −0.383; p = 0.10). No voxels within the amygdala–vmPFC pathway were negatively correlated with amygdala BOLD responses to fearful versus neutral faces. Valence and arousal ratings of the faces did not show significant relationships with FA values. FA values, which range from 0 to 1.0, were >0.3, mitigating concern that these effects were related to any gray matter differences as these are typically observed with FA values <0.2.
There was some evidence of correlation with white matter located within the right hemisphere, but this effect was confined to a right ventral striatum (MNI 25, 6, −12; t(18) = 3.69) region similar to that on the left (Fig. 1, coronal slices, Y = 0–9) and did not show evidence of consistent correlation throughout a pathway to the prefrontal cortex.
Correlations with anxiety
Individual differences in trait anxiety scores were negatively correlated with the mean FA value of the entirety of the identified amygdala–vmPFC pathway (r = −0.515; p < 0.05) (Fig. 2). The magnitude of this effect was strongest in the posterior region (MNI −10, 12, −11; t(18) = 4.12) of the amygdala–vmPFC pathway (supplemental Fig. 1, available at www.jneurosci.org as supplemental material).
To determine the anatomical specificity of this relationship between FA and trait anxiety, an additional whole-brain voxelwise regression analysis was performed on the FA images, using trait anxiety as a regressor. Results showed that the relationship between FA and trait anxiety is specific to white matter fibers within the frontal and temporal areas, primarily overlapping with the originally defined amygdala–vmPFC pathway (supplemental Fig. 2, available at www.jneurosci.org as supplemental material).
Trait anxiety was not correlated with amygdala BOLD signal increases to fearful versus neutral faces (r = −0.222; p = 0.35).
Discussion
In this study, we used individual differences in fMRI responses in the amygdala to identify a white matter pathway between the amygdala and vmPFC using DTI. We then independently demonstrated that individual differences in the structural integrity of this pathway predicted differences in trait anxiety levels.
A compelling aspect of the present data set is that the FA correlation map—which simply shows the spatial locations of voxels in which FA values were associated with amygdala reactivity to fearful faces—was almost exclusively allocated within the probabilistic tractography results that identified a similar amygdala–vmPFC pathway in an independent study (Johansen-Berg et al., 2008). We observed white matter circuitry that connected the lateral ventral amygdala with the dorsal amygdala extending in a superior direction into the substantia innominata region of the ventral basal forebrain, ventral striatum/nucleus accumbens, and vmPFC. These data converge with a previous DTI study of both human and macaque brains that identified comparable patterns of anatomical connectivity across both species, including extensive connections between the amygdala and vmPFC (Croxson et al., 2005). Our data extend these findings by showing that individual differences in the strength of this pathway can be used to predict behavioral outcomes, in this case, trait anxiety.
We are careful to refer to the identified white matter as a “putative” pathway because there are a number of candidate white matter sources for the FA values observed here. First, the more lateral FA voxels could be consistent with the location of the uncinate fasciculus (UF). The UF is the major fiber tract connecting the anterior temporal lobe to the inferior frontal cortex (Ebeling and von Cramon, 1992; Highley et al., 2002). Critically, the UF interconnects the amygdala and the orbitofrontal lobe (Ebeling and von Cramon, 1992). In Figure 1, the most lateral FA voxels, just lateral and ventral to the putamen (see coronal slices, Y = 6–18), are consistent with the location of the UF. But both the present study and the previous findings of Johansen-Berg et al. (2008) also reveal more medial white matter. These could potentially represent shorter, multisynaptic tracts between the amygdala and the ventral basal forebrain (Aggleton et al., 1980) and striatum (Fudge et al., 2002) and then connections from these regions and the orbitofrontal cortex (Ferry et al., 2000). However, they could represent fiber tracts projecting in a lateral direction that serve to connect more medial gray matter structures (e.g., amygdala, nucleus basalis) with the UF [Selden et al. (1998), their Fig. 3C].
Since the anatomical connections between the amygdala and vmPFC are bidirectional (Amaral et al., 1992; Ghashghaei et al., 2007), the present data may suggest that more coherent structural connectivity between the amygdala and vmPFC, implying the possibility of more efficient communication and information exchange between the two brain regions, predicts more favorable outcomes in terms of anxiety. Indeed, a positive correlation between structural and functional connectivity between the amygdala and vmPFC has been demonstrated, such that individuals with stronger anatomical connections between these two brain regions demonstrate stronger functional connectivity (Wang et al., 2009). Furthermore, a recent study showed an inverse relationship between low-expressing serotonin transporter gene (5-HTTLPR) alleles, which have been known to be associated with anxiety or depression, and FA values of the left frontal uncinate fasciculus, which overlaps with the amygdala–vmPFC pathway in our study (Pacheco et al., 2009). Together with the present results, these studies collectively suggest that the coherence of white matter between the amygdala and vmPFC is strongly linked with beneficial functional outcomes in terms of anxiety.
These structural results complement the functional findings of Pezawas et al. (2005), who showed that increased functional coupling between the amygdala and vmPFC predicted lower levels of temperamental anxiety. This study also reported that only the functional connectivity of the amygdala-vmPFC was associated with anxiety—whereas other measures, including amygdala reactivity to fearful faces, similar to our findings, were not. The lack of a strong relationship between amygdala reactivity to fearful faces and anxiety suggests that, in psychiatrically normal subjects, assessing the strength of the connectivity between the amygdala and vmPFC, be it structural or functional, could be a better predictor of anxiety.
That said, it should be noted that some studies have reported a significant positive correlation between amygdala reactivity and measures of anxiety (Bishop et al., 2004; Etkin et al., 2004; Haas et al., 2007; Dickie and Armony, 2008). Interestingly, we note that the majority of these findings have linked amygdala reactivity to unaware or unattended fearful faces and anxiety (Bishop et al., 2004; Etkin et al., 2004; Dickie and Armony, 2008), and thus, the present effect may be related to our use of overtly presented fearful faces.
In the present study, a positive correlation was observed between the strength of this amygdala–vmPFC pathway and amygdala reactivity to fearful faces. Consistent with this finding, Pezawas et al. (2005) demonstrated that amygdala activity was positively correlated with vmPFC activity. These findings might at first seem counterintuitive since numerous studies have shown an inverse relationship between the amygdala and vmPFC (Kim et al., 2003; Shin et al., 2005; Urry et al., 2006; Hare et al., 2008). First, these inverse effects have been observed during tasks in which subjects regulated their responses to presented stimuli. Thus, whether an inverse relationship is observed between the amygdala and the vmPFC may depend heavily on the nature of the task used. Perhaps the current task using passive viewing of fearful faces was the basis of the observed identified pathway and future DTI studies could use paradigms that call for more active regulation for comparison. Second, other studies showing an inverse relationship between the amygdala and vmPFC involved clinically pathological groups (Rauch et al., 2000; Shin et al., 2005; Etkin and Wager, 2007). For example, Shin et al. (2005) observed such an effect during the passive viewing of fearful faces in a group of subjects with posttraumatic stress disorder. Thus, perhaps the exaggerated amygdala responses observed in patient groups call for attempts by the vmPFC to regulate such responses (Shin et al., 2005; Etkin and Wager, 2007). DTI studies of these patient populations combined with these behavioral tests during fMRI could help clarify these issues.
Finally, we note that there are numerous constructs that can be used to measure the propensity for negative affect, or negative temperament. For example, although conceptually similar, trait anxiety as measured by the STAI is not identical with neuroticism (Costa and McCrae, 1997) nor harm avoidance (Cloninger, 1986). Harm avoidance is defined as a tendency to respond intensely to aversive stimuli and as an aversion to punishment, novelty, and nonreward (Cloninger, 1986), whereas neuroticism is a tendency to experience negative emotions and interpret situations as threatening (Costa and McCrae, 1997). To date, brain imaging studies of prefrontal–amygdala interactions have used only one of these available measures within a given study [e.g., STAI (Etkin et al., 2004; present study); harm avoidance (Pezawas et al., 2005); neuroticism (Haas et al., 2007)]. Future studies using multiple measures will be needed to dissect the relationship between these related but ultimately unique measures of negative affect and prefrontal–amygdala connectivity.
In conclusion, our data provide evidence that the strength of an anatomical amygdala–prefrontal pathway predicted lower levels of normal trait anxiety. Given that one of the easiest and most convenient ways to elicit human amygdala activation is to overtly present fearful faces, our data offer a simple strategy for combining fMRI and DTI to identify structural pathways that can then predict behavioral outcomes.
Footnotes
-
This work was supported by National Institute of Mental Health Grant MH080716. We thank Erika J. Ruberry, Alison L. Duffy, and Rebecca A. Loucks for technical assistance.
- Correspondence should be addressed to M. Justin Kim, Department of Psychological and Brain Sciences, Dartmouth College, 6207 Moore Hall, Hanover, NH 03755. justin.m.kim{at}dartmouth.edu