Ictal and interictal brain activation in episodic migraine: Neural basis for extent of allodynia

Abstract

In some patients, migraine attacks are associated with symptoms of allodynia which can be localized (cephalic) or generalized (extracephalic). Using functional neuroimaging and cutaneous thermal stimulation, we aimed to investigate the differences in brain activation of patients with episodic migraine (n = 19) based on their allodynic status defined by changes between ictal and interictal pain tolerance threshold for each subject at the time of imaging. In this prospective imaging study, differences were found in brain activity between the ictal and interictal visits in the brainstem/pons, thalamus, insula, cerebellum and cingulate cortex. Significant differences were also observed in the pattern of activation along the trigeminal pathway to noxious heat stimuli in no allodynia vs. generalized allodynia in the thalamus and the trigeminal nucleus but there were no activation differences in the trigeminal ganglion. The functional magnetic resonance imaging (fMRI) findings provide direct evidence for the view that in migraine patients who are allodynic during the ictal phase of their attacks, the spinal trigeminal nucleus and posterior thalamus become hyper-responsive (sensitized)–to the extent that they mediate cephalic and extracephalic allodynia, respectively. In addition, descending analgesic systems seem as “switched off” in generalized allodynia.

About the Authors:

Nasim Maleki

Contributed equally to this work with: Nasim Maleki, Edina Szabo

Roles Data curation, Formal analysis, Writing – original draft, Writing – review & editing

Current address: Psychiatric Neuroimaging Division, Department of Psychiatry, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States of America

Affiliations Center for Pain and the Brain (PAIN Research Group), Boston, MA, United States of America, Department of Anesthesiology, Critical Care and Pain Medicine, Boston Children’s Hospital, Harvard Medical School, Boston, MA, United States of America

Edina Szabo

Contributed equally to this work with: Nasim Maleki, Edina Szabo

Roles Data curation, Formal analysis, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliations Center for Pain and the Brain (PAIN Research Group), Boston, MA, United States of America, Department of Anesthesiology, Critical Care and Pain Medicine, Boston Children’s Hospital, Harvard Medical School, Boston, MA, United States of America

ORCID logo https://orcid.org/0000-0003-4087-5373

Lino Becerra

Roles Data curation, Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing

Affiliations Center for Pain and the Brain (PAIN Research Group), Boston, MA, United States of America, Department of Anesthesiology, Critical Care and Pain Medicine, Boston Children’s Hospital, Harvard Medical School, Boston, MA, United States of America

Eric Moulton

Roles Methodology, Writing – review & editing

Affiliations Center for Pain and the Brain (PAIN Research Group), Boston, MA, United States of America, Department of Anesthesiology, Critical Care and Pain Medicine, Boston Children’s Hospital, Harvard Medical School, Boston, MA, United States of America, Department of Ophthalmology, Boston Children's Hospital, Harvard Medical School, Boston, MA, United States of America

Steven J. Scrivani

Roles Data curation, Writing – review & editing

Affiliations Center for Pain and the Brain (PAIN Research Group), Boston, MA, United States of America, Department of Diagnostic Sciences, Craniofacial Pain Center, Tufts University School of Dental Medicine, Boston, MA, United States of America, Department of Public Health and Community Medicine, Pain Research, Education and Policy Program, Tufts University School of Medicine, Boston, MA, United States of America

Rami Burstein

Roles Conceptualization, Data curation, Writing – original draft, Writing – review & editing

Affiliation: Department of Anesthesia and Critical Care, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA, United States of America

David Borsook

Roles Conceptualization, Data curation, Visualization, Writing – original draft, Writing – review & editing

Affiliations Center for Pain and the Brain (PAIN Research Group), Boston, MA, United States of America, Department of Anesthesiology, Critical Care and Pain Medicine, Boston Children’s Hospital, Harvard Medical School, Boston, MA, United States of America, Department of Psychiatry, Brain Imaging Center, McLean Hospital, Harvard Medical School, Belmont, MA, United States of America

Abstract

In some patients, migraine attacks are associated with symptoms of allodynia which can be localized (cephalic) or generalized (extracephalic). Using functional neuroimaging and cutaneous thermal stimulation, we aimed to investigate the differences in brain activation of patients with episodic migraine (n = 19) based on their allodynic status defined by changes between ictal and interictal pain tolerance threshold for each subject at the time of imaging. In this prospective imaging study, differences were found in brain activity between the ictal and interictal visits in the brainstem/pons, thalamus, insula, cerebellum and cingulate cortex. Significant differences were also observed in the pattern of activation along the trigeminal pathway to noxious heat stimuli in no allodynia vs. generalized allodynia in the thalamus and the trigeminal nucleus but there were no activation differences in the trigeminal ganglion. The functional magnetic resonance imaging (fMRI) findings provide direct evidence for the view that in migraine patients who are allodynic during the ictal phase of their attacks, the spinal trigeminal nucleus and posterior thalamus become hyper-responsive (sensitized)–to the extent that they mediate cephalic and extracephalic allodynia, respectively. In addition, descending analgesic systems seem as “switched off” in generalized allodynia.

Figures

Fig 6

Table 2

Fig 1

Table 1

Fig 2

Fig 3

Fig 4

Fig 5

Fig 6

Table 2

Fig 1

Table 1

Fig 2

Citation: Maleki N, Szabo E, Becerra L, Moulton E, Scrivani SJ, Burstein R, et al. (2021) Ictal and interictal brain activation in episodic migraine: Neural basis for extent of allodynia. PLoS ONE 16(1): e0244320. https://doi.org/10.1371/journal.pone.0244320

Editor: Andrea Antal, University Medical Center Goettingen, GERMANY

Received: May 21, 2020; Accepted: December 7, 2020; Published: January 4, 2021

Copyright: © 2021 Maleki et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: At the time the imaging datasets were generated, there was no paragraph in the participant consent forms with regard to open access data availability (the study was approved by the McLean Hospital Institutional Review Board). Since participants were not able to give or refuse their assent for public data sharing, the authors are not allowed to share MRI acquisitions publicly. However, data access requests can be made to to the corresponding author ([email protected]) or the Department of Anesthesiology, Critical Care and Pain Medicine at Boston Children’s Hospital ([email protected]).

Funding: The work was supported by NIH Grants NS073997 (NINDS) and K24NS064050 (NINDS) to David Borsook.

Competing interests: The authors have declared that no competing interests exist.

1. Introduction

