Skip to main content
  • Research article
  • Open access
  • Published:

Individual pain sensitivity is associated with resting-state cortical activities in healthy individuals but not in patients with migraine: a magnetoencephalography study

Abstract

Background

Pain sensitivity may determine the risk, severity, prognosis, and efficacy of treatment of clinical pain. Magnetic resonance imaging studies have linked thermal pain sensitivity to changes in brain structure. However, the neural correlates of mechanical pain sensitivity remain to be clarified through investigation of direct neural activities on the resting-state cortical oscillation and synchrony.

Methods

We recorded the resting-state magnetoencephalographic (MEG) activities of 27 healthy individuals and 30 patients with episodic migraine (EM) and analyzed the source-based oscillatory powers and functional connectivity at 2 to 59 Hz in pain-related cortical regions, which are the bilateral anterior cingulate cortex (ACC), medial orbitofrontal (MOF) cortex, lateral orbitofrontal (LOF) cortex, insula cortex, primary somatosensory cortex (SI), primary motor cortex (MI), and posterior cingulate cortex (PCC). The mechanical punctate pain threshold (MPPT) was obtained at the supraorbital area (the first branch of the trigeminal nerve dermatome, V1) and the forearm (the first thoracic nerve dermatome, T1) and further correlated with MEG measures.

Results

The MPPT is inversely correlated with the resting-state relative powers of gamma oscillation in healthy individuals (all corrected P < 0.05). Specifically, inverse correlation was noted between the MPPT at V1 and gamma powers in the bilateral insula (r = − 0.592 [left] and − 0.529 [right]), PCC (r = − 0.619 and − 0.541) and MI (r = − 0.497 and − 0.549) and between the MPPT at T1 and powers in the left PCC (r = − 0.561) and bilateral MI (r = − 0.509 and − 0.520). Furthermore, resting-state functional connectivity at the delta to beta bands, especially between frontal (MOF, ACC, LOF, and MI), parietal (PCC), and sensorimotor (bilateral SI and MI) regions, showed a positive correlation with the MPPT at V1 and T1 (all corrected P < 0.05). By contrast, in patients with EM, the MPPT was not associated with resting-state cortical activities.

Conclusions

Pain sensitivity in healthy individuals is associated with the resting-state gamma oscillation and functional connectivity in pain-related cortical regions. Further studies must be conducted in a large population to confirm whether resting-state cortical activities can be an objective measurement of pain sensitivity in individuals without clinical pain.

Introduction

Pain sensitivity may determine the risk, severity, prognosis, and treatment efficacy of clinical pain [1, 2]. Reduced pain sensitivity may delay recognition and undermine treatment efficacy in acute pain [3], whereas elevated pain sensitivity may increase health care costs and susceptibility to chronic pain conditions [4]. In the era of precision medicine, objective assessment of pain sensitivity at the individual level is an unmet need. Pain sensitivity is highly variable across individuals. Notably, genetic [5, 6], environmental [7, 8], and psychological [9, 10] factors influence individual pain perception [11].

Some studies have explored the neural correlates of pain sensitivity. Studies involving structural magnetic resonance imaging (MRI) have identified that cortical thickness or volume density change in the pain-related cortical network may contribute to varying pain sensitivity across healthy individuals. More specifically, high pain sensitivity was associated with cortical thickening in the primary somatosensory cortex (SI), posterior cingulate cortex (PCC), and orbitofrontal cortex [12]. Additionally, the intensity rating of thermal pain was negatively correlated with grey matter density in the SI, PCC, precuneus, intraparietal sulcus, and inferior parietal cortex [13]. In line with the aforementioned MRI studies, neurophysiological studies obtaining direct neural signals showed that noxious stimuli activated a widely distributed brain network related to pain processing, including the SI, primary motor cortex (MI), insula, anterior cingulate cortex (ACC), medial frontal cortex, and PCC [14]. Moreover, activations among some of these regions were greater in pain-sensitive individuals than in pain-insensitive ones [15], and the magnitude increased with increasing stimulus intensity or perceived magnitude of pain intensity [16]. The aforementioned findings suggest the pain-related cortical network as the structural correlate of individual thermal pain sensitivity, and noxious-evoked neural oscillations and synchronizations in these regions may cause the intersubject variability in pain perception. However, how the neural correlates of spontaneous cortical activities with the mechanical pain sensitivity and in clinical pain remains elusive.

The rating of individual pain in response to noxious stimulus has been coded from cortical activation [16,17,18,19], and pain intensity during heat stimulation was related to gamma oscillation in the medial frontal cortex [19]. In another study, during capsaicin-heat pain, peak alpha frequency over the sensorimotor region was inversely correlated with individual pain intensity [18]. The inconsistent findings regarding cortical regions and oscillations may be attributable to the influence of the pain modalities and the involvement of cognitive processes, such as the salience or attention effect. Therefore, some studies have investigated prestimulus functional connectivity instead of stimulus-evoked responses and found that the connectivity of the anterior insula cortex [20] and frontocentral network [21] determined the pain perception of the subsequent noxious stimulus. Thus, brain activities or synchrony might involve the neurophysiological mechanisms for individual pain perception. Moreover, brain oscillations and synchrony serve integrative functions through flexibly regulating information flow among the cortical regions [22,23,24]. Thus, exploring oscillation and synchrony during the resting-state condition might yield promising insights into how functionally diverse processes relevantly reflect the intersubject variability of pain sensitivity.

This study investigates the hypothesis that resting-state cortical activities at the pain-related cortical network underpin interindividual pain variability. To characterize temporal–spatial features of cortical oscillations and cortico-cortical synchronization within this pain-related network, the present study used magnetoencephalography (MEG) to record brain activities during the resting-state condition. Moreover, this study assessed individual pain sensitivity with the mechanical punctate pain threshold (MPPT) instead of the thermal pain threshold used in most of the earlier pertinent studies. Furthermore, we recruited patients with episodic migraine (EM) to determine whether the study findings are exclusive to patients with pain disorder. We selected EM here because patients with EM have been characterized as having an aberrant pain sensitivity threshold [25], heightened cortical excitability [26, 27], and altered resting-state cortical oscillations and connectivity in pain-related regions [28,29,30]. Moreover, whether the underlying pain sensitivity mechanism is reshaped for pain disorder remains uncertain. The specific aims of this study were to (1) elucidate the relationship between cortical oscillations and pain sensitivity, (2) investigate the correlation of cortical synchronization with pain sensitivity, and (3) examine the effect of pain disorder in patients with EM on the cortical mechanism of pain sensitivity.

