Skip to main content

Altered effective connectivity of posterior thalamus in migraine with cutaneous allodynia: a resting-state fMRI study with granger causality analysis



Most migraineurs develop cutaneous allodynia (CA) during migraine, and the underlying mechanism of CA in migraine is thought to be sensitization of the third-order trigeminovascular neurons in the posterior thalamic nuclei. This study aimed to investigate whether the ascending/descending pathway associated with the thalamus is disturbed in migraineurs with CA (MWCA) using effective connectivity analysis of resting-state functional magnetic resonance imaging.


Thirty four migraineurs without aura (14 MWCA and 20 migraineurs without CA (MWoCA)) and 25 matched healthy controls (HC) were recruited in the study. The effective connectivity pathways associated with the posterior thalamus (PTH) were investigated using the Granger causality analysis. We chose bilateral PTH as two individual seeds, and compared MWCA with MWoCA and HC, respectively. Spearman correlation analysis was performed to test the correlation between the abnormal effective connectivity and the allodynia severity of MWCA.


Compared with MWoCA, MWCA showed decreased inflows from the left limbic regions and dorsal medial prefrontal cortex (dmPFC) to the ipsilateral PTH, as well as increased inflow from the right ventral medial prefrontal cortex (vmPFC) to the ipsilateral PTH; no significantly different outflows from the bilateral PTH to other regions were found. Compared with HC, MWCA showed increased outflows from the left PTH to the bilateral vmPFC, decreased outflows from the right PTH to the bilateral temporoparietal areas, decreased inflow from the left parietooccipital area to the ipsilateral PTH, and increased inflows from the right dorsolateral prefrontal cortex and the bilateral temporoparietal areas to the right PTH. Correlation analyses revealed that the disturbed connectivities between PTH and cuneus, as well as PTH and middle frontal gyrus were associated with the allodynia severity of MWCA.


MWCA demonstrated disrupted effective connection pathways between the PTH and other cortical or subcortical regions that participated in multi-dimentional pain processing. Our findings highlight the dysfunctional ascending and descending pain network at the thalamic-level and may help to illuminate the possible pathophysiologic mechanisms of CA.


Migraine is a disabling and common neurological disorder, typically characterized by unilateral, throbbing, pulsating headache [1]. It affects about 12 % of general population and causes substantial personal and social burden [2]. During headache attacks, more than 60 % of migraineurs develop cutaneous allodynia (CA) [3, 4], which is defined as a painful perception of an innocuous stimulation of the skin [5]. CA is associated with high frequency of attacks, and has been proposed as a risk factor for progression to chronic migraine [6, 7]. In addition, it has been also proposed as a predictor of poor response to triptan therapy [8]. Therefore, probing the mechanism of CA in migraine may have significant implications for our understanding of the pathophysiology of migraine and its prognosis.

The hypothetical theory regarding the pathophysiological mechanism of CA is central sensitization of trigeminal neurons that take inputs from the intracranial and extracranial structures [911]. Central sensitization is a condition in which neurons have lower activation thresholds, increased responsiveness to afferent inputs, increased spontaneous activity, and enlarged receptive fields [12]. It is hypothesized that sensitization of the second-order trigeminovascular neurons in the spinal trigeminal nucleus propels the development of cephalic allodynia, and sensitization of the third-order trigeminovascular neurons in the posterior thalamic nuclei mediates the development of extracephalic allodynia [13, 14]. Functional neuroimaging studies in humans and animals have shown activation of the trigeminal neurons in the posterior/dorsal thalamic nuclei in migraine [9, 15, 16]. In addition, activation of the trigeminal system has been found to be associated with the ascending nociceptive transmission via the trigemino-thalamo-cortical pathway [17]. As we know that thalamus is a relay region subserving both sensory and motor processing, with axonal fibers not only projecting out to the cerebral cortex in all directions but also receiving feedback informations from those multiple cortical areas [18]. For pain perception, thalamus is believed to modulate ascending nociceptive information before transmitting nociceptive inputs to cortical structures [19, 20]. Furthermore, thalamus is also a part of the descending pain inhibition pathways that modulate nociceptive inputs from spinal trigeminal nucleus [21]. However, how the posterior thalamus affects the brain systems and whether it is affected by those brain regions in migraineurs with CA are poorly understood.

Effective connectivity and functional connectivity are two effective techniques to address the integration of functionally specialized areas in human brain [22]. They both allow the inference of regarding communication between spatially remote brain regions, whereas effective connectivity, utilizing information about time-lagged relationships between brain regions, allows additional inferring about the directionality of information transfer within functionally connected networks [23]. Previous resting-state functional magnetic resonance imaging (rs-fMRI) studies of migraine have found abnormal functional connectivity between thalamus and other pain matrix regions, such as periaqueductal grey (PAG), nucleus cuneiformis (NCF), or the affective regions (including insula, anterior cingulate cortex, and amygdala) [2426]. However, these studies are neither directly investigating the CA in migraine, nor are applying the thalamic-level effective connectivity analysis. Therefore, we hypothesize that studying the direct influence (i.e., effective connectivity) of thalamus with other structures can strengthen our understanding of the neural circuitry of CA in migraine.

In the present study, we aim to investigate the effective connectivity patterns of bilateral posterior thalamus with the rest of the brain in migraineurs with CA. Granger causality analysis (GCA) is a special approach to explore such effective connectivity among regions, and can be used to obtain the deductive network based on certain hypothetical seed region without requiring prior knowledge [27]. So by employing the GCA method on rs-fMRI data, we can reconstruct an effective connectivity network associated with thalamus, which would provide us with fresh insights into the descending or ascending pain pathways at the thalamic-level in migraine with CA. To the best of our knowledge, this is the first study to examine the causal interactions between the posterior thalamus and the rest brain regions in migraineurs with CA using rs-fMRI.