During a migraine attack, most patients develop cutaneous allodynia, defined as pain perception evoked by normally non-noxious stimuli [1, 2]. That is, patients with migraine may experience pain from common daily activities, such as wearing glasses, or taking a shower [2, 3]. The pain in migraine is thought to contain an early activation and sensitization of first-order neurons in the trigeminal ganglion characterized by limited or no skin sensitivity to the head (no allodynia), causing activation in second-order neurons in the spinal trigeminal nucleus. Whereas sensitization of the first-order neurons explains the throbbing perception of migraine headaches and their aggravation by physical activity, second-order neurons in the spinal trigeminal nucleus appear to be responsible for mediating the development of allodynia on the ipsilateral side of head (localized or cephalic allodynia), and subsequent sensitization of third-order thalamic trigeminovascular neurons (particularly in the pulvinar region of the thalamus) has been related to the development of allodynia on the contralateral side of the head (localized or cephalic allodynia) and body (generalized or extracephalic allodynia) [4].

The associated manifestations of allodynia have clinical implications including (1) period of time for effective treatment in episodic migraine (in patients with allodynia, the effectiveness of triptans is much higher when they are administered early in the attack) [1, 5], (2) differences in triggers in allodynic vs. non-allodynic migraine (patients with allodynia have more triggers than those with no allodynia) [6], and (3) predicting migraine progression (allodynia is a risk factor for migraine chronification) [7].

Most imaging research has been done between migraine attacks (i.e., interictally) [8], there are fewer studies focusing on patients during an acute migraine attack (i.e., ictally) [9]. Previous findings have shown changes in the trigeminal system in healthy subjects after noxious stimulation [10] and in patients with migraine [11]. Nevertheless, there are no studies investigating the differences between allodynic states in the pathway. Furthermore, generalized allodynia may require other brain systems (i.e., cognitive, sensory) which could provide a basis for differences in underlying brain mechanisms that may relate to treatment resistance or migraine progression [7, 12]. Altered brain systems in migraine hyperexcitability have been well-documented but the issue of how brain responses may relate to the extent of allodynia has not been assessed in imaging studies before.

Our aim was to investigate brain activation in allodynic and non-allodynic migraine patients during their ictal vs. interictal state to assess cortical and subcortical differences underlying the phenomenon of migraine-associated allodynia as defined according to changes between ictal and interictal pain tolerance threshold in response to thermal stimuli using quantitative sensory testing (QST). We focused on two questions: (1) Are there any differences in response to noxious thermal stimulation in the migraine brain between ictal vs. interictal states along the trigeminovascular pathway using functional magnetic resonance imaging (fMRI)?; and (2) How are the brain response changes manifested in the generalized allodynic state different from those in the non-allodynic state? The observations and insights gathered from the brain changes in allodynia may help understand if the progression from no allodynia to generalized allodynia engages brain areas involved in the sensory, emotional and other aspects of migraine and how each is involved.

2. Methods

2.1. Patient recruitment and consent

