A Multisensory fMRI Investigation of Nociceptive-Preferential Cortical Regions and Responses

The existence of nociceptive-specific brain regions has been a controversial issue for decades. Multisensory fMRI studies, which examine fMRI activities in response to various types of sensory stimulation, could help identify nociceptive-specific brain regions, but previous studies are limited by sample size and they did not differentiate nociceptive-specific regions and nociceptive-preferential regions, which have significantly larger responses to nociceptive input. In this study, we conducted a multisensory fMRI experiment on 80 healthy participants, with the aim to determine whether there are certain brain regions that specifically or preferentially respond to nociceptive stimulation. By comparing the evoked fMRI responses across four sensory modalities, we found a series of brain regions specifically or preferentially involved in nociceptive sensory input. Particularly, we found different parts of some cortical regions, such as insula and cingulate gyrus, play different functional roles in the processing of nociceptive stimulation. Hence, this multisensory study improves our understanding of the functional integrations and segregations of the nociceptive-related regions.

Introduction

Pain is an unpleasant sensory and emotional experience associated with actual or potential tissue damage or described in terms of such damage (IASP, International Association for the Study of Pain) (Loeser and Treede, 2008). By using neuroimaging methods, such as functional magnetic resonance imaging (fMRI), electroencephalography (EEG), and positron emission tomography (PET), we can detect brain activities and cortical regions that are induced by pain perception (Tölle et al., 1999; Iannetti et al., 2005; Moisset and Bouhassira, 2007; Wager et al., 2013; Jackson et al., 2020). More specifically, it has been well-documented that nociceptive pain (Mischkowski et al., 2018), which is caused by noxious stimuli in intact tissue, can trigger brain activities in a network of cortical regions, including the thalamus, insula, the primary and secondary somatosensory cortices, and the cingulate cortex, which are often called “pain matrix” (Ingvar, 1999; Peyron et al., 2000, 2002; Apkarian et al., 2005).

However, whether the pain matrix is specific to pain (i.e., activated by pain stimulation only) has long been a controversial issue. Some pain researchers believe that at least part of the pain matrix is specific to nociceptive pain (Stern et al., 2006; Tracey and Mantyh, 2007; Boly et al., 2008; Borsook et al., 2010), but more evidence have shown that almost all cortical regions in the pain matrix can be activated by other sensory modalities (Downar et al., 2000; Macaluso and Driver, 2005; Lui et al., 2008; Mouraux and Iannetti, 2009; Mouraux et al., 2011). Some multisensory studies, which used a random sequence of different types of sensory stimulation to examine activated cortical regions, found that non-nociceptive somatosensory, auditory, and visual stimulation can activate cortical regions largely similar to those regions activated nociceptive somatosensory stimulation (Downar et al., 2000; Mouraux and Iannetti, 2009; Mouraux et al., 2011). Based on this phenomenon, many researchers put forward another possible point of view: the so-called “pain matrix” may be related to the significant input of sensory stimulation and the ability to respond to it, regardless of the type of sensory stimulation input (Iannetti and Mouraux, 2010; Legrain et al., 2011; Mouraux et al., 2011; Iannetti et al., 2013). In another word, the pain matrix is not specific to pain.

However, these multisensory studies were based on a relatively small sample size (the number of participants was smaller than 20), and the statistical power of the results was weak. Another important limitation of existing multisensory studies concerning the nociceptive pain is that, they were focused on the identification of nociceptive-specific regions, but overlooked the possibility that some cortical regions may be preferential to nociceptive pain, even they are not specific to nociceptive pain. Nociceptive-preferential regions are cortical regions that can be evoked by stimulation of various types of sensory modalities but are more preferentially activated by nociceptive stimulation. In another word, nociceptive-preferential regions may have significantly greater response to nociceptive pain rather than to other sensory modalities. However, the existence of nociceptive-preferential regions has been seldom investigated by a multisensory study. It is important to identify and differentiate nociceptive-specific and nociceptive-preferential cortical regions to improve our understanding of the functional integrations and segregations of these regions.

In this study, we hypothesized that at least some cortical regions are preferentially involved in the processing of nociceptive input and validated our hypothesis based on fMRI data recorded from a relatively large number of participants (80 healthy participants) in a multisensory experiment. We applied a random sequence of stimuli with four different sensory modalities, including nociceptive somatosensory, non-nociceptive somatosensory, auditory, and visual, on these participants, and recorded the fMRI responses. We used the conjunction analysis on the evoked fMRI activities to identify cortical regions that are specifically activated by one type of sensory modality (i.e., modality-specific brain regions). We further compared evoked fMRI magnitudes among four sensory modalities to identify cortical regions that are preferentially activated by one specific sensory modality (i.e., modality-preferential brain regions). These identified nociceptive-specific and -preferential regions were further compared to check their overlap and difference.

Materials and Methods

Participants

Eighty healthy right-handed participants aged 21.0 ± 3.1 years, with 36 men and 44 women, were recruited. All participants were healthy and had no neurological or cardiovascular diseases or neurological disorders or acute/chronic pain. The study was approved by the local Ethics Committee and all participants gave written informed consent.

Experiment

