Neural interactions in the human frontal cortex dissociate reward and punishment learning

Abstract
eLife assessment
Introduction
Results
Discussion
Methods
Data availability
References
Article and author information
Metrics

Abstract

How human prefrontal and insular regions interact while maximizing rewards and minimizing punishments is unknown. Capitalizing on human intracranial recordings, we demonstrate that the functional specificity toward reward or punishment learning is better disentangled by interactions compared to local representations. Prefrontal and insular cortices display non-selective neural populations to rewards and punishments. Non-selective responses, however, give rise to context-specific interareal interactions. We identify a reward subsystem with redundant interactions between the orbitofrontal and ventromedial prefrontal cortices, with a driving role of the latter. In addition, we find a punishment subsystem with redundant interactions between the insular and dorsolateral cortices, with a driving role of the insula. Finally, switching between reward and punishment learning is mediated by synergistic interactions between the two subsystems. These results provide a unifying explanation of distributed cortical representations and interactions supporting reward and punishment learning.

eLife assessment

This is an important information-theoretic re-analysis of human intracranial recordings during reward and punishment learning. It provides convincing evidence that reward and punishment learning is represented in overlapping regions of the brain while relying on specific inter-regional interactions. This preprint will be interesting to researchers in systems and cognitive neuroscience.

https://doi.org/10.7554/eLife.92938.3.sa0

Introduction

Reward and punishment learning are two key facets of human and animal behavior, because they grant successful adaptation to changes in the environment and avoidance of potential harm. These learning abilities are formalized by the law of effect (Thorndike, 1898; Bouton, 2007) and they pertain the goal-directed system, which supports the acquisition of action-outcome contingencies and the selection of actions according to expected outcomes, as well as current goal and motivational state (Dickinson and Balleine, 1994; Balleine and Dickinson, 1998; Balleine and O’Doherty, 2010; Dolan and Dayan, 2013; Balleine, 2019).

At the neural level, the first hypothesis suggests that these abilities are supported by distinct frontal areas (Pessiglione and Delgado, 2015; Palminteri and Pessiglione, 2017). Indeed, an anatomical dissociation between neural correlates of reward and punishment prediction error (PE) signals has been observed. PE signals are formalized by associative models (Rescorla et al., 1972) and reinforcement learning theory (Sutton and Barto, 2018) as the difference between actual and expected action outcomes. Reward prediction error (RPE) signals have been observed in the midbrain, ventral striatum and ventromedial prefrontal cortex (vmPFC) (Schultz et al., 1997; O’Doherty et al., 2004; O’Doherty et al., 2001; Pessiglione et al., 2006; D’Ardenne et al., 2008; Steinberg et al., 2013; Palminteri et al., 2015; Gueguen et al., 2021). Punishment prediction error (PPE) signals have been found in the anterior insula (aINS), dorsolateral prefrontal cortex (dlPFC), lateral orbitofrontal cortex (lOFC), and amygdala (O’Doherty et al., 2001; Seymour et al., 2005; Pessiglione et al., 2006; Yacubian et al., 2006; Gueguen et al., 2021). Evidence from pharmacological manipulations and lesion studies also indicates that reward and punishment learning can be selectively affected (Frank et al., 2004; Bódi et al., 2009; Palminteri et al., 2009; Palminteri et al., 2012). Complementary evidence, however, suggests that reward and punishment learning may instead share common neural substrates. Indeed, hubs of the reward circuit, such as the midbrain dopamine systems and vmPFC, contain neural populations encoding also punishments (Tom et al., 2007; Matsumoto and Hikosaka, 2009; Plassmann et al., 2010; Monosov and Hikosaka, 2012). Taken together, it is still unclear whether reward and punishment learning recruit complementary cortical circuits and whether differential interactions between frontal regions support the encoding of RPE and PPE.

To address this issue, we repose on recent literature proposing that learning reflects a network phenomenon emerging from neural interactions distributed over cortical-subcortical circuits (Bassett and Mattar, 2017; Hunt and Hayden, 2017; Averbeck and Murray, 2020; Averbeck and O’Doherty, 2022). Indeed, cognitive functions emerge from the dynamic coordination over large-scale and hierarchically organized networks (Varela et al., 2001; Bressler and Menon, 2010; Reid et al., 2019; Panzeri et al., 2022; Thiebaut de Schotten and Forkel, 2022; Miller et al., 2024; Noble et al., 2024) and accumulating evidence supports that information about task variables is widely distributed across brain circuits, rather than anatomically localized (Parras et al., 2017; Saleem et al., 2018; Steinmetz et al., 2019; Urai et al., 2022; Voitov and Mrsic-Flogel, 2022).

Accordingly, we investigated whether reward and punishment learning arise from complementary cortico-cortical functional interactions, defined as statistical relationships between the activity of different cortical regions (Panzeri et al., 2022), within and/or between brain regions of the frontal cortex. In particular, we investigated whether reward and punishment prediction errors are encoded by redundancy- and/or synergy-dominated functional interactions in the frontal cortex. The search for synergy- and redundancy-dominated interactions is motivated by recent hypotheses suggesting that a trade-off between redundancy for robust sensory and motor functions and synergistic interaction may be important for flexible higher cognition (Luppi et al., 2024). On one hand, we reasoned that redundancy-dominated brain networks may be associated with neural interactions subserving similar functions. Redundant interactions may appear in collective states dominated by oscillatory synchronization (Engel et al., 2001; Varela et al., 2001; Buzsáki and Draguhn, 2004; Fries, 2015) or resonance phenomena (Vinck et al., 2023). Such collective states may give rise to selective patterns of information flow (Buehlmann and Deco, 2010; Kirst et al., 2016; Battaglia and Brovelli, 2020). On the other, synergy-dominated brain networks may be associated with functionally-complementary interactions. Indeed, synergistic interactions have been reported between distant transmodal regions during high-level cognition (Luppi et al., 2022) and, at the microscale, in populations of neurons within a cortical column of the visual cortex and across areas of the visuomotor network (Nigam et al., 2019; Varley et al., 2023). The notion of redundant and synergistic interactions resonates with the hypothesis that brain interactions regulate segregation and integration processes to support cognitive functions (Wang et al., 2021; Deco et al., 2015; Sporns, 2013; Finc et al., 2020; Cohen and D’Esposito, 2016; Braun et al., 2015; Shine et al., 2016).

In order to study redundancy- and synergy-dominated interactions, we used formal definitions from Partial Information Decomposition (PID; Williams and Beer, 2010; Wibral et al., 2017; Lizier et al., 2018). The PID decomposes the total information that a set of source variables (i.e. pairs of brain signals) encodes about a specific target variable (i.e. prediction errors) into components representing shared (redundant) encoding between the variables, unique encoding by some of the variables, or synergistic encoding in the combination of different variables. Within this framework, we used a metric known as interaction information (McGill, 1954; Ince et al., 2017), which quantifies whether a three-variable interaction (i.e. pairs of brain regions and the PE variable) is either synergy- or redundancy-dominated. We predicted that redundancy-dominated functional interactions would engage areas with similar functional properties (e.g. those encoding RPE), whereas synergy-dominated relations would be observed between areas performing complementary functions (e.g. the encoding of RPE and PPE).

We investigated neural interactions within and between four cortical regions, namely the aINS, dlPFC, lOFC, and vmPFC, by means of intracerebral EEG (iEEG) data collected from epileptic patients while performing a reinforcement learning task (Gueguen et al., 2021). We found various proportions of intracranial recordings encoding uniquely RPE or PPE signals or both, suggesting a local mixed representation of PEs. We then identified two distinct learning-related subsystems dominated by redundant interactions. A first subsystem with RPE-only interactions between the vmPFC and lOFC, and a second subsystem with PPE-only interactions between the aINS and dlPFC. Within each redundant-dominated subsystem, we demonstrated differential patterns of directional interactions, with the vmPFC and aINS playing a driving role in the reward and punishment learning circuits, respectively. Finally, these two subsystems interacted during the encoding of PE signals irrespectively of the context (reward or punishment), through synergistic collaboration between the dlPFC and vmPFC. We concluded that the functional specificity toward reward or punishment learning is better disentangled by interactions compared to local representations. Overall, our results provide a unifying explanation of distributed cortical representations and interactions supporting reward and punishment learning.

Results

iEEG data, behavioral task, and computational modeling

We analyzed iEEG data from sixteen pharmacoresistant epileptic patients implanted with intracranial electrodes (Gueguen et al., 2021). A total of 248 iEEG bipolar derivations located in the aINS, dlPFC, vmPFC, and lOFC regions (Figure 1A) and 1788 pairs of iEEG signals, both within and across brain regions (Figure 1B) were selected for further analysis. Single subject anatomical repartition is shown in Figure 1—figure supplement 1. Participants performed a probabilistic instrumental learning task and had to choose between two cues to either maximize monetary gains (for reward cues) or minimize monetary losses (for punishment cues) (Figure 1C). Overall, they selected more monetary gains and avoided monetary losses but the task structure was designed so that the number of trials was balanced between reward and punishment conditions (Figure 1D).

Figure 1 with 2 supplements see all

Download asset Open asset

intracerebral EEG (iEEG) implantation, behavioral task, and computational modeling.

(A) Anatomical location of intracerebral electrodes across the 16 epileptic patients. Anterior insula (aINS, n=75), dorsolateral prefrontal cortex (dlPFC, n=70), lateral orbitofrontal cortex (lOFC, n=59), ventromedial prefrontal cortex (vmPFC, n=44), (B) Number of pairwise connectivity links (i.e. within patients) within and across regions, (C) Example of a typical trial in the reward (top) and punishment (bottom) conditions. Participants had to select one abstract visual cue among the two presented on each side of a central visual fixation cross and subsequently observed the outcome. Duration is given in milliseconds, (D) Number of trials where participants received outcomes +1€ (142±44, mean ± std) vs. 0€ (93±33) in the rewarding condition (blue) and outcomes 0€ (141±42) to –1€ (93±27) in the punishment condition (red), (E) Across participants trial-wise reward prediction error (PE) (Reward prediction error, RPE - blue) and punishment PE (PPE - red), ±95% confidence interval.

We estimated trial-wise prediction errors by fitting a Q-learning model to behavioral data. Fitting the model consisted in adjusting the constant parameters to maximize the likelihood of observed choices. We used three constant parameters: (i) the learning rate α accounting for how fast participants learned new pairs of cues; (ii) the choice temperature β to model different levels of exploration and exploitation; (iii) Ө parameter to account for the tendency to repeat the choice made in the previous trial. The RPE and PPE were obtained by taking the PE for rewarding and punishing pairs of cues, respectively. RPE and PPE showed high absolute values early during learning and tended toward zero as participants learned to predict the outcome (Figure 1E). The convergence toward zero of RPE and PPE was stable at the single subject level (Figure 1—figure supplement 2).

Local mixed encoding of PE signals

At the neural level, we first investigated local correlates of prediction error signals by studying whether RPEs and PPEs are differentially encoded in prefrontal and insular regions. To this end, we performed model-based information theoretical analyses of iEEG gamma activities by computing the mutual information (MI) between the across-trials modulations in RPE or PPE signals and the gamma band power in the aINS, dlPFC, lOFC and vmPFC. The MI allowed us to detect both linear and non-linear relationships between the gamma activity and the PE. Preliminary spectrally-resolved analyses showed that the frequency range significantly encoding prediction errors was between 50 and 100 Hz (Figure 2—figure supplement 1). We thus extracted for each trial time-resolved gamma power within the 50–100 Hz range using a multi-taper approach for further analyses. MI analysis between gamma power and prediction error signals displayed significant group-level effects in all four cortical regions (Figure 2A) and globally reproduced previous findings based on general linear model analyses (Gueguen et al., 2021). Interestingly, we observed a clear spatial dissociation between reward and punishment PE signaling. Whereas the vmPFC and dlPFC displayed complementary functional preferences for RPE and PPE, respectively, the aINS and the lOFC carried similar amounts of information about both R/PPE (Figure 2A).