Thirty-four (N = 34) right-handed migraineurs (aged 19–44 years) without aura as well as 25 age-, sex-, and handness-matched healthy controls were recruited for this study. The diagnosis of migraine without aura was made according to the Second Edition of the International Classification of Headache Disorder (ICHD- II) [28]. The protocol of this study was approved by the Medical Ethics Committee of West China Hospital, Sichuan University, and all participants provided written informed consent. To avoid any possible pharmacological interference, all subjects had to be weaned off analgesic drugs 1 week or longer, not in preventive treatment and had not used any other drugs for at least 1 month prior to study participation. Subjects did not having a migraine attack at least 72 h prior to scanning. In addition, patients would be excluded if they had a migraine precipitated during the 2-day following-up. Potential subjects were excluded if they had any contraindication to MRI, had a prior brain injury, had a neurologic disorder other than migraine, had a psychiatric disorder other than anxiety or depression, or if they had any acute or chronic pain disorder other than migraine.

Clinical parameters

The Allodynia Symptom Checklist-12 (ASC-12) scale was administrated to divide all enrolled migraineurs into two groups: migraineurs with or without cutaneous allodynia (MWCA or MWoCA). ASC-12 is a 12-item questionnaire presented interictally to migraine patients, estimating the prevalence and severity of CA in the migraine population. The yielding score of each subject was placed into one of 4 categories: 0–2 = no allodynia; 3–5 = mild allodynia; 6–8 = moderate allodynia; 9 or more = severe allodynia. Then all participants were interviewed regarding demographic features (e.g., age, sex, and education). Additional information including the migraine history (e.g., onset age, frequency and duration of attacks, headache characteristic, and the pain intensity (evaluated with the visual analogue scale (VAS))) and the impact of headache (evaluated with the Migraine Disability Assessment Scale (MIDAS) [29] and Headache Impact Test (HIT)-6 [30]) was collected in those ones with migraine. For all subjects, psychiatric assessment, including 24-Hamilton Depression Scale (24- HAMD) [31] and 14-Hamilton Anxiety Scale (14-HAMA) [32] was also administered to assess the depression and anxiety state.

Statistical analysis of clinical data

Demographic and clinical features of MWCA compared to MWoCA and HC respectively were conducted using an independent sample t-test or chi square test, as appropriate.

Data acquisition

The experiment was performed on a 3.0 Tesla scanner (Trio Tim, Siemens, Erlangen, Germany) with a 16-channel birdcage head coil, and tightly padded clamps were used to minimize head motion. A routine T1 weighted imaging was first performed, and then the resting-state fMRIs were obtained by using an echo-planar imaging sequence with the following parameters: voxel size 3.75 × 3.75 × 5 mm3, TR 2000 ms, TE 30 ms, FOV 240 × 240 mm2, matrix 64 × 64, and slice thickness 5 mm with no gap, number of slices 30, number of time points 180. Throughout the scanning, subjects were instructed to lay in the scanner supine, relax, close their eyes but be awake.

Data spatial processing

