Skip to main content

Increased connectivity of pain matrix in chronic migraine: a resting-state functional MRI study



To investigate the whole-brain resting-state functional connectivity in patients with chronic migraine (CM) using a data-driven method.


We prospectively recruited patients with either episodic migraine (EM) or CM aged 18–60 years who visited the headache clinic of the Samsung Medical Center from July 2016 to December 2017. All patients underwent 3 T MRI using an identical scanner. Patients were considered interictal if they did not have a migraine headache at the day and ± 1 days of functional MRI acquisition. Using the group-independent component analysis (ICA), connectivity analysis with a weighted and undirected network model was performed. The between-group differences in degree centrality (DC) values were assessed using 5000 permutation tests corrected with false discovery rate (FDR).


A total of 62 patients (44 EM and 18 CM) were enrolled in this study. Among the seven functionally interpretable spatially independent components (ICs) identified, only one IC, interpreted as the pain matrix, showed a significant between-group difference in DC (CM > EM, p = 0.046). This association remained significant after adjustment for age, sex, migraine with aura (MWA), allodynia, depression, and anxiety (p = 0.038). The pain matrix was functionally correlated with the hypothalamus (p = 0.040, EM > CM) and dorsal raphe nucleus (p = 0.039, CM > EM) with different levels of strength in EM and CM.


CM patients have a stronger connectivity in the pain matrix than do EM patients. Functional alteration of the pain network might play a role in migraine chronification.


Migraine is a neurological disorder characterized by episodic headaches associated with nausea, vomiting, and increased sensitivity to external stimuli. Chronic migraine (CM) is a devastating subtype of migraine, which is defined as headache days of 15 or more per month and migraine days of eight or more per month for > 3 months [1]. CM has an estimated prevalence of 1.5% worldwide [2]. CM is more disabling and results in a much higher disease burden than episodic migraine (EM) [3].

About 3% of EMs progress to CMs annually [4]. However, the pathophysiology of migraine chronification is still unknown. Although epidemiological studies revealed risk factors that promote the conversion of EM to CM, biological mechanism of migraine chronification has not been fully elucidated yet, particularly in the absence of medication overuse. Patients can develop CM with or without triggers such as stressful life events, weight gain, and caffeine overuse, which are also common in the lives of healthy people [4, 5]. Therefore, it is likely that a predisposition to migraine chronification exists. To elucidate this, researchers have investigated functional features of the CM brain [6, 7].

To date, studies on functional neuroimaging features of CM are relatively scarce and focused only on predefined brain areas. Schwedt et al. revealed that affective pain regions (anterior insula, amygdala, and PAG) are functionally connected to other brain regions differently in CM patients and normal controls [8]. In studies involving experimental fMRI, CM patients have an enhanced activation of brain regions such as the anterior hypothalamus in response to nociceptive stimuli and nociceptive trigeminal nucleus in response to visual stimuli, compared to normal controls [6, 7]. These studies suggest that distinct functional characteristics of CM exist. However, no study has investigated the whole-brain functional features of CM in comparison to those of EM.

In this study, we aimed to investigate the whole brain resting-state functional connectivity in patients with CM compared to those with EM using a data-driven method. A thorough clinical evaluation and a strict correction for multiple comparisons were performed.


Study subjects

We prospectively recruited new migraine patients who visited the headache clinic of the Samsung Medical Center from July 2016 to December 2017. Patients who were 1) aged 18–60 years, 2) diagnosed with EM or CM, and 3) currently not on migraine preventive medications were included in the study. Patients were excluded if they had 1) medication-overuse headache, 2) chronic pain disorders other than migraine, 3) an alleged diagnosis of major psychiatric disorders such as bipolar affective disorder and schizophrenia, or 4) were currently undergoing treatment for depression or anxiety. Migraine diagnosis was made by two headache specialists (M.J.L. and C.S.C.) according to the International Classification of Headache Disorders, 3rd edition beta version [9].

Standard protocol approvals, registrations, and patient consents

The Samsung Medical Center Institutional Review Board approved this study. All patients and controls provided written informed consent prior to participation.

Clinical evaluation

All patients completed a structured questionnaire designed to characterize their headaches. Subsequently, an investigator (M.J.L.) interviewed all patients to verify their responses on the questionnaires. The presence of allodynia was confirmed during the interview. Patients also completed the Allodynia Severity Checklist-12, Patient Health Questionnaire-9 (PHQ-9), and Hospital Anxiety and Depression Scale (HADS) [10,11,12].

Functional magnetic resonance image acquisition