Each participant received four types of sensory modality: nociceptive somatosensory, non-nociceptive somatosensory, auditory, and visual, presented in a random order with similar attentional level. Nociceptive somatosensory stimuli were infrared neodymium yttrium aluminum perovskite (Nd: YAP) laser pulses delivered to the left-hand dorsum. The laser pulse duration was approximately 4 ms, and the spot size was 7 mm in diameter. Non-nociceptive somatosensory stimuli were bipolar square wave electrical pulses stimulation (100 Hz frequency, 1,000 μs pulse width, single pulse, forward current, and 1 ms duration). The stimuli were delivered through a pair of skin electrodes (1 cm inter-electrode distance) placed at the inside of the left arm, over the median nerve. Auditory stimuli adopted binaural monophonic stimulation (800 Hz auditory tones lasting 50 ms, 5 ms rise and fall times) and delivered through pneumatic earphones. Visual stimuli consisted of a bright white disk displayed on the projection screen for 50 ms. For each sensory modality, two stimulus intensity levels were used: 3 and 3.5 J for nociceptive somatosensory, 2 and 4 mA for non-nociceptive somatosensory, 76 dB SPL and 88 dB SPL for auditory stimuli, two different grayscale images for visual stimuli.

The experiment consisted of two runs. Each run consisted of 4 types of modality stimuli, and each was delivered 20 times in a pseudo-random order, 10 times for each intensity level. A white fixation cross was presented for 600 s at the center of the black screen, and the stimulus was delivered in the meantime. Ten seconds after one stimulus, the participant was required to rate the intensity of the stimulus received using a visual analog scale (VAS) ranging between 0 (not perceived) and 10 (maximum intensity) within 5 s. The interval between each stimulus is 22–23 s, and each run lasted for about 15 min.

fMRI data were acquired using a standard gradient echo-planar imaging sequence by a GE 3.0T MRI scanner with the following parameters: 43 axial slices, thickness/gap = 3/0 mm, time of echo (TE) = 30 ms, time of repetition (TR) = 2,000 ms, acquisition matrix = 64 × 64, field of view (FOV) = 192 × 192 mm², flip angle (FA) = 90°, acquisition time = 30:12 min. A total of two sessions with 454 functional volumes were collected in the multisensory stimulation task. At the end of the task experiment, a high resolution T1-weighted structural image was collected for spatial registration with the following parameters: 176 sagittal slices, TE = 2.992 ms, TR = 6.896 ms, inversion time T1 = 450 ms, 1 mm slice thickness with no gap, 1 × 1 mm in-plane resolution, acquisition matrix = 256 × 256, acquisition time = 4:36 min.

fMRI Data Preprocessing

fMRI data were preprocessed using Statistical Parametric Mapping (SPM12)¹. The first three volumes were removed for signal equilibration, and 451 volumes were remained. Functional images were slice-time corrected, motion-corrected (realigned), coregistered to the structural images, spatially normalized to the Montreal Neurological Institute (MNI) template, resampled to 3 × 3 × 3 mm voxel size, spatial smoothed using a full-width at half-maximum (FWHM) Gaussian kernel of 8 mm, and then high-pass temporally filtered (1/128 Hz cutoff).

General Linear Model

The task activated fMRI patterns were analyzed using the General Linear Model (GLM) in SPM. GLM was performed in two levels. The first-level analysis was performed- for each participant with regressors modeling the occurrence of each condition (nociceptive somatosensory, non-nociceptive somatosensory, auditory, visual, and the rating period) and six additional regressors of head motion parameters. Four contrast analyses corresponding to the four stimulus types were performed to assess the BOLD responses of each participant and to produce a statistical parametric map for each participant. The first-level statistical parametric maps of all participants were then entered into a second-level one-sample t-test analysis for each stimulation. The second-level statistical results were further performed by the Benjamini-Hochberg procedure to decrease the False Discovery Rate (FDR) and the corrected significance level was set to 0.0001. We also used other significance thresholds to examine the robustness of the results (for example, results obtained with P < 0.05 with FDR correction can be found in Supplementary Figure 1).

Identification of Modality-Specific Regions (Conjunction Analysis)

To identify brain regions specific to each type of stimulation (i.e., activities only elicited by one sensory modality, but not by other three), we performed a conjunction analysis based on the Minimum Statistic compared to the Conjunction Null (MS/CN) method (Nichols et al., 2005). We calculated the intersection of four modalities’ statistical maps, which were obtained by GLM to identify multimodal regions, i.e., those brain regions co-activated by all four modalities. And the brain regions specific to one specific sensory modality can be identified as those regions that were activated by only this sensory modality but not by others. All the identified voxels were defined using the automatic anatomical labeling (AAL) atlas (Collins et al., 1998), and the threshold of cluster size was set to 20 (i.e., only a cluster with more than 20 voxels was retained).

Identification of Modality-Preferential Regions (PSC Analysis)

Further, in order to identify modality-preferential regions (i.e., those regions are preferentially activated by one specific sensory modality), we extracted percentage signal change (PSC) at each voxel (Pernet, 2014). PSC features were used here because PSC can more accurately fit the fMRI data and compensate temporal shift than GLM-based coefficients (Calhoun et al., 2004). PSC was calculated as the ratio between the magnitude of the BOLD response accounting for the temporal shift and the mean of the adjusted time series. Before extracting PSC features, preprocessed fMRI data (see section “fMRI Data Preprocessing”) were further processed with the following extra steps: (1) linearly detrending the time-series data to eliminate signal drift caused by factors such as machine heat and pulse, (2) regressing out six motion parameters and the rating regressor modeling the occurrence of the rating period.