Figure 2 with 2 supplements see all

Download asset Open asset

Local mixed encoding of reward and punishment prediction error signals.

(A) Time-courses of mutual information (MI in bits) estimated between the gamma power and the reward (blue) and punishment (red) prediction error (PE) signals. The solid line and the shaded area represent the mean and SEM of the across-contacts MI. Significant clusters of MI at the group level are plotted with horizontal bold lines (p<0.05, cluster-based correction, non-parametric randomization across epochs), (B) Instantaneous proportions of task-irrelevant (gray) and task-relevant bipolar derivations presenting a significant relation with either the reward prediction error (RPE) (blue), the punishment prediction error (PPE) (red) or with both RPE and PPE (purple). Data is aligned to the outcome presentation (vertical line at 0 s).

To better characterize the spatial granularity of PE encoding, we further studied the specificity of individual brain regions by categorizing bipolar derivations as either: (i) RPE-specific; (ii) PPE-specific; (iii) PE-unspecific responding to both R/PPE; (iv) PE-irrelevant (i.e. non-significant ones) (Figure 2B). All regions displayed a local mixed encoding of prediction errors with temporal dynamics peaking around 500 ms after outcome presentation. The vmPFC and dlPFC differentially responded to reward and punishment PEs, and contained approximately 30% of RPE- and PPE-specific contacts, respectively. In both regions, the proportion of RPE- and PPE-specific bipolar derivations was elevated for approximately 1 s after outcome presentation. The lOFC also contained a large proportion of PPE-specific bipolar derivations, but displayed more transient dynamics lasting approximately 0.5 s. The aINS had similar proportions of bipolar derivations specific for the RPE and PPE (20%), with temporal dynamics lasting approximately 0.75 s. Importantly, all regions contained approximately 10% of PE-unspecific bipolar derivations that responded to both RPE and PPE, especially in the aINS and dlPFC. The remaining bipolar derivations were categorized as PE-irrelevant. A complementary analysis, conducted to evaluate inter-subject reproducibility, revealed that local encoding in the lOFC and vmPFC was represented in 30 to 50% of the subjects. In contrast, this encoding was found in 50 to 100% of the subjects in the aINS and dlPFC (Figure 2—figure supplement 2).

Taken together, our results demonstrate that reward and avoidance learning are not supported by highly selective brain activations, but rather from a mixed or mixed encoding of RPE or PPE signals distributed over the prefrontal and insular cortices. Nevertheless, such distributed encoding seems to involve two complementary systems primarily centered over the vmPC and dlPFC, respectively.

Encoding of PE signals occurs with redundancy-dominated subsystems

To better understand the observed complex encoding of reward and punishment PEs, we tested the hypothesis that functional dissociations occur with differential and distributed interactions between prefrontal and insular cortices. To address this question, we performed model-based network-level analyses based on the PID framework (Williams and Beer, 2010; Wibral et al., 2017; Lizier et al., 2018). We particularly used the interaction information (McGill, 1954; Ince et al., 2017) to quantify whether a three-variable interaction (i.e. pairs of brain regions, and the PE variable) is either synergy- and redundancy-dominated (Williams and Beer, 2010). Indeed, interaction information (II) can be either positive or negative. A negative value indicates a net redundancy (i.e. a pair of recordings are carrying similar information about the PE), whereas a positive value indicates a net synergistic effect (i.e. a pair of recordings are carrying complementary information about the PE). We computed the time-resolved II across trials between the gamma activity of pairs of iEEG signals and PEs. To differentiate cortico-cortical interactions for reward and punishment learning, we first calculated the II separately for RPEs and PPEs. RPE- and PPE-specific analyses exclusively showed negative modulations of II, therefore, indicating the presence of redundancy-dominated local and long-range interactions (Figure 3).

Figure 3 with 1 supplement see all

Download asset Open asset

Encoding of prediction error (PE) signals occurs with redundancy-dominated subsystems.

Dynamic interaction information (II in bits) within- (A) and between-regions (B) about the RPE (II_RPE) and PPE (II_PPE) are plotted in blue and red. Significant clusters of II_RPE and II_PPE are displayed with horizontal bold blue and red lines (p<0.05, cluster-based correction, non-parametric randomization across epochs). Significant differences between II_RPE and II_PPE are displayed in green. Shaded areas represent the SEM. The vertical gray line at 0 s represents the outcome presentation.

To better characterize the local interactions encoding reward and punishment PEs, we computed the II between pairs of gamma band signals recorded within the aINS, dlPFC, lOFC, and vmPFC. Within-region II analyses showed that significant RPE-specific interactions were exclusively observed in the vmPFC and lOFC, whereas PPE-specific interactions were present only in the dlPFC. In addition, the aINS was found to display both RPE- and PPE-specific interactions (Figure 3A). A relevant sign of high specificity for either reward or punishment PE signals was the presence of a significant cluster dissociating RPE and PPE in the vmPFC and dlPFC only (green clusters in Figure 3A).

To investigate the nature of long-range interactions, we next computed the II for RPE and PPE between signals from different brain regions (Figure 3B). Similarly, results exclusively showed redundancy-dominated interactions (i.e. negative modulations). RPE-specific interactions were observed between the lOFC and vmPFC, whereas PPE-specific interactions were observed between the aINS and dlPFC and to a smaller extent between the dlPFC and lOFC, peaking at 500 ms after outcome presentation. A significant difference between RPE and PPE was exclusively observed in the lOFC-vmPFC and aINS-dlPFC interactions, but not between dlPFC and lOFC (green clusters in Figure 3B). The analysis of inter-subject reproducibility revealed that both within-area and across-area significant redundant interactions were carried by 30 to 60% of the subjects (Figure 3—figure supplement 1). Taken together, we conclude that the encoding of RPE and PPE signals occurs with redundancy-dominated subsystems that differentially engage prefronto-insular regions.

Contextual directional interactions within redundant subsystems

Previous analyses of II are blind to the direction of information flows. To address this issue, we estimated the transfer entropy (TE) (Schreiber, 2000) on the gamma power during the rewarding (TE_Rew) and punishment conditions (TE_Pun), between all possible pairs of contacts. As a reminder, the TE is an information-theoretic measure that quantifies the degree of directed statistical dependence or ‘information flow’ between time series, as defined by the Wiener-Granger principle Wiener, 1956; Granger, 1969. Delay-specific analyses of TE showed that a maximum delay of information transfer between pairs of signals comprised an interval between 116 and 236 ms (Figure 4—figure supplement 1). We thus computed the TE for all pairs of brain regions within this range of delays and detected temporal clusters where the TE significantly differed between conditions (TE_Rew >TE_Pun or TE_Pun >TE_Rew). Only two pairs of brain regions displayed statistically-significant modulations in TE (Figure 4). We observed that the TE from the aINS to the dlPFC (TE_aINS→dlPFC) peaked at approximately 400 ms after outcome onset and was significantly stronger during the punishment condition compared to the rewarding condition. By contrast, the information flow around ~800 ms from the vmPFC to the lOFC (TE_vmPFC→lOFC) was significantly stronger during the rewarding condition. No other brain interactions were found significant (Figure 4—figure supplement 2). Overall, these results demonstrate that the two redundancy-dominated RPE- and PPE-specific networks (Figure 3B) are characterized by differential directional interactions. The vmPFC and aINS act as drivers in the two systems, whereas the dlPFC and lOFC play the role of receivers, thus suggesting a flow of PE-specific information within the network.

Figure 4 with 2 supplements see all

Download asset Open asset

Contextual modulation of information transfer.

Time courses of transfer entropy (TE, in bits) from the anterior insula (aINS) to the dorsolateral prefrontal cortex (dlPFC) (aINS→dlPFC) and from the vmPFC to the lateral orbitofrontal cortex (lOFC) (vmPFC→lOFC), estimated during the rewarding condition (TE_Rew in blue) and punishing condition (TE_Pun in red). Significant differences (p<0.05, cluster-based correction, non-parametric randomization across epochs) of TE between conditions are displayed with horizontal bold lines (blue for TE_Rew >TE_Pun and red for TE_Pun >TE_Rew). Shaded areas represent the SEM. The vertical gray line at 0 s represents the outcome presentation.

Integration of PE signals occurs with synergy-dominated interactions between segregated sub-systems

Since learning required participants to concurrently explore rewarding and punishment outcomes, we finally investigated the nature of cortico-cortical interactions encoding both RPE and PPE signals. We estimated the II about the full PEs, i.e., the information carried by co-modulation of gamma power between all pairs of contacts about PE signals (Figure 5—figure supplement 1). Encoding of PEs was specifically associated with significantly positive II between the dlPFC and vmPFC (II_dlPFC-vmPFC Figure 5A). Such between-regions synergy-dominated interaction occurred approximately between 250 and 600 ms after outcome onset.

Figure 5 with 3 supplements see all

Download asset Open asset

Synergistic interactions about the full prediction error (PE) signals between recordings of the dlPFC and vmPFC.

(A) Dynamic interaction information (II in bits) between the dorsolateral prefrontal cortex (dlPFC) and vmPFC about the full prediction error (II_dlPFC-vmPFC). Hot and cold colors indicate synergy- and redundancy-dominated II about the full PE. Significant clusters of II are displayed with a horizontal bold green line (p<0.05, cluster-based correction, non-parametric randomization across epochs). Shaded areas represent the SEM. The vertical gray line at 0 s represents the outcome presentation. (B) Dynamic II_dlPFC-vmPFC binned according to the local specificity PPE-RPE (II_PPE-RPE in pink) or mixed (II_Mixed in purple) (C) Distributions of the mean of the II_PPE-RPE and II_Mixed for each pair of recordings (II_PPE-RPE: one-sample t-test against 0; dof = 34; P fdr-corrected=0.015*; T=2.86; CI(95%)=[6.5e-5, 3.9e-4]; II_Mixed: dof = 33; P fdr-corrected=0.015*; T=2.84; CI(95%)=[5.4e-5, 3.3e-4]).

We then investigated if the synergy between the dlPFC and vmPFC encoding global PEs could be explained by their respective local specificity. Indeed, we previously reported larger proportions of recordings encoding the PPE in the dlPFC and the RPE in the vmPFC (Figure 2B). Therefore, it is possible that the positive II_dlPFC-vmPFC could be mainly due to complementary roles where the dlPFC brings information about the PPE only and the vmPFC brings information to the RPE only. To test this possibility, we computed the II_dlPFC-vmPFC for groups of bipolar derivations with different local specificities. As a reminder, bipolar derivations were previously categorized as RPE or PPE specific if their gamma activity were modulated according to the RPE only, to the PPE only, or to both (Figure 2B). We obtained four categories of II. The first two categories, named II_RPE-RPE and II_PPE-PPE, reflect the II estimated between RPE- and PPE- bipolar derivations from the dlPFC and vmPFC. The third category (II_PPE-RPE) refers to the II estimated between PPE-specific bipolar recordings from the dlPFC and RPE-specific bipolar recordings from the vmPFC. Finally, the fourth category, named II_Mixed, includes the remaining possibilities (i.e. RPE-Both, PPE-Both, and Both-Both) (Figure 5—figure supplement 2). Interestingly, we found significant synergistic interactions between recordings with mixed specificity i.e., II_PPE-RPE and II_Mixed between 250 and 600ms after outcome onset (Figure 5B and C). Consequently, the II_dlPFC-vmPFC is partly explained by the dlPFC and vmPFC carrying PPE- and RPE-specific information (II_PPE-RPE) together with interactions between non-specific recordings (II_Mixed). In addition, we simulated data to demonstrate that synergistic interactions can emerge between regions with the same local specificity (Figure 5—figure supplement 3). Taken together, the integration of the global PE signals occurred with a synergistic interaction between recordings with mixed specificity from the dlPFC and vmPFC.

Discussion