All the study subjects underwent magnetic resonance (MR) imaging using a 3 T MR scanner (Achieva; Philips Medical Systems, Best, the Netherlands). T1-weighted structural images with the following imaging parameters were obtained: repetition time (TR) = 9.87 ms; echo time (TE) = 4.59 ms; field of view (FOV) = 240 × 240 mm2; number of slices = 360; slice thickness = 0.5 mm; and pixel resolution = 0.5 mm2. The imaging parameters for resting-state functional magnetic resonance image (rs-fMRI) are as follows: TR = 3000 ms; TE = 35 ms; flip angle = 90o; FOV = 220 × 220 mm2; number of slices = 35; slice thickness = 4 mm; pixel resolution = 1.7 mm2; and number of volumes = 100.

For all included patients, we assessed the presence and characteristics of headache and the use of acute medications at the day and ± 1 days of fMRI acquisition. Patients were considered interictal if they did not have migraine headache, defined as any headache of moderate to severe intensity, headaches with nausea, vomiting, photophobia, or phonophobia, or headaches that led to the taking of acute migraine medications, at the day and ± 1 days of fMRI acquisition.

Data preprocessing

Imaging data were preprocessed using fusion of neuroimaging preprocessing (FuNP) pipeline that integrated the AFNI and FSL software [13]. T1-weighted structural MR images were processed by correcting the magnetic field inhomogeneity and removing non-brain tissues. The rs-fMRI data were also processed. The volumes of data obtained during the first 12 s (i.e., four volumes) were discarded to allow the magnetic field to be saturated. The frame-wise displacement (FD) between time series volumes was calculated and the volumes with FD exceeding 0.5 mm were removed [14]. Head motion correction was performed on the remaining time series volumes. Slice timing correction was performed and intensity normalization with a mean value of 10,000 was applied to all the volumes. Nuisance variables such as contributions from white matter, cerebrospinal fluid, head motion, heart, breathing, and the large vein were removed using the FIX software [15]. The low-resolution fMRI data were registered onto the high-resolution T1-weighted data and subsequently onto the Montreal Neurological Institute (MNI) standard space. A band-pass filter with frequency between 0.009 and 0.08 Hz and spatial smoothing with a full width at half maximum of 6 mm was applied.

Group ICA