Methods

Participants

This study recruited healthy individuals who did not have medical or family histories of pain disorders and had not experienced any significant pain condition during the past year. Furthermore, patients with EM (monthly headache days: 1–14) were enrolled from the Headache Clinic of Taipei Veterans General Hospital (VGHTPE), and the diagnosis (code 1.1: migraine without aura) was according to the criteria of the International Classification of Headache Disorders, Third Edition [31]. All participants were right-handed, denied having any history of systemic or major neuropsychiatric disease, and had normal physical and neurological examination results as well as normal brain MRI results. Participants who were taking any medication (eg, migraine preventive medications) on a daily basis were excluded. Furthermore, patients with migraine who overused medication were excluded. The hospital’s institutional review board (IRB) approved the study protocol (VGHTPE-IRB, 2015–10-001 BC), and each participant provided written informed consent.

All participants underwent scheduled MPPT and MEG recordings (both detailed below) and were instructed not to take any analgesics or other medications within 3 days before the recordings. Patients with migraine were assessed during their interictal phase, which was defined as the absence of acute migraine within 2 days before (days − 1 and − 2) and after (days + 1 and + 2) the aforementioned assessment. We evaluated the severity of depression and anxiety using the Hospital Anxiety and Depression Scale (HADS) and functional disability in migraine using the Migraine Disability Assessment Scale (MIDAS).

Pain sensitivity measurement

Pain sensitivity measurement was performed in the chronic pain examination room (constant 20 °C room temperature and no windows) in the Department of Neurology at VGHTPE. Quantitative measurements of the MPPT were defined as the lowest intensity perceived as painful for participants, and MPPT values were obtained as follows. MPTT was determined using the standard rigid electronic von Frey device (ALMEMO 2450, AHLBORN, Germany), which consists of a 1000-g internal load cell connected to a probe and a rigid tip (diameter: 0.8 mm; IITC Life Science Inc., USA). The system digitally measures and displays the values in grams, with a resolution of 0.1 g. During measurement, the tip was applied perpendicularly to the skin surface, with ascending stimulus intensity applied at 25 g/s [32]. Participants were instructed to inform the examiner immediately upon experiencing a painful sensation. Assessment target areas included the left supraorbital (ie, the first branch of the trigeminal nerve dermatome, V1) and proximal medio-ventral forearm (ie, the first thoracic nerve dermatome, T1). Breaks were taken between each stimulus, and the order of stimulation sites was randomized. Sensory stimuli on each target area were delivered 5 times; the 2 extreme values were excluded, and the average values of the remaining 3 recordings were calculated for further analyses.

Resting-state MEG recording

A 5-min resting-state MEG recording [33, 34] was obtained for each participant, each of whom sat comfortably with eyes closed but remained awake and relaxed. If the participant fell asleep or had excessive within-run head movement, the recording was stopped and then rerun. Furthermore, a 3-min empty-room recording was conducted to capture sensor and environmental noises, which were applied to calculate the noise covariance for further source model analysis. MEG data were recorded with the digitization rate of 600 Hz using a whole-scalp 306-channel neuromagnetometer (Vectorview; Elekta Neuromag, Helsinki, Finland) composed of 102 identical triple sensor elements (1 magnetometer and 2 orthogonal planar gradiometers). In total, 4 coils representing the head position were placed on the participant’s scalp, specified by the nasion and 2 preauricular points using Cartesian coordinates and measured with a 3-dimensional (3D) digitizer. For accurate registration, approximately 50 additional scalp points were digitized. These head landmarks and points enabled further coordinate alignment between the MEG and MRI data. Additionally, electrooculography (EOG) and electrocardiography (ECG) activities during MEG recording were simultaneously recorded for data preprocessing. MRI images of individual brain structures were acquired using a 3 T MR system (Siemens Magnetom Tim Trio), with a TR of 9.4 ms, TE of 4 ms, recording matrix of 256 × 256 pixels, field of view of 256 mm, and slice thickness of 1 mm.

Data preprocessing

To obtain intrinsic spontaneous cortical activities and reject the contaminations of nonbrain or environmental artifacts, MEG data were preprocessed as follows (Fig. 1): (1) MaxFilter from the Neuromag software system was applied to remove external noise from MEG recordings [35, 36], (2) the data segments containing artifacts from head movement, muscle activities or environmental noise were manually rejected, (3) notch filters (60 Hz and its harmonics) were used to remove powerline contaminations, (4) identified heartbeat and eye blinking events from ECG and EOG data were used to define the projectors through principal component analysis separately. The principal components meeting the artifact’s sensor topography were then manually excluded through orthogonal projection [37].

Fig. 1
figure 1

Pipeline of source-based resting-state magnetoencephalographic spectral and functional connectivity analysis

