Functional Connectivity of Reward Processing in the Brain

Controversial results have been reported concerning the neural mechanisms involved in the processing of rewards and punishments. On the one hand, there is evidence suggesting that monetary gains and losses activate a similar fronto-subcortical network. On the other hand, results of recent studies imply that reward and punishment may engage distinct neural mechanisms. Using functional magnetic resonance imaging (fMRI) we investigated both regional and interregional functional connectivity patterns while participants performed a gambling task featuring unexpectedly high monetary gains and losses. Classical univariate statistical analysis showed that monetary gains and losses activated a similar fronto-striatal-limbic network, in which main activation peaks were observed bilaterally in the ventral striatum. Functional connectivity analysis showed similar responses for gain and loss conditions in the insular cortex, the amygdala, and the hippocampus that correlated with the activity observed in the seed region ventral striatum, with the connectivity to the amygdala appearing more pronounced after losses. Larger functional connectivity was found to the medial orbitofrontal cortex for negative outcomes. The fact that different functional patterns were obtained with both analyses suggests that the brain activations observed in the classical univariate approach identifies the involvement of different functional networks in the current task. These results stress the importance of studying functional connectivity in addition to standard fMRI analysis in reward-related studies.


INTRODUCTION
Bribing someone with cash or, alternatively, threatening him with a high penalty powerfully infl uences behavior. Expectations of both punishment and reward also have an impact on future decisions, since we try to increase the probability of the occurrence of positive reinforcement and to minimize the recurrence of adverse event (Daw et al., 2006).
The delineation of the neural circuits subserving the processing of rewards and punishments and their translation into action is therefore of great importance, but controversial fi ndings have been reported. On the one hand, there is ample evidence indicating that monetary gains and losses activate a similar fronto-subcortical network but to a differential degree (Dreher, 2007;Gottfried et al., 2003;Marco-Pallares et al., 2007;Nieuwenhuis et al., 2005;Tom et al., 2007;van Veen et al., 2004). On the other hand, recent studies suggest that reward and punishment outcomes may be processed by different neural circuits (Frank et al., 2004;Wrase et al., 2007;Yacubian et al., 2006). For example, modeling work of Frank et al. (2004) distinguishes between two excitatory/inhibitory pathways in the basal ganglia, which show differential modulation during positive and negative reinforcement processing. Dopamine (DA) release is typically evoked by positive outcomes, and in turn increases the activity in the excitatory pathway and suppresses the activity in the inhibitory connection. In contrast, negative events are associated with dips in DA levels, and thus show opposite effects. Drawing on neuroimaging results, Yacubian et al. (2006) distinguished two turn might mediate various aspects of behavior. For instance, Cohen et al. (2005) demonstrated an increase of the estimated functional connectivity between the anterior cingulate cortex (ACC) and the nucleus accumbens (NAcc), when comparing high-risk vs. low-risk gambling decisions. Klein et al. (2007) in a feedback-based learning task have shown a dynamically changing functional connectivity pattern between the dorsal ACC, the hippocampal formation, and the ventral striatum that enables the learning of stimulus-reward associations across the task. Cohen et al. (2008) also reported recently that the microstructural properties of white matter tracts connecting the amygdala to the hippocampus, orbitofrontal cortex (OFC), and the ventral striatum predicted functional connectivity derived from functional magnetic resonance imaging (fMRI) time series and participants' behavior following both positive and negative feedback in a reversal learning task.
Which regions might interact in the processing of reward? The ventral striatum has been proposed to be involved in the selection of appropriate responses and the modulation of goal-directed behavior (Berridge and Robinson, 1998;Kelley and Berridge, 2002). This region has shown increased activation in the presence of positive reward outcomes (monetary gains) when compared to negative outcomes (monetary losses) in several studies (Delgado et al., 2000(Delgado et al., , 2003May et al., 2004;Riba et al., 2008). More importantly, studies by Tom et al. (2007) have recently shown that activation in the ventral striatum decreased as the size of the potential loss increased. The ventral striatum receives synaptic inputs from the OFC, dorsolateral prefrontal cortex, and limbic structures such as hippocampus and amygdala (Groenewegen et al., 1999). It is also the target of dense dopaminergic projections originating in the ventral tegmental area (VTA) (Voorn et al., 1986). The VTA has been repeatedly associated to learning and motivation (Schultz, 2007). Therefore the ventral striatum is a key candidate in the study of functional connectivity in the context of reward processing.
With this aim, healthy volunteers were involved in a simple gambling task modifi ed from previous studies (Gehring and Willoughby, 2002;Marco-Pallares et al., 2008;Riba et al., 2008). In this task unexpected and large monetary gains and losses (henceforth boost trials) occurred infrequently in addition to frequent gain and loss trials of smaller magnitudes. We used the ''beta series correlation'' method proposed by  to examine event-related changes in whole-brain functional connectivity with the ventral striatum as a seed region and to compare gain and loss functional connectivity patterns. If two regions are functionally interacting in response to a specifi c event, the fl uctuation of blood oxygen level dependent (BOLD) activity of both regions across trials should be correlated. We predict that the role of reward and loss outcomes might be functionally differentiated by examining patterns of brain connectivity. To our knowledge, no other study has directly compared functional connectivity between reward and loss outcomes.

PARTICIPANTS
Seventeen young adult students [10 women, 21.6 ± 2.6 (SD) mean age] from the University of Barcelona participated in the study. All participants were healthy, right-handed native Spanish speakers with no history of neurological or psychiatric episodes. They all gave written informed consent to a protocol approved by the University of Barcelona ethics committee.

TASK DESIGN
Several important modifi cations were made to the monetary gambling task designed by Gehring and Willoughby (2002). In the standard trials (80%) a warning signal was presented ("*"; 500 ms duration) followed by the presentation of two numbers (5 and 25) displayed in white against a black background in the two possible combinations, [5 25] or [25 5]. Participants had to select one of the two numbers by pressing a spatially corresponding button with the left or right index fi nger. One second after the choice, one of the numbers turned green and the other turned red. If the number selected by the participant changed to red, the participant incurred a loss of the corresponding amount of money in Euro cent. In contrast, if the number turned into green, this indicated a gain.
In addition to the standard trials described above, two additional conditions were created in order to assess brain responses to unexpected rewards and losses. In 10 % of the trials ("boost unexpected trials"), an unexpectedly large gain or loss occurred. Independently of the chosen item (either 5 or 25) the feedback turned into 125 (125 € cents), again having either green or red color indicating wins or losses. To control for the fact that boost trials were both, large and unexpected, in an additional 10% of the trials ("similar unexpected") the chosen number turned to either 7 (instead of 5) or 27 (instead of 25). While these trials were unexpected, the magnitude of the gain or loss was virtually unchanged. The four possible outcomes for the standard trials ([25 5] [5 25] [5 25] [25 5]), and the respective versions for the unexpected similar trials and for the unexpected boost trials were presented in random order. Additionally, each run included 12 randomized fi xation trials that lasted 20 s.
Participants were provided with an initial 10 € sum and were encouraged to win as much as possible. They had performed in a similar event-related potential task (without unexpected trials) several weeks to months earlier and thus were familiar with the task in general. They were informed about the potential occurrence of unexpected trials. The experiment comprised four blocks of 140 trials each. These combinations were counterbalanced by condition, making the statistically expected outcome 0 on each trial in order to avoid confounds of differential probability of gains or losses. At the end of each run, participants were informed about the accumulated amount of money. At the end of the experiment, participants were paid the fi nal amount obtained.
One of the attributes that boosts the impact of reward is the degree of uncertainty that exists in the estimation of the action's value. Reward information becomes more relevant as the uncertainty of the reward prediction and the magnitude of the outcome increases (Fiorillo et al., 2003;Schultz et al., 1997). Indeed, two recent fMRI studies have provided evidence that the NAcc response is larger to less likely rewards (Abler et al., 2007;Yacubian et al., 2006). Moreover, dynamic changes in functional connectivity patterns have been shown over the time course of an experiment indicating learning of stimulus-reward associations (Klein et al., 2007). We therefore assumed that the reward-related activations would be more pronounced and constant across the experiment for the unexpected boost trials which is why we primarily focus on this class of trials. To check whether the results obtained for the boost trials can be generalized to the standard condition [Gain/Loss (5/25)] or instead refl ect only surprise, functional and connectivity gain and loss-related patterns were assessed for standard trials as well. To compare these trials, we created four analogous functional contrasts: Gain (125) vs. fi xation, Gain (5 + 25) vs. fi xation, Loss (125) vs. fi xation and Loss (5 + 25) vs. fi xation. Second, brain regions responding selectively to gains and losses were defi ned by the gain vs. loss contrast in the boost and the standard trials refl ecting the effect of valence in the expected and unexpected condition.

fMRI ACQUISITION
fMRI data was collected using a 3T whole-body MRI scanner (Siemens Magnetom Trio, Erlangen, Germany). Visual images were back-projected onto a screen using an LED-projector and participants viewed the images through a mirror on the head coil. Magnet compatible response buttons were used. Conventional high-resolution structural images [magnetization-prepared, rapidacquired gradient echoes sequence, 192 slice sagittal, TR = 2500 ms, TE = 4.77 ms, TI = 1100 ms, fl ip angle = 7°, 1 mm thickness (isotropic voxels)] were followed by functional images sensitive to BOLD contrast (echo planar T2*-weighted gradient echo sequence, TR = 2000 ms, TE = 30 ms, fl ip angle = 80°). Each functional run consisted of 336 sequential whole-brain volumes comprising 32 axial slices aligned to the plane intersecting the anterior and posterior commissures, 3.5 mm in-plane resolution, 4 mm thickness, no gap, positioned to cover all but the most superior region of the brain and the cerebellum.

MRI ANALYSIS
Functional images were analyzed using standard procedures implemented in the Statistical Parameter Mapping software (SPM2, http://www.fi l.ion.ucl.ac.uk/spm). First, functional volumes were phase shifted in time with reference to the fi rst slice to minimize purely acquisition-dependent signal-variations across slices. Head-movement artefacts were corrected based on an affi ne rigid body transformation, where the reference volume was the fi rst image of the fi rst run (e.g., Friston et al., 1996). Functional data was then spatially smoothed with an 8 mm full-width half-maximum (FWHM) isotropic Gaussian Kernel.
For the group-level analyses, realigned functional data was averaged and the mean functional image was normalized to a standard stereotactic space using the EPI derived MNI template (ICBM 152, Montreal Neurological Institute) provided by SPM2. After an initial 12-parameter affi ne transformation, an iterative non-linear normalization was applied using discrete cosine basis functions by which brain warps are expanded in SPM2 (Ashburner and Friston, 1999). Resulting normalization parameters derived for the mean image were applied to the whole functional set. Finally, functional EPI volumes were resampled into 4 mm cubic voxels and then spatially smoothed with an 8 mm FWHM isotropic Gaussian Kernel to minimize effects of inter-subject anatomical differences. Notice, that all statistical analyses, with the exception of the group statistics in the functional connectivity analysis (see below) were performed in native space (i.e., without spatial normalization).

Univariate fMRI analysis
The statistical evaluation of our data was based on a least-square estimation using the general linear model by modeling the different conditions with a regressor waveform convolved with a canonical hemodynamic response function (Friston et al., 1998). Thus, an event-related design matrix was created including the conditions of interest: Gain 5, Gain 25, Gain 7/27, Gain 125, Loss 5, Loss 25, Loss 7/27, Loss 125, and fi xation. Both native and normalized data was high-pass fi ltered (to a maximum of 1/90 Hz), and serial autocorrelations were estimated using an autoregressive model [AR(1) model]. Resulting estimates were used for non-sphericity correction during the model estimation. Confounding effects in global mean were removed by proportional scaling, and signal-correlated motion effects were minimized by including the estimated movement parameters.

Region of interest analysis
Our approach for identifying functional networks requires the defi nition of an a priori region of interest (ROI), which is then used to determine which voxels throughout the whole-brain are functionally interacting with the selected ROI. Given the strong evidence implicating the ventral striatum in the processing rewards and losses, functional ROIs were defi ned as the largest cluster of connected voxels which were more strongly activated for gains than losses (P < 0.01 uncorrected, 20 voxels spatial extent, restricted to the ventral striatum). In those participants for whom the activation extended beyond the ventral striatum the threshold was increased to P < 0.001 or P < 0.0001 in order to restrict the seed region to the ventral striatum (Figure 3).

Functional connectivity analysis
Functional connectivity analysis was performed using the method proposed by  using the parameter estimates obtained in the context of the general linear model. The beta series bivariate method uses a standard general linear model approach but adapts the model in such a way that separate beta values (general linear model parameter estimates) are obtained in order to modulate activation changes for each component of each individual trial. With this approach, a series of parameter estimates can be extracted from a seed region, correlated across the brain to identify specifi c networks.
The analysis was based on the hypothesis that if different regions are involved in a network, strongly correlated activity patterns should be observed among them and functional brain connec tions can thus be inferred. In particular, each specifi c trial was modeled as an independent covariate in the study design matrix, which allowed the assessment of a trial-to-trial parameter estimate for each condition. For each participant, the parameter estimates series obtained were sorted by conditions. For the conditions of interest (standard loss, boost loss, standard gain, and boost gain), the corresponding parameter estimates were averaged across all voxels of the left and right NAcc ROI. Thereafter, individual native space correlation maps were generated for each participant and condition by correlating the seed region with the beta series of each voxel in the whole-brain. To allow statistical inferences to be made based on the correlation magnitude, an arc-hyperbolic tangent transform was applied to normalize the correlation coeffi cients, which were then z-transformed. The z-transformed correlation maps were spatially smoothed with an 8 mm FWHM) isotropic Gaussian Kernel to minimize effects of inter-subject anatomical differences. The resulting condition-specifi c connectivity maps were then normalized by applying the corresponding normalization parameters, which had been computed earlier. The individual contrast images were entered into a second-level analysis employing a one-sample t-test using a random effects analysis within the general linear model in order to characterize the networks involved in processing gains and losses.
Maps thresholded at P < 0.00001, uncorrected for multiple comparisons, were used for further discussion. Thereafter, correlation maps were compared between the gain and loss contrast at the group-level applying a paired samples t-test, using a P < 0.001 uncorrected for multiple comparison threshold.
Finally, under the assumption that functional interactions between brain regions should refl ect differential functional brain activations, we investigated how the differential connectivity profi le observed between the gain and loss condition interacted with the functional activation pattern: group-level standard parameter estimates were averaged separately by condition and compared with each other in the statistically signifi cant connectivity cluster. The connectivity cluster applied for testing the expected interaction was defi ned under the differential connectivity pattern, using a P < 0.05 (uncorrected for multiple comparisons) threshold.

UNIVARIATE ANALYSIS FOR GAIN AND LOSS TRIALS
In the present gambling task, we focus our analysis on the unexpected gain and loss boost trials. Standard functional univariate analysis was performed in order to compare the overall pattern of activity for gains and losses. This analysis revealed a very similar fronto-subcortical-parietal network when monetary gains or losses were incurred (see Tables 1 and 2). To summarize the tabulated . Additionally, no signifi cant differences were found in the inverse contrast (loss vs. gain trials), even after lowering the threshold to P < 0.05 uncorrected. Overall, the analysis shows that positive and negative outcomes evoked very similar brain activity (see Figures 1 and 2).

FUNCTIONAL CONNECTIVITY ANALYSIS
Concerning functional connectivity, we used the ventral striatum (NAcc) identifi ed above as a seed region (see Figure 3) to contrast gains vs. losses and to determine which brain regions signifi cantly correlated with activity in the ventral striatum. First, functional connectivity was examined separately for gains and losses and then tested for signifi cant differences in connectivity between both conditions in boost and standard trials. The list of the reported functional connectivity interactions is presented in Tables 3 and 4. Similar connectivity patterns were identifi ed for the boost and the standard condition involving an extensive network of regions including the hippocampus, insular cortex, and OFC, which exhibited activity that correlated signifi cantly with the activity seen in the ventral striatum in gain trials as well as in loss trials (see Figures 4 and 5).
We also investigated which regions showed a signifi cant difference in correlation with the NAcc for gain > loss trials. Whereas this contrast revealed signifi cant differences in connectivity neither for the boost nor for the standard condition, we found signifi cant differences in the inverse contrast (i.e., loss > gain) for the boost trials in the medial OFC (peak activity, −8, 24, −16). To assess the functional consequences of the differential connectivity pattern group-level standard parameter estimates were compared between conditions in the statistically signifi cant connectivity cluster. Indeed, the differential connectivity pattern also revealed a differential functional activation pattern when the gain and loss conditions were directly compared [F (1,16) = 5.3, P < 0.036, cluster size = 44 voxels extent].
A clear increase in BOLD signal was observed for gains, whereas a decrease was seen for losses (see Figure 6).
As visual inspection of Figure 4 clearly suggested a difference in connectivity pattern between gain and losses in particular in the amygdala area we performed a further exploratory ROI analysis. However, no signifi cant differences were observed between both conditions in the right or left amygdala region.

DISCUSSION
Applying a recently developed functional connectivity procedure  in an event-related fMRI experiment featuring monetary rewards and losses, we examined whether the interregional interactions maintained by the ventral striatum (NAcc) can be used to characterize and, possibly, dissociate the processing of gains and losses. Standard univariate fMRI analysis revealed a very similar neural network for processing of gains and losses, except for a larger activation observed in the ventral striatum in gain trials. This pattern agrees with previous fMRI studies that underscored the fundamental role of this neural network in the encoding, updating, and maintenance of rewards and punishments (Breiter et al., 2001;Delgado et al., 2000Delgado et al., , 2003Knutson et al., 2001Knutson et al., , 2003McClure et al., 2004;Riba et al., 2008;Tom et al., 2007;Yacubian et al., 2006).

FIGURE 1 | Univariate fMRI analysis in the boost trials.
Coronal and sagittal views of group average whole-brain univariate functional analysis superimposed on a group-averaged structural MRI image in standard stereotactic space (T-score overlays). Shown are the gain vs. fi xation contrast (A) and the loss vs.
Brain networks processing reward FIGURE 2 | Univariate fMRI analysis in the standard trials. Coronal and sagittal views of group average whole-brain univariate functional analysis superimposed on a group-averaged structural MRI image in standard stereotactic space (T-score overlays). Shown are the gain vs. fi xation contrast (A) and the loss vs. fi xation contrast (B) using different statistical thresholds (P < 0.0001; P < 0.0005; P < 0.001). Gain and loss patterns are simultaneously represented in (C): gain (green, P < 0.001), loss (red, P < 0.001), and conjunction gain പ loss (yellow, P < 0.001 and P < 0.001).

FIGURE 3 | Mean region of interest for the boost and standard condition.
The left and right NAcc ROI was functionally defi ned for each participant in native space by identifying the statistically signifi cant activation cluster in the Gain (5 + 25) vs. Loss (5 + 25) and in the Gain (125) vs. Loss (125) contrasts. Displayed is the mean region of interest after normalizing and superimposing the individual ROIs in the standard (A) and boost (B) condition on a group-averaged structural MRI image in standard stereotactic space. The functional connectivity analysis using the ventral striatum as a seed region revealed a topographically distinct subcortical-limbicanterior prefrontal network when compared to the previous standard fMRI analysis. Whereas the general connectivity patterns for gain and loss trials showed similarities, there were also remarkable differences. First, the correlation between ventral striatum and the OFC was stronger in loss trials than in gain trials for the boost condition. Second, the amygdala cluster that showed signifi cant correlations to the ventral striatum seed region was considerably larger for losses compared to gains. As this result was not signifi cant, it awaits confi rmation before a differential functional connectivity between both regions during the processing of gains and losses can be concluded. In a more general way, the different neural network patterns observed between the standard analysis and the connectivity analysis stress the importance of studying functional connectivity as a complementary tool, as it has been successfully demonstrated in previous studies  (Buchsbaum et al., 2005;Fiebach and Schubotz, 2006;Gazzaley et al., 2004;Ranganath et al., 2005;.

FUNCTIONAL CONNECTIVITY ANALYSIS OF GAINS AND LOSSES
Connectivity analysis revealed a network including the OFC, the insular cortex, the amygdala, and the hippocampus that correlated with the activity observed in the seed region (ventral striatum) in the processing of gains and losses. Following seminal work by Alexander et al. (1986) a number of frontal-basal ganglia circuits dedicated to learning, motivation and motor action selection have been described which modulate cortical processing (Kelley, 2004;Lisman and Grace, 2005;Münte et al., 2008). It has to be pointed out that the functional loops that have been described for these different domains overlap and share some important processing stations. Thus, a more realistic view might be that of a network, as depicted in Figure 7 following earlier work by Kelley (2004). In this fi gure we highlighted those areas of the learning/motivation/motor preparation network that comprised the network processing gains and losses as revealed by our connectivity analysis. In the following we discuss the role of the different players within this network for the processing of gains and losses.

Amygdala and hippocampus
A main result of the experience of positive (monetary gains in our experiment) and negative (losses) events is the modulation of future behavior on the basis of learned stimulus-outcome contingencies. Interestingly, the activation of the substantia nigra (SN)/VTA and the hippocampus has been associated with novelty processing and facilitation of memory formation (Schott et al., 2004(Schott et al., , 2006Wittmann et al., 2008). This recent work in humans dovetails nicely with a proposal by Lisman and Grace (2005) that envisions the hippocampus and the midbrain dopaminergic neurons as parts of a functional loop. Novel information that is not yet stored in long-term memory activates the hippocampus which projects to the subiculum, the NAcc, ventral pallidum and to the midbrain, leading to novelty-dependent fi ring of these cells. This fi ring evokes a release of DA that enhances long-term potentiation in the hippocampus and as a result the system is learning. In light of the Lisman and Grace account the interregional connectivity between the NAcc and the hippocampus probably refl ects the integration of contextual aspects related to reward processing (Adcock et al., 2006;Moore and Price, 1999;Wittmann et al., 2005). As no signifi cant differences were observed in the amygdala's connectivity to the striatum between gain and loss trials, the numerically greater proportion of activated voxels in the loss trials should be interpreted with caution. The amygdala has traditionally been associated with the processing of aversive states (Baxter and Murray, 2001;Yacubian et al., 2006) but others have also proposed a central role for this structure in processing both aversive and pleasant stimuli (Hamann et al., 1999;Phan et al., 2004;Phelps  and LeDoux, 2005;Salinas and White, 1998). It is interesting to note that Yacubian et al. (2006) on the basis of standard activation contrasts distinguished two systems involved in predictions of gains on the one hand and losses on the other hand. In their account the amygdala is involved in the prediction of possible losses whereas the ventral striatum is involved in the prediction of gains. While we examined reward/loss delivery rather than prediction, our results likewise hint at a greater role of the amygdala in the processing of losses than in the processing of gains. These differences appear to be qualitative rather than quantitative in the present task, however. Previous studies have indicated that the OFC-amygdala circuit is critical for the adaptation to changes in stimulus-reward or response-reward mappings (Cools et al., 2004;Goto and Grace, 2005;Kesner and Rogers, 2004;Kringelbach, 2005). Such adaptation processes seem to be more crucial after loss-events than gain-events, because the former indicate the need for a change in behavior.

The insular cortex
Our results suggest that interregional functional connectivity between the insula and the ventral striatum is equally strong during the processing of gains and losses. This disagrees with earlier fi ndings that the insula is primarily involved in the processing of negative events (Morris et al., 1999;Phillips et al., 1998). Within the context of reward processing, the anterior insula has been related to reward-related uncertainty (Contreras et al., 2007;Ernst et al., 2002;Hsu et al., 2005;Paulus et al., 2003;Preuschoff et al., 2008;Reynolds and Zahm, 2005;Sell et al., 2000) and risk aversion (Kuhnen and Knutson, 2005;Preuschoff and Bossaerts, 2007). In fact, in order to determine the correct value of an uncertain outcome, it is necessary to evaluate the risk, which is then used to estimate the probability of a possible reward (Preuschoff and Bossaerts, 2007). The latter function may explain why this region has also been implicated in drug craving by several imaging studies (Contreras et al., 2007;Garavan et al., 2000;Sell et al., 2000;Wang et al., 1999). Anatomically, the insular cortex has bidirectional connections with many structures implicated in reward and decision making, including the OFC, the ACC, the NAcc, and the amygdala (Reynolds and Zahm, 2005). Because of its extensive interconnectivity, the insular cortex might be crucial for integrating emotion-related and interoceptive information and to feed this information forward to the OFC and ACC, thus infl uencing decision making, as well as directly affecting other reward-related limbic regions like the amygdala and NAcc.

The role of the medial OFC in reward processing
The medial OFC was differentially involved in the processing of loss compared to gain outcomes. The role of the OFC in the FIGURE 4 | Connectivity analysis for the boost trials. Coronal and axial views of the group average whole-brain interregional interactions with the ventral striatum superimposed on a group-averaged structural MRI image in standard stereotactic space (T-score overlays). Functional connectivity is examined in the reward (A), and punishment (B) condition using different statistical thresholds (P < 0.0001; P < 0.0005; P < 0.001). Gains and losses connectivity patterns are simultaneously depicted in (C): gain (green, P < 0.001), loss (red, P < 0.001), and conjunction gain പ loss (yellow, P < 0.001 and P < 0.001).
processing of rewards and punishment has been extensively documented (Rolls, 1996(Rolls, , 2000. O'Doherty et al. (2001) reported an increase in the activity of the medial OFC as a function of the magnitude of the reward or loss incurred. They showed a graded increase in the activation of this region in relation to reward, but also a decrease relative to baseline when a punishment was delivered, a pattern similar to the one documented in Figure 6 for the present study. The deactivation observed in the OFC for loss trials might refl ect diminished presynaptic input from NAcc neurons. This idea is supported by reports that NAcc activity is suppressed after reward outcomes (Breiter et al., 2001;Delgado et al., 2000;May et al., 2004) and when anticipated rewards are not obtained (Breiter et al., 2001;Delgado et al., 2000;Knutson et al., 2001Knutson et al., , 2003O'Doherty et al., 2002). Moreover, Yacubian et al. (2006) found higher deactivation in the NAcc during loss trials when the loss condition was less likely to occur and in Tom et al. (2007) activation in both the ventral striatum and the medial OFC decreased as the size of the potential loss increased. This suggests that activation decreases in the NAcc and the medial OFC are related to the impact of a negative outcome. The present study adds the important fi nding of a stronger functional connectivity observed in losses compared to gains between the OFC and the NAcc. Interestingly, the OFC has recently been studied by a combined functional/structural connectivity approach using a feedback-guided reversal learning task . In this study tractography of amygdala-OFC connections predicted an individual's learning behavior following rule reversals, stressing the OFC's role in learning from negative events in concert with the amygdala.

CONCLUSIONS
In sum, using connectivity analysis with the ventral striatum as a seed region a brain network comprising the amygdala, hippocampus, insular cortex, and OFC was revealed for both boost and standard trials. While the network was similar for both, gain and loss trials, both conditions were only functionally differentiated for the boost trials in the medial OFC for which greater functional connectivity to the ventral striatum was shown in loss trials. Interestingly, a more extended connectivity from the NAcc to the amygdala was observed for loss trials and may refl ect that amygdala's role in the implementation of behavioral changes after loss-events.
In light of models mainly based on anatomical and physiological work in animals (c.f. Figure 7) it is interesting to speculate about the information fl ow within the brain network revealed for the processing of positive and negative feedback events. It appears that, FIGURE 5 | Connectivity analysis for the standard trials. Coronal views and axial of the group average whole-brain interregional interactions with the ventral striatum superimposed on a group-averaged structural MRI image in standard stereotactic space (T-score overlays). Functional connectivity is examined in the reward (A), and punishment (B) condition using different statistical thresholds (P < 0.0001; P < 0.0005; P < 0.001). Gains and losses connectivity patterns are simultaneously depicted in (C): gain (green, P < 0.001), loss (red, P < 0.001), and conjunction gain പ loss (yellow, P < 0.001 and P < 0.001).
depending on the particular processing requirements, nodes of the wider motivation/learning/motor action network work together in a fl exible manner, leading to situation-specifi c functional connectivity patterns. Functional connectivity analyses have the unique power to highlight these situation-specifi c patterns. For example, in a recent study Kahnt et al. (2008) showed that the dorsal and ventral striatum are differentially connected to different midbrain regions tentatively identifi ed as the SN and VTA. Moreover, in this study functional connectivity between the dorsal striatum and the substantia was related to the impact of different feedback stimuli on a participant's future behavior. It is important to bear in mind that spatial overlap between the gain and loss-related networks does not rule out the possibility that different neural populations within the same regions are involved in these different conditions. For example, there is evidence for excitatory and inhibitory cells in the ventral striatum that have functionally different properties. Indeed, Frank et al., (2004) have proposed a differential modulation of these populations during positive and negative reinforcement. These differences can not be disentangled by comparing the functional patterns, since excitatory and inhibitory contributions can not be differentiated in terms of the evoked BOLD response (Logothetis, 2008).
The importance of the dorsal ACC and the ventral striatum in the adjustment of behavior based on the reinforcement learning signals emanating from midbrain regions has been stressed by recent theories (e.g., Holroyd and Coles, 2002). Surprisingly, our functional connectivity analyses with the ventral striatum as a seed did not reveal the dorsal ACC, although anatomical connections do exist and previous studies have reported functional connectivity between them Klein et al., 2007).
Finally, we would like to point out that connectivity analyses are inherently correlational and, thus, without a good anatomical and neurocognitive model no strong statements can be made about the direction of information fl ow. Coactivity shows that different regions are functionally related, but further experimental manipulations are needed to address the nature of this functional relationship.

FIGURE 7 | Reward processing network in light of the current results.
Illustration of the reward processing network revealed by the current connectivity analysis (black boxes, black arrows) embedded in a wider motivation/learning network (gray boxes and arrows). The wider network is slightly modifi ed (omitting unspecifi c hypothalamic/thalamic projections) from Kelley (2004).