Our study revealed the presence of specific functional interactions between prefrontal and insular cortices about reward and punishment prediction error signals. We first provided evidence for a mixed encoding of reward and punishment prediction error signals in each cortical region. We then identified a first subsystem specifically encoding RPEs with emerging redundancy-dominated interactions within and between the vmPFC and lOFC, with a driving role of the vmPFC. A second subsystem specifically encoding PPEs occurred with redundancy-dominated interactions within and between the aINS and dlPFC, with a driving role of the aINS. Switching between the encoding of reward and punishment PEs involved a synergy-dominated interaction between these two systems mediated by interactions between the dlPFC and vmPFC (Figure 6).

Figure 6

Download asset Open asset

Summary of findings.

The four nodes represent the investigated regions, namely the anterior insula (aINS), the dorsolateral and ventromedial parts of the prefrontal cortex (dlPFC and vmPFC, and the lateral orbitofrontal cortex lOFC). The outer disc represents the local mixed encoding i.e., the different proportions of contacts over time having a significant relationship between the gamma power and PE signals. In blue, is the proportion of contacts with a significant relation with the PE across rewarding trials (RPE-specific). Conversely, in red for punishment trials (PPE-specific). In purple, the proportion of contacts with a significant relationship with both the reward prediction error (RPE) and punishment prediction error (PPE). In gray, is the remaining proportion of non-significant contacts. Regarding interactions, we found that information transfer between aINS and dlPFC carried redundant information about PPE only and information transfer between vmPFC and lOFC about RPE only. This information transfer occurred with a leading role of the aINS in the punishment context and the vmPFC in the rewarding context. Finally, we found synergistic interactions between the dlPFC and the vmPFC about the full PE, without splitting into rewarding and punishing conditions.

Local mixed representations of prediction errors

Amongst the four investigated core-learning regions, the vmPFC was the only region to show a higher group-level preference for RPEs. This supports the notion that the vmPFC is functionally more specialized for the processing outcomes in reward learning, as previously put forward by human fMRI meta-analyses (Yacubian et al., 2006; Diekhof et al., 2012; Bartra et al., 2013; Garrison et al., 2013; Fouragnan et al., 2018). The dlPFC, instead, showed a stronger selectivity for punishment PE, thus supporting results from fMRI studies showing selective activations for aversive outcomes (Liu et al., 2011; Garrison et al., 2013; Fouragnan et al., 2018). On the contrary, the aINS and lOFC did not show clear selectivity for either reward or punishment PEs. The aINS carried a comparable amount of information about the RPE and PPE, thus suggesting that the insula is part of the surprise-encoding network (Fouragnan et al., 2018; Loued-Khenissi et al., 2020). Previous study reported a stronger link between the gamma activity of the aINS and the PPE compared to the RPE (Gueguen et al., 2021). This discrepancy in the results could be explained by the measures of information we are using here that are able to detect both linear and non-linear relationships between gamma activity and PE signals (Ince et al., 2017). The lOFC showed an initial temporal selectivity for PPE followed by a delayed one about the RPE. This is in accordance with fMRI and human intracranial studies which revealed that the lOFC was activated when receiving punishing outcomes, but also contains reward-related information (O’Doherty et al., 2001; Saez et al., 2018; Gueguen et al., 2021).

By taking advantage of the multi-site sampling of iEEG recordings, we quantified the heterogeneity in functional selectivity within each area and showed that the region-specific tendency toward either RPE or PPEs (Figure 2A) could be explained by the largest domain-specific proportion of contacts (Figure 2B). In other words, if a region showed a larger proportion of contacts being RPE-specific, the amount of information about the RPE at the group-level was also larger. Interestingly, we observed that 5 to 20% of contacts within a given region encoded both the RPE and PPE, thus revealing local mixed representations. Consequently, a strict dichotomous classification of learning-related areas as either reward, and punishment may fail to capture important properties of the individual nodes of the learning circuit, such as the functional heterogeneity in the encoding of PEs. These results suggest that the human prefrontal cortex exhibits a mixed local selectivity for prediction error signals at the mesoscopic scale. This view is in line with recent literature showing that the prefrontal cortex contains single neurons exhibiting mixed selectivity for multiple task variables (Meyers et al., 2008; Rigotti et al., 2013; Stokes et al., 2013; Panzeri et al., 2015; Parthasarathy et al., 2017; Bernardi et al., 2020). In the learning domain, single-unit studies have reported neurons encoding both rewarding and aversive outcomes in the OFC of the primate (Morrison and Salzman, 2009; Monosov and Hikosaka, 2012; Hirokawa et al., 2019). Mixed selectivity provides computational benefits, such as increasing the number of binary classifications, improving cognitive flexibility, and simplifying readout by downstream neurons (Fusi et al., 2016; Helfrich and Knight, 2019; Ohnuki et al., 2021; Panzeri et al., 2022). We suggest that the encoding of cognitive variables such as prediction error signals is supported by similar principles based on mixed selectivity at the meso- and macroscopic level, and may provide a natural substrate for cognitive flexibility and goal-directed learning (Rigotti et al., 2013).

Redundancy-dominated interactions segregate reward and punishment learning subsystems

We then tested whether the encoding of RPE and PPE signals could be supported by differential cortico-cortical interactions within and between frontal brain regions. To do so, we exploited the interaction information (II) (McGill, 1954; Ince et al., 2017) to quantify whether the amount of information bound up in a pair of gamma responses and PE signals is dominated by redundant or synergistic interactions (Williams and Beer, 2010). The II revealed redundancy-dominated interactions specific for RPE and PPE in the vmPFC and the dlPFC, respectively (Figure 3A). The aINS was the only region for which the between-contacts II did not increase the functional selectivity, with large redundant interactions for both RPE and PPE signals. This suggests that within-area redundant interactions can potentially amplify the functional specificity, despite the presence of local mixed selectivity (Figure 2A). Such ‘winner-take-all’ competition could be implemented by mutual inhibition mechanisms, which have been suggested to be essential in reward-guided choice (Hunt et al., 2012; Jocham et al., 2012; Strait et al., 2014; Hunt and Hayden, 2017).

Across-areas interaction information revealed two subsystems with redundancy-dominated interactions. A reward subsystem with RPE-specific interactions between the lOFC and vmPFC, and a punishment subsystem with PPE-specific interactions between the aINS and dlPFC (Figure 3B). Although a significant modulation selective for RPE was also present in the interaction between dlPFC and lOFC peaking around 500 ms after outcome presentation, a significant difference between the encoding of RPE and PPE was exclusively observed in the lOFC-vmPFC and aINS-dlPFC interactions (green clusters in Figure 3B). This result suggests that the observed functionally-distinct learning circuits for RPE and PPEs are associated with differential cortico-cortical interactions, rather than distinct local properties. More generally, our results suggest that redundancy-based network-level interactions are related to the functional specificity observed in neuroimaging and lesion studies (Pessiglione and Delgado, 2015; Palminteri and Pessiglione, 2017).

We then investigated differential communication patterns and directional relations within the two redundancy-dominated circuits (Kirst et al., 2016; Palmigiano et al., 2017). We identified significant information routing patterns, and dissociating reward and punishment learning (Figure 4). Within the reward subsystem, the vmPFC played a driving role toward the lOFC only during the rewarding condition. Conversely, within the punishment subsystem, the aINS played a driving role toward the dlPFC only during the punishment condition. These results support the notion that redundancy-dominated cognitive networks are associated with the occurrence of information-routing capabilities, where signals are communicated on top of collective reference states (Battaglia and Brovelli, 2020).

Here, we quantified directional relationships between regions using the transfer entropy (Schreiber, 2000), which is a functional connectivity measure based on the Granger-Wiener causality principle. Tract tracing studies in the macaque have revealed strong interconnections between the lOFC and vmPFC in the macaque (Carmichael and Price, 1996; Ongür and Price, 2000). In humans, cortico-cortical anatomical connections have mainly been investigated using diffusion magnetic resonance imaging (dMRI). Several studies found strong probabilities of structural connectivity between the anterior insula with the orbitofrontal cortex and the dorsolateral part of the prefrontal cortex (Cloutman et al., 2012; Ghaziri et al., 2017), and between the lOFC and vmPFC (Heather Hsu et al., 2020). In addition, the statistical dependency (e.g. coherence) between the LFP of distant areas could be potentially explained by direct anatomical connections (Schneider et al., 2021; Vinck et al., 2023). Taken together, the existence of an information transfer might rely on both direct or indirect structural connectivity. However, here we also reported differences in TE between rewarding and punishing trials given the same backbone anatomical connectivity (Figure 4). Our results are further supported by a recent study involving drug-resistant epileptic patients with resected insula who showed poorer performance than healthy controls in case of risky loss compared to risky gains (Von Siebenthal et al., 2017).

Encoding the full PE is supported by synergistic interactions between subsystems

Humans can flexibly switch between learning strategies that allow the acquisition of stimulus-action-outcomes associations in changing contexts. We investigated how RPE and PPE subsystems coordinated to allow such behavioral flexibility. To do so, we searched for neural correlates of PEs irrespectively of the context (reward or punishment learning) in between-regions interactions. We found that the encoding of global PE signals was associated with synergy-dominated interactions between the two subsystems, mediated by the interactions between the dlPFC and the vmPFC (Figure 5). Importantly, such synergy-dominated interaction reveals that the joint representation of the dlPFC and vmPFC is greater than the sum of their individual contributions to the encoding of global PE signals. Thus, it suggests that successful adaptation in varying contexts requires both the vmPFC and dlPFC for the encoding of global PE signals.

Role of redundant and synergistic interactions in brain network coordination

At the macroscopic level, few studies investigated the potential role of redundant and synergistic interactions. By combining functional and diffusion MRI, recent work suggested that redundant interactions are predominantly associated with structurally coupled and functionally segregated processing. In contrast, synergistic interactions preferentially support functional integrative processes and complex cognition across higher-order brain networks (Luppi et al., 2022). Triadic synergistic interactions between the continuous spike counts recorded within and across areas of the visuomotor network have been shown to carry behaviorally-relevant information and to display the strongest modulations during the processing of visual information and movement execution (Varley et al., 2023). Finally, cortical representations of prediction error signals in the acoustic domain observed tone-related and instantaneous redundant interactions, such as time-lagged synergistic interactions within and across temporal and frontal regions of the auditory system (Gelens et al., 2023).

At the microscopic level, the amount of information encoded by a population of neurons can be modulated by pairwise and higher-order interactions, producing varying fractions of redundancy and synergy (Averbeck et al., 2006; Panzeri et al., 2015; Panzeri et al., 2022). Synergistic and redundant pairs of neurons can be identified by estimating the amount of information contained in the joint representation minus the sum of the information carried by individual neurons (Schneidman et al., 2003). Redundant coding is intricately linked to correlated activity (Gutnisky and Dragoi, 2008) and can spontaneously emerge due to the spatial correlations present in natural scenes by triggering neurons with overlapping receptive fields. Correlations between the trial-by-trial variations of neuronal responses could limit the amount of information encoded by a population (Bartolo et al., 2020; Kafashan et al., 2021) and facilitate readout by downstream neurons (Salinas and Sejnowski, 2001). While redundancy has been at the heart of heated debates and influential theories, such as efficient coding and redundancy compression in sensory areas (Barlow, 2001), synergy phenomena have been described to a lesser extent. Recently, a study reported synergistic coding in a V1 cortical column together with structured correlations between synergistic and redundant hubs (Nigam et al., 2019). Taken together, we suggest that population codes with balancing proportions of redundancy and synergy offer a good compromise between system robustness and resilience to cell loss and the creation of new information (Panzeri et al., 2022). We suggest that redundancy-dominated interactions confer robustness and network-level selectivity for complementary learning processes, which may lead to functional integration processes. On the other hand, synergy-dominated interactions seem to support neural interactions between redundancy-dominated networks, thus supporting functional integrative processes in the brain. In addition, our study suggests that redundant and synergistic interactions occur across multiple spatial scales from local to large-scale.

Conclusion