The preprocessed rs-fMRI data of all subjects were temporally concatenated, and group independent component analysis (ICA) was performed to automatically generate spatially independent components (ICs) using the FSL MELODIC software [16]. The generated ICs were classified into signal and noise components with two criteria. First, the cross-correlation between the generated ICs and known resting state networks was calculated, and ICs with correlation values less than 0.25 were considered as noise components [17]. Second, signal and noise components were classified by visual inspection based on their spatial map, time series, and frequency spectrum [18, 19]. To identify the functional characteristics of each IC, we performed ‘cognitive decoding’ using Neurosynth software ( [20]. Neurosynth is an open-source software platform for meta-analyses that enables us to identify relevant specific terms relevant to given activation maps by searching large-scale studies (3228 terms in 14,371 studies, as of March 2019). We applied the cognitive decoding process to the z-statistic map of each IC to identify relevant terms. The cognitive decoding process resulted in correlation values between the z-statistic map of each IC and the activation map of specific terms.

Main analysis: whole-brain functional connectivity analysis

Connectivity analysis with a weighted and undirected network model was performed. Graph nodes were defined using functionally interpretable ICs (i.e. node = IC) and graph edges were defined as the partial correlation with L2-norm between the time series of different nodes [21]. The correlation values were soft-thresholded to avoid binarizing edge weights [22, 23]. The soft-thresholded correlation values were transformed to z-values using Fisher’s r-to-z transformation. Degree centrality (DC), which measures the importance of a given node, was calculated for each node by summing all edge weights connected to a given node [24]. DC values of each node were used for identifying differences between the EM and CM groups.

Secondary analysis: connectivity with the hypothalamus, dorsal raphe nucleus, and periaqueductal gray

After identifying brain networks which showed significant between-group differences, we tested whether the identified network was relevant to migraine pathophysiology. We defined three regions of interest (ROIs): the hypothalamus, dorsal raphe nucleus (DRN), and periaqueductal gray (PAG). The hypothalamus and PAG were manually drawn, while the DRN was defined using the Harvard ascending arousal network atlas via image co-registration (Fig. 1). [25] The centroid coordinates of ROIs in the MNI standard space were consistent with previous studies (hypothalamus: x = 0, y = − 4, z = − 9; PAG: x = 1, y = − 31, z = − 9; DRN: x = 1, y = − 32, z = − 17) [26,27,28,29]. The time series were extracted from each ROI and their correlation with the time series of the identified functional network in the main analysis were computed between all possible pairs. The correlation values were transformed to z-values using Fisher’s r-to-z transformation. The z-transformed correlation values were used for identifying differences between the EM and CM groups.

Fig. 1
figure 1

Region of interest segmentation results. The hypothalamus, dorsal raphe nuclei (DRN), and periaqueductal gray (PAG) were segmented on a three-dimensional brain atlas

Statistical analysis

Clinical variables were compared between EM and CM groups using the Chi-square test, Fisher’s exact test, or Mann-Whitney tests. The differences in DC values between EM and CM groups were assessed using permutation tests followed by false discovery rate (FDR) correction [30, 31]. Subjects were randomly assigned to the EM and CM groups 5000 times, and a null distribution was constructed. The ICs with DC values outside 95% of the null distribution were considered significant ICs with significant between-group differences. The p-values were further corrected using FDR (p < 0.05, corrected) [30]. Multivariable linear regression analysis with adjustment for age, sex, presence of aura, allodynia, depression (PHQ-9 scores ≥8), anxiety (HADS-A scores of ≥8), disease duration, headache intensity, and acute antimigraine drug use/month was performed. Pearson’s correlation analysis between the strength of identified brain networks and monthly headache days was performed to assess if the group difference was the consequence of frequent headaches. The correlation analysis was also performed between the strength of identified brain networks and clinical variables such as patients’ disease duration, allodynia, anxiety, and depression scores. Interaction analysis was performed to determine a possible modifying effect of the presence of depression and anxiety. Statistical analysis was performed using MATLAB 2017a (Mathworks Inc., Natick, MA, USA) and SPSS software (IBM-SPSS. Chicago, IL, USA).


Study subjects

Among 86 eligible patients, 64 (45 EM and 19 CM) underwent the interictal study. Among them, two (1 EM and 1 CM) were excluded from the analysis because of poor image quality. Finally, data from 62 patients (44 EM and 18 CM) were used for the analysis. Demographics and characteristics of patients are summarized in Table 1.

Table 1 Demographics and characteristics of study participants

Functional network identification

The group-ICA approach automatically generated nine ICs (Fig. 2). Two ICs (white matter and noise components) were excluded from further analyses. Finally, seven functionally interpretable ICs were identified. All but IC 3 were compatible with known resting-state functional networks: IC 1 (visual network) comprises the lingual gyrus, and superior- and inferior- occipital cortices; IC 2 (default mode network) comprises the posterior cingulate cortex and precuneus; IC 4 (executive control network) comprises the medial prefrontal cortex, orbitofrontal cortex, and anterior cingulate cortex (ACC); IC 5 (frontoparietal network) comprises the superior frontal cortex, angular gyrus, and posterior cingulate cortex; IC 6 (frontoparietal network) comprises the orbitofrontal cortex, and superior and inferior parietal gyri; and IC 7 (sensorimotor network) comprises the pre-central and post-central gyri and paracentral lobule. Major components of IC 3 included the dorsolateral prefrontal cortex (DLPFC), anterior insula, ACC, thalamus, and precuneus (Fig. 3). IC 3 also included supramarginal gyrus, planum temporale, premotor cortex, and cerebellum. Based on the cognitive decoding process using Neurosynth, we identified IC 3 as the pain matrix (Table 2).

Fig. 2
figure 2

Resting-state networks identified using independent component analysis. Nine automatically generated independent components (ICs). ICs 8 and 9 were considered as noise components and were therefore excluded. ICs 1 to 7 are functionally interpretable ICs

Fig. 3
figure 3

Pain matrix. The functional network (IC 3) included the dorsolateral prefrontal cortex, anterior insular cortex, anterior cingulate cortex, thalamus, and precuneus, suggesting a pain matrix

Table 2 The results of the cognitive decoding process of the ICs using Neurosynth software

Differences in functional connectivity

Figure 4 shows between-group differences in DC values among identified ICs. A significant between-group difference was found only in IC 3 (pain matrix). Patients with CM showed stronger connectivity in terms of DC in the pain matrix than those with EM (uncorrected p = 0.0066 and FDR-corrected p = 0.0462). This between-group difference remained significant after adjustment for covariates such as age, sex, migraine with aura (MWA), allodynia, depression, anxiety, disease duration, headache intensity, and acute antimigraine drug use/month (Table 3).

Fig. 4
figure 4

Group comparison among identified resting-state networks. DC values were compared between EM and CM by using permutation tests with FDR. Each bar and error bar represent the mean and standard errors of mean, respectively

Table 3 Univariable and multivariable analysis of between-group difference in the connectivity of the pain matrix

To investigate the relative importance of the sub-regions of the pain matrix, we calculated mean z-statistic values of the ICA weights from the sub-regions within the pain matrix. The sub-region with the highest z-statistic value was ACC (= 3.693) followed by precuneus (= 3.650), DLPFC (= 3.548), premotor cortex (= 3.004), supramarginal gyrus (= 2.956), planum temporale (= 2.937), cerebellum (= 2.825), anterior insula (= 2.789), and thalamus (= 2.419). The results showed that ACC was the most important sub-region in the IC 3 and might be the largest contributing factor to explain the between-group differences between EM and CM groups.

Clinical correlates of the pain matrix connectivity

We performed correlation analyses between clinical variables and the connectivity (i.e. DC values) of pain matrix. No significant correlation was found between the pain matrix connectivity and headache days (r = 0.0444, p = 0.7321), HADS-D score (r = 0.1080, p = 0.4239), HADS-A score (r = 0.0322, p = 0.8119), PHQ-9 score (r = − 0.0400, p = 0.7638), ASC-12 score (r = 0.0306, p = 0.8136), and disease duration (r = − 0.0910, p = 0.4818).

Neither depression nor anxiety modified the association between CM and the pain matrix connectivity (P for interaction = 0.479 and 0.425, respectively). The presence of mild non-migrainous headache on the day of fMRI acquisition also did not modify this association (P for interaction = 0.372).

Functional correlates of the pain matrix connectivity

Figure 5 shows connectivity between key regions involved in migraine pathophysiology and the pain matrix. The strength of the functional connection between the pain matrix and the hypothalamus (CM > EM, FDR-corrected p = 0.0399) and DRN (EM > CM, FDR-corrected p = 0.0390) was different between groups. No significant between-group difference was found in the connectivity between the pain matrix and PAG (FDR-corrected p = 0.2738).

Fig. 5
figure 5

Functional connectivity of IC 3 (pain matrix) with the hypothalamus, DRN, and PAG. a Functional connectivity of IC 3 with the hypothalamus, dorsal raphe nuclei (DRN), and periaqueductal gray (PAG) are illustrated. The mean (SD) of edge values are shown and the width of the lines indicate the magnitude of edge weights. b Between-group analysis of connectivity between IC 3 and each region are summarized. Patients with CM showed a stronger connectivity between IC 3 and the hypothalamus (FDR-corrected p = 0.0399), whereas the connectivity between the DRN and IC 3 was stronger in patients with EM (FDR-corrected p = 0.0390). No between-group difference in the connectivity between IC 3 and the PAG were noted (FDR-corrected p = 0.2738)


In this study, we found that 1) the connectivity in the pain matrix differed between EM and CM patients; 2) the pain matrix connectivity was not correlated with headache frequency or psychiatric comorbidities; and 3) the strength of the functional connection between the pain matrix and the hypothalamus and DRN was different between EM and CM groups. An enhanced connectivity of the pain matrix may play a role in migraine chronification.