The data were pre-processed using SPM8 (The Wellcome Department of Cognitive Neurology, London, UK, The first five time points of the resting-state data were discarded due to instability of the initial MRI signal, leaving 175 time points remaining for further processing. Then the functional images were slice-timing corrected and realigned to the first volume using a six-parameter rigid body transformation. Subjects with maximum head translation exceeded 2 mm or maximum rotation exceeded 2° were excluded. The mean image generated was then spatially normalized into standard stereotactic space, using the Montreal Neurological Institute (MNI) echo planar imaging (EPI) template. Computed transformation parameters were applied to all functional images, interpolated to isotropic voxels of 2 mm3 and the resulting images were smoothed using an 4-mm full-width half-maximum (FWHM), isotropic Gaussian kernel. Then, using the Data Processing Assistant for resting-state fMRI (DPARSF) package [33], linear drift was removed. A band-pass frequency filter (0.01 < f < 0.08 Hz) was then applied to reduce physiological high frequency noise [34]. To further reduce the effects of confounding factors unlikely to be involved in specific regional correlation, we also removed several sources of spurious variance by linear regression, including six head motion parameters, and average signals from cerebrospinal fluid, white matter according to previous fMRI studies [35, 36].

Granger causality analysis and statistical analysis

We used Granger causality to describe the effective connectivity between the reference time series of the seed regions (left and right posterior thalamus (PTH)) and the time series of each voxel within the whole brain. The coordinates (peak MNI coordinates: the right PTH = 3, −14, 6; the left PTH = −6, −21, −3) were obtained based on previous studies that showed structural or functional alterations in PTH in migraineurs [37, 38]. Two 6 mm-radius sphere seeds based on the peak coordinates were then designed for GCA. Bivariate first-order coefficient-based voxelwise GCA was performed using the REST-GCA in the REST toolbox ( [39]. Granger causality estimates the causal effect of the seed region on every other voxel in the brain (X to Y effect), as well as the Y to X effect, the causal effect of every voxel in the brain (Y) on the seed region (X). A positive coefficient from X to Y indicates that activity in region X exerts a causal influence on the activity region Y in the same direction (i.e., positive influence). Similarly, a negative coefficient from X to Y suggests that the activity of region X exerts an opposing directional influence on the activity of region Y (i.e., negative influence). Using this approach, we are able to build a Granger-causal model based on the temporal elements of regional BOLD activity. The generated voxel-wise GCA maps were than transformed to z scores to improve the normality [40]. Two-tailed two-sample t-tests were conducted on the causal effects to compare MWCA to MWoCA and HC respectively with a Gaussian random field (GRF)-corrected significance (voxel-level of P < 0.01 and joint cluster-level of P < 0.05). Age, sex ratio, and education were applied as covariates in the two-sample t- tests.

Clinical correlation analysis

To investigate the association between allodynia severity and abnormal effective connectivity, the regions showing significantly different (increased or decreased) Granger influences between the MWCA and MWoCA, or MWCA and HC comparisons were extracted as regions of interest (ROIs). Mean Granger causality values within these ROIs were correlated against the ASC scores of MWCA using Spearman correlation analysis. Statistical analyses were performed in SPSS 17.0 (SPSS Inc., Chicago, Illinois, USA) and threshold was set at P < 0.05.


Demographics and clinical data

Demographic and psychiatric characteristics of all subjects were given in Table 1. No subject was excluded due to head movement. Of the 34 migraineurs without aura, 14 were classified as MWCA and 20 as MWoCA according to their ASC scores. All the enrolled MWCA were with mild (3–5 scores) or moderate (6–8 scores) allodynia [4]. All subjects were right-handed. No significant difference was found for gender, age, or education either between MWCA and MWoCA, or MWCA and HC groups. Both the mean 24-HAMD and 14-HAMA scores of MWCA were significantly higher than those of HC (P value was 0.000 and 0.004, respectively), while they were not significantly different between the two migraine groups.

Table 1 Demographic and psychiatric characteristics of the study participants

Table 2 showed the headache features of the two migraine groups. In respect to the ASC-12, the MWCA group got a significantly higher score comparing to the MWoCA group (P = 0.012). There were no significant difference between the two groups in course of migraine disease and average headache intensity assessed by VAS. HIT-6 and MIDAS also achieved balance between the two migraine groups. However, MWoCA showed a significantly higher headache frequency and attack duration than MWCA (P = 0.018 and 0.028, respectively).

Table 2 Headache features of the patients with migraine

Abnormality for causal influence to and from the left PTH

Two-sample t-test results of resting-state effective connectivity to and from the left PTH were illustrated in Table 3, Figs. 1, and 3a. Compared to HC, MWCA showed significantly decreased causal inflow to the left PTH from the ipsilateral parietooccipital regions (including the cuneus, posterior cingulate cortex (PCC), and precuneus), and significantly increased causal outflow to the bilateral ventromedial prefrontal cortex (vmPFC, including the left rectus gyrus (RG), medial orbitofrontal cortex (mOFC), bilateral mSFG, and right anterior cingulate cortex (ACC)) (Fig. 1a). Both the limbic regions (including the left amygdala, hippocampus, basal ganglia (BG, including the putamen and pallidum)) and dorsomedial prefrontal cortex (dmPFC, including the left medial superior frontal gyrus (mSFG), supplementary motor area (SMA), and bilateral middle cingulate cortex (MCC)) had decreased influence to the ipsilateral PTH in MWCA comparing to those in MWoCA (Fig. 1b). No significantly different causal outflow from the left PTH to rest of brain was found between two migraine comparisons.

Table 3 Two sample t-test (voxel-level P < 0.01 and cluster-level P < 0.05 Gaussian random field corrected) of difference for causal influence to and from the left posterior thalamus
Fig. 1
figure 1

Altered effective connectivity to and from the left PTH in MWCA compared with HC (a) and MWoCA (b) (voxel-level P < 0.01 and cluster-level P < 0.05, Gaussian random field corrected) respectively. Compared with HC, MWCA shows decreased inflow from the L_PTH to the left cuneus, precuneus, and PCC, and increased outflow to the left RG, mOFC, SFG, and ACC (a). Compared with MWoCA, MWCA shows decreased inflow from ipsilateral amygdala, hippocampus, putamen, pallidum, mSFG, SMA, and MCC (b). ACC, anterior cingulate cortex; HC, healthy control; L_PTH, right posterior thalamus; MCC, middle cingulate cortex; mOFC, medial orbitofrontal cortex; mSFG, medial superior frontal gyrus; MWCA, migraineurs with cutaneous allodynia; MWoCA, migraineurs without cutaneous allodynia; PCC, posterior cingulate cortex; PTH, posterior thalamus; RG, rectus gyrus; SFG, superior frontal gyrus; SMA, supplementary motor area

Abnormality for causal influence to and from the right PTH

As shown in Table 4, Figs. 2a, and 3b, compared with HC, MWCA group demonstrated significantly increased causal inflow to the right PTH from the bilateral temporoparietal areas (including the superior temporal gyrus (STG) and inferior parietal lobule (IPL, including the supramarginal gyrus and angular gyrus)), as well as decreased causal outflow to the bilateral temporoparietal areas; MWCA also exhibited increased causal inflow to the right PTH from the ipsilateral dorsolateral prefrontal cortex (DLPFC, including the opercular, triangular, and orbital part of right inferior frontal gyrus (IFG) and the orbital part of middle frontal gyrus (MFG)). Compared to MWoCA, MWCA had increased inflow to the right PTH from the ipsilateral vmPFC (including the olfactory cortex and caudate nucleus) (Fig. 2b) and no significantly different causal outflow from the right PTH to rest of brain was found.

Table 4 Two sample t-test (voxel-level P < 0.01 and cluster-level P < 0.05 Gaussian random field corrected) of difference for causal influence to and from the right posterior thalamus
Fig. 2
figure 2

Altered effective connectivity to and from the right PTH in MWCA compared with HC (a) and MWoCA (b) (voxel-level P < 0.01 and cluster-level P < 0.05, Gaussian random field corrected) respectively. Compared with HC, MWCA showed increased inflow to R_PTH from the bilateral STG, SMG, IPL, IFG, and MFG, and decreased outflow to the same regions (a). Compared with MWoCA, MWCA showed increased inflow from the OC and caudate to the R_PTH (b). HC, healthy control; IFG, inferior frontal gyrus; IPL, inferior parietal lobule; MFG, middle frontal gyrus; MWCA, migraineurs with cutaneous allodynia; MWoCA, migraineurs without cutaneous allodynia; OC, olfactory cortex; PTH, posterior thalamus; R_PTH, right posterior thalamus; SMG, supramarginal gyrus; STG, superior temporal gyrus

Fig. 3
figure 3

Illustration of the abnormal thalamic-level Granger-causal model in MWCA compared with the two control groups. a Abnormal effective connectivity pathways associated with the Left PTH. b Abnormal effective connectivity pathways associated with the right PTH. The node size is represented by absolute T-value of corresponding cluster. The red and blue flags represent the positive and negative correlation between the strength of abnormal effective connectivity pathway and ASC (index of allodynia severity), respectively. DLPFC, dorsolateral prefrontal cortex; dmPFC, dorsal medial prefrontal cortex; HC, healthy control; L, left; MWCA, migraineurs with cutaneous allodynia; MWoCA, migraineurs without cutaneous allodynia; PTH, posterior thalamus; R, right; vmPFC, ventral medial prefrontal cortex

Correlations results

As characterized in Fig. 3, negative correlation was found between the ASC and causal inflow from the left cuneus to the left PTH in MWCA vs. HC (r = −0.72, P = 0.004). Positive correlation was found between the ASC and causal inflow from the right MFG to the right PTH (r = 0.62, P = 0.019) in MWCA vs. HC.


The present rs-fMRI study with Granger causality analysis demonstrated effective connectivity alterations from and to the bilateral PTH in MWCA comparing to MWoCA and HC respectively. Specifically, the abnormal effective connectivity pathways included: 1) one thalamo-cortico-thalamic circuit and one ipsilateral feedback inflow to the PTH in MWCA compared with MWoCA, i.e., bidirectional connectivity between the bilateral PTH and the mPFC (including the ventral and dorsal part), and inflow from the left limbic regions (including the amygdala, hippocampus, and BG) to the ipsilateral PTH; 2) one bidirectional thalamic-cortical circuit and another two ipsilateral feedback inflows to the PTH in MWCA compared with HC, i.e., bidirectional connectivity between the right PTH and the bilateral temporoparietal regions (including the STG and IPL), inflows from the left parietooccipital regions (including the precuneus, cuneus, and PCC) to the ipsilateral PTH and from right DLPFC (including the IFG and MFG) to the ipsilateral PTH. Moreover, the disturbed effective connectivities between PTH and cuneus, as well as PTH and MFG were associated with the severity of allodynia. These results confirmed the previous notion that central sensitisation of the third-order trigeminal neurons in the thalamus was implicated in the mechanism of CA in migraine [14, 41]. Additionally, some regions in those abnormal effective connectivity pathways were involved in affective, cognitive, and sensory discriminative dimensions of pain, which suggested that CA implicated a multi-dimensional pain processing dysfunction.