Specifically, for each sensory modality (i.e., nociceptive somatosensory, non-nociceptive somatosensory, auditory, and visual), we extracted time-series data according to the stimulus onset for each trial at each voxel, and then we calculated the PSC for each stimulus, as described previously. We then averaged the PSC cross trials for each sensory modality. Then, we obtained the average PSC values of each voxel for each sensory modality and each subject.

To identify voxels that are preferentially activated by one certain sensory modality, we compared extracted PSC features among sensory modalities using a voxel-wise one-way repeated measures analysis of variance (ANOVA), followed by post-hoc Tukey-Kramer comparisons. To address the multiple comparisons problem, we computed the Benjamini-Hochberg correction of the FDR and set the corrected significance level to 0.0001. To examine the robustness of the results, brain regions identified using P < 0.05 with FDR correction are shown in Supplementary Figure 2. For each voxel, if its average PSC of one sensory modality was significantly larger than the average PSCs of other three sensory modalities, then this voxel was considered as “preferentially” activated by this sensory modality. Above analysis was performed for each sensory modality and we finally had four modality-preferential maps for four types of sensory stimulation. All the identified voxels were defined using the AAL atlas (Iannetti et al., 2013), and the threshold of cluster size was set to 20.

Further, we examined the across-trial correlation between subjective rating and PSCs in the modality-preferential regions for each sensory modality. For each subject, the PSCs of each trial were averaged across all voxels within the modality-preferential regions. We then calculated Pearson’s correlations between PSCs and VAS ratings across trials for each subject. The obtained correlation coefficients for each subject were then normalized by Fisher’s z-transformation, and one-sample t-test was performed to determine whether there was a significant cross-subject correlation between the VAS rating and PSCs for each sensory modality.

Results

Behavioral Results

Figure 1 shows the VAS ratings of each type of sensory modality. The average VAS rating for each sensory modality was (mean ± SD): nociceptive somatosensory: 4.858 ± 1.472; non-nociceptive somatosensory, 4.715 ± 1.335; auditory: 4.550 ± 1.241; visual: 4.771 ± 1.061. There was no significant difference among the VAS ratings of four sensory modalities (F (2.655, 209.752) = 1.936, P = 0.132, one-way repeated measures ANOVA).

FIGURE 1

Figure 1. Visual analog scale (VAS) ratings of each type of sensory modality. Each dot represents the average rating of one participant to one type of sensory modality.

General Linear Model Analysis

The GLM-based activated regions are shown in Figure 2. Nociceptive somatosensory, non-nociceptive somatosensory, auditory, and visual stimuli elicited widely distributed fMRI brain activities throughout the brain regions. Besides, the considerable overlaps among the brain regions activated by different sensory modalities were noticeable.

FIGURE 2

Figure 2. Brain regions activated by each type of sensory stimulation, as revealed by General Linear Model (GLM) (P < 0.0001, FDR-corrected).

Multimodal regions (i.e., those brain regions co-activated by nociceptive somatosensory, non-nociceptive somatosensory, auditory, and visual stimuli) on wide-field included the prefrontal cortex, sub-lobar (mainly including the thalamus, caudate, and putamen), temporal lobe (mainly including superior temporal gyrus, middle temporal gyrus), insula, cingulate gyrus, precuneus, inferior parietal lobule, precentral gyrus, and cuneus (Figure 3).

FIGURE 3

Figure 3. Binary multimodal regions defined by conjunction analysis.

Modality-Specific Regions Identified by Conjunction Analysis

Conjunction analysis in specific sensory modality showed that there were some brain regions specifically elicited by nociceptive somatosensory, non-nociceptive somatosensory, auditory, or visual stimuli at P < 0.0001 with FDR correction (Figure 4 and Table 1), while P < 0.05 for the specific brain regions see Supplementary Figure 1. Nociceptive somatosensory-specific regions were mainly in the ipsilateral and contralateral cingulate gyrus, the ipsilateral and contralateral precuneus, and the contralateral frontal lobe (including the supplementary motor area, dorsolateral and medial part of superior frontal gyrus). Brain responses activated by stimulation of non-nociceptive somatosensory were mainly in the contralateral postcentral and inferior parietal lobule. Auditory-specific regions were mainly in the temporal gyrus, and there was also a small number of activations in the frontal gyrus, paracentral lobule, rolandic operculum, postcentral gyrus and supramarginal gyrus. Visual-specific regions were mainly in the occipital lobe, extending into the temporal and parietal lobe, and also in the precentral and postcentral gyrus.

FIGURE 4

Figure 4. Sensory modality-specific regions defined by conjunction analysis of General Linear Model (GLM) maps. Significant results were identified at P < 0.0001 with FDR correction. ACG, anterior cingulate gyrus; ANG, angular gyrus; CAL, calcarine sulcus; MCG, middle cingulate gyrus; FFG, fusiform gyrus; IPL, inferior parietal lobule; ITG, inferior temporal gyrus; MFG, middle frontal gyrus; MOG, middle occipital gyrus; MTG, middle temporal gyrus; PCG, posterior cingulate gyrus; PCL, paracentral lobule; PCUN, precuneus; PoCG, postcentral gyrus; ROL, rolandic operculum; PreCG, precentral gyrus; SFGdor, dorsolateral part of superior frontal gyrus; SFGmed, medial part of superior frontal gyrus; SMA, supramarginal gyrus; SOG, superior occipital gyrus; SPL, superior parietal lobule; STG, superior temporal gyrus.