Data-driven vs. ROI-based method

To date, the neural mechanism underlying migraine chronification is still unknown. To unveil functional characteristics of CM, functional neuroimaging is used for research. Using a resting-state or task-specific functional MRI, specific brain regions were tested with a priori hypothesis using ROI-based approaches [6,7,8]. However, no study has compared whole-brain connectivity features between CM and EM. In this study, the group ICA approach was adopted to define large-scale brain networks. The major advantage of using group ICA over pre-defined atlases is that it is a data-driven approach. There are many existing pre-defined atlases including Brodmann areas, automated anatomical labeling (AAL), and Harvard-Oxford atlases [32, 33]. However, the atlases were constructed using different pools of subjects and each atlas provides a different number of brain regions, which might reduce the reproducibility of the neuroimaging studies. In addition, the pre-defined atlases might not reflect the functional characteristics adequately since they were derived from a different set of subjects. The group-ICA approach yields brain networks that share similar activity patterns among the patients and thus, reflects the functional characteristics of the data more robustly than the atlas-based approach.

Pain matrix in CM

In this study, a functional network (pain matrix), which comprised the ACC, anterior insulae, thalami, DLPFC, precuneus, supramarginal gyri, and cerebellum, differed between EM and CM in its functional connectivity. The concept of a pain matrix has been challenged by studies which reported that similar areas are activated in response to non-nociceptive stimuli [34]. In addition, the pain matrix overlaps with the salience circuit, which is implicated in chronic pain processing [35]. However, in our study, the pain matrix additionally involves areas such as the DLPFC and supramarginal gyri, which are considered as major components of the central executive network. In addition to the salience circuit, these regions overlap with areas involved in pain experience (the ACC, anterior insular, and thalamus) [36], cognitive modulation of pain sensitivity (precuneus) [37], pain expectation (DLPFC, insula, ACC, globus pallidus, putamen, thalamus, and cerebellum) [38], and pain catastrophizing on mild pain (ACC, insula, DLPFC, precuneus, thalamus, putamen, inferior parietal lobule, parahippocampal gyrus) [39]. In concordance with the results of a previous study using an experimental fMRI paradigm which reported that migraineurs have enhanced pain-induced activity of the pain matrix [40], we observed that the functional connectivity of pain matrix was more greater in patients with CM.