One remarkable finding of the present study was the abnormally causal connectivity between the bilateral PTH and the mPFC (including the dorsal and ventral mPFC) in both MWCA vs. MWoCA and MWCA vs. HC. Results of the comparison between the two migraine groups suggested that CA may be associated with the decreased inflow from the left dmPFC to the ipsilateral PTH and increased inflow from the right vmPFC to the ipsilateral PTH. mPFC has widespread anatomical and functional connections with thalamus [42], and is thought to play a critical role in emotional-cognitive processing [43]. Furthermore, the two subregions of mPFC (i.e., vmPFC and dmPFC) have functional differentces,: the dmPFC is involved in evaluation of the perceived pain sensation [44], while the vmPFC has an important role in regulation of sensory and affective pain [45]. Therefore, the present results suggested that patients with CA had an impaired ability of evaluation of the perceived pain sensation. This was consistent with the definition of CA which is defined as pain (an exaggeration of the sensory inputs) resulting from the application of a non-noxious stimulus such as heat, cold or pressure to normal skin. Results also suggested that MWCA showed increased descending pain modulation which might be considered a response to the exaggerated pain information. Previous study has suggested that the medial thalamus and the mesial cortex (including the mPFC and ACC) involve the affective-motivational component of pain [46]. Other resting-state fMRI studies also found abnormal mPFC activity to be implicated in the abnormal affective pain perception, cognition, and modulation in migraineurs [25, 36, 47]. However, whether the dysfunctional effective connectivities between PTH and mPFC contribute to CA or are just cumulative effective results of CA still merits further study.

Another abnormal effective connectivity associated with CA was the decreased influence from the amygdala, hippocampus, putamen, and pallidum to the left PTH. All these subcortical regions have dense connections to thalamus [18, 48], and have been implicated in the cognitive/affective aspect of pain processing in migraine [4952]. The involvement of these subcortical regions may correlate with the strongly relevant individual concern, like fears and anxiety [53] resulting from recurring and repetitive attacks of migraine [54]. Furthermore, the result is supported by a previous fMRI study which found that migraine patients with very high frequency of migraine attacks per month in response to the thermal stimuli demonstrated significantly lower activation throughout the caudate, putamen, and pallidum nuclei of the BG than those with low frequency of migraine attacks [51]. Therefore, based on these findings, we speculate that present abnormal cognitive and emotional pain processes may result from the accumulating damages from recurring and repetitive attacks, which may involve the cognitive-emotional impairment in migraineurs with CA.