TABLE 1

Table 1. List of Modality-specific brain regions as revealed by General Linear Model (GLM).

Modality-Preferential Regions Identified by PSC Analysis

Figure 5 and Table 2 showed the PSC features that are preferentially modulated by sensory modality, as revealed by voxel-wise one-way repeated measures ANOVA at P < 0.0001 with FDR correction, while P < 0.05 for the preferential brain regions see Supplementary Figure 2. The nociceptive somatosensory-preferential regions include the cingulate gyrus, frontal lobe (including the supplementary motor area, the opercular part of inferior frontal gyrus, the dorsolateral part of superior frontal gyrus), insula, and rolandic operculum. The auditory-preferential regions include the ipsilateral and contralateral temporal gyrus, and visual-preferential regions were mainly in and around the occipital lobe. However, we did not find any non-nociceptive somatosensory-preferential cortical regions.

FIGURE 5

Figure 5. Sensory modality-preferential regions defined by voxel-wise one-way repeated measures ANOVA. Significant results were identified at P < 0.0001 with FDR correction. ACG, anterior cingulate gyrus; ANG, angular gyrus; MCG, middle cingulate gyrus; FFG, fusiform gyrus; IFGoperc, opercular part of inferior frontal gyrus; INS, insula; IOG, inferior occipital gyrus; IPL, inferior parietal lobule; ITG, inferior temporal gyrus; MOG, middle occipital gyrus; MTG, middle temporal gyrus; ROL, rolandic operculum; SFGdor, dorsolateral part of superior frontal gyrus; SMA, supramarginal gyrus; SOG, superior occipital gyrus; SPL, superior parietal lobule; STG, superior temporal gyrus.

TABLE 2

Table 2. List of modality-preferential brain regions as revealed by percentage signal change (PSC).

Further, we list the overlaps between modality-preferential regions and modality-specific regions in Table 3. It can be seen that PSC-based modality-preferential regions and GLM-based modality-specific regions were partially similar, though some differences existed. For example, we can find non-nociceptive-specific regions but cannot find non-nociceptive-preferential regions. According to the definitions of nociceptive-specific and -preferential regions, nociceptive-specific regions should be a subset of nociceptive- preferential regions.

TABLE 3

Table 3. Common regions of modality-specific regions and modality-preferential regions.

We tested for correlation between VAS ratings and PSCs for the modality-preferential regions. The correlation analysis revealed a significant correlation for visual (P = 0.0042) and nociceptive somatosensory (P = 1.34268 × 10^–9), and a marginally significant correlation for auditory (P = 0.0536).

Discussion

In this study, we conducted a multisensory fMRI experiment on a relatively large number of subjects (80 healthy participants), with the aim to identify nociceptive-specific regions and nociceptive-preferential brain regions. By using GLM and conjunction analysis, we identified a series of modality-specific regions, which were in response to one certain type of sensory modality only. In addition, by using PSC features and ANOVA, we identified a series of modality-preferential regions, which preferentially responded to a certain sensory modality, instead of simply following the conventional concept of sensory-specific regions.

Sensory Modalities Used in This Study

In this work, we used four types of sensory stimuli, nociceptive somatosensory, non-nociceptive somatosensory, auditory, and visual, aiming to identify nociceptive-specific or -preferential cortical regions. Previous studies (Baliki et al., 2009; Mouraux and Iannetti, 2009; Mouraux et al., 2011; Liberati et al., 2016) have repeatedly found that, nociceptive-somatosensory-evoked brain regions are not only overlapped with non-nociceptive-somatosensory-evoked regions, but also overlapped with visual- or auditory-evoked regions. Importantly, there may exist some regions that are evoked by nociceptive-somatosensory and visual (or auditory), but not by non-nociceptive-somatosensory stimuli (see Supplementary Figures 3, 4, which were obtained from our data). Therefore, we used visual or auditory stimuli to rule out more false-positive results that were not real nociceptive-specific cortical regions. Actually, a perfect multisensory study to investigate modality-specific regions should include all types of sensory modalities, which is impossible in practice. Therefore, to achieve a tradeoff, we followed other classical multisensory studies (Downar et al., 2000; Mouraux et al., 2011; Morrow et al., 2019; Sanchez et al., 2020) to include visual and auditory stimuli, which are common and easy to implement, in our experiment. In addition, because the brain response patterns of visual and auditory stimuli are well-documented, the results could be used to benchmark the data analysis method used.

Comparison Between Modality-Specific and Modality-Preferential Regions

By using the conventional GLM and conjunction analysis, we identified a set of nociceptive-specific regions, including the dorsolateral part of superior frontal gyrus, supplementary motor area, medial part of superior frontal gyrus, cingulate gyrus and precuneus. On the other hand, by based on PSC features, the identified nociceptive-preferential regions include the rolandic operculum, dorsolateral part of superior frontal gyrus, opercular part of inferior frontal gyrus, supplementary motor area, insula, anterior, and middle cingulate gyrus.

