Skip to main content

Involvement of the cerebellum in structural connectivity enhancement in episodic migraine

Abstract

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.

Peer Review reports

Introduction

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,3,4,5,6,7,8,9,10,11] and functional [2, 10,11,12,13,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,9,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,5,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
figure 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:

$$\:{a}_{ij}=\frac{2}{{V}_{i}+{V}_{j}}{\:n}_{ij}$$
(1)

where \(\:{a}_{ij}\) is the entry of the connectivity matrix concerning ROI i and ROI j, \(\:{n}_{ij}\) is the number of streamlines between the two ROIs given by the tractogram, and \(\:{V}_{i}\) and \(\:{V}_{j}\) 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

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
figure 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
figure 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
figure 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
figure 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
figure 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,39,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,9,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.

Data availability

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

References

  1. 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. https://doi.org/10.1186/s10194-022-01402-2

    Article  Google Scholar 

  2. Ashina S, Bentivegna E, Martelletti P, Eikermann-Haerter K (2021) Structural and functional brain changes in Migraine. Pain Ther 10:211–223. https://doi.org/10.1007/s40122-021-00240-5

    Article  PubMed  PubMed Central  Google Scholar 

  3. 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. https://doi.org/10.1002/hbm.23729

    Article  PubMed  PubMed Central  Google Scholar 

  4. 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. https://doi.org/10.1177/0333102419885392

    Article  PubMed  Google Scholar 

  5. 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. https://doi.org/10.1002/hbm.25267

    Article  Google Scholar 

  6. Chong CD, Schwedt TJ (2015) Migraine affects white-matter tract integrity: a diffusion-tensor imaging study. Cephalalgia 35:1162–1171. https://doi.org/10.1177/0333102415573513

    Article  PubMed  Google Scholar 

  7. 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. https://doi.org/10.1186/s10194-021-01216-8

    Article  PubMed  PubMed Central  Google Scholar 

  8. Silvestro M, Tessitore A, Caiazzo G et al (2021) Disconnectome of the migraine brain: a connectopathy model. J Headache Pain 22:102. https://doi.org/10.1186/s10194-021-01315-6

    Article  PubMed  PubMed Central  Google Scholar 

  9. 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. https://doi.org/10.1007/s00234-021-02716-9

    Article  PubMed  Google Scholar 

  10. 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. https://doi.org/10.1007/s11682-016-9533-6

    Article  PubMed  Google Scholar 

  11. 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. https://doi.org/10.3390/brainsci11020234

    Article  PubMed  PubMed Central  Google Scholar 

  12. Cobos KL, Long X, Lebel C et al (2023) Increased hippocampal efficiency is associated with greater headache frequency in adolescents with chronic headache. https://doi.org/10.1093/texcom/tgad013. Cerebral Cortex Communications tgad013

  13. 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. https://doi.org/10.3389/fnbeh.2020.567588

    Article  PubMed  PubMed Central  Google Scholar 

  14. Wang M, Tutt JO, Dorricott NO et al (2022) Involvement of the cerebellum in migraine. Front Syst Neurosci 16:984406. https://doi.org/10.3389/fnsys.2022.984406

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Mehnert J, May A (2019) Functional and structural alterations in the migraine cerebellum. J Cereb Blood Flow Metab 39:730–739. https://doi.org/10.1177/0271678X17722109

    Article  PubMed  Google Scholar 

  16. 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. https://doi.org/10.1186/s10194-019-1045-5

    Article  PubMed  PubMed Central  Google Scholar 

  17. 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. https://doi.org/10.1177/0333102416643527

    Article  PubMed  Google Scholar 

  18. 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. https://doi.org/10.1007/s11682-022-00690-1

    Article  PubMed  PubMed Central  Google Scholar 

  19. Lucas C (2021) Migraine with aura. Rev Neurol (Paris) 177:779–784. https://doi.org/10.1016/j.neurol.2021.07.010

    Article  PubMed  CAS  Google Scholar 

  20. 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. https://doi.org/10.1007/s11916-011-0233-z

    Article  PubMed  Google Scholar 

  21. (2018) Headache Classification Committee of the International Headache Society (IHS) The International Classification of Headache Disorders, 3rd edition. Cephalalgia 38:1–211. https://doi.org/10.1177/0333102417738202

  22. Woolrich MW, Jbabdi S, Patenaude B et al (2009) Bayesian analysis of neuroimaging data in FSL. NeuroImage 45:S173–S186. https://doi.org/10.1016/j.neuroimage.2008.10.055

    Article  PubMed  Google Scholar 

  23. 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. https://doi.org/10.1016/j.neuroimage.2019.116137

    Article  PubMed  Google Scholar 

  24. 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. https://doi.org/10.1016/j.neuroimage.2018.07.066

    Article  PubMed  Google Scholar 

  25. Koay CG, Basser PJ (2006) Analytically exact correction scheme for signal extraction from noisy magnitude MR signals. J Magn Reson 179:317–322. https://doi.org/10.1016/j.jmr.2006.01.016

    Article  PubMed  CAS  Google Scholar 

  26. 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. https://doi.org/10.1016/j.neuroimage.2014.07.061

    Article  PubMed  Google Scholar 

  27. 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. https://doi.org/10.1016/j.neuroimage.2012.06.005

    Article  PubMed  Google Scholar 

  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 (2015) SIFT2: enabling dense quantitative assessment of brain white matter connectivity using streamlines tractography. NeuroImage 119:338–351. https://doi.org/10.1016/j.neuroimage.2015.06.092

    Article  PubMed  Google Scholar 

  30. 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. https://doi.org/10.1093/cercor/bhx179

    Article  PubMed  Google Scholar 

  31. 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. https://doi.org/10.1006/nimg.2001.0978

    Article  PubMed  CAS  Google Scholar 

  32. Zalesky A, Fornito A, Bullmore ET (2010) Network-based statistic: identifying differences in brain networks. NeuroImage 53:1197–1207. https://doi.org/10.1016/j.neuroimage.2010.06.041

    Article  PubMed  Google Scholar 

  33. 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. https://doi.org/10.1016/j.neuroimage.2013.10.046

    Article  PubMed  Google Scholar 

  34. Rubinov M, Sporns O (2010) Complex network measures of brain connectivity: uses and interpretations. NeuroImage 52:1059–1069. https://doi.org/10.1016/j.neuroimage.2009.10.003

    Article  PubMed  Google Scholar 

  35. 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. https://doi.org/10.1016/j.neuroimage.2007.07.053

    Article  PubMed  Google Scholar 

  36. 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. https://doi.org/10.1523/JNEUROSCI.6709-10.2011

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Stoodley CJ, Schmahmann JD (2010) Evidence for topographic organization in the cerebellum of motor control versus cognitive and affective processing. Cortex 46:831–844. https://doi.org/10.1016/j.cortex.2009.11.008

    Article  PubMed  PubMed Central  Google Scholar 

  38. Gu L, Wang Y, Shu H (2022) Association between migraine and cognitive impairment. J Headache Pain 23:88. https://doi.org/10.1186/s10194-022-01462-4

    Article  PubMed  PubMed Central  Google Scholar 

  39. Giffin NJ, Ruggiero L, Lipton RB et al (2003) Premonitory symptoms in migraine: an electronic diary study. Neurology 60:935–940. https://doi.org/10.1212/01.WNL.0000052998.58526.A9

    Article  PubMed  CAS  Google Scholar 

  40. Gil-Gouveia R, Martins IP (2019) Cognition and cognitive impairment in Migraine. Curr Pain Headache Rep 23:84. https://doi.org/10.1007/s11916-019-0824-7

    Article  PubMed  Google Scholar 

  41. C. S, W. W (2001) Nociceptive visceral stimulation modulates the activity of cerebellar Purkinje cells. Exp Brain Res 140:122–126. https://doi.org/10.1007/s002210100824

  42. Younis S, Hougaard A, Noseda R, Ashina M (2019) Current understanding of thalamic structure and function in migraine. Cephalalgia 39:1675–1682. https://doi.org/10.1177/0333102418791595

    Article  PubMed  Google Scholar 

  43. Middleton FA, Strick PL (2001) Cerebellar projections to the Prefrontal cortex of the Primate. J Neurosci 21:700–712. https://doi.org/10.1523/JNEUROSCI.21-02-00700.2001

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. 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. https://doi.org/10.1002/cne.903540104

    Article  PubMed  CAS  Google Scholar 

  45. 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. https://doi.org/10.1111/j.1749-6632.1992.tb22800.x

    Article  PubMed  CAS  Google Scholar 

  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. https://doi.org/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. https://doi.org/10.1212/WNL.0000000000012029

  48. 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. https://doi.org/10.1111/j.1526-4610.2008.01133.x

    Article  PubMed  Google Scholar 

  49. 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. https://doi.org/10.3389/fnana.2017.00138

    Article  PubMed  PubMed Central  Google Scholar 

  50. Veréb D, Szabó N, Tuka B et al (2020) Temporal instability of salience network activity in migraine with aura. Pain 161:856. https://doi.org/10.1097/j.pain.0000000000001770

    Article  PubMed  Google Scholar 

  51. 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. https://doi.org/10.1186/s10194-016-0610-4

    Article  PubMed Central  Google Scholar 

  52. Schmahmann JD (2019) The cerebellum and cognition. Neurosci Lett 688:62–75. https://doi.org/10.1016/j.neulet.2018.07.005

    Article  PubMed  CAS  Google Scholar 

  53. Waxman SG (2017) Clinical neuroanatomy, 28e. McGraw-Hill Education, New York, NY

    Google Scholar 

  54. Marciszewski KK, Meylakh N, Di Pietro F et al (2019) Fluctuating Regional Brainstem Diffusion Imaging measures of microstructure across the Migraine cycle. https://doi.org/10.1523/ENEURO.0005-19.2019. eNeuro 6:ENEURO.0005-19.2019

  55. Yeh C, Jones DK, Liang X et al (2021) Mapping structural connectivity using Diffusion MRI: challenges and opportunities. Magn Reson Imaging 53:1666–1682. https://doi.org/10.1002/jmri.27188

    Article  Google Scholar 

  56. Calamante F (2019) The Seven Deadly sins of measuring Brain Structural Connectivity using Diffusion MRI streamlines Fibre-Tracking. Diagnostics 9:115. https://doi.org/10.3390/diagnostics9030115

    Article  PubMed  PubMed Central  Google Scholar 

  57. Bassett DS, Bullmore ET (2017) Small-world brain networks revisited. Neuroscientist 23:499–516. https://doi.org/10.1177/1073858416667720

    Article  PubMed  Google Scholar 

Download references

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.

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: https://doi.org/10.54499/LA/P/0083/2020, https://doi.org/10.54499/UIDP/50009/2020, and https://doi.org/10.54499/UIDB/50009/2020).

Author information

Authors and Affiliations

Authors

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.

Corresponding author

Correspondence to Ana Matoso.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

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/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Matoso, A., Fouto, A.R., Esteves, I. et al. Involvement of the cerebellum in structural connectivity enhancement in episodic migraine. J Headache Pain 25, 154 (2024). https://doi.org/10.1186/s10194-024-01854-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s10194-024-01854-8

Keywords