Reduced influence from the left parietooccipital area (including cuneus, precuneus, and PCC) to the ipsilateral PTH was also found in MWCA comparing with HC, and was negatively correlated with the ASC (allodynia assessment index), which suggested the abnormal feedback flow through the parietooccipital-thalamic circuit to be implicated in the underlying mechanism of CA. Most of these regions are important parts of the default mode network (DMN), which is the most consistently reported intrinsic connectivity networks during resting-state and has been implicated in sensation integration, self-relevant, and internally cognitive-attentional dimensions of pain in migraineurs [55, 56]. Decreased connectivity between the parietooccipital regions and PTH reported herein may suggest an abnormal feedback from DMN regions to PTH, and certainly affects the following pain processing as relevant to the DMN function in migraine. In addition, the results also found increased influence from the right DLPFC (including IFG and MFG) to ipsilateral PTH in MWCA vs. HC, as well as its positive relationship with ASC. Concerning pain perception, the prefrontal regions including IFG have an important role in modulating pain processing during attentional manipulation paradigms [57]. DLPFC, along with medial thalamus, has also been reported to engage in the pain processing during heat allodynia [58]. Thus we speculate that the DLPFC exerts abnormal active control on sensory integration involving in the CA through the cortico-thalamus pathway.

The result revealed an abnormal bidirectional interaction between the right PTH and the bilateral temporoparietal areas (including the STG and IPL) in MWCA vs. HC. Specifically, there were increased inflow from the bilateral temporoparietal area to the right PTH and decreased outflow from the right PTH to the bilateral temporoparietal areas. Previous morphological and functional studies in migraine showed that migraineurs have reduced gray matter in STG and IPL [59, 60], as well as less pain-induced activation of STG than HC in the interictal state [61]. STG is an essential structure involved in auditory processing and social cognition [62], and is generally associated with the IPL to serve the so-called ventral stream of the visual pathway [63], which acts as a link between auditory and visual processing, perception and memory [64]. When taken together with the functional properties of these brain regions, the present findings suggested that migraineurs suffered from a condition of global dysfunction of sensory integration and memory processes, as well as abnormal descending modulation through the cortico-thalamic pathway.

Study limitations

There are several limitations in the present study. First, the sample size of patients was small. The severity of allodynia symptoms, the length of disease course, the history of medication use and other characteristics were possible confounding variables, but we could not make any subgroup comparisons due to the few participants. Second, the current study did not present the large results of the MWoCA vs. HC, because we exclusively focused on the results that was directly associated with CA comparing MWCA with either MWoCA or HC. Finally, we have not performed a long-term evaluation or rescanned the subjects yet, so whether the functional alteration is a predictor or a consequence of central sensitization cannot be determined by the present study.


The present study provides evidence of the involvement of the third-order trigeminovascular neurons in the posterior thalamic nuclei in the pathophysiological mechanism of CA in migraine. The disrupted effective connectivity pathways between PTH and other cortical or subcortical regions (such as the prefrontal cortex, limbic regions, parietooccipital and temporoparietal areas) exhibit a dysfunctional multi-dimensional (including the sensory-discriminative, cognitive, and affective domains) pain processing pattern in MWCA. These findings highlight the dysfunctional ascending and descending pain network at the thalamic-level. In a future study, the implicated regions in the present study will be chosen as seed regions for further effective connectivity analysis, in order to construct a comprehensive pathway model specifically accounting for CA in migraine.



anterior cingulate cortex


Allodynia symptom checklist


basal ganglia


cutaneous allodynia


dorsolateral prefrontal cortex


default mode network


dorsal medial prefrontal cortex


granger causality analysis


healthy control


inferior frontal gyrus


inferior parietal lobule


middle cingulate cortex


middle frontal gyrus


medial orbitofrontal cortex


medial prefrontal cortex


medial superior frontal cortex


migraineurs with cutaneous allodynia


migraineurs without cutaneous allodynia


olfactory cortex


posterior cingulate cortex


posterior thalamus


rectus gyrus


resting-state functional magnetic resonance imaging


supplementary motor area


supramarginal gyrus


superior temporal gyrus


