- Research article
- Open Access
Structural brain network characteristics in patients with episodic and chronic migraine
The Journal of Headache and Pain volume 22, Article number: 8 (2021)
Migraine is a primary headache disorder that can be classified into an episodic (EM) and a chronic form (CM). Network analysis within the graph-theoretical framework based on connectivity patterns provides an approach to observe large-scale structural integrity. We test the hypothesis that migraineurs are characterized by a segregated network.
19 healthy controls (HC), 17 EM patients and 12 CM patients were included. Cortical thickness and subcortical volumes were computed, and topology was analyzed using a graph theory analytical framework and network-based statistics. We further used support vector machines regression (SVR) to identify whether these network measures were able to predict clinical parameters.
Network based statistics revealed significantly lower interregional connectivity strength between anatomical compartments including the fronto-temporal, parietal and visual areas in EM and CM when compared to HC. Higher assortativity was seen in both patients’ group, with higher modularity for CM and higher transitivity for EM compared to HC. For subcortical networks, higher assortativity and transitivity were observed for both patients’ group with higher modularity for CM. SVR revealed that network measures could robustly predict clinical parameters for migraineurs.
We found global network disruption for EM and CM indicated by highly segregated network in migraine patients compared to HC. Higher modularity but lower clustering coefficient in CM is suggestive of more segregation in this group compared to EM. The presence of a segregated network could be a sign of maladaptive reorganization of headache related brain circuits, leading to migraine attacks or secondary alterations to pain.
Migraine is a multifactorial neurovascular disorder which affects about 12% of the general population  and rates among the most disabling diseases [2, 3]. In the episodic form of migraine (EM), headache occurs less than 15 days per month, whereas in the chronic form (CM), it occurs on 15 or more per month for at least three consecutive months . Apart from functional brain alterations [5,6,7,8], several [9,10,11,12,13,14,15] but not all studies  reported multi-regional volumetric alterations in white mater and gray matter. Grey matter volume (GMV) alterations can get worse (i.e. seen as decreases in GMV) over time, e.g. in sensory-discriminative brain regions , although the existences of longitudinal changes in GMV are still under debate [17, 18]. Regarding the spatial locations of migraine-related GMV changes, previous studies reported alterations in the occipital , frontal [8, 11, 14, 19,20,21], temporal cortex , somatosensory , parietal cortex [19, 20] and cerebellar regions . Similarly, other studies demonstrated altered cortical thickness (CT) in migraineurs compared to controls, seen in the visual cortex [23,24,25], somatosensory cortex , frontal cortex , and temporo-parietal cortex .
However, a different analysis framework is required - moving away from regional GMV or CT differences – to assess the inter-relation between the described brain regions on a network level. This approach would allow relating GMV connectedness and integrity to symptom severity (i.e., EM and CM) and would thus allow a systematic and integrative way to analyze structural abnormalities in patients with migraine. In particular, a network analysis using a graph theoretical framework has been extensively used for observing effect of various disorders in brain network integrity. This framework considers brain regions as nodes, and the interrelations between them as edges to form a network [28, 29]. The network formation begins with the collection of relational data among elements of a neurobiological system, which may vary from anatomical networks of associations between morphometry of cortical regions, inter-regional white matter projections, or multi-dimensional time series and their statistical dependencies or causal relations in behavior in social interactions. Once this data is corrected, normalized and assembled into the mathematical form of a graph or network, the common mathematical framework of graph theory is applied to obtain a set of measures to observe different alterations in the network . Using the anatomical networks reconstructed using the gray matter volumes (GMV) and functional network using resting-state functional magnetic resonance imaging (rs-fMRI)-correlations, Liu and colleagues demonstrated that network characteristics were disrupted in females with EM . From the networks from diffusion tensor imaging (DTI) and rs-fMRI it has been shown that CM exhibited altered rich club organization (higher connection density, abnormal small-world organization with increased global efficiency) compared to healthy controls (HC). It was further concluded that the higher ‘bridgeness’ in patients with non-rich club regions might increase the integration among pain-related brain circuits with more excitability but less inhibition for the modulation of migraine . In summary, the graph theory findings indicate that migraineurs lose structural network integrity, most likely seen as maladaptive integration among pain-related brain circuits, resulting in a disturbed balance of neuronal excitation and inhibition.
Yet, a systematic comparison of structural morphometric measures complimented with the brain network analysis has not been performed between patients with low (EM) and high (CM) occurrence of monthly migraine attacks. We thus examined cortical and subcortical morphometric changes leading to brain network reorganization in both patients with EM and CM relative to HC. Based on extant findings, we hypothesize to see the strongest global structural network alterations, seen as dis-integrated networks in patients with CM compared to HC and EM.
Design and study duration
This is the primary analysis of the reported data using a cross-sectional design. Other imaging data (MR spectroscopy (MRS) and Arterial Spin Labeling (ASL) MR imaging) have been collected for all participants and results are presented elsewhere . No statistical power calculation was conducted prior to the study. The sample size was based on the available data (during the study interval) and was similar to a recent ASL study in episodic migraine patients . All data was collected between December 2013 and July 2015.
19 right-handed HC, 17 right-handed patients with EM and 12 patients with CM were included for the study. The detailed demographic data are listed in Table 1. All patients fulfilled the modified ICHD-III-beta diagnostic criteria for EM or CM . None of the HC demonstrated signs of EM or CM according to these criteria (family history of migraine was allowed). During the enrollment process, we excluded all patients, which suffered from comorbid tension-type headache. Six of 12 CM also fulfilled the criteria for medication overuse headache (MOH), which is line to the literature [36, 37]. MOH is defined as headache that develops or significantly worsens during overuse of acute pain medication . For all participants, exclusion criteria were severe psychiatric disorders, cardiac problems (e.g. severe hypertension), other headache disorders or other neurologic disorders such as epilepsy, stroke, traumatic brain injury, neck injury or cerebrovascular disease. All participants completed prospective headache diaries, the Migraine Disability Assessment (MIDAS)  and Hamilton Anxiety (HADS-A) and Depression (HADS-D) Score  questionnaires. Acute and prophylactic medication was recorded prior to the study interval for each patient (see Table 2). We assessed the attacks/month based on the MIDAS questionnaire. Here, the label “headache attacks/month rate” (Table 1) reflects the average number of migraine headache days in the last three months prior to the MRI (i.e., an attack frequency of 4.0 in EM means that the average number of headache days was four per month across in this group). We recorded aura occurrence in all patients electronically in a table. Patients were free from migraine attacks at least 48 h before and after the scan. The study was approved by the ethics committee of canton Zurich (KEK number E-37/2007), Switzerland. All subjects provided written informed consent prior to study enrolment. Both groups received 50 Swiss Francs reimbursement for their study participation. Patients were recruited by advertisement (Intranet of the Hospital and mailing lists) and word-of-mouth.
Whole-brain magnetic resonance imaging (MRI) was performed on a 3 T scanner (Philips Ingenia, Netherlands) with a 32-channel receive-only head coil at the Neuroimaging Center of the University Hospital Zurich. 3D T1-weighted magnetization prepared rapid gradient echo (MPRAGE) sequence was acquired for each subject with TE/TI/TR = 2.52/900/1900 ms, flip angle = x°, field of view (FOV) = 256 × 256 mm2, matrix size = 256 × 256, slab thickness = 192 mm, voxel size = 1 × 1 × 1 mm3. Subjects’ scans were examined for any major anatomical abnormalities by an experienced neuroradiologist.
Cortical and subcortical morphometric analysis
Data from all subjects were analyzed using FreeSurfer version 5.3.0 (http://surfer.nmr.mgh.harvard.edu). This automated anatomic parcellation procedure enables one to extract reliable estimates of various cortical and subcortical measures including thickness, volume, area, curvature etc. . The procedure includes several steps: intensity normalization, skull stripping, Talairach transformation, and atlas-based assignment of neuro-anatomical labels, which are described in detail in previous studies [41, 42]. To describe here briefly, all subjects were run using the “recon-all” processing stream with default parameters to create a cortical surface model. This process includes above-mentioned procedures along with motion correction, averaging of multiple T1 volumes, removal of non-brain tissue and grey matter white matter boundary tessellation to create the surface model. This obtained model is then further used with its intensity and continuity information from the entire 3D volume in segmentation and deformation procedures to generate cortical thickness, calculated as the closest distance from the gray/white boundary to the gray/CSF boundary at each vertex on the tessellated surface . This process of obtained morphometric measures have been validated using histological  and manual measurements  and have been demonstrated to have very high reliability across different scanner and field strengths . These cortical and subcortical morphometric measures are very effective in depicting the regional alterations. However, the effect of these regional changes to other associated regions could have a significant impact in overall information transfer leading to various functional modifications. Hence, to observe these network level differences between the groups the computed cortical thickness (CT) and subcortical volumes (SCV) from all the subjects were further processed using brain network analysis.
Brain network analysis
Graph theoretical measures of network modularity, distance, and local information transfer was computed using the CT and subcortical volumes obtained from FreeSurfer using Brain Connectivity Toolbox  (https://sites.google.com/site/bctnet/). The group level correlations between the cortical regions and subcortical volumes and the differences between them were then computed in different network densities for observing the steady topological changes [47, 48]. The details of the analysis have been explained elsewhere [49, 50].
Among different measures computed in the study, below is the brief overview of those relevant for the study, with simplistic illustration in Fig. 1.
Modularity is a measure of the degree, to which the network is subdivided into densely interconnected nodes (modules) with sparse connections to other network or modules. Louvain algorithm was used for computing the modularity which is a hierarchical clustering algorithm, that recursively merges communities into a single node and executes the modularity .
Transitivity is the ratio between the number of triangles and the number of triplets in the graph.
Assortativity is correlation coefficient between the degrees of all nodes on two opposite ends of a link, higher (positive) assortativity indicating the nodes tend to link to other nodes with the same or similar degree.
Clustering coefficient is the fraction of triangles around a node representing the node’s neighbors that are also neighbors of each other.
Betweenness centrality of a node is the fraction of all shortest paths in the network that contain a given node. A node with higher edge-betweenness centrality participates in a large number of shortest paths.
In addition, network-based statistic (NBS) was used to assess differences in the inter-regional connectivity between the groups. NBS analysis performs the mass-univariate testing at every connection comprising the graph controlling for multiple comparisons through evaluating the null hypothesis at the level of interconnected subnetworks rather than individual connections . Here, the connectivity matrices obtained from the association of CT and SCV between the regions across a range of network densities were subjected to NBS analysis. The analysis primary goal was to identify the sub-network with the regions shown to have significant difference in various network properties using graph theoretical measures. Further details regarding the procedure are mentioned elsewhere [53, 54].
For the graph theoretical framework analysis, we used the CT values of a network comprising 34 cortical regions in each hemisphere based on Desikan atlas . Subsequently, we used the sub-cortical volumes from nine regions in each hemisphere and brain stem for the sub-cortical network. For assessing the statistical significance of graph metrics between patients and HC, a nonparametric permutation tests with 5000 iterations were applied [56, 57]. Given that, CT is sensitive to age and sex, they were further used as covariates for the analysis. In each repetition, the regional data for each subject were randomly reassigned to one of the two groups and an association matrix was obtained. The network measures were then calculated for all the networks at each density. Here, density represents cost of the network computed by fraction of present connections to all possible connections. Hence, the network measures derived at each density would specify the alterations in network behavior at different levels of fragmentation (from full, partial to discontinuous connectivity). This method of thresholding ensures that all the regions (nodes) of the network are connected while discarding spurious connections (edges) [47, 58]. The actual between-group difference in network measures was then placed in the corresponding permutation distribution and a two-tailed p-value (at 5% significance level, false discovery rate (FDR) corrected) was calculated based on its percentile position .
To assess the statistical significance for CT correlations with different clinical parameters, QDEC – a FreeSurfer statistical toolbox was used. Here, surface maps depicting regions with significant differences in the correlation with CT at each vertex were determined with general linear models (GLMs) using p < 0.001 as the threshold for a significant cluster. In addition, we further performed the GLM analysis to observe the association between the CT change and different clinical scores including HADS-A, HADS-D, hours of sleep and attacks per month.
To validate the significance of these network measures, we further applied support vector machine analysis to predict the clinical scores used in the diagnostic criteria for migraineurs. Here, we performed a support vector regressor (SVR) analysis – representing a machine-learning-based multiple regression method - that could associate the observed and trained values and present the regression coefficient for the accuracy of the prediction . The regression coefficient of 0.5 obtained after 10-fold cross validation is considered borderline significant result.
Structural network analysis
We found no significant difference (p > 0.05) in age and sex between the groups: EM - HC, CM - HC and EM - CM. Comparing networks obtained using graph theoretical framework between EM and HC, we found significantly (p < 0.05, FDR corrected) higher transitivity and assortativity in EM (Fig. 2). We further obtained the centrality measures, namely mean node and edge-betweenness to be higher in EM. For the contrast ‘CM – HC’, we found significantly higher modularity and assortativity in CM. We further obtained higher centrality measures (mean node and edge-betweenness) in CM. For the comparison of the networks between two patient groups, we found significantly higher modularity but lower clustering coefficient and transitivity in CM compared to EM.
When comparing similar networks obtained using values of subcortical volumes between EM and HC, we found higher transitivity and assortativity significantly in EM. The comparison CM and HC revealed significantly higher modularity, transitivity and assortativity in CM. All sub-cortical network results are shown in Fig. 3.
Among the cortical regions showing significant regional network difference in terms of clustering, degree and nodal edge-betweenness obtained using CT, NBS further revealed distinct networks with lower interregional connectivity for EM and CM when compared to HC. For EM, all 20 nodes (regions) showing the graph theory differences formed a network of significantly (p < 0.05, corrected) reduced connectivity in comparison to HC (Fig. 4a). However, for CM, out of 22 nodes (regions) only 19 formed a network of significantly (p < 0.05, corrected) reduced connectivity when compared to HC (Fig. 4b). For the comparison between two migraine groups EM and CM, NBS exhibited two distinct subnetworks for the contrast EM < CM and CM < EM (Fig. 5a, b). As expected, the contrast CM < EM showed lower structural connectivity between greater number of nodes than in contrast EM < CM (16 ROIs compared to 11 ROIs), indicating higher structural connectivity loss for chronic migraineurs than episodic.
Similarly, among the subcortical regions showing the graph theoretical difference in measures - clustering, degree and nodal edge-betweenness obtained using SCV, NBS revealed distinct networks for EM and CM with lower interregional connectivity when compared to HC. For EM, we revealed that out of 18 subcortical regions analyzed, 12 showed graph theoretical differences and formed a network with significantly reduced connectivity compared to HC (Fig. 4c). However, for CM, only 10 subcortical regions showed graph theoretical differences and exhibit a subnetwork with significantly reduced connectivity compared to HC (Fig. 4d). For the comparison of EM and CM, interestingly only the contrast EM < CM yielded marginally significant subnetwork with only four ROIs (Fig. 5c) and the contrast CM < EM was not significant.
Finally, the network measures (assortativity, transitivity and modularity) obtained using both cortical thickness and subcortical volumes showing a significant difference between the groups, additionally yielded a significant interrelation for MIDAS and attacks per month for the migraineurs (Fig. 6). Considering all network measures together, the SVR yield MIDAS (across EM and CM) with a regression coefficient of 0.792 for CT, 0.712 for SCV and attacks with 0.798 for CT, 0.715 for SCV. Using only assortativity as network measure, the association for MIDAS and attacks was with a regression coefficient of 0.779 for CT, 0.690 for SCV and 0.812 for CT, 0.798 for SCV respectively. Similarly, for transitivity – MIDAS, it was 0.655 for CT, 0.647 for SCV and for transitivity – attacks, it was 0.649 for CT and 0.589 for SCV. Modularity alone could reveal the interrelation to the MIDAS and attacks with regression coefficient of 0.656 for CT, 0.698 for SCV and 0.649 for CT, 0.735 for SCV respectively.
All reported CT group-differences and correlations were observed with age and sex as nuisance variables (p < 0.001, uncorrected as none of them survived multiple comparison (FDR) correction). Comparing EM to HC (supp. figure 1), average CT was significantly increased in EM in left lateral occipital cortex, supramarginal gyrus as well as in the right insula, lingual gyrus and precuneus. Additionally, the correlation of the average CT with HADS-A was significantly different in right precentral and inferior parietal cortex between EM and HC. Similarly, the correlation of the average CT with HADS-D was significantly different in left lingual gyrus and right supramarginal gyrus as well as the correlation of the average CT with hours of sleep in left superior parietal and right supramarginal gyrus and caudal middle frontal cortex between EM and HC.
The average CT was significantly increased in CM compared to HC in the left insula and posterior cingulate cortex (PCC) and significantly decreased in the bilateral inferior parietal and right lateral occipital cortex (supp. figure 2). Additionally, the correlation of the average CT with HADS-A was significantly different in right caudal anterior cingulate and precentral cortex between CM and HC. Similarly, the correlation of the average CT with HADS-D was significantly different in left lingual and PCC as well as the right supramarginal gyrus and precuneus between CM and HC. The correlation of the average CT with the hours of sleep differ significantly in left pars opercularis, superior parietal, insula and right lateral occipital and supramarginal gyrus.
The contrast ‘EM – CM’ (supp. figure 3) revealed significantly decreased average CT in EM in the left insula and significantly increased in right insula-, supramarginal- and postcentral gyrus. Additionally, the correlation of the average CT with HADS-A was significantly different in the bilateral inferior parietal and left superior parietal lobule. Similarly, the correlation of the average CT with HADS-D was significantly dissimilar in left insula between EM and CM. The correlation of the average CT with hours of sleep varied significantly in left superior parietal, insula and right supramarginal, postcentral and insula between EM and CM. Moreover, the correlation of the average CT with headache attacks per month was significantly different in left insula between these groups.
This study reported group differences in structural networks based on CT between HC, EM, and CM. Significant group differences were found in various graph theory measures and these measures were able to predict the clinical scores MIDAS and migraine attacks per month. Based on the results, we conclude that the impact of migraine severity leads to strong structural impairments and dysfunctional neural network configurations.
Structural network alterations in migraineurs
Using both GMV and rs-fMRI data for reconstructing structural and functional connectivity, Liu et al. demonstrated that EM showed abnormal global topology in both structural and functional networks, characterized by higher mean clustering coefficients  . Even though these findings are noteworthy in aiding to comprehend the pathology, studies investigating structural topological changes, which might have led to these functional alterations and provide a probable explanation of symptoms in different types of migraineurs, are still distant. In our study, EM displayed higher transitivity and assortativity along with the centrality measures than HC, indicative of a shift in the hubs for information transfer in this group. Similarly, higher modularity, assortativity and centrality measures (mean node and edge-betweenness) in CM than in HC furthers indicates the network being more assortative and segregated in CM. This finding is additionally supported by the fact that CM were found to have significantly higher modularity but lower clustering coefficient and transitivity compared to EM. As the NBS analysis yielded lower interregional connectivity in migraineurs, we conclude that migraineurs display disturbed connections not only in a localized brain area but also between regions of different anatomical compartments including the frontal, temporal, parietal and visual areas.
Importantly, as clinically relevant, we demonstrated the importance of these network measures by being able to interrelate the MIDAS scores and attacks per month for the migraineurs. The interdependence between attack frequency and network measures was highest when using assortativity followed by all network measures (assortativity, transitivity and modularity) together. For MIDAS, using all three measures or the assortativity lead to a similar association. Hence, from these and previous findings we could speculate that there might be as well a transition from episodic to chronic with the intensification of network segregation and assortativity in migraineurs.
Out of 18 analyzed subcortical regions, 12 and 10 for EM and CM respectively, revealed graph theoretical differences and formed a network with significantly reduced connectivity compared to HC. Strikingly, the right caudate, left thalamus, left accumbens area and right hippocampus showed stronger alterations in CM compared to EM. This finding indicates that chronic migraine impairs structural network integrity not only on the cortical but also on the subcortical level.
Brain morphometric alterations in migraineurs
Some studies reported CT increases in the somatosensory cortex  or higher visual brain regions, including V3A and MT+ , whereas others reported decreased CT in migraine patents with medication overuse headache  or in EM [25, 62, 63]. Our results indicate especially the CM showed a thicker CT compared to HC but that both EM and CM showed rather increases in CT, especially in visual brain regions (left occipital cortex, lingual gyrus). This could be related to the high presence of aura in our sample. Yet, Granziera et al., (2006) reported a CT increase in patients with and without aura in higher visual (V3A and MT+) brain regions. The CT increase could index altered excitability of the cortex in EM and even more so in CM. It has been shown that transcranial direct stimulation in EM can lead to reduced number of migraine days, indicating that neuronal excitation can be re-normalized resulting in lower migraine occurrence [64, 65]. Our novel finding of differentiating EM and CM based on association of CT and headache attacks observed in left insula underlines (a) the validity of CT as a structural marker to differentiate migraine subgroups and (b) the role of the salience network in migraine.
Functional alterations in CM and MOH
Alterations in rs-fMRI have been demonstrated in several studies [66,67,68,69,70,71,72]. For example, Lee et al. (2019) applied rs-fMRI and observed increased functional connectivity of the pain matrix in CM . Comparable to our findings, inter-regional coupling was altered in CM in a particular (pain) network comprising anterior insula, thalamus, prefrontal cortex, precuneus and anterior cingulate cortex. Using a seed-based approach, Schwedt et al. (2013) reported rs-fMRI connectivity differences between CM and controls in the anterior insula, amygdala, pulvinar, thalamus, middle temporal cortex, and periaqueductal gray (PAG) . Our study adds that not functional connectivity shows a correlation to migraine frequency  but also graph theoretical measures. Recently, our group applied a network and a seed based (PAG) approach and found abnormal rs-fMRI connectivity in patients with MOH compared to patients with myofascial pain and healthy controls . The alterations in MOH were seen as hyperconnectivity of the salience network (bilateral anterior insular cortex, dorsal anterior cingulate cortex, supplementary motor area), which correlated to white matter alterations in parts of this network (i.e. insular cortex), but also in the para-hippocampus, cerebellum and visual regions. As we observed alterations in structural connectivity in EM to HM in some of these regions (insula and anterior cingulate cortex), we cannot exclude that some of the observed effects related to CM might partially be mediated by the presence of MOH in some of the CM.
Some of the published rs-fMRI studies compared patients without aura to controls, and this sub-group comparison is not possible in our study, as the number of EM or CM without (and with) aura is too low. Yet, our SCV analysis extends rs-fMRI studies in CM as this group showed lower structural connectivity in EM and CM patients in comparison to HC in a network comprising amygdala, caudate nucleus, pallidum, thalamus, hippocampus and bilateral cerebellum. This observation is line to our previous study, which demonstrated structural (GMV) alterations in the basal ganglia (e.g., caudate nucleus and pallidum), hippocampus as well as in the thalamus, cerebellum, brainstem (PAG), visual and frontal cortex in patients with MOH . Another group reported increases in GMV in the right amygdala and right putamen in CM compared to controls, and headache frequency correlated positively with GMV in the putamen as well as in frontal and temporal regions .
In our study, CM (compared to EM), showed structural alterations with a right-hemispheric dominance. In addition, the study by Chen et al. (2016) reported right-dominant alterations in rs-fMRI connectivity, i.e. between the marginal division of neostriatum and the right middle temporal or right middle frontal gyrus . However, this was the case for EM, CM, and CM with MOH (compared to HC). Yet, a right dominance of GMV increased was also reported by a previous study in CM (some with MOH) .
It would have been interesting to examine the link between network alterations and patients’ disease duration, as done in a previous publication using resting-state fMRI . However, this parameter could not be reliably assessed as many of the patients could just remember that they had migraine for “some years” or “since their adolescence or early adulthood”. Our group size was moderate, especially for the group of CM patients. Yet, we found a systematic increase in impairment on structural integrity with the presence of migraine attacks, indicating that CM caused the strongest structural abnormalities compared to HC and EM. The findings from CT and subcortical volume analysis were uncorrected which might be because of the moderate sample size. Nonetheless, the results present a significant basis for the network measures observed and gives a more enhanced tool for understanding the migraine pathophysiology.
In conclusion, the level of impairment (migraine days per month) was associated with altered GVM but additionally with disturbed structural network integrity. The observation of an under-segregated network, especially in patients with CM, could be a sign of a maladaptive, elevated integration among pain-related brain circuits, leading to more excitability but less inhibition for the modulation of migraine.
Availability of data and materials
The datasets generated and/or analyzed during the current study are not publicly available due to patient consent but could be available from the corresponding author on reasonable request and would be decided upon individual basis.
Support Vector Regressions
Functional magnetic resonance imaging
Grey matter Volume
Diffusion Tensor Imaging
Migraine Disability Assessment
Hamilton Anxiety Score
Hamilton Depression Score
Magnetization Prepared Rapid Gradient Echo
Network Based Statistics
Field of View
Default Mode Network
Manzoni GC, Stovner LJ (2010) Epidemiology of headache. Handb Clin Neurol 97:3–22
Lipton RB, Bigal ME, Diamond M, Freitag F, Reed ML, Stewart WF et al (2007) Migraine prevalence, disease burden, and the need for preventive therapy. Neurology. 68(5):343–349
Vos T, Flaxman AD, Naghavi M, Lozano R, Michaud C, Ezzati M et al (2012) Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990-2010: a systematic analysis for the global burden of disease study 2010. Lancet. 380(9859):2163–2196
Headache Classification Committee of the International Headache Society (IHS) (2018) The International Classification of Headache Disorders, 3rd edition. Cephalalgia 38(1):1–211
Maleki N, Gollub RL (2016) What have we learned from brain functional connectivity studies in migraine headache? Headache. 56(3):453–461
Skorobogatykh K, van Hoogstraten WS, Degan D, Prischepa A, Savitskaya A, Ileen BM et al (2019) Functional connectivity studies in migraine: what have we learned? J Headache Pain 20(1):108
Xu G, Cheng S, Qu Y, Cheng Y, Zhou J, Li Z et al (2020) The functional alterations in primary migraine: a systematic review and meta-analysis protocol. Medicine (Baltimore) 99(10):e19019
Maleki N, Becerra L, Brawn J, Bigal M, Burstein R, Borsook D (2012) Concurrent functional and structural cortical alterations in migraine. Cephalalgia. 32(8):607–620
Rocca MA, Ceccarelli A, Falini A, Colombo B, Tortorella P, Bernasconi L et al (2006) Brain gray matter changes in migraine patients with T2-visible lesions: a 3-T MRI study. Stroke 37(7):1765–1770
Liu J, Lan L, Li G, Yan X, Nan J, Xiong S et al (2013) Migraine-related gray matter and white matter changes at a 1-year follow-up evaluation. J Pain 14(12):1703–1708
Schmidt-Wilcke T, Ganssbauer S, Neuner T, Bogdahn U, May A (2008) Subtle grey matter changes between migraine patients and healthy controls. Cephalalgia. 28(1):1–4
Lakhan SE, Avramut M, Tepper SJ (2013) Structural and functional neuroimaging in migraine: insights from 3 decades of research. Headache. 53(1):46–66
Valfre W, Rainero I, Bergui M, Pinessi L (2008) Voxel-based morphometry reveals gray matter abnormalities in migraine. Headache. 48(1):109–117
Jin C, Yuan K, Zhao L, Zhao L, Yu D, von Deneen KM et al (2013) Structural and functional abnormalities in migraine patients without aura. NMR Biomed 26(1):58–64
Mehnert J, May A. Functional and structural alterations in the migraine cerebellum. J Cereb Blood Flow Metab. 2017:271678X17722109
Matharu MS, Good CD, May A, Bahra A, Goadsby PJ (2003) No change in the structure of the brain in migraine: a voxel-based morphometric study. Eur J Neurol 10(1):53–57
Mehnert J, Schulte L, May A (2020) No grey matter alterations in longitudinal data of migraine patients. Brain 143(11):e93
Sheng L, Zhao P, Ma H, Yuan C, Zhong J, Dai Z et al (2020) A lack of consistent brain grey matter alterations in migraine. Brain 143(6):e45
Kim J, Suh SI, Seol H, Oh K, Seo WK, Yu SW et al (2008) Regional grey matter changes in patients with migraine: a voxel-based morphometry study. Cephalalgia. 28(6):598–604
Schmitz N, Admiraal-Behloul F, Arkink EB, Kruit MC, Schoonman GG, Ferrari MD et al (2008) Attack frequency and disease duration as indicators for brain damage in migraine. Headache. 48(7):1044–1055
Schmitz N, Arkink EB, Mulder M, Rubia K, Admiraal-Behloul F, Schoonman GG et al (2008) Frontal lobe structure and executive function in migraine patients. Neurosci Lett 440(2):92–96
DaSilva AF, Granziera C, Snyder J, Hadjikhani N (2007) Thickening in the somatosensory cortex of patients with migraine. Neurology. 69(21):1990–1995
Gaist D, Hougaard A, Garde E, Reislev NL, Wiwie R, Iversen P et al (2018) Migraine with visual aura associated with thicker visual cortex. Brain
Granziera C, DaSilva AF, Snyder J, Tuch DS, Hadjikhani N (2006) Anatomical alterations of the visual motion processing network in migraine with and without aura. PLoS Med 3(10):e402
Magon S, May A, Stankewitz A, Goadsby PJ, Schankin C, Ashina M et al (2019) Cortical abnormalities in episodic migraine: a multi-center 3T MRI study. Cephalalgia. 39(5):665–673
Hougaard A, Amin FM, Hoffmann MB, Larsson HB, Magon S, Sprenger T et al (2015) Structural gray matter abnormalities in migraine relate to headache lateralization, but not aura. Cephalalgia. 35(1):3–9
Schwedt TJ, Berisha V, Chong CD (2015) Temporal lobe cortical thickness correlations differentiate the migraine brain from the healthy brain. PLoS One 10(2):e0116687
Sporns O (2013) Structure and function of complex brain networks. Dialogues Clin Neurosci 15(3):247–262
Koirala N, Fleischer V, Glaser M, Zeuner KE, Deuschl G, Volkmann J et al (2018) Frontal lobe connectivity and network community characteristics are associated with the outcome of subthalamic nucleus deep brain stimulation in patients with Parkinson’s disease. Brain Topogr 31(2):311–321
Bassett DS, Sporns O (2017) Network neuroscience. Nat Neurosci 20(3):353–364
Liu J, Zhao L, Li G, Xiong S, Nan J, Li J et al (2012) Hierarchical alteration of brain structural and functional networks in female migraine sufferers. PLoS One 7(12):e51250
Li K, Liu L, Yin Q, Dun W, Xu X, Liu J et al (2017) Abnormal rich club organization and impaired correlation between structural and functional connectivity in migraine sufferers. Brain Imaging Behav 11(2):526–540
Michels L, Villanueva J, O'Gorman RL, Muthuramam M, Koirala N, Buechler R et al (2019) Interictal hyperperfusion in the higher visual cortex in patients with episodic migraine. Headache 59(10):1808–1820 Accepted
Hodkinson DJ, Veggeberg R, Wilcox SL, Scrivani S, Burstein R, Becerra L et al (2015) Primary somatosensory cortices contain altered patterns of regional cerebral blood flow in the Interictal phase of migraine. PLoS One 10(9):e0137971
Headache Classification Committee of the International Headache S (2013) The international classification of headache disorders, 3rd edition (beta version). Cephalalgia. 33(9):629–808
May A, Schulte LH (2016) Chronic migraine: risk factors, mechanisms and treatment. Nat Rev Neurol 12(8):455–464
Dodick DW, Turkel CC, DeGryse RE, Aurora SK, Silberstein SD, Lipton RB et al (2010) OnabotulinumtoxinA for treatment of chronic migraine: pooled results from the double-blind, randomized, placebo-controlled phases of the PREEMPT clinical program. Headache. 50(6):921–936
Stewart WF, Lipton RB, Dowson AJ, Sawyer J (2001) Development and testing of the migraine disability assessment (MIDAS) questionnaire to assess headache-related disability. Neurology. 56(6 Suppl 1):S20–S28
Zigmond AS, Snaith RP (1983) The hospital anxiety and depression scale. Acta Psychiatr Scand 67(6):361–370
Dale AM, Fischl B, Sereno MI (1999) Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage. 9(2):179–194
Fischl B, Dale AM (2000) Measuring the thickness of the human cerebral cortex from magnetic resonance images. Proc Natl Acad Sci 97(20):11050–11055
Fischl B, Salat DH, Busa E, Albert M, Dieterich M, Haselgrove C et al (2002) Whole brain segmentation: automated labeling of neuroanatomical structures in the human brain. Neuron. 33(3):341–355
Rosas HD, Liu AK, Hersch S, Glessner M, Ferrante RJ, Salat DH et al (2002) Regional and progressive thinning of the cortical ribbon in Huntington's disease. Neurology. 58(5):695–701
Kuperberg GR, Broome MR, McGuire PK, David AS, Eddy M, Ozawa F et al (2003) Regionally localized thinning of the cerebral cortex in schizophrenia. Arch Gen Psychiatry 60(9):878–888
Reuter M, Schmansky NJ, Rosas HD, Fischl B (2012) Within-subject template estimation for unbiased longitudinal image analysis. Neuroimage. 61(4):1402–1418
Rubinov M, Sporns O (2010) Complex network measures of brain connectivity: uses and interpretations. Neuroimage. 52(3):1059–1069
Hosseini SM, Hoeft F, Kesler SR (2012) GAT: a graph-theoretical analysis toolbox for analyzing between-group differences in large-scale structural and functional brain networks. PLoS One 7(7):e40709
Fleischer V, Radetz A, Ciolac D, Muthuraman M, Gonzalez-Escamilla G, Zipp F et al (2019) Graph theoretical framework of brain networks in multiple sclerosis: a review of concepts. Neuroscience 403:35–53. https://doi.org/10.1016/j.neuroscience.2017.10.033
Fleischer V, Koirala N, Droby A, Gracien R-M, Deichmann R, Ziemann U et al (2019) Longitudinal cortical network reorganization in early relapsing–remitting multiple sclerosis. Ther Adv Neurol Disord 12:1756286419838673
Fleischer V, Groger A, Koirala N, Droby A, Muthuraman M, Kolber P et al (2017) Increased structural white and grey matter network connectivity compensates for functional decline in early multiple sclerosis. Mult Scler 23(3):432–441
Blondel VDG, Jean-Loup, Lambiotte R, Lefebvre E (2008) Fast unfolding of communities in large networks Journal of Statistical Mechanics. Theory Exp 10:10008
Zalesky A, Fornito A, Bullmore ET (2010) Network-based statistic: identifying differences in brain networks. Neuroimage. 53(4):1197–1207
Koirala N, Anwar AR, Ciolac D, Glaser M, Pintea B, Deuschl G et al (2019) Alterations in white matter network and microstructural integrity differentiate Parkinson's disease patients and healthy subjects. Front Aging Neurosci 11:191
Korgaonkar MS, Fornito A, Williams LM, Grieve SM (2014) Abnormal structural networks characterize major depressive disorder: a connectome analysis. Biol Psychiatry 76(7):567–574
Desikan RS, Ségonne F, Fischl B, Quinn BT, Dickerson BC, Blacker D et al (2006) An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage. 31(3):968–980
Bassett DS, Bullmore E, Verchinski BA, Mattay VS, Weinberger DR, Meyer-Lindenberg A (2008) Hierarchical organization of human cortical networks in health and schizophrenia. J Neurosci 28(37):9239–9248
He Y, Chen Z, Evans A (2008) Structural insights into aberrant topological patterns of large-scale cortical networks in Alzheimer's disease. J Neurosci 28(18):4756–4766
Achard S, Bullmore E (2007) Efficiency and cost of economical brain functional networks. PLoS Comput Biol 3(2):e17
Bernhardt BC, Chen Z, He Y, Evans AC, Bernasconi N (2011) Graph-theoretical analysis reveals disrupted small-world organization of cortical thickness correlation networks in temporal lobe epilepsy. Cereb Cortex 21(9):2147–2157
Drucker H, Burges CJC, Kaufman L, Smola A, Vapnik V (1996) Support vector regression machines. Proceedings of the 9th International Conference on Neural Information Processing Systems. MIT Press, Denver, Colorado, pp 155–161
Riederer F, Schaer M, Gantenbein AR, Luechinger R, Michels L, Kaya M et al (2017) Cortical Alterations in Medication-Overuse Headache. Headache 57(2):255–265. https://doi.org/10.1111/head.12993. Epub 2016 Dec 2
Messina R, Rocca MA, Colombo B, Valsasina P, Horsfield MA, Copetti M et al (2013) Cortical abnormalities in patients with migraine: a surface-based analysis. Radiology. 268(1):170–180
Seifert CL, Magon S, Staehle K, Zimmer C, Foerschler A, Radue EW et al (2012) A case-control study on cortical thickness in episodic cluster headache. Headache. 52(9):1362–1368
Vigano A, D'Elia TS, Sava SL, Auve M, De Pasqua V, Colosimo A et al (2013) Transcranial direct current stimulation (tDCS) of the visual cortex: a proof-of-concept study based on interictal electrophysiological abnormalities in migraine. J Headache Pain 14:23
Pohl H, Moisa M, Jung HH, Brenner K, Aschmann J, Riederer F et al (2020) Long-term effects of self-administered Transcranial direct current stimulation in episodic migraine prevention: results of a randomized controlled trial. Neuromodulation. https://doi.org/10.1111/ner.13292. Online ahead of print
Coppola G, Di Renzo A, Petolicchio B, Tinelli E, Di Lorenzo C, Serrao M et al (2020) Increased neural connectivity between the hypothalamus and cortical resting-state functional networks in chronic migraine. J Neurol 267(1):185–191. https://doi.org/10.1007/s00415-019-09571-y. Epub 2019 Oct 12
Chen Z, Chen X, Liu M, Liu S, Shu S, Ma L et al (2016) Altered functional connectivity of the marginal division in migraine: a resting-state fMRI study. J Headache Pain 17(1):89
Liu J, Zhao L, Lei F, Zhang Y, Yuan K, Gong Q et al (2015) Disrupted resting-state functional connectivity and its changing trend in migraine suffers. Hum Brain Mapp 36(5):1892–1907
Hubbard CS, Khan SA, Keaser ML, Mathur VA, Goyal M, Seminowicz DA (2014) Altered Brain Structure and Function Correlate with Disease Severity and Pain Catastrophizing in Migraine Patients. eNeuro 1(1):e20 14
Schwedt TJ, Schlaggar BL, Mar S, Nolan T, Coalson RS, Nardos B et al (2013) Atypical resting-state functional connectivity of affective pain regions in chronic migraine. Headache. 53(5):737–751
Hadjikhani N, Ward N, Boshyan J, Napadow V, Maeda Y, Truini A et al (2013) The missing link: enhanced functional connectivity between amygdala and visceroceptive cortex in migraine. Cephalalgia. 33(15):1264–1268
Michels L, Christidi F, Steiger VR, Sandor PS, Gantenbein AR, Landmann G et al (2017) Pain modulation is affected differently in medication-overuse headache and chronic myofascial pain - a multimodal MRI study. Cephalalgia 37(8):764–779. https://doi.org/10.1177/0333102416652625. Epub 2016 June 1
Lee MJ, Park BY, Cho S, Kim ST, Park H, Chung CS (2019) Increased connectivity of pain matrix in chronic migraine: a resting-state functional MRI study. J Headache Pain 20(1):29
Riederer F, Marti M, Luechinger R, Lanzenberger R, von Meyenburg J, Gantenbein AR et al (2012) Grey matter changes associated with medication-overuse headache: correlations with disease related disability and anxiety. World J Biol Psychiatry 13(7):517–525
Neeb L, Bastian K, Villringer K, Israel H, Reuter U, Fiebach JB (2017) Structural Gray matter alterations in chronic migraine: implications for a progressive disease? Headache. 57(3):400–416
We thank all subjects for the participation in this study. We thank Catharina Fritz-Rochner for help with patient recruitment and data analysis.
We greatly appreciate the financial support by the Swiss Headache Society (Hansruedi Isler Forschungsstipendium). SG and MM acknowledge the support from SFB-TR-128.
Ethics approval and consent to participate
The study was approved by the ethics committee of canton Zurich (KEK number E-37/2007), Switzerland. All subjects provided written informed consent prior to study enrolment.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Figure S1. Illustration of between-group differences for CT. (A) Regions depicting significant difference in the CT between EM and HC. Regions depicting significant difference between EM and HC in the correlation of average CT with HADS-A (B), HADS-D (C), and hours of sleep (D). A summary of the significant regions is reported in supplementary Table 1. All results are shown at p < 0.001 (uncorrected).
: Figure S2. Chronic migraine patients (CM) and healthy controls (HC): Regions depicting significant difference in the CT between CM and HC (A). Regions depicting significant difference between CM and HC in the correlation of average CT with HADS-A (B), HADS-D (C), and hours of sleep (D). A summary of the significant regions is reported in supplementary Table 2. All results are shown at p < 0.001 (uncorrected).
: Figure S3. Episodic migraine patients (EM) and chronic migraine patients (CM): Regions depicting significant difference in the CT between EM and CM (A). Regions depicting significant difference between EM and CM in the regression of average CT with HADS-A (B), HADS-D (C), hours of sleep (D) and number of headache attacks per month (E). A summary of the significant regions is reported in Supplementary Table 3. All results are shown at p < 0.001 (uncorrected).
About this article
Cite this article
Michels, L., Koirala, N., Groppa, S. et al. Structural brain network characteristics in patients with episodic and chronic migraine. J Headache Pain 22, 8 (2021). https://doi.org/10.1186/s10194-021-01216-8
- Graph theory