We can see from Tables 1–3 that, although the identified modality-specific and modality-preferential regions shared some common regions (as listed in Table 3), they still have remarkable difference. Theoretically, modality-specific regions should be a subset of modality-preferential regions, because one region only evoked by one type of sensory modality will have greater responses to this sensory stimulation than to other types of sensory stimulation. But our results did not agree with the above ideal and expected results. For example, we found non-nociceptive-specific regions but could not find any non-nociceptive-specific regions. It is because that the methods used to identify modality-specific and modality-preferential regions are different and it is not fair to directly compare modality-specific and modality-preferential regions identified in this work. Actually, we followed the experimental design and data analysis procedure in Mouraux et al. (2011) to check the existence of modality-specific regions. Our findings are in general consistent with the results reported in Mouraux et al. (2011). However, because this study has much more participants, we can identify some nociceptive-specific regions but Mouraux et al. (2011) did not. To further identify modality-preferential regions, we used a new data analysis pipeline (PSC analysis). As a consequence, the modality-specific and modality-preferential regions identified in the present study are not directly comparable.

Functional Roles of Nociceptive-Preferential Regions

In the following, we mainly discussed a set of important brain regions identified as nociceptive-preferential regions. For other sensory modalities, their specific or preferential regions are highly consistent with literature, so these regions are not discussed here.

Posterior Insula

We observed that in the insula, multimodal regions (i.e., voxels activated by nociceptive somatosensory, non-nociceptive somatosensory, auditory, and visual regardless of the sensory modality) and nociceptive somatosensory-specific regions included the anterior insula, while nociceptive somatosensory-preferential regions included the anterior and posterior insula. This finding is consistent with previous results supporting that the posterior insula is involved in the processing of pain perception and nociceptive sensory input (Segerdahl et al., 2015; Wu et al., 2017; Liberati et al., 2019), and this region can cause pain when it is directly electrically stimulated (Mazzola et al., 2009, 2012, 2019). Our results showed that the posterior insula was nociceptive-preferential, but not specific for nociceptive somatosensory, suggesting subregions of insula may play different roles in sensory perception.

Supplementary Motor Area and Dorsolateral Superior Frontal Gyrus

Our results indicate that the BOLD response in the supplementary motor area and dorsolateral superior frontal gyrus was nociceptive somatosensory-specific and nociceptive somatosensory-preferential. Both supplementary motor area and dorsolateral superior frontal gyrus belong to the frontal lobe. The dorsolateral superior frontal gyrus is found to be correlated with the control of complex movements (Martino et al., 2011), and the response in the supplementary motor area is often considered to play an important role in motor control (Eccles, 1982; Tanji and Shima, 1994; Tanji, 1996). Besides, imaging movements of a hand without performing, it can also activate the supplementary motor area (Lotze et al., 1999; Al-Wasity et al., 2021). Therefore, these brain regions related to motor perception were activated may be interpreted as the participants have a strong desire to inhibit the movement and take away their hands because of the unpleasant feeling of nociceptive somatosensory stimuli.

Cingulate Cortex

Nociceptive somatosensory-specific and nociceptive somatosensory-preferential regions included the cingulate gyrus, but they did not overlap. Nociceptive somatosensory-specific regions in the cingulum included the anterior and middle cingulate gyrus, while nociceptive somatosensory-preferential regions included the anterior, middle and posterior cingulate gyrus. Numerous studies repeatedly observe that the cingulate gyrus is involved in nociception (Davis et al., 1995; Dum et al., 2009; Jensen et al., 2016). The anterior cingulate cortex is involved in pain processing (Devinsky et al., 1995; Hutchison et al., 1999; Fulbright et al., 2001), the middle cingulate cortex is correlated with the attentional orienting to pain (Peyron et al., 2000), and the posterior cingulate cortex is thought to be correlated with pain memory (Bromm, 2001). These results and our results showed different roles of different parts of cingulate gyrus in the processing of sensory inputs.

Limitation and Future Work

Our results are different from the previous study (Mouraux et al., 2011), which claimed that nociceptive-specific regions do not exist. Such a difference may be due to the difference of the number of participants [14 participants in the previous study (Mouraux et al., 2011) and 80 participants in the present study] and the difference of experimental design. On the other hand, this study is still not sufficient to ascertain the existence of nociceptive-specific regions because of the following three reasons. First, the number of participants was still limited. Second, the experimental design did not cover all possible types of simulation and all possible stimulation intensity. We only used four types of stimuli (visual, auditory, non- nociceptive somatosensory, and nociceptive somatosensory) in the experiments. It is almost impossible to exhaust the full range of sensory modalities, so this work is a limited but feasible attempt. Also, we believe that the real nociceptive-specific cortical regions might be a subset of the nociceptive -specific regions identified in this study, because (1) both real and identified specific regions can respond to one sensory modality, (2) real specific regions do not respond to all other possible sensory modalities while identified specific regions do not respond to a small number of other sensory modalities included in the experiment. Third, these identified nociceptive-specific regions were sometimes reported to be respond to various types of stimulation and tasks in literature. It is always difficult to claim one brain region is specific to one type of stimulation or task, because normally a brain region has multiple functions. Anyhow, we believe this study is still an important step to advance our understanding of human brain’s processing of different types of sensory modalities.

Conclusion

