
==== Front
J Headache Pain
J Headache Pain
The Journal of Headache and Pain
1129-2369
1129-2377
Springer Milan Milan

1854
10.1186/s10194-024-01854-8
Research
Involvement of the cerebellum in structural connectivity enhancement in episodic migraine
Matoso Ana anamatoso@tecnico.ulisboa.pt

1
Fouto Ana R. 1
Esteves Inês 1
Ruiz-Tagle Amparo 1
Caetano Gina 1
da Silva Nuno A. 2
Vilela Pedro 3
Gil-Gouveia Raquel 45
Nunes Rita G. 1
Figueiredo Patrícia 1
1 grid.9983.b 0000 0001 2181 4263 Institute for Systems and Robotics – Lisboa and Department of Bioengineering, Instituto Superior Técnico, University of Lisbon, Lisbon, Portugal
2 https://ror.org/03jpm9j23 grid.414429.e 0000 0001 0163 5700 Learning Health, Hospital da Luz, Lisbon, Portugal
3 https://ror.org/03jpm9j23 grid.414429.e 0000 0001 0163 5700 Imaging Department, Hospital da Luz, Lisbon, Portugal
4 https://ror.org/03jpm9j23 grid.414429.e 0000 0001 0163 5700 Neurology Department, Hospital da Luz, Lisbon, Portugal
5 https://ror.org/03b9snr86 grid.7831.d 0000 0001 0410 653X Center for Interdisciplinary Research in Health, Universidade Católica Portuguesa, Lisbon, Portugal
18 9 2024
18 9 2024
2024
25 1 15415 7 2024
31 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Background

The pathophysiology of migraine remains poorly understood, yet a growing number of studies have shown structural connectivity disruptions across large-scale brain networks. Although both structural and functional changes have been found in the cerebellum of migraine patients, the cerebellum has barely been assessed in previous structural connectivity studies of migraine. Our objective is to investigate the structural connectivity of the entire brain, including the cerebellum, in individuals diagnosed with episodic migraine without aura during the interictal phase, compared with healthy controls.

Methods

To that end, 14 migraine patients and 15 healthy controls were recruited (all female), and diffusion-weighted and T1-weighted MRI data were acquired. The structural connectome was estimated for each participant based on two different whole-brain parcellations, including cortical and subcortical regions as well as the cerebellum. The structural connectivity patterns, as well as global and local graph theory metrics, were compared between patients and controls, for each of the two parcellations, using network-based statistics and a generalized linear model (GLM), respectively. We also compared the number of connectome streamlines within specific white matter tracts using a GLM.

Results

We found increased structural connectivity in migraine patients relative to healthy controls with a distinct involvement of cerebellar regions, using both parcellations. Specifically, the node degree of the posterior lobe of the cerebellum was greater in patients than in controls and patients presented a higher number of streamlines within the anterior limb of the internal capsule. Moreover, the connectomes of patients exhibited greater global efficiency and shorter characteristic path length, which correlated with the age onset of migraine.

Conclusions

A distinctive pattern of heightened structural connectivity and enhanced global efficiency in migraine patients compared to controls was identified, which distinctively involves the cerebellum. These findings provide evidence for increased integration within structural brain networks in migraine and underscore the significance of the cerebellum in migraine pathophysiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s10194-024-01854-8.

Keywords

Migraine
MRI
Structural connectivity
Cerebellum
Tractography
http://dx.doi.org/10.13039/501100001871 Fundação para a Ciência e a Tecnologia 2023.03810.BDANA SFRH/BD/139561/2018 PTDC/EMD-EMD/29675/2017 LARSyS - FCT10.54499/LA/P/0083/2020 issue-copyright-statement© Springer-Verlag Italia S.r.l., part of Springer Nature 2024
==== Body
pmcIntroduction

About 14% of people suffer from migraine, making it one of the most prevalent neurological diseases worldwide [1]. Although its pathophysiology is still not fully understood, multiple studies have found structural [2–11] and functional [2, 10–14] abnormalities throughout extensive brain networks. In particular, the cerebellum has an important role in pain processing, and it has been shown to be altered in numerous neuroimaging studies of migraine [14]. Mehnert et al. [15] reported gray matter volume changes in the crus (I and II), while Qin et al. [16] found microstructural changes in the vermis VI and in the bilateral lobules V and VI of the cerebellum. Evidence also suggests that migraine patients are more prone to ischemic lesions in the posterior lobe of the cerebellum [17]. Importantly, two studies [12, 13] have reported functional connectivity changes involving the crus, as well as the insula and the hippocampus, namely in terms of the network measures of centrality and efficiency. Taken together, these studies highlight the significant role of the cerebellum in the structural and functional brain changes observed in patients with episodic migraine.

Regarding structural connectivity, several studies have reported alterations in migraine based on diffusion Magnetic Resonance Imaging (dMRI) parameters [4, 5, 8–10, 18]. Some studies found increased structural connectivity between cortical regions involved in pain perception and processing [8], and between subcortical regions such as the thalamus and the caudate nucleus [4], suggesting that subcortical networks are strengthened in migraine [4, 5]. Other studies investigated connectome differences from a graph theory perspective using the connectivity between regions as edges and the regions as nodes. Li et al. [10] found a decreased characteristic path length in the structural connectome of migraine patients, which points to an increased number of short-distance connections and hence an increased integration. Similarly, increased integration was found in subcortical regions such as the putamen, the pallidum, and the thalamus as well as in the parahippocampal gyrus, and the anterior cingulate gyrus [3, 10]. Although the cerebellum has an important role in pain processing, and it has been shown to be altered in numerous neuroimaging studies of migraine [14], to our knowledge no study of structural connectivity using tractography in migraine has included it.

Here, we aim to investigate structural connectivity changes in migraine, considering the whole brain including cortical and subcortical regions as well as the cerebellum. It is important to note that some of the studies mentioned above have heterogeneous cohorts, including male and female patients, with and without aura, and with chronic and episodic migraine [4–6, 15]. This makes it difficult to disentangle the specific connectivity patterns for each subtype of migraine, as it is very likely that different subtypes have different brain functional and connectivity patterns. In our study, we focused solely on low-frequency episodic migraine without aura in female patients, which is by far the most prevalent group of migraine patients [19, 20].

Methods

Participants

This study analyzed data from the MigN2treat cohort (https://welcome.isr.tecnico.ulisboa.pt/projects/multimodal-neuroimaging-biomarkers-throughout-the-migraine-cycle-towards-neurofeedback-training-for-personalized-anti-migraine-treatment/) cohort which comprised 14 female patients diagnosed with episodic migraine without aura according to the criteria of the 3rd edition of the International Classification of Headache Disorders (ICHD-III) [21], and 15 healthy controls. Both healthy controls and patients verified the following inclusion criteria: (a) age between 18 and 55 years; (b) had at least 9 years of education; (c) had Portuguese as their first language; and the following exclusion criteria: (a) diagnosis of a neurological condition (other than migraine for patients); (b) diagnosis of any psychiatric disorder (severe anxiety or depressive symptoms were excluded by the use of standardized self-report scales State-Trait Anxiety Inventory and the Zung Self-Rating Depression scale); (c) daily use of psychoactive medication including migraine prophylaxis for the patient group; (d) being pregnant, breastfeeding, post-menopause, or using of contraception precluding cyclic menses; (e) contraindications for performing an MRI (e.g. being claustrophobic); and (f) evidence of incident brain lesion or structural abnormalities on the T1-weighted MRI.

Healthy controls were matched for gender, contraceptive use, and menstrual phase at the time of scanning and balanced for age. All participants provided written informed consent and the study was carried out following the Declaration of Helsinki upon approval by the local Ethics Committee.

Imaging data was acquired during the post-ovulation period, corresponding to the interictal phase, as participants in the cohort experienced menstrual-related migraine attacks. Prior to scanning, participants were required to be free from pain for at least 48 h, with confirmation of the absence of a migraine attack obtained 72 h post-scan.

The following demographic and clinical parameters were collected: age (years), frequency (migraine attacks per month), age of onset of migraine (years), and disease duration (years).

Data acquisition

The MRI data was acquired using a 3 Tesla system (Siemens Vida), equipped with a 64-channel receiver head coil from June 2019 to November 2022. Regarding the dMRI data, an echo planar imaging (EPI) sequence was used with the following parameters: Time of Repetition (TR) = 6800 ms, Time to Echo (TE) = 89 ms, 66 axial slices with an in-plane Generalized Autocalibrating Partially Parallel Acquisitions (GRAPPA) factor of 2, Simultaneous Multi-slice (SMS) factor of 3, flip angle = 90º, and 2 mm isotropic voxel size. The sampling scheme consisted of 3 diffusion shells (multi-shell) with b = 400, 1000, 2000 s/mm2 along 32, 32, 60 unique gradient directions, respectively, 8 non-diffusion-weighted volumes (b0) and an additional 3 b0 images with opposite phase encoding (Posterior-Anterior). The total imaging time was 15 min and 47 s. One T1-weighted image was for each subject also acquired with a 3D magnetization-prepared rapid gradient echo (MPRAGE) sequence with the following parameters: 1 mm isotropic voxel size, TR = 2300 ms, TE = 2.98 ms, Time of Inversion (TI) = 900 ms, flip angle = 9º, FOV = 256 × 240 × 128 mm3 having lasted 5 min and 12 s. A Fluid-Attenuated Inversion Recovery (FLAIR) image was also acquired and assessed for white matter lesions using the following parameters: 0.9 mm isotropic resolution, TR = 5000 ms, TE = 386 ms, TI = 1800 ms, (FOV) = 240 × 240 × 151 mm3, GRAPPA = 2, and lasted 5 min and 57 s, with 0.9 mm isotropic resolution. Both anatomical scans were evaluated by an experienced neuroradiologist.

Data preprocessing

An overview of the image processing and subsequent analyses can be found in Fig. 1.

Fig. 1 Overview of the data processing pipeline, from the acquisition of the images to the connectome generation, including the different analyses performed

Image processing was conducted using a combination of FSL [22] (version 6.0.5) and MRtrix [23] (version 3.0.3) tools. The following preprocessing procedures were applied to the dMRI data based on the DESIGNER (Diffusion parameter EStImation with Gibbs and NoisE Removal) pipeline [24]: denoising (dwidenoise), Gibbs ringing correction (mrdegibbs), and Rician bias correction were performed according to [25], followed by geometric and eddy-current distortions and motion correction with FSL (eddy and topup). After this, bias field correction was done with MRtrix (dwibiascorrect using the -ants option).

Tractography

To calculate the streamlines (tracts) of white matter fibers, MRtrix was used. Firstly, the basis functions for each tissue type were estimated from the subject’s dMRI data. Then, the fiber orientation density functions (FODf) were calculated using multi-shell multi-tissue constrained spherical deconvolution [26] followed by their normalization to enable the comparison between subjects. Afterward, Anatomically-Constrained Tractography (ACT) [27] was performed which makes use of prior anatomical knowledge to restrict the generated streamlines by, for instance, preventing streamlines from ending in implausible places such as in the middle of white matter or in the cerebrospinal fluid (CSF). This is done by segmenting the brain into five tissue types (5TT) using the 5ttgen function of MRtrix and then using the grey and white matter masks to create a seed region along the interface of these two regions from which seed points were randomly generated. From these points, the streamlines were reconstructed thus creating a tractogram. Ten million streamlines were seeded with a maximum length of 25 cm using the iFOD2 tracking algorithm (probabilistic) [28]. Subsequently, spherical deconvolution-informed filtering of tractograms (SIFT2) [29] was performed. Finally, the tractogram was utilized to determine the connectivity between different regions of interest (ROI) and to generate a connectome.

Connectome generation

Two sets of ROIs covering the whole brain were defined using two different parcellations: (i) Schaefer + SC + CB parcellation: 100 cortical regions from the Schaefer atlas [30] combined with 12 subcortical regions (SC) and 26 cerebellum (CB) regions from the Automated Anatomical Labeling atlas with 116 nodes (AAL116) [31]; and (ii) AAL116 parcellation: the 116 regions of the AAL116 atlas, which includes cortical regions as well as the same SC and CB regions as in i).

For each of these parcellations, a connectivity matrix was calculated based on the streamlines that connect each pair of ROIs and the volumes of those ROIs, according to:

1 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{a}_{ij}=\frac{2}{{V}_{i}+{V}_{j}}{\:n}_{ij}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{a}_{ij}$$\end{document} is the entry of the connectivity matrix concerning ROI i and ROI j, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{n}_{ij}$$\end{document} is the number of streamlines between the two ROIs given by the tractogram, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{V}_{i}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{V}_{j}$$\end{document} are the volumes of ROI i and ROI j, respectively. Note that, due to the absence of directionality in dMRI, the connectivity matrix is symmetric. Also, self-connections (streamlines that begin and end in the same ROI) were not considered, and as such the connectivity matrix has a zero diagonal.

Connectome analysis

The analyses present in this section were all performed using MATLAB version: 9.14.0 (R2023a).

Connectivity analysis

To detect significant connectivity changes between patients and controls, the network-based statistics toolbox (NBS version 1.2) [32] was used, which performs mass-univariate testing at every connection evaluating the null hypothesis at the level of interconnected subnetworks instead of individual connections. It leverages permutation testing to assess the statistical significance across subnetworks of connections. The default value of 3.1 was used as the test statistic’s lower threshold to apply to the network so that only the more significant connections survive. NBS then examines whether these significant connections form larger clusters/subnetworks that are unlikely to have occurred by chance by performing a 5000 permutations (default value) testing. Clusters whose family-wise error rate (FWER) corrected p-value is lower than 0.05 were considered statistically significant. Age was used as a covariate.

For a more straightforward interpretation of the results, the regions (nodes) of each parcellation were grouped into larger regions (node aggregates) according to the larger-scale organization of each atlas: (i) Schaefer + SC + CB parcellation: the 7 canonical resting-state functional networks identified by Yeo [33] + subcortical regions + crus + posterior lobe of the cerebellum (PLC) + vermis; and (ii) AAL116 parcellation: 4 lobes (occipital, parietal, frontal and temporal) + subcortical regions + crus + PLC + vermis. For all regions but the vermis, left and right regions were considered separately.

Graph metrics analysis

The following graph metrics were calculated on the node aggregates described above using the Brain Connectivity Toolbox (version 2019) [34]: global metrics (characteristic path length, global efficiency, average degree, clustering coefficient, and small-worldness) and a local metric (node degree). For the comparisons between patients and controls, a Generalized Linear Model (GLM) was fitted to the data, including group as the main effect and age as a covariate. In the patients, the correlation between the graph metrics and the clinical parameters was calculated using the Spearman correlation coefficient. For both the group comparisons and the correlations, a p-value lower than 0.05 was considered significant for the global metrics, whilst a p-value lower than 0.05 corrected for the number of nodes using Bonferroni correction was considered significant for the nodal metrics (p < 0.0024 and p < 0.0033, for the Schaefer + SC + CB parcellation and AAL116 parcellation, respectively).

White matter tract analysis

Additionally, using the tractogram of each subject, the number of reconstructed streamlines was computed within each ROI of a white matter atlas (ICBM-DTI-81 white-matter labels atlas [35]). A GLM was fitted to the data, which included the group as the main effect and age as a covariate. A p-value lower than 0.001 (corrected for the number of ROIs of the atlas using Bonferroni correction) was considered significant.

Results

Study population and MRI data

Concerning the clinical and demographic data collected about the subjects, Table 1 presents the statistical analysis performed. No significant changes were found between patients and controls concerning their age (p = 0.08).

Table 1 Clinical and demographic characteristics of migraine patients (M) and healthy controls (HC). No significant changes were found between patients and controls concerning their age (p = 0.08). ✝Mann-Whitney U test

Parameter	Group	Mean ± SD	P-value	
Age (years)	M	35.5 ± 8.7	0.08✝	
HC	30.8 ± 6.8	
Frequency (migraine attacks/month)	M	2.4 ± 1.9	-	
Age Onset (Years)	M	15.7 ± 7.3	-	
Pain Intensity (0–10)	M	6.62 ± 1.31	-	
Disease Duration (Years)	M	20.1 ± 11.0	-	

Furthermore, no relevant white matter lesions were found in the anatomical scans.

An example of a tractogram and the corresponding connectivity matrices, obtained from a representative patient with each parcellation, are shown in Figure S1 of the Supplementary Material.

Connectivity analysis

Several connections exhibited significantly increased connectivity in patients relative to controls, for both parcellations (listed in Table S1 and Table S2 in the Supplementary Material). Specifically, a network with increased connectivity in patients was found in both parcellations with p = 0.03 (corrected for FWER) for the Schaefer + SC + CB parcellation and p = 0.04 (corrected for FWER) for the AAL116 parcellation.

For the visualization of significant connectome differences, we display the total number of connections exhibiting significant differences between patients and controls within each node aggregate in Fig. 2, both for the Schaefer + SC + CB parcellation (left) and the AAL116 parcellation (right).

Fig. 2 Connectograms of connections exhibiting significantly increased connectivity in patients (M) vs. controls (HC), detected by NBS using both the Schaefer + SC + CB parcellation (M > HC, p = 0.03) and the AAL116 parcellation (M > HC, p = 0.04). The edges represent the sum of connections with significant differences within each region (node aggregate). L = Left, R = Right, VN = Visual Network, SMN = Somatossensory Network, FPN = Fronto-Parietal Network, DAN = Dorsal Attention Network, VAN = Ventral Attention Network, DMN = Default Mode Network, LN = Limbic Network, Occ = Occipital, Frt = Frontal, Par = Parietal, Tmp = Temporal, SC = Subcortical, PLC = Posterior Lobe of the Cerebellum

Despite some differences between the results obtained with the two parcellations, common changes in connectivity patterns exist. Specifically, increased connectivity in patients was consistently observed with both parcellations between the right crus and the left posterior lobe of the cerebellum, as well as between the frontal and parietal lobe network and the vermis. Additionally, increased connectivity was also evident between various cerebellar regions and occipital regions (vision network). In summary, a distinctive pattern of increased cerebellar connectivity was consistently observed in patients relative to controls, using both parcellations.

No network was found to have decreased connectivity in migraine patients.

Graph metrics analysis

Concerning local graph metrics, the only significant between-groups difference was found in the node degree of the right and left posterior lobe of the cerebellum, which was increased in patients relative to controls using both parcellations, as shown in Fig. 3. The results for the global graph metrics are presented in Fig. 4. A significant decrease of the characteristic path length and a significant increase of the global efficiency were observed in patients relative to controls. The average degree was significantly higher in patients compared to controls, using both parcellations, while the clustering coefficient and the small-worldness showed no significant differences between groups. The p-values obtained for each metric are presented in Table S3 of the Supplementary Material.

Fig. 3 Local graph metrics that show significant differences between patients and controls for both parcellations. The boxplots represent distributions of the metrics across subjects. Significant differences between groups are indicated (*p < 0.05, **p < 0.01 corrected for multiple comparisons). The node degree of the right and left Posterior Lobes of the Cerebellum (PLC) is increased in migraine patients relative to healthy controls (M > HC)