ventral medial prefrontal cortex


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

    Article  Google Scholar 

  2. Lipton RB, Bigal ME, Diamond M, Freitag F, Reed ML, Stewart WF, Group AA (2007) Migraine prevalence, disease burden, and the need for preventive therapy. Neurology 68(5):343–349. doi:10.1212/01.wnl.0000252808.97649.21

    Article  Google Scholar 

  3. de Tommaso M, Guido M, Libro G, Losito L, Sciruicchio V, Monetti C, Puca F (2002) Abnormal brain processing of cutaneous pain in migraine patients during the attack. Neurosci Lett 333(1):29–32

    Article  PubMed  Google Scholar 

  4. Lipton RB, Bigal ME, Ashina S, Burstein R, Silberstein S, Reed ML, Serrano D, Stewart WF, American Migraine Prevalence Prevention Advisory G (2008) Cutaneous allodynia in the migraine population. Ann Neurol 63(2):148–158. doi:10.1002/ana.21211

    Article  Google Scholar 

  5. Aguggia M (2012) Allodynia and migraine. Neurol Sci 33(Suppl 1):S9–S11. doi:10.1007/s10072-012-1034-9

    Article  PubMed  Google Scholar 

  6. Bigal ME, Lipton RB (2006) Modifiable risk factors for migraine progression. Headache 46(9):1334–1343. doi:10.1111/j.1526-4610.2006.00577.x

    Article  PubMed  Google Scholar 

  7. Louter MA, Bosker JE, van Oosterhout WPJ, van Zwet EW, Zitman FG, Ferrari MD, Terwindt GM (2013) Cutaneous allodynia as a predictor of migraine chronification. Brain 136:3489–3496. doi:10.1093/brain/awt251

    Article  PubMed  Google Scholar 

  8. Burstein R, Collins B, Jakubowski M (2004) Defeating migraine pain with triptans: a race against the development of cutaneous allodynia. Ann Neurol 55(1):19–26. doi:10.1002/ana.10786

    Article  CAS  PubMed  Google Scholar 

  9. Burstein R, Jakubowski M, Garcia-Nicas E, Kainz V, Bajwa Z, Hargreaves R, Becerra L, Borsook D (2010) Thalamic sensitization transforms localized pain into widespread allodynia. Ann Neurol 68(1):81–91. doi:10.1002/ana.21994

    Article  PubMed Central  PubMed  Google Scholar 

  10. Burstein R, Yarnitsky D, Goor-Aryeh I, Ransil BJ, Bajwa ZH (2000) An association between migraine and cutaneous allodynia. Ann Neurol 47(5):614–624

    Article  CAS  PubMed  Google Scholar 

  11. Tietjen GE, Brandes JL, Peterlin BL, Eloff A, Dafer RM, Stein MR, Drexler E, Martin VT, Hutchinson S, Aurora SK, Recober A, Herial NA, Utley C, White L, Khuder SA (2009) Allodynia in migraine: association with comorbid pain conditions. Headache 49(9):1333–1344. doi:10.1111/j.1526-4610.2009.01521.x

    Article  PubMed  Google Scholar 

  12. Lee MC, Zambreanu L, Menon DK, Tracey I (2008) Identifying brain activity specifically related to the maintenance and perceptual consequence of central sensitization in humans. J Neurosci 28(45):11642–11649. doi:10.1523/JNEUROSCI.2638-08.2008

    Article  CAS  PubMed  Google Scholar 

  13. Bernstein C, Burstein R (2012) Sensitization of the trigeminovascular pathway: perspective and implications to migraine pathophysiology. J Clin Neurol 8(2):89–99. doi:10.3988/jcn.2012.8.2.89

    Article  PubMed Central  PubMed  Google Scholar 

  14. Landy S, Rice K, Lobo B (2004) Central sensitisation and cutaneous allodynia in migraine: implications for treatment. CNS Drugs 18(6):337–342

    Article  PubMed  Google Scholar 

  15. Shields KG, Goadsby PJ (2005) Propranolol modulates trigeminovascular responses in thalamic ventroposteromedial nucleus: a role in migraine? Brain 128(Pt 1):86–97. doi:10.1093/brain/awh298

    PubMed  Google Scholar 

  16. Zagami AS, Lambert GA (1990) Stimulation of cranial vessels excites nociceptive neurones in several thalamic nuclei of the cat. Exp Brain Res 81(3):552–566

    Article  CAS  PubMed  Google Scholar 

  17. Sokolov AY, Lyubashina OA, Sivachenko IB, Panteleev SS (2014) Effects of intravenous metamizole on ongoing and evoked activity of dura-sensitive thalamic neurons in rats. Eur J Pharmacol 731:58–64. doi:10.1016/j.ejphar.2014.03.006

    Article  CAS  PubMed  Google Scholar 

  18. Herrero MT, Barcia C, Navarro JM (2002) Functional anatomy of thalamus and basal ganglia. Childs Nerv Syst 18(8):386–404. doi:10.1007/s00381-002-0604-1

    Article  PubMed  Google Scholar 

  19. Kupers R, Kehlet H (2006) Brain imaging of clinical pain states: a critical review and strategies for future studies. Lancet Neurol 5(12):1033–1044. doi:10.1016/S1474-4422(06)70624-X

    Article  PubMed  Google Scholar 

  20. Noseda R, Jakubowski M, Kainz V, Borsook D, Burstein R (2011) Cortical projections of functionally identified thalamic trigeminovascular neurons: implications for migraine headache and its associated symptoms. J Neurosci 31(40):14204–14217. doi:10.1523/JNEUROSCI.3285-11.2011

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Ab Aziz CB, Ahmad AH (2006) The role of the thalamus in modulating pain. MJMS 13(2):11–18

    PubMed Central  PubMed  Google Scholar 

  22. Friston KJ (1994) Functional and effective connectivity in neuroimaging: a synthesis. Hum Brain Mapp 2(1–2):56–78. doi:10.1002/hbm.460020107

    Article  Google Scholar 

  23. Leavitt VM, Wylie G, Genova HM, Chiaravalloti ND, DeLuca J (2012) Altered effective connectivity during performance of an information processing speed task in multiple sclerosis. Mult Scler 18(4):409–417. doi:10.1177/1352458511423651

    Article  PubMed  Google Scholar 

  24. Schwedt TJ, Larson-Prior L, Coalson RS, Nolan T, Mar S, Ances BM, Benzinger T, Schlaggar BL (2014) Allodynia and descending pain modulation in migraine: a resting state functional connectivity analysis. Pain Med 15(1):154–165. doi:10.1111/pme.12267

    Article  PubMed Central  PubMed  Google Scholar 

  25. Schwedt TJ, Schlaggar BL, Mar S, Nolan T, Coalson RS, Nardos B, Benzinger T, Larson-Prior LJ (2013) Atypical resting-state functional connectivity of affective pain regions in chronic migraine. Headache 53(5):737–751. doi:10.1111/head.12081

    Article  PubMed Central  PubMed  Google Scholar 

  26. Xue T, Yuan K, Cheng P, Zhao L, Zhao L, Yu D, Dong T, von Deneen KM, Gong Q, Qin W, Tian J (2013) Alterations of regional spontaneous neuronal activity and corresponding brain circuit changes during resting state in migraine without aura. NMR Biomed 26(9):1051–1058. doi:10.1002/nbm.2917

    Article  PubMed  Google Scholar 

  27. Friston K (2009) Causal modelling and brain connectivity in functional magnetic resonance imaging. PLoS Biol 7(2), e33. doi:10.1371/journal.pbio.1000033

    Article  PubMed  Google Scholar 

  28. Silberstein SD, Olesen J, Bousser MG, Diener HC, Dodick D, First M, Goadsby PJ, Gobel H, Lainez MJ, Lance JW, Lipton RB, Nappi G, Sakai F, Schoenen J, Steiner TJ, International Headache S (2005) The international classification of headache disorders, 2nd edition (ICHD-II)--revision of criteria for 8.2 medication-overuse headache. Cephalalgia 25(6):460–465. doi:10.1111/j.1468-2982.2005.00878.x

    Article  CAS  PubMed  Google Scholar 

  29. Stewart WF, Lipton RB, Kolodner KB, Sawyer J, Lee C, Liberman JN (2000) Validity of the migraine disability assessment (MIDAS) score in comparison to a diary-based measure in a population sample of migraine sufferers. Pain 88(1):41–52

    Article  CAS  PubMed  Google Scholar 

  30. Shin HE, Park JW, Kim YI, Lee KS (2008) Headache impact test-6 (HIT-6) scores for migraine patients: their relation to disability as measured from a headache diary. J Clin Neurol 4(4):158–163. doi:10.3988/jcn.2008.4.4.158

    Article  PubMed Central  PubMed  Google Scholar 

  31. Hamilton M (1967) Development of a rating scale for primary depressive illness. Br J Soc Clin Psychol 6(4):278–296

    Article  CAS  PubMed  Google Scholar 

  32. Hamilton M (1959) The assessment of anxiety states by rating. Br J Med Psychol 32(1):50–55

    Article  CAS  PubMed  Google Scholar 

  33. Chao-Gan Y, Yu-Feng Z (2010) DPARSF: A MATLAB Toolbox for “Pipeline” Data Analysis of Resting-State fMRI. Front Syst Neurosci 4:13. doi:10.3389/fnsys.2010.00013

    PubMed Central  PubMed  Google Scholar 

  34. Biswal B, Yetkin FZ, Haughton VM, Hyde JS (1995) Functional Connectivity in the Motor Cortex of Resting Human Brain Using Echo-Planar Mri. Magnet Reson Med 34 (4):537-541. doi: 10.1002/mrm.1910340409

  35. Fox MD, Snyder AZ, Vincent JL, Corbetta M, Van Essen DC, Raichle ME (2005) The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci U S A 102(27):9673–9678. doi:10.1073/pnas.0504136102

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Liao W, Qiu C, Gentili C, Walter M, Pan Z, Ding J, Zhang W, Gong Q, Chen H (2010) Altered effective connectivity network of the amygdala in social anxiety disorder: a resting-state FMRI study. PLoS One 5(12), e15238. doi:10.1371/journal.pone.0015238

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Zhao L, Liu JX, Dong XL, Peng YL, Yuan K, Wu FM, Sun JB, Gong QY, Qin W, Liang FR (2013) Alterations in regional homogeneity assessed by fMRI in patients with migraine without aura stratified by disease duration. J Headache Pain 14:85. doi:10.1186/1129-2377-14-85

    Article  PubMed Central  PubMed  Google Scholar 

  38. Stankewitz A, Schulz E, May A (2013) Neuronal correlates of impaired habituation in response to repeated trigemino-nociceptive but not to olfactory input in migraineurs: an fMRI study. Cephalalgia 33(4):256–265. doi:10.1177/0333102412470215

    Article  CAS  PubMed  Google Scholar 

  39. Song XW, Dong ZY, Long XY, Li SF, Zuo XN, Zhu CZ, He Y, Yan CG, Zang YF (2011) REST: a toolkit for resting-state functional magnetic resonance imaging data processing. PLoS One 6(9), e25031. doi:10.1371/journal.pone.0025031

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Zang ZX, Yan CG, Dong ZY, Huang J, Zang YF (2012) Granger causality analysis implementation on MATLAB: a graphic user interface toolkit for fMRI data processing. J Neurosci Methods 203(2):418–426. doi:10.1016/j.jneumeth.2011.10.006

    Article  PubMed  Google Scholar 

  41. Aguggia M, Saracco MG, Cavallini M, Bussone G, Cortelli P (2013) Sensitization and pain. Neurol Sci 34(1):37–40. doi:10.1007/s10072-013-1382-0

    Article  Google Scholar 

  42. Mufson EJ, Mesulam MM (1984) Thalamic connections of the insula in the rhesus monkey and comments on the paralimbic connectivity of the medial pulvinar nucleus. J Comp Neurol 227(1):109–120. doi:10.1002/cne.902270112

    Article  CAS  PubMed  Google Scholar 

  43. Etkin A, Egner T, Kalisch R (2011) Emotional processing in anterior cingulate and medial prefrontal cortex. Trends Cogn Sci 15(2):85–93. doi:10.1016/j.tics.2010.11.004

    Article  PubMed Central  PubMed  Google Scholar 

  44. Napadow V, Dhond RP, Kim J, LaCount L, Vangel M, Harris RE, Kettner N, Park K (2009) Brain encoding of acupuncture sensation--coupling on-line rating with fMRI. NeuroImage 47(3):1055–1065. doi:10.1016/j.neuroimage.2009.05.079

    Article  PubMed Central  PubMed  Google Scholar 

  45. Wang GQ, Cen C, Li C, Cao S, Wang N, Zhou Z, Liu XM, Xu Y, Tian NX, Zhang Y, Wang J, Wang LP, Wang Y (2015) Deactivation of excitatory neurons in the prelimbic cortex via Cdk5 promotes pain sensation and anxiety. Nat Commun 6:7660. doi:10.1038/ncomms8660

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Sewards TV, Sewards MA (2003) Representations of motivational drives in mesial cortex, medial thalamus, hypothalamus and midbrain. Brain Res Bull 61(1):25–49

    Article  PubMed  Google Scholar 

  47. Jin C, Yuan K, Zhao L, Zhao L, Yu D, von Deneen KM, Zhang M, Qin W, Sun W, Tian J (2013) Structural and functional abnormalities in migraine patients without aura. NMR Biomed 26(1):58–64. doi:10.1002/nbm.2819

    Article  PubMed  Google Scholar 

  48. Aggleton JP, Mishkin M (1984) Projections of the amygdala to the thalamus in the cynomolgus monkey. J Comp Neurol 222(1):56–68. doi:10.1002/cne.902220106

    Article  CAS  PubMed  Google Scholar 

  49. Faria V, Erpelding N, Lebel A, Johnson A, Wolff R, Fair D, Burstein R, Becerra L, Borsook D (2015) The migraine brain in transition: girls versus Boys. Pain. doi:10.1097/j.pain.0000000000000292

    PubMed  Google Scholar 

  50. Maleki N, Becerra L, Brawn J, McEwen B, Burstein R, Borsook D (2013) Common hippocampal structural and functional changes in migraine. Brain Struct Funct 218(4):903–912. doi:10.1007/s00429-012-0437-y

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Maleki N, Becerra L, Nutile L, Pendse G, Brawn J, Bigal M, Burstein R, Borsook D (2011) Migraine attacks the Basal Ganglia. Mol Pain 7:71. doi:10.1186/1744-8069-7-71

    Article  PubMed Central  PubMed  Google Scholar 

  52. Neugebauer V (2015) Amygdala pain mechanisms. Handb Exp Pharmacol 227:261–284. doi:10.1007/978-3-662-46450-2_13

    Article  PubMed Central  PubMed  Google Scholar 

  53. Brosschot JF (2002) Cognitive-emotional sensitization and somatic health complaints. Scand J Psychol 43(2):113–121

    Article  PubMed  Google Scholar 

  54. Smitherman TA, Kolivas ED, Bailey JR (2013) Panic disorder and migraine: comorbidity, mechanisms, and clinical implications. Headache 53(1):23–45. doi:10.1111/head.12004

    Article  PubMed  Google Scholar 

  55. Xue T, Yuan K, Zhao L, Yu DH, Zhao LM, Dong T, Cheng P, von Deneen KM, Qin W, Tian J (2012) Intrinsic brain network abnormalities in migraines without aura revealed in resting-state fMRI. PLoS One 7(12)), e52927. doi:10.1371/journal.pone.0052927

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  56. Tessitore A, Russo A, Giordano A, Conte F, Corbo D, De Stefano M, Cirillo S, Cirillo M, Esposito F, Tedeschi G (2013) Disrupted default mode network connectivity in migraine without aura. J Headache Pain 14:89. doi:10.1186/1129-2377-14-89

    Article  PubMed Central  PubMed  Google Scholar 

  57. Petrovic P, Petersson KM, Ghatan PH, Stone-Elander S, Ingvar M (2000) Pain-related cerebral activation is altered by a distracting cognitive task. Pain 85(1–2):19–30. doi:10.1016/S0304-3959(99)00232-8

    Article  CAS  PubMed  Google Scholar 

  58. 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. doi:10.1093/brain/awg102

    Article  CAS  PubMed  Google Scholar 

  59. Valfre W, Rainero I, Bergui M, Pinessi L (2008) Voxel-based morphometry reveals gray matter abnormalities in migraine. Headache 48(1):109–117. doi:10.1111/j.1526-4610.2007.00723.x

    Article  PubMed  Google Scholar 

  60. Coppola G, Di Renzo A, Tinelli E, Iacovelli E, Lepre C, Di Lorenzo C, Di Lorenzo G, Di Lenola D, Parisi V, Serrao M, Pauri F, Fiermonte G, Bianco F, Pierelli F (2015) Evidence for brain morphometric changes during the migraine cycle: a magnetic resonance-based morphometry study. Cephalalgia 35(9):783–791. doi:10.1177/0333102414559732

    Article  PubMed  Google Scholar 

  61. 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(12):947–958. doi:10.1177/0333102414526069

    Article  PubMed Central  PubMed  Google Scholar 

  62. Bigler ED, Mortensen S, Neeley ES, Ozonoff S, Krasny L, Johnson M, Lu J, Provencal SL, McMahon W, Lainhart JE (2007) Superior temporal gyrus, language function, and autism. Dev Neuropsychol 31(2):217–238. doi:10.1080/87565640701190841

    Article  PubMed  Google Scholar 

  63. Kravitz DJ, Saleem KS, Baker CI, Ungerleider LG, Mishkin M (2013) The ventral visual pathway: an expanded neural framework for the processing of object quality. Trends Cogn Sci 17(1):26–49. doi:10.1016/j.tics.2012.10.011

    Article  PubMed Central  PubMed  Google Scholar 

  64. Shapiro K, Hillstrom AP, Husain M (2002) Control of visuotemporal attention by inferior parietal and superior temporal cortex. Curr Biol 12(15):1320–1325

    Article  CAS  PubMed  Google Scholar 

Download references


This study was supported by National Natural Science Foundation (Grant Nos.81371536, 81000605, 81500959, and 81220108013), Natural Science Foundation of Guangdong (Grant No.S20120200-10867), Sichuan province science and technology plan project (Grant No.2015HH0036), and the Program for Changjiang Scholars and Innovative Research Team in University of China.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Junran Zhang or Qiyong Gong.

Additional information

Competing interest

The authors declare that they have no competing interests.

Authors’ contribution

TW: experimental design, image data analysis, results interpretation, manuscript drafting; WZ: manuscript revision, Grammatical error correction; NC: data acquisition, results interpretation, manuscript drafting; JL: results interpretation, manuscript drafting; JPZ, QL, HH, and LH: manuscript revision; JRZ and QYG: conception and design of the study, results interpretation, manuscript revision; All authors read and approved the final manuscript.

Ting Wang and Ning Chen are co-first authors

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, T., Chen, N., Zhan, W. et al. Altered effective connectivity of posterior thalamus in migraine with cutaneous allodynia: a resting-state fMRI study with granger causality analysis. J Headache Pain 17, 17 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: