1. Introduction {#s0005}
===============

Hypersensitivity to sensory input is one of the core symptoms of the migraine attack ([@bb0255]; [@bb0245]) but in principle exists already interictally in migraineurs ([@bb0270]). This irritability to external input becomes pronounced in the preictal phase of the migraine cycle ([@bb0100]; [@bb0150]; [@bb0180]) which may explain the perception of prodromal symptoms. The migraine brain processes sensory input such as auditory ([@bb0010]; [@bb0300]), olfactory ([@bb0080]; [@bb0275]), somatosensory ([@bb0040]), visual ([@bb0030]; [@bb0095]; [@bb0175]) or nociceptive ([@bb0205]) stimuli different than a non-migraine population, which has been discussed to be due to hyperresponsivity of primary sensory areas ([@bb0030]) or a lack of habituation in response to repetitive visual stimulation ([@bb0020]; [@bb0055]; [@bb0070]; [@bb0160]).

Recently, it has been suggested that not only the individual primary sensory systems are altered but also multisensory processes ([@bb0030]; [@bb0215]) raising the question whether multisensory integration ([@bb0185]; [@bb0250]) and appraisal itself are also affected in migraine patients ([@bb0035]; [@bb0270]).

To test our hypotheses of enhanced functional coupling in migraineurs during nociception and visual stimulation, we adapted a robust trigeminal pain model used in fMRI studies ([@bb0145]; [@bb0245]; [@bb0280]) for electroencephalographic studies and compared interictal episodic migraineurs with healthy controls using concurrent source localization methods and measures of functional coupling.

2. Material and methods {#s0010}
=======================

2.1. Subjects and experimental design {#s0015}
-------------------------------------

Thirty-one female episodic interictal migraineurs (minimum of 24 h before and after the pervious/next attack) and 31 pairwise age matched female healthy volunteers without a history of headache participated. Male patients were not excluded, but since migraine is more common in females, by chance no male patient volunteered to take part in the study. The definition of migraine followed IHS classification ([@bb0115]) and all migraineurs filled out a migraine diary. Twelve migraineurs suffered from migraine with visual aura (fortifications). All patients reported increased sensitivity to light and nausea accompanying their attacks. 27 patients further report increased sensitivity to sound. All migraineurs were called after the day of the experiment until the next attack occurred to define the number of days until next attack. The established paradigm ([@bb0195]; [@bb0260]; [@bb0280]) on nociception and visual excitation consists of separately applied 15 nociceptive, trigeminal stimuli (gaseous ammonia), 15 olfactory stimuli (gaseous rose odor) and 15 control stimuli (air puffs) to stimulate the trigeminal system. These three conditions were implemented as single boluses into a constant stream of air which flowed through a Teflon-tube into the participants left nostril. Additionally, all participants received 15 trials of binocular checkerboard stimulation. Stimuli were presented (jittered from 8 to 10 s) in a randomized order with the only constrain that two ammonia pulses where not allowed to follow each other directly to prevent habituation. After each stimulation, participants were asked to rate the intensity (in case of rose odor, air puffs and visual stimulation) or level of pain (in case of ammonia) on a scale from 0 to 100 (where 0 means no pain at all and 100 worst imaginable pain) and also unpleasantness on a scale from −50 to 50 (where −50 means very pleasant and 50 very unpleasant). The exact timing is shown in [Fig. 1](#f0005){ref-type="fig"}.Fig. 1Timeline of the experimental design. The red cross marks the onset of a reaction task to monitor the subjects attention. Actual stimuli are delivered either through a Teflon-tube with a constant stream of air into the subjects left nostril (for gaseous ammonia, rose odor, air puffs) or via a monitor (visual stimulation). After each trial the intensity and unpleasantness of the stimulus was rated. Each condition was repeated 15 times and trials were presented in randomized order.Fig. 1

The repetitive visual stimulation used a circular flickering checkerboard pattern consisting of 5 rings each with 12 white or black squares. The patterns color (black/white) was reversed at a rate of 8 Hz to elicit steady state visually evoked potential (SSVEP) ([@bb0210]). White was defined as highest and black as lowest possible luminance of the monitor. The stimuli were presented on a TFT Monitor (Samsung Syncmaster P2370(G), refreshing rate: 75 Hz, resolution: 1920 × 1080 pixel, length of diagonal: 23 in., width: 51 cm, height: 28.8 cm, reaction time of the monitor: 2 ms, max. luminance: 250 cd/m^2^) controlled by a NVIDIA (Geforce GTX 1080) video card. The visual stimulus had a diameter of 14.8 cm and was in 75 cm distance from the subjects head. The horizontal and vertical visual angles are accordingly 10.9°. The background was throughout the experiment set to 50% gray (125 cd/m^2^). Patients were seated vertically and horizontally centered relative to the screen\'s center.

During the experiment we measured electroencephalography (EEG) at 60 locations of the 10--20 system ([@bb0045]) at the head using a BrainAmp device (Brain Products, Germany) with active electrodes sampled at 1000 Hz. Four additional electrodes were used to obtain vertical and horizontal electrooculogram (EOG) to attenuate for eye movement and blinking artifacts. Two subjects (1 migraine, 1 control) were removed from the analysis due to artifacts. The final group therefore consisted of 30 migraineurs and 30 healthy controls for which we report the demographics in [Table 1](#t0005){ref-type="table"}. Ethical approval was obtained from the local ethics authority (PV4896). Patients gave written informed consent.Table 1Demographic and clinical characteristics.Table 1PatientsHealthy ControlsN3030Age (mean ± SD \[min, max\])27,2 ± 6,6 (20, 45)27.8 ± 8.0 (19,52)Without Aura/with Aura18/12--Days with headache per month (mean ± SD \[min, max\])4,3 ± 2,7 (1, 12)--Attacks per Month2,8 ± 2,1 (1, 8)--Duration of disease in years (≤11/12--17/≥ 18)9/10/7--Days until next attack16,29 ± 17,01 (2, 71)--Days after last attack15,32 ± 11,33 (3, 55)--

2.2. Preprocessing of EEG data {#s0020}
------------------------------

EEG raw data were re-referenced to a common average, cut into epochs between −500 and 3000 ms regarding stimulus onset, and high-pass filtered at 0.5 Hz using the FieldTrip toolbox ([@bb0230]). Powerline artifacts were reduced by a notch filter at 50 Hz. Eye movement and blinking artifacts were automatically eliminated by regressing out the electrode differences from the two vertical as well as from the two horizontal EOG channels using the procedure described by Parra and colleagues ([@bb0235]). Supplementary Fig. S1 demonstrates the (missing) correlation between EOG and EEG channels after artifact removal. Thereafter, all trials passed the automated muscle detection routine of FieldTrip and an overall z-score higher than 15. All identified artifact loaded trials were completely rejected from further analysis, leaving 94.9% of the trials (87.8% in ammonia condition, 95.9% in rose, 98.7% in checker and 97.3% in air condition) for the analysis. There were no significant differences in the number of remaining trials between groups over all trials as well as in the individual conditions (two-sample, two-tailed *t*-tests, all *p* \> .2).

Time frequency transformation of the individual trials were calculated using the multi-taper method ([@bb0155]; [@bb0200]; [@bb0290]) for frequencies window of 2 to 100 Hz with frequency step of 1 Hz and frequency resolution from 1 to 10 Hz, depending on the frequency under observation (higher resolution for higher frequencies) using the implementation of the SPM12 toolbox (<http://www.fil.ion.ucl.ac.uk/spm/>). Temporal resolution was set to 800 ms and the temporal steps to 50 ms. The resulting time-frequency spectra were - on a single trial level - recalculated as relative changes to baseline (defined as 500 to 0 ms before stimulus onset) by division, logarithmically transformed, and then averaged within the individuals showing the induced responses ([@bb0065]) following the robust averaging protocol within SPM12. Afterwards the individual averages were cropped to a temporal window from 0 to 2500 ms regarding stimulus onset.

2.3. Analysis of time frequency dynamics {#s0025}
----------------------------------------

From the conditions "rose odor" and "ammonia" were subtracted the control condition "air puffs". We then calculated a one sampled *t*-test for the main effect of nociception\>control, rose\>control and visual stimulation across all subjects and a further two sample t-test to show differences between healthy controls and migraineurs on smoothed (1 Hz, 500 ms Gaussian Kernel) image projections in the frequency range of 1 to 25 Hz at the electrodes Pz and Oz. These electrode positions were chosen from the literature regarding time-frequency responses after pain and trigeminal input ([@bb0110]), or visual input, respectively. Results are reported positive if they passed the FWE-corrected threshold of *p* \< .05. Frequency bands are defined as: beta 12.5--28 Hz, alpha 8--12 Hz, theta 3.5--7.5 Hz, and delta 0.5--3 Hz.

2.4. Source localization and phase coupling {#s0030}
-------------------------------------------

We run source localization for the main effect of repetitive visual input with exact low resolution brain electromagnetic tomography (eLORETA) implemented into the MATLAB toolbox METH by Guido Nolte (<https://www.uke.de/dateien/institute/neurophysiologie-und-pathophysiologie/downloads/meg_methods.zip>) and using a lead-field with 5003 voxel covering the whole brain (cortex, midbrain and cerebellum) to locate power values from the cross-frequency-spectrum. We further used a linearly constrained minimum variance (LCMV) beamformer to show the consistency between both methods as suggested by [@bb0170]. All voxel showing power change (dB) relative to baseline had to further pass FDR correction (p \< .05) to deal with multiple comparison e.g. number of voxel. The peak voxel of the main effect of repetitive visual stimulation was used as a seed for the analysis on coupling differences by Multivariate Interaction Measure (MIM) in source space between migraineurs and healthy controls ([@bb0085]). This method is a multivariate interaction measure, in fact a multivariate version of imaginary coherence, and is calculated here on source level to include all three source orientations for each voxel and does not require the specification of a source orientation. It is robust to artifacts of volume conduction and invariant to linear and static transformations of data within each voxel. Single subject results for each voxel entered a multivariate permutation test to assess group differences, since MIM-values are not normally distributed, and *p*-values were corrected with the *tmax*-method for multiple comparisons ([@bb0025]; [@bb0105]; [@bb0310]). This procedure was implemented in MATLAB by David Groppe (Version 1.1, <https://de.mathworks.com/matlabcentral/fileexchange/54585-mult-comp-perm-t2>). For visual inspection and to further test relations with clinical markers we also included results at an exploratory statistical threshold of *p* \< .01 uncorrected.

2.5. Co-modulation with clinical markers of migraine {#s0035}
----------------------------------------------------

Behavioral data as well as results from the EEG (power in individual frequency bands, and coupling strength) were additionally tested for differences between patients with and without aura. Since we noticed quite long interictal periods in our patients with a broad variance (i.e. standard deviation), we also correlated the power of the individual frequency bands with the length of the interictal period as well as with the severity of the disease ([@bb0055]). Results from the coupling analysis, i.e. the Multivariate Interaction Measure, were furthermore (non-parametrically) correlated with the patients\' clinical parameters, (i) headache days per month, (ii) number of headache free days before, and (iii) after the EEG measurement by means of Spearman\'s rho.

2.6. Data availability {#s0040}
----------------------

Unconditional access to anonymized data is available to qualified investigators on request to the corresponding author.

3. Results {#s0045}
==========

3.1. Behavior {#s0050}
-------------

While two-sample, two-tailed *t*-tests did not reveal any significant differences in the intensity ratings between healthy controls and migraineurs, the unpleasantness ratings of the visual and rose odor were significantly (*p* = .002 and *p* = .04) higher in migraineurs compared to healthy controls. This effect was not seen in unpleasantness ratings of ammonia or air puff ratings (compare e.g. [@bb0190]). The behavior was not significantly different between patients with and without aura.

3.2. Time-frequency dynamics {#s0055}
----------------------------

The main effect of nociceptive stimulation showed significant increases in theta (for both electrodes of interest) and beta band (electrode Oz). Decreases could be shown in the alpha band after roughly 1.5 s regarding stimulus onset at electrode Pz. The comparison of migraineurs and healthy controls did not reveal significant differences in the nociceptive domain at the conservative threshold (FWE, *p* \< .05) chosen. The main effect of visual stimulation was significant in three time-frequency clusters and showed a robust effect of SSVEP at 8 Hz in Oz (see [Table 2](#t0010){ref-type="table"} and [Fig. 2](#f0010){ref-type="fig"}). Contrasting migraineurs and healthy controls revealed increased power of the SSVEP (8 Hz) in a temporal window of 350 to 1150 ms regarding stimulus onset in patients only with an effect size of 1.46 (Cohen\'s d).Table 2Results for time-frequency analysis of the EEG, FWE corrected (*p* \< .05) effects in the time-frequency domain.Table 2Electrode positionPower increase (+) or decrease (−)Frequency bandFrequency \[Hz\]Time \[ms\]T-valuep-Value (FWE-corrected)Degree of freedomNociceptionPz+Delta/theta3--6350--11505.580.00059−Alpha9--101450--20003.750.01459Oz+Delta/theta2--5100--18007.970.00059+Beta13--23750--20004.400.00359  Main effect visualPz+Delta/theta/alpha2--8100--60010.020.00059−Delta/theta2--7700--200010.370.00059−Beta13--21300--20008.950.00059Oz+Delta/theta/alpha2--9100--50012.480.00059+Flickering (SSVEP)8100--20005.560.00059−Alpha13200--20006.200.00059−Delta/theta2--6700--20008.130.00059  Group effect visual (controls \< patients)Oz+Flickering (SSVEP)8350--11505.580.03858Fig. 2Effect of visual repetitive (top) and trigeminal nociceptive stimulation (bottom) in the time-frequency domain from 0 to 2.5 s and 0 to 30 Hz for the electrodes Pz (top) and Oz (bottom). Significant effects (FWE-corrected, p \< .05) are signified by black and - for the SSVEP - white rectangles. Displayed is the logarithmic transformation of power P relative to baseline power P~base~.Fig. 2

3.3. Source localization and coupling {#s0060}
-------------------------------------

The main effect of the SSVEP showed highest power in early visual areas (Cuneus, MNI coordinates 0, −83, 8 \[xyz/mm\]) (see [Fig. 3](#f0015){ref-type="fig"}, top) which was confirmed using the eLORETA and the LCVM beamformer as source localization tools (see Supplementary Fig. S2). The analysis of differences in coupling of the SSVEP showed that migraineurs had a significant (corrected for multiple comparison) increased coupling from early visual to two peak voxel located in the left temporal pole, more specifically at the coordinates \[−68, −8 −23\] and \[−60, −8, −23\] (xyz, mm in MNI-space, *p* = .0404 and *p* = .0344 adjusted for multiple comparisons, respectively). Visual inspection at lower statistical threshold (uncorrected, *p* \< .01) further reveals increased coupling to two clusters namely left anterior cingulate cortex (ACC, peak voxel coordinate \[−8, 38, 30\]) and right premotor cortex (PMC) (peak voxel coordinate \[30, 0, 68\]) as shown in [Fig. 3](#f0015){ref-type="fig"} (bottom).Fig. 3Source localization of the FDR-corrected (*p* \< .05) main effect (in logarithmic transformation of power P relative to baseline power P~base~) of repetitive visual stimulation (top) and increased coupling in the migraineurs (bottom) to premotor cortex (PMC), anterior cingulate cortex (ACC) and temporal pole (TP) by means of multivariate interaction measure (MIM) taking the peak voxel of the main effect as seed (marked in blue) for the phase (noted as negative decadic logarithmic transformation of uncorrected p-values lower than 0.01). Sagittal slices are presented for 4 x-coordinates in MNI space.Fig. 3

3.4. Co-modulation with clinical markers of migraine {#s0065}
----------------------------------------------------

The MIM between Cuneus and the two aforementioned significant clusters in the middle temporal gyrus (MTG) (*r* = 0.340, *p* = .035 and *r* = 0.370, *p* = .024, respectively) and the ACC (*r* = 0.467, *p* = .005) correlated positively with the number of headache free days before the experiment (Spearman\'s Rho, one sided test) see [Fig. 4](#f0020){ref-type="fig"}. This means that the coupling was stronger the longer the duration of the last attack. No correlation (positive or negative) was found regarding the coupling of the cuneus and any cortical areas with the days until the next attack. We found no significant differences between patients with and without aura. Furthermore the severity and length of the interictal period was not correlated to the power of the individual frequency bands in neither of the two conditions.Fig. 4Significant correlation (Spearman) of headache free days since the last migraine attack and coupling strength from cuneus to both locations in the temporal pole (TP 1 and TP 2) and anterior cingulate cortex (ACC).Fig. 4

4. Discussion {#s0070}
=============

Although our nociceptive stimulation paradigm showed robust results in fMRI ([@bb0265]; [@bb0245]; [@bb0275]) and reproduced pain related electrophysiological effects such as increases in theta- and decreases in alpha band ([@bb0140]; [@bb0240]; [@bb0285]) we did not find significant differences in nociceptive processing between migraineurs and controls despite investigating rather large groups. Such differences are therefore mainly either sub-cortical (brainstem ([@bb0305]), thalamus ([@bb0220]), hypothalamus ([@bb0245])) and therefore not accessible by EEG, or rather subtle alterations in the cortical transmission of nociceptive input in migraineurs. Another possibility is that such differences exist but perhaps more in the ictal phase. Although there are some reports about visual stimulation and EEG changes in migraine patients ([@bb0020]; [@bb0055]; [@bb0070]), we note a surprising lack of such data regarding trigeminal pain and nociceptive input ([@bb0075]).

In contrast to the nociceptive system, the neuronal response following visual input showed robust differences between groups: Interictally, episodic migraineurs showed a significant higher activation of the primary visual area (cuneus) in response to repetitive visual stimulation. This increased response was notably in the same frequency band as the visual input flicker, corresponding to steady-state visually evoked potentials (SSVEP). This finding can be seen in the light of a lack of habituation to slowly flickering light in migraineurs. Of note, fast flickering frequencies around 20 Hz, known as photic driving, show inconsistent results ([@bb0050]; [@bb0090]). Nevertheless, similar findings of a lack of habituation are known to happen in migraine patients also in other sensory domains ([@bb0165], [@bb0160]) and have been discussed as a hyperresponsive cortical reaction to sensory (here: visual) input. Hyperresponsivity to visual stimulation in interictal migraine patients has also been reported in a previous PET study ([@bb0030]) for two different light intensities, and also in fMRI ([@bb0135]). More importantly, using concurrent analyzing methods of inter-regional coupling invariant to volume conduction we could link this finding to other cortical hubs, especially the temporal pole ([@bb0060]; [@bb0205]), the premotor cortex and the ACC. Although altered functional connectivities have been reported in healthy subjects following stressful visual stimulation ([@bb0130]), our finding of enhanced functional connectivity of the aforementioned cortical areas is remarkable as it shows harmonized oscillations of widespread cortical areas as a response to visual input. This finding prompts the question whether the hypersensitive appearance of migraineurs, known to have lower detection thresholds in numerous sensory dimensions (pain, heat, vision, auditory, touch ([@bb0120])) and altered thresholds in contrast detection and orientation ([@bb0225]; [@bb0295]), may indeed prove to be due to a generally altered stimulus-response pattern of the migraine brain. This migraine specific stimulus-response pattern would meet and possibly drive oscillations in specific single brain areas such as the thalamus or hypothalamus which are currently thought to define thresholds and generate impulses for attack generation ([@bb0005]; [@bb0015]; [@bb0190]; [@bb0245]; [@bb0275]).

Our data show not only generally elevated levels of coupling between the visual system, the motor system, the temporal pole and ACC but link them to the cycling of the disease, as the activity correlation within these systems and the days since the last attack increases.

Perception changes of sensory input during a migraine cycle are a continuum across different phases with a different magnitude of sensitivity ([@bb0150]; [@bb0180]; [@bb0255]), but also a different connectivity strength of cortical coupling. This is further evidence for a migraine specific disturbance in cortical processing and underlines the importance to focus on migraine cyclic behavior when investigating this disease. We note that we could not show a correlation towards the next attack which may be due to the higher variance of days towards the next attack in our cohort. It may be necessary to investigate patients with a higher attack frequency to prove such correlations. Nevertheless, our finding suggest migraine specific alterations of cortical pathways which is in line with recent results gained from resting-state fMRI ([@bb0125]).

Appendix A. Supplementary data {#s0075}
==============================

Supplementary materialImage 1

This work was supported by the German Research Foundation, SFB936/A5 to A.M., and SFB936/Z3 to G.N.

Supplementary data to this article can be found online at <https://doi.org/10.1016/j.nicl.2019.101815>.