In this study, we conducted a multisensory fMRI experiment on 80 healthy participants, with the aim to determine whether there are certain brain regions that specifically or preferentially respond to nociceptive stimulation. We found a series of brain responses involved in nociceptive sensory input. Specifically, we identified nociceptive-specific regions including the dorsolateral part of superior frontal gyrus, supplementary motor area, medial part of superior frontal gyrus, cingulate gyrus and precuneus, and nociceptive-preferential regions including the rolandic operculum, dorsolateral part of superior frontal gyrus, opercular part of inferior frontal gyrus, supplementary motor area, insula, anterior and middle cingulate gyrus. The nociceptive somatosensory specific and preferential regions may play a vital role in the processing of nociceptive sensory input.

Data Availability Statement

The datasets presented in this article are not readily available due to the nature of this research, participants of this study did not agree for their data to be shared publicly, so supporting data is not available. Requests to access the datasets should be directed to ZZ, zgzhang@szu.edu.cn.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of Institute of Psychology, Chinese Academy Sciences, the Ethics Committee of Liaoning Normal University, and the Ethics Committee of Shenzhen University. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin. Written informed consent was obtained from the individual(s), and minor(s)’ legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.

Author Contributions

ZZ, LL, and GH conceived of the presented idea. XZ developed the theory and performed the computations. ZL, LS, and LZ supervised the findings of this experiment. All authors discussed the results and contributed to the final manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant 81871443, in part by the Science, Technology and Innovation Commission of Shenzhen Municipality Technology Fund under Grant JCYJ20170818093322718, and in part by Shenzhen Peacock Plan under Grant KQTD2016053112051497.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2021.635733/full#supplementary-material

Footnotes

1. ^ www.fil.ion.ucl.ac.uk/spm/software/spm12

References

Al-Wasity, S., Vogt, S., Vuckovic, A., and Pollick, F. E. (2021). Upregulation of Supplementary Motor Area Activation with fMRI Neurofeedback during Motor Imagery. Eneuro 2021:8.

Google Scholar

Apkarian, A. V., Bushnell, M. C., Treede, R.-D., and Zubieta, J.-K. (2005). Human brain mechanisms of pain perception and regulation in health and disease. Eur. J. Pain 9, 463–484. doi: 10.1016/j.ejpain.2004.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Baliki, M. N., Geha, P. Y., and Apkarian, A. V. (2009). Parsing pain perception between nociceptive representation and magnitude estimation. J. Neurophysiol. 101, 875–887. doi: 10.1152/jn.91100.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Boly, M., Faymonville, M.-E., Schnakers, C., Peigneux, P., Lambermont, B., Phillips, C., et al. (2008). Perception of pain in the minimally conscious state with PET activation: an observational study. Lancet Neurol. 7, 1013–1020. doi: 10.1016/s1474-4422(08)70219-9

CrossRef Full Text | Google Scholar

Borsook, D., Sava, S., and Becerra, L. (2010). The pain imaging revolution: advancing pain into the 21st century. Neurosci. 16, 171–185. doi: 10.1177/1073858409349902

PubMed Abstract | CrossRef Full Text | Google Scholar

Bromm, B. (2001). Brain images of pain. Physiology 16, 244–249. doi: 10.1152/physiologyonline.2001.16.5.244

PubMed Abstract | CrossRef Full Text | Google Scholar

Calhoun, V., Stevens, M., Pearlson, G., and Kiehl, K. (2004). fMRI analysis with the general linear model: removal of latency-induced amplitude bias by incorporation of hemodynamic derivative terms. Neuroimage 22, 252–257. doi: 10.1016/j.neuroimage.2003.12.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Collins, D. L., Zijdenbos, A. P., Kollokian, V., Sled, J. G., Kabani, N. J., Holmes, C. J., et al. (1998). Design and construction of a realistic digital brain phantom. IEEE Transac. Med. Imag. 17, 463–468. doi: 10.1109/42.712135

CrossRef Full Text | Google Scholar

Davis, K. D., Wood, M. L., Crawley, A. P., and Mikulis, D. J. (1995). fMRI of human somatosensory and cingulate cortex during painful electrical nerve stimulation. Neuroreport 7, 321–325. doi: 10.1097/00001756-199512290-00077

CrossRef Full Text | Google Scholar

Devinsky, O., Morrell, M. J., and Vogt, B. A. (1995). Contributions of anterior cingulate cortex to behaviour. Brain 118, 279–306. doi: 10.1093/brain/118.1.279

PubMed Abstract | CrossRef Full Text | Google Scholar

Downar, J., Crawley, A. P., Mikulis, D. J., and Davis, K. D. (2000). A multimodal cortical network for the detection of changes in the sensory environment. Nat. Neurosci. 3, 277–283. doi: 10.1038/72991

PubMed Abstract | CrossRef Full Text | Google Scholar

Dum, R. P., Levinthal, D. J., and Strick, P. L. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Eccles, J. C. (1982). The initiation of voluntary movements by the supplementary motor area. Archiv für Psychiatr. Nervenkrankheiten 231, 423–441. doi: 10.1007/bf00342722

PubMed Abstract | CrossRef Full Text | Google Scholar

Fulbright, R. K., Troche, C. J., Skudlarski, P., Gore, J. C., and Wexler, B. E. (2001). Functional MR imaging of regional brain activation associated with the affective experience of pain. Am. J. Roentgenol. 177, 1205–1210. doi: 10.2214/ajr.177.5.1771205