Our report of mixed representation of reward and punishment prediction error signals explains the discrepancy in the attribution of a functional specificity to the core learning cortical regions. Instead, we propose that functional specialization for reward and punishment PE signals occurs with redundancy-dominated interactions within the two subsystems formed by the vmPFC-lOFC and aINS-dlPFC, respectively. Within each subsystem, we observed asymmetric and directional interactions with the vmPFC and aINS playing a driving role in the reward and punishment learning circuits. Finally, switching between reward and punishment learning was supported by synergistic collaboration between subsystems. This supports the idea that higher-order integration between functionally-distinct subsystems are mediated by synergistic interactions. Taken together, our results provide a unifying view reconciling distributed cortical representations with interactions supporting reward and punishment learning. They highlight the relevance of considering learning as a network-level phenomenon by linking distributed and functionally redundant subnetworks through synergistic interactions hence supporting flexible cognition (Fedorenko and Thompson-Schill, 2014; Petersen and Sporns, 2015; Bassett and Mattar, 2017; Hunt and Hayden, 2017; Averbeck and Murray, 2020; Averbeck and O’Doherty, 2022).

Methods

Data acquisition and experimental procedure

Intracranial EEG recordings

Intracranial electroencephalography (iEEG) recordings were collected from sixteen patients presenting pharmaco-resistant focal epilepsy and undergoing presurgical evaluation (33.5±12.4 years old, 10 females). As the location of the epileptic foci could not be identified through noninvasive methods, neural activity was monitored using intracranial stereotactic electroencephalography. Multi-lead and semi-rigid depth electrodes were stereotactically implanted according to the suspected origin of seizures. The selection of implantation sites was based solely on clinical aspects. iEEG recordings were performed at the clinical neurophysiology epilepsy departments of Grenoble and Lyon Hospitals (France). iEEG electrodes had a diameter of 0.8 mm, 2 mm wide, 1.5 mm apart, and contained 8–18 contact leads (Dixi, Besançon, France). For each patient, 5–17 electrodes were implanted. Recordings were conducted using an audio–video-EEG monitoring system (Micromed, Treviso, Italy), which allowed simultaneous recording of depth iEEG channels sampled at 512 Hz (six patients), or 1024 Hz (12 patients) [0.1–200 Hz bandwidth]. One of the contacts located in the white matter was used as a reference. Anatomical localizations of iEEG contacts were determined based on post-implant computed tomography scans or post-implant MRI scans coregistered with pre-implantation scans (Lachaux et al., 2003; Chouairi et al., 2022). All patients gave written informed consent and the study received approval from the ethics committee (CPP 09-CHUG-12, study 0907) and from a competent authority (ANSM no: 2009-A00239-48).

Limitations

iEEG have been collected from pharmacoresistant epileptic patients who underwent deep electrode probing for preoperative evaluation. However, we interpreted these data as if collected from healthy subjects and assumed that epileptic activity does not affect the neural realization of prediction error. To best address this question, we excluded electrodes contaminated with pathological activity and focused on task-related changes and multi-trial analysis to reduce the impact of incorrect or task-independent neural activations. Therefore, our results may benefit from future replication in healthy controls using non-invasive recordings. Despite the aforementioned limitations, we believe that access to deep intracerebral EEG recordings of human subjects can provide privileged insight into the neural dynamics that regulate human cognition, with outstanding spatial, temporal, and spectral precision. In the long run, this type of data could help bridge the gap between neuroimaging studies and electrophysiological recordings in nonhuman primates.

Preprocessing of iEEG data

Bipolar derivations were computed between adjacent electrode contacts to diminish contributions of distant electric sources through volume conduction, reduce artifacts, and increase the spatial specificity of the neural data. Bipolar iEEG signals can approximately be considered as originating from a cortical volume centered within two contacts (Brovelli et al., 2005; Bastin et al., 2016; Combrisson et al., 2017), thus providing a spatial resolution of approximately 1.5–3 mm (Lachaux et al., 2003; Jerbi et al., 2009; Chouairi et al., 2022). Recording sites with artifacts and pathological activity (e.g. epileptic spikes) were removed using visual inspection of all of the traces of each site and each participant.

Definition of anatomical regions of interest

Anatomical labeling of bipolar derivations was performed using the IntrAnat software (Deman et al., 2018). The 3D T1 pre-implantation MRI gray/white matter was segmented and spatially normalized to obtain a series of cortical parcels using MarsAtlas (Auzias et al., 2016) and the Destrieux atlas (Destrieux et al., 2010). 3D coordinates of electrode contacts were then coregistered on post-implantation images (MRI or CT). Each recording site (i.e. bipolar derivation) was labeled according to its position in a parcellation scheme in the participant’s native space. Thus, the analyzed dataset only included electrodes identified to be in the gray matter. Four regions of interest (ROIs) were defined for further analysis: (1) the ventromedial prefrontal cortex (vmPFC) ROI was created by merging six (three per hemisphere) parcels in MarsAlas (labeled PFCvm, OFCv, and OFCvm in MarsAtlas) corresponding to the ventromedial prefrontal cortex and fronto-medial part of the orbitofrontal cortex, respectively; (2) the lateral orbitofrontal cortex (lOFC) ROI included four (two per hemisphere) MarsAtlas parcels (MarsAtlas labels: OFCvl and the OFCv); (3) the dorsolateral prefrontal cortex (dlPFC) ROI was defined as the inferior and superior bilateral dorsal prefrontal cortex (MarsAtlas labels: PFrdli and PFrdls); (4) the anterior insula (aINS) ROI was defined as the bilateral anterior part of the insula (Destrieux atlas labels: Short insular gyri, anterior circular insular sulcus and anterior portion of the superior circular insular sulcus). The total number of bipolar iEEG derivations for the four ROIS was 44, 59, 70, and 75 for the vmPFC, lOFC, dlPFC, and aINS, respectively (Figure 1A). As channels with artifacts or epileptic activities were removed here, the number of recordings differs from a previous study (Gueguen et al., 2021).

Behavioral task and set-up

Participants were asked to participate in a probabilistic instrumental learning task adapted from previous studies (Pessiglione et al., 2006; Palminteri et al., 2012). Participants received written instructions that the goal of the task was to maximize their financial payoff by considering reward-seeking and punishment avoidance as equally important. Instructions were reformulated orally if necessary. Participants started with a short session, with only two pairs of cues presented on 16 trials, followed by 2–3 short sessions of 5 min. At the end of this short training, all participants were familiar with the timing of events, with the response buttons and all reached a threshold of at least 70% of correct choices during both reward and punishment conditions. Participants then performed three to six sessions on a single testing occurrence, with short breaks between sessions. Each session was an independent task, with four new pairs of cues to be learned. Cues were abstract visual stimuli taken from the Agathodaimon alphabet. The two cues of a pair were always presented together on the left and right of a central fixation cross and their relative position was counterbalanced across trials. On each trial, one pair was randomly presented. Each pair of cues was presented 24 times for a total of 96 trials per session. The four pairs of cues were divided into two conditions. A rewarding condition where the two pairs could either lead the participants to win one euro or nothing (+1€ vs. 0€) and a symmetric punishment condition where the participants could either lose one euro or nothing (–1€ vs. 0€). Rewarding and punishing pairs of cues were presented in an intermingled random manner and participants had to learn the four pairs at once. Within each pair, the two cues were associated with the two possible outcomes with reciprocal probabilities (0.75/0.25 and 0.25/0.75). To choose between the left or right cues, participants used their left or right index to press the corresponding button on a joystick (Logitech Dual Action). Since the position on the screen was counterbalanced, response (left versus right) and value (good vs. bad cue) were orthogonal. The chosen cue was colored in red for 250 ms and then the outcome was displayed on the screen after 1000 ms. To win money, participants had to learn by trial and error which cue-outcome association was the most rewarding in the rewarding condition and the least penalizing in the punishment condition. Visual stimuli were delivered on a 19-inch TFT monitor with a refresh rate of 60 Hz, controlled by a PC with Presentation 16.5 (Neurobehavioral Systems, Albany, CA).

Computational model of learning

To model choice behavior and estimate prediction error signals, we used a standard Q-learning model (Watkins and Dayan, 1992) from reinforcement learning theory (Sutton and Barto, 2018). For a pair of cues A and B, the model estimates the expected value of choosing A (Qa) or B (Qb), given previous choices and received outcomes. Q-values were initiated to 0, corresponding to the average of all possible outcome values. After each trial t, the expected value of choosing a stimulus (e.g. A) was updated according to the following update rule:

Q_{a_{t + 1}} = Q_{a_{t}} + α δ_{t}

with α the learning rate weighting the importance given to new experiences and δ, the outcome prediction error signals at a trial t defined as the difference between the obtained and expected outcomes:

δ_{t} = R_{t} - Q a_{t}

with R_t the reinforcement value among –1€, 0€, and 1€. The probability of choosing a cue was then estimated by transforming the expected values associated with each cue using a softmax rule with a Gibbs distribution. An additional Ө parameter was added in the softmax function to the expected value of the chosen option on the previous trial of the same cue to account for the tendency to repeat the choice made on the previous trial. For example, if a participant chose option A on trial t, the probability of choosing A at trial t+1 was obtained using:

P a_{t + 1} = \frac{e^{Q a_{t} + θ / β}}{e^{Q a_{t} + θ / β} + e^{Q b_{t} / β}}

with β the choice temperature for controlling the ratio between exploration and exploitation. The three free parameters α, β, and Ө were fitted per participant and optimized by minimizing the negative log-likelihood of choice using the MATLAB fmincon function, initialized at multiple starting points of the parameter space (Palminteri et al., 2015). Estimates of the free parameters, the goodness of fit and the comparison between modeled and observed data can be seen in Table 1 and Figure 1 in Gueguen et al., 2021.

Table 1

Results of the one-sample t-test performed against 0.

	T-value	p-value	p-value(FDR corrected)	dof	CI 95%
II_PPE-RPE	2859	0.007**	0.015*	34	[6.5e-05, 3.9e-04]
II_Mixed	2841	0.008**	0.015*	33	[5.4e-05, 3.3e-04]
II_PPE-PPE	1,25	0.2667	0.3556	5	[–7.1e-05, 2.1e-04]
II_RPE-RPE	0733	0.4912	0.4912	6	[–3.1e-05, 5.8e-05]

iEEG data analysis

Estimate of single-trial gamma-band activity

Here, we focused solely on broadband gamma for three main reasons. First, it has been shown that the gamma band activity correlates with both spiking activity and the BOLD fMRI signals (Mukamel et al., 2005; Niessing et al., 2005; Lachaux et al., 2007; Nir et al., 2007), and it is commonly used in MEG and iEEG studies to map task-related brain regions (Brovelli et al., 2005; Crone et al., 2006; Vidal et al., 2006; Ball et al., 2008; Jerbi et al., 2009; Lachaux et al., 2012; Cheyne and Ferrari, 2013). Therefore, focusing on the gamma band facilitates linking our results with the fMRI and spiking literature on probabilistic learning. Second, single-trial and time-resolved high-gamma activity can be exploited for the analysis of cortico-cortical interactions in humans using MEG and iEEG techniques (Brovelli et al., 2015; Brovelli et al., 2017; Combrisson et al., 2022a). Finally, while previous analyses of the current dataset (Gueguen et al., 2021) reported an encoding of PE signals at different frequency bands, the power in lower frequency bands were shown to carry redundant information compared to the gamma band power. In the current study, we thus estimated the power in the gamma band using a multitaper time-frequency transform based on Slepian tapers (Percival and Walden, 1993; Mitra and Pesaran, 1999). To extract gamma-band activity from 50 to 100 Hz, the iEEG time series were multiplied by 9 orthogonal tapers (15 cycles for a duration of 200ms and with a time-bandwidth for frequency smoothing of 10 Hz), centered at 75 Hz and Fourier-transformed. To limit false negative proportions due to multiple testings, we down-sampled the gamma power to 256 Hz. Finally, we smoothed the gamma power using a 10-point Savitzky-Golay filter. We used MNE-Python (Gramfort et al., 2013) to inspect the time series, reject contacts contaminated with pathological activity, and estimate the power spectrum density (mne.time_frequency.psd_multitaper) and the gamma power (mne.time_frequency.tfr_multitaper).

