Abstract
Among younger adults, the ability to willfully regulate negative affect, enabling effective responses to stressful experiences, engages regions of prefrontal cortex (PFC) and the amygdala. Because regions of PFC and the amygdala are known to influence the hypothalamic-pituitary-adrenal axis, here we test whether PFC and amygdala responses during emotion regulation predict the diurnal pattern of salivary cortisol secretion. We also test whether PFC and amygdala regions are engaged during emotion regulation in older (62- to 64-year-old) rather than younger individuals. We measured brain activity using functional magnetic resonance imaging as participants regulated (increased or decreased) their affective responses or attended to negative picture stimuli. We also collected saliva samples for 1 week at home for cortisol assay. Consistent with previous work in younger samples, increasing negative affect resulted in ventral lateral, dorsolateral, and dorsomedial regions of PFC and amygdala activation. In contrast to previous work, decreasing negative affect did not produce the predicted robust pattern of higher PFC and lower amygdala activation. Individuals demonstrating the predicted effect (decrease < attend in the amygdala), however, exhibited higher signal in ventromedial prefrontal cortex (VMPFC) for the same contrast. Furthermore, participants displaying higher VMPFC and lower amygdala signal when decreasing compared with the attention control condition evidenced steeper, more normative declines in cortisol over the course of the day. Individual differences yielded the predicted link between brain function while reducing negative affect in the laboratory and diurnal regulation of endocrine activity in the home environment.
Introduction
Responding adaptively to stressful experiences is facilitated by the ability to regulate negative emotion. Previous work in younger samples indicates that deliberately diminishing negative affective responses yields increased activation in lateral and dorsal regions of prefrontal cortex (PFC) and/or decreased activation in the amygdala (Beauregard et al., 2001; Ochsner et al., 2002, 2004; Schaefer et al., 2002; Levesque et al., 2003, 2004; Kalisch et al., 2005; Phan et al., 2005; Harenski and Hamann, 2006; Ohira et al., 2006), a pattern suggesting that PFC projections to the amygdala exert a top-down, inhibitory influence. The absence of direct connections between lateral/dorsal PFC and the amygdala, however, indicates that this apparently inhibitory effect must occur indirectly. One plausible route is via medial PFC (MPFC), which has direct connections with lateral/dorsal PFC and also with the amygdala (Price, 2005), and has been implicated in extinction of conditioned fear (Quirk et al., 2003; Phelps et al., 2004).
The amygdala projects to the hypothalamus, a set of nuclei that induces secretion of cortisol from the adrenal cortex by virtue of a hormonal cascade along the hypothalamic-pituitary-adrenal (HPA) axis (Sullivan et al., 2004). Cortisol is a hormone that stimulates glucose production and is therefore secreted in abundance in situations requiring mobilization of energy. Transient increases in circulating cortisol are superimposed on a normatively steep morning-to-evening diurnal rhythm in which the highest levels of cortisol are observed in the early morning hours, reaching their nadir in the evening (Baum and Grunberg, 1995). Flattened cortisol slopes are associated with stressful life circumstances (Polk et al., 2005; Ranjit et al., 2005) and with larger waist circumference, lower levels of social support, higher levels of perceived stress, and poorer explicit memory performance in healthy individuals (Abercrombie et al., 2004). In women with metastatic breast cancer, flatter slopes are associated with higher mortality rates (Sephton et al., 2000). Synthesizing these findings, one pathway leading to flattened cortisol slopes may be the failure to effectively regulate one’s emotional responses, a possibility we investigate in the current report.
Older individuals were recruited for this effort from an ongoing longitudinal study focused on determining whether regulation-induced changes in brain function predict health and well-being in subsequent longitudinal follow-ups. They used cognitive reappraisal to amplify or attenuate their responses to unpleasant photographs in an event-related functional magnetic resonance imaging (fMRI) protocol during which pupil diameter was recorded to index cognitive resource allocation. In addition, participants collected saliva samples for 1 week at home to characterize each person’s diurnal slope of cortisol secretion. We made four primary predictions: (1) replicating previous studies, deliberate attempts to increase and decrease negative affect experience should produce increased lateral/dorsal PFC activation; (2) whereas decreasing negative affect should produce decreased amygdala activation, increasing negative affect should produce increased amygdala activation; (3) MPFC and amygdala responses should be inversely associated during negative affect reduction; and (4) individuals exhibiting the steepest declines in cortisol from morning to night should also exhibit the biggest changes in PFC and the amygdala, a pattern suggesting adaptive HPA function in the face of effective emotion regulation.
Materials and Methods
Participants
Our sample of n = 19 participants (11 females) was recruited from a larger sample of 84 right-handed respondents of the Wisconsin Longitudinal Study (WLS) who participated in a previous study in our laboratory (Jackson et al., 2003; Rosenkranz et al., 2003; Urry et al., 2004). The WLS is a longstanding effort to comprehensively describe the sociodemographic and psychological trajectories of ∼10,000 men and women who were Wisconsin high school seniors in 1957 (for more information, see http://www.ssc.wisc.edu/~wls). Participants in the WLS, including participants in the present study, will continue to be followed in subsequent years, thus providing the opportunity to determine whether emotion regulation-related brain function at this stage predicts health and well-being later in life. Ages ranged from 62 to 64 years.
Participants reported no history of neurological issues or contraindications for MRI scanning (e.g., implanted cardiac or ferrous metal devices). In addition, using a brief screen of cognitive status, the Neurobehavioral Cognitive Status Examination (Northern California Neurobehavioral Group, Fairfax, CA), participants displayed intact performance relative to norms for healthy individuals ages 60–64 years across scales assessing attention, memory, comprehension, reasoning, and judgment (mean ± SD throughout, 80.05 ± 2.33) (Macaulay et al., 2003). Although psychiatric interviews were not conducted, participants fell well within normal limits on questionnaire measures of depressive [Beck Depression Inventory-II (BDI) (Beck et al., 1996)] (3.09 ± 4.11) and anxious [Trait version of the State-Trait Anxiety Inventory, Form X (STAI) (Spielberger et al., 1983)] (31.41 ± 7.43) symptomatology, the two most common forms of psychopathology. [In geriatric samples, a cutoff score of 8 is recommended for the BDI (Kim et al., 2002), whereas a cutoff of 54–55 (determined using the State form) is recommended for the STAI (Kvaal et al., 2005).] Data from two participants (females) were excluded because of excessive movement during the scan, leaving us a final analytic sample size of n = 17 (9 females) for the MR data, with a mean ± SD age of 62.9 ± 0.4 years. Because of technical problems, pupil diameter data were unavailable for an additional five subjects, thus n = 14 (8 females) for pupil analyses (n = 12 subjects providing both pupil and fMRI data). Informed consent was obtained and the rights of participants were protected according to the procedures set forth by the University of Wisconsin–Madison Health Sciences Institutional Review Board. Participants received $70.00 for their participation in the MRI session.
Stimulus materials
We used a randomized event-related paradigm in which we presented one of two sets of digital color photographs (800 × 600 resolution) counterbalanced across participants in the scanner. The photographs were selected from the International Affective Picture System (IAPS) (Center for the Study of Emotion and Attention, 1999). Across the two sets, 72 negative photos were selected according to the IAPS norms to be both highly unpleasant [1, most unpleasant to 9, most pleasant (2.35 ± 0.57)] and arousing [1, least arousing to 9, most arousing (5.82 ± 0.80)], whereas 24 neutral photos were selected to be neither pleasant nor unpleasant (5.04 ± 0.29) and nonarousing (3.14 ± 0.79). In both sets, the negative pictures were significantly more unpleasant and arousing than the neutral pictures according to the normative ratings (all p values <0.05). Picture presentation order was pseudorandomized with the constraint that no more than five pictures of the same valence were shown consecutively. Stimulus presentation was accomplished using the E-Prime software suite (Psychology Software Tools, Pittsburgh, PA) while visual stimulation was delivered via a fiber-optic goggle system (Avotec, Stuart, FL) interfaced with an eye-tracking system (SensoMotoric Instruments, Needham, MA) to measure pupil dilation during the scan session.
Emotion regulation task
Using a variant of the paradigm validated in our previous published electrophysiological work (Jackson et al., 2000), participants were trained to follow one of three instructions during each picture trial: “increase,” “decrease,” or “attend.” The instruction to increase, presented via pneumatic headphones (Avotec), served as a cue to amplify the intensity of their negative affect (subjects heard the word “enhance”). For this task, participants were trained to either (1) imagine themselves or a loved one experiencing the situation being depicted, or (2) imagine a more extreme outcome than the one depicted (e.g., in response to a picture of a ferocious dog, a participant might imagine that the dog’s leash broke and so the dog is about to bite them). Conversely, the instruction to decrease signified the cue to reduce the intensity of their negative affect, for which they were trained to either (1) view the situation as fake or unreal, or (2) imagine that the situation being depicted had a different outcome than the one suggested (e.g., victims of a car accident survived and healed well) (subjects heard the word “suppress”). Alternatively, on attend trials, participants were instructed to maintain their attention to the picture without changing their negative affective experience (subjects heard the word “maintain”). Presentation of regulation conditions was pseudorandomized such that no more than three trials of the same regulation condition were shown consecutively.
Procedure
Participants completed the laboratory portion of this experiment at the University of Wisconsin–Madison in two sessions separated by 1 d. At the first session, participants underwent a “mock” scan to become acclimated to the sights, sounds, and procedures of the scanning environment and to train them to complete the regulation task using a standard set of instructions presented via an E-Prime stimulus presentation program. Participants were positioned inside the bore of an inactive MRI scanner shell, complete with bed, head coil, response box, and goggles, so they could practice completing a block of 24 trials of the regulation task. Simulated scanner sounds and task instructions were piped through ear-bud headphones. The practice block was succeeded by follow-up queries to determine whether participants were using the strategies as instructed and to shape as necessary.
The actual MRI scan session always occurred from 10:00 A.M. to 12:00 P.M. on the morning after the mock session. The timing was held constant in this manner so that associations between fMRI signal and salivary cortisol could not be attributed to variation in the time of day. This scan session was divided into four blocks of 24 trials. Before the real scan, the experimenter presented the task instructions again to consolidate participant understanding of their duties in the scanner. Anatomical images were collected after the functional runs.
Trial structure
At the beginning of each 17 s trial, a white fixation cross was presented in the center of a black screen for 1 s. This was coupled with a simultaneous tone to ensure visual attention to the impending trial. The fixation cross/tone combination was followed by presentation of a picture for 11 s. The regulation instruction (increase, attend, or decrease) was delivered at 3, 4, 5, or 6 s after picture onset, a design feature that provided time for unregulated affect to be elicited before applying voluntary regulatory strategies. The variation in time of instruction onset allowed us to uniquely partition variance associated with regulation of negative emotion at instruction onset separately from that associated with generation of emotion at picture onset. Participants were instructed to follow the task instruction until cued to rate their success at doing so. Success ratings were provided on a four-point scale that was depicted visually for 5 s (1, no emotion; 2, not at all successful; 3, somewhat successful; and 4, very successful).
Data acquisition and reduction
Functional and anatomical MRI.
Whole-brain functional and anatomical images were acquired using a 3.0 tesla GE SIGNA MRI scanner (GE Medical Systems, Waukesha, WI) with LX software (version VH2), a transmit–receive quadrature birdcage coil, and Nvi (40 mT/m; slew rate, 150 mT · m−1 · ms−1) gradients. Functional images were acquired using a T2*-weighted gradient-echo echo-planar imaging (EPI) pulse sequence with 4 mm sagittal slices [interslice gap, 1 mm; echo time (TE), 30.0 ms; repetition time (TR), 2000.0 ms; flip angle, 60°; field of view (FOV), 240 × 240 mm; 64 × 64 data acquisition matrix]. The number of slices required to achieve whole-brain coverage ranged from 30 to 33. With a TR of 2 s and a 17 s trial length, our effective temporal resolution was 1 s. High-resolution T1-weighted spoiled gradient-echo three-dimensional anatomical images were collected in 124 contiguous 1.2 mm slices (TE, 8.0 ms; TR, 21.0 ms; flip angle, 30°; FOV, 240 × 240 mm; 256 × 256 data acquisition matrix), thus enabling localization of functional effects.
The functional MRI data were processed using the Analysis of Functional NeuroImages (AFNI version 2.40e) suite of image analysis routines. During off-line reconstruction, the first five whole-brain images were discarded from each of four scan runs to ensure steady-state magnetization. Because our interest centered on structures whose ventral locations are subject to susceptibility-related signal dropout, we limited our single-subject time series analyses to voxels meeting a signal threshold of at least ∼10% of the whole-brain maximum signal value for the first time point (after discarding the first five images). The raw EPI data were converted to AFNI format with correction for slice timing differences attributable to interleaved acquisition. Using six rigid body motion parameters, each time point was then registered to the one closest in time to collection of our anatomical images. We applied a high-pass filter (0.02 Hz) to remove slow drift in the signal over time and then concatenated the four scan runs into one single-subject dataset for each participant.
Single-subject general linear models (GLMs) using full deconvolution (B. D. Ward, Deconvolution Analysis of fMRI Time Series Data, http://afni.nimh.nih.gov/afni) were performed to estimate the impulse response function for each of six conditions of interest (16 time points). In these omnibus analyses, stick-function predictors were used to indicate picture onset (negative, neutral), as well as regulation instruction onset (negative–increase, negative–decrease, negative–attend, and neutral–attend). Responses to picture onset versus regulation onset were statistically separable because of the 3–6 s jitter in instruction onset. In addition, to account for variability in stimulus unpleasantness across conditions, we modeled normative ratings of unpleasantness gleaned from the IAPS dataset. Ratings were centered separately for the negative and neutral pictures so as to be orthogonal to the valence effects. [There was a small but reliable difference in stimulus unpleasantness such that, across sets, the negative pictures in the increase (2.12 ± 0.50) and attend conditions (2.20 ± 0.42) were more unpleasant than those in the decrease condition (2.72 ± 0.80). By modeling the normative ratings of unpleasantness (centered separately within the negative and neutral valence categories), we were able to statistically control for any small differences in blood oxygenation level-dependent (BOLD) signal that might be accounted for by unpleasantness rather than regulation instruction.] To account for motion-related change in BOLD signal, we also modeled the six motion parameters obtained during volume registration. Trials for which the pupil diameter data were missing for >50% of the total trial length, likely indicating suboptimal viewing of the picture, were eliminated, as were trials that contained time points in which motion estimates exceeded 1.5 mm.
A truncated area-under-the-curve (AUC) metric was calculated by summing estimates of percentage signal change (100 × β coefficient/baseline) across an 8 s window after delivery of the regulation instruction (5–12 s, inclusive). This 8 s summary statistic maximizes chances of capturing the peak of the hemodynamic response even if the peak of the response varies temporally as a function of brain region and/or condition. Additionally, it eliminates noise associated with variation in the rise and fall of the response because those time points are excluded from the AUC summary statistic. Voxels were resampled to 2 mm3 and transformed into Talairach space. Estimates of AUC for percentage signal change were spatially blurred with a 5 mm fill-width half-maximal Gaussian filter.
To test for main effects of regulation condition, we conducted a voxelwise repeated-measures ANOVA in AFNI. We focused on negative trials only because no regulation instructions were provided during neutral trials. Statistically defined clusters of activation were identified using whole-brain Monte Carlo simulations (Alpha Sim program by AFNI) to achieve a corrected cluster threshold of p < 0.05. Specifically, clusters reaching a contiguous volume of at least 296 mm3 at a voxelwise threshold of p < 0.001 or 696 mm3 at a voxelwise threshold of p < 0.005 were considered significant at p < 0.05. Mean AUC estimates across all voxels in each cluster were extracted for each subject and submitted to repeated-measures GLM analyses and pairwise comparisons in SPSS version 12 (SPSS, Chicago, IL). To test a priori hypotheses about amygdala involvement in regulation, mean AUC estimates across all voxels in two (left and right) amygdala regions of interest (ROIs), defined using the Talairach Atlas provided with AFNI, were extracted for each subject and similarly submitted to follow-up GLM and pairwise analyses.
Salivary cortisol.
Using the Salivette saliva collection device (Sarstedt, Newton, NC), a cotton roll housed within a clear plastic tube, participants collected six saliva samples per day for 7 consecutive days at home during the week immediately after their MRI scan session. On each day, the first sample was collected 0.5 h after awakening. Four additional samples were collected at random times within each of the following 3 h blocks of time: 9:00 A.M. to 12:00 P.M., 12:00 P.M. to 3:00 P.M., 3:00 P.M. to 6:00 P.M., and 6:00 P.M. to 9:00 P.M. The final sample was collected at bedtime. Using this method, we ensured that we could adequately resolve the diurnal pattern while minimizing the chance that recurring daily events would systematically bias our estimates.
Samples were centrifuged for 10 min at 5000 rpm to remove cellular and bacterial debris that are inherent in saliva samples. We then conducted a cortisol enzyme immunoassay (Salimetrics, State College, PA). Assay results were considered acceptable only if the coefficient of variation (CV) for the duplicate measurement of a sample was ≤20%. The mean interassay CV was 7.4%, and the mean intraassay CV was 3.8%. The detection limit for this assay is 0.007 μg/dl. Diurnal slope was calculated by regressing log-transformed cortisol values on time for each subject and day. Estimates of slope were then averaged across days and submitted to correlational and regression analyses to assess associations with the aggregate estimates of PFC and amygdala fMRI signal while regulating negative affect.
Pupil diameter.
Horizontal pupil diameter was acquired continuously (sampling rate, 60 Hz) during the scan session. The pupil dilation data were cleaned and processed using algorithms designed by Siegle, Granholm, and Steinhauer (2002, unpublished Matlab code) with Matlab software (MathWorks, Natick, MA) and adapted in our laboratory (L. L. Greischar, 2003, unpublished Matlab code). Blinks were identified and eliminated using local regression slopes and amplitude thresholds. Missing data points were then estimated using linear interpolation. A five-sample rolling average was calculated to smooth the signal, and slow, irrelevant drifts in the pupil diameter data over the course of the scan session were removed via linear detrending over blocks of 24 trials. Trials requiring >50% interpolation during the picture presentation period as well as trials whose value fell greater than ±4 SD from the within-subjects mean were eliminated. Pupil diameter was aggregated into 0.5 s bins, and then means of the unstandardized pupil diameter residuals (removing variance associated with normative ratings of unpleasantness, the same strategy followed in our fMRI analysis) were computed for each epoch during the picture presentation period. These values were range-corrected within subjects [(current pupil diameter − minimum pupil diameter) ÷ (maximum pupil diameter − minimum pupil diameter)] and then aggregated within conditions across trials. Mean diameter for the 0.5 s picture period immediately before the regulation instruction was then subtracted from the mean pupil diameter during each of 10 0.5-s picture periods after instruction up until picture offset within each condition. Proportional change in diameter was then computed for each of the 10 time points after the instruction [i.e., (post − pre)/pre].
We next computed across-subjects GLM analyses to assess the effect of regulation instructions on pupil diameter Although insufficient power precluded treating time as a factor in one omnibus test, we assessed temporal specificity of the increase, attend, and decrease instructions by separately evaluating activity averaged across time points within each of three post-instruction time bins [early (0.5–2.0 s post-instruction average), middle (2.5–3.5 s post-instruction), and late (4.0–5.0 s post-instruction)].
Behavioral data.
Ratings of perceived success were averaged across trials within each regulation condition on a subject-by-subject basis. They were then submitted to a groupwise repeated-measures GLM.
Results
Overview
We first present the main effects of emotion regulation condition (increase, attend, decrease) on fMRI signal. These analyses address our first two predictions concerning PFC and amygdala involvement in voluntary emotion regulation. To address the third and fourth predictions concerning associations between the amygdala and MPFC, and between regions implicated in emotion regulation and salivary cortisol, we next present correlational and regression analyses. Finally, we present the main effects of emotion regulation condition on pupil diameter and ratings of perceived success.
Amygdalar tracking of negative affect regulation
We hypothesized that voluntary regulation of negative affect would linearly modulate amygdalar activation, with increasing yielding higher activation compared to attending, and attending, in turn, yielding higher activation compared to decreasing. Consistent with this pattern, a multivariate GLM testing the effect of condition (increase, attend, decrease) and hemisphere (left, right) on amygdalar response to negative pictures revealed a significant main effect of condition (F(2,14) = 4.85; p = 0.025) (Fig. 1) and a significant within-subjects condition linear contrast (F(1,15) = 9.63; p = 0.007). Mean AUC when increasing (0.77 ± 2.00) was higher than when both attending (0.05 ± 1.92; p = 0.017) and decreasing (−0.16 ± 2.27; p = 0.007) in response to negative trials. Although mean AUC was higher for attend than decrease trials, this predicted effect was not reliable (p = 0.410). There was no effect of hemisphere (F(1,15) = 0.01; p = 0.919), nor was there a significant hemisphere × condition interaction (F(2,14) = 1.18; p = 0.338).
PFC involvement in regulation of negative affect
As depicted in Figure 2, the voxelwise repeated-measures ANOVA testing the main effect of condition (increase, attend, decrease) on AUC in the BOLD response to negative pictures revealed frontal clusters (corrected cluster p < 0.05 based on Monte Carlo simulations) in the left inferior frontal gyrus (IFG) [Brodmann’s area (BA) 44], which includes portions of Broca’s area, a superior dorsal region of the left medial frontal gyrus (MFG) (BA 6), which includes portions of supplementary motor area, a more anterior dorsal region of the left MFG (BA 10), and the left superior frontal gyrus (SFG) (BA 9), part of the dorsolateral PFC. Mean AUC was computed across all voxels within each of these frontal regions and submitted to repeated-measures GLMs followed by pairwise comparisons. Mean AUC in the left IFG cluster when increasing (1.70 ± 0.96) was higher than when attending (0.86 ± 0.85; p < 0.001) or decreasing (0.95 ± 0.84; p < 0.001) (Fig. 2a). The difference between decreasing and attending in the left IFG was not significant (p = 0.423). The same pattern (increase > both attend and decrease; attend = decrease) was also evident for the left SFG (increase, 0.78 ± 0.79; attend, −0.01 ± 0.84; decrease, −0.08 ± 1.00) (Fig. 2b). A linear pattern of effects, similar to that shown in the amygdala, was found in the anterior dorsal MFG (BA 10) (Fig. 2c). In this region, increasing (0.86 ± 1.52) generated higher mean AUC compared to attending (0.07 ± 1.32), and attending was higher in turn than decreasing (−0.18 ± 1.40). Finally, in the superior dorsal MFG/BA 6, both increasing (2.91 ± 1.66) and decreasing (2.34 ± 1.78) resulted in significantly higher activation relative to attending (1.83 ± 1.88), but increasing was still higher than decreasing in this cluster (p < 0.001) (Fig. 2d).
Regulation effects in other brain regions
Table 1 presents a comprehensive list of findings that emerged in the voxelwise ANOVA, including the four regions described above that are also depicted in Figure 2 (noted in Table 1). The cluster with the greatest spatial extent, spanning multiple anatomical regions, was a left posterior medial cluster that included posterior cingulate gyrus, precuneus, and culmen. Also spanning multiple regions, the cluster with the second largest extent included left middle temporal and occipital gyri and the superior temporal gyrus (BA 19 and 39). Additional areas of activation were noted in the right and left fusiform gyrus, right and left thalamus, right and left caudate, and left hippocampus, parahippocampal gyrus, and insula. For each cluster, mean AUC estimates were calculated across all voxels identified in each of these clusters and submitted to a repeated-measures GLM, and significant results were followed up with pairwise comparisons. These results are organized in Table 1 according to the pattern of differences between averaged responses in each cluster for the three regulation conditions (increase, attend, decrease). With the exception of the two dorsal MFG clusters (anterior and superior) described previously (Table 1, bottom), the effect of regulation condition was the same: increasing resulted in higher mean AUC compared to both attending and decreasing, which were not significantly different from one another. There were no main effects of regulation condition evident in VMPFC.
Inverse associations between amygdala and VMPFC
To determine whether MPFC regions show reciprocal coupling with amygdala activation, particularly when decreasing negative affect as predicted, we first calculated the MR signal difference in the amygdala during instructions to regulate relative to the attend control condition. Negative scores in the amygdala for the decrease–attend contrast (i.e., decrease < attend) indicate more effective attenuation of negative affect. Positive scores in the amygdala for the increase–attend contrast (i.e., increase > attend) indicate more effective amplification of negative affect. We then used this amygdala signal difference score as a groupwise predictor of the same signal difference in the remainder of the brain volume on a voxelwise basis.
Decrease–attend
As shown in Figure 3a, participants who most successfully reduced MR signal in the left amygdala showed the highest activation in both right and left VMPFC (BA 11 and 32) (maximum, t(16) = −5.28 at Talairach coordinates x = 5, y = 37, z = −12; and maximum, t(16) = −4.79 at x = −23, y = 43, z = −10, respectively). Although a similar pattern was observed using the right amygdala as a predictor, it did not survive cluster thresholding.
Increase–attend
No significant associations were evident in parallel analyses using the increase–attend signal difference in the amygdala (left or right) as a predictor of the voxelwise increase–attend signal difference.
Specificity
The associations reported above between the left and right VMPFC and the amygdala appear to be specific to the decrease–attend contrast as opposed to representing a more general effect of engaging in regulatory processes. However, to test this more rigorously, we next performed two hierarchical multiple regression analyses. In the first analysis, we entered increase–attend activation in the left amygdala on the first step and increase–attend activation in the left VMPFC on the second step. (Increase–attend activation for VMPFC was extracted from the clusters identified in the voxelwise decrease–attend regression analysis described in the preceding paragraph.) On the third and final step, we entered decrease–attend activation in the left VMPFC. In the second analysis, we observed the same structure, this time entering right rather than left VMPFC on the second and third steps. Supporting a unique decrease–attend effect, results indicate that, even after removing variance explained by increase–attend activation in the left VMPFC and left amygdala, decrease–attend activation in the left VMPFC still explains 44% of the variance in decrease–attend activation in the left amygdala (ΔF(1,13) = 20.71; p = 0.001). Similarly, the right VMPFC still explains 26% of the variance in decrease–attend activation in the left amygdala (ΔF(1,13) = 11.15; p = 0.005). In both analyses, increase–attend activation in neither the left nor the right VMPFC nor in the left amygdala explains unique variance in decrease–attend activation in the left amygdala (all four p values >0.40).
Does VMPFC mediate the association between dorsal PFC and the amygdala?
The inverse association between the amygdala and VMPFC suggested that VMPFC might be exerting a top-down inhibitory effect on the amygdala. In our next set of analyses, we tested whether VMPFC might also be serving the role of intermediary between the amygdala and other regions of PFC that were implicated in emotion regulation (even those for which their involvement was apparently limited to the increase condition). Interestingly, higher activation in this anterior dorsal MFG (BA 10) region was associated with higher activation in the left and right VMPFC when decreasing relative to attending (r(16) = 0.70, p = 0.002; and r(16) = 0.61, p = 0.009, respectively). In addition, subjects demonstrating higher decrease–attend activation in anterior dorsal MFG (BA 10) showed lower decrease–attend activation in the left amygdala (r(15) = −0.54; p = 0.03). Decrease–attend activations in the left superior dorsal MFG (BA 6), left SFG, and left IFG were uncorrelated with one another, nor were they correlated with left amygdalar activation for decrease–attend (p values >0.10). These results suggest that the only candidate region possibly acting on VMPFC was the anterior dorsal MFG (BA 10) because it is the only one to be associated with both VMPFC and the amygdala. [For the sake of completeness, across-subject correlations between the statistically defined frontal regions of interest identified in our tests for regulation main effects (anterior and superior dorsal MFG, SFG, and IFG) and between those regions and the a priori amygdala ROI are presented in Table 2.]
We next conducted the Sobel test of mediation as described and implemented by Preacher and Hayes (2004). In brief, this test provides a means of calculating and testing the significance of the indirect effect of a predictor variable (BA 10) on a criterion variable (amygdala) by way of a hypothesized mediator (VMPFC). This test supported a mediating role for the VMPFC between the anterior dorsal MFG (BA 10) and the amygdala (Z = −2.86, p = 0.0042 for left VMPFC; and Z = −2.33, p = 0.012 for right VMPFC). Moreover, as a test of specificity, reversing the order (i.e., treating VMPFC as the predictor and BA 10 as the mediator) suggested that BA 10 does not account for the association between VMPFC and the amygdala (Z = 0.48, p = 0.63 for left VMPFC; and Z = −0.30, p = 0.76 for right VMPFC). As suggested by the absence of relevant correlations, the VMPFC did not exhibit a mediating relationship with any of the other frontal regions identified earlier [i.e., SFG, IFG, or superior dorsal MFG (BA 6)] and the amygdala.
Patterns of cortisol secretion
Salivary cortisol, expressed in log-transformed micrograms per deciliters exhibited the expected pattern of declining values over the course of the day, as indicated by negative values of diurnal slope averaged across study days in all participants (−0.17 ± 0.05). Levels of cortisol secretion were significantly higher in the first morning sample, taken 30 min after awakening (−0.35 ± 0.17) compared with the last sample of the day, taken just before bedtime (−1.25 ± 0.20; t(15) = 12.28; p < 0.001). These morning and bedtime samples were not significantly correlated in this sample (r(15) = −0.23; p = 0.40), but together they explained 94% of the variance in diurnal slope (F(2,13) = 107.36; p < 0.001), with higher values in the morning being associated with steeper (more negative) slopes (β = −0.49; p < 0.001) and higher values in the evening being associated with flatter (less negative) slopes (β = 0.73; p < 0.001).
Cortisol secretion and its association with emotion regulatory brain responses
We next examined groupwise correlations between brain responses during emotion regulation and the diurnal slope of salivary cortisol. Regions of interest for these analyses were the left and right amygdala ROIs described previously, and the PFC clusters identified either in the voxelwise regression analysis treating left amygdala as a predictor (left and right VMPFC) and in the voxelwise test of regulation main effects (SFG, IFG, anterior dorsal MFG, and superior dorsal MFG).
Decrease–attend
As depicted in Figure 3b, we found that individuals who showed more activation in the left and right VMPFC when decreasing relative to attending (r = −0.72, n = 16, p = 0.002 (bottom left); and r = −0.64, n = 16, p = 0.007 (bottom middle), respectively) and smaller activation in the left amygdala when decreasing relative to attending (r = 0.72; n = 16; p = 0.002 (bottom right)] exhibited sharper declines in cortisol over the course of the day. Although activation in the right amygdala was also positively associated with cortisol slope, this correlation was not significant. Of the four frontal clusters identified in the test of regulation main effects (SFG, IFG, and anterior and superior dorsal MFG), only the anterior dorsal MFG (BA 10) exhibited any significant associations with cortisol. Interestingly, the pattern of associations mirrored those uncovered between VMPFC and cortisol. Specifically, subjects with higher activation in anterior dorsal MFG when decreasing relative to attending exhibited steeper declines in cortisol over the course of the day (r = −0.57; n = 16; p = 0.02).
Increase–attend
Two significant correlations were detected between increase–attend brain responses for the same regions tested above and salivary cortisol slope. Lower activation in the left amygdala when increasing relative to attending was associated with steeper declines in cortisol over the course of the day (r = 0.58; n = 16; p = 0.017). Conversely, lower increase–attend activation in the right VMPFC was associated with flatter declines in cortisol (r = −0.53; n = 16; p = 0.037). These findings are in the same direction as those reported above for the decrease–attend contrast.
Specificity
Because the direction of associations between cortisol slope and brain responses in the left amygdala and right VMPFC while increasing (relative to the attend control condition) is the same as those between cortisol slope and decreasing negative affect, it is possible that processes in these regions related to regulatory efforts in general are the key to exhibiting adaptive patterns of diurnal cortisol secretion. To address this question, we next conducted two hierarchical multiple regression analyses to determine whether decrease–attend activation in these two regions predicts unique variance in cortisol slope after removing variance associated with increase–attend activation. In the first analysis, we entered increase–attend and decrease–attend activation in the left amygdala as predictors of cortisol slope. Results indicate that decrease–attend activation in the left amygdala uniquely explains 27% of the variance in cortisol slope (ΔF(1,13) = 8.74; p = 0.011), even after removing variance explained by increase–attend activation in the left amygdala, which itself is not significant (β = 0.33; p = 0.12). In the second analysis, we entered increase–attend and decrease–attend activation in the right VMPFC as predictors of cortisol slope. Results indicate that decrease–attend activation in the right VMPFC uniquely explains 20% of the variance in cortisol slope (ΔF(1,13) = 4.94; p = 0.045), even after removing variance explained by increase–attend activation in the right VMPFC, which itself is not significant (β = −0.28; p = 0.24).
Unique contributions of morning versus evening cortisol values
Simultaneous regression analyses, in which we treated the first morning and bedtime cortisol samples as predictors of activation in VMPFC and the amygdala, suggested that bedtime cortisol levels may be more important than morning cortisol levels in associations between these regions and diurnal cortisol slope. Specifically, bedtime (β = −0.58; p = 0.013) but not morning (β = 0.28; p = 0.20) cortisol levels explained unique variance in left VMPFC when decreasing relative to attending to negative affect. A similar pattern of associations was evident for right VMPFC for bedtime (β = −0.45; p = 0.065) and for morning (β = 0.33; p = 0.17). Finally, bedtime (β = 0.60; p = 0.012) but not morning (β = −0.22; p = 0.30) cortisol levels explained unique variance in the left amygdala when decreasing relative to attending. The same direction of associations as noted above for VMPFC was discovered for anterior dorsal MFG, although neither morning (β = 0.32; p = 0.21) nor bedtime (β = −0.38; p = 0.14) cortisol levels reached significance in this region. We did not conduct parallel analyses for the increase–attend contrast because no unique brain–cortisol associations were uncovered above.
Signal quality
The ventral nature of the VMPFC and amygdala regions described above increases the risk of susceptibility-related signal drop-out relative to their dorsal counterparts. Moreover, our sample is older; thus, although the age range is restricted to just a few years in the present sample, one might expect more morphometric variation than would be evident in a younger sample. Signal drop-out and morphometric variation would both be evident in variation across subjects in baseline BOLD signal, in which case it was possible that baseline signal might explain the inverse between-subjects associations between the VMPFC and left amygdalar regions for the decrease–attend signal difference. We adopted a number of strategies to rule out this possibility. For one, as noted in Materials and Methods, we limited our single-subject time series analyses to voxels meeting a signal threshold of at least ∼10% of the whole-brain maximum signal value for the first time point (after discarding the first five images to ensure steady-state magnetization). Second, on a groupwise basis, we extracted estimates of baseline signal for the right and left VMPFC and for the left amygdala and computed zero-order correlations between baseline signal and decrease–attend signal in each region. No significant correlations emerged (p values >0.30). We also used the baseline values as control variables in hierarchical regression analyses. Activation in left and right VMPFC clusters for the decrease–attend difference continued to predict unique variance in left amygdala signal for the decrease–attend comparison even when separately controlling for variation in baseline signal in these regions (β = −0.86, ΔR2 = 0.72, ΔF(1,13) = 35.58, p < 0.001 for right VMPFC; and β = −0.90, ΔR2 = 0.76, ΔF(1,13) = 44.82, p < 0.001 for left VMPFC). Variation in baseline signal as a function of susceptibility artifact and/or volumetric variation in this aged sample does not therefore explain the inverse association between VMPFC and amygdala reported above. However, this does not rule out the possibility that morphometric variation might still be at work because morphometric changes may be manifest in task-related functional activations and not in the baseline state.
Medication usage
All participants reported taking at least one form of medication within 2 weeks of the scan session, not unsurprising given their advanced age. Medications reported include the following: anti-arthritis (n = 5), anti-cholesterols (n = 6), antidepressants (n = 2), anti-hypertensives (n = 5), thyroid hormones (n = 7), estrogen hormones (n = 6), insulin (n = 5), and/or antibiotics (n = 2). This raises the possibility that medication usage may account for some of the observed covariation between VMPFC and the amygdala and/or between these regions and cortisol. However, inspection of the zero-order scatter plots suggested that these subjects did not carry the effects reported. In addition, excluding the individuals reporting use of these medications in the past 2 weeks does not affect the pattern of results. In all cases, the direction of the associations is preserved, and, in most cases, the correlations remain significant or nearly so despite the loss of statistical power.
Regulation-induced changes in pupil diameter
Figure 4 presents proportional change in pupil diameter over time after delivery of increase, attend, and decrease instructions on negative picture trials. We predicted that increasing and decreasing negative affect would result in greater proportional increases in pupil diameter compared to attending. Suggesting sustained processing during the two active forms of regulation, both increasing (0.10 ± 0.072) and decreasing (0.07 ± 0.047), resulted in larger pupil diameter compared to attending (0.03 ± 0.053; p = 0.002 and 0.043, respectively), during the late regulation period (F(2,12) = 7.31; p = 0.008). Moreover, the difference in pupil diameter between increasing and decreasing during the late regulation period only bordered on significant (p = 0.085). In the early (F(2,12) = 6.46; p = 0.012) and middle (F(2,12) = 4.91; p = 0.028) (Fig. 4, gray shading) regulation periods, increasing negative affect (0.10 ± 0.049 and 0.12 ± 0.064, respectively) resulted in larger pupil diameter compared to both attending (0.07 ± 0.043 and 0.06 ± 0.054) and decreasing (0.07 ± 0.037 and 0.08 ± 0.047; p values <0.05). There was no difference in pupil diameter for decreasing versus attending during the early and middle time periods (p values >0.40).
Associations between brain activation and pupil diameter
The fMRI analyses identified several regions of PFC and the amygdala that are more active when participants increase or decrease negative affect. Moreover, having just shown that increasing and decreasing negative affect produce significant changes in pupil diameter, we next examined an exploratory question: to what extent is the role of these regions in affect regulation associated with the degree of effort required by the task? We turned to our pupil diameter data to address this question. First, to index individual differences in the extent to which participants expended effort in the two active regulation conditions controlling for attention to negative information, we computed two pupil diameter difference scores: increase–attend and decrease–attend. A positive number indicates greater effort for increase or decrease relative to the control condition. We then calculated correlations between each of these difference metrics and its isomorphic contrast in brain activation for the frontal and amygdalar clusters identified above, regions of a priori interest given their role in regulation of negative affect in this study. Table 3 summarizes these associations. Notably, subjects showing higher activation in the left IFG and left SFG clusters for the decrease–attend contrast also showed higher decrease–attend pupil diameter difference scores in early, middle, and late regulation periods. Conversely, subjects showing higher activation in the left amygdala for the decrease–attend contrast showed lower scores on the decrease–attend pupil diameter difference in the late regulation period, an effect that nearly surpassed the significance threshold. There were no significant correlations between increase–attend activation in the frontal and amygdalar clusters identified in our main effect analyses and the increase–attend pupil diameter difference scores.
Behavioral data
A multivariate GLM testing the effect of condition (increase, attend, or decrease in response to negative stimuli) on ratings of task success, in which ratings ranged from 1 to 4, revealed no significant effect (p > 0.20; increase, 2.91 ± 0.61; attend, 2.76 ± 0.73; decrease, 2.73 ± 0.57).
Discussion
From this work, we draw three conclusions about the nature of the neural and endocrine correlates of emotion regulation, which we operationalized as the use of cognitive reappraisal to increase and decrease felt emotion in response to unpleasant visual stimuli. First, changes in fMRI signal in the amygdala and regions of PFC when attempting to reduce the experience of negative affect, specifically in an anterior region of BA 10 within the medial frontal gyrus and bilateral VMPFC, are lawfully related to diurnal cortisol secretion in the home environment. Second, experimenter-cued (but participant-generated) reframing of negative visual stimuli produces effects in lateral/dorsal PFC and the amygdala in older adults that are similar to, but less robust than, those reported in some previous studies with younger adults. Third, inverse functional coupling between VMPFC and amygdala regions during negative affect reduction is consistent with the conclusion that VMPFC exerts a top-down inhibitory influence on the activity of the amygdala. Importantly, BA 10 exerts an indirect effect on the amygdala that is mediated by bilateral regions of VMPFC.
Individuals who were good at diminishing negative affect, as evidenced by higher VMPFC and lower amygdala activation, exhibited a more normative (i.e., steeper) decline in salivary cortisol over the course of the day, with bedtime cortisol levels being more important than morning levels in this regard. These findings likely do not reflect regulatory skill or general cognitive task performance because accounting for variance associated with intensifying negative affect, itself not a significant predictor, did not obliterate these effects. Animal studies suggest that chronic stress is associated with elevated glucocorticoid secretion in the evening but not in the morning (Marti et al., 1993; Janssens et al., 1995). Additionally, dysregulation of cortisol secretion is evident in psychopathological conditions such as posttraumatic stress disorder and major depression (Burke et al., 2005; Yehuda et al., 2006). As reviewed by Sullivan and Gratton (1999), MPFC is rich in glucocorticoid receptors and has been implicated in negative feedback inhibition of the HPA axis. Disrupted negative feedback inhibition may produce chronic dysregulation of cortisol secretion (Diorio et al., 1993; McEwen, 1999; Sapolsky, 2003), which itself may contribute to disease states (Sephton et al., 2000). This pattern of findings is consistent with the idea that (1) a circuit involving the VMPFC and the amygdala underlies effective emotion regulation ability, and (2) this circuit mediates HPA disturbances associated with chronic stress and psychopathology. In subsequent waves of data collection within the Wisconsin Longitudinal Study, participants in the current study will complete measures of health and well-being that will enable us to use the present brain data to test whether activity in these circuits also predicts longer-term physical and mental health outcomes.
In previous studies, deliberate attenuation of negative affect yielded activation in dorsolateral, ventral lateral, and dorsomedial (including dorsal anterior cingulate cortex) regions of PFC, a finding that has been interpreted to reflect the execution of cognitive control for regulatory purposes (Ochsner et al., 2004; Phan et al., 2005). Interestingly, our analyses of condition main effects yielded only one region in which decreasing negative affect yielded stronger activation than unregulated attention (i.e., supplementary motor area), and even here higher activation was observed when increasing compared with decreasing negative affect. The most striking difference between this study and those reported previously is the advanced age of our sample. Advanced age is associated with changes in brain structure and/or function as well as with differences in cognitive and emotional processing. Aging may therefore be an important factor in understanding the paucity of frontal main effects with decreasing negative affect in this report. The inconsistency with previous work might also reflect the mixed gender nature of our sample and differences in design and methodology. Space considerations preclude an extended discussion of these issues here, but interested readers may refer to the supplemental material (available at www.jneurosci.org).
Whereas we find few PFC effects for the decrease condition, correlational analyses indicate that subjects who allocated more cognitive resources to decreasing than attending to negative affect, as indexed by larger pupil diameters, exhibited more activation in the SFG and IFG and less activation in the amygdala when decreasing than attending. Thus, whereas the involvement of these regions in negative affect reduction was not supported in tests of main effects, lawful associations with cognitive processing demands support the effectiveness of the manipulation. Similarly, lawful associations were also identified between the amygdala and VMPFC when participants deliberately diminished their affective response. Relative to their responses when passively experiencing negative affect, individuals who effectively reduced amygdalar responses when decreasing negative affect exhibited activation in bilateral VMPFC, an effect that was not evident anywhere else in PFC nor when increasing negative affect. Ventromedial PFC has been implicated in human and animal studies of conditioned fear extinction (Milad and Quirk, 2002; Phelps et al., 2004). In addition, Kim et al. (2003) demonstrated that VMPFC and amygdala responses are inversely correlated during processing of surprised facial expressions. Importantly, stimulation of the rodent analog of VMPFC reduces responses in the central nucleus of the amygdala (Quirk et al., 2003).
Based on the literature reviewed above, it should not be surprising that the VMPFC was implicated in the ability of subjects to transform their experience of negative stimuli to be less negative. Notably, the VMPFC regulatory effect is not observed as a mean difference; rather, it was only observed in analyses that harnessed variation across individuals. Whereas significant mean differences between conditions are less likely in the face of such variability, we discovered that those who are good at attenuating negative affect (i.e., they display less activation in the amygdala; decrease < attend) are the same ones exhibiting greater bilateral activation in VMPFC (decrease > attend). Because we control for emotional responses that occur spontaneously by subtracting out responses in the attend condition, this inverse association reflects an emotion regulatory as opposed to reactivity effect in this particular experimental context.
Interestingly, VMPFC appeared to serve a bridging function between the anterior dorsal MFG (BA 10) and the amygdala. Recent theories of BA 10 function suggest that this region is involved in coordination of information processing and information transfer when multiple processes are engaged in the service of a behavioral goal (Ramnani and Owen, 2004) and in coordinating attention between external stimuli and internal thoughts (Burgess et al., 2005). These functions match those required to regulate negative affect in our emotion regulation task, a task that likely involves multiple streams of information processing that would require such coordination, including maintaining representations of the behavioral goal to decrease negative affect, monitoring and modifying spontaneous appraisals arising in response to the images, and evaluating behavioral success. Although promising, the triangulation of associations between BA 10 and the amygdala, with VMPFC playing a mediating role, was not predicted on an a priori basis; thus, cautious interpretation is warranted. Moreover, although the pattern described above is consistent with a top-down process, the correlational nature of these associations and the bidirectional anatomical connectivity necessarily precludes such causal inference.
Conclusion
Regions of prefrontal cortex and the amygdala are responsible for the ability to regulate negative affect, findings that are bolstered here by the use of pupil diameter to verify the sustained active nature of the prescribed cognitive reappraisal strategies (increase, decrease) compared with an attention control condition. These data further indicate that there are individual differences in the ability to regulate, particularly in the ability to attenuate, negative affect as it unfolds over time. Individuals who show bigger reductions in the response of the amygdala when decreasing negative affect compared with the attention control are the same individuals who evidence greater responses in bilateral VMPFC for decreasing negative affect in response to aversive pictures compared with the attention control condition. Furthermore, linking the responses of the brain in a laboratory setting to the exigencies of daily life, higher VMPFC and lower amygdala signal during the downregulation of negative affect relative to the attention control condition, are associated with a more adaptive diurnal rhythm of circulating free cortisol secretion over the course of 1 week at home. It remains an issue for future studies to determine whether these associations generalize to other age groups because only older individuals were studied here. The neural mechanisms of emotion regulation may be different in older versus younger individuals. In summary, the pattern of findings herein suggests that functional coupling between the PFC and the amygdala enables effective regulation of negative emotion. Furthermore, the functioning of PFC–amygdala circuitry during regulation of negative affect predicts longer-term regulation of endocrine activity that may be important for health and well-being.
Footnotes
-
This work was supported by National Institutes of Mental Health Grants P50-MH61083, P50-MH52354, P50-MH069315, and R37-MH43454, National Institute on Aging Grant P01-AG021079, and General Clinical Research Centers Program of the National Center for Research Resources, National Institutes of Health Grant M01-RR03186. We gratefully acknowledge the following individuals: Daren C. Jackson and Alexander J. Shackman for helpful consultation, Michael Anderle and Ron Fisher for assistance with collection of the imaging data, Kristin Rae Evans for conducting the salivary cortisol assays, Joseph A. Hampel and Ann K. Chamberlain for assistance with cortisol preprocessing, and two anonymous reviewers for their very insightful comments.
- Correspondence should be addressed to Heather L. Urry at her present address: Department of Psychology, Tufts University, 490 Boston Avenue, Medford, MA 02155. Email: heather.urry{at}tufts.edu