Abstract
Background
Previous functional MRI studies have revealed that ongoing clinical pain in different chronic pain syndromes is directly correlated to the connectivity strength of the resting default mode network (DMN) with the insula. Here, we investigated seed-based resting state DMN-insula connectivity during acute migraine headaches.
Methods
Thirteen migraine without aura patients (MI) underwent 3 T MRI scans during the initial six hours of a spontaneous migraine attack, and were compared to a group of 19 healthy volunteers (HV). We evaluated headache intensity with a visual analogue scale and collected seed-based MRI resting state data in the four core regions of the DMN: Medial prefrontal cortex (MPFC), posterior cingulate cortex (PCC), and left and right inferior parietal lobules (IPLs), as well as in bilateral insula.
Results
Compared to HV, MI patients showed stronger functional connectivity between MPFC and PCC, and between MPFC and bilateral insula. During migraine attacks, the strength of MPFC-to-insula connectivity was negatively correlated with pain intensity.
Conclusion
We show that greater subjective intensity of pain during a migraine attack is associated with proportionally weaker DMN-insula connectivity. This is at variance with other chronic extra-cephalic pain disorders where the opposite was found, and may thus be a hallmark of acute migraine head pain.
Introduction
During the last decades, neuroimaging methods have greatly contributed to the understanding of the central nervous system mechanisms involved in migraine pathophysiology, showing that the disorder is associated with macro- and microstructural, neurochemical, and functional alterations in various brain areas (1). However, despite enormous progress in the field, the brain circuitry subtending and encoding the perception of head pain is still partly elusive and presents several methodological obstacles. Most studies have assessed functional brain changes in response to experimentally induced noxious stimuli or during induced migraine attacks by nitric oxide donors, various other vasodilating chemicals or hypoxia (2–5). PET recordings during spontaneous attacks have focused on brain stem and photophobia-related activations (6–10), but the central nervous system underpinnings of headache intensity were not assessed.
Among the various neuroimaging procedures, resting state functional MRI (RS-fMRI) analyses spontaneous BOLD signal fluctuations not attributable to an explicit in- or output, and determines between which brain areas resting activity is related (11). Seed-based or region of interest (ROI) analysis is one method of identifying networks and connectivity in the brain. In this case, the signal from only a certain voxel or a cluster of voxels, known as the seed or ROI, is used to calculate its correlations with other voxels of the brain, which provides a more precise and detailed look at specific connectivity in brain areas of interest. Seed-based analysis of RS-fMRI data is thus suitable for examining the brain activity and connections associated with spontaneous pain, i.e. in the absence of any experimental intervention.
As in certain chronic pain conditions (12–14), seed-based RS-fMRI studies in migraine between attacks have provided evidence of a disruption of the default mode network (DMN) (15–20), a set of brain regions active at rest and involved in self-referential orientation and monitoring (21). Intriguingly, in chronic pain disorders, such as fibromyalgia and low-back pain, such studies have revealed a positive correlation between the intensity of pain and resting DMN connectivity with the insula (INS) (22,23) and this correlation persists with manoeuvre-induced exacerbation of pain in chronic low back pain patients (24). To the best of our knowledge, there are no RS-fMRI studies using seed-based analysis to determine the temporal dynamics of the functional connectivity between DMN and insula during spontaneous migraine attacks, or a fortiori its correlation with headache intensity.
We decided therefore to use a seed-based approach to search for abnormal functional connectivity patterns of the DMN and between the DMN, insula, and the severity of head pain in episodic migraine without aura patients scanned during an attack. We hypothesized that, akin to other pain disorders, migraine headache would be associated with increased connectivity of the DMN with insula and the strength of this connectivity would be proportional to headache intensity.
Materials and methods
Participants
Clinical and demographic characteristics of healthy volunteers (HV) and migraine without aura patients recorded during (MI) attacks. Data are expressed as means ± SD.
All participants received a complete description of the study and granted written informed consent. The institutional review board of the Faculty of Pharmacy and Medicine, University of Rome, Polo Pontino, Italy, approved the project. The study complied with the principles set out in the WMA Declaration of Helsinki – Ethical principles for medical research involving human subjects.
Imaging protocols
MRI data were obtained on a Siemens 3 T Verio scanner using a 12-channel head coil. Structural anatomic scans were performed using a T1-weighted sagittal magnetization-prepared rapid gradient echo (MP-RAGE) series (TR: 1900 ms, TE: 2.93 ms, 176 sagittal slices, 0.508 × 0.508 × 1 mm3 voxels). We acquired an interleaved double-echo Turbo Spin Echo sequence proton density and T2-weighted images (repetition time: 3320 ms, echo time: 10/103 ms, matrix: 384 ×384, field of view: 220 mm, slice thickness: 4 mm, gap: 1.2 mm, 50 axial slices). Functional MRI data were obtained using T2*-weighted, echo-planar imaging (TR: 3000 ms, TE: 30 ms, 40 axial slices, 3.906 ×3.906 × 3 mm, 150 volumes). To accelerate fMRI acquisitions and minimize distortions, the actual EPI uses multiband sequences with simultaneous echo refocusing and parallel imaging (called Generalized Autocalibrating Partially Parallel Acquisition, GRAPPA, by Siemens). Functional resting scans lasted seven minutes and 30 seconds, during which participants were instructed to relax, avoid motion, and keep their eyes closed.
Data processing and analyses
All images were processed using SPM12 (Statistical Parametric Mapping, version 12). Functional connectivity analyses were performed using the CONN v16b connectivity toolbox (28) in MATLAB (MathWorks, Inc).
Pre-processing involved the following steps based on SPM12 algorithms. All images from a single subject were realigned using a six-parameter rigid body process, resliced by a cubic spline interpolation; outcomes, motion, and rotations were recorded in a text file for every subject and used as first level covariates. The structural (T1 – MPRAGE) and functional data were co-registered for each subject dataset; data were transformed into a common stereotactic space based on Talairach and Tournoux (29) and resampled isotropically at 3 × 3 × 3 mm. After that, the spatially normalized functional images were smoothed by 8 mm in each direction.
The processing steps listed below were performed using CONN, following the procedure reported elsewhere (30). First, each participant’s structural dataset was segmented to grey matter, white matter, and cerebrospinal fluid (CSF). The pre-process step removes sources of possible confounders: BOLD signals from white matter and CSF, realignment parameters of each subject as covariates. Moreover, a scrubbing pre-processing step based on the CONN artefact detector tool (ART) detected the offending scans for every subject; therefore the output text files were used in the denoising step for scrubbing. None of the participants’ volumes was removed from the analysis due to excessive motion. The rest effect condition convolved with hemodynamic response function and the bandpass filter values were 0.01 and 0.1 Hz.
First level functional connectivity analysis was performed using the seed-based approach. We collected seed-based resting state data in the four core regions of the DMN: Medial prefrontal cortex (MPFC), posterior cingulate cortex (PCC), and bilateral inferior parietal lobules (IPLs). In addition, we collected seed-based resting state data from the insula bilaterally. All seeds – or region of interests, ROIs – were independent from our data and derived from the Automated Anatomical Labeling (AAL) atlas (31). A total of six ROIs were defined. Thereafter, functional connectivity of the different seeds was determined in each single subject, based on a bivariate correlation method.
The images resulting from the first level analysis provided the seed-to-voxel connectivity maps for each source, for each subject and for the rest condition. We used these data to compute the Pearson’s correlation maximum coefficient value between the time course of each of the four DMN areas and that of bilateral insula for each subject.
The second level connectivity analysis used a general linear model analysis to investigate within-group and between-group (HV vs. MI) combinations of effects, using the DMN areas and insulae as contrasts and the subjects as regressors.
We also calculated t-scores (uncorrected p-values) between each pair of seed areas (ROI-to-ROI analysis).
Statistical analysis
Group differences for demographic data were assessed using ANOVA. To explore the within-group connectivity patterns, a one-sided one sample t-test was performed on the individual normalized connectivity maps in each group (HV and MI).
The DMN functional connectivity maps in each group were obtained using the independent component analysis (ICA) analysis toolbox in GIFT v4.0 (32) without correction at a significance level of 0.001 and a minimum cluster size of five voxels. These functional connectivity maps were used to visually represent the spatial characteristics of the DMN areas in Figure 1.
Spatial characteristics of the areas of the Default Mode Network included in the seed-based resting state analysis for MI patients (in blue) and HV (in yellow/red): Posterior cingulate cortex (PCC), inferior parietal lobule (IPL), medial prefrontal cortex (MPFC). We also collected seed-based resting state data from the insula (INS) bilaterally.
The difference between groups of connectivity between the four DMN areas and bilateral insula was evaluated using an unpaired two sample t-test with an uncorrected threshold of p < 0.05.
We used Pearson’s test to search for correlations between each of the four DMN ROIs and the insula bilaterally in each subject, and to relate these values (separately for right and left insula) to individual VAS scores, as well as other clinical variables. Correlations with a p value ≤ 0.05 were considered significant.
Results
All 32 subjects completed the scanning sessions. Demographic variables are summarized in Table 1. Patients reported no premonitory symptoms prior to the headache phase. Mean headache intensity immediately before the scans was 7.7 ± 0.9 (range: 6–10) on the 10 cm visual analogue scale. None of the patients withdrew from the study due to migraine symptoms or to MRI-induced discomfort. Structural brain MRIs were normal in all participants.
Seed-based resting state fMRI
To visualize the strength of the DMN connectivity, we created group-level t-maps for healthy volunteers (HV) and ictal migraine patients (MI).
The two sample t-test disclosed a significantly greater correlation in MI between MPFC and PCC (t = 2.33, p = 0.013) than in HV (Figure 1). MI patients also exhibited stronger connectivity than HV between MPFC or right IPL and right (t = 1.71, p = 0.049; t = 1.82, p = 0.039, respectively) or left insula (t = 1.91, p = 0.033; t = 1.76, p = 0.044, respectively) (Figure 2). We failed to find a significant correlation between PCC and insula bilaterally.
Default Mode Network (DMN) inter-regional resting state-functional connectivity in migraine patients during an attack (MI) and healthy volunteers (HV). Left panel: Illustrative image from an MI patient showing activated DMN areas overlaid on the “ch2better” brain template image for visualization. Right panel: Histograms of correlation coefficients (r) between middle prefrontal cortex (MPFC) and posterior cingulate cortex (PCC) (upper), between MPFC and right and left insula (INS) (middle), and between right inferior parietal lobule (IPL) and right and left INS (lower). Correlations (mean + standard errors) are significantly greater in MI than in HV for all connectivities.
Clinical correlations
In MI, the strength of MPFC-to-INS connectivity was negatively correlated with subjective perception of migraine headache severity as measured by the VAS score (right F = 5.59, p = 0.038, r = −0.585, R-sq(adj) = 27.7%; left F = 5.43, p = 0.040 r = −0.575, R-sq(adj) = 26.95%) (Figure 3). We found no significant correlation between other clinical features of migraine and MPFC-to-INS connectivity or any other ROI-to-ROI connectivity.
In patients, MPFC–to-insula connectivity is negatively correlated with the intensity of the migraine headache as measured on a visual analogue scale.
Discussion
The main findings of our study can be summarized as follows: 1) compared with HV, migraine patients had stronger intrinsic connectivity between MPFC and PCC regions of the DMN; 2) in MI, seed-based analyses disclosed increased connectivity of MPFC and right IPL DMN areas to bilateral insula; and 3) contrary to our initial hypothesis, the strength of MPFC-to-insula connectivity was negatively correlated with the severity of perceived migraine headache. Altogether, our findings identify the insula bilaterally and MPFC as key areas of maladaptive intrinsic connectivity during migraine attacks.
Interregional DMN connectivity
To the best of our knowledge, abnormal connectivity between areas of DMN integrity have not been reported before during migraine attacks. Disruption of intra-DMN connectivity was, however, found in migraine patients who were scanned between attacks (15–20). Interictal spontaneous activity was significantly decreased in the major regions of the DMN (15–18), including the PCC, IPL and MPFC, and the latter activity was negatively correlated with the duration of the migraine disorder (16,17). In an independent component RS-fMRI analysis we recently found, in migraine between attacks, that the DMN was less connected with a network composed of visuo-spatial and medial visual cortical areas, and that the strength of the latter network was in turn associated with microstructural thalamic abnormalities (19). Here, with a seed-based analysis, migraine patients scanned during an attack had a significant increase in connectivity strength between MPFC and PCC, ROIs belonging to the DMN, as compared with HV, instead of a decrease as previously reported interictally.
Among the set of DMN regions, the midline MPFC and of PCC have distinct roles (21,33). The PCC is a tonically active brain region that may continuously gather extrinsic, and possibly intrinsic information (21). The PCC responds to both emotional and non-emotional stimuli (34). A large body of data assigns to the MPFC a role in the integrative processing of emotional and cognitive processes by incorporating emotional biasing signals or markers into decision-making processes (21). It appears that MPFC and PCC normally act synergistically. In fact, Kucyi et al. (35) showed that when healthy subjects’ minds wander away from incoming nociceptive stimulations, there is increased connectivity between the MPFC and PCC regions of the DMN, as well as areas belonging to the salience network (SN) such as the insula. Concordantly, we show that, in patients who are at rest during the headache phase, connectivity between PCC and MPFC is increased, and that MPFC is highly connected with bilateral insula. This might be the consequence of a spontaneous tendency of some migraineurs to try to mind-wander away from their pain because of the absence of other immediately available coping strategies.
DMN-to-insula connectivity
In animal studies, the insular cortex, anchored with the SN, is reciprocally connected to multiple sensory, motor, limbic, and association areas of the brain (36). The insular cortex and areas within the DMN share significant topographic reciprocal connectivity and form an anatomically tightly coupled network ideally placed to integrate information from several brain regions (36). This emphasizes that the insula has numerous functions; for instance, it takes part in processing subjective emotional experiences, in sensorimotor integration, in goal-directed cognition, and in autonomic response. Moreover, the insula may channel inputs into the SN and regulate spontaneous switching between the DMN and the executive control network (ECN) to optimize performance (37). Previous neuroimaging studies during migraine attacks have repeatedly provided evidence of activation in bilateral insula, which was considered as a non-specific response to pain (2,4,9,10,38). Between attacks, migraineurs also had abnormal insular cortex connectivity (39,40) and habituation in response to repeated sensory stimuli (41,42). In the present study, we show that migraine patients during an attack have greater bilateral insula connectivity with MPFC and right IPL areas of the DMN than healthy volunteers, in the absence of a specific task. Taken together with our previous results obtained in the same group of subjects showing reduced functional connectivity between the ECN and the attentive networks (27), the present findings suggest that, during a migraine attack, the control by the insula of the spontaneous switching between the DMN and the ECN to optimize performance is abnormal, favouring activation of the former over the latter. Additionally, since the prefrontal cortex is deemed to mediate attenuation of pain via cognitive control mechanisms (43), increased MPFC-to-insula connectivity during the migraine attack could be a physiological mechanism to cognitively attenuate pain perception.
DMN-to-INS correlation with severity of migraine headache
Another striking finding in our study is the correlation between DMN-to-insula connectivity and severity of the migraine headache. During migraine attacks, individual r-scores of MPFC-to-bilateral insula connectivity correlated negatively with subjective perception of headache intensity. Remarkably, this finding is at variance with previous studies on the association between clinical pain and DMN connectivity in chronic pain disorders like fibromyalgia and chronic low back pain. As a matter of fact, in these studies greater clinical pain was associated with greater MPFC connectivity with the insula (22–24) and not lesser, as shown here. This raises the intriguing possibility that the anti-correlation between MPFC-to-insula connectivity and head pain severity is a distinct feature of episodic migraine attacks. It may also be related to the finding, in patients between attacks, that connectivity between insula and the DMN is positively correlated with the usual duration of migraine attacks (15). Moreover, when taken separately, the regional metabolism of the insular cortex (44) and the regional homogeneity within the MPFC (17) are significantly anti-correlated with disease duration, and insular metabolism is anti-correlated with lifetime headache frequency (44). Further studies are needed to verify whether the aberrant ictal/interictal connectivity patterns in migraine are permanent functional abnormalities of a progressive brain disorder or, more likely, plastic changes of brain functions relative to the time of an attack, as previously clearly shown with clinical neurophysiology methods (45).
Interestingly, besides pain-induced DMN and insula activation, the above-mentioned study by Kucyi et al. (35) has shown stronger dynamic functional connectivity between MPFC and periaqueductal gray (PAG), an opiate-rich region mediating pain suppression, suggesting that sustained DMN activity during pain, for instance during head pain as in our study, may be a physiological strategy of the brain to try to reduce or limit pain (46). In this context, further studies are worthwhile to determine whether the microstructural changes in the PAG and its abnormal connectivity with the prefrontal cortex (47,48), particularly with worsening of attack frequency (49) and intensity (50), may be primary dysfunctions or secondary maladaptive consequences of the abnormal connectivity between networks during the migraine cycle.
It also remains to be determined if the abnormal DMN-to-INS connectivity could be a consequence of a ‘thalamocortical dysrhythmia’, a model for the cyclic functional brain abnormalities found in migraine (51,52). For instance, Wang et al. (53) found that inflow/outflow signals between the posterior thalamic nuclei and the limbic system, MPFC, and temporo-parieto-occipital areas were decreased in relation to the severity of migraine clinical features. This might support the hypothesis that a disrupted thalamic control of brain networks, such as that found by us previously in patients during spontaneous attacks (27), could underlie the distinct resting state connectivity patterns found in the present study.
We acknowledge some limitation of our study. Firstly, cross-sectional fMRI studies limit inferences about causality. Secondly, as all patients were ictal at the time of scanning, the results provide no information about brain functions between attacks. It is nevertheless striking that disruptions of the DMN and insula connectivity have been consistently observed in episodic migraine also between attacks. Future work using the same methods is necessary to follow patients longitudinally, i.e. during attacks as well as at different time points during the pain-free periods. Since our study sample is rather small, larger studies are also needed to replicate our findings and to include high-frequency and chronic migraine patients.
Conclusions
In patients recorded during a full-blown migraine attack, we report increased connectivity between areas belonging to the DMN as well between the latter and bilateral insula. Moreover, we show that the strength of the DMN-to-INS connectivity is anti-correlated with intensity of the migraine headache, which contrasts with the opposite finding in chronic pain disorders. Altogether, these results support, on the one hand, resting DMN-to-INS connectivity as a state-specific neuroimaging marker for the encoding of the clinical characteristics of pain, and, on the other hand, significantly extend these data by providing evidence that a distinct connectivity pattern may encode intensity of acute migraine headache. Hence, comparative studies of other acute painful conditions are necessary to verify the specificity of these findings for migraine.
Article highlights
How the brain encodes the severity of acute migraine headache is unknown. We found that migraineurs during attacks had stronger intrinsic default-mode network (DMN) connectivity and that the strength of DMN-to-insula connectivity was anticorrelated with the severity of headache, which is at variance with publications on chronic pain disorders.
Footnotes
Acknowledgment
The contribution of the G.B. Bietti Foundation in this paper was supported by the Ministry of Health and Fondazione Roma.
Declaration of conflicting interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Fondazione Roma, Italian Ministry of Health.