The matrix identified in this study is different from the functional MRI markers of acute pain which include somatosensory areas (S1 and S2) and PAG [41]. CM can be either a predisposition to or state of frequent headaches. However, CM brains did not show markers of continuing acute pain. Instead, the insulae and ACC, which play a major role in chronic pain [42], as well as other pain-related and cognitive areas, had a stronger functional connectivity in CM. Our data suggest that a stronger connectivity of the pain matrix is a characteristic of the CM brain, which might play a major role in migraine chronification.

Different functional features between CM and EM

Whether EM and CM are different disorders or in a single continuum has been debated for a long time. In this resting-state fMRI study, we suggest that CM has functional characteristics distinct from EM. Previously, only a limited number of studies investigated resting-state fMRI features of CM and most of them focused on specific structures of interest (e.g. amygdala, insula, and ACC) [8, 43]. Our findings are in line with previous study results on involvement of limbic structures in CM, especially ACC and insular cortex. Although our cross-sectional study is not suitable for proving any causal relationship, the functional connectivity did not correlate with headache frequencies or psychologic comorbidities, suggesting that our findings are not a consequence of them but a predisposition to migraine chronification. We are currently conducting a prospective fMRI study to test the change of brain functional characteristics in association with disease courses in patients with migraine ( Identifier: NCT03487978).

Comparison with functional neuroimaging of chronic pain

fMRI has been used in research on chronic pain disorders in several previous studies. The prefrontal cortices, insulae, and ACCs are reported to be activated in most chronic pain disorders. However, conflicting results on the resting-state connectivity of networks comprising the aforementioned areas exist. Specifically, a greater connectivity was reported between the default mode network and ACC in patients with diabetic neuropathic pain [44]; between the default mode network and insular cortices in patients with fibromyalgia [45]; in the salience network, central executive network, and default mode network in pediatric patients with complex regional pain syndrome (CRPS) [46]; and between the right insular cortices and cingulate gyri in patients with fibromyalgia [47]; while reduced connectivity among the medial prefrontal cortex, insular cortex, and ACC were found in patients with chronic pain disorders such as CRPS, knee osteoarthritis, and chronic back pain [48].

Clinically, CM differs from other chronic pain disorders because the “chronicity” in migraine does not imply persistent pain but increased days of headache, which is a combination of frequent episodic attacks (a function of the headache generator) and lower tendency toward clear remission (a function of either an enhanced pain signaling or decreased pain modulation). Our main analysis showed that an increased connectivity within the pain matrix may play a role in migraine chronification. However, the role of the migraine generator should not be overlooked because small structures such as the hypothalamus and PAG might have been missed by the whole-brain, data-driven approach used in our study. Thus, we further tested three structures involved in migraine pathophysiology: the hypothalamus, a migraine generator; PAG, a pain modulatory center; and DRN, the main serotonergic center activated during migraine attack [49,50,51]. As a result, patients with CM showed an increased connectivity between the pain matrix and hypothalamus compared to those with EM, while the connectivity between the pain matrix and DRN was weaker in CM patients. Taken together with a recent report of increased hypothalamic activation to painful stimuli in CM [6], we suggest that the hypothalamus is more easily activated by external stimuli and strongly connected to the pain matrix in patients with CM, while the brain connectivity between the pain matrix and serotonergic system in patients with CM is relatively weak. Our data suggest that the pain matrix is another key player in CM pathophysiology.

Our study has several strengths. First, we carefully defined CM and performed interictal fMRI imaging in patients with CM. Second, a data-driven method was used, and a strict statistical correction was performed to avoid pitfalls during multiple comparisons. The limitations of our study include 1) the small number of participants, 2) the lack of normal controls, and 3) the cross-sectional nature of the study. To overcome these limitations, we are currently conducting a longitudinal fMRI study in migraineurs and controls to study the effects of the disease and time on fMRI changes. In addition, we could not investigate some important structures (e.g. dorsal pons) and nuclei involved in migraine pathophysiology (e.g., the nucleus raphe magnus and superior salivatory nucleus) because of non-availability of reliable atlases.


In conclusion, CM has an enhanced functional connectivity of the pain matrix which has a different functional connection to hypothalamus and DRN compared to EM patients. Functional alteration of the pain network might play a role in migraine chronification.



Automated anatomical labeling


Chronic migraine


Complex regional pain syndrome


Degree centrality


Dorsolateral prefrontal cortex


Dorsal raphe nucleus


Episodic migraine


Frame-wise displacement


False discovery rate


Field of view


Hospital Anxiety and Depression Scale


Independent component


Independent component analysis


Montreal Neurological Institute


Magnetic resonance


Migraine with aura


Periaqueductal gray


Patient Health Questionnaire-9


Region of interest


Resting-state functional magnetic resonance image


Echo time


