Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Research Articles, Behavioral/Cognitive

Pain-Related Expectation and Prediction Error Signals in the Anterior Insula Are Not Related to Aversiveness

Sepideh Fazeli and Christian Büchel
Journal of Neuroscience 18 July 2018, 38 (29) 6461-6474; DOI: https://doi.org/10.1523/JNEUROSCI.0671-18.2018
Sepideh Fazeli
Department of Systems Neuroscience, University Medical Center Hamburg-Eppendorf, Hamburg, Germany 20246
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christian Büchel
Department of Systems Neuroscience, University Medical Center Hamburg-Eppendorf, Hamburg, Germany 20246
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Christian Büchel
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

The anterior insula has repeatedly been linked to the experience of aversive stimuli, such as pain. Previously, we showed that the anterior insula is involved in the integration of pain intensity and its prior expectation. However, it is unclear whether this integration occurs by a pain-specific expectation or a more general expectation of an aversive event. To dissociate these possibilities, we conducted an experiment using painful stimuli and aversive pictures with three levels of aversiveness on human male volunteers. Stimuli were preceded by a probabilistic, combined modality and intensity cue in a full factorial design. Subjective ratings of pain intensity and skin conductance responses were best explained by a combination of actual pain intensity and expected pain intensity. In addition, using fMRI, we investigated the neuronal implementation of the integration of prior expectation and pain intensity. Similar to subjective ratings and autonomic responses, the dorsal anterior insula represented pain intensity and expectations. The ventral anterior insula additionally represented the absolute difference of the two terms (i.e., the prediction error). The posterior insula only represented pain intensity. Importantly, the pattern observed in the anterior insula was only present if the cued modality was correct (i.e., expect pain); in case of an incorrect modality cue (i.e., expect aversive picture), the ventral anterior insula simply represented pain intensity. The stimulus expectation and prediction error specificity in the ventral anterior insula indicates the integration of expectation with painful stimuli in this area. Importantly, this pattern cannot be explained by aversiveness.

SIGNIFICANCE STATEMENT The anterior insula has been shown to integrate pain intensity and their expectation. However, it is unclear whether this integration is pain-specific or related more generally to an aversive event. To address this, we combined painful stimuli and aversive pictures with three levels of aversiveness. The ventral anterior insula represented pain intensity, expectation, and their absolute difference (i.e., the prediction error). Importantly, this pattern was only observed if the cued modality was correct. In case of an incorrect modality cue, this area simply represented as pain intensity. The stimulus expectation and prediction error specificity in the ventral anterior insula indicates the integration of expectation with painful stimuli in this area. Importantly, this pattern cannot be explained by aversiveness.

  • expectation
  • human fMRI
  • pain
  • predictive coding
  • somatosensory perception

Introduction

Accurate and fast perception of pain is crucial for an organism to avoid tissue damage. In many cases, sensory cues precede pain and allow the organism to generate an expectation or prediction. The influence of such expectations on pain perception has been established in basic cued pain paradigms (Koyama et al., 2005; Keltner et al., 2006; Wiech et al., 2008; Atlas et al., 2010) and in the context of placebo and nocebo effects (Tracey, 2010; Häuser et al., 2012). It has previously been speculated that predictive coding, which integrates predictions and violation of predictions (i.e., prediction errors), could be a suitable mechanism mediating the effect of expectations on pain (Büchel et al., 2014; Geuter et al., 2017). In this framework, expectations (i.e., predictions) and the mismatch between these expectations and the sensory inputs (i.e., prediction errors) are assumed to be combined to minimize future prediction errors by updating the expectations at multiple levels of the neuronal hierarchy. So far, this framework has mainly been applied to visual and auditory processing (Rao and Ballard, 1999; Sterzer et al., 2008; Egner et al., 2010; Okada et al., 2018) and has only recently been extended to pain perception (Geuter et al., 2017). However, pain is not unimodal but comprises other aspects, such as aversiveness, unpleasantness, and salience, which are shared by negative stimuli in other sensory modalities. Therefore, the observed expectation and prediction error signals in the anterior insula and the intensity representation in the posterior insula (Geuter et al., 2017) might not be related to pain per se but could be related to aversiveness in general. This alternative hypothesis can be tested by adding a second aversive stimulus modality. Consequently, we adapted our cued pain paradigm (Geuter et al., 2017) by adding aversive pictures as a second sensory modality. In both modalities, we presented aversive stimuli at three levels of intensity. Each stimulus was preceded by a combined modality and intensity cue. Cues and stimuli were combined in a full factorial design with factors of cue modality, cue intensity, stimulus modality, and stimulus intensity. This design allowed us to investigate whether expecting an aversive stimulus with a specific intensity only affects the neuronal activity in response to a stimulus in the cued modality or whether this influence also generalizes to a stimulus in the other modality.

Within the framework of predictive coding, the prediction error represents the “surprise” induced by a mismatch between the sensory signals encountered and those predicted. If the previously observed prediction error signal in the anterior insula is related to aversiveness in general, we should observe such a signal if the cued intensity differs from the perceived intensity, regardless of the cued modality: a prediction error for a high-intensity painful stimulus, when a low-intensity aversive picture or a low-intensity painful stimulus was expected. On the contrary, if the anterior insula does not represent aversiveness in general, we should only observe prediction errors if the intensity cue is incongruent with the experienced intensity, but only if the modality cue matches: if a high intensity painful stimulus is perceived, when a low-intensity pain stimulus is cued, but not when a low-intensity aversive picture cue has been presented. Therefore, investigating the modality cue specificity of expectations and prediction errors for pain can potentially shed light on the characteristic of the pain processing system and a possible implementation of predictive coding therein.

Materials and Methods

Participants

Thirty-two healthy human male volunteers (average age: 26 years; range: 19–38 years) successfully participated in the study and were paid as compensation for their participation. Exclusion criteria were neurological, psychiatric, dermatological diseases, pain conditions, current medication, or substance abuse. All volunteers gave their informed consent. The study was approved by the Ethics board of the Hamburg Medical Association.

Stimuli and task

The stimuli were either thermal or visual stimuli. Thermal stimulation was performed using a 30 × 30 mm2 Peltier thermode (CHEPS Pathway, Medoc) at three different intensities: warm (42°C), medium painful (46°C), and highly painful (48°C) to the left radial forearm at different skin patches in each experimental block to avoid sensitization. The stimulation pattern was chosen based on the distribution of average human pain thresholds to have three distinguishable levels of thermal stimulation. The visual stimuli were pictures chosen from the International Affective Picture System (IAPS) (Lang et al., 2005) database at three levels of aversiveness. The classification of IAPS pictures was performed manually to divide the reported valance spectrum into six separate groups, where three of them were assigned to the training set and the other three were assigned to the fMRI set. The pictures in the fMRI set had three levels of aversiveness of which the low aversive category had aversiveness values of 2.02 ± 0.05 (mean ± SE), the medium aversive category had aversiveness values of 4.06 ± 0.02 (mean ± SE), and the high aversive category had aversiveness values of 5.23 ± 0.01 (mean ± SE). Before each picture or heat stimulus, a visual cue was presented. The color of the cue (triangle) indicated the modality of the stimulus (yellow for picture and red for heat). The size of the cue indicated its intensity (small, medium, and large triangle for low, medium, and high intensities, respectively).

At the beginning of each trial (Fig. 1), the cue was presented for 500 ms, indicating the modality and intensity of the following stimulus. The cue validity for the stimulus modality was 70%. Cue validity with respect to intensity was 60% (for details of all cue stimulus combinations and frequency, see Fig. 1F). After a jittered blank period of 300–700 ms, volunteers received the visual or heat stimulus with a duration of 2 s. After stimulus presentation, volunteers had to rate the intensity of the aversiveness of the picture or the heat stimulus within a time window of 2 s on a 4 point rating scale (1 labeled as “neutral” and 4 as “strong”). A fixation cross was presented throughout the trial. In addition, we added 4 catch trials to each block, in which volunteers were asked to report the preceding cue within 8 s. No stimulation was given in these trials. Volunteers were rewarded with 50 cents for a correct answer in each catch trial. Finally, a variable intertrial interval (2.5 ± 0.5 s) followed each trial, during which only a fixation cross was displayed.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Experimental design and predicted response patterns for INT, INT+EXP, INT+PE, and INT+EXP+PE models. A, In the INT model, the response to the stimulus is only driven by its intensity. B, In the INT+EXP model, the stimulus response is the weighted sum of the stimulus intensity (left) and the stimulus expectation (middle). C, In the INT+PE model, the stimulus response is the weighted sum of the stimulus intensity (left) and the stimulus prediction error (middle). D, For the INT+EXP+PE model, the stimulus response is the weighted sum of the stimulus intensity (left), the stimulus expectation (second from left), and the prediction error (third from left). Right, Examples shown here assume an equal contribution of the expectation and the prediction error signals to the stimulus response. Different contributions of the expectation or prediction error can lead to a different pattern for the stimulus response. E, Trial structure: At the beginning of each trial, a central cue indicated the probability of the intensity and the modality of the upcoming stimulus. Then after a blank screen with variable duration, volunteers received the stimulus (2 s), followed by the rating (2 s). The intertrial intervals were randomly jittered between 2 and 3 s. Volunteers were asked to fixate a central-fixation-dot (except for the rating duration). F, Table illustrating the cue-stimulus contingencies.