Local correlates of PE signals

To quantify the local encoding of prediction error (PE) signals in the four ROIs, we used information-theoretic metrics. To this end, we computed the time-resolved mutual information (MI) between the single-trial gamma-band responses and the outcome-related PE signals. As a reminder, mutual information is defined as:

I (X, Y) = H (X) - H (X | Y)

In this equation, the variables X and Y represent the across-trials gamma-band power and the PE variables, respectively. H(X) is the entropy of X, and H(X|Y) is the conditional entropy of X given Y. In the current study, we used a semi-parametric binning-free technique to calculate MI, called Gaussian-Copula Mutual Information (GCMI) (Ince et al., 2017). The GCMI is a robust rank-based approach that allows the detection of any type of monotonic relation between variables and it has been successfully applied to brain signals analysis (Colenbier et al., 2020; Michelmann et al., 2021; Ten Oever et al., 2021). Mathematically, the GCMI is a lower-bound estimation of the true MI and it does not depend on the marginal distributions of the variables, but only on the copula function that encapsulates their dependence. The rank-based copula-normalization preserves the relationship between variables as long as this relation is strictly increasing or decreasing. As a consequence, the GCMI can only detect monotonic relationships. Nevertheless, the GCMI is of practical importance for brain signal analysis for several reasons. It allows to estimate the MI on a limited number of samples and it contains a parametric bias correction to compensate for the bias due to the estimation on smaller datasets. It allows to compute the MI on uni- and multivariate variables that can either be continuous or discrete see Table 1 in Ince et al., 2017. Finally, it is computationally efficient, which is a desired property when dealing with a large number of iEEG contacts recording at a high sampling rate. Here, the GCMI was computed across trials and it was used to estimate the instantaneous amount of information shared between the gamma power of iEEG contacts and RPE (MI_RPE = I(γ; RPE)) and PPE signals (M_PPE = I(γ; PPE)).

Network-level interactions and PE signals

The goal of network-level analyses was to characterize the nature of cortico-cortical interactions encoding reward and punishment PE signals. In particular, we aimed to quantify: (1) the nature of the interdependence between pairs of brain ROIs in the encoding of PE signals; (2) the information flow between ROIs encoding PE signals. These two questions were addressed using Interaction Information and Transfer Entropy analyses, respectively.

Interaction Information analysis

In classical information theory, interaction information (II) provides a generalization of mutual information for more than two variables (McGill, 1954; Ince et al., 2017). For the three-variables case, the II can be defined as the difference between the total, or joint, mutual information between ROIs (R₁ and R₂) and the third behavioral variable (S), minus the two individual mutual information between each ROI and the behavioral variable. For a three variables multivariate system composed of two sources R₁, R₂, and a target S, the II is defined as:

\begin{array}{ll} I I (R_{1}; R_{2}; S) & = I (S; R_{1} ∣ R_{2}) - I (R_{1}; S) \\ = I (R_{1}, R_{2}; S) - I (R_{1}; S) - (R_{2}; S) \end{array}

Unlike mutual information, the interaction information can be either positive or negative. A negative value of interaction information indicates a net redundant effect between variables, whereas positive values indicate a net synergistic effect (Williams and Beer, 2010). Here, we used the II to investigate the amount of information and the nature of the interactions between the gamma power of pairs of contacts (γ₁, γ₂) about the RPE (II_RPE = II(γ₁, γ₂; RPE)) and PPE signals (II_PPE = II(γ₁, γ₂; PPE)). The II was computed by estimating the MI quantities of equation (5) using the GCMI between contacts within the same brain region or across different regions.

Transfer entropy analysis

To quantify the degree of communication between neural signals, the most successful model-free methods rely on the Wiener-Granger principle (Wiener, 1956; Granger, 1969). This principle identifies information flow between time series when future values of a given signal can be predicted from the past values of another, above and beyond what can be achieved from its autocorrelation. One of the most general information theoretic measures based on the Wiener-Granger principle is Transfer Entropy (TE) (Schreiber, 2000). The TE can be formulated in terms of conditional mutual information (Schreiber, 2000; Kaiser and Schreiber, 2002):

T E (X \to Y) = I (X_{P a s t}; Y_{t} ∣ Y_{p a s t})

Here, we computed the TE on the gamma activity time courses of pairs of iEEG contacts. We used the GCMI to estimate conditional mutual information. For an interval [d₁, d₂] of n_delays, the final TE estimation was defined as the mean over the TE estimated at each delay:

T E (X \to Y)_{[d_{1}, d_{2}]} = \frac{1}{n_{d e l a y s}} . \sum_{d = d_{1}}^{d_{2}} I (X_{d}; Y_{t} ∣ Y_{d})

Statistical analysis

We used a group-level approach based on non-parametric permutations, encompassing non-negative measures of information (Combrisson et al., 2022a). The same framework was used at the local level (i.e. the information carried by a single contact) or at the network level (i.e. the information carried by pairs of contacts for the II and TE). To take into account the inherent variability existing at the local and network levels, we used a random-effect model. To generate the distribution of permutations at the local level, we shuffled the PE variable across trials 1000 times and computed the MI between the gamma power and the shuffled version of the PE. The shuffling led to a distribution of MI reachable by chance, for each contact and at each time point (Combrisson and Jerbi, 2015). To form the group-level effect, we computed a one-sample t-test against the permutation mean across the MI computed on individual contacts taken from the same brain region, at each time point. The same procedure was used on the permutation distribution to form the group-level effect reachable by chance. We used cluster-based statistics to correct for multiple comparisons (Maris and Oostenveld, 2007). The cluster-forming threshold was defined as the 95^th percentile of the distribution of t-values obtained from the permutations. We used this threshold to form the temporal clusters within each brain region. We obtained cluster masses on both the true t-values and the t-values computed on the permutations. To correct for multiple comparisons, we built a distribution made of the largest 1000 cluster masses estimated on the permuted data. The final corrected p-values were inferred as the proportion of permutations exceeding the t-values. To generate the distributions of II and TE reachable by chance, we respectively shuffled the PE variable across trials for the II and the gamma power across trials of the source for the TE (Vicente et al., 2011). The rest of the significance testing procedure at the network level is similar to the local level, except that it is not applied within brain regions but within pairs of brain regions.

Software

Information-theoretic metrics and group-level statistics, are implemented in a homemade Python software called Frites (Combrisson et al., 2022b). The interaction information can be computed using the frites.conn.conn_ii function and the transfer entropy using the frites.conn.conn_te function.

Data availability

The Python scripts and notebooks to reproduce the results presented here are hosted on GitHub, copy archived at Combrisson, 2024. The preprocessed data used here can be downloaded from Dryad.

The following data sets were generated

1. Combrisson E
2. Basanisi R
3. Gueguen MCM
4. Rheims S
5. Kahane P
6. Bastin J
7. Brovelli A
(2024) Dryad Digital Repository
Neural interactions in the human frontal cortex dissociate reward and punishment learning.
https://doi.org/10.5061/dryad.jdfn2z3k4

References

(2016) MarsAtlas: a cortical parcellation atlas for functional mapping
Human Brain Mapping 37:1573–1592.
https://doi.org/10.1002/hbm.23121
- PubMed
- Google Scholar
(2006) Neural correlations, population coding and computation
Nature Reviews. Neuroscience 7:358–366.
https://doi.org/10.1038/nrn1888
- PubMed
- Google Scholar
1. Averbeck BB
2. Murray EA
(2020) Hypothalamic interactions with large-scale neural circuits underlying reinforcement learning and motivated behavior
Trends in Neurosciences 43:681–694.
https://doi.org/10.1016/j.tins.2020.06.006
- PubMed
- Google Scholar
1. Averbeck B
2. O’Doherty JP
(2022) Reinforcement-learning in fronto-striatal circuits
Neuropsychopharmacology 47:147–162.
https://doi.org/10.1038/s41386-021-01108-0
- PubMed
- Google Scholar
1. Ball T
2. Demandt E
3. Mutschler I
4. Neitzel E
5. Mehring C
6. Vogt K
7. Aertsen A
8. Schulze-Bonhage A
(2008) Movement related activity in the high gamma range of the human EEG
NeuroImage 41:302–310.
https://doi.org/10.1016/j.neuroimage.2008.02.032
- PubMed
- Google Scholar
1. Balleine BW
2. Dickinson A
(1998) Goal-directed instrumental action: contingency and incentive learning and their cortical substrates
Neuropharmacology 37:407–419.
https://doi.org/10.1016/s0028-3908(98)00033-1
- PubMed
- Google Scholar
1. Balleine BW
2. O’Doherty JP
(2010) Human and rodent homologies in action control: corticostriatal determinants of goal-directed and habitual action
Neuropsychopharmacology 35:48–69.
https://doi.org/10.1038/npp.2009.131
- PubMed
- Google Scholar
1. Balleine BW
(2019) The meaning of behavior: discriminating reflex and volition in the brain
Neuron 104:47–62.
https://doi.org/10.1016/j.neuron.2019.09.024
- PubMed
- Google Scholar
1. Barlow H
(2001) Redundancy reduction revisited
Network 12:241–253.
https://doi.org/10.1080/net.12.3.241.253
- PubMed
- Google Scholar
(2020) Information-limiting correlations in large neural populations
The Journal of Neuroscience 40:1668–1678.
https://doi.org/10.1523/JNEUROSCI.2072-19.2019
- PubMed
- Google Scholar
(2013) The valuation system: a coordinate-based meta-analysis of BOLD fMRI experiments examining neural correlates of subjective value
NeuroImage 76:412–427.
https://doi.org/10.1016/j.neuroimage.2013.02.063
- PubMed
- Google Scholar
1. Bassett DS
2. Mattar MG
(2017) A Network neuroscience of human learning: potential to inform quantitative theories of brain and behavior
Trends in Cognitive Sciences 21:250–264.
https://doi.org/10.1016/j.tics.2017.01.010
- PubMed
- Google Scholar
1. Bastin J
2. Deman P
3. David O
4. Gueguen M
5. Benis D
6. Minotti L
7. Hoffman D
8. Combrisson E
9. Kujala J
10. Perrone-Bertolotti M
11. Kahane P
12. Lachaux JP
13. Jerbi K
(2016) Direct recordings from human anterior insula reveal its leading role within the error-monitoring network
Cerebral Cortex 01:bhv352.
https://doi.org/10.1093/cercor/bhv352
- Google Scholar
Book
1. Battaglia D
2. Brovelli A
(2020) Functional Connectivity and Neuronal Dynamics: Insights from Computational Methods
MIT Press.
https://doi.org/10.7551/mitpress/11442.001.0001
- Google Scholar
1. Bernardi S
2. Benna MK
3. Rigotti M
4. Munuera J
5. Fusi S
6. Salzman CD
(2020) The geometry of abstraction in the hippocampus and prefrontal cortex
Cell 183:954–967.
https://doi.org/10.1016/j.cell.2020.09.031
- PubMed
- Google Scholar
1. Bódi N
2. Kéri S
3. Nagy H
4. Moustafa A
5. Myers CE
6. Daw N
7. Dibó G
8. Takáts A
9. Bereczki D
10. Gluck MA
(2009) Reward-learning and the novelty-seeking personality: a between- and within-subjects study of the effects of dopamine agonists on young parkinson’s patients
Brain 132:2385–2395.
https://doi.org/10.1093/brain/awp094
- PubMed
- Google Scholar
Book
1. Bouton ME
(2007)
Learning and Behavior: A Contemporary Synthesis