PubMed Abstract | CrossRef Full Text | Google Scholar

Hutchison, W. D., Davis, K. D., Lozano, A. M., Tasker, R. R., and Dostrovsky, J. O. (1999). Pain-related neurons in the human cingulate cortex. Nat. Neurosci. 2, 403–405. doi: 10.1038/8065

PubMed Abstract | CrossRef Full Text | Google Scholar

Iannetti, G., Zambreanu, L., Cruccu, G., and Tracey, I. (2005). Operculoinsular cortex encodes pain intensity at the earliest stages of cortical processing as indicated by amplitude of laser-evoked potentials in humans. Neuroscience 131, 199–208. doi: 10.1016/j.neuroscience.2004.10.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Iannetti, G. D., and Mouraux, A. (2010). From the neuromatrix to the pain matrix (and back). Experimental brain research 205, 1–12. doi: 10.1007/s00221-010-2340-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Iannetti, G. D., Salomons, T. V., Moayedi, M., Mouraux, A., and Davis, K. D. (2013). Beyond metaphor: contrasting mechanisms of social and physical pain. Trends Cogn. Sci. 17, 371–378. doi: 10.1016/j.tics.2013.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Ingvar, M. (1999). Pain and functional imaging. Philosophical Transactions of the Royal Society of London. Series B: Biol. Sci. 354, 1347–1358. doi: 10.1098/rstb.1999.0483

PubMed Abstract | CrossRef Full Text | Google Scholar

Jackson, J. B., O’Daly, O., Makovac, E., Medina, S., de Lara Rubio, A., McMahon, S. B., et al. (2020). Noxious pressure stimulation demonstrates robust, reliable estimates of brain activity and self-reported pain. NeuroImage 221:117178. doi: 10.1016/j.neuroimage.2020.117178

PubMed Abstract | CrossRef Full Text | Google Scholar

Jensen, K. B., Regenbogen, C., Ohse, M. C., Frasnelli, J., Freiherr, J., and Lundström, J. N. (2016). Brain activations during pain: a neuroimaging meta-analysis of patients with pain and healthy controls. Pain 157, 1279–1286. doi: 10.1097/j.pain.0000000000000517

PubMed Abstract | CrossRef Full Text | Google Scholar

Legrain, V., Iannetti, G. D., Plaghki, L., and Mouraux, A. (2011). The pain matrix reloaded: a salience detection system for the body. Prog. Neurobiol. 93, 111–124. doi: 10.1016/j.pneurobio.2010.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Liberati, G., Algoet, M., Santos, S. F., Ribeiro-Vaz, J. G., Raftopoulos, C., and Mouraux, A. (2019). Tonic thermonociceptive stimulation selectively modulates ongoing neural oscillations in the human posterior insula: Evidence from intracerebral EEG. NeuroImage 188, 70–83. doi: 10.1016/j.neuroimage.2018.11.059

PubMed Abstract | CrossRef Full Text | Google Scholar

Liberati, G., Klöcker, A., Safronova, M. M., Ferrao Santos, S., Ribeiro Vaz, J.-G., Raftopoulos, C., et al. (2016). Nociceptive local field potentials recorded from the human insula are not specific for nociception. PLoS Biol. 14:e1002345. doi: 10.1371/journal.pbio.1002345

PubMed Abstract | CrossRef Full Text | Google Scholar

Loeser, J. D., and Treede, R.-D. (2008). The Kyoto protocol of IASP basic pain Terminology. Pain 137, 473–477. doi: 10.1016/j.pain.2008.04.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Lotze, M., Montoya, P., Erb, M., Hülsmann, E., Flor, H., Klose, U., et al. (1999). Activation of cortical and cerebellar motor areas during executed and imagined hand movements: an fMRI study. J. Cogn. Neurosci. 11, 491–501. doi: 10.1162/089892999563553

PubMed Abstract | CrossRef Full Text | Google Scholar

Lui, F., Duzzi, D., Corradini, M., Serafini, M., Baraldi, P., and Porro, C. A. (2008). Touch or pain? Spatio-temporal patterns of cortical fMRI activity following brief mechanical stimuli. Pain 138, 362–374. doi: 10.1016/j.pain.2008.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Macaluso, E., and Driver, J. (2005). Multisensory spatial interactions: a window onto functional integration in the human brain. Trends Neurosci. 28, 264–271. doi: 10.1016/j.tins.2005.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Martino, J., Gabarrós, A., Deus, J., Juncadella, M., Acebes, J., Torres, A., et al. (2011). Intrasurgical mapping of complex motor function in the superior frontal gyrus. Neuroscience 179, 131–142. doi: 10.1016/j.neuroscience.2011.01.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Mazzola, L., Isnard, J., Peyron, R., Guenot, M., and Mauguiere, F. (2009). Somatotopic organization of pain responses to direct electrical stimulation of the human insular cortex. Pain 146, 99–104. doi: 10.1016/j.pain.2009.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Mazzola, L., Isnard, J., Peyron, R., and Mauguière, F. (2012). Stimulation of the human cortex and the experience of pain: Wilder Penfield’s observations revisited. Brain 135, 631–640. doi: 10.1093/brain/awr265

PubMed Abstract | CrossRef Full Text | Google Scholar