Procedure

The experiment consisted of two sessions: a behavioral training session and an fMRI session. Before the training session, volunteers were informed about the procedure and gave their written informed consent. They were informed about the experimental paradigm and the nature of the cues and the contingencies between the cues and the stimuli to avoid learning effects during the main fMRI session. The training session contained up to three blocks and was terminated once volunteers had learned the cue regularities. The termination criterion was defined as detecting >75% catch trials and a significant discrimination (repeated-measures ANOVA [rmANOVA]: p < 0.05) of the three levels of pain and aversive picture intensities. The fMRI session consisted of four blocks. Each block consisted of 130 trials (including four catch trials) and lasted ∼15 min. The trial order within each block was pseudo-randomized, and the order of blocks was randomized across volunteers. During each block, fMRI signal and skin conductance responses (SCRs) were measured. A high-resolution anatomical image for each volunteer was acquired at the end of the fMRI session. The whole fMRI experiment lasted for ∼2 h.

Data acquisition

fMRI data were recorded on a Trio 3 tesla system equipped with a 32-channel head coil (Siemens). Forty transversal slices (voxel size 2 × 2 × 2 mm, 1 mm interslice gap) were obtained within each volume using a T2*-sensitive single-shot EPI sequence (TR = 1.24 s, TE = 26 ms, flip angle: 60°, FOV: 220 × 220 mm) with parallel imaging (GRAPPA, in-plane acceleration factor 2) (Griswold et al., 2002) and simultaneous multislice acquisitions (“multiband,” slice acceleration factor 2) (Feinberg et al., 2010; Moeller et al., 2010; Xu et al., 2013) as described previously (Setsompop et al., 2012). The corresponding image reconstruction algorithm was provided by the University of Minnesota Center for Magnetic Resonance Research. Slices were tilted ∼30° relative to the AC-PC line to improve coverage in the brainstem. Additionally, T1-weighted structural images (1 × 1 × 1 mm resolution) were acquired using an MPRAGE sequence (TR = 2300 ms, TE = 2.98 ms, flip-angle = 9°). Stimulus presentation, response recording, and timing relative to MR data acquisition were done using Psychophysics Toolbox 3 (http://www.psychtoolbox.org). Skin conductance was recorded using an EDA100C MRI system (Biopac Systems), amplified using an analog amplifier (MP150, Biopac Systems), and sampled at 100 Hz using a CED 1401 analog-digital converter (Cambridge Electronic Design). SCR electrodes were connected to the palm (thenar and hypothenar eminences) of the left hand. During the study, lights in the MRI room were dimmed, and luminance was kept constant across volunteers.

Data analyses

SCRs.

The phasic skin conductance drive locked to the onset of each stimulus was estimated from the raw SCRs in each trial using an established deconvolution method as implemented in ledalab (Benedek and Kaernbach, 2010). The resulting phasic responses were then separately averaged for different conditions within a temporal window from 2 to 5 s after stimulus onset for the painful stimuli and within a window from 0 to 3 s for the aversive pictures, which was selected to cover the peak of the averaged SCR signal in each modality across volunteers.

fMRI data preprocessing.

Functional imaging data were analyzed using MATLAB (version 2014b, The MathWorks) and SPM12 (Wellcome Trust Centre for Neuroimaging, London). The first five volumes of each functional block were discarded. The remaining images were spatially realigned to correct for motion artifacts and were then spatially normalized using DARTEL, a high-dimensional warping algorithm implemented in SPM (Ashburner, 2007). The functional images were spatially smoothed using a Gaussian kernel with an FWHM of 6 mm. First-level models included separate regressors using a boxcar with a duration of 2 s for each of the 36 experimental conditions: 6 stimuli (3 levels of pain + 3 levels of aversive images) × 3 intensity cues × 2 modality cues. Additional regressors with a duration of 2 s modeled the rating period. Another regressor (boxcar, duration 9 s) modeled the catch trials. In addition, the model also included six motion parameters estimated during realignment. All regressors were convolved with a canonical hemodynamic response function as implemented in SPM. The first-level contrast estimates for the 18 pain conditions were the basis for our computational models. Model estimation used a linear mixed-effect model as implemented in MATLAB (version 2014b, The MathWorks). Models were estimated for all gray-matter voxels. Gray-matter voxels were determined using the segmentation routine in SMP12 at the single volunteer level, then averaged across volunteers, and finally thresholded at an 80% gray matter probability.

Neurological pain signature (NPS) analyses.

The NPS is a multivariate pattern of brain activity with the aim of distinguishing experimental pain from other conditions (Wager et al., 2013; Krishnan et al., 2016). The NPS response mirrors subjective pain and has been suggested as a surrogate marker for heat pain intensity (Wager et al., 2013). We computed NPS expression values for each of the experimental conditions based on the regressors described above separately for each participant (Wager et al., 2011, 2013; Krishnan et al., 2016). The resulting NPS values were then the basis for our model comparison.

ROI masks and analysis.

According to the finding reported by Geuter et al. (2017), insular cortex is the area among pain-responsive regions representing expectation and prediction error signals. Based on these findings, we mainly focused our analysis on the insular cortex and its anatomically defined subdivisions as defined by the Destrieux atlas (Fischl et al., 2004; Destrieux et al., 2010). Using insular subdivisions, we tested our four models on the average parameter estimates extracted from voxels within these ROIs (for the detailed model comparison routine, see below). Even though the insular cortex is the focus of this study because of its key role in the integration of pain expectation as the internal model, and sensory signals as the physical input (Ostrowsky et al., 2002; Nieuwenhuys, 2012; Geuter et al., 2017), the processing of nociceptive information is not limited to the insular cortex. We therefore conducted additional ROI-based analyses in several other pain-related regions, including the parietal operculum (SII), postcentral sulcus (SI), frontal medial cortex, anterior cingulate (ACC), posterior cingulate (PCC), amygdala, hippocampus, thalamus, cerebellum, and periaqueductal gray (PAG). Anatomical masks were the same as used by Geuter et al. (2017), except for the PCC, frontal medial cortex, hippocampus, and cerebellum, which were generated by thresholding the probability map of the corresponding areas in the Harvard-Oxford Atlas (Desikan et al., 2006) at 50% probability (https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/Atlases). Similar to insular subdivisions, we tested our four models on the average parameter estimates from voxels within these masks.

Modeling behavioral, SCR, NPS, and voxelwise fMRI data.

Crossing our three levels of stimulus intensity, three levels of intensity cues, and two levels of modality cues resulted in 18 (3 × 3 × 2) different pain conditions, in which nine stimuli were preceded by a valid modality cue and nine stimuli were preceded by an invalid modality cue. For these valid and invalid cue conditions, subjective ratings, SCR, NPS, and voxelwise fMRI data were fitted separately to four different models. All models included a stimulus intensity term. The first model (intensity model [INT]) only contained the stimulus intensity term as the physical characteristic (intensity) of the stimulus (Fig. 1A). The second model (expectation model [INT+EXP]) additionally comprised an expectation term represented by the intensity cue (Figure 1B). The third model (prediction error model [INT+PE]), comprised two terms: the prediction error (i.e., the absolute difference between intensity cue and stimulus intensity) and stimulus intensity (Fig. 1C). This model was inspired by several studies reporting reduced neuronal activity in response to expected versus unexpected stimuli (Egner et al., 2010; Jiang et al., 2013). The final full model (INT+EXP+PE) contained both the expectation and the prediction error term, as well as the intensity term. Conceptually, this model assumes that the response to a painful stimulus is a weighted sum of its intensity value, the prediction error as well as the expectation (Fig. 1D). This model was inspired by studies showing the existence of both expectation and prediction error signals in different brain regions (Egner et al., 2010; Geuter et al., 2017). Prediction errors were coded in absolute differences, motivated by previous studies in the visual and auditory system (den Ouden et al., 2010; Todorovic et al., 2011; Todorovic and de Lange, 2012; Boll et al., 2013; Geuter et al., 2017). Furthermore, to account for the uneven “perceptual distance” between three levels of intensity, we added a second-order intensity term to all models. Model comparisons were performed on behavioral ratings, autonomous (SCR) data, as well as voxelwise BOLD signals and the NPS score. Model fit and parameter estimation were conducted using the FITLME function in MATLAB (version 2014b, The MathWorks).

Model comparison and statistical analysis.

To find the best of four alternative models, we used a nested-model comparison in a backward selection fashion. We started with the INT+EXP+PE model as the most comprehensive model and compared it with the INT+EXP and INT+PE models using likelihood-ratio tests. The likelihood-ratio test constructs a test statistic using the log-likelihood objective function evaluated under the unrestricted model parameter estimates (INT+EXP+PE) and the restricted model parameter estimates (INT+EXP or INT+PE). The test statistic follows a χ2 distribution with 1 df. Comparing the unrestricted INT+EXP+PE model with the two alternative-restricted INT+EXP and INT+PE models resulted in two p values. If both p values were below our threshold (e.g., p < 0.05), then INT+EXP+PE (unrestricted model) was selected as the best model. Otherwise, the model associated with the biggest p value was discarded, and the fit to the other restricted model was compared with the fit to the INT model as the simplest model using likelihood-ratio test. If the ensuing p value was below the threshold, then the intermediate model (in our case either INT+EXP or INT+PE) was the winning model. Otherwise, the INT model as the simplest model was determined as the best model. In summary, model selection proceeded according to the following algorithm:

  • 1: Fit the data to the full model (INT+EXP+PE).

  • 2: Compare the model with the two simpler models (INT+EXP and INT+ PE).

  • 3: Find the largest associated p value from the previous tests.

  • ≫ If the largest p value is above threshold, go to Step 4.

  • ≫ Otherwise, determine the full unrestricted model as the best model and stop.

  • 4: Compare the model with smallest associated p value with the INT model.

  • ≫ If the p value is above threshold, determine the INT model as the best model and stop.

  • ≫ Otherwise, determine the intermediate unrestricted model as the best model and stop.

The criteria used to select each model as the winning model are also illustrated in Figure 2. In addition, we used Bayesian model weights to verify the model comparison results of the likelihood-ratio test. The Bayesian model weight wM for the model M is given by the following: Embedded Image where BIC is the Bayesian Information Criterions computed using FITLME function in MATLAB (version 2014b, The MathWorks). Each wM is between 0 and 1, and the model with the biggest wM is the best model. Furthermore, given the computed BIC, the Bayes factor comparing two alternative model was obtained as follows: Embedded Image Therefore, the ratio of any two wM is equal to the corresponding Bayes factor comparing those models. For behavioral, SCR, and NPS data, the p value threshold was set at 0.01. At the fMRI voxel level, we used an ROI (insula) based correction for multiple comparisons using false discovery rate (FDR) correction at p = 0.05. FDR correction was performed using the FDR_BH function in MATLAB (version 2014b, The MathWorks). Moreover, we also tested for main effects and interactions across stimulus intensity, intensity cue, as well as modality cue using rmANOVA as implemented in MATLAB (version 2014b, The MathWorks). We also conducted paired t test, where all tests were two-sided unless otherwise stated. All statistical tests were conducted in MATLAB (version 2014b, The MathWorks).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Flowchart of the model selection procedure. A, The EXP+PE+INT model will be considered the winning model if the associated p values of the likelihood ratio tests are <0.05. B, The EXP+INT is considered the best model, if the p value for the likelihood ratio test of EXP+PE+INT against PE+INT is <0.05, but the p value for EXP+PE+INT versus EXP+INT is >0.05. C, The PE +INT model is considered the best model, if the p value for the likelihood ratio test of EXP+PE+INT versus EXP+INT is <0.05, and the p value for EXP+PE+INT versus EXP+INT is >0.05. D, If both p values of the likelihood ratio tests are >0.05, the INT model is considered the best model.

Results

Behavioral results

To validate our intensity manipulation for heat pain and to verify the discriminability between different levels of pain, we first tested for the main effect of stimulus intensity (Fig. 3A). Results regarding the aversive pictures are beyond the scope of this report but are depicted in Figure 3A for the sake of comparison. An rmANOVA, based on the behavioral data, revealed significant effects for the three levels of pain (F(2,62) = 736.4, p = 6 × 10−44, rmANOVA). Similarly, the main effects of intensity cue (F(2,62) = 30.7, p = 5 × 10−5, rmANOVA) and modality cue (F(1,31) = 10.85, p = 0.0025, rmANOVA) were significant. The two-way interaction between stimulus intensity and intensity cue (F(4,124) = 7.09, p = 3 × 10−5, rmANOVA) and the interaction between modality cue and intensity cue (F(2,62) = 11.45, p = 0.0001, rmANOVA) were also significant (Table 1). The rmANOVA results regarding the aversive pictures are summarized in Table 1.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Behavioral, autonomous responses and NPS to thermal pain and aversive pictures. A, Color bars represent the subjective ratings of warm (38°), medium-painful (46°), and highly painful (48°) stimuli (left) and subjective ratings for aversiveness of IAPS pictures (right). B, SCRs to the pain and aversive pictures. C, NPS responses to painful stimuli and aversive pictures. Error bars indicate SEM across volunteers.

View this table:
  • View inline
  • View popup
Table 1.

Main effect of stimulus intensity, intensity cue, modality cue, as well as their corresponding two-way and three-way interactions: subjective ratings, SCR, and NPS in both pain and picture conditions

In addition, we fitted the behavioral ratings to our four alternative models (INT, INT+EXP, INT+PE, and INT+EXP+PE). Backward model selection revealed that the INT+EXP model best explained the data in the conditions with a valid modality cue. In contrast, the INT model best explained the data in the invalid modality cue conditions. In addition to the likelihood-ratio test, BIC and Bayesian model weights were computed for each model to validate the consistency of our model selection routine. Comparing the four models in the valid modality cue conditions based on their corresponding BIC, INT+EXP+PE model showed the biggest weight (BIC = 144.2 and BICw = 0.6), even though the Bayes Factor comparing INT+EXP+PE and INT+EXP model was BFPC>EXP = 1.52, which indicates their difference is not substantial; and as the expectation model (INT+EXP) has less degree of freedom, it can be determined as the best model describing the behavioral rating data in the valid modality cue conditions (Lewandowsky and Farrell, 2010). For the invalid modality cue conditions, again the INT+EXP+PE model had the biggest weight (BIC = 144.2 and BICw = 0.29), but there was no substantial difference between INT+EXP+PE model and the other three models (BFPC>EXP = 1.18, BFPC>PE = 1.21, BFPC>IC = 1.34). Therefore, the INT model as the simplest model was identified as the best model describing the behavioral ratings in the invalid modality cue conditions. In summary, both model selection strategies identified the INT model as the best model describing the behavioral ratings in the invalid modality cue conditions and the INT+EXP model as the best model for the data in the conditions with a valid modality cue.

SCR

Analogous to the behavioral data, the main effects of stimulus intensity (F(2,62) = 53.2, p = 3 × 10−14, rmANOVA), intensity cue (F(2,62) = 8.48, p = 0.006, rmANOVA), and modality cue (F(1,31) = 12.71, p = 0.001, rmANOVA) were significant. Similar to the ratings, the two-way interaction between stimulus intensity and intensity cue (F(4,124) = 12.33, p = 2 × 10−8, rmANOVA) as well as the interaction between the modality cue and the intensity cue (F(2,62) = 8.4, p = 0.0006, rmANOVA) were significant (Table 1; Fig. 3B). The best model describing the SCR data were identified using likelihood-ratio test in both valid and invalid modality cue conditions. As in the behavioral ratings, INT+EXP was the winning model in explaining the SCR data in the valid cue conditions. The INT model was the best model in the invalid modality cue conditions. This result was verified using Bayesian model weights for the model comparisons for the valid and invalid modality cue conditions. For valid modality cue conditions, the INT+EXP model had the biggest BICw of 0.78, whereas the BICw of INT+EXP+PE model was 0.21, the BICw of INT+PE model was 4 × 10−4, and the INT model had a BICw of 0.002. For the invalid modality cue conditions, INT+PE had the biggest BICw of 0.47, which was only 1.25 times bigger than the BICw of the INT model with the value of 0.37. INT+EXP+PE and INT+EXP models had BICws of 0.08 and 0.06, respectively. Therefore, the INT model can be considered as the best model in the invalid modality cue conditions.

fMRI data

In a first step, we analyzed fMRI data from all conditions using the NPS (Wager et al., 2013). For the NPS data, an rmANOVA showed a significant effect for the three levels of pain (F(2,62) = 122.65, p = 3 × 10−22, rmANOVA). Neither the main effect of intensity cue nor modality cue was significant. The two-way interaction between modality cue and intensity cue was significant (F(2,62) = 5.1, p = 0.008, rmANOVA); however, the other two-way interactions were not significant (Fig. 3C; Table 1). As for the behavioral data, we fitted the NPS-derived scores to all our four models. Analogous to the behavioral ratings, the INT+EXP model best explained the data for the valid modality cue conditions using hierarchical model comparison based on the likelihood-ratio test. In the invalid modality cue conditions, the INT model was best in explaining the NPS data based on the nested likelihood-ratio test. This was confirmed by Bayesian model weights in the valid and invalid modality cue conditions. The INT+EXP model had the biggest BICw of 0.88 (with the Bayes Factor of 14.83 compared with the INT+EXP+PE model as the second wining model) for valid modality cue conditions. The Bayesian model weights for INT+EXP+PE, INT+PE, and INT models were 0.06, 0.004, and 0.05, respectively. In the invalid modality cue conditions, the INT model had the biggest BICw of 0.89 (with the Bayes Factor of 16.15 compared with the INT+PE model as the second wining model). The other alternative models of INT+EXP+PE, INT+EXP, and INT+PE had the BICw of 0.003, 0.053, and 0.055, respectively. As control, we also tested the NPS on picture stimuli (Table 1; Fig. 3C).

We next investigated the response patterns of voxels in subregions of insular cortex according to a surface atlas (Fischl et al., 2004; Destrieux et al., 2010). Regions in this atlas are defined according to sulcal depth (Fischl et al., 2004; Destrieux et al., 2010). To obtain an intuitive visualization of the cortical representation for each model, we color-coded surface vertices according to the best explanatory model (four different colors) on the inflated cortical surface. As shown in Figure 4, the spatial distribution of the insular subregions according to this atlas closely resembled the distribution of the winning models. Activity patterns in the anterior division of the circular insula (ventral part of anterior insula [vAI]) were best explained by the INT+EXP+PE model, activity patterns in the superior division of the circular insula (dorsal part of anterior Insula) were best explained by the INT+EXP model, and the inferior division of the circular insula showed activation patterns compatible with the INT model. To verify this, we averaged the contrast estimates in each anatomical subregion of the insula and performed a formal model comparison on this average activation pattern. For the valid modality cue conditions, the INT+EXP+PE model best accounted for the fMRI pattern in the anterior division of circular insula, the INT+EXP model best described the response pattern in the superior division of circular insula, and the INT model was the best model explaining fMRI activation patterns in the inferior division of circular insula (Fig. 4, inserted). For all invalid modality cue conditions, the INT model was the best model across all insular subdivisions. We also compared the individual contribution of the expectation and the prediction error terms in voxels best explained by the INT+EXP+PE model within the insula. For this purpose, we mapped the corresponding parameter estimates within the insular cortex on the inflated cortical surfaces (Fig. 5A,B). This analysis revealed that the prediction error parameter estimate and the expectation parameter estimate are not significantly different in both hemispheres. We also observed that the expectation weight was not significantly different in INT+EXP representative cluster compared with INT+EXP+PE voxels in both hemispheres (Fig. 5C).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Voxelwise representation of the best model explaining responses to thermal pain within insular cortex. A, Contrast estimates were fitted to four alternative models of INT, INT+EXP, INT+PE, and INT+EXP+PE (valid modality cue). The INT+EXP+PE model (red) best explains the data in the anterior division of the circular insula (vAI). Color coding is based on likelihood-ratio test model comparison (FDR-corrected, q < 0.05). q values are the adjusted p values found using an optimized FDR approach in voxels within the insula. The INT+EXP model (yellow) best explains the fMRI data in the superior division of the circular insula (dAI). Within the inferior division of circular insula, the INT model (blue) is the winning model. B, When the modality cue was invalid, the contrast estimates were best explained by the INT model (blue) in all regions. Error bars indicate SEM across volunteers.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Voxelwise map of the respective weights of the individual components of the INT+EXP+PE and INT+EXP models within insular cortex for the valid modality cue conditions. A, Weights for the expectation term are plotted for all INT+EXP+PE representative voxels in insular cortex (i.e., Fig. 4, red voxels). Bar graph represents the average of the expectation weight within voxels best explained by the INT+EXP+PE model. Error bars indicate SEM across volunteers. B, Voxelwise map for the weight of the prediction error term as well as its average for all voxels best explained by the INT+EXP+PE model. C, Map of the expectation term and its average across voxels best explained by the INT+EXP model (i.e., Fig. 4, yellow voxels).

In an exploratory fashion, we also performed formal model comparisons with the threshold of 0.01 outside the insular ROI (Fig. 6) on the surface and in all voxels, including subcortical areas. Here we additionally observed groups of voxels representing the INT+EXP model as the best model in the ACC, but only in valid modality cue conditions (Fig. 6A). Among the subcortical areas, the PAG showed activation patterns that were best explained by the INT+EXP model (Fig. 6C). As in cortical regions, this was only true if the modality cue was valid. Otherwise, intensity coding was identified as the best model (Fig. 6B,C).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Color map representing the best model for each voxel at an uncorrected threshold of 0.01 projected on flattened cortex and inflated brain surfaces. A, In addition to the insular clusters depicted in Figure 4, an INT+EXP cluster is also evident within the ACC and visual cortex. There is also a small cluster of the INT+PE model within the cingulate cortex (green). B, For the invalid modality cue, the INT model is the best model across all pain-responsive areas. C, Color-coded model selection for subcortical areas. In parts of the PAG, the INT+EXP model best described the data. For the invalid modality cue, the INT model best explains PAG responses.

In an explorative fashion, we also examined the averaged parameter estimates in other pain-related anatomically defined ROIs for the valid modality cue conditions. Based on rmANOVA, we found a significant main effect of stimulus intensity in all ROIs, including ACC, PCC, parietal operculum, postcentral, frontal medial cortex, hippocampus, amygdala, thalamus, cerebellum, as well as the PAG (for detailed results, see Table 2). Moreover, among all ROIs, a significant main effect of cue intensity was observed in the PAG (F(2,62) = 3.88, p = 0.02, rmANOVA). The two-way interaction between modality cue and intensity cue was significant in the ACC (LH: F(2,62) = 4.62, p = 0.01; RH: F(2,62) = 4.63, p = 0.01, rmANOVA). With regards to model selection, we used condition-specific parameter estimates for each voxel and subjected those parameter estimates to our formal model comparison. As for the behavioral data, the model best explaining the data were identified using hierarchical likelihood-ratio tests. The INT+EXP model was determined as the best model (for detailed results, see Fig. 7; Table 3) in the ACC and PAG. In all other regions, the INT model was the winning model.

View this table:
  • View inline
  • View popup
Table 2.

Main effect of stimulus intensity, intensity cue, modality cue, as well as their corresponding two-way and three-way interactions for all ROIs, including the SII, SI, frontal medial cortex, ACC, PCC, amygdala, hippocampus, thalamus, cerebellum, and PAGa

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Parameter estimates of ROIs. Mean of the parameter estimates (± SE) at three pain intensities of high, medium, and low in different intensity cue conditions for valid modality cue trials are plotted for left (LH) and right (RH) hemispheres (except for the midline PAG) in each ROI. The fit of the parameter estimates to the best model among the four nested models is also depicted. In the ACC and PAG, the best model was determined as the INT+EXP model. Otherwise, the INT model was the best explanatory model. Bayesian model weights (BICw) for all different models are illustrated in each ROI. Error bars indicate SEM across volunteers.

View this table:
  • View inline
  • View popup
Table 3.

Bayesian model weights for all ROIs, including the SII, SI, frontal medial cortex, ACC, PCC, amygdala, hippocampus, thalamus, cerebellum, and PAGa

Discussion

Our data show a clear discriminability of different levels of pain intensity based on behavioral ratings, autonomic measures, and a neural compound marker (NPS). Furthermore, we observed a strong effect of the intensity cue on pain perception, in agreement with previous findings (Koyama et al., 2005; Keltner et al., 2006; Atlas et al., 2010; Atlas and Wager, 2012; Wiech, 2016). Using computational modeling, activity patterns in the posterior insula were best explained by pure intensity coding, whereas activity in the dorsal anterior insula (dAI) was best explained by a model comprising the expected intensity. Finally, activity in the vAI was best explained by a model comprising stimulus intensity, cued intensity, and the prediction error (i.e., the difference between perceived and cued intensity) (Geuter et al., 2017). Finally, we observed a significant two-way interaction between modality cue and intensity cue, indicating a modality specific effect of expectation (Langner et al., 2011). Using formal modeling, we could reveal that the contribution of expectation in dAI and expectation plus prediction errors in vAI were only present for valid modality cues. In cases of an invalid modality cue, insular regions all showed simple intensity coding patterns.

Insular cortex

The insular cortex is crucial for pain processing (Ostrowsky et al., 2002; for review, see Nieuwenhuys (2012). It is ideally suited to represent pain intensity and related features because of its structural and functional organization (Barrett and Simmons, 2015; Nomi et al., 2017). While the traditional assumption (Cerliani et al., 2012; Cloutman et al., 2012; Jakab et al., 2012; Dennis et al., 2014) refers to an anterior–posterior dichotomy within insular cortex, more recent neuroimaging studies suggest a tripartite organization (Deen et al., 2011; Nomi et al., 2017).

In this tripartite organization, the posterior insula with its direct spinothalamocortical projections (originating in lamina I dorsal horn neurons) (Craig, 2002; Dum et al., 2009) as well as its intense connectivity to the somatosensory cortices (Wiech et al., 2014) makes it an ideal candidate for encoding the intensity of nociceptive thermal stimuli. Such sensory coding within the PI has previously been observed for both tonic and phasic pain (Segerdahl et al., 2015; Geuter et al., 2017) and even innocuous thermal stimuli (Davis et al., 1998). Our observation of a pure intensity coding pattern in the posterior insula confirms these findings.

The dAI is functionally connected to a cognitive control network (Dosenbach et al., 2007) and has stronger connections to posterior regions of the cingulate cortex (Taylor et al., 2009). Therefore, the dAI is presumably involved in expectation and attentional mechanisms (Deen et al., 2011). Moreover, the dAI as a part of cortical limbic areas with an agranular laminar structure is considered to encode expectation signals (Barrett and Simmons, 2015; Chanes and Barrett, 2016). This is in accordance with our findings showing a representation of stimulus intensity and expectation (INT+EXP model) in the dAI. Previous studies have also identified a role of this area in prestimulus effects and attentional modulation (Ploner et al., 2010; Wiech et al., 2010).

In contrast, the vAI exhibits strong connectivity to the posterior insula, ACC, thalamus, amygdala, and basal ganglia, which are also contributing to pain processing (DuPont et al., 2003; Baliki et al., 2009; Mutschler et al., 2009; Uddin et al., 2014). Such a strong connection of the vAI to these limbic structures makes it an ideal mediator between areas representing pain intensity and expectation to evaluate the expectation against the sensory input and estimating a prediction error (Downar et al., 2002; Seymour et al., 2004; Geuter et al., 2017). Our observed representation of stimulus intensity, expectation, and prediction error (INT+EXP+PE model) in the vAI is in agreement with these findings.

Expectation and prediction error in the valid modality

Predictions and prediction errors are key ingredients for perceptual theories, such as predictive coding (Knill and Pouget, 2004), and theories of brain function, such as the free energy model, including active inference (Knill and Pouget, 2004; Friston, 2005, 2008). In the predictive coding framework, the brain is relying on its internal model of the physical reality but tries to make it an adequate model by updating it based on the incoming data. Consequently, the integration of “prediction errors” into this model will change it over time to better reflect the physical reality. These concepts can account for extra-classical-receptive-field effects (Rao and Ballard, 1999), repetition suppression (Summerfield et al., 2008; Todorovic and de Lange, 2012), illusory contours formation (Dura-Bernal et al., 2011), and mismatch negativity (Ylinen et al., 2016). Such a mechanism was also proposed as a model to account for expectation effects in pain (Geuter et al., 2017) and placebo hypoalgesia (Büchel et al., 2014; Grahl et al., 2018). According to predictive coding in the case of pain, we hypothesize the presence of sensory specific representations of the expectation and prediction error terms at the hierarchical level where the intensity of pain is encoded similar to other modalities (Rao and Ballard, 1999; Friston, 2005). Crucially, this predicts that these representations are restricted to the case in which the modality cue is valid. As our second stimulus category comprised aversive pictures, this allows to test whether aversiveness, which is shared by pain and picture stimuli, can account for insula activation.

Importantly, our data show that the prediction error signal in the anterior insula is not generally related to aversiveness, which is shared by the pictures and the painful stimuli, but is limited to the pain expectation. This is indicated by the fact that we could only observe the “expectation” effect in the valid modality cue conditions. In case the modality cue was invalid (when volunteers expected an aversive picture but received pain), the intensity cue had no effect on the response patterns in the insula. In other words, an expectation that shares an aversiveness component (i.e., aversive pictures) is not encoded in these areas. This is further corroborated by the absence of a prediction error in the insula in this case. As a consequence, our data suggest a more specific role for the vAI, representing expectations and prediction errors related to pain. This indicates that this area is part of the hierarchy where pain intensity integration occurs and can be seen in analogy to other stimulus modalities, in which feature-selective predictive coding within the fusiform face area in response to faces, and within the parahippocampal place area in response to houses (den Ouden et al., 2010; Egner et al., 2010; Jiang et al., 2013) has been shown. Moreover, the observed functional dissociation between insular subregions resonates with previous accounts of insular subspecialization for different features of aversive stimuli (Baliki et al., 2009; Deen et al., 2011; Nomi et al., 2017).

Signed versus unsigned prediction errors

Prediction errors could either be represented as absolute deviations (i.e., absolute deviation of outcome from expectation) or signed deviations (i.e., more or less pain than expected). Unfortunately, due to the design of our experiment, distinguishing between unsigned and signed prediction errors cannot be dissociated in our model. Our data indicate the representation of unsigned prediction error in the anterior insula. Similar observations of unsigned prediction errors have been made in the anterior insula, amygdala, and visual and auditory cortices (den Ouden et al., 2009, 2010; Todorovic et al., 2011; Todorovic and de Lange, 2012; Boll et al., 2013; Geuter et al., 2017). Unlike the signed prediction error, which can be used to guide approach and avoidance behavior, this type of error (i.e., absolute prediction error) determines the captured attention, which is related to stimulus saliency rather than its actual value (Kahnt et al., 2014).

In conclusion, the vAI exhibits all features necessary for the implementation of a predictive coding type integration of pain. Importantly, the functional organization revealed by our model comparison is remarkably similar to the structurally defined tripartite organization of the insular cortex (Fischl et al., 2004; Destrieux et al., 2010). Our results corroborate the previously observed gradient from stimulus intensity coding to perceptual coding from posterior to anterior insula (Geuter et al., 2017). This might also reflect a gradient from nociception to pain, where somatosensory information from posterior somatosensory regions is more and more integrated with cognitive factors such as expectation, which finally allows the estimation of prediction errors in the anterior insula.

Footnotes

  • This work was supported by European Research Council Advanced Grant ERC-2010-AdG_20100407 and Deutsche Forschungsgemeinschaft Grant SFB 936 Project A06. We thank the University of Minnesota Center for Magnetic Resonance Research for providing the image reconstruction algorithm for the simultaneous multislice acquisitions; Tor D. Wager for sharing the NPS pattern; and Jürgen Finsterbusch, Katrin Bergholz, Kathrin Wendt, Waldemar Schwarz, and Matthias Kerkemeyer for valuable assistance with data collection. S.F. thanks Bjoern Horing for help on the NPS analysis, discussions, and comments on the manuscript draft.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Dr. Sepideh Fazeli, Department of Systems Neuroscience, University Medical Center Hamburg-Eppendorf, Martinistrasse 52, Hamburg, Germany 20246. sfazeli{at}uke.de

References

  1. ↵
    1. Ashburner J
    (2007) A fast diffeomorphic image registration algorithm. Neuroimage 38:95–113. doi:10.1016/j.neuroimage.2007.07.007 pmid:17761438
    OpenUrlCrossRefPubMed
  2. ↵
    1. Atlas LY,
    2. Wager TD
    (2012) How expectations shape pain. Neurosci Lett 520:140–148. doi:10.1016/j.neulet.2012.03.039 pmid:22465136
    OpenUrlCrossRefPubMed
  3. ↵
    1. Atlas LY,
    2. Bolger N,
    3. Lindquist MA,
    4. Wager TD
    (2010) Brain mediators of predictive cue effects on perceived pain. J Neurosci 30:12964–12977. doi:10.1523/JNEUROSCI.0057-10.2010 pmid:20881115
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Baliki MN,
    2. Geha PY,
    3. Apkarian AV
    (2009) Parsing pain perception between nociceptive representation and magnitude estimation. J Neurophysiol 101:875–887. doi:10.1152/jn.91100.2008 pmid:19073802
    OpenUrlCrossRefPubMed
  5. ↵
    1. Barrett LF,
    2. Simmons WK
    (2015) Interoceptive predictions in the brain. Nat Rev Neurosci 16:419–429. doi:10.1038/nrn3950 pmid:26016744
    OpenUrlCrossRefPubMed
  6. ↵
    1. Benedek M,
    2. Kaernbach C
    (2010) Decomposition of skin conductance data by means of nonnegative deconvolution. Psychophysiology 47:647–658. doi:10.1111/j.1469-8986.2009.00972.x pmid:20230512
    OpenUrlCrossRefPubMed
  7. ↵
    1. Boll S,
    2. Gamer M,
    3. Gluth S,
    4. Finsterbusch J,
    5. Büchel C
    (2013) Separate amygdala subregions signal surprise and predictiveness during associative fear learning in humans. Eur J Neurosci 37:758–767. doi:10.1111/ejn.12094 pmid:23278978
    OpenUrlCrossRefPubMed
  8. ↵
    1. Büchel C,
    2. Geuter S,
    3. Sprenger C,
    4. Eippert F
    (2014) Placebo analgesia: a predictive coding perspective. Neuron 81:1223–1239. doi:10.1016/j.neuron.2014.02.042 pmid:24656247
    OpenUrlCrossRefPubMed
  9. ↵
    1. Cerliani L,
    2. Thomas RM,
    3. Jbabdi S,
    4. Siero JC,
    5. Nanetti L,
    6. Crippa A,
    7. Gazzola V,
    8. D'Arceuil H,
    9. Keysers C
    (2012) Probabilistic tractography recovers a rostrocaudal trajectory of connectivity variability in the human insular cortex. Hum Brain Mapp 33:2005–2034. doi:10.1002/hbm.21338 pmid:21761507
    OpenUrlCrossRefPubMed
  10. ↵
    1. Chanes L,
    2. Barrett LF
    (2016) Redefining the role of limbic areas in cortical processing. Trends Cogn Sci 20:96–106. doi:10.1016/j.tics.2015.11.005 pmid:26704857
    OpenUrlCrossRefPubMed
  11. ↵
    1. Cloutman LL,
    2. Binney RJ,
    3. Drakesmith M,
    4. Parker GJ,
    5. Lambon Ralph MA
    (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. doi:10.1016/j.neuroimage.2011.11.016 pmid:22100771
    OpenUrlCrossRefPubMed
  12. ↵
    1. Craig AD
    (2002) How do you feel? Interoception: the sense of the physiological condition of the body. Nat Rev Neurosci 3:655–666. doi:10.1038/nrn894 pmid:12154366
    OpenUrlCrossRefPubMed
  13. ↵
    1. Davis KD,
    2. Kwan CL,
    3. Crawley AP,
    4. Mikulis DJ
    (1998) Functional MRI study of thalamic and cortical activations evoked by cutaneous heat, cold, and tactile stimuli. J Neurophysiol 80:1533–1546. doi:10.1152/jn.1998.80.3.1533 pmid:9744957
    OpenUrlCrossRefPubMed
  14. ↵
    1. Deen B,
    2. Pitskel NB,
    3. Pelphrey KA
    (2011) Three systems of insular functional connectivity identified with cluster analysis. Cereb Cortex 21:1498–1506. doi:10.1093/cercor/bhq186 pmid:21097516
    OpenUrlCrossRefPubMed
  15. ↵
    1. den Ouden HE,
    2. Friston KJ,
    3. Daw ND,
    4. McIntosh AR,
    5. Stephan KE
    (2009) A dual role for prediction error in associative learning. Cereb Cortex 19:1175–1185. doi:10.1093/cercor/bhn161 pmid:18820290
    OpenUrlCrossRefPubMed
  16. ↵
    1. den Ouden HE,
    2. Daunizeau J,
    3. Roiser J,
    4. Friston KJ,
    5. Stephan KE
    (2010) Striatal prediction error modulates cortical coupling. J Neurosci 30:3210–3219. doi:10.1523/JNEUROSCI.4458-09.2010 pmid:20203180
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Dennis EL,
    2. Jahanshad N,
    3. McMahon KL,
    4. de Zubicaray GI,
    5. Martin NG,
    6. Hickie IB,
    7. Toga AW,
    8. Wright MJ,
    9. Thompson PM
    (2014) Development of insula connectivity between ages 12 and 30 revealed by high angular resolution diffusion imaging. Hum Brain Mapp 35:1790–1800. doi:10.1002/hbm.22292 pmid:23836455
    OpenUrlCrossRefPubMed
  18. ↵
    1. Desikan RS,
    2. Ségonne F,
    3. Fischl B,
    4. Quinn BT,
    5. Dickerson BC,
    6. Blacker D,
    7. Buckner RL,
    8. Dale AM,
    9. Maguire RP,
    10. Hyman BT,
    11. Albert MS,
    12. Killiany RJ
    (2006) An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31:968–980. doi:10.1016/j.neuroimage.2006.01.021 pmid:16530430
    OpenUrlCrossRefPubMed
  19. ↵
    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. doi:10.1016/j.neuroimage.2010.06.010 pmid:20547229
    OpenUrlCrossRefPubMed
  20. ↵
    1. Dosenbach NU,
    2. Fair DA,
    3. Miezin FM,
    4. Cohen AL,
    5. Wenger KK,
    6. Dosenbach RA,
    7. Fox MD,
    8. Snyder AZ,
    9. Vincent JL,
    10. Raichle ME,
    11. Schlaggar BL,
    12. Petersen SE
    (2007) Distinct brain networks for adaptive and stable task control in humans. Proc Natl Acad Sci U S A 104:11073–11078. doi:10.1073/pnas.0704320104 pmid:17576922
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Downar J,
    2. Crawley AP,
    3. Mikulis DJ,
    4. Davis KD
    (2002) A cortical network sensitive to stimulus salience in a neutral behavioral context across multiple sensory modalities. J Neurophysiol 87:615–620. doi:10.1152/jn.00636.2001 pmid:11784775
    OpenUrlCrossRefPubMed
  22. ↵
    1. Dum RP,
    2. Levinthal DJ,
    3. Strick PL
    (2009) The spinothalamic system targets motor and sensory areas in the cerebral cortex of monkeys. J Neurosci 29:14223–14235. doi:10.1523/JNEUROSCI.3398-09.2009 pmid:19906970
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. DuPont S,
    2. Bouilleret V,
    3. Hasboun D,
    4. Semah F,
    5. Baulac M
    (2003) Functional anatomy of the insula: new insights from imaging. Surg Radiol Anat 25:113–119. doi:10.1007/s00276-003-0103-4 pmid:12819943
    OpenUrlCrossRefPubMed
  24. ↵
    1. Dura-Bernal S,
    2. Wennekers T,
    3. Denham SL
    (2011) The role of feedback in a hierarchical model of object perception. Adv Exp Med Biol 718:165–179. doi:10.1007/978-1-4614-0164-3_14 pmid:21744218
    OpenUrlCrossRefPubMed
  25. ↵
    1. Egner T,
    2. Monti JM,
    3. Summerfield C
    (2010) Expectation and surprise determine neural population responses in the ventral visual stream. J Neurosci 30:16601–16608. doi:10.1523/JNEUROSCI.2770-10.2010 pmid:21147999
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Feinberg DA,
    2. Moeller S,
    3. Smith SM,
    4. Auerbach E,
    5. Ramanna S,
    6. Gunther M,
    7. Glasser MF,
    8. Miller KL,
    9. Ugurbil K,
    10. Yacoub E
    (2010) Multiplexed echo planar imaging for sub-second whole brain FMRI and fast diffusion imaging. PLoS One 5:e15710. doi:10.1371/journal.pone.0015710 pmid:21187930
    OpenUrlCrossRefPubMed
  27. ↵
    1. Fischl B,
    2. van der Kouwe A,
    3. Destrieux C,
    4. Halgren E,
    5. Ségonne F,
    6. Salat DH,
    7. Busa E,
    8. Seidman LJ,
    9. Goldstein J,
    10. Kennedy D,
    11. Caviness V,
    12. Makris N,
    13. Rosen B,
    14. Dale AM
    (2004) Automatically parcellating the human cerebral cortex. Cereb Cortex 14:11–22. doi:10.1093/cercor/bhg087 pmid:14654453
    OpenUrlCrossRefPubMed
  28. ↵
    1. Friston K
    (2005) A theory of cortical responses. Philos Trans R Soc Lond B Biol Sci 360:815–836. doi:10.1098/rstb.2005.1622 pmid:15937014
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Friston K
    (2008) Hierarchical models in the brain. PLoS Comput Biol 4:e1000211. doi:10.1371/journal.pcbi.1000211 pmid:18989391
    OpenUrlCrossRefPubMed
  30. ↵
    1. Geuter S,
    2. Boll S,
    3. Eippert F,
    4. Büchel C
    (2017) Functional dissociation of stimulus intensity encoding and predictive coding of pain in the insula. eLife 6:e24770. doi:10.7554/eLife.24770 pmid:28524817
    OpenUrlCrossRefPubMed
  31. ↵
    1. Grahl A,
    2. Onat S,
    3. Büchel C
    (2018) The periaqueductal gray and Bayesian integration in placebo analgesia. eLife 7:e32930. doi:10.7554/eLife.32930 pmid:29555019
    OpenUrlCrossRefPubMed
  32. ↵
    1. Griswold MA,
    2. Jakob PM,
    3. Heidemann RM,
    4. Nittka M,
    5. Jellus V,
    6. Wang J,
    7. Kiefer B,
    8. Haase A
    (2002) Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn Reson Med 47:1202–1210. doi:10.1002/mrm.10171 pmid:12111967
    OpenUrlCrossRefPubMed
  33. ↵
    1. Häuser W,
    2. Hansen E,
    3. Enck P
    (2012) Nocebo phenomena in medicine: their relevance in everyday clinical practice. Dtsch Arztebl Int 109:459–465. doi:10.3238/arztebl.2012.0459 pmid:22833756
    OpenUrlCrossRefPubMed
  34. ↵
    1. Jakab A,
    2. Molnár PP,
    3. Bogner P,
    4. Béres M,
    5. Berényi EL
    (2012) Connectivity-based parcellation reveals interhemispheric differences in the insula. Brain Topogr 25:264–271. doi:10.1007/s10548-011-0205-y pmid:22002490
    OpenUrlCrossRefPubMed
  35. ↵
    1. Jiang J,
    2. Summerfield C,
    3. Egner T
    (2013) Attention sharpens the distinction between expected and unexpected percepts in the visual brain. J Neurosci 33:18438–18447. doi:10.1523/JNEUROSCI.3308-13.2013 pmid:24259568
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Kahnt T,
    2. Park SQ,
    3. Haynes JD,
    4. Tobler PN
    (2014) Disentangling neural representations of value and salience in the human brain. Proc Natl Acad Sci U S A 111:5000–5005. doi:10.1073/pnas.1320189111 pmid:24639493
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Keltner JR,
    2. Furst A,
    3. Fan C,
    4. Redfern R,
    5. Inglis B,
    6. Fields HL
    (2006) Isolating the modulatory effect of expectation on pain transmission: a functional magnetic resonance imaging study. J Neurosci 26:4437–4443. doi:10.1523/JNEUROSCI.4463-05.2006 pmid:16624963
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Knill DC,
    2. Pouget A
    (2004) The Bayesian brain: the role of uncertainty in neural coding and computation. Trends Neurosci 27:712–719. doi:10.1016/j.tins.2004.10.007 pmid:15541511
    OpenUrlCrossRefPubMed
  39. ↵
    1. Koyama T,
    2. McHaffie JG,
    3. Laurienti PJ,
    4. Coghill RC
    (2005) The subjective experience of pain: where expectations become reality. Proc Natl Acad Sci U S A 102:12950–12955. doi:10.1073/pnas.0408576102 pmid:16150703
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Krishnan A,
    2. Woo CW,
    3. Chang LJ,
    4. Ruzic L,
    5. Gu X,
    6. López-Solà M,
    7. Jackson PL,
    8. Pujol J,
    9. Fan J,
    10. Wager TD
    (2016) Somatic and vicarious pain are represented by dissociable multivariate brain patterns. eLife 5:e15166. doi:10.7554/eLife.15166 pmid:27296895
    OpenUrlCrossRefPubMed
  41. ↵
    1. Lang PJ,
    2. Bradley MM,
    3. Cuthbert BN
    (2005) International affective picture system (IAPS): affective ratings of pictures and instruction manual. Gainesville, FL: National Institute of Mental Health, Center for the Study of Emotion and Attention.
  42. ↵
    1. Langner R,
    2. Kellermann T,
    3. Boers F,
    4. Sturm W,
    5. Willmes K,
    6. Eickhoff SB
    (2011) Modality-specific perceptual expectations selectively modulate baseline activity in auditory, somatosensory, and visual cortices. Cereb Cortex 21:2850–2862. doi:10.1093/cercor/bhr083 pmid:21527785
    OpenUrlCrossRefPubMed
  43. ↵
    1. Lewandowsky S,
    2. Farrell S
    (2010) Computational Modeling in Cognition: Principles and Practice. California, CA: SAGE Publications, Inc.
  44. ↵
    1. Moeller S,
    2. Yacoub E,
    3. Olman CA,
    4. Auerbach E,
    5. Strupp J,
    6. Harel N,
    7. Uğurbil K
    (2010) Multiband multislice GE-EPI at 7 tesla, with 16-fold acceleration using partial parallel imaging with application to high spatial and temporal whole-brain fMRI. Magn Reson Med 63:1144–1153. doi:10.1002/mrm.22361 pmid:20432285
    OpenUrlCrossRefPubMed
  45. ↵
    1. Mutschler I,
    2. Wieckhorst B,
    3. Kowalevski S,
    4. Derix J,
    5. Wentlandt J,
    6. Schulze-Bonhage A,
    7. Ball T
    (2009) Functional organization of the human anterior insular cortex. Neurosci Lett 457:66–70. doi:10.1016/j.neulet.2009.03.101 pmid:19429164
    OpenUrlCrossRefPubMed
  46. ↵
    1. Nieuwenhuys R
    (2012) The insular cortex: a review. Prog Brain Res 195:123–163. doi:10.1016/B978-0-444-53860-4.00007-6 pmid:22230626
    OpenUrlCrossRefPubMed
  47. ↵
    1. Nomi JS,
    2. Schettini E,
    3. Broce I,
    4. Dick AS,
    5. Uddin LQ
    (2017) Structural connections of functionally defined human insular subdivisions. Cereb Cortex 2017:1–12. doi:10.1093/cercor/bhx211 pmid:28968768
    OpenUrlCrossRefPubMed
  48. ↵
    1. Okada K,
    2. Matchin W,
    3. Hickok G
    (2018) Neural evidence for predictive coding in auditory cortex during speech production. Psychon Bull Rev 25:423–430. doi:10.3758/s13423-017-1284-x pmid:28397076
    OpenUrlCrossRefPubMed
  49. ↵
    1. Ostrowsky K,
    2. Magnin M,
    3. Ryvlin P,
    4. Isnard J,
    5. Guenot M,
    6. Mauguière F
    (2002) Representation of pain and somatic sensation in the human insula: a study of responses to direct electrical cortical stimulation. Cereb Cortex 12:376–385. doi:10.1093/cercor/12.4.376 pmid:11884353
    OpenUrlCrossRefPubMed
  50. ↵
    1. Ploner M,
    2. Lee MC,
    3. Wiech K,
    4. Bingel U,
    5. Tracey I
    (2010) Prestimulus functional connectivity determines pain perception in humans. Proc Natl Acad Sci U S A 107:355–360. doi:10.1073/pnas.0906186106 pmid:19948949
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Rao RP,
    2. Ballard DH
    (1999) Predictive coding in the visual cortex: a functional interpretation of some extra-classical receptive-field effects. Nat Neurosci 2:79–87. doi:10.1038/4580 pmid:10195184
    OpenUrlCrossRefPubMed
  52. ↵
    1. Segerdahl AR,
    2. Mezue M,
    3. Okell TW,
    4. Farrar JT,
    5. Tracey I
    (2015) The dorsal posterior insula subserves a fundamental role in human pain. Nat Neurosci 18:499–500. doi:10.1038/nn.3969 pmid:25751532
    OpenUrlCrossRefPubMed
  53. ↵
    1. Setsompop K,
    2. Gagoski BA,
    3. Polimeni JR,
    4. Witzel T,
    5. Wedeen VJ,
    6. Wald LL
    (2012) Blipped-controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g-factor penalty. Magn Reson Med 67:1210–1224. doi:10.1002/mrm.23097 pmid:21858868
    OpenUrlCrossRefPubMed
  54. ↵
    1. Seymour B,
    2. O'Doherty JP,
    3. Dayan P,
    4. Koltzenburg M,
    5. Jones AK,
    6. Dolan RJ,
    7. Friston KJ,
    8. Frackowiak RS
    (2004) Temporal difference models describe higher-order learning in humans. Nature 429:664–667. doi:10.1038/nature02581 pmid:15190354
    OpenUrlCrossRefPubMed
  55. ↵
    1. Sterzer P,
    2. Frith C,
    3. Petrovic P
    (2008) Believing is seeing: expectations alter visual awareness. Curr Biol 18:R697–R698. doi:10.1016/j.cub.2008.06.021 pmid:18727901
    OpenUrlCrossRefPubMed
  56. ↵
    1. Summerfield C,
    2. Trittschuh EH,
    3. Monti JM,
    4. Mesulam MM,
    5. Egner T
    (2008) Neural repetition suppression reflects fulfilled perceptual expectations. Nat Neurosci 11:1004–1006. doi:10.1038/nn.2163 pmid:19160497
    OpenUrlCrossRefPubMed
  57. ↵
    1. Taylor KS,
    2. Seminowicz DA,
    3. Davis KD
    (2009) Two systems of resting state connectivity between the insula and cingulate cortex. Hum Brain Mapp 30:2731–2745. doi:10.1002/hbm.20705 pmid:19072897
    OpenUrlCrossRefPubMed
  58. ↵
    1. Todorovic A,
    2. de Lange FP
    (2012) Repetition suppression and expectation suppression are dissociable in time in early auditory evoked fields. J Neurosci 32:13389–13395. doi:10.1523/JNEUROSCI.2227-12.2012 pmid:23015429
    OpenUrlAbstract/FREE Full Text
  59. ↵
    1. Todorovic A,
    2. van Ede F,
    3. Maris E,
    4. de Lange FP
    (2011) Prior expectation mediates neural adaptation to repeated sounds in the auditory cortex: an MEG study. J Neurosci 31:9118–9123. doi:10.1523/JNEUROSCI.1425-11.2011 pmid:21697363
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Tracey I
    (2010) Getting the pain you expect: mechanisms of placebo, nocebo and reappraisal effects in humans. Nat Med 16:1277–1283. doi:10.1038/nm.2229 pmid:20948533
    OpenUrlCrossRefPubMed
  61. ↵
    1. Uddin LQ,
    2. Kinnison J,
    3. Pessoa L,
    4. Anderson ML
    (2014) Beyond the tripartite cognition-emotion-interoception model of the human insular cortex. J Cogn Neurosci 26:16–27. doi:10.1162/jocn_a_00462 pmid:23937691
    OpenUrlCrossRefPubMed
  62. ↵
    1. Wager TD,
    2. Atlas LY,
    3. Leotti LA,
    4. Rilling JK
    (2011) Predicting individual differences in placebo analgesia: contributions of brain activity during anticipation and pain experience. J Neurosci 31:439–452. doi:10.1523/JNEUROSCI.3420-10.2011 pmid:21228154
    OpenUrlAbstract/FREE Full Text
  63. ↵
    1. Wager TD,
    2. Atlas LY,
    3. Lindquist MA,
    4. Roy M,
    5. Woo CW,
    6. Kross E
    (2013) An fMRI-based neurologic signature of physical pain. N Engl J Med 368:1388–1397. doi:10.1056/NEJMoa1204471 pmid:23574118
    OpenUrlCrossRefPubMed
  64. ↵
    1. Wiech K
    (2016) Deconstructing the sensation of pain: the influence of cognitive processes on pain perception. Science 354:584–587. doi:10.1126/science.aaf8934 pmid:27811269
    OpenUrlAbstract/FREE Full Text
  65. ↵
    1. Wiech K,
    2. Ploner M,
    3. Tracey I
    (2008) Neurocognitive aspects of pain perception. Trends Cogn Sci 12:306–313. doi:10.1016/j.tics.2008.05.005 pmid:18606561
    OpenUrlCrossRefPubMed
  66. ↵
    1. Wiech K,
    2. Lin CS,
    3. Brodersen KH,
    4. Bingel U,
    5. Ploner M,
    6. Tracey I
    (2010) Anterior insula integrates information about salience into perceptual decisions about pain. J Neurosci 30:16324–16331. doi:10.1523/JNEUROSCI.2087-10.2010 pmid:21123578
    OpenUrlAbstract/FREE Full Text
  67. ↵
    1. Wiech K,
    2. Jbabdi S,
    3. Lin CS,
    4. Andersson J,
    5. Tracey I
    (2014) Differential structural and resting state connectivity between insular subdivisions and other pain-related brain regions. Pain 155:2047–2055. doi:10.1016/j.pain.2014.07.009 pmid:25047781
    OpenUrlCrossRefPubMed
  68. ↵
    1. Xu J,
    2. Moeller S,
    3. Auerbach EJ,
    4. Strupp J,
    5. Smith SM,
    6. Feinberg DA,
    7. Yacoub E,
    8. Uğurbil K
    (2013) Evaluation of slice accelerations using multiband echo planar imaging at 3 T. Neuroimage 83:991–1001. doi:10.1016/j.neuroimage.2013.07.055 pmid:23899722
    OpenUrlCrossRefPubMed
  69. ↵
    1. Ylinen S,
    2. Huuskonen M,
    3. Mikkola K,
    4. Saure E,
    5. Sinkkonen T,
    6. Paavilainen P
    (2016) Predictive coding of phonological rules in auditory cortex: a mismatch negativity study. Brain Lang 162:72–80. doi:10.1016/j.bandl.2016.08.007 pmid:27588355
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 38 (29)
Journal of Neuroscience
Vol. 38, Issue 29
18 Jul 2018
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Pain-Related Expectation and Prediction Error Signals in the Anterior Insula Are Not Related to Aversiveness
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Pain-Related Expectation and Prediction Error Signals in the Anterior Insula Are Not Related to Aversiveness
Sepideh Fazeli, Christian Büchel
Journal of Neuroscience 18 July 2018, 38 (29) 6461-6474; DOI: 10.1523/JNEUROSCI.0671-18.2018

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Pain-Related Expectation and Prediction Error Signals in the Anterior Insula Are Not Related to Aversiveness
Sepideh Fazeli, Christian Büchel
Journal of Neuroscience 18 July 2018, 38 (29) 6461-6474; DOI: 10.1523/JNEUROSCI.0671-18.2018
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • expectation
  • human fMRI
  • pain
  • predictive coding
  • somatosensory perception

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Prostaglandin E2 induces long-lasting inhibition of noradrenergic neurons in the locus coeruleus and moderates the behavioral response to stressors
  • Detection of spatially-localized sounds is robust to saccades and concurrent eye movement-related eardrum oscillations (EMREOs)
  • Structural and functional development of inhibitory connections from the medial nucleus of the trapezoid body to the superior paraolivary nucleus
Show more Research Articles

Behavioral/Cognitive

  • Prostaglandin E2 induces long-lasting inhibition of noradrenergic neurons in the locus coeruleus and moderates the behavioral response to stressors
  • Detection of spatially-localized sounds is robust to saccades and concurrent eye movement-related eardrum oscillations (EMREOs)
  • Rewarding capacity of optogenetically activating a giant GABAergic central-brain interneuron in larval Drosophila
Show more Behavioral/Cognitive
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.