OUP USA.
- Google Scholar
1. Braun U
2. Schäfer A
3. Walter H
4. Erk S
5. Romanczuk-Seiferth N
6. Haddad L
7. Schweiger JI
8. Grimm O
9. Heinz A
10. Tost H
11. Meyer-Lindenberg A
12. Bassett DS
(2015) Dynamic reconfiguration of frontal brain networks during executive cognition in humans
PNAS 112:11678–11683.
https://doi.org/10.1073/pnas.1422487112
- PubMed
- Google Scholar
1. Bressler SL
2. Menon V
(2010) Large-scale brain networks in cognition: emerging methods and principles
Trends in Cognitive Sciences 14:277–290.
https://doi.org/10.1016/j.tics.2010.04.004
- PubMed
- Google Scholar
(2005) High gamma frequency oscillatory activity dissociates attention from intention in the human premotor cortex
NeuroImage 28:154–164.
https://doi.org/10.1016/j.neuroimage.2005.05.045
- PubMed
- Google Scholar
1. Brovelli XA
2. Chicharro D
3. Badier JM
4. Wang H
5. Jirsa V
(2015) Characterization of cortical networks and corticocortical functional connectivity mediating arbitrary visuomotor mapping
The Journal of Neuroscience 35:12643–12658.
https://doi.org/10.1523/JNEUROSCI.4892-14.2015
- PubMed
- Google Scholar
1. Brovelli A
2. Badier JM
3. Bonini F
4. Bartolomei F
5. Coulon O
6. Auzias G
(2017) Dynamic reconfiguration of visuomotor-related functional connectivity networks
The Journal of Neuroscience 37:839–853.
https://doi.org/10.1523/JNEUROSCI.1672-16.2016
- PubMed
- Google Scholar
1. Buehlmann A
2. Deco G
(2010) Optimal information transfer in the cortex through synchronization
PLOS Computational Biology 6:e1000934.
https://doi.org/10.1371/journal.pcbi.1000934
- PubMed
- Google Scholar
1. Buzsáki G
2. Draguhn A
(2004) Neuronal oscillations in cortical networks
Science 304:1926–1929.
https://doi.org/10.1126/science.1099745
- PubMed
- Google Scholar
1. Carmichael ST
2. Price JL
(1996) Connectional networks within the orbital and medial prefrontal cortex of macaque monkeys
The Journal of Comparative Neurology 371:179–207.
https://doi.org/10.1002/(SICI)1096-9861(19960722)371:2<179::AID-CNE1>3.0.CO;2-#
- Google Scholar
1. Cheyne D
2. Ferrari P
(2013) MEG studies of motor cortex gamma oscillations: evidence for a gamma “fingerprint” in the brain?
Frontiers in Human Neuroscience 7:575.
https://doi.org/10.3389/fnhum.2013.00575
- PubMed
- Google Scholar
(2022) Preoperative deficiency anemia in digital replantation: a marker of disparities, increased length of stay, and hospital cost
Journal of Hand and Microsurgery 14:147–152.
https://doi.org/10.1055/s-0040-1714152
- PubMed
- Google Scholar
(2012) The variation of function across the human insula mirrors its patterns of structural connectivity: evidence from in vivo probabilistic tractography
NeuroImage 59:3514–3521.
https://doi.org/10.1016/j.neuroimage.2011.11.016
- PubMed
- Google Scholar
1. Cohen JR
2. D’Esposito M
(2016) The Segregation and Integration of Distinct Brain Networks and Their Relationship to Cognition
The Journal of Neuroscience 36:12083–12094.
https://doi.org/10.1523/JNEUROSCI.2965-15.2016
- PubMed
- Google Scholar
(2020) Disambiguating the role of blood flow and global signal with partial information decomposition
NeuroImage 213:116699.
https://doi.org/10.1016/j.neuroimage.2020.116699
- PubMed
- Google Scholar
1. Combrisson E
2. Jerbi K
(2015) Exceeding chance level by chance: the caveat of theoretical chance levels in brain signal classification and statistical assessment of decoding accuracy
Journal of Neuroscience Methods 250:126–136.
https://doi.org/10.1016/j.jneumeth.2015.01.010
- PubMed
- Google Scholar
(2017) From intentions to actions: neural oscillations encode motor processes through phase, amplitude and phase-amplitude coupling
NeuroImage 147:473–487.
https://doi.org/10.1016/j.neuroimage.2016.11.042
- PubMed
- Google Scholar
1. Combrisson E
2. Allegra M
3. Basanisi R
4. Ince RAA
5. Giordano BL
6. Bastin J
7. Brovelli A
(2022a) Group-level inference of information-based measures for the analyses of cognitive brain networks from neurophysiological data
NeuroImage 258:119347.
https://doi.org/10.1016/j.neuroimage.2022.119347
- PubMed
- Google Scholar
(2022b) Frites: a python package for functional connectivityanalysis and group-level statistics of neurophysiological data
Journal of Open Source Software 7:3842.
https://doi.org/10.21105/joss.03842
- Google Scholar
Software
1. Combrisson E
(2024) Papercode, version swh:1:rev:7772b6216b89bd783eb6895fc9199d1e1f97462c
Software Heritage.