Mazzola, L., Mauguiere, F., and Isnard, J. (2019). Functional mapping of the human insula: data from electrical stimulations. Revue neurologique 175, 150–156. doi: 10.1016/j.neurol.2018.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Mischkowski, D., Palacios-Barrios, E. E., Banker, L., Dildine, T. C., and Atlas, L. Y. (2018). Pain or nociception? Subjective experience mediates the effects of acute noxious heat on autonomic responses. Pain 159:699. doi: 10.1097/j.pain.0000000000001132

PubMed Abstract | CrossRef Full Text | Google Scholar

Moisset, X., and Bouhassira, D. (2007). Brain imaging of neuropathic pain. Neuroimage 37, S80–S88.

Google Scholar

Morrow, J., Mosher, C., and Gothard, K. (2019). Multisensory neurons in the primate amygdala. J. Neurosci. 39, 3663–3675. doi: 10.1523/jneurosci.2903-18.2019

PubMed Abstract | CrossRef Full Text | Google Scholar

Mouraux, A., Diukova, A., Lee, M. C., Wise, R. G., and Iannetti, G. D. (2011). A multisensory investigation of the functional significance of the “pain matrix”. Neuroimage 54, 2237–2249. doi: 10.1016/j.neuroimage.2010.09.084

PubMed Abstract | CrossRef Full Text | Google Scholar

Mouraux, A., and Iannetti, G. D. (2009). Nociceptive laser-evoked brain potentials do not reflect nociceptive-specific neural activity. J. Neurophysiol. 101, 3258–3269. doi: 10.1152/jn.91181.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Nichols, T., Brett, M., Andersson, J., Wager, T., and Poline, J.-B. (2005). Valid conjunction inference with the minimum statistic. Neuroimage 25, 653–660. doi: 10.1016/j.neuroimage.2004.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Pernet, C. R. (2014). Misconceptions in the use of the General Linear Model applied to functional MRI: a tutorial for junior neuro-imagers. Front. Neurosci. 8:1. doi: 10.15728/bbr.2004.1.1.1

CrossRef Full Text | Google Scholar

Peyron, R., Frot, M., Schneider, F., Garcia-Larrea, L., Mertens, P., Barral, F., et al. (2002). Role of operculoinsular cortices in human pain processing: converging evidence from PET, fMRI, dipole modeling, and intracerebral recordings of evoked potentials. Neuroimage 17, 1336–1346. doi: 10.1006/nimg.2002.1315

PubMed Abstract | CrossRef Full Text | Google Scholar

Peyron, R., Laurent, B., and Garcia-Larrea, L. (2000). Functional imaging of brain responses to pain. A review and meta-analysis (2000). Clin. Neurophysiol. 30, 263–288. doi: 10.1016/s0987-7053(00)00227-6

CrossRef Full Text | Google Scholar

Sanchez, G., Hartmann, T., Fuscà, M., Demarchi, G., and Weisz, N. (2020). Decoding across sensory modalities reveals common supramodal sil signatures of conscious perception. Proc. Natl. Acad. Sci. 117, 7437–7446. doi: 10.1073/pnas.1912584117

PubMed Abstract | CrossRef Full Text | Google Scholar

Segerdahl, A. R., Mezue, M., Okell, T. W., Farrar, J. T., and Tracey, I. (2015). The dorsal posterior insula subserves a fundamental role in human pain. Nat. Neurosci. 18, 499–500. doi: 10.1038/nn.3969

PubMed Abstract | CrossRef Full Text | Google Scholar

Stern, J., Jeanmonod, D., and Sarnthein, J. (2006). Persistent EEG overactivation in the cortical pain matrix of neurogenic pain patients. Neuroimage 31, 721–731. doi: 10.1016/j.neuroimage.2005.12.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanji, J. (1996). New concepts of the supplementary motor area. Curr. Opin. Neurobiol. 6, 782–787. doi: 10.1016/s0959-4388(96)80028-6

CrossRef Full Text | Google Scholar

Tanji, J., and Shima, K. (1994). Role for supplementary motor area cells in planning several movements ahead. Nature 371, 413–416. doi: 10.1038/371413a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Tölle, T. R., Kaufmann, T., Siessmeier, T., Lautenbacher, S., Berthele, A., Munz, F., et al. (1999). Region−specific encoding of sensory and affective components of pain in the human brain: a positron emission tomography correlation analysis. Ann. Neurol. 45, 40–47. doi: 10.1002/1531-8249(199901)45:1<40::aid-art8>3.0.co;2-l

CrossRef Full Text | Google Scholar

Tracey, I., and Mantyh, P. W. (2007). The cerebral signature for pain perception and its modulation. Neuron 55, 377–391. doi: 10.1016/j.neuron.2007.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Wager, T. D., Atlas, L. Y., Lindquist, M. A., Roy, M., Woo, C.-W., and Kross, E. (2013). An fMRI-based neurologic signature of physical pain. N. Engl. J. Med. 368, 1388–1397. doi: 10.1056/nejmoa1204471

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, R., Wang, F., Yang, P.-F., and Chen, L. M. (2017). High-resolution functional MRI identified distinct global intrinsic functional networks of nociceptive posterior insula and S2 regions in squirrel monkey brain. Neuroimage 155, 147–158. doi: 10.1016/j.neuroimage.2017.04.067

PubMed Abstract | CrossRef Full Text | Google Scholar