Furthermore, for further source modeling analysis, T1-weighted structural volumetric images were automatically reconstructed into a surface model using BrainVISA (4.5.0, http://brainvisa.info). The detailed geometric reconstruction of the scalp, brain gray and white matter, and tessellations provided a topographical 3D representation of the brain surface and was used to estimate the gray and white matter border.

Source-based spectral power and functional connectivity analysis

In combination with a structural surface model, the distributed source model of resting-state MEG data was estimated using depth-weighted minimum norm estimation analysis (Fig. 1). A forward model was established by overlapping sphere method for rapid calculation of realistic head shapes [38], which presented each cortical vertex as a current dipole and included approximately 15,000 vertices in the whole brain model. The inverse operator estimated the distribution of current sources that account for data recorded at the MEG sensors. For group analysis, the cortical source model of each participant was then morphed into a common source space defined by the Colin27 anatomy [39]. Through dynamic source analysis, the current density of cortical activation in each individual could be obtained at any cortical region and any time point.

In this study, the regions of interest (ROIs) were defined in the T1 template volume using Mindboggle cortical parcellation [40] and were composed of 14 brain regions, including the bilateral ACC, medial orbitofrontal cortex (MOF), lateral orbitofrontal cortex (LOF), insula cortex, SI, MI, and PCC. These ROIs were involved in the sensory and affective aspects of cortical pain representation and were related to pain sensitivity in terms of the structural volume or thickness [12,13,14]. Besides, the bilateral auditory cortices, apart from the pain regions, were also selected as negative control. The time-varying source density of each ROI was individually derived from the averaged source density of each vertex within the ROI. To obtain the resting-state oscillations and functional connectivity, the dynamic source density of each ROI was further analyzed as follows. First, oscillatory power at each ROI was estimated using the Welch method (window duration: 5 s with 50% overlap) and defined as absolute power. Moreover, the oscillatory power was normalized to relative power through the division of the power at each frequency band by the total power, which has been reported to adequately reduce the interindividual variability of the oscillatory magnitude [41]. Second, the amplitude envelope correlation analysis [42], which orthogonalized the signals to remove zero-lag interactions [43], was used to calculated the oscillatory functional connectivity between ROIs and then constructed the full 14 × 14 adjacency matrix. Oscillatory power and functional connectivity were categorized according to frequency bands: delta (2–4 Hz), theta (5–7 Hz), alpha (8–13 Hz), beta (14–29 Hz), and gamma (30–59 Hz), and averaged in each frequency range. The MEG data preprocessing and analysis were performed using Brainstorm [44].

Statistical analysis

The demographics and MPPT of the healthy individuals and EM were compared using Student’s t test or a chi-square test as appropriate. Pearson’s correlation was used to determine the correlation between the MPPT (V1 and T1) and cortical oscillations (absolute and relative) and synchronization (14 × 14 ROIs) at the delta to gamma frequency bands, as well as the oscillations and synchronization in the bilateral auditory cortex. Throughout the statistical analyses, false discovery rate correction (FDR) was used for multiple comparisons, and a P value of < 0.05 was considered statistically significant.

Results

Demographics and pain sensitivity responses

This study evaluated 57 participants (27 healthy individuals and 30 patients with EM). The groups did not differ in age and sex (both P > 0.05; Table 1). The MPPT at V1 was higher in patients with EM than in healthy individuals (P = 0.041); however, no difference was observed at T1. Moreover, the anxiety and depression scores were higher in patients with EM than in healthy individuals (all P < 0.01).

Table 1 Subject Characteristics

Correlations between resting-state cortical oscillations and the MPPT

In healthy individuals, absolute cortical powers within 14 pain-related regions did not correlate with the MPPT at V1 and T1 at 5 frequency bands (all corrected P > 0.05; Fig. 2a). However, relative power at the gamma band was significantly linked to the MPPT (Fig. 2b), but not at other frequency bands (delta to beta bands). Notably, significant negative correlations between gamma power and the MPPT were observed between MPPTs at V1 in the bilateral insula (left, r = − 0.592; right, r = − 0.529), PCC (left, r = − 0.619; right, r = − 0.541), and MI regions (left, r = − 0.497; right, r = − 0.549) and between MPPT at T1 and power in the left PCC (r = − 0.561) and bilateral MI regions (left, r = − 0.509; right, r = − 0.52; Fig. 3). These relationships indicated that the high gamma power of resting-state cortical activities represented a low pain threshold, that is, high pain sensitivity. Remarkably, in the auditory cortex, no clear correlation was noted between cortical oscillations and MPPT (Supplementary Table 1).

Fig. 2
figure 2

Relationship of the mechanical punctate pain threshold with the (a) absolute power and (b) relative power at the delta to gamma bands in the pain-related regions in healthy individuals. The P-value of the correlation analysis is color-coded, with a significant value denoted with yellow. *P < .05 with false discovery rate correction. MPPT, mechanical punctate pain threshold; ACC, anterior cingulate cortex; LOF, lateral orbitofrontal cortex; MOF, medial orbitofrontal cortex; SI, primary somatosensory cortex; PCC, posterior cingulate cortex; MI, primary motor cortex; V1, left supraorbital (the first branch of the trigeminal nerve dermatome); T1, proximal medio-ventral forearm (the first thoracic nerve dermatome)

Fig. 3
figure 3

Significant correlation of the mechanical punctate pain threshold with resting-state gamma activities at the (a) left insula, (b) right insula, (c) left PCC, (d) right PCC, (e) left MI, and (f) right MI in healthy individuals. *P < .05 with false discovery rate correction. MPPT, mechanical punctate pain threshold; PCC, posterior cingulate cortex; MI, primary motor cortex; V1, left supraorbital (the first branch of the trigeminal nerve dermatome); T1, proximal medio-ventral forearm (the first thoracic nerve dermatome)

Correlations of the MPPT with resting-state functional connectivity

Significant positive correlations were observed between functional connectivity and the MPPT (P < 0.05 with FDR correction), indicating that strong cortical synchronization exhibited a high pain threshold. For the MPPT at V1, strong correlations were noted at the delta, theta, alpha, or beta bands within the connectivity between pain-related regions, which mainly included the cortical areas of salience network (insula and ACC), sensorimotor network (SI and MI) and some areas in the default mode network (PCC and MOF) (predominant connections are detailed in Fig. 4 (a)). Regarding the MPPT at T1, significant correlations were observed at the alpha and beta bands mainly within the connectivity at ACC and sensorimotor network (Fig. 4 (b)). Notably, the connectivity between auditory cortex at each frequency band did not correlate with the MPPT (Supplementary Table 2).

Fig. 4
figure 4

In healthy individuals, a The relationship between mechanical punctate pain threshold at V1 and functional connectivity within pain-related regions. Significant correlations are noted at the delta to beta bands; b The relationship between mechanical punctate pain threshold at T1 and functional connectivity within pain-related regions. Significant correlations are observed at the alpha and beta bands. Each node presents the centroid of each region of interest. The lower part shows the correlation coefficients of the significant connectivity. MPPT, mechanical punctate pain threshold; ACC, anterior cingulate cortex; LOF, lateral orbitofrontal cortex; MOF, medial orbitofrontal cortex; SI, primary somatosensory cortex; PCC, posterior cingulate cortex; MI, primary motor cortex; V1, left supraorbital (the first branch of the trigeminal nerve dermatome)

Correlation of pain sensitivity with resting-state activities in EM

By contrast, in patients with EM, the power of cortical oscillations did not correlate with the MPPT at all frequency bands within the 14 ROIs (all corrected P > 0.05; Fig. 5). Similarly, resting-state functional connectivity between all ROIs at all frequency bands did not correlate with the MPPT (all P > 0.05).

Fig. 5
figure 5

In patients with episodic migraine, no correlation between the mechanical punctate pain threshold and relative resting-state power. MPPT, mechanical punctate pain threshold; ACC, anterior cingulate cortex; LOF, lateral orbitofrontal cortex; MOF, medial orbitofrontal cortex; SI, primary somatosensory cortex; PCC, posterior cingulate cortex; MI, primary motor cortex; V1, left supraorbital (the first branch of the trigeminal nerve dermatome); T1, proximal medio-ventral forearm (the first thoracic nerve dermatome)

Discussion

This MEG study demonstrated that the MPPT in healthy individuals was associated with resting-state cortical oscillation and synchronization in the pain-related cortical regions, but not in the auditory cortex. Specifically, the MPPT was inversely correlated with the relative power of gamma oscillation in the bilateral insula, PCC, and MI regions and positively correlated with cortical synchronization within the pain-related cortical network. By contrast, patients with EM did not exhibit such a relationship between pain sensitivity and cortical activities.

Association between resting-state gamma oscillation and pain sensitivity

In the present study, large gamma activities were linked to low MPPT (high pain sensitivity), implying that the resting-state gamma oscillation might be a signature of pain sensitivity. The association is in agreement with the functional roles of gamma oscillation in pain perception [45], attentional effect of pain [46], and subjective pain intensity [47]. Gamma power activities have been used to reliably encoded subjective pain ratings to ongoing experimental pain in healthy participants [48, 49] and the ongoing pain intensity in patients with chronic back pain [50]. Moreover, enhanced resting-state gamma oscillatory activity in the prefrontal cortex and cerebellum was suggested to be a characteristic marker in patients with chronic neuropathic pain [51]. In a recent animal study, gamma oscillation could increase the recruitment of neural activation in a large network between cortical and subcortical structures, and thus enhance pain perception [52]. Taken together, pain perception may be characterized by resting-state gamma oscillation, and elevated gamma activities may facilitate the neural coupling of nociceptive input.

This study further found that gamma activities underlying mechanical pain perception were observed predominantly in the bilateral insula, PCC, and MI. Moreover, a multivariate analysis of laser-evoked pain-related EEG activities suggested the objective neural correlate of individual pain sensitivity from FCz and Cz electrodes [49]. Another electrophysiological approach in humans and rodents demonstrated that the gamma event-related synchronizations sampled using Cz electrodes reliably predict pain sensitivity across individuals [48]. The findings that used EEG FCz and Cz electrodes to coded pain sensitivity [48, 49] could imply that the nociceptive-related insula and MI regions may involve the pain sensitivity mechanism. This partly agreed with the fMRI findings that indicated that insula activation coded pain intensity [53]. Additionally, regarding the difference in brain structure in the PCC, pain sensitivity was related to cortical thickness [12] and gray matter density [13]. Because of discrepancies in the recording (EEG, MRI, or fMRI vs MEG) and experimental design (task-evoked vs spontaneous activities), cortical involvement in pain sensitivity remains debatable; however, spontaneous gamma activities at the insula, MI, and PCC are suggested to be neural signatures underlying the pain sensitivity mechanism.

Brain synchrony coupling pain sensitivity

The MPPT was coupled with cortical synchronization within the pain-related network. This echoed previous findings that prestimulus frontoparietal connectivity [21] or insular-brainstem connectivity [20] was related to subsequent pain perception. Moreover, a recent resting-state fMRI study determined that the medial-frontal and frontal-parietal networks were linked to the pain threshold [54]. Regarding the cortical function in pain stimulation, induced hemodynamic responses in the SI, ACC, and prefrontal regions reflected the subjective experience of pain [15]. Moreover, the functional connectivity in the prefrontal cortex and PCC, which are engaged in executive control and the default mode network, respectively, was associated with change in pain sensitivity for sleep disruption [55]. Similarly, the relationship between differences in brain structure and pain sensitivity was observed in the PCC, SI, and frontal regions [12, 13]. The interaction between prefrontal or ACC and sensorimotor regions has been suggested to characterize the affective or cognitive modulation of pain perception [56, 57]. Moreover, the nociceptive process was shaped from the interaction within the sensorimotor network [58, 59]. Accordingly, individual differences in pain sensitivity might result from the integration of sensory, affective, and cognitive states of the individual, indicating that oscillatory connectivity in the pain-related network could underpin pain sensitivity.

Aberrant neural mechanism for pain sensitivity in EM

The elevated MPPT in EM echoed our recent published findings [60], suggesting the altered pain sensitivity for the neuropathology of migraine, which was further evidenced for the disassociation between the MPPT and resting-state cortical oscillations or cortical synchronizations in the present study. Notably, in consistent with our findings, Schwedt and Chong [61] examined the correlations between heat pain thresholds and cortical thickness using functional MRI, and revealed a significant negative correlation in healthy subjects but not in patients with migraine. In part, this dissociation may reflect the decreased cortical excitability [26, 27] during the interictal period and abnormal resting-state cortical oscillations and connectivity in pain-related regions [28,29,30] in patients with migraine. Specifically, in a resting-state condition, patients with migraine were associated with reorganized functional connectivity in the default mode, salience, or sensorimotor network [29, 62], implying that the functions of the resting brain were reshaped by clinical pain.

Study limitations

One limitation to this study is the generalizability of the present findings to pain sensitivity with respect to other sensory modalities (eg, cold pain, hot pain, or pressure pain). Evidence suggests a close relationship between thermal and mechanical pain sensitivity [63]. However, a comparative study across sensory modalities is necessary, and this issue is beyond the scope of the present study. In addition, the processing of pain sensitivity may differ between bodily parts, as suggested by earlier psychophysical studies [64]. The varying spatial acuity of pain, innervation density, and receptive field size of nociceptors across different regions may be potential reasons for this phenomenon. Nevertheless, functional interactions within the pain-related regions considerably reflect pain sensitivity at these bodily parts.

While previous migraine studies [29, 65, 66] studying the relationship of high frequency oscillations with aberrant pain processing, the sampling rate (600 Hz) in this study limits the calculation of the high frequency oscillation (> 200 Hz). However, the power and functional connectivity at high gamma band (60–100 Hz) was further analyzed for the correlations with the MPPT. The results show no significant correlation at this frequency band (all corrected p > 0.05; Supplementary Tables 3–5). Furthermore, the analysis of brain activity of radial sources from cingulate cortex may be difficult with MEG. In previous MEG studies, task-related [67] and resting-state [68] ACC connectivity were clearly observed using depth-weighting MNE analysis and 306-channel whole-head Elekta Neuromag system (the same device as the present study). The reasons that we could clearly measure the MEG activity from cingulate cortex may be: (1) the source orientation in the cingulate cortex is not all radial; (2) depth-weighting algorithm alleviates the reduction of the radial source amplitude with increasing depth [69]. Herein, we delineated the spectral magnitude of the resting-state cortical activities from one healthy subject as shown in the Supplementary Fig. 1. The ACC activities are clearly discernible and the activation strengths are ranked in the middle across all cortical regions.

Finally, the dissociation between pain sensitivity and cortical activities in migraine is assumed to be due to the dysfunction of the resting-state pain network; however, the current study design was incapable of deciphering whether this dissociation in migraine is due to disease per se or its heterogeneous phenotype (headache profiles, psychiatric comorbidities, etc.).

Conclusion

Mechanical pain sensitivity in healthy individuals is associated with the resting-state gamma oscillation and functional connectivity within pain-related brain networks. The relationship between cortical activation and pain sensitivity is not noted in EM for cortical reorganization. Resting-state cortical activities may be a brain signature of pain sensitivity in healthy individuals, and confirmatory studies in a large population are necessary.

Availability of data and materials

The datasets analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

MEG:

Magnetoencephalography

EM:

Episodic migraine

ACC:

Anterior cingulate cortex

MOF:

Medial orbitofrontal

LOF:

Lateral orbitofrontal

SI:

Primary somatosensory cortex

MI:

Primary motor cortex

PCC:

Posterior cingulate cortex

MPPT:

Mechanical punctate pain threshold

V1:

First branch of the trigeminal nerve dermatome

T1:

First thoracic nerve dermatome

MIDAS:

Migraine Disability Assessment

HADS-A:

The anxiety subscale of the Hospital Anxiety and Depression Scale

HADS-D:

The depression subscale of the Hospital Anxiety and Depression Scale

MRI:

Magnetic resonance imaging

EOG:

Electrooculography

ECG:

Electrocardiography

ROIs:

Regions of interest

References

  1. Meints SM, Mawla I, Napadow V, Kong J, Gerber J, Chan ST, Wasan AD, Kaptchuk TJ, McDonnell C, Carriere J, Rosen B, Gollub RL, Edwards RR (2019) The relationship between catastrophizing and altered pain sensitivity in patients with chronic low-back pain. Pain 160:833–843

    Article  PubMed  PubMed Central  Google Scholar 

  2. Nahman-Averbuch H, Leon E, Hunter BM, Ding L, Hershey AD, Powers SW, King CD, Coghill RC (2019) Increased pain sensitivity but normal pain modulation in adolescents with migraine. Pain 160:1019–1028

    Article  PubMed  Google Scholar 

  3. Sheifer SE, Gersh BJ, Yanez ND 3rd, Ades PA, Burke GL, Manolio TA (2000) Prevalence, predisposing factors, and prognosis of clinically unrecognized myocardial infarction in the elderly. J Am Coll Cardiol 35:119–126

    Article  CAS  PubMed  Google Scholar 

  4. Diatchenko L, Slade GD, Nackley AG, Bhalang K, Sigurdsson A, Belfer I, Goldman D, Xu K, Shabalina SA, Shagin D, Max MB, Makarov SS, Maixner W (2005) Genetic basis for individual variations in pain perception and the development of a chronic pain condition. Hum Mol Genet 14:135–143

    Article  CAS  PubMed  Google Scholar 

  5. Kim H, Neubert JK, San Miguel A, Xu K, Krishnaraju RK, Iadarola MJ, Goldman D, Dionne RA (2004) Genetic influence on variability in human acute experimental pain sensitivity associated with gender, ethnicity and psychological temperament. Pain 109:488–496

    Article  PubMed  Google Scholar 

  6. Norbury TA, MacGregor AJ, Urwin J, Spector TD, McMahon SB (2007) Heritability of responses to painful stimuli in women: a classical twin study. Brain 130:3041–3049

    Article  PubMed  Google Scholar 

  7. Chesler EJ, Wilson SG, Lariviere WR, Rodriguez-Zas SL, Mogil JS (2002) Identification and ranking of genetic and laboratory environment factors influencing a behavioral trait, thermal nociception, via computational analysis of a large data archive. Neurosci Biobehav Rev 26:907–923

    Article  PubMed  Google Scholar 

  8. MacGregor AJ, Griffiths GO, Baker J, Spector TD (1997) Determinants of pressure pain threshold in adult twins: evidence that shared environmental influences predominate. Pain 73:253–257

    Article  CAS  PubMed  Google Scholar 

  9. Ochsner KN, Ludlow DH, Knierim K, Hanelin J, Ramachandran T, Glover GC, Mackey SC (2006) Neural correlates of individual differences in pain-related fear and anxiety. Pain 120:69–77

    Article  PubMed  Google Scholar 

  10. Seminowicz DA, Davis KD (2006) Cortical responses to pain in healthy individuals depends on pain catastrophizing. Pain 120:297–306

    Article  PubMed  Google Scholar 

  11. Nielsen CS, Staud R, Price DD (2009) Individual differences in pain sensitivity: measurement, causation, and consequences. J Pain 10:231–237

    Article  PubMed  Google Scholar 

  12. Erpelding N, Moayedi M, Davis KD (2012) Cortical thickness correlates of pain and temperature sensitivity. Pain 153:1602–1609

    Article  PubMed  Google Scholar 

  13. Emerson NM, Zeidan F, Lobanov OV, Hadsel MS, Martucci KT, Quevedo AS, Starr CJ, Nahman-Averbuch H, Weissman-Fogel I, Granovsky Y, Yarnitsky D, Coghill RC (2014) Pain sensitivity is inversely related to regional grey matter density in the brain. Pain 155:566–573

    Article  PubMed  Google Scholar 

  14. Apkarian AV, Bushnell MC, Treede RD, Zubieta JK (2005) Human brain mechanisms of pain perception and regulation in health and disease. Eur J Pain 9:463–484

    Article  PubMed  Google Scholar 

  15. Coghill RC, McHaffie JG, Yen YF (2003) Neural correlates of interindividual differences in the subjective experience of pain. Proc Natl Acad Sci U S A 100:8538–8542

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Coghill RC, Sang CN, Maisog JM, Iadarola MJ (1999) Pain intensity processing within the human brain: a bilateral, distributed mechanism. J Neurophysiol 82:1934–1943

    Article  CAS  PubMed  Google Scholar 

  17. Hsiao FJ, Chen WT, Liao KK, Wu ZA, Ho LT, Lin YY (2008) Oscillatory characteristics of nociceptive responses in the SII cortex. Can J Neurol Sci 35:630–637

    Article  PubMed  Google Scholar 

  18. Furman AJ, Meeker TJ, Rietschel JC, Yoo S, Muthulingam J, Prokhorenko M, Keaser ML, Goodman RN, Mazaheri A, Seminowicz DA (2018) Cerebral peak alpha frequency predicts individual differences in pain sensitivity. Neuroimage 167:203–210

    Article  PubMed  Google Scholar 

  19. Nickel MM, May ES, Tiemann L, Schmidt P, Postorino M, Ta Dinh S, Gross J, Ploner M (2017) Brain oscillations differentially encode noxious stimulus intensity and pain intensity. Neuroimage 148:141–147

    Article  PubMed  Google Scholar 

  20. Ploner M, Lee MC, Wiech K, Bingel U, Tracey I (2010) Prestimulus functional connectivity determines pain perception in humans. Proc Natl Acad Sci U S A 107:355–360

    Article  CAS  PubMed  Google Scholar 

  21. Taesler P, Rose M (2016) Prestimulus theta oscillations and connectivity modulate pain perception. J Neurosci 36:5026–5033

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Akam T, Kullmann DM (2014) Oscillatory multiplexing of population codes for selective communication in the mammalian brain. Nat Rev Neurosci 15:111–122

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Fries P (2015) Rhythms for cognition: communication through coherence. Neuron 88:220–235

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Schnitzler A, Gross J (2005) Normal and pathological oscillatory communication in the brain. Nat Rev Neurosci 6:285–296

    Article  CAS  PubMed  Google Scholar 

  25. Nahman-Averbuch H, Shefi T, Schneider VJ 2nd, Li D, Ding L, King CD, Coghill RC (2018) Quantitative sensory testing in patients with migraine: a systematic review and meta-analysis. Pain 159:1202–1223

    Article  PubMed  Google Scholar 

  26. Chen WT, Hsiao FJ, Ko YC, Liu HY, Wang PN, Fuh JL, Lin YY, Wang SJ (2018) Comparison of somatosensory cortex excitability between migraine and "strict-criteria" tension-type headache: a magnetoencephalographic study. Pain 159:793–803

    Article  PubMed  Google Scholar 

  27. Hsiao FJ, Wang SJ, Lin YY, Fuh JL, Ko YC, Wang PN, Chen WT (2018) Somatosensory gating is altered and associated with migraine chronification: a magnetoencephalographic study. Cephalalgia 38:744–753

    Article  PubMed  Google Scholar 

  28. Lee MJ, Park BY, Cho S, Park H, Kim ST, Chung CS (2019) Dynamic functional connectivity of the migraine brain: a resting-state functional magnetic resonance imaging study. Pain 160:2776–2786

    Article  PubMed  Google Scholar 

  29. Wu D, Zhou Y, Xiang J, Tang L, Liu H, Huang S, Wu T, Chen Q, Wang X (2016) Multi-frequency analysis of brain connectivity networks in migraineurs: a magnetoencephalography study. J Headache Pain 17:38

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Yang FC, Chou KH, Hsu AL, Fuh JL, Lirng JF, Kao HW, Lin CP, Wang SJ (2018) Altered brain functional Connectome in migraine with and without restless legs syndrome: a resting-state functional MRI study. Front Neurol 9:25

    Article  PubMed  PubMed Central  Google Scholar 

  31. Headache Classification Committee of the International Headache Society (IHS) (2018) The international classification of headache disorders, 3rd edition. Cephalalgia 38:1–211

    Google Scholar 

  32. Reitz MC, Hrncic D, Treede RD, Caspani O (2016) A comparative behavioural study of mechanical hypersensitivity in 2 pain models in rats and humans. Pain 157:1248–1258

    Article  PubMed  Google Scholar 

  33. Hsiao FJ, Yu HY, Chen WT, Kwan SY, Chen C, Yen DJ, Yiu CH, Shih YH, Lin YY (2015) Increased intrinsic connectivity of the default mode network in temporal lobe epilepsy: evidence from resting-state MEG recordings. PLoS One 10:e0128787

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Hsiao FJ, Wang SJ, Lin YY, Fuh JL, Ko YC, Wang PN, Chen WT (2017) Altered insula-default mode network connectivity in fibromyalgia: a resting-state magnetoencephalographic study. J Headache Pain 18:89

    Article  PubMed  PubMed Central  Google Scholar 

  35. Taulu S, Kajola M, Simola J (2004) Suppression of interference and artifacts by the signal space separation method. Brain Topogr 16:269–275

    Article  PubMed  Google Scholar 

  36. Taulu S, Simola J (2006) Spatiotemporal signal space separation method for rejecting nearby interference in MEG measurements. Phys Med Biol 51:1759–1768

    Article  CAS  PubMed  Google Scholar 

  37. Florin E, Baillet S (2015) The brain's resting-state activity is shaped by synchronized cross-frequency coupling of neural oscillations. Neuroimage 111:26–35

    Article  PubMed  Google Scholar 

  38. Huang MX, Mosher JC, Leahy RM (1999) A sensor-weighted overlapping-sphere head model and exhaustive head model comparison for MEG. Phys Med Biol 44:423–440

    Article  CAS  PubMed  Google Scholar 

  39. Holmes CJ, Hoge R, Collins L, Woods R, Toga AW, Evans AC (1998) Enhancement of MR images using registration for signal averaging. J Comput Assist Tomogr 22:324–333

    Article  CAS  PubMed  Google Scholar 

  40. Klein A, Tourville J (2012) 101 labeled brain images and a consistent human cortical labeling protocol. Front Neurosci 6:171

    Article  PubMed  PubMed Central  Google Scholar 

  41. Babiloni C, Ferri R, Binetti G, Cassarino A, Dal Forno G, Ercolani M, Ferreri F, Frisoni GB, Lanuzza B, Miniussi C, Nobili F, Rodriguez G, Rundo F, Stam CJ, Musha T, Vecchio F, Rossini PM (2006) Fronto-parietal coupling of brain rhythms in mild cognitive impairment: a multicentric EEG study. Brain Res Bull 69:63–73

    Article  PubMed  Google Scholar 

  42. Colclough GL, Woolrich MW, Tewarie PK, Brookes MJ, Quinn AJ, Smith SM (2016) How reliable are MEG resting-state connectivity metrics? Neuroimage 138:284–293

    Article  CAS  PubMed  Google Scholar 

  43. Brookes MJ, Woolrich MW, Barnes GR (2012) Measuring functional connectivity in MEG: a multivariate approach insensitive to linear source leakage. Neuroimage 63:910–920

    Article  CAS  PubMed  Google Scholar 

  44. Tadel F, Baillet S, Mosher JC, Pantazis D, Leahy RM (2011) Brainstorm: a user-friendly application for MEG/EEG analysis. Comput Intell Neurosci 2011:879716

    Article  PubMed  PubMed Central  Google Scholar 

  45. Gross J, Schnitzler A, Timmermann L, Ploner M (2007) Gamma oscillations in human primary somatosensory cortex reflect pain perception. PLoS Biol 5:e133

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Hauck M, Lorenz J, Engel AK (2007) Attention to painful stimulation enhances gamma-band activity and synchronization in human sensorimotor cortex. J Neurosci 27:9270–9277

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Zhang ZG, Hu L, Hung YS, Mouraux A, Iannetti GD (2012) Gamma-band oscillations in the primary somatosensory cortex--a direct and obligatory correlate of subjective pain intensity. J Neurosci 32:7429–7438

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Hu L, Iannetti GD (2019) Neural indicators of perceptual variability of pain across species. Proc Natl Acad Sci U S A 116:1782–1791

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Schulz E, Zherdin A, Tiemann L, Plant C, Ploner M (2012) Decoding an individual's sensitivity to pain from the multivariate analysis of EEG data. Cereb Cortex 22:1118–1123

    Article  PubMed  Google Scholar 

  50. May ES, Nickel MM, Ta Dinh S, Tiemann L, Heitmann H, Voth I, Tolle TR, Gross J, Ploner M (2019) Prefrontal gamma oscillations reflect ongoing pain intensity in chronic back pain patients. Hum Brain Mapp 40:293–305

    Article  PubMed  Google Scholar 

  51. Zhou R, Wang J, Qi W, Liu FY, Yi M, Guo H, Wan Y (2018) Elevated resting state gamma oscillatory activities in electroencephalogram of patients with post-herpetic neuralgia. Front Neurosci 12:750

    Article  PubMed  PubMed Central  Google Scholar 

  52. Tan LL, Oswald MJ, Heinl C, Retana Romero OA, Kaushalya SK, Monyer H, Kuner R (2019) Gamma oscillations in somatosensory cortex recruit prefrontal and descending serotonergic pathways in aversion and nociception. Nat Commun 10:983

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Geuter S, Boll S, Eippert F, Buchel C (2017) Functional dissociation of stimulus intensity encoding and predictive coding of pain in the insula. Elife 6:e24770

  54. Tu Y, Zhang B, Cao J, Wilson G, Zhang Z, Kong J (2019) Identifying inter-individual differences in pain threshold using brain connectome: a test-retest reproducible study. Neuroimage 202:116049

    Article  PubMed  Google Scholar 

  55. Letzen JE, Remeniuk B, Smith MT, Irwin MR, Finan PH, Seminowicz DA (2020) Individual differences in pain sensitivity are associated with cognitive network functional connectivity following one night of experimental sleep disruption. Hum Brain Mapp 41:581–593

    Article  PubMed  Google Scholar 

  56. Garcia-Larrea L, Peyron R (2013) Pain matrices and neuropathic pain matrices: a review. Pain 154(Suppl 1):S29–S43

    Article  PubMed  Google Scholar 

  57. Lorenz J, Minoshima S, Casey KL (2003) Keeping pain out of mind: the role of the dorsolateral prefrontal cortex in pain modulation. Brain 126:1079–1091

    Article  CAS  PubMed  Google Scholar 

  58. Sliwinski C, Nees TA, Puttagunta R, Weidner N, Blesch A (2018) Sensorimotor activity partially ameliorates pain and reduces nociceptive fiber density in the chronically injured spinal cord. J Neurotrauma 35:2222–2238

    Article  PubMed  PubMed Central  Google Scholar 

  59. Senkowski D, Hofle M, Engel AK (2014) Crossmodal shaping of pain: a multisensory approach to nociception. Trends Cogn Sci 18:319–327

    Article  PubMed  Google Scholar 

  60. Pan LH, Wang YF, Lai KL, Chen WT, Chen SP, Ling YH, Chou LW, Treede RD, Wang SJ (2020) Mechanical punctate pain threshold is associated with headache frequency and phase in patients with migraine. In press, Cephalalgia

    Google Scholar 

  61. Schwedt TJ, Chong CD (2014) Correlations between brain cortical thickness and cutaneous pain thresholds are atypical in adults with migraine. PLoS One 9:e99791

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  62. Chong CD, Schwedt TJ, Hougaard A (2019) Brain functional connectivity in headache disorders: a narrative review of MRI investigations. J Cereb Blood Flow Metab 39:650–669

    Article  PubMed  Google Scholar 

  63. Grundstrom H, Larsson B, Arendt-Nielsen L, Gerdle B, Kjolhede P (2019) Associations between pain thresholds for heat, cold and pressure, and pain sensitivity questionnaire scores in healthy women and in women with persistent pelvic pain. Eur J Pain 23:1631–1639

    Article  PubMed  Google Scholar 

  64. Mancini F, Bauleo A, Cole J, Lui F, Porro CA, Haggard P, Iannetti GD (2014) Whole-body mapping of spatial acuity for pain and touch. Ann Neurol 75:917–924

    Article  PubMed  PubMed Central  Google Scholar 

  65. Li F, Xiang J, Wu T, Zhu D, Shi J (2016) Abnormal resting-state brain activity in headache-free migraine patients: a magnetoencephalography study. Clin Neurophysiol 127:2855–2861

    Article  PubMed  Google Scholar 

  66. Liu H, Ge H, Xiang J, Miao A, Tang L, Wu T, Chen Q, Yang L, Wang X (2015) Resting state brain activity in patients with migraine: a magnetoencephalography study. J Headache Pain 16:525

    Article  PubMed  Google Scholar 

  67. Khan S, Gramfort A, Shetty NR, Kitzbichler MG, Ganesan S, Moran JM, Lee SM, Gabrieli JD, Tager-Flusberg HB, Joseph RM, Herbert MR, Hamalainen MS, Kenet T (2013) Local and long-range functional connectivity is reduced in concert in autism spectrum disorders. Proc Natl Acad Sci U S A 110:3107–3112

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Tang W, Liu H, Douw L, Kramer MA, Eden UT, Hamalainen MS, Stufflebeam SM (2017) Dynamic connectivity modulates local activity in the core regions of the default-mode network. Proc Natl Acad Sci U S A 114:9713–9718

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Tadel F, Bock E, Niso G, Mosher JC, Cousineau M, Pantazis D, Leahy RM, Baillet S (2019) MEG/EEG group analysis with brainstorm. Front Neurosci 13:76

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank the study participants for actively participating.

Funding

The authors disclose receipt of the following financial support for the research, authorship, and article publication: Ministry of Science and Technology of Taiwan (MOST 107–2314-B-075-015-MY2–2 to WT Chen, 108–2321-B-010-001 and 108–2321-B-010-014-MY2 to SJ Wang, and 107–2221-E-010-007 and 108–2221-E-010-004 to FJ Hsiao), Yen Tjing Ling Medical Foundation (CI-109-1 to FJ Hsiao), and Taipei-Veterans General Hospital (V108C-129 and V107C-091 to WT Chen). This work was supported by the Brain Research Center, National Yang-Ming University from The Featured Areas Research Center Program within the framework of the Higher Education Sprout Project by the Ministry of Education in Taiwan. The funders had no role in the study design, data collection and analysis, decision to publish, or manuscript preparation.

Author information

Authors and Affiliations

Authors

Contributions

FJ Hsiao, WT Chen, and SJ Wang conceived and design the work and wrote the article. FJ Hsiao, WT Chen, HY Liu, YF Wang, SP Chen, KL Lai, and SJ Wang acquired the data. FJ Hsiao, WT Chen, HY Liu, YF Wang, SP Chen, KL Lai, LLH Pan, SJ Wang analyzed the data and participated in the discussion and provided the comments. The authors have read and approved the manuscript.

Corresponding authors

Correspondence to Fu-Jung Hsiao or Shuu-Jiun Wang.

Ethics declarations

Ethics approval and consent to participate

All participants completed informed consent forms after receiving a complete explanation of the study. The Institutional Review Board of Taipei Veterans General Hospital approved the study protocol (VGHTPE-IRB, 2015–10-001 BC).

Consent for publication

Not applicable.

Competing interests

FJ Hsiao, WT Chen, HY Liu, YF Wang, SP Chen, KL Lai, and LLH Pan declare no potential conflicts of interest. SJ Wang has served on the advisory boards of Eli Lilly, Daiichi-Sankyo, Taiwan Pfizer and Taiwan Norvatis. He has received honoraria as a moderator from Allergan, Pfizer, Eli Lilly, Bayer, and Eisai.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hsiao, FJ., Chen, WT., Liu, HY. et al. Individual pain sensitivity is associated with resting-state cortical activities in healthy individuals but not in patients with migraine: a magnetoencephalography study. J Headache Pain 21, 133 (2020). https://doi.org/10.1186/s10194-020-01200-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s10194-020-01200-8

Keywords