Repetition time


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

    Google Scholar 

  2. Natoli JL, Manack A, Dean B, Butler Q, Turkel CC, Stovner L et al (2010) Global prevalence of chronic migraine: a systematic review. Cephalalgia 30:599–609

    CAS  Article  Google Scholar 

  3. Lanteri-Minet M, Duru G, Mudge M, Cottrell S (2011) Quality of life impairment, disability and economic burden associated with chronic daily headache, focusing on chronic migraine with or without medication overuse: a systematic review. Cephalalgia 31:837–850

    Article  Google Scholar 

  4. Scher AI, Stewart WF, Ricci JA, Lipton RB (2003) Factors associated with the onset and remission of chronic daily headache in a population-based study. Pain 106:81–89

    CAS  Article  Google Scholar 

  5. Scher AI, Stewart WF, Lipton RB (2004) Caffeine as a risk factor for chronic daily headache: a population-based study. Neurology 63:2022–2027

    CAS  Article  Google Scholar 

  6. Schulte LH, Allers A, May A (2017) Hypothalamus as a mediator of chronic migraine: evidence from high-resolution fMRI. Neurology 88:2011–2016

    Article  Google Scholar 

  7. Schulte LH, Allers A, May A (2018) Visual stimulation leads to activation of the nociceptive trigeminal nucleus in chronic migraine. Neurology 90:e1973–e1978

    Article  Google Scholar 

  8. Schwedt TJ, Schlaggar BL, Mar S, Nolan T, Coalson RS, Nardos B et al (2013) Atypical resting-state functional connectivity of affective pain regions in chronic migraine. Headache 53:737–751

    Article  Google Scholar 

  9. Headache Classification Committee of the International Headache S (2013) The international classification of headache disorders, 3rd edition (beta version). Cephalalgia 33:629–808

    Article  Google Scholar 

  10. Lipton RB, Bigal ME, Ashina S, Burstein R, Silberstein S, Reed ML et al (2008) Cutaneous allodynia in the migraine population. Ann Neurol 63:148–158

    Article  Google Scholar 

  11. Seo JG, Park SP (2015) Validation of the patient health Questionnaire-9 (PHQ-9) and PHQ-2 in patients with migraine. J Headache Pain 16:65

    Article  Google Scholar 

  12. Mykletun A, Stordal E, Dahl AA (2001) Hospital anxiety and depression (HAD) scale: factor structure, item analyses and internal consistency in a large population. Br J Psychiatry 179:540–544

    CAS  Article  Google Scholar 

  13. Park BY, Byeon K, Park H (2019) FuNP (fusion of neuroimaging preprocessing) pipelines: a fully automated preprocessing software for functional magnetic resonance imaging. Front Neuroinform 13:5

    Article  Google Scholar 

  14. Power JD, Barnes KA, Snyder AZ, Schlaggar BL, Petersen SE (2012) Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59:2142–2154

    Article  Google Scholar 

  15. Salimi-Khorshidi G, Douaud G, Beckmann CF, Glasser MF, Griffanti L, Smith SM (2014) Automatic denoising of functional MRI data: combining independent component analysis and hierarchical fusion of classifiers. Neuroimage 90:449–468

    Article  Google Scholar 

  16. Beckmann CF, DeLuca M, Devlin JT, Smith SM (2005) Investigations into resting-state connectivity using independent component analysis. Philos Trans R Soc Lond Ser B Biol Sci 360:1001–1013

    Article  Google Scholar 

  17. Smith SM, Fox PT, Miller KL, Glahn DC, Fox PM, Mackay CE et al (2009) Correspondence of the brain's functional architecture during activation and rest. Proc Natl Acad Sci U S A 106:13040–13045

    CAS  Article  Google Scholar 

  18. Griffanti L, Douaud G, Bijsterbosch J, Evangelisti S, Alfaro-Almagro F, Glasser MF et al (2017) Hand classification of fMRI ICA noise components. Neuroimage 154:188–205

    Article  Google Scholar 

  19. Kelly RE Jr, Alexopoulos GS, Wang Z, Gunning FM, Murphy CF, Morimoto SS et al (2010) Visual inspection of independent components: defining a procedure for artifact removal from fMRI data. J Neurosci Methods 189:233–245

    Article  Google Scholar 

  20. Yarkoni T, Poldrack RA, Nichols TE, Van Essen DC, Wager TD (2011) Large-scale automated synthesis of human functional neuroimaging data. Nat Methods 8:665–670

    CAS  Article  Google Scholar 

  21. Smith SM, Miller KL, Salimi-Khorshidi G, Webster M, Beckmann CF, Nichols TE et al (2011) Network modelling methods for FMRI. Neuroimage 54:875–891

    Article  Google Scholar 

  22. Mumford JA, Horvath S, Oldham MC, Langfelder P, Geschwind DH, Poldrack RA (2010) Detecting network modules in fMRI time series: a weighted network analysis approach. Neuroimage 52:1465–1476

    Article  Google Scholar 

  23. Schwarz AJ, McGonigle J (2011) Negative edges and soft thresholding in complex network analysis of resting state functional connectivity data. Neuroimage 55:1132–1146

    Article  Google Scholar 

  24. Rubinov M, Sporns O (2010) Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52:1059–1069

    Article  Google Scholar 

  25. Edlow BL, Takahashi E, Wu O, Benner T, Dai G, Bu L et al (2012) Neuroanatomic connectivity of the human ascending arousal system critical to consciousness and its disorders. J Neuropathol Exp Neurol 71:531–546

    Article  Google Scholar 

  26. Denuelle M, Fabre N, Payoux P, Chollet F, Geraud G (2007) Hypothalamic activation in spontaneous migraine attacks. Headache 47:1418–1426

    PubMed  Google Scholar 

  27. Linnman C, Moulton EA, Barmettler G, Becerra L, Borsook D (2012) Neuroimaging of the periaqueductal gray: state of the field. Neuroimage 60:505–522

    Article  Google Scholar 

  28. Ulrich M, Keller J, Gron G (2016) Dorsal raphe nucleus Down-regulates medial prefrontal cortex during experience of flow. Front Behav Neurosci 10:169

    Article  Google Scholar 

  29. Beliveau V, Svarer C, Frokjaer VG, Knudsen GM, Greve DN, Fisher PM (2015) Functional connectivity of the dorsal and median raphe nuclei at rest. Neuroimage 116:187–195

    Article  Google Scholar 

  30. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol:289–300

  31. Chen C, Witte M, Heemsbergen W, van Herk M (2013) Multiple comparisons permutation test for image based data mining in radiotherapy. Radiat Oncol 8:293

    Article  Google Scholar 

  32. Brodmann K, Vergleichende Lokalisationslehre der Grosshirnrinde in ihren Prinzipien dargestellt auf Grund des Zellenbaues. 1909: Barth

  33. Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N et al (2002) Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15:273–289

    CAS  Article  Google Scholar 

  34. Legrain V, Iannetti GD, Plaghki L, Mouraux A (2011) The pain matrix reloaded: a salience detection system for the body. Prog Neurobiol 93:111–124

    Article  Google Scholar 

  35. Borsook D, Edwards R, Elman I, Becerra L, Levine J (2013) Pain and analgesia: the value of salience circuits. Prog Neurobiol 104:93–105

    Article  Google Scholar 

  36. Bantick SJ, Wise RG, Ploghaus A, Clare S, Smith SM, Tracey I (2002) Imaging how attention modulates pain in humans using functional MRI. Brain 125:310–319

    Article  Google Scholar 

  37. Goffaux P, Girard-Tremblay L, Marchand S, Daigle K, Whittingstall K (2014) Individual differences in pain sensitivity vary as a function of precuneus reactivity. Brain Topogr 27:366–374

    Article  Google Scholar 

  38. Koyama T, McHaffie JG, Laurienti PJ, Coghill RC (2005) The subjective experience of pain: where expectations become reality. Proc Natl Acad Sci U S A 102:12950–12955

    CAS  Article  Google Scholar 

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

    Article  Google Scholar 

  40. Schwedt TJ, Chong CD, Chiang CC, Baxter L, Schlaggar BL, Dodick DW (2014) Enhanced pain-induced activity of pain-processing regions in a case-control study of episodic migraine. Cephalalgia 34:947–958

    Article  Google Scholar 

  41. Wager TD, Atlas LY, Lindquist MA, Roy M, Woo CW, Kross E (2013) An fMRI-based neurologic signature of physical pain. N Engl J Med 368:1388–1397

    CAS  Article  Google Scholar 

  42. Malinen S, Vartiainen N, Hlushchuk Y, Koskinen M, Ramkumar P, Forss N et al (2010) Aberrant temporal and spatial brain activity during rest in patients with chronic pain. Proc Natl Acad Sci U S A 107:6493–6497

    CAS  Article  Google Scholar 

  43. Chen Z, Chen X, Liu M, Dong Z, Ma L, Yu S (2017) Altered functional connectivity of amygdala underlying the neuromechanism of migraine pathogenesis. J Headache Pain 18:7

    Article  Google Scholar 

  44. Cauda F, Sacco K, Duca S, Cocito D, D'Agata F, Geminiani GC et al (2009) Altered resting state in diabetic neuropathic pain. PLoS One 4:e4542

    Article  Google Scholar 

  45. Napadow V, LaCount L, Park K, As-Sanie S, Clauw DJ, Harris RE (2010) Intrinsic brain connectivity in fibromyalgia is associated with chronic pain intensity. Arthritis Rheum 62:2545–2555

    Article  Google Scholar 

  46. Becerra L, Sava S, Simons LE, Drosos AM, Sethna N, Berde C et al (2014) Intrinsic brain networks normalize with treatment in pediatric complex regional pain syndrome. Neuroimage Clin 6:347–369

    Article  Google Scholar 

  47. Ichesco E, Schmidt-Wilcke T, Bhavsar R, Clauw DJ, Peltier SJ, Kim J et al (2014) Altered resting state connectivity of the insular cortex in individuals with fibromyalgia. J Pain 15:815–826 e1

    CAS  Article  Google Scholar 

  48. Baliki MN, Mansour AR, Baria AT, Apkarian AV (2014) Functional reorganization of the default mode network across chronic pain conditions. PLoS One 9:e106133

    Article  Google Scholar 

  49. Hamel E (2007) Serotonin and migraine: biology and clinical implications. Cephalalgia 27:1293–1300

    CAS  Article  Google Scholar 

  50. Weiller C, May A, Limmroth V, Juptner M, Kaube H, Schayck RV et al (1995) Brain stem activation in spontaneous human migraine attacks. Nat Med 1:658–660

    CAS  Article  Google Scholar 

  51. Goadsby PJ, Holland PR, Martins-Oliveira M, Hoffmann J, Schankin C, Akerman S (2017) Pathophysiology of migraine: a disorder of sensory processing. Physiol Rev 97:553–622

    Article  Google Scholar 

Download references


We thank Ms. Miran Jung and Ms. Jinju Lee for aiding in data management.


This study was supported by the National Research Foundation of Korea (NRF) grants funded by the Korean government (MSIP) (Nos. 2017R1A2B2009086 and 2017R1A2B4007254). The SK chemical supported the data management. No financial arrangement has a possible conflict of interest.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations



MJL and C-SC designed the study, MJL, SC, STK, and C-SC acquired the data, MJL, B-yP, and HP analyzed the data, all authors interpreted the data. MJL and B-yP drafted the manuscript, SC, STK, HP, and C-SC revised it crtically for important intellectual content. All authors approved the final version of manuscript. All authors agreed to be accountable for all aspects of the work.

Corresponding author

Correspondence to Chin-Sang Chung.

Ethics declarations

Ethics approval and consent to participate

The Samsung Medical Center Institutional Review Board approved this study. All patients and controls provided written informed consent prior to participation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lee, M.J., Park, By., Cho, S. et al. Increased connectivity of pain matrix in chronic migraine: a resting-state functional MRI study. J Headache Pain 20, 29 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Migraine
  • Chronic migraine
  • Functional MRI
  • Neuroimaging
  • Pain matrix