Subject recruitment and the experiment protocol were approved by the McLean Hospital Institutional Review Board and met the scientific and ethical guidelines for human pain research of the Helsinki Accord (http://ohsr.od.nih.gov/guidelines/helsinki.html) and the International Association for the Study of Pain. Subjects were recruited through participating physicians and online advertisements. After a telephone interview, participants were scheduled for a screening session during which all aspects of the study were described (including study procedures, time commitment, and the types of pain stimuli to be used). To ensure compliance with the inclusion and exclusion criteria prior to subject enrollment, a study physician performed a full medical examination, and a detailed medical history was taken from the patients (including information on the age of migraine onset, duration and frequency of the migraine attack, presence of aura, type of medication the patients used, and their history of allodynic status). Written informed consent was received from all participants by the study staff after answering any study-related questions.

2.2. Inclusion/exclusion criteria

The inclusion criteria were as follows: (i) episodic migraine as defined by the International Classification for Headache Disorders (ICHD-2) [13]; (ii) Beck Depression Inventory II (BDI-II; [14]) scores ≤ 25; (iii) had episodic migraine for at least three years or longer. The exclusion criteria were as follows: (i) other significant medical condition (except migraine); (ii) depression measured by BDI-II [14] with a cut-off score of ≥ 25 (implicating moderate to severe depression); (iii) on preventatives; (iv) positive on pregnancy or drug screening tests; (v) having MRI contraindications (such as having a heart pacemaker or unremovable body piercing etc.). Participants who took medications for their migraine headaches were asked to discontinue their pain medication for one dosing period before their scan sessions. They were allowed to continue taking these medications after the scan visit. Subjects were also asked to refrain from eating or drinking caffeinated beverages for 12 hours before their scan and from eating in general up to 6 hours before the scan session.

2.3. Study visits

The study consisted of two visits, one during the interictal period (baseline scan) and one during an ictal event (migraine scan). The initial interictal visit was scheduled when the patients had been migraine free for 72 hours and did not have any symptoms or signs of developing a migraine attack. For the ictal visit, arrangements were made so that the research team was available between 6am to 8pm on weekdays and weekends to scan the patients within 2.5 to 5 hours of the migraine attack onset.

2.4. Quantitative sensory testing (QST)

In order to determine individual sensitivity and pain thresholds to thermal stimuli, each subject participated in QST prior to each scan session. Thermal pain thresholds were measured on the face and dorsum of the hand, and the side of the body was determined based on which side the pain is felt more often during a migraine attack. The same area on the back of the hand was tested for all patients. Participants received between 3–5 thermal stimuli and the average of these values +1˚C was the threshold temperature that the participant would receive during the scan.

2.5. Procedure

The experimental paradigm included two visits to the McLean Hospital Brain Imaging Center. The scans during both visits consisted of a series of anatomical scans followed by a functional imaging scan. For both functional scans, a painful heat stimulus (threshold +1°C) was applied to the hand. After the ictal scan, patients were allowed to take their pain medications.

2.6. Stimulation paradigm

Patients rated their pain while lying in the MRI scanner during each functional scan using a rating dial. They were instructed to assess the pain intensity evoked from each stimulus using a visual analog scale (VAS; 0 = no pain, 10 = maximum pain).

2.7. Thermal stimulation

A Federal Drug Administration-approved Thermal Sensory Analyzer (TSA-II, Medoc Advanced Medical Systems, Ramat Yishai, Israel) was used to deliver a thermal stimulus through a probe that was adapted to rest on the hand. The probe is 1.6 x 1.6 cm, or about one-half the size of the thumb. For warm heat, stimuli were given at 41˚C on the hand (a standard temperature probe that is commonly used as a low-innocuous stimulus in studies of episodic migraine; e.g., [15–18]). For painful heat, stimuli were given at +1˚C above the subject’s individualized pain threshold value (as determined by QST) on the hand. Stimuli were repeated three times, each for a length of 25 seconds and separated by 30 seconds of no stimuli.

2.8. Imaging

2.8.1. Image acquisition.

Imaging was performed in a 3.0 T Siemens (Erlangen, Germany) Trio scanner using a circularly polarized (CP) coil. A Gradient Echo (GE) echo planar imaging (EPI) sequence (TR = 2500 ms, TE = 30 ms, resolution = 3.5 x 3.5 x 3.5 mm3, matrix = 64 x 64, 74 volumes, 41 slices) was used to acquire the functional images. High resolution, T1-weighted structural data were also collected from each patient using a 3D MPRAGE pulse sequence (TR = 2100 ms, TE = 2.74 ms, flip angle = 12°, 128 sagittal slices, resolution = 1.33 × 1.0 × 1.0 mm3).

2.8.2. Image analysis.

Functional imaging data processing and analyses were carried out using FMRIB Software Library (FSL) (www.fmrib.ax.ac.uk/fsl), version 4.1.9. The preprocessing steps included skull stripping using the Brain Extraction Tool (BET) with bias field correction and neck removal. The volumes were motion corrected using FMRIB's Linear Image Registration Tool (MCFLIRT) and B0 unwarping was performed using FMRIB's Utility for Geometrically Unwarping EPIs (FUGUE) to correct for some field inhomogeneities. The data were spatially smoothed with a Gaussian kernel of 5 mm full width at half maximum (FWHM) and a high-pass temporal filter with cutoff of 60s was applied. First-level analysis on functional imaging data for each subject was carried out using FMRI Expert Analysis Tool (FEAT) Version 5.98. The recorded temperature traces for each participant were entered as the explanatory variables (EVs) for heat stimuli. EVs were convolved with a gamma-variate hemodynamic response function. Functional volumes were spatially normalized to the Montreal Neurological Institute (MNI152) template for group analysis, and individual data of right-sided migraine patients (who were stimulated on the right side) were flipped along the y-axis so that they correspond with the patients with left-sided migraine. Thus, in the group results the left side is the ipsilateral stimulation side, while the right side is the contralateral stimulation side. Patients with no dominant migraine side were also stimulated on the left. Contrasts were defined for the non-noxious and noxious thermal stimuli for each subject, and in the group level analysis contrasts were generated to compare allodynic patients to non-allodynic patients, and to investigate patients during their ictal and interictal state. Higher level group contrast analysis was carried out to calculate the intragroup and intergroup comparison activation maps using FMRIB’s Local Analysis of Mixed Effects (FLAME1) [19]. All statistical parametric maps were thresholded using a Gaussian Mixture Model (GMM) technique [20] which is a generalized false discovery rate (FDR) correction for multiple comparisons (p < 0.05).

2.8.3. Region of interest (ROI) analyses.

In addition to the whole-brain approach, two region of interest (ROI) analyses were conducted: (1) ROIs were determined for the trigeminal ganglion and trigeminal nucleus based upon previously implicated areas in the trigeminal nociceptive pathway [21–23], and the thalamus in the standard MNI space using the MNI structural atlas (available in FSL) [24]; (2) ROIs were defined in the thalamus based on a probabilistic atlas of thalamic sub-regions segmented according to their white-matter connectivity to cortical areas [25, 26]. More specifically, the Oxford thalamic connectivity atlas was used (included in FSL), which reports probabilistic connectivity values based on diffusion tractography. This comprises bilateral thalamic subregions that have been identified to show specific anatomical connections to different cortices. For both ROI analyses, the group average contrast of parameter estimates (COPE) values were converted to percentage (%) signal change in ictal vs. interictal scans for noxious stimulation using FSL’s Featquery for which the contrast maps were thresholded according to the GMM determined thresholds.

3. Results

3.1. Patients' characteristics

Fig 1 shows the flowchart of patient inclusions in the study. From the 51 episodic migraine patients, 26 patients completed both scan visits for this study. However, six patients were excluded from the analysis due to the long time separation between the interictal scan and the scan obtained during the migraine attack (> 1 year). Twenty patients completed both visits within 9 months (13 females, 7 males; 44.6 ± 11.7 years old, age range: 24–61 years). One female subject was excluded later from the imaging analysis due to the concerns of residual effects of a beta-blocker. Nineteen episodic migraine patients (12 females) were included in the final imaging analysis. A detailed description of demographic and migraine attack characteristics is shown in Table 1. There were no significant differences between males and females in age and migraine characteristics, but females used more preventative medications than males. In addition, age, sex, migraine characteristics and medication history were not significantly different between the allodynia groups defined by QST testing (see below). All of the participants had migraine for at least more than four years and experienced five to six migraine attacks on average each month with a range of one to twelve migraine attacks per month.

[Figure omitted. See PDF.]

Fig 1. Flowchart showing the sequence of patient inclusion.

https://doi.org/10.1371/journal.pone.0244320.g001

[Figure omitted. See PDF.]

Table 1. Demographic and clinical characteristics of the migraine patients.

https://doi.org/10.1371/journal.pone.0244320.t001

3.2. Timing of the scan visits

Time elapsed between the first and second visit ranged from 1 month to 9 months with an average of 3 months (SD = 2.5). Since we were interested in assessing the differences in brain activation among migraine patients based on their allodynic status defined by changes between ictal and interictal pain tolerance threshold at the time of imaging, migraine scan was delayed relative to the time of onset for an average of 4 hours [1]. Once the patients were at the imaging center at McLean Hospital the scans happened shortly after they were admitted within 34.5 minutes on average (SD = 19.9). During this time QST was performed to measure the pain thresholds from which the allodynic status was determined (see below). Patients were experiencing comparable levels of pain (M = 6, SD = 1, on a scale of 0 (no pain) to 10 (maximum imaginable pain)).

3.3. Medication history

Subjects were on a wide range of medications, including over the counter medicines (e.g., Advil, Excedrin) and prescribed medications (e.g., Imitrex, Maxalt). S1 Table in S1 File lists the medication history for each participant in the study. While a number of patients had a history of being on preventative medication, none of them were currently on such medication (except for one subject that was later excluded from the analysis). Subjects abstained from taking their migraine medications for one dosing interval prior to their interictal scan session. They also had to refrain from taking their medication on the day of their ictal visit until after the scan was completed. Medication did not change between the first and second visit.

3.4. Pain tolerance threshold and allodynia state

Pain tolerance thresholds were determined on the hand and the face for all of the patients during both interictal (baseline (B)) and ictal (migraine (M)) visits using QST. The results for the allodynia status are presented in Fig 2. Allodynia state was determined based on the changes in the pain tolerance thresholds (>3°C) between face and hand when stimulated by the heat [1, 27, 28]. In total there were seven patients who experienced no allodynia (Face: 46.2°C (B) vs. 47.3°C (M), Hand: 45.6°C (B) vs. 47.4°C (M)), four patients who were categorized as having localized allodynia (Face: 46.7°C (B) vs. 42.1°C (M), Hand: 45.7°C (B) vs. 46.5°C (M)) and the remaining eight patients experienced generalized allodynia (Face: 47.3°C (B) vs. 41.7°C (M), Hand: 47.1°C (B) vs. 43.2°C (M)).

[Figure omitted. See PDF.]

Fig 2. Change in the pain tolerance thresholds.

Allodynia state was determined according to the changes in the pain tolerance thresholds between interictal and ictal visits on face and hand. If the pain tolerance change was >3°C on both face and hand the patient was categorized as having generalized allodynia. Localized allodynia was defined if >3°C change was observed on face but not hand. If there was no change or if the change was < 3°C on face and hand that was considered no allodynia. In total there were 7 patients who experienced generalized allodynia and 8 who were categorized as having no allodynia with the remaining 4 experiencing localized allodynia.

https://doi.org/10.1371/journal.pone.0244320.g002

3.5. Visual analog scale (VAS) rating of pain intensity

The noxious heat stimulus (threshold +1°C) applied to the dorsum of the hand evoked on average a pain intensity of 6.1 ± 3 (SD) during the interictal scans and 7.4 ± 2.6 (SD) during the ictal scans on the 0–10 scale. In spite of the fact that participants were pain-free in the baseline (B) visit and experiencing head pain in the migraine (M) visit, there was no significant difference between the VAS ratings between the ictal vs. interictal visits (t(18) = 1.84, p = 0.1). The VAS data were compared within the allodynia groups as well and no significant differences were found between the ictal vs. interictal VAS scores in the allodynic groups: no allodynia (Face: 6.6 (B) vs. 6.3 (M), Hand: 5.4 (B), vs. 7.1 (M)), localized allodynia (Face: 4.6 (B) vs. 6 (M), Hand: 4.9 (B) vs. 6(M)), generalized allodynia (Face: 7.34 (B) vs. 7.1 (M), Hand: 7.9 (B) vs. 7.2 (M)).

3.6. Imaging results

Imaging results are presented in Figs 3 to 6. Group average findings for brain activation in response to noxious heat (threshold +1°C) applied to the dorsum of the hand during the interictal and ictal scans are presented in S2A and S2B Tables in S1 File.

[Figure omitted. See PDF.]

Fig 3. Brain activation for ictal vs. interictal state.

Contrast analysis of the ictal visit vs. interictal visit in response to noxious thermal stimulation on the left hand is presented. Hot colors (yellow to red colors) indicate stronger activation during the ictal state. The results are based on longitudinal data of 19 migraine patients, and non-noxious heat vs. noxious heat contrast maps. Statistical maps were thresholded and corrected for multiple comparisons using Gaussian mixture modeling (GMM) approach. Coordinates are in Montreal Neurological Institute (MNI) space.

https://doi.org/10.1371/journal.pone.0244320.g003

[Figure omitted. See PDF.]

Fig 4. Brain activation in generalized allodynia vs. no allodynia.

Contrast analysis of the ictal vs. interictal visits in response to noxious thermal stimulation on the left hand is presented for patients with generalized allodynia vs. patients with no allodynia. Hot colors (yellow to red colors) indicate stronger activation during the ictal state. The results are based on longitudinal data of 7 patients with generalized allodynia and 8 patients with no allodynia. Statistical maps were thresholded and corrected for multiple comparisons using Gaussian mixture modeling (GMM) approach. Coordinates are in Montreal Neurological Institute (MNI) space.

https://doi.org/10.1371/journal.pone.0244320.g004

[Figure omitted. See PDF.]

Fig 5. Signal change differences along the trigeminal pathway based on Region of Interest (ROI) analysis.

Signal change values were extracted using Featquery. For those patients who received thermal stimulation on the right side, the brains were flipped along the y-axis so that all images were in accordance regarding the site of the migraine (and the stimulation side). The left side is the ipsilateral stimulation side and the right is the contralateral stimulation side. The group average results for percentage signal change between generalized and no allodynia groups are shown for the left thalamus and left trigeminal nucleus ipsilateral to the side of stimulation. For the trigeminal ganglion, the signal change is reported for the contralateral region (right side). Error bars represent standard deviations. The * symbol indicates group differences at a significant level of p < 0.05.

https://doi.org/10.1371/journal.pone.0244320.g005

[Figure omitted. See PDF.]

Fig 6. Region of Interest (ROI) based comparison results for percentage (%) signal change in various thalamic regions according to the cortical connectivity in ictal vs. interictal scans.

Signal change values were extracted using Featquery. Thalamic sub-regions (indicated by different colors) were segmented based on the white-matter connectivity between the thalamus and cortex. The label for each color is listed in the x-axis. Bars represent changes in the percentage of the average signal change between ictal vs. interictal activation in response to noxious heat applied to the dorsum of the hand for patients with generalized allodynia and with no allodynia. The main differences are in the thalamic regions that have connections with temporal and prefrontal cortical regions with an increased signal change in generalzied allodynia patients. In contrast, thalamic regions with connectivity to the posterior parietal regions showed stronger signal change in the no allodynia group. Error bars represent standard deviations. The * symbol indicates group differences at a significant level of p < 0.05.

https://doi.org/10.1371/journal.pone.0244320.g006

3.6.1. Interictal (baseline) state.

Ipsilateral activation in response to noxious heat stimulation was observed in the thalamus, middle cingulate, precentral, and superior frontal gyrus. Contralateral activation was found in the anterior insula, caudate, palladium, temporal pole, middle temporal, and postcentral gyrus. And bilateral activation was seen in response to stimulation in the supramarginal, inferior parietal, inferior operculum gyrus and cerebellum, although activation in the cerebellum was mainly contralateral.

3.6.2. Ictal (migraine) state.

Ipsilateral activation in response to noxious heat stimulation was observed in the middle cingulate, precentral, and superior frontal gyrus. Contralateral activation was seen in the pons, anterior insula, caudate, palladium, temporal pole, middle temporal, and postcentral gyrus, as well as in the cerebellum. Bilateral activation was observed in thalamus, rolandic operculum in the occipital cortex, supramarginal, inferior parietal and inferior operculum gyrus.

3.6.3. Ictal vs. interictal states.

Contrast analysis results for brain activation in ictal vs. interictal scans in response to noxious heat (threshold +1°C) applied to the hand are shown in Fig 3 and the corresponding activations are shown in Table 2A and 2B. The main differences were significantly stronger ipsilateral activation during the migraine scan in the anterior and posterior insula, supramarginal gyrus and cerebellum, as well as stronger contralateral activation in the anterior and middle cingulate, trigeminal nucleus and brainstem (pons). Significantly stronger bilateral activation during the migraine attack was found in the occipital cortex, and posterior cingulate. There were also brain regions that showed increased activation during the interictal scan compared to the ictal scan. The main differences were in the frontal cortical regions and temporal regions. There was an increased activation in the ipsilateral precentral gyrus, postcentral gyrus, as well as putamen. Contralateral stronger responses were found in the inferior, superior and middle temporal gyrus.

[Figure omitted. See PDF.]

Table 2. a. Whole brain contrast results in ictal vs. interictal scans (ictal > interictal). b. Whole brain contrast results in ictal vs. interictal scans (interictal > ictal).

https://doi.org/10.1371/journal.pone.0244320.t002

3.6.4. Patients with generalized allodynia vs. no allodynia.

Additional analyses were performed to compare brain activation in ictal vs. interictal scans in response to noxious heat in patients who showed changes in pain threshold in both face and hand (generalized allodynia, n = 7), and in the patients who showed no change in the pain thresholds between the two visits (no allodynia, n = 8). The results are presented in Fig 4. The main differences between the two groups were in the activation of the contralateral brainstem (pons), insula and thalamus, which were present only in the generalized allodynia group and not in the no allodynia group. On the other hand, the no allodynia group showed increased activity during a migraine attack in the posterior cingulate and periaqueductal gray (PAG) that was not observed in the generalized allodynia group.

In Fig 5 ROI-based comparison results for % signal change between the interictal and ictal scans for noxious stimulation of the hand are shown. The activations are presented along the trigeminal nociceptive pathway in patients with generalzied vs. no allodynia for the left thalamus and left trigeminal nucleus, ipsilateral to the side of stimulation and for the trigeminal ganglion on the contralateral side. The results show a higher level of signal change, suggestive of an increased hypersensitivity, for generalized allodynia vs. no allodynia in the thalamus and the trigeminal nucleus. The level activity for the trigeminal ganglion was comparable between the two groups.

The results in Fig 6 present % signal changes in ictal vs. interictal scans for thalamic sub-regions according to their cortical connectivity. For this analysis, thalamic segmentation was performed using a probabilistic atlas of sub-thalamic regions that were generated based on white-matter connectivity between the thalamus and cortical areas [25, 26]. The bars represent changes in the percentage of the average signal change between ictal vs. interictal brain activation in response to noxious heat (threshold +1°C) for subjects with generalized allodynia and subjects with no allodynia. For both groups the biggest changes were in the thalamic region with connectivity to temporal and prefrontal regions in the cortex. The main differences were found in the thalamic regions that have connections with temporal and prefrontal cortical regions with a stronger signal change in generalzied allodynia patients. In the no allodynia group, there was a stronger signal change in the posterior parietal region.

4. Discussion

4.1. Overall findings

We designed a study to assess the brain activation in migraine patients in response to noxious thermal stimulation during both a migraine attack (ictal) and between attacks (interictal) in the context of the changes in the thermal pain tolerance threshold as an index of allodynia. Aside from the allodynic measures, there are a number of features that distinguish this study from others: (1) the migraine happened naturally and was not induced [29–32]; (2) the timing of the ictal scan relative to the onset of the attack was considerably consistent and similar among the patients [4]; and (3) the ictal scan happened at a timepoint at which allodynic symptoms would be present [1, 4, 9].

4.2. Thermal thresholds differentiate allodynic state

Differences in thermal thresholds between ictal vs. interictal visits for migraine patients were used to categorized patients’ allodynic state which allowed defining three distinct groups of no-allodynia, localized allodynia, and generalized allodynia with distinct pain tolerance profiles as shown in Fig 2. Questions about experiencing allodynic symptoms were also asked before enrolling patients to ensure compliance with the inclusion criteria. QST is sometimes described as the gold standard for detecting whether a patient with migraine has allodynia at a particular time [3, 27], although it is also prone to temporal sampling bias. Nevertheless, these patterns are consistent with our prior reports on the development of allodynia in migraine patients [1, 4, 9] and were used to classify patients as shown in Fig 1.

4.3. Brain activation differences between ictal and interictal states

As shown in Fig 3 and Table 2A and 2B we observed differences between these two states with increased activation to thermal stimuli in trigemino-thalamic pathways (brainstem, including the trigeminal nucleus, pons and thalamus) and other brain regions (including the cingulate, temporal, anterior and posterior insula cortices, and cerebellum). Activation in these regions to thermal stimuli was previously reported in interictal migraine patients by our group [8, 33–37] and others [38–40]. These areas seem to be enhanced during the ictal state vs. interictal state by thermal stimulation.

4.4. Differences between allodynic states–trigemino-thalamic pathway

Our major findings were the observed differences between the no and generalized allodynia groups including increased activations in ictal vs. interictal state of the contralateral brainstem (pons), insula and thalamus, which were observed only in the generalized allodynia group and not in the no allodynia group, and that the no allodynia group showed increased activation during a migraine attack in the posterior cingulate and periaqueductal gray (PAG) that was not observed in the generalized allodynia group.

These findings suggest that in the no allodynia group, the PAG does what it supposes to do which is inserting inhibitory effects on ascending nociceptive neurons in spinal and trigeminal dorsal horn [41, 42]. Allodynia may be spatially distributed as defined in abnormal spatial distribution of pain in fibromyalgia [43, 44]. On the other hand, larger activation in the thalamus, pons and insula of generalized allodynic patients may suggest that the development of generalized allodynia is more complicated than previously thought. The development of hyper-responsiveness in the spinal trigeminal nucleus and posterior thalamus may depend not only on the continuous arrival of pain signals from the meninges, but also on some sort of a “breakdown” in the proper function of brainstem and cortical areas that regulate thalamocortical inputs.

4.5. Differences between allodynic states and ictal vs. interictal phases–non-trigemino-thalamic pathway

As noted in the results, regions showing increased activation in ictal vs. interictal state outside the traditional trigemino-thalamic pathways included a number of cortical regions: the cingulate cortex, orbitofrontal regions, the anterior and posterior insula, the calcarine, and the supramarginal cortices. The cingulate displayed greater increases in the contrast for ictal vs. non-ictal state for the anterior (ACC), middle (MCC) and posterior (PCC). These regions reflect possible complex involvement in pain emotion [45] and other complex processes, such as outcome evaluation [46], and the ACC is activated in migraine using experimental heat stimuli [16, 17]. The MCC is putatively involved in pain/avoidance and reward approach functions amongst others [47]. The PCC is considered to be involved in functions that include arousal and attention [48]. The orbitofrontal cortices also showed a difference between the two states. This region is involved in emotional and executive functions [49] as well as reward value for taste [50] and may thus be implicated in alterations in both taste and smell in patients with migraine. The supramarginal cortex, a sensory association cortex is involved in spatial perception and spatial contrast which is abnormal in migraine patients during their attack [51] is also activated by thermal stimuli in the interictal period. Other studies do not shed light on functionality but confirm changes in this region–i.e., thinned cortex in the supramarginal temporal region reported in episodic migraine [52]. Some of these regions were strongly activated in the allodynic state as well during the ictal scan (compared to the interictal scan). The insula is well known to participate in migraine [34], and while the anterior component considered to be involved in autonomic changes and salience [53], the posterior insula is associated with pain related changes [54]. In the generalized allodynia group increased subcortical activations were also observed in regions (i.e., amygdala, caudate) associated with altered pain and emotional processing in migraine [37, 55–57]. Furthermore, the stronger activation of the brainstem in the allodynic state is in line with the notion that the brainstem is likely to contribute to the sensory processing alterations in migraine and might be related to the intensity of allodynia and pain hypersensitivity [15, 58].

Thus, multiple cortical and subcortical regions display significantly increased responsivity to noxious heat during the ictal (and allodynic) state. A few findings contribute to these enhanced profiles. As noted in our prior report, cortical connectivity may also be involved in cortico-cortical connections noted for sensory areas in migraine [35] suggesting integration of information across sensory modalities. The possible implication of this is that areas that are involved in sensory allodynia may contribute or be driven in part by other sensory manifestations of supersensitized primary cortical areas based on the concept of central sensitization. These increased activations suggest that normal function is in abeyance and thus multiple cortical behavioral manifestations may be present in patients with allodynia. These have not been evaluated in migraine patients.

4.6. Caveats

There are a number of caveats that are worth consideration: (1) Numbers of subjects: Getting scans of migraine patients during their migraine attack is a real challenge. Although 51 patients were recruited to the study, and our data for ictal vs. interictal comparison were reasonable (26 patients completed both visits, 19 patients were included in the final analysis), the patient numbers for comparisons for the allodynic vs. non-allodynic states were small. However, our psychophysics measures helped “phenotyping” these subgroups. Considering that this was a prospective clinical imaging study, and the first visit was the interictal scan, it might be understandable that less participants came back to the ictal visit where brain activations were measured during the migraine headache. Since we were interested in the neural responses to noxious stimulation in the ictal vs. interictal states and how the changes manifest in the generalized allodynic state are different from those in the non-allodynic state, we could only focus on patients who completed both visits. Also, no healthy controls were included in the study. (2) Migraine cycle: Although all interictal scans were performed at least 72 hours after the last migraine attack and the patients did not have migraine during the baseline visit, their headaches were not monitored following the interictal scan. This information might have been useful for determining the specific phase of the migraine cycle. (3) Allodynic state: In patients who experience mild forms of allodynia, it is likely that the symptoms of allodynia may not always be present or experienced similarly with every single attack. As such, symptoms experienced on ‘average’ by a patient (as it is also assessed by allodynia symptoms questionnaires) and the QST-based determinations at a given time may not always align. (4) Order effects: There is an order effect due to subjects not being randomized between interictal and ictal-phase scan sessions which was not possible to avoid in this study. (5) Laterality: When the patients came for their ictal scan the stimulation was applied on the side that patients were feeling the most pain, or if it was equal on both sides it was applied to the left side. For all of the subjects consistently the sites were the same for both scans. However, since for the analysis the data from the right-sided migraine patients were flipped and analyzed together with the data from the left-sided migraine patients, it is difficult to make any inferences about the laterality effects of the migraine itself in the results. Flipping the brain images could be also a limitation because of the structural and functional differences between the right and left hemispheres [59]. (6) Medication: The issue of medication (noted in S1 Table in S1 File) used by the patients is always a confound in differentiating between the disease and medication effects on the brain. The patients were asked to refrain from taking any medications for 48 hours before their scan. In this study, most patients were on either triptans or non-steroidal anti-inflammatory drugs (NSAIDs) or no medication. However, since this was a paired-wise comparison and each subject served as his/her own control, at an individual subject level the residual effects of medications were canceled out in the analysis of differences between interictal and ictal scans.

4.7. Conclusions

Allodynia during a migraine attack is quite common and present in over 60% of patients. Based on the present findings, there are distinct patterns and activations in the no vs. generalized allodynia patients and a potential trigeminovascular pathway was detected that may facilitate the development of generalized allodynia. In particular, it seems that increased functional activity of the spinal trigeminal nucleus and posterior thalamus during an attack distinguishes patients who are allodynic vs. those who are not during their attacks. Furthermore, the thalamic subregion that connects to the temporal and prefrontal cortices shows higher activation during a migraine attack. Previous studies have demonstrated activation of the temporal pole during an attack and changes in plasticity of the temporal cortex in patients with high migraine frequency [34, 36].

Supporting information

[Figure omitted. See PDF.]

S1 File.

https://doi.org/10.1371/journal.pone.0244320.s001

(DOCX)

Acknowledgments

The authors would like to thank Gabi Barmettler and Lauren Nutile for their assistance.

Citation: Maleki N, Szabo E, Becerra L, Moulton E, Scrivani SJ, Burstein R, et al. (2021) Ictal and interictal brain activation in episodic migraine: Neural basis for extent of allodynia. PLoS ONE 16(1): e0244320. https://doi.org/10.1371/journal.pone.0244320

References

1. Burstein R, Yarnitsky D, Goor-Aryeh I, Ransil BJ, Bajwa ZH. An association between migraine and cutaneous allodynia. Ann Neurol. 2000;47: 614–624. pmid:10805332

2. Bigal ME, Ashina S, Burstein R, Reed ML, Buse D, Serrano D, et al. Prevalence and characteristics of allodynia in headache sufferers. Neurology. 2008;70: 1525–1533. pmid:18427069

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

4. Burstein R, Cutrer MF, Yarnitsky D. The development of cutaneous allodynia during a migraine attack: Clinical evidence for the sequential recruitment of spinal and supraspinal nociceptive neurons in migraine. Brain J Neurol. 2000;123 (Pt 8): 1703–1709. pmid:10908199

5. Burstein R, Collins B, Jakubowski M. Defeating migraine pain with triptans: A race against the development of cutaneous allodynia. Ann Neurol. 2004;55: 19–26. pmid:14705108

6. Baldacci F, Vedovello M, Ulivi M, Vergallo A, Poletti M, Borelli P, et al. Triggers in allodynic and non-allodynic migraineurs. A clinic setting study. Headache. 2013;53: 152–160. pmid:23278553

7. Louter MA, Bosker JE, van Oosterhout WPJ, van Zwet EW, Zitman FG, Ferrari MD, et al. Cutaneous allodynia as a predictor of migraine chronification. Brain J Neurol. 2013;136: 3489–3496. pmid:24080152

8. Sprenger T, Borsook D. Migraine changes the brain: Neuroimaging makes its mark. Curr Opin Neurol. 2012;25: 252–262. pmid:22487570

9. Burstein R, Jakubowski M, Garcia-Nicas E, Kainz V, Bajwa Z, Hargreaves R, et al. Thalamic sensitization transforms localized pain into widespread allodynia. Ann Neurol. 2010;68: 81–91. pmid:20582997

10. DaSilva AFM, Becerra L, Makris N, Strassman AM, Gonzalez RG, Geatrakis N, et al. Somatotopic activation in the human trigeminal pain pathway. J Neurosci Off J Soc Neurosci. 2002;22: 8183–8192. pmid:12223572

11. Stankewitz A, Aderjan D, Eippert F, May A. Trigeminal nociceptive transmission in migraineurs predicts migraine attacks. J Neurosci Off J Soc Neurosci. 2011;31: 1937–1943. pmid:21307231

12. Young WB. Allodynia as a complication of migraine: Background and management. Curr Treat Options Neurol. 2009;11: 3–9. pmid:19094830

13. ICHD-II. The International Classification of Headache Disorders. Cephalalgia 24. 2004;24: 9–160. pmid:14979299

14. Beck AT. Beck Depression Inventory Manual. 2nd edition. San Antonio: TX: Psychological Corporation.; 1996.

15. Moulton EA, Burstein R, Tully S, Hargreaves R, Becerra L, Borsook D. Interictal dysfunction of a brainstem descending modulatory center in migraine patients. PLoS ONE. 2008;3. pmid:19030105

16. Russo A, Esposito F, Conte F, Fratello M, Caiazzo G, Marcuccio L, et al. Functional interictal changes of pain processing in migraine with ictal cutaneous allodynia. Cephalalgia. 2017;37: 305–314. pmid:27084886

17. Russo A, Tessitore A, Esposito F, Di Nardo F, Silvestro M, Trojsi F, et al. Functional changes of the perigenual part of the anterior cingulate cortex after external trigeminal neurostimulation in migraine patients. Front Neurol. 2017;8: 282. pmid:28663737

18. Tessitore A, Russo A, Esposito F, Giordano A, Taglialatela G, De Micco R, et al. Interictal cortical reorganization in episodic migraine without aura: an event-related fMRI study during parametric trigeminal nociceptive stimulation. Neurol Sci Off J Ital Neurol Soc Ital Soc Clin Neurophysiol. 2011;32 Suppl 1: S165–167. pmid:21533737

19. Woolrich MW, Jbabdi S, Patenaude B, Chappell M, Makni S, Behrens T, et al. Bayesian analysis of neuroimaging data in FSL. NeuroImage. 2009;45: S173–S186. pmid:19059349

20. Pendse G, Borsook D, Becerra L. Enhanced false discovery rate using Gaussian mixture models for thresholding fMRI statistical maps. NeuroImage. 2009;47: 231–261. pmid:19269334

21. Becerra L, Morris S, Bazes S, Gostic R, Sherman S, Gostic J, et al. Trigeminal neuropathic pain alters responses in CNS circuits to mechanical (brush) and thermal (cold and heat) stimuli. J Neurosci Off J Soc Neurosci. 2006;26: 10646–10657. pmid:17050704

22. Noseda R, Burstein R. Migraine pathophysiology: Anatomy of the trigeminovascular pathway and associated neurological symptoms, cortical spreading depression, sensitization, and modulation of pain. Pain. 2013;154 Suppl 1: S44–53. pmid:24347803

23. Moulton E, Pendse G, Morris S, Strassman A, Aiello-Lammens M, Becerra L, et al. Capsaicin-induced thermal hyperalgesia and sensitization in the human trigeminal nociceptive pathway: An fMRI study. NeuroImage. 2007;35: 1586–1600. pmid:17407825

24. Mazziotta J, Toga A, Evans A, Fox P, Lancaster J, Zilles K, et al. A probabilistic atlas and reference system for the human brain: International Consortium for Brain Mapping (ICBM). Philos Trans R Soc Lond Ser B. 2001;356: 1293–1322. pmid:11545704

25. Behrens TEJ, Johansen-Berg H, Woolrich MW, Smith SM, Wheeler-Kingshott C a. M, Boulby PA, et al. Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging. Nat Neurosci. 2003;6: 750–757. pmid:12808459

26. Johansen-Berg H, Behrens TEJ, Sillery E, Ciccarelli O, Thompson AJ, Smith SM, et al. Functional-anatomical validation and individual variation of diffusion tractography-based segmentation of the human thalamus. Cereb Cortex N Y N 1991. 2005;15: 31–39. pmid:15238447

27. Nahman-Averbuch H, Shefi T, Schneider VJ, Li D, Ding L, King CD, et al. Quantitative sensory testing in patients with migraine: A systematic review and meta-analysis. Pain. 2018;159: 1202–1223. pmid:29781957

28. Kunz M, Peter J, Huster S, Lautenbacher S. Pain and disgust: The facial signaling of two aversive bodily experiences. PLOS ONE. 2013;8: e83277. pmid:24349479

29. Butt JH, Rostrup E, Hansen AS, Lambertsen KL, Kruuse C. Induction of migraine-like headache, but not aura, by cilostazol in patients with migraine with aura. Brain J Neurol. 2018;141: 2943–2951. pmid:30137217

30. Christensen CE, Younis S, Deen M, Khan S, Ghanizada H, Ashina M. Migraine induction with calcitonin gene-related peptide in patients from erenumab trials. J Headache Pain. 2018;19: 105. pmid:30409109

31. Ghanizada H, Al-Karagholi MA-M, Arngrim N, Olesen J, Ashina M. PACAP27 induces migraine-like attacks in migraine patients. Cephalalgia. 2020;40: 57–67. pmid:31299857

32. Maniyar FH, Sprenger T, Monteith T, Schankin C, Goadsby PJ. Brain activations in the premonitory phase of nitroglycerin-triggered migraine attacks. Brain J Neurol. 2014;137: 232–241. pmid:24277718

33. Maleki N, Linnman C, Brawn J, Burstein R, Becerra L, Borsook D. Her versus his migraine: Multiple sex differences in brain function and structure. Brain. 2012;135: 2546–2559. pmid:22843414

34. Maleki N, Becerra L, Brawn J, Bigal M, Burstein R, Borsook D. Concurrent functional and structural cortical alterations in migraine. Cephalalgia. 2012;32: 607–620. pmid:22623760

35. Hodkinson DJ, Veggeberg R, Kucyi A, van Dijk KRA, Wilcox SL, Scrivani SJ, et al. Cortico-cortical connections of primary sensory areas and associated symptoms in migraine. eNeuro. 2017;3. pmid:28101529

36. Moulton EA, Becerra L, Maleki N, Pendse G, Tully S, Hargreaves R, et al. Painful heat reveals hyperexcitability of the temporal pole in interictal and ictal migraine states. Cereb Cortex N Y N 1991. 2011;21: 435–448. pmid:20562317

37. Wilcox SL, Veggeberg R, Lemme J, Hodkinson DJ, Scrivani S, Burstein R, et al. Increased functional activation of limbic brain regions during negative emotional processing in migraine. Front Hum Neurosci. 2016;10: 366. pmid:27507939

38. Schwedt TJ, Chiang C-C, Chong CD, Dodick DW. Functional MRI of migraine. Lancet Neurol. 2015;14: 81–91. pmid:25496899

39. Tedeschi G, Russo A, Tessitore A. Relevance of functional neuroimaging studies for understanding migraine mechanisms. Expert Rev Neurother. 2013;13: 275–285. pmid:23448217

40. Mehnert J, May A. Functional and structural alterations in the migraine cerebellum. J Cereb Blood Flow Metab. 2019;39: 730–739. pmid:28737061

41. Edwards RR. Individual differences in endogenous pain modulation as a risk factor for chronic pain. Neurology. 2005;65: 437–443. pmid:16087910

42. Staud R, Weyl EE, Price DD, Robinson ME. Mechanical and heat hyperalgesia highly predict clinical pain intensity in patients with chronic musculoskeletal pain syndromes. J Pain Off J Am Pain Soc. 2012;13: 725–735. pmid:22739051

43. Staud R. Abnormal pain modulation in patients with spatially distributed chronic pain: fibromyalgia. Rheum Dis Clin North Am. 2009;35: 263–274. pmid:19647141

44. Vierck CJ. Mechanisms underlying development of spatially distributed chronic pain (fibromyalgia). Pain. 2006;124: 242–263. pmid:16842915

45. Vogt BA. Pain and emotion interactions in subregions of the cingulate gyrus. Nat Rev Neurosci. 2005;6: 533–544. pmid:15995724

46. Jahn A, Nee DE, Alexander WH, Brown JW. Distinct regions of anterior cingulate cortex signal prediction and outcome evaluation. NeuroImage. 2014;95: 80–89. pmid:24667454

47. Vogt BA. Midcingulate cortex: Structure, connections, homologies, functions and diseases. J Chem Neuroanat. 2016;74: 28–46. pmid:26993424

48. Leech R, Sharp DJ. The role of the posterior cingulate cortex in cognition and disease. Brain J Neurol. 2014;137: 12–32. pmid:23869106

49. Rudebeck PH, Rich EL. Orbitofrontal cortex. Curr Biol CB. 2018;28: R1083–R1088. pmid:30253144

50. Rolls ET. Convergence of sensory systems in the orbitofrontal cortex in primates and brain design for emotion. Anat Rec A Discov Mol Cell Evol Biol. 2004;281A: 1212–1225. pmid:15470678

51. Benedek K, Tajti J, Janáky M, Vécsei L, Benedek G. Spatial contrast sensitivity of migraine patients without aura. Cephalalgia. 2002;22: 142–145. pmid:11972583

52. Magon S, May A, Stankewitz A, Goadsby PJ, Schankin C, Ashina M, et al. Cortical abnormalities in episodic migraine: A multi-center 3T MRI study. Cephalalgia. 2019;39: 665–673. pmid:30525946

53. Borsook D, Veggeberg R, Erpelding N, Borra R, Linnman C, Burstein R, et al. The Insula: A “Hub of Activity” in Migraine. Neurosci Rev J Bringing Neurobiol Neurol Psychiatry. 2016;22: 632–652. pmid:26290446

54. Lu C, Yang T, Zhao H, Zhang M, Meng F, Fu H, et al. Insular cortex is critical for the perception, modulation, and chronification of pain. Neurosci Bull. 2016;32: 191–201. pmid:26898298

55. Szabó E, Galambos A, Kocsel N, Édes AE, Pap D, Zsombók T, et al. Association between migraine frequency and neural response to emotional faces: An fMRI study. NeuroImage Clin. 2019;22. pmid:31146320

56. Maleki N, Becerra L, Nutile L, Pendse G, Brawn J, Bigal M, et al. Migraine attacks the Basal Ganglia. Mol Pain. 2011;7: 71. pmid:21936901

57. Borsook D, Upadhyay J, Chudler EH, Becerra L. A key role of the basal ganglia in pain and analgesia—insights gained through human functional imaging. Mol Pain. 2010;6: 27. pmid:20465845

58. Chong CD, Plasencia JD, Frakes DH, Schwedt TJ. Structural alterations of the brainstem in migraine. NeuroImage Clin. 2016;13: 223–227. pmid:28003961

59. Hervé P-Y, Zago L, Petit L, Mazoyer B, Tzourio-Mazoyer N. Revisiting human hemispheric specialization with neuroimaging. Trends Cogn Sci. 2013;17: 69–80. pmid:23317751