https://archive.softwareheritage.org/swh:1:dir:d0f6a1bc4776dce6390c104511e78c8e30f51a89;origin=https://github.com/brainets/papercode;visit=swh:1:snp:f0cbea21b2baf19ef23042c111ebd0df79deab3e;anchor=swh:1:rev:7772b6216b89bd783eb6895fc9199d1e1f97462c
(2006) High-frequency gamma oscillations and human brain mapping with electrocorticography
Progress in Brain Research 159:275–295.
https://doi.org/10.1016/S0079-6123(06)59019-3
- PubMed
- Google Scholar
(2008) BOLD responses reflecting dopaminergic signals in the human ventral tegmental area
Science 319:1264–1267.
https://doi.org/10.1126/science.1150605
- PubMed
- Google Scholar
1. Deco G
2. Tononi G
3. Boly M
4. Kringelbach ML
(2015) Rethinking segregation and integration: contributions of whole-brain modelling
Nature Reviews. Neuroscience 16:430–439.
https://doi.org/10.1038/nrn3963
- PubMed
- Google Scholar
1. Deman P
2. Bhattacharjee M
3. Tadel F
4. Job AS
5. Rivière D
6. Cointepas Y
7. Kahane P
8. David O
(2018) Intranat electrodes: a free database and visualization software for intracranial electroencephalographic data processed for case and group studies
Frontiers in Neuroinformatics 12:40.
https://doi.org/10.3389/fninf.2018.00040
- PubMed
- Google Scholar
1. Destrieux C
2. Fischl B
3. Dale A
4. Halgren E
(2010) Automatic parcellation of human cortical gyri and sulci using standard anatomical nomenclature
NeuroImage 53:1–15.
https://doi.org/10.1016/j.neuroimage.2010.06.010
- PubMed
- Google Scholar
1. Dickinson A
2. Balleine B
(1994) Motivational control of goal-directed action
Animal Learning & Behavior 22:1–18.
https://doi.org/10.3758/BF03199951
- Google Scholar
1. Diekhof EK
2. Kaps L
3. Falkai P
4. Gruber O
(2012) The role of the human ventral striatum and the medial orbitofrontal cortex in the representation of reward magnitude - an activation likelihood estimation meta-analysis of neuroimaging studies of passive reward expectancy and outcome processing
Neuropsychologia 50:1252–1266.
https://doi.org/10.1016/j.neuropsychologia.2012.02.007
- PubMed
- Google Scholar
1. Dolan RJ
2. Dayan P
(2013) Goals and habits in the brain
Neuron 80:312–325.
https://doi.org/10.1016/j.neuron.2013.09.007
- PubMed
- Google Scholar
(2001) Dynamic predictions: oscillations and synchrony in top-down processing
Nature Reviews. Neuroscience 2:704–716.
https://doi.org/10.1038/35094565
- PubMed
- Google Scholar
1. Fedorenko E
2. Thompson-Schill SL
(2014) Reworking the language network
Trends in Cognitive Sciences 18:120–126.
https://doi.org/10.1016/j.tics.2013.12.006
- PubMed
- Google Scholar
1. Finc K
2. Bonna K
3. He X
4. Lydon-Staley DM
5. Kühn S
6. Duch W
7. Bassett DS
(2020) Dynamic reconfiguration of functional brain networks during working memory training
Nature Communications 11:2435.
https://doi.org/10.1038/s41467-020-15631-z
- PubMed
- Google Scholar
(2018) Separate neural representations of prediction error valence and surprise: evidence from an fMRI meta-analysis
Human Brain Mapping 39:2887–2906.
https://doi.org/10.1002/hbm.24047
- PubMed
- Google Scholar
(2004) By carrot or by stick: cognitive reinforcement learning in parkinsonism
Science 306:1940–1943.
https://doi.org/10.1126/science.1102941
- PubMed
- Google Scholar
1. Fries P
(2015) Rhythms for cognition: communication through coherence
Neuron 88:220–235.
https://doi.org/10.1016/j.neuron.2015.09.034
- PubMed
- Google Scholar
(2016) Why neurons mix: high dimensionality for higher cognition
Current Opinion in Neurobiology 37:66–74.
https://doi.org/10.1016/j.conb.2016.01.010
- PubMed
- Google Scholar
(2013) Prediction error in reinforcement learning: a meta-analysis of neuroimaging studies
Neuroscience and Biobehavioral Reviews 37:1297–1310.
https://doi.org/10.1016/j.neubiorev.2013.03.023
- PubMed
- Google Scholar
1. Gelens F
2. Äijälä J
3. Roberts L
4. Komatsu M
5. Uran C
6. Jensen MA
7. Miller KJ
8. Ince RAA
9. Garagnani M
10. Vinck M
11. Canales-Johnson A
(2023) Distributed representations of prediction error signals across the cortical hierarchy are synergistic
Neuroscience 01:e3735.
https://doi.org/10.1101/2023.01.12.523735
- Google Scholar
1. Ghaziri J
2. Tucholka A
3. Girard G
4. Houde JC
5. Boucher O
6. Gilbert G
7. Descoteaux M
8. Lippé S
9. Rainville P
10. Nguyen DK
(2017) The corticocortical structural connectivity of the human insula
Cerebral Cortex 27:1216–1228.
https://doi.org/10.1093/cercor/bhv308
- PubMed
- Google Scholar
1. Gramfort A
2. Luessi M
3. Larson E
4. Engemann DA
5. Strohmeier D
6. Brodbeck C
7. Goj R
8. Jas M
9. Brooks T
10. Parkkonen L
11. Hämäläinen M
(2013) MEG and EEG data analysis with MNE-Python
Frontiers in Neuroscience 7:267.
https://doi.org/10.3389/fnins.2013.00267
- PubMed
- Google Scholar
1. Granger CWJ
(1969) Investigating causal relations by econometric models and cross-spectral methods
Econometrica 37:424.
https://doi.org/10.2307/1912791
- Google Scholar
1. Gueguen MCM
2. Lopez-Persem A
3. Billeke P
4. Lachaux JP
5. Rheims S
6. Kahane P
7. Minotti L
8. David O
9. Pessiglione M
10. Bastin J
(2021) Anatomical dissociation of intracerebral signals for reward and punishment prediction errors in humans
Nature Communications 12:3344.
https://doi.org/10.1038/s41467-021-23704-w
- PubMed
- Google Scholar
1. Gutnisky DA
2. Dragoi V
(2008) Adaptive coding of visual information in neural populations
Nature 452:220–224.
https://doi.org/10.1038/nature06563
- PubMed
- Google Scholar
1. Heather Hsu CC
2. Rolls ET
3. Huang CC
4. Chong ST
5. Zac Lo CY
6. Feng J
7. Lin CP
(2020) Connections of the human orbitofrontal cortex and inferior frontal gyrus
Cerebral Cortex 30:5830–5843.
https://doi.org/10.1093/cercor/bhaa160
- PubMed
- Google Scholar
Book
1. Helfrich RF
2. Knight RT
(2019) Cognitive Neurophysiology of the Prefrontal cortex
In: D’Esposito M, Grafman JH, editors. Handbook of Clinical Neurology. Elsevier. pp. 35–59.
https://doi.org/10.1016/B978-0-12-804281-6.00003-3
- Google Scholar
1. Hirokawa J
2. Vaughan A
3. Masset P
4. Ott T
5. Kepecs A
(2019) Frontal cortex neuron types categorically encode single decision variables
Nature 576:446–451.
https://doi.org/10.1038/s41586-019-1816-9
- PubMed
- Google Scholar
(2012) Mechanisms underlying cortical activity during value-guided choice
Nature Neuroscience 15:470–476.
https://doi.org/10.1038/nn.3017
- PubMed
- Google Scholar
1. Hunt LT
2. Hayden BY
(2017) A distributed, hierarchical and recurrent framework for reward-based choice
Nature Reviews. Neuroscience 18:172–182.
https://doi.org/10.1038/nrn.2017.7
- PubMed
- Google Scholar
1. Ince RAA
2. Giordano BL
3. Kayser C
4. Rousselet GA
5. Gross J
6. Schyns PG
(2017) A statistical framework for neuroimaging data analysis based on mutual information estimated via A gaussian copula
Human Brain Mapping 38:1541–1573.
https://doi.org/10.1002/hbm.23471
- PubMed
- Google Scholar
1. Jerbi K
2. Ossandón T
3. Hamamé CM
4. Senova S
5. Dalal SS
6. Jung J
7. Minotti L
8. Bertrand O
9. Berthoz A
10. Kahane P
11. Lachaux JP
(2009) Task-related gamma-band dynamics from an intracerebral perspective: review and implications for surface EEG and MEG
Human Brain Mapping 30:1758–1771.
https://doi.org/10.1002/hbm.20750
- PubMed
- Google Scholar
1. Jocham G
2. Hunt LT
3. Near J
4. Behrens TEJ
(2012) A mechanism for value-guided choice based on the excitation-inhibition balance in prefrontal cortex
Nature Neuroscience 15:960–961.
https://doi.org/10.1038/nn.3140
- PubMed
- Google Scholar
(2021) Scaling of sensory information in large neural populations shows signatures of information-limiting correlations
Nature Communications 12:473.
https://doi.org/10.1038/s41467-020-20722-y
- PubMed
- Google Scholar
1. Kaiser A
2. Schreiber T
(2002) Information transfer in continuous processes
Physica D 166:43–62.
https://doi.org/10.1016/S0167-2789(02)00432-3
- Google Scholar
(2016) Dynamic information routing in complex networks
Nature Communications 7:11061.
https://doi.org/10.1038/ncomms11061
- PubMed
- Google Scholar
(2003) Intracranial EEG and human brain mapping
Journal of Physiology, Paris 97:613–628.
https://doi.org/10.1016/j.jphysparis.2004.01.018
- PubMed
- Google Scholar
1. Lachaux JP
2. Fonlupt P
3. Kahane P
4. Minotti L
5. Hoffmann D
6. Bertrand O
7. Baciu M
(2007) Relationship between task-related gamma oscillations and BOLD signal: new insights from combined fMRI and intracranial EEG
Human Brain Mapping 28:1368–1375.
https://doi.org/10.1002/hbm.20352
- PubMed
- Google Scholar
(2012) High-frequency neural activity and human cognition: past, present and possible future of intracranial EEG research
Progress in Neurobiology 98:279–301.
https://doi.org/10.1016/j.pneurobio.2012.06.008
- PubMed
- Google Scholar
1. Liu X
2. Hairston J
3. Schrier M
4. Fan J
(2011) Common and distinct networks underlying reward valence and processing stages: a meta-analysis of functional neuroimaging studies
Neuroscience and Biobehavioral Reviews 35:1219–1236.
https://doi.org/10.1016/j.neubiorev.2010.12.012
- PubMed
- Google Scholar
(2018) Information decomposition of target effects from multi-source interactions: perspectives on previous, current and future work
Entropy 20:307.
https://doi.org/10.3390/e20040307
- PubMed
- Google Scholar
(2020) Anterior insula reflects surprise in value-based decision-making and perception
NeuroImage 210:116549.
https://doi.org/10.1016/j.neuroimage.2020.116549
- PubMed
- Google Scholar
1. Luppi AI
2. Mediano PAM
3. Rosas FE
4. Holland N
5. Fryer TD
6. O’Brien JT
7. Rowe JB
8. Menon DK
9. Bor D
10. Stamatakis EA
(2022) A synergistic core for human brain evolution and cognition
Nature Neuroscience 25:771–782.
https://doi.org/10.1038/s41593-022-01070-0
- PubMed
- Google Scholar
(2024) Information decomposition and the informational architecture of the brain
Trends in Cognitive Sciences 28:352–368.
https://doi.org/10.1016/j.tics.2023.11.005
- Google Scholar
1. Maris E
2. Oostenveld R
(2007) Nonparametric statistical testing of EEG- and MEG-data
Journal of Neuroscience Methods 164:177–190.
https://doi.org/10.1016/j.jneumeth.2007.03.024
- PubMed
- Google Scholar
1. Matsumoto M
2. Hikosaka O
(2009) Two types of dopamine neuron distinctly convey positive and negative motivational signals
Nature 459:837–841.
https://doi.org/10.1038/nature08028
- PubMed
- Google Scholar
1. McGill W
(1954) Multivariate information transmission
Transactions of the IRE Professional Group on Information Theory 4:93–111.
https://doi.org/10.1109/TIT.1954.1057469
- Google Scholar
(2008) Dynamic population coding of category information in inferior temporal and prefrontal cortex
Journal of Neurophysiology 100:1407–1419.
https://doi.org/10.1152/jn.90248.2008
- PubMed
- Google Scholar
1. Michelmann S
2. Price AR
3. Aubrey B
4. Strauss CK
5. Doyle WK
6. Friedman D
7. Dugan PC
8. Devinsky O
9. Devore S
10. Flinker A
11. Hasson U
12. Norman KA
(2021) Moment-by-moment tracking of naturalistic learning and its underlying hippocampo-cortical interactions
Nature Communications 12:5394.
https://doi.org/10.1038/s41467-021-25376-y
- PubMed
- Google Scholar
(2024) Cognition is an emergent property
Current Opinion in Behavioral Sciences 57:101388.
https://doi.org/10.1016/j.cobeha.2024.101388
- Google Scholar
1. Mitra PP
2. Pesaran B
(1999) Analysis of dynamic brain imaging data
Biophysical Journal 76:691–708.
https://doi.org/10.1016/S0006-3495(99)77236-X
- PubMed
- Google Scholar
1. Monosov IE
2. Hikosaka O
(2012) Regionally distinct processing of rewards and punishments by the primate ventromedial prefrontal cortex
The Journal of Neuroscience 32:10318–10330.
https://doi.org/10.1523/JNEUROSCI.1801-12.2012
- PubMed
- Google Scholar
1. Morrison SE
2. Salzman CD
(2009) The convergence of information about rewarding and aversive stimuli in single neurons
The Journal of Neuroscience 29:11471–11483.
https://doi.org/10.1523/JNEUROSCI.1815-09.2009
- PubMed
- Google Scholar
1. Mukamel R
2. Gelbard H
3. Arieli A
4. Hasson U
5. Fried I
6. Malach R
(2005) Coupling between neuronal firing, field potentials, and FMRI in human auditory cortex
Science 309:951–954.
https://doi.org/10.1126/science.1110913
- PubMed
- Google Scholar
(2005) Hemodynamic signals correlate tightly with synchronized gamma oscillations
Science 309:948–951.
https://doi.org/10.1126/science.1110948
- PubMed
- Google Scholar
(2019) Synergistic coding of visual information in columnar networks
Neuron 104:402–411.
https://doi.org/10.1016/j.neuron.2019.07.006
- PubMed
- Google Scholar
1. Nir Y
2. Fisch L
3. Mukamel R
4. Gelbard-Sagiv H
5. Arieli A
6. Fried I
7. Malach R
(2007) Coupling between neuronal firing rate, gamma LFP, and BOLD fMRI is related to interneuronal correlations
Current Biology 17:1275–1285.
https://doi.org/10.1016/j.cub.2007.06.066
- PubMed
- Google Scholar
(2024) The tip of the iceberg: a call to embrace anti-localizationism in human neuroscience research
Imaging Neuroscience 2:1–10.
https://doi.org/10.1162/imag_a_00138
- Google Scholar
(2001) Abstract reward and punishment representations in the human orbitofrontal cortex
Nature Neuroscience 4:95–102.
https://doi.org/10.1038/82959
- PubMed
- Google Scholar
(2004) Dissociable roles of ventral and dorsal striatum in instrumental conditioning
Science 304:452–454.
https://doi.org/10.1126/science.1094285
- PubMed
- Google Scholar
1. Ohnuki T
2. Osako Y
3. Manabe H
4. Sakurai Y
5. Hirokawa J
(2021) Over-representation of fundamental decision variables in the prefrontal cortex underlies decision bias
Neuroscience Research 173:1–13.
https://doi.org/10.1016/j.neures.2021.07.002
- PubMed
- Google Scholar
1. Ongür D
2. Price JL
(2000) The organization of networks within the orbital and medial prefrontal cortex of rats, monkeys and humans
Cerebral Cortex 10:206–219.
https://doi.org/10.1093/cercor/10.3.206
- PubMed
- Google Scholar
(2017) Flexible information routing by transient synchrony
Nature Neuroscience 20:1014–1022.
https://doi.org/10.1038/nn.4569
- PubMed
- Google Scholar
(2009) Pharmacological modulation of subliminal learning in Parkinson’s and Tourette’s syndromes
PNAS 106:19179–19184.
https://doi.org/10.1073/pnas.0904035106
- PubMed
- Google Scholar
1. Palminteri S
2. Justo D
3. Jauffret C
4. Pavlicek B
5. Dauta A
6. Delmaire C
7. Czernecki V
8. Karachi C
9. Capelle L
10. Durr A
11. Pessiglione M
(2012) Critical roles for anterior insula and dorsal striatum in punishment-based avoidance learning
Neuron 76:998–1009.
https://doi.org/10.1016/j.neuron.2012.10.017
- PubMed
- Google Scholar
Website
(2015) Contextua modulation of value signals in reward and punishment learning. Nat Commun
Accessed January 2, 2019.

http://www.nature.com/articles/ncomms9096
1. Palminteri S
2. Pessiglione M
(2017) Opponent brain systems for reward and punishment learning
Decision Neuroscience 2017:291–303.
https://doi.org/10.1016/B978-0-12-805308-9.00023-3
- Google Scholar
1. Panzeri S
2. Macke JH
3. Gross J
4. Kayser C
(2015) Neural population coding: combining insights from microscopic and mass signals
Trends in Cognitive Sciences 19:162–172.
https://doi.org/10.1016/j.tics.2015.01.002
- PubMed
- Google Scholar
1. Panzeri S
2. Moroni M
3. Safaai H
4. Harvey CD
(2022) The structures and functions of correlations in neural population codes
Nature Reviews. Neuroscience 23:551–567.
https://doi.org/10.1038/s41583-022-00606-4
- PubMed
- Google Scholar
(2017) Neurons along the auditory pathway exhibit a hierarchical organization of prediction error
Nature Communications 8:2148.
https://doi.org/10.1038/s41467-017-02038-6
- PubMed
- Google Scholar
(2017) Mixed selectivity morphs population codes in prefrontal cortex
Nature Neuroscience 20:1770–1779.
https://doi.org/10.1038/s41593-017-0003-2
- PubMed
- Google Scholar
Book
1. Percival DB
2. Walden AT
(1993) Spectral Analysis for Physical Applications
cambridge university press.
https://doi.org/10.1017/CBO9780511622762
- Google Scholar
(2006) Dopamine-dependent prediction errors underpin reward-seeking behaviour in humans
Nature 442:1042–1045.
https://doi.org/10.1038/nature05051
- PubMed
- Google Scholar
1. Pessiglione M
2. Delgado MR
(2015) The good, the bad and the brain: neural correlates of appetitive and aversive values underlying decision making
Current Opinion in Behavioral Sciences 5:78–84.
https://doi.org/10.1016/j.cobeha.2015.08.006
- PubMed
- Google Scholar
1. Petersen SE
2. Sporns O
(2015) Brain networks and cognitive architectures
Neuron 88:207–219.
https://doi.org/10.1016/j.neuron.2015.09.027
- PubMed
- Google Scholar
(2010) Appetitive and aversive goal values are encoded in the medial orbitofrontal cortex at the time of decision making
The Journal of Neuroscience 30:10799–10808.
https://doi.org/10.1523/JNEUROSCI.0788-10.2010
- PubMed
- Google Scholar
1. Reid AT
2. Headley DB
3. Mill RD
4. Sanchez-Romero R
5. Uddin LQ
6. Marinazzo D
7. Lurie DJ
8. Valdés-Sosa PA
9. Hanson SJ
10. Biswal BB
11. Calhoun V
12. Poldrack RA
13. Cole MW
(2019) Advancing functional connectivity research from association to causation
Nature Neuroscience 22:1751–1760.
https://doi.org/10.1038/s41593-019-0510-4
- PubMed
- Google Scholar
Book
(1972)
Classical Conditioning II: Current Research and Theory

Appleton-Century-Crofts.
- Google Scholar
1. Rigotti M
2. Barak O
3. Warden MR
4. Wang XJ
5. Daw ND
6. Miller EK
7. Fusi S
(2013) The importance of mixed selectivity in complex cognitive tasks
Nature 497:585–590.
https://doi.org/10.1038/nature12160
- PubMed
- Google Scholar
1. Saez I
2. Lin J
3. Stolk A
4. Chang E
5. Parvizi J
6. Schalk G
7. Knight RT
8. Hsu M
(2018) Encoding of multiple reward-related computations in transient and sustained high-frequency activity in human oFC
Current Biology 28:2889–2899.
https://doi.org/10.1016/j.cub.2018.07.045
- PubMed
- Google Scholar
(2018) Coherent encoding of subjective spatial position in visual cortex and hippocampus
Nature 562:124–127.
https://doi.org/10.1038/s41586-018-0516-1
- PubMed
- Google Scholar
1. Salinas E
2. Sejnowski TJ
(2001) Correlated neuronal activity and the flow of neural information
Nature Reviews. Neuroscience 2:539–550.
https://doi.org/10.1038/35086012
- PubMed
- Google Scholar
1. Schneider M
2. Broggini AC
3. Dann B
4. Tzanou A
5. Uran C
6. Sheshadri S
7. Scherberger H
8. Vinck M
(2021) A mechanism for inter-areal coherence through communication based on connectivity and oscillatory power
Neuron 109:4050–4067.
https://doi.org/10.1016/j.neuron.2021.09.037
- PubMed
- Google Scholar
(2003) Synergy, redundancy, and independence in population codes
The Journal of Neuroscience 23:11539–11553.
https://doi.org/10.1523/JNEUROSCI.23-37-11539.2003
- PubMed
- Google Scholar
1. Schreiber T
(2000) Measuring information transfer
Physical Review Letters 85:461–464.
https://doi.org/10.1103/PhysRevLett.85.461
- PubMed
- Google Scholar
(1997) A neural substrate of prediction and reward
Science 275:1593–1599.
https://doi.org/10.1126/science.275.5306.1593
- PubMed
- Google Scholar
(2005) Opponent appetitive-aversive neural processes underlie predictive learning of pain relief
Nature Neuroscience 8:1234–1240.
https://doi.org/10.1038/nn1527
- PubMed
- Google Scholar
1. Shine JM
2. Bissett PG
3. Bell PT
4. Koyejo O
5. Balsters JH
6. Gorgolewski KJ
7. Moodie CA
8. Poldrack RA
(2016) The Dynamics of Functional Brain Networks: Integrated Network States during Cognitive Task Performance
Neuron 92:544–554.
https://doi.org/10.1016/j.neuron.2016.09.018
- PubMed
- Google Scholar
1. Sporns O
(2013) Network attributes for segregation and integration in the human brain
Current Opinion in Neurobiology 23:162–171.
https://doi.org/10.1016/j.conb.2012.11.015
- PubMed
- Google Scholar
(2013) A causal link between prediction errors, dopamine neurons and learning
Nature Neuroscience 16:966–973.
https://doi.org/10.1038/nn.3413
- PubMed
- Google Scholar
(2019) Distributed coding of choice, action and engagement across the mouse brain
Nature 576:266–273.
https://doi.org/10.1038/s41586-019-1787-x
- PubMed
- Google Scholar
1. Stokes MG
2. Kusunoki M
3. Sigala N
4. Nili H
5. Gaffan D
6. Duncan J
(2013) Dynamic coding for cognitive control in prefrontal cortex
Neuron 78:364–375.
https://doi.org/10.1016/j.neuron.2013.01.039
- PubMed
- Google Scholar
(2014) Reward value comparison via mutual inhibition in ventromedial prefrontal cortex
Neuron 82:1357–1366.
https://doi.org/10.1016/j.neuron.2014.04.032
- PubMed
- Google Scholar
Book
1. Sutton RS
2. Barto AG
(2018)
Reinforcement Learning: An Introduction

MIT press.
- Google Scholar
(2021) An engram of intentionally forgotten information
Nature Communications 12:6443.
https://doi.org/10.1038/s41467-021-26713-x
- PubMed
- Google Scholar
1. Thiebaut de Schotten M
2. Forkel SJ
(2022) The emergent properties of the connected brain
Science 378:505–510.
https://doi.org/10.1126/science.abq2591
- PubMed
- Google Scholar
1. Thorndike EL
(1898) Animal intelligence: an experimental study of the associative processes in animals
The Psychological Review 2:i–109.
https://doi.org/10.1037/h0092987
- Google Scholar
1. Tom SM
2. Fox CR
3. Trepel C
4. Poldrack RA
(2007) The neural basis of loss aversion in decision-making under risk
Science 315:515–518.
https://doi.org/10.1126/science.1134239
- PubMed
- Google Scholar
(2022) Large-scale neural recordings call for new insights to link brain and behavior
Nature Neuroscience 25:11–19.
https://doi.org/10.1038/s41593-021-00980-9
- PubMed
- Google Scholar
(2001) The brainweb: phase synchronization and large-scale integration
Nature Reviews. Neuroscience 2:229–239.
https://doi.org/10.1038/35067550
- PubMed
- Google Scholar
(2023) Information-processing dynamics in neural networks of macaque cerebral cortex reflect cognitive state and behavior
PNAS 120:e2207677120.
https://doi.org/10.1073/pnas.2207677120
- PubMed
- Google Scholar
1. Vicente R
2. Wibral M
3. Lindner M
4. Pipa G
(2011) Transfer entropy--a model-free measure of effective connectivity for the neurosciences
Journal of Computational Neuroscience 30:45–67.
https://doi.org/10.1007/s10827-010-0262-3
- PubMed
- Google Scholar
(2006) Visual grouping and the focusing of attention induce gamma-band oscillations at different frequencies in human magnetoencephalogram signals
Journal of Cognitive Neuroscience 18:1850–1862.
https://doi.org/10.1162/jocn.2006.18.11.1850
- PubMed
- Google Scholar
(2023) Principles of large-scale neural interactions
Neuron 111:987–1002.
https://doi.org/10.1016/j.neuron.2023.03.015
- PubMed
- Google Scholar
1. Voitov I
2. Mrsic-Flogel TD
(2022) Cortical feedback loops bind distributed representations of working memory
Nature 608:381–389.
https://doi.org/10.1038/s41586-022-05014-3
- PubMed
- Google Scholar
(2017) Decision-making impairments following insular and medial temporal lobe resection for drug-resistant epilepsy
Social Cognitive and Affective Neuroscience 12:128–137.
https://doi.org/10.1093/scan/nsw152
- PubMed
- Google Scholar
1. Wang R
2. Liu M
3. Cheng X
4. Wu Y
5. Hildebrandt A
6. Zhou C
(2021) Segregation, integration, and balance of large-scale resting brain networks configure different cognitive abilities
PNAS 118:e2022288118.
https://doi.org/10.1073/pnas.2022288118
- PubMed
- Google Scholar
1. Watkins CJCH
2. Dayan P
(1992) Q-learning
Machine Learning 8:279–292.
https://doi.org/10.1007/BF00992698
- Google Scholar
(2017) Partial information decomposition as a unified approach to the specification of neural goal functions
Brain and Cognition 112:25–38.
https://doi.org/10.1016/j.bandc.2015.09.004
- PubMed
- Google Scholar
Book
1. Wiener N
(1956)
The Theory of Prediction

Mod Math Eng.
- Google Scholar
Preprint
1. Williams PL
2. Beer RD
(2010) Nonnegative decomposition of multivariate information
arXiv.

http://arxiv.org/abs/1004.2515
- Google Scholar
(2006) Dissociable systems for gain- and loss-related value predictions and errors of prediction in the human brain
The Journal of Neuroscience 26:9530–9537.
https://doi.org/10.1523/JNEUROSCI.2915-06.2006
- PubMed
- Google Scholar

Article and author information

Author details

Etienne Combrisson

Institut de Neurosciences de la Timone, Aix Marseille Université, Marseille, France

Contribution
Conceptualization, Data curation, Software, Formal analysis, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing

For correspondence
e.combrisson@gmail.com

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-7362-3247
Ruggero Basanisi

Institut de Neurosciences de la Timone, Aix Marseille Université, Marseille, France

Contribution
Software, Methodology

Competing interests
No competing interests declared
Maelle CM Gueguen

Univ. Grenoble Alpes, Inserm, U1216, Grenoble Institut Neurosciences, Grenoble, France

Contribution
Conceptualization, Data curation

Competing interests
No competing interests declared
Sylvain Rheims

Department of Functional Neurology and Epileptology, Hospices Civils de Lyon and University of Lyon, Lyon, France

Contribution
Data curation

Competing interests
No competing interests declared
Philippe Kahane

Univ. Grenoble Alpes, Inserm, U1216, CHU Grenoble Alpes, Grenoble Institut Neurosciences, Grenoble, France

Contribution
Data curation

Competing interests
No competing interests declared
Julien Bastin

Univ. Grenoble Alpes, Inserm, U1216, Grenoble Institut Neurosciences, Grenoble, France

Contribution
Data curation, Supervision, Funding acquisition, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-0533-7564
Andrea Brovelli

Institut de Neurosciences de la Timone, Aix Marseille Université, Marseille, France

Contribution
Conceptualization, Resources, Software, Supervision, Funding acquisition, Investigation, Methodology, Writing - original draft, Project administration, Writing - review and editing

For correspondence
andrea.brovelli@univ-amu.fr

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-5342-1330

Funding

Agence Nationale de la Recherche (ANR-18-CE28-0016)

Etienne Combrisson
Ruggero Basanisi
Julien Bastin
Andrea Brovelli

Agence Nationale de la Recherche (ANR-17-CE37-0018)

Maelle CM Gueguen
Sylvain Rheims
Philippe Kahane
Julien Bastin

Agence Nationale de la Recherche (ANR- 13-TECS-0013)

Maelle CM Gueguen
Sylvain Rheims
Philippe Kahane
Julien Bastin

HORIZON EUROPE Framework Programme (604102)

Julien Bastin

HORIZON EUROPE Framework Programme (945539)

Etienne Combrisson
Andrea Brovelli

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We would like to express our gratitude to Benjamin Morillon, Manuel R Mercier, and Stefano Palminteri for their valuable comments on an earlier draft of this manuscript and for their feedback on our responses to reviewer comments.

Ethics

All patients gave written informed consent and the study received approval from the ethics committee (CPP 09-CHUG-12, study 0907) and from a competent authority (ANSM no: 2009-A00239-48).

Version history

Preprint posted: May 2, 2023 (view preprint)
Sent for peer review: October 2, 2023
Preprint posted: November 9, 2023 (view preprint)
Preprint posted: June 13, 2024 (view preprint)
Version of Record published: June 28, 2024 (version 1)

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.92938. This DOI represents all versions, and will always resolve to the latest one.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.