Fig. 4 Global graph metrics in patients and controls for both parcellations. The boxplots represent distributions across subjects. Significant differences between groups are indicated (*p < 0.05, **p < 0.01, ***p < 0.001). The characteristic path length was decreased in patients relative to controls in both parcellations, whilst the global efficiency and the average degree were increased. No significant changes were found in the clustering coefficient and the small-worldness

Moreover, a significant correlation (Spearman rank correlation r2 = 0.26, p = 0.038) was found between the characteristic path length obtained using the Schaefer + SC + CB parcellation in patients and their age onset of migraine (Fig. 5). No other significant correlations between graph metrics and clinical parameters were found.

Fig. 5 Correlation between the characteristic path length (from the connectome obtained with the Schaefer + SC + CB parcellation) and the age onset of migraine for the patients (Spearman rank correlation r2 = 0.26 with p = 0.038)

White matter tract analysis

Regarding the change in the number of reconstructed white matter streamlines, although not statistically significant, there is a slight bilateral increase in the anterior limb of the internal capsule (ALIC) in the patient group relative to controls, as shown in Fig. 6.

Fig. 6 Number of reconstructed streamlines in the left and right anterior limb of the internal capsule (ALIC) (from the ICBM-DTI-81 white-matter labels atlas) for patients (M) and controls (HC). The boxplots represent distributions of the metrics across subjects. No significant differences between groups are indicated. In both ROIs, there was a slight increase in the number of reconstructed streamlines in patients

Discussion

In this study, the structural connectome of a group of migraine patients with episodic migraine without aura was studied using dMRI and tractography considering two different whole-brain parcellations. We identified increased connectivity of cerebellar regions in migraine patients, which was reflected in the increased node degree of the posterior lobe of the cerebellum. Consistently with such increased connectivity, we also find an increased number of streamlines in the anterior limb of the internal capsule. From a global connectome perspective, we found increased global efficiency and decreased characteristic path length in patients, which was correlated with the age onset of migraine.

Connectivity analysis

The overall patterns of increased connectivity in patients relative to controls were characterized by a strong involvement of cerebellar regions, namely the crus, the posterior lobe of the cerebellum, and the vermis.

The connectivity analysis showed increased connectivity between the crus and the posterior lobe of the cerebellum. It is believed that the crus exhibits increased activity in response to painful stimuli and has cognitive and emotional representations [36]. Since the crus and the posterior lobe of the cerebellum are involved in cognitive processes [37], alterations in their connectivity may be associated with cognitive impairments reported in migraine patients [38–40]. In addition, stimulation of the posterior lobe of the cerebellum has been shown to modulate the response to noxious visceral stimuli in mice [41], which could be related to migraine pathology.

We also found increased connectivity between the vermis and other regions. The posterior vermis is the anatomical substrate of the limbic cerebellum [37] and thus there can be a disruption in the cerebro-cerebellar-limbic loops of cerebellar input that could be indicative of the emergence of emotional disturbances like mood changes often observed in migraine patients [11, 39].

The cerebellum has an inhibitory role in pain processing [42], having several connections to the prefrontal cortex (via the thalamus) [43]. Therefore, the increased connectivity between the cerebellum and other regions could be indicative of a dysfunctional negative feedback loop in which the thalamus is not sensing the inhibitory signal [42]. These connectivity results could be indicative of the importance of the cerebellum in the modulation of migraine possibly through the modulation of calcitonin gene-related peptide (CGRP) [14, 15]. CGRP modulates nociception and assists in the onset of migraine, and thus it might play a role in reshaping connectivity and stabilizing synapses in the cerebellar circuitry as shown in early preclinical studies [44, 45]. In fact, many migraine-specific high-efficacy treatments target CGRP [46, 47], hence changes in the connectivity pattern in the cerebellum could be potentially explored both as potential biomarkers for treatment response or as a predictor of treatment efficacy.

Our results not only contribute further evidence that the cerebellum plays a central role in migraine pathophysiology but also emphasize the importance of including it as a region of interest in connectome studies of migraine, with the robustness of our analysis being supported by the observation of similar patterns across two different brain parcellations.

Graph metrics analysis

At the global connectome level, decreased characteristic path length and increased global efficiency were found in patients compared to controls, which can be associated with a dysfunctional modulation of pain networks in migraine and a faster dissemination of pain-related information [8–10]. Moreover, the negative correlation found between the characteristic path length of patients and their age onset could potentially imply a plastic adaptation to migraine over time (since the patient developed this disorder), which has been reported by several diffusion tensor imaging studies [4, 6, 48, 49] in migraine patients. These changes could also contribute to an increased pain perception as has been reported in functional studies [50, 51]. Additionally, some studies have also reported correlations between graph theory metrics or diffusion metrics and patients’ clinical data. For instance, Dai et al. [9] reported a correlation between local efficiency and the visual analog scale in episodic migraine without aura, while Planchuelo-Gómez et al. [4] and Chong et al. [6] reported a relationship between diffusion metrics such as fractional anisotropy (FA) and mean diffusivity and the years lived with migraine (disease duration), using less homogeneous cohorts than ours, suggesting that chronicity might be exacerbating neural abnormalities.

In terms of local graph metrics, we found a bilateral increase in the node degree of the posterior lobe of the cerebellum, which is consistent with the increased connectivity of this region with the rest of the brain. As the posterior lobe of the cerebellum has both cognitive and limbic functions [52], the increased node degree of this region might be a result of a compensatory mechanism leading to effects on cognitive processing and emotional state.

White matter tract analysis

Complementarily to the analysis of the structural connections between gray matter regions, the corresponding white matter tracts were also evaluated considering the number of streamlines that passed through each tract and found a slight increase in the number of streamlines bilaterally in the anterior limb of the internal capsule of patients. The internal capsule contains ascending and descending axons including fibers to and from the thalamus, and to and from the cerebellum [53]. This increase is therefore consistent with the increased connectivity of the cerebellum, and could also be indicative of the dysfunctional negative feedback loop mentioned before. On the other hand, the increased number of reconstructed fibers could also be a consequence of an increased FA in the cerebellum of patients, which other researchers [54] have found. Future studies should clarify these different possible reasons.

Limitations and innovative aspects

This study has a few limitations, the main one being the small sample size. Despite that, we focused on a homogeneous patient group with controlled conditions and selected a well-matched control group in terms of gender, and menstrual cycle phase, minimizing individual variability.

Furthermore, our connectivity metrics are based on the number of reconstructed streamlines obtained by using a specific tractography algorithm. Although this methodology is well-established and commonly used in the literature, it comes with several challenges in terms of quantification and interpretation. However, throughout the tractography pipeline, decisions were made to tackle these challenges. For instance, two of the main biases of this type of analysis are streamline termination biases, where the reconstructed streamlines end in unreasonable voxels (such as in the middle of CSF), and streamline quantification biases where the number of reconstructed fibers does not represent the actual fiber count for a given voxel [55, 56]. In this work, we employed the ACT framework that tackles the first bias by seeding the streamlines from the interface between grey and white matter, and we used a streamline filtering technique (SIFT2) to deal with the quantification of fiber density in each voxel. Currently, using the number of streamlines as a proxy for fiber count allied to these advanced post-processing techniques is the state-of-the-art for structural connectivity analysis generation [55]. Finally calculating local and global graph theory metrics provides a complete assessment of structural connectome. By using weighted connectivity matrices to represent the structural connectome instead of binary matrices, there is a better representation of biological properties [57].

This is the first study that uses tractography based on multi-shell data to estimate the structural connectome in migraine. Moreover, the inclusion of the cerebellum in two different brain parcellations is also novel. Finally, our study is unique by focusing on a homogeneous cohort of female-only patients with menstrual and menstrual-related low-frequency episodic migraine, and including healthy controls matched for the menstrual phase.

Conclusion

In conclusion, our findings reveal structural connectome alterations in a cohort of patients diagnosed with menstruation-related low-frequency episodic migraine without aura, prominently implicating the cerebellum. These results substantiate the notion of heightened integration across whole brain networks underscoring the previously acknowledged yet often overlooked role of the cerebellum in migraine pathophysiology.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Acknowledgements

We acknowledge the Portuguese Science Foundation and LARSyS for supporting our study. We also express our gratitude to the participants of this study as well as the collaborators from Hospital da Luz Lisboa.

Author contributions

AM performed the data analysis, created the figures and wrote the manuscript. AF, IE, AR and GC were responsible for the data acquisition and preprocessing and reviewed the manuscript. NS and PV provided the resources for the data acquisition and reviewed the manuscript. RG was responsible for funding acquisition, for the supervision of data acquisition and reviewed the manuscript. PF and RN were responsible for funding acquisition, they supervised the project and revised the manuscript. All authors read and approved the final manuscript.

Funding

This work was supported by the Portuguese Science Foundation through grants 2023.03810.BDANA, SFRH/BD/139561/2018, COVID/BD/153268/2023, PTDC/EMD-EMD/29675/2017, LISBOA-01-0145-FEDER-029675. This work was also supported by LARSyS funding (DOI: 10.54499/LA/P/0083/2020, 10.54499/UIDP/50009/2020, and 10.54499/UIDB/50009/2020).

Data availability

Data is available upon request. The code can be found in the following repository: https://github.com/anamatoso/connectivity-analysis-diffusion.

Declarations

Competing interests

The authors declare no competing interests.

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Stovner LJ Hagen K Linde M Steiner TJ The global prevalence of headache: an update, with analysis of the influences of methodological factors on prevalence estimates J Headache Pain 2022 23 1 17 10.1186/s10194-022-01402-2 34979902
Stovner LJ, Hagen K, Linde M, Steiner TJ (2022) The global prevalence of headache: an update, with analysis of the influences of methodological factors on prevalence estimates. J Headache Pain 23:1–17. 10.1186/s10194-022-01402-234979902
2. Ashina S Bentivegna E Martelletti P Eikermann-Haerter K Structural and functional brain changes in Migraine Pain Ther 2021 10 211 223 10.1007/s40122-021-00240-5 33594593
Ashina S, Bentivegna E, Martelletti P, Eikermann-Haerter K (2021) Structural and functional brain changes in Migraine. Pain Ther 10:211–223. 10.1007/s40122-021-00240-533594593
3. Liu J Ma S Mu J Integration of white matter network is associated with interindividual differences in psychologically mediated placebo response in migraine patients Hum Brain Mapp 2017 38 5250 5259 10.1002/hbm.23729 28731567
Liu J, Ma S, Mu J et al (2017) Integration of white matter network is associated with interindividual differences in psychologically mediated placebo response in migraine patients. Hum Brain Mapp 38:5250–5259. 10.1002/hbm.2372928731567
4. Planchuelo-Gómez Á García-Azorín D Guerrero ÁL Structural connectivity alterations in chronic and episodic migraine: a diffusion magnetic resonance imaging connectomics study Cephalalgia 2020 40 367 383 10.1177/0333102419885392 31674222
Planchuelo-Gómez Á, García-Azorín D, Guerrero ÁL et al (2020) Structural connectivity alterations in chronic and episodic migraine: a diffusion magnetic resonance imaging connectomics study. Cephalalgia 40:367–383. 10.1177/033310241988539231674222
5. Planchuelo-Gómez Á García‐Azorín D Guerrero ÁL Multimodal fusion analysis of structural connectivity and gray matter morphology in migraine Hum Brain Mapp 2021 42 908 921 10.1002/hbm.25267
Planchuelo-Gómez Á, García‐Azorín D, Guerrero ÁL et al (2021) Multimodal fusion analysis of structural connectivity and gray matter morphology in migraine. Hum Brain Mapp 42:908–921. 10.1002/hbm.25267
6. Chong CD Schwedt TJ Migraine affects white-matter tract integrity: a diffusion-tensor imaging study Cephalalgia 2015 35 1162 1171 10.1177/0333102415573513 25712998
Chong CD, Schwedt TJ (2015) Migraine affects white-matter tract integrity: a diffusion-tensor imaging study. Cephalalgia 35:1162–1171. 10.1177/033310241557351325712998
7. Michels L Koirala N Groppa S Structural brain network characteristics in patients with episodic and chronic migraine J Headache Pain 2021 22 8 10.1186/s10194-021-01216-8 33657996
Michels L, Koirala N, Groppa S et al (2021) Structural brain network characteristics in patients with episodic and chronic migraine. J Headache Pain 22:8. 10.1186/s10194-021-01216-833657996
8. Silvestro M Tessitore A Caiazzo G Disconnectome of the migraine brain: a connectopathy model J Headache Pain 2021 22 102 10.1186/s10194-021-01315-6 34454429
Silvestro M, Tessitore A, Caiazzo G et al (2021) Disconnectome of the migraine brain: a connectopathy model. J Headache Pain 22:102. 10.1186/s10194-021-01315-634454429
9. Dai L Zheng Q Chen X Altered brain structural topological properties and its correlations with clinical characteristics in episodic migraine without aura Neuroradiology 2021 63 2099 2109 10.1007/s00234-021-02716-9 34212221
Dai L, Zheng Q, Chen X et al (2021) Altered brain structural topological properties and its correlations with clinical characteristics in episodic migraine without aura. Neuroradiology 63:2099–2109. 10.1007/s00234-021-02716-934212221
10. Li K Liu L Yin Q Abnormal rich club organization and impaired correlation between structural and functional connectivity in migraine sufferers Brain Imaging Behav 2017 11 526 540 10.1007/s11682-016-9533-6 26922054
Li K, Liu L, Yin Q et al (2017) Abnormal rich club organization and impaired correlation between structural and functional connectivity in migraine sufferers. Brain Imaging Behav 11:526–540. 10.1007/s11682-016-9533-626922054
11. De Tommaso M Vecchio E Quitadamo SG Pain-related brain connectivity changes in migraine: a narrative review and proof of Concept about possible novel treatments interference Brain Sci 2021 11 234 10.3390/brainsci11020234 33668449
De Tommaso M, Vecchio E, Quitadamo SG et al (2021) Pain-related brain connectivity changes in migraine: a narrative review and proof of Concept about possible novel treatments interference. Brain Sci 11:234. 10.3390/brainsci1102023433668449
12. Cobos KL, Long X, Lebel C et al (2023) Increased hippocampal efficiency is associated with greater headache frequency in adolescents with chronic headache. 10.1093/texcom/tgad013. Cerebral Cortex Communications tgad013
13. Ke J Yu Y Zhang X Functional alterations in the posterior insula and cerebellum in Migraine without Aura: a resting-state MRI study Front Behav Neurosci 2020 14 567588 10.3389/fnbeh.2020.567588 33132860
Ke J, Yu Y, Zhang X et al (2020) Functional alterations in the posterior insula and cerebellum in Migraine without Aura: a resting-state MRI study. Front Behav Neurosci 14:567588. 10.3389/fnbeh.2020.56758833132860
14. Wang M Tutt JO Dorricott NO Involvement of the cerebellum in migraine Front Syst Neurosci 2022 16 984406 10.3389/fnsys.2022.984406 36313527
Wang M, Tutt JO, Dorricott NO et al (2022) Involvement of the cerebellum in migraine. Front Syst Neurosci 16:984406. 10.3389/fnsys.2022.98440636313527
15. Mehnert J May A Functional and structural alterations in the migraine cerebellum J Cereb Blood Flow Metab 2019 39 730 739 10.1177/0271678X17722109 28737061
Mehnert J, May A (2019) Functional and structural alterations in the migraine cerebellum. J Cereb Blood Flow Metab 39:730–739. 10.1177/0271678X1772210928737061
16. Qin Z He X-W Zhang J Structural changes of cerebellum and brainstem in migraine without aura J Headache Pain 2019 20 93 10.1186/s10194-019-1045-5 31477012
Qin Z, He X-W, Zhang J et al (2019) Structural changes of cerebellum and brainstem in migraine without aura. J Headache Pain 20:93. 10.1186/s10194-019-1045-531477012
17. Koppen H Boele H-J Palm-Meinders IH Cerebellar function and ischemic brain lesions in migraine patients from the general population Cephalalgia 2017 37 177 190 10.1177/0333102416643527 27059879
Koppen H, Boele H-J, Palm-Meinders IH et al (2017) Cerebellar function and ischemic brain lesions in migraine patients from the general population. Cephalalgia 37:177–190. 10.1177/033310241664352727059879
18. Rahimi R Dolatshahi M Abbasi-Feijani F Microstructural white matter alterations associated with migraine headaches: a systematic review of diffusion tensor imaging studies Brain Imaging Behav 2022 16 2375 2401 10.1007/s11682-022-00690-1 35710680
Rahimi R, Dolatshahi M, Abbasi-Feijani F et al (2022) Microstructural white matter alterations associated with migraine headaches: a systematic review of diffusion tensor imaging studies. Brain Imaging Behav 16:2375–2401. 10.1007/s11682-022-00690-135710680
19. Lucas C Migraine with aura Rev Neurol (Paris) 2021 177 779 784 10.1016/j.neurol.2021.07.010 34384631
Lucas C (2021) Migraine with aura. Rev Neurol (Paris) 177:779–784. 10.1016/j.neurol.2021.07.01034384631
20. Katsarava Z Buse DC Manack AN Lipton RB Defining the differences between episodic migraine and chronic migraine Curr Pain Headache Rep 2012 16 86 92 10.1007/s11916-011-0233-z 22083262
Katsarava Z, Buse DC, Manack AN, Lipton RB (2012) Defining the differences between episodic migraine and chronic migraine. Curr Pain Headache Rep 16:86–92. 10.1007/s11916-011-0233-z22083262
21. (2018) Headache Classification Committee of the International Headache Society (IHS) The International Classification of Headache Disorders, 3rd edition. Cephalalgia 38:1–211. 10.1177/0333102417738202
22. Woolrich MW Jbabdi S Patenaude B Bayesian analysis of neuroimaging data in FSL NeuroImage 2009 45 S173 S186 10.1016/j.neuroimage.2008.10.055 19059349
Woolrich MW, Jbabdi S, Patenaude B et al (2009) Bayesian analysis of neuroimaging data in FSL. NeuroImage 45:S173–S186. 10.1016/j.neuroimage.2008.10.05519059349
23. Tournier J-D Smith R Raffelt D MRtrix3: a fast, flexible and open software framework for medical image processing and visualisation NeuroImage 2019 202 116137 10.1016/j.neuroimage.2019.116137 31473352
Tournier J-D, Smith R, Raffelt D et al (2019) MRtrix3: a fast, flexible and open software framework for medical image processing and visualisation. NeuroImage 202:116137. 10.1016/j.neuroimage.2019.11613731473352
24. Ades-Aron B Veraart J Kochunov P Evaluation of the accuracy and precision of the diffusion parameter EStImation with Gibbs and NoisE removal pipeline NeuroImage 2018 183 532 543 10.1016/j.neuroimage.2018.07.066 30077743
Ades-Aron B, Veraart J, Kochunov P et al (2018) Evaluation of the accuracy and precision of the diffusion parameter EStImation with Gibbs and NoisE removal pipeline. NeuroImage 183:532–543. 10.1016/j.neuroimage.2018.07.06630077743
25. Koay CG Basser PJ Analytically exact correction scheme for signal extraction from noisy magnitude MR signals J Magn Reson 2006 179 317 322 10.1016/j.jmr.2006.01.016 16488635
Koay CG, Basser PJ (2006) Analytically exact correction scheme for signal extraction from noisy magnitude MR signals. J Magn Reson 179:317–322. 10.1016/j.jmr.2006.01.01616488635
26. Jeurissen B Tournier J-D Dhollander T Multi-tissue constrained spherical deconvolution for improved analysis of multi-shell diffusion MRI data NeuroImage 2014 103 411 426 10.1016/j.neuroimage.2014.07.061 25109526
Jeurissen B, Tournier J-D, Dhollander T et al (2014) Multi-tissue constrained spherical deconvolution for improved analysis of multi-shell diffusion MRI data. NeuroImage 103:411–426. 10.1016/j.neuroimage.2014.07.06125109526
27. Smith RE Tournier J-D Calamante F Connelly A Anatomically-constrained tractography: improved diffusion MRI streamlines tractography through effective use of anatomical information NeuroImage 2012 62 1924 1938 10.1016/j.neuroimage.2012.06.005 22705374
Smith RE, Tournier J-D, Calamante F, Connelly A (2012) Anatomically-constrained tractography: improved diffusion MRI streamlines tractography through effective use of anatomical information. NeuroImage 62:1924–1938. 10.1016/j.neuroimage.2012.06.00522705374
28. Tournier J-D, Calamante F, Connelly A (2010) Improved probabilistic streamlines tractography by 2nd order integration over fibre orientation distributions
29. Smith RE Tournier J-D Calamante F Connelly A SIFT2: enabling dense quantitative assessment of brain white matter connectivity using streamlines tractography NeuroImage 2015 119 338 351 10.1016/j.neuroimage.2015.06.092 26163802
Smith RE, Tournier J-D, Calamante F, Connelly A (2015) SIFT2: enabling dense quantitative assessment of brain white matter connectivity using streamlines tractography. NeuroImage 119:338–351. 10.1016/j.neuroimage.2015.06.09226163802
30. Schaefer A Kong R Gordon EM Local-global parcellation of the Human Cerebral Cortex from intrinsic functional connectivity MRI Cereb Cortex 2018 28 3095 3114 10.1093/cercor/bhx179 28981612
Schaefer A, Kong R, Gordon EM et al (2018) Local-global parcellation of the Human Cerebral Cortex from intrinsic functional connectivity MRI. Cereb Cortex 28:3095–3114. 10.1093/cercor/bhx17928981612
31. Tzourio-Mazoyer N Landeau B Papathanassiou D Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain NeuroImage 2002 15 273 289 10.1006/nimg.2001.0978 11771995
Tzourio-Mazoyer N, Landeau B, Papathanassiou D et al (2002) Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. NeuroImage 15:273–289. 10.1006/nimg.2001.097811771995
32. Zalesky A Fornito A Bullmore ET Network-based statistic: identifying differences in brain networks NeuroImage 2010 53 1197 1207 10.1016/j.neuroimage.2010.06.041 20600983
Zalesky A, Fornito A, Bullmore ET (2010) Network-based statistic: identifying differences in brain networks. NeuroImage 53:1197–1207. 10.1016/j.neuroimage.2010.06.04120600983
33. Yeo BTT Krienen FM Chee MWL Buckner RL Estimates of segregation and overlap of functional connectivity networks in the human cerebral cortex NeuroImage 2014 88 212 227 10.1016/j.neuroimage.2013.10.046 24185018
Yeo BTT, Krienen FM, Chee MWL, Buckner RL (2014) Estimates of segregation and overlap of functional connectivity networks in the human cerebral cortex. NeuroImage 88:212–227. 10.1016/j.neuroimage.2013.10.04624185018
34. Rubinov M Sporns O Complex network measures of brain connectivity: uses and interpretations NeuroImage 2010 52 1059 1069 10.1016/j.neuroimage.2009.10.003 19819337
Rubinov M, Sporns O (2010) Complex network measures of brain connectivity: uses and interpretations. NeuroImage 52:1059–1069. 10.1016/j.neuroimage.2009.10.00319819337
35. Hua K Zhang J Wakana S Tract probability maps in stereotaxic spaces: analyses of white matter anatomy and tract-specific quantification NeuroImage 2008 39 336 347 10.1016/j.neuroimage.2007.07.053 17931890
Hua K, Zhang J, Wakana S et al (2008) Tract probability maps in stereotaxic spaces: analyses of white matter anatomy and tract-specific quantification. NeuroImage 39:336–347. 10.1016/j.neuroimage.2007.07.05317931890
36. Moulton EA Elman I Pendse G Aversion-related circuitry in the Cerebellum: responses to noxious heat and unpleasant images J Neurosci 2011 31 3795 3804 10.1523/JNEUROSCI.6709-10.2011 21389234
Moulton EA, Elman I, Pendse G et al (2011) Aversion-related circuitry in the Cerebellum: responses to noxious heat and unpleasant images. J Neurosci 31:3795–3804. 10.1523/JNEUROSCI.6709-10.201121389234
37. Stoodley CJ Schmahmann JD Evidence for topographic organization in the cerebellum of motor control versus cognitive and affective processing Cortex 2010 46 831 844 10.1016/j.cortex.2009.11.008 20152963
Stoodley CJ, Schmahmann JD (2010) Evidence for topographic organization in the cerebellum of motor control versus cognitive and affective processing. Cortex 46:831–844. 10.1016/j.cortex.2009.11.00820152963
38. Gu L Wang Y Shu H Association between migraine and cognitive impairment J Headache Pain 2022 23 88 10.1186/s10194-022-01462-4 35883043
Gu L, Wang Y, Shu H (2022) Association between migraine and cognitive impairment. J Headache Pain 23:88. 10.1186/s10194-022-01462-435883043
39. Giffin NJ Ruggiero L Lipton RB Premonitory symptoms in migraine: an electronic diary study Neurology 2003 60 935 940 10.1212/01.WNL.0000052998.58526.A9 12654956
Giffin NJ, Ruggiero L, Lipton RB et al (2003) Premonitory symptoms in migraine: an electronic diary study. Neurology 60:935–940. 10.1212/01.WNL.0000052998.58526.A912654956
40. Gil-Gouveia R Martins IP Cognition and cognitive impairment in Migraine Curr Pain Headache Rep 2019 23 84 10.1007/s11916-019-0824-7 31511992
Gil-Gouveia R, Martins IP (2019) Cognition and cognitive impairment in Migraine. Curr Pain Headache Rep 23:84. 10.1007/s11916-019-0824-731511992
41. C. S, W. W (2001) Nociceptive visceral stimulation modulates the activity of cerebellar Purkinje cells. Exp Brain Res 140:122–126. 10.1007/s002210100824
42. Younis S Hougaard A Noseda R Ashina M Current understanding of thalamic structure and function in migraine Cephalalgia 2019 39 1675 1682 10.1177/0333102418791595 30079744
Younis S, Hougaard A, Noseda R, Ashina M (2019) Current understanding of thalamic structure and function in migraine. Cephalalgia 39:1675–1682. 10.1177/033310241879159530079744
43. Middleton FA Strick PL Cerebellar projections to the Prefrontal cortex of the Primate J Neurosci 2001 21 700 712 10.1523/JNEUROSCI.21-02-00700.2001 11160449
Middleton FA, Strick PL (2001) Cerebellar projections to the Prefrontal cortex of the Primate. J Neurosci 21:700–712. 10.1523/JNEUROSCI.21-02-00700.200111160449
44. Morara S Sternini C Provini L Rosina A Developmentally regulated expression of α- and β‐calcitonin gene‐related peptide mRNA and calcitonin gene‐related peptide immunoreactivity in the rat inferior olive J Comp Neurol 1995 354 27 38 10.1002/cne.903540104 7615873
Morara S, Sternini C, Provini L, Rosina A (1995) Developmentally regulated expression of α- and β‐calcitonin gene‐related peptide mRNA and calcitonin gene‐related peptide immunoreactivity in the rat inferior olive. J Comp Neurol 354:27–38. 10.1002/cne.9035401047615873
45. Morara S Rosina A Provini L CGRP as a marker of the climbing fibers during the development of the Cerebellum in the rat Ann N Y Acad Sci 1992 657 461 463 10.1111/j.1749-6632.1992.tb22800.x 1637100
Morara S, Rosina A, Provini L (1992) CGRP as a marker of the climbing fibers during the development of the Cerebellum in the rat. Ann N Y Acad Sci 657:461–463. 10.1111/j.1749-6632.1992.tb22800.x1637100
46. Goadsby PJ, Silberstein SD, Yeung PP et al (2020) Long-term safety, tolerability, and efficacy of fremanezumab in migraine: a randomized study. Neurology 95. 10.1212/WNL.0000000000010600
47. Goadsby PJ, Reuter U, Lanteri-Minet M et al (2021) Long-term efficacy and safety of Erenumab: results from 64 weeks of the LIBERTY Study. Neurology 96. 10.1212/WNL.0000000000012029
48. Schmitz N Admiraal-Behloul F Arkink EB Attack frequency and disease duration as indicators for brain damage in Migraine Headache 2008 48 1044 1055 10.1111/j.1526-4610.2008.01133.x 18479421
Schmitz N, Admiraal-Behloul F, Arkink EB et al (2008) Attack frequency and disease duration as indicators for brain damage in Migraine. Headache 48:1044–1055. 10.1111/j.1526-4610.2008.01133.x18479421
49. Szabó N Faragó P Király A Evidence for plastic processes in migraine with aura: a Diffusion Weighted MRI Study Front Neuroanat 2018 11 138 10.3389/fnana.2017.00138 29387002
Szabó N, Faragó P, Király A et al (2018) Evidence for plastic processes in migraine with aura: a Diffusion Weighted MRI Study. Front Neuroanat 11:138. 10.3389/fnana.2017.0013829387002
50. Veréb D Szabó N Tuka B Temporal instability of salience network activity in migraine with aura Pain 2020 161 856 10.1097/j.pain.0000000000001770 31815918
Veréb D, Szabó N, Tuka B et al (2020) Temporal instability of salience network activity in migraine with aura. Pain 161:856. 10.1097/j.pain.000000000000177031815918
51. Wang T Chen N Zhan W Altered effective connectivity of posterior thalamus in migraine with cutaneous allodynia: a resting-state fMRI study with granger causality analysis J Headache Pain 2016 17 17 10.1186/s10194-016-0610-4
Wang T, Chen N, Zhan W et al (2016) Altered effective connectivity of posterior thalamus in migraine with cutaneous allodynia: a resting-state fMRI study with granger causality analysis. J Headache Pain 17:17. 10.1186/s10194-016-0610-4
52. Schmahmann JD The cerebellum and cognition Neurosci Lett 2019 688 62 75 10.1016/j.neulet.2018.07.005 29997061
Schmahmann JD (2019) The cerebellum and cognition. Neurosci Lett 688:62–75. 10.1016/j.neulet.2018.07.00529997061
53. Waxman SG Clinical neuroanatomy, 28e 2017 New York, NY McGraw-Hill Education
Waxman SG (2017) Clinical neuroanatomy, 28e. McGraw-Hill Education, New York, NY
54. Marciszewski KK, Meylakh N, Di Pietro F et al (2019) Fluctuating Regional Brainstem Diffusion Imaging measures of microstructure across the Migraine cycle. 10.1523/ENEURO.0005-19.2019. eNeuro 6:ENEURO.0005-19.2019
55. Yeh C Jones DK Liang X Mapping structural connectivity using Diffusion MRI: challenges and opportunities Magn Reson Imaging 2021 53 1666 1682 10.1002/jmri.27188
Yeh C, Jones DK, Liang X et al (2021) Mapping structural connectivity using Diffusion MRI: challenges and opportunities. Magn Reson Imaging 53:1666–1682. 10.1002/jmri.27188
56. Calamante F The Seven Deadly sins of measuring Brain Structural Connectivity using Diffusion MRI streamlines Fibre-Tracking Diagnostics 2019 9 115 10.3390/diagnostics9030115 31500098
Calamante F (2019) The Seven Deadly sins of measuring Brain Structural Connectivity using Diffusion MRI streamlines Fibre-Tracking. Diagnostics 9:115. 10.3390/diagnostics903011531500098
57. Bassett DS Bullmore ET Small-world brain networks revisited Neuroscientist 2017 23 499 516 10.1177/1073858416667720 27655008
Bassett DS, Bullmore ET (2017) Small-world brain networks revisited. Neuroscientist 23:499–516. 10.1177/107385841666772027655008
