Skip to main content

Shifts in structural connectome organization in the limbic and sensory systems of patients with episodic migraine

Abstract

Migraine is a complex neurological condition characterized by recurrent headaches, which is often accompanied by various neurological symptoms. Magnetic resonance imaging (MRI) is a powerful tool for investigating whole-brain connectivity patterns; however, systematic assessment of structural connectome organization has rarely been performed. In the present study, we aimed to examine the changes in structural connectivity in patients with episodic migraines using diffusion MRI. First, we computed structural connectivity using diffusion MRI tractography, after which we applied dimensionality reduction techniques to the structural connectivity and generated three low-dimensional eigenvectors. We subsequently calculated the manifold eccentricity, defined as the Euclidean distance between each data point and the center of the data in the manifold space. We then compared the manifold eccentricity between patients with migraines and healthy controls, revealing significant between-group differences in the orbitofrontal cortex, temporal pole, and sensory/motor regions. Between-group differences in subcortico-cortical connectivity further revealed significant changes in the amygdala, accumbens, and caudate nuclei. Finally, supervised machine learning effectively classified patients with migraines and healthy controls using cortical and subcortical structural connectivity features, highlighting the importance of the orbitofrontal and sensory cortices, in addition to the caudate, in distinguishing between the groups. Our findings confirmed that episodic migraine is related to the structural connectome changes in the limbic and sensory systems, suggesting its potential utility as a diagnostic marker for migraine.

Peer Review reports

Introduction

Migraine is a neurological disorder characterized by recurrent headaches accompanied by nausea, vomiting, and an increased sensitivity to sound or light [1]. Magnetic resonance imaging (MRI) is widely used as a technique to uncover the neurological underpinnings of migraines [2]. In particular, structural connectivity, defined using diffusion tractography, has previously been adopted to assess changes in inter-regional relationships in tract strengths between different brain regions, beyond the simple investigation of regional abnormalities in patients with migraine [3, 4]. By calculating graph-theoretical parameters such as nodal centrality and efficiency measures, researchers have observed changes in the tract strengths in pain-processing-related brain regions, particularly in the sensory and limbic networks, which were found to be associated with prolonged disease duration and higher pain severity [3, 4].

Several challenges have arisen in these studies. Indeed, previous studies of structural connectivity have reported varying results [3,4,5,6]. For example, one study found changes in the structural connectome hubs in the limbic and emotional networks that were not related to disease severity [3]. In contrast, another study observed network differences in the occipital, temporal, and parietal lobes associated with disease severity and duration [4]. This inconsistency may be due to methodological issues in diffusion tractography, such as the choice of features for constructing structural connectivity (e.g., streamline count, density, or cross-section) and tractography methods (e.g., deterministic or probabilistic). As such, methods that can extract major information from structural connectivity data without considering the information that may vary across methodological choices must be considered.

Manifold learning, also called gradient analysis, is one technique that can be used to explore connectome organization [7,8,9,10]. The essence of this analytical method lies in dimensionality reduction, which reduces high-dimensional connectome data to multiple low-dimensional eigenvectors (i.e., gradients) [8, 10], allowing for an effective explanation of the entire dataset using a small set of gradients, while retaining sufficient information [11,12,13]. The reliability of the gradient technique has previously been demonstrated for several neurological disorders [11,12,13,14]. For example, our previous study successfully revealed alterations in functional connectome organization in patients with migraines using functional MRI-based gradient analysis [14]. The only limitation of gradient analysis lies in the interpretation of the results, as only multivariate analysis can be used to assess between-group differences across multiple gradients. Thus, the directionality of association cannot be obtained from univariate analysis.

In the present study, we opted to perform gradient analysis to investigate whole-brain structural connectome changes in patients with episodic migraine. Specifically, we used “manifold eccentricity,” which consolidates multiple gradients into a single vector [15]. This enabled us to assess the directionality of the statistics by quantifying the expansion and contraction of the data points in the manifold space. We further hypothesized that patients with migraine would show changes in manifold eccentricity compared to healthy controls.

Methods

Study participants

Patients with episodic migraines were recruited from an academic headache clinic between October 2020 and November 2021. The inclusion criteria for patients were: (1) age 18–50 years; (2) not taking preventive medications; and (3) premenopausal status in female patients. The exclusion criteria included: (1) chronic migraine, medication-overuse headaches, chronic pain disorders other than migraine, and psychiatric disorders such as bipolar affective disorder or schizophrenia; (2) contraindications for 3T MRI, including use of a tissue expander, pacemaker, non-detachable metal objects, orthodontic devices, or electrical leads or implants in the body; (3) pregnancy; (4) claustrophobia requiring sedation during MR scanning; (5) inability to report their headache or complete the headache diary due to cognitive decline; and (6) disagreement with the study procedures. Migraine was diagnosed according to the ICHD-3 criteria by a single investigator (MJL), a neurologist specializing in headaches. This study was approved by the Samsung Medical Center Institutional Review Board, and all participants provided written informed consent to participate. This study was part of an ongoing longitudinal project registered at ClinicalTrials.gov (NCT03487978).

MRI data acquisition

T1-weighted structural and diffusion MRI data were obtained using a Philips Ingenia 3T scanner (Philips, Amsterdam, Netherlands). T1-weighted data were acquired using a turbo field echo sequence in the sagittal plane (repetition time [TR] = 8.1 ms; echo time [TE] = 3.7 ms; field of view [FOV] = 256 × 256 mm2; voxel size = 1 mm isotropic; and number of slices = 180). Diffusion MRI data were acquired using a spin-echo echo-planar imaging sequence in the axial plane (TR = 7,062 ms; TE = 91 ms; FOV = 220 × 220 mm2; voxel size = 1.719 × 1.719 × 3 mm3; number of slices = 47; b-value = 1,000 s/mm2; number of diffusion directions = 64; and number of b0 images = 1). T1-weighted MRI scanning was performed for 3 min 20 s, and diffusion MRI for 7 min 50 s.

MRI data preprocessing

The T1-weighted MRI data were preprocessed using the fusion of neuroimaging preprocessing (FuNP) surface-based pipeline [16], which included gradient nonlinearity correction, non-brain tissue removal, and intensity normalization. The cortical surfaces were generated using FreeSurfer v7.1.1 (Boston, MA, USA) [17] by following the boundaries between the white and pial surfaces [18,19,20]. The mid-thickness surface was generated by averaging the white and pial surfaces, and was used to generate an inflated surface. Diffusion MRI data were preprocessed using MRtrix3 v3.0.2 [21], including denoising, Gibbs ringing artifact removal, susceptibility distortion correction, head motion correction, and eddy current correction. Anatomically-constrained tractography was further performed using different tissue types of the cortical and subcortical gray matter, white matter, and cerebrospinal fluid, as defined using T1-weighted MRI [22]. T1-weighted and diffusion MRI data were co-registered using FSL v6.0 [23], and different tissue types were transformed into the native diffusion MRI space. After estimating the single-shell response functions [24], constrained spherical deconvolution and intensity normalization were subsequently performed [25]. A tractogram was further generated using a probabilistic approach with 40 million streamlines [26, 27]. Options were set with a maximum tract length of 250 and a fractional anisotropy cutoff of 0.06. Spherical deconvolution-informed filtering of tractograms (SIFT2) was further applied to reconstruct the whole-brain streamlines weighted by cross-section multipliers [28]. The structural connectivity matrix was constructed by mapping the reconstructed cross-section streamlines onto a sub-parcellation of the Desikan–Killiany atlas with 300 parcels [29]. The sub-parcellation of the Desikan–Killiany atlas was defined in MICAPIPE (https://github.com/MICA-MNI/micapipe/tree/master/parcellations) [30], which subdivided the FreeSurfer segmentation into several sub-regions.

Structural connectivity gradient generation

We generated low-dimensional representations of structural connectivity (i.e., gradients) by applying manifold learning, a nonlinear dimensionality reduction technique [10]. The application of manifold learning to high-dimensional data allows the generation of multiple low-dimensional eigenvectors to construct a newly defined low-dimensional space (i.e., a manifold space). First, we constructed a group representative matrix computed using distance-dependent thresholding to preserve long-range connections [31]. Using the BrainSpace toolbox (https://github.com/MICA-MNI/BrainSpace [10]), we then applied diffusion map embedding, a robust and computationally efficient nonlinear dimensionality reduction technique [32, 33], to the group-representative structural connectivity matrix. Individual gradients were estimated by applying diffusion map embedding to the individual structural connectivity matrix and were aligned to group representative gradients using Procrustes alignment [10, 34].

Manifold eccentricity and between-group differences

We subsequently generated three structural connectivity gradients that sufficiently explained total connectivity information and showed biologically interpretable axes [15, 35, 36]. Multiple gradients were further summarized into a single feature termed manifold eccentricity, defined as the Euclidean distance between each data point and the center of the data in the manifold space [15]. Differences in manifold eccentricity between patients with migraines and healthy controls were assessed using non-parametric permutation tests. The subject indices were randomly shuffled, and an analysis of covariance (ANCOVA) was performed by controlling for age and sex. This process was repeated 10,000 times. A null distribution of the between-group differences was constructed, and the p-value was calculated by dividing the number of absolute permuted statistical values larger than absolute value of the real statistic by the number of permutations. Multiple comparisons across brain regions were corrected using a false discovery rate (FDR) < 0.05 [37]. To quantify the between-group difference effects on manifold eccentricity according to brain networks, we further summarized the statistical values according to seven intrinsic functional communities, including the visual, somatomotor, dorsal attention, ventral attention, limbic, frontoparietal, and default mode cortices [38], and cortical hierarchical levels, including idiotypic, unimodal association, heteromodal association, and paralimbic areas [39].

Changes in subcortical structures

In addition to investigating the changes in cortical manifold eccentricity, we further assessed the between-group differences in subcortical structural connectivity by analyzing the degree of subcortico-cortical structural connectivity, defined as the sum of edge weights connecting each subcortical region to all cortical regions. The subcortical structures of the thalamus, caudate, putamen, pallidum, hippocampus, amygdala, and accumbens were defined from the T1-weighted data using FSL FIRST [40]. We further conducted 10,000 permutation tests, and multiple comparisons across the subcortical structures were corrected using an FDR threshold of < 0.05.

Classification between patients with migraine and healthy controls

To validate the above features, we adopted supervised machine learning to classify patients with migraines and healthy controls using the cortical manifold eccentricity and degree values of subcortico-cortical connectivity. Specifically, following the regression of age and sex from the features, we applied the least absolute shrinkage and selection operator (LASSO) method to select the imaging features [41], and entered the selected features into a linear regression model. We further performed the classification task using a five-fold cross-validation framework by dividing the data into training (4/5 partitions) and testing (1/5 partitions) datasets. The procedure was repeated 100 times using different training and testing datasets to avoid subject selection bias. Classification performance was assessed using precision, recall, and area under the receiver operating characteristic (ROC) curves (AUC), and the mean scores with standard deviation (SD) across 100 repetitions were reported.

Sensitivity analysis

A) Excluding patients with migraine with aura

To evaluate the impact of migraine with or without aura on between-group differences between patients with migraine and healthy controls, we further assessed the differences in the manifold eccentricity and degree values of subcortico-cortical connectivity after excluding patients with migraine who had an aura (n = 7).

B) Different parcellations

To assess the robustness of the findings across different spatial scales, we repeated the analyses using a sub-parcellation of the Desikan–Killiany atlas with 200 and 400 parcels [29]. We further performed the same analyses using the Schaefer atlas with 300 parcels to represent a functional parcellation scheme [42].

C) Different migraine phases

We obtained the participants’ headache status within ± 2 days of MRI scanning. Patients were considered ictal if they had headaches of any intensity on the day of scanning, interictal if they were headache-free on ± 2 days of scanning, and peri-ictal if they were headache-free on the day of scanning, but developed headaches within two days of scanning. To investigate whether the migraine phase affects the structural connectome organization, we conducted separate analyses for patients in the interictal, peri-ictal, and ictal phases.

D) Effects of depression and anxiety

To evaluate the impact of anxiety or depression on migraines, we assessed the differences in manifold eccentricity and degree values of subcortico-cortical structural connectivity between healthy controls and migraine patients without depression (n = 39), as well as between healthy controls and migraine patients with depression (n = 8). We further compared two between-group differences by calculating the linear correlation, and conducted the same analysis for patients with migraine without anxiety (n = 39) and with anxiety (n = 8).

E) Low- vs. high-frequency episodic migraine

A sensitivity analysis between low- and high-frequency episodic migraines was not conducted because of the small sample size of patients with high-frequency episodic migraines (n = 4). Instead, we repeated the analysis to assess the between-group differences using healthy controls and patients with migraines by considering only low-frequency episodic migraines (n = 43).

F) Manifold eccentricity calculation using multiple eigenvectors

We generated the manifold eccentricity using three eigenvectors that explained approximately 38% of the connectome information in the primary analysis. To assess the effects of the number of eigenvectors, we further performed independent analyses by calculating manifold eccentricity using multiple eigenvectors that explained approximately 50, 70, and 90% of the connectome information.

G) Between-group differences without controlling for age and sex

In the primary analysis, we controlled for age and sex, while assessing between-group differences between patients with migraines and healthy controls. In addition, we performed the analysis without controlling for age and sex.

Results

Study participants

We recruited 108 participants, including 65 patients with migraines and 43 healthy controls, and excluded 18 patients and two healthy controls. Specifically, we excluded seven patients who did not undergo T1-weighted or diffusion MRI scanning, as well as 11 patients and two healthy controls in whom cortical surface reconstruction failed due to severe head movements. Finally, 47 patients with migraine (mean ± SD age = 34.3 ± 8.3 years; sex = 74.5% female) and 41 healthy controls (age = 35.2 ± 7.7 years; sex = 75.6% female) were enrolled in this study. Among the patients with migraine, 26 were in the interictal phase, six in the peri-ictal phase, and 15 in the ictal phase. Seven patients had migraine with aura, while the remaining 40 did not. Patients showed a mean of 6.9 ± 4.84 days of headache per month and 4.2 ± 2.96 days of migraine per month. In particular, four patients had high-frequency episodic migraines (8 or more migraine days/month), eight had depression, and eight had anxiety. The detailed demographic and clinical information of the participants are presented in Table 1.

Table 1 Demographic and clinical information of study participants

Between-group differences in the manifold eccentricity

We used the first three gradients that accounted for approximately 38.23% of the structural connectivity matrix. Consistent with previous studies [15, 35, 43], the first gradient (G1) showed left-right, the second (G2) demonstrated anterior-posterior, and the third (G3) exhibited sensory-transmodal axes (Fig. 1A). After calculating the manifold eccentricity (Fig. 1B) [15], we assessed between-group differences, which revealed significant expansions in the orbitofrontal cortex and temporal pole and contractions in the sensorimotor regions in patients with migraines (Fig. 1C). When summarizing the between-group differences according to the seven intrinsic functional communities [38] and cortical hierarchical levels [39], the limbic network showed the strongest effects, followed by the heteromodal association areas (Fig. 1C).

Fig. 1
figure 1

Between-group differences in manifold eccentricity between patients with migraines and healthy controls. (A) The structural connectome was estimated using diffusion MRI tractography (top). The group representative matrix was computed using distance-dependent thresholding, and the group-level gradients (G1, G2, and G3) were generated using a diffusion map embedding algorithm. Individual gradients were aligned to the group gradients using Procrustes alignment. (B) The schema of manifold eccentricity is shown with a three-dimensional scatter plot (top). The average manifold eccentricity maps of healthy controls and patients with migraines are shown on the brain surfaces (bottom). (C) T-statistics of between-group differences in the manifold eccentricity are shown (left top). The effects were stratified according to seven intrinsic function communities and four cortical hierarchical levels using spider plots (left bottom). Expansion or contraction in the manifold eccentricity of each brain region is shown as dots using arrows to differentiate healthy controls (blue) from patients with migraines (red). Gray dots indicate the regions that did not show significant effects (right). Abbreviations: FDR, false discovery rate; MRI, magnetic resonance imaging; SC, structural connectivity

Between-group differences in subcortico-cortical connectivity

Using the degree values of the subcortico-cortical structural connectivity (Fig. 2A), we identified significant between-group differences in the caudate, amygdala, and accumbens, characterized by increased degree values in patients with migraines (Fig. 2B).

Fig. 2
figure 2

Between-group differences in the subcortico-cortical structural connectivity. (A) The degree values of each subcortical region were calculated with a row-wise sum of the subcortico-cortical structural connectivity. (B) The t-statistics of between-group differences in the degree values are reported on brain surfaces (left), while the effects of each subcortical region are summarized using a spider plot (right). Abbreviation: FDR, false discovery rate

Classification between patients with migraine and healthy controls

Across 100 repetitions, the sensory/motor cortex, inferior parietal lobule, temporal cortex, and caudate largely contributed to the group classification (Fig. 3A). The mean ± SD values for precision, recall, and AUC were 0.84 ± 0.05, 0.84 ± 0.05, and 0.90 ± 0.03, respectively (Fig. 3B), indicating that these structural connectivity features could be considered as markers for migraine diagnosis.

Fig. 3
figure 3

Classification between patients with migraines and healthy controls.  (A) Selected probabilities based on the LASSO coefficient are shown on brain surfaces. (B) Prediction performance is shown as ROC curves. The red, gray, and blue lines indicate the mean ROC curve across 100 iterations, the ROC curves of all iterations, and the baseline, respectively. Abbreviations: AUC, area under the ROC curve; ROC, receiver operating characteristic; LASSO, least absolute shrinkage and selection operator

Sensitivity analyses

Multiple sensitivity analyses demonstrated the robustness of our findings, as summarized below:

A) Excluding patients with migraine with aura

We repeated the analysis, excluding patients with migraine with aura, and observed that the spatial patterns of the t-statistics aligned with the main findings (linear correlation coefficient r = 0.92, p < 0.001; Supplementary Fig. 1).

B) Different parcellations

We further repeated the analysis using different spatial granularities based on the sub-parcellation of the Desikan–Killiany atlas with 200 and 400 parcels, and observed similar results (r = 0.78, p < 0.001 for 200 parcels; r = 0.91, p < 0.001 for 400 parcels; Supplementary Fig. 2A–B). We also assessed between-group differences using the Schaefer atlas with 300 parcels, which yielded consistent findings (r = 0.76, p < 0.001; Supplementary Fig. 2C). Overall, the results demonstrate the robustness of the findings across different parcellation schemes.

C) Different migraine phases

We examined the between-group differences in manifold eccentricity and degree values of subcortico-cortical structural connectivity in patients with migraines in the interictal, peri-ictal, or ictal phases, finding consistent results across different migraine phases: interictal (r = 0.90, p < 0.001), ictal (r = 0.92, p < 0.001), and peri-ictal (r = 0.52, p < 0.001; Supplementary Fig. 3).

D) Effects of depression and anxiety

We found consistent findings when comparing the statistics of migraine patients with depression or anxiety with those without depression or anxiety (r = 0.83, p < 0.001 for depression; r = 0.72, p < 0.001 for anxiety; Supplementary Fig. 4).

E) Low- vs. high-frequency episodic migraine

We assessed the between-group differences by considering only patients with low-frequency episodic migraines, and observed findings consistent with the main results (r = 0.93, p < 0.001; Supplementary Fig. 5).

F) Manifold eccentricity calculation using multiple eigenvectors

We repeated the assessment of between-group differences in manifold eccentricity using four, seven, and 11 eigenvectors that showed approximately 50%, 70%, and 90% of the explanations, respectively. These findings were largely consistent with the main results derived using the three eigenvectors (r = 0.95, p < 0.001 for four eigenvectors; r = 0.94, p < 0.001 for seven eigenvectors; and r = 0.91, p < 0.001 for 11 eigenvectors; Supplementary Fig. 6). These results indicate that these three eigenvectors may contain sufficient information to explain the connectome data.

G) Between-group differences without controlling for age and sex

When we assessed the between-group differences without controlling for age and sex, we observed similar findings to the initial analysis (r = 1.00, p < 0.001; Supplementary Fig. 7), indicating that age and sex did not significantly affect the results.

Discussion

Migraines are a form of headache characterized by atypical brain structures; however, little is known about the structural connectome organization in patients with episodic migraines. Using manifold learning techniques, in the present study, we identified atypical connectome organization in the somatomotor and limbic networks in patients with migraines. We further demonstrated the utility of manifold features in distinguishing patients with migraines from healthy controls. Our findings provide insights into the understanding of macroscale structural connectome organization in patients with episodic migraines, further shedding light on the neural underpinnings of this condition.

Manifold learning is an approach used to delineate brain organization. In the present study, we applied this approach to structural connectivity data to define structural manifold eccentricity. Previous functional MRI studies have revealed a distinct cortical hierarchy along the sensory-transmodal axis [8, 10]. Similarly, another study based on microstructural profiles revealed a similar pattern known as the sensory-fugal axis [44]. In the present study, we applied this approach to structural connectivity data and defined structural manifold eccentricity. The manifold eccentricity depicts the distance between each brain region and the center of all data points in the manifold space, thereby quantifying changes in connectome organization in terms of expansion and contraction [15]. Prior research has shown that the expansion of manifold eccentricity within the manifold space may be interpreted as increases in nodal and within-module connectivity and segregation with other brain regions [15]. In the present study, we found that the limbic regions, particularly the orbitofrontal cortex and temporal pole, showed expansion of the connectome manifold. These findings may indicate greater differentiation between connectivity with limbic regions compared to connectivity with the rest of the brain. In contrast, the sensorimotor regions exhibited contractions in manifold eccentricity, indicating more integrated patterns in the brain.

The orbitofrontal cortex is involved in pain-related brain networks that mediate changes in the value of pain [45], in which altered serotonergic and dopaminergic systems may disrupt sensory integration and decision-making processes [46, 47]. The temporal pole participates in pain processing by mediating affective responses to stimuli [48,49,50]. In the present study, we found that these limbic cortices showed an expansion of manifold eccentricity in patients with migraines, suggesting more dissimilar connectome patterns compared with other brain regions. Our findings further indicate that expansions in manifold eccentricity in the limbic cortex are related to a larger differentiation in connectivity from other brain regions, which could lead to reduced efficiency in processing sensory and cognitive information. Sensory processing abnormalities have been well described in patients with migraines [51,52,53,54,55]. In particular, the pain matrix, which consists of the thalamic, sensorimotor, and prefrontal regions, integrates pain-related sensory and cognitive responses, and contributes to active nociceptive processing [56, 57]. Although the exact biological underpinnings of the pain matrix remain debatable, it is clear that altered nodal connectivity in sensory regions is associated with pain processing in patients with migraine.

Subcortical regions are also crucial for understanding the brain of patients with migraines. For example, the amygdala receives sensory inputs from the limbic cortices of the orbitofrontal, temporal, and cingulate regions, delivering this information to subcortical structures [58]. These subcortico-cortical circuits play a role in pain modulation, while alterations in these circuits are known to increase pain [58]. In one previous study, patients with episodic migraines showed that altered connectivity in the reward system, specifically between the accumbens and amygdala, might have been affected by endogenous dopaminergic signaling, which could be the underlying neural mechanism for pain processing in patients with migraines [59]. As such, increased degree values in the subcortico-cortical connectivity in the accumbens and amygdala may indicate increased reception of pain information from the cortical areas.

Multiple sensitivity analyses considering different spatial scales, types of parcellation scheme, the number of eigenvectors used to construct manifold eccentricity, and removal of potential confounding effects by considering migraine with/without aura and ictal/peri-ictal/interictal phases have demonstrated the reliability and robustness of our findings. We found consistent between-group differences across different migraine phases, indicating the role of the structural connectome as a reliable diagnostic marker for migraine. Furthermore, a machine learning framework for classifying patients with migraines and healthy controls validated the clinical utility of our approach. Although the generalizability of the findings should be tested using independent external datasets in future studies, the present results enhance our understanding of the macroscale perspectives of the brains of migraine patients.

In the present study, we employed manifold learning techniques to assess the changes in the structural connectome organization using diffusion MRI-based tractography data. We thus provide insights into the structural connectivity in patients with migraines, extending beyond the functional alterations explored in our previous study [14]. This study contributes to a deeper understanding of whole-brain structural connectivity in patients with migraines, potentially facilitating the discovery of robust neuroimaging biomarkers.

This study has some limitations which should be considered. First, we did not assess the effect of the white matter hyperintensity, which has been shown to be associated with migraine pathology [60, 61]. This parameter can be measured using fluid-attenuated inversion recovery data, and future studies are required to investigate the relationship between white matter hyperintensity burden and structural connectome changes. Second, our analysis focused on identifying structural connectome changes in patients with migraines and their implications for group classification. In future studies, we will further explore how these changes are related to headaches. Third, several methodological choices can be considered for constructing structural connectivity; these may include: tractography approaches, either deterministic or probabilistic; filtering methods, including spherical deconvolution informed filtering of tractograms [28, 62] and convex optimization modeling for microstructure informed tractography [63, 64]; and quantification of structural connectivity, such as fiber count, density, or cross-section streamlines. Fourth, we only studied the structural connectivity changes in patients with migraines in the current study. Further investigation of the interactions between the structural and functional networks may provide new insights into migraine connectopathy. Finally, the current study explored connectivity changes in patients with episodic migraine. Comparisons of connectome changes between patients with episodic and chronic migraines, as well as the effects of psychological disorders, should be investigated further.

Data availability

All data supporting the findings of this study will be made available upon request from the corresponding author, Mi Ji Lee. The data are not publicly available due to Institutional Review Board restrictions. The codes for eigenvector generation are available at https://github.com/MICA-MNI/BrainSpace; the codes for calculating network communication measures are available at https://sites.google.com/site/bctnet/; and the codes for statistical analyses are available at https://github.com/MICA-MNI/ENIGMA.

References

  1. Walter K (2022) What Is Migraine? JAMA 327:93. https://doi.org/10.1001/jama.2021.21857

  2. Hasan H, Irfan Waheed R, Bin Arif T et al (2020) Gray and White Matter Changes in migraineurs: a review of literature. SN Compr Clin Med 2:2185–2196. https://doi.org/10.1007/s42399-020-00576-7

    Article  Google Scholar 

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

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

  5. Planchuelo-Gómez Á, García-Azorín D, Guerrero ÁL et al (2019) 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 

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

  7. Bernhardt BC, Smallwood J, Keilholz S, Margulies DS (2022) Gradients in brain organization. NeuroImage 251:118987. https://doi.org/10.1016/j.neuroimage.2022.118987

    Article  PubMed  Google Scholar 

  8. Margulies DS, Ghosh SS, Goulas A et al (2016) Situating the default-mode network along a principal gradient of macroscale cortical organization. Proceedings of the National Academy of Sciences 113:12574–12579

  9. Haak KV, Marquand AF, Beckmann CF (2018) Connectopic mapping with resting-state fMRI. NeuroImage 170:83–94. https://doi.org/10.1016/j.neuroimage.2017.06.075

    Article  PubMed  Google Scholar 

  10. Vos de Wael R, Benkarim O, Paquola C et al (2020) BrainSpace: a toolbox for the analysis of macroscale gradients in neuroimaging and connectomics datasets. Commun Biol 3:103. https://doi.org/10.1038/s42003-020-0794-7

    Article  PubMed  PubMed Central  Google Scholar 

  11. Cross N, Paquola C, Pomares FB et al (2021) Cortical gradients of functional connectivity are robust to state-dependent changes following sleep deprivation. NeuroImage 226:117547. https://doi.org/10.1016/j.neuroimage.2020.117547

    Article  PubMed  Google Scholar 

  12. Vos de Wael R, Royer J, Tavakol S et al (2021) Structural connectivity gradients of the temporal lobe serve as Multiscale axes of Brain Organization and cortical evolution. Cereb Cortex 31:5151–5164. https://doi.org/10.1093/cercor/bhab149

    Article  PubMed  PubMed Central  Google Scholar 

  13. Royer J, Larivière S, Rodriguez-Cruces R et al (2023) Cortical microstructural gradients capture memory network reorganization in temporal lobe epilepsy. Brain 146:3923–3937. https://doi.org/10.1093/brain/awad125

    Article  PubMed  PubMed Central  Google Scholar 

  14. Lee CH, Park H, Lee MJ, Park B (2023) Whole-brain functional gradients reveal cortical and subcortical alterations in patients with episodic migraine. Hum Brain Mapp 44:2224–2233. https://doi.org/10.1002/hbm.26204

    Article  PubMed  PubMed Central  Google Scholar 

  15. Park B, Bethlehem RAI, Paquola C et al (2021) An expanding manifold in transmodal regions characterizes adolescent reconfiguration of structural connectome organization. Elife 10:e64694. https://doi.org/10.7554/eLife.64694

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Park B, Byeon K, Park H (2019) FuNP (Fusion of Neuroimaging Preprocessing) Pipelines: a fully automated Preprocessing Software for functional magnetic resonance imaging. Front Neuroinform 13

  17. Fischl B (2012) FreeSurfer Neuroimage 62:774–781. https://doi.org/10.1016/j.neuroimage.2012.01.021

    Article  PubMed  Google Scholar 

  18. Dale AM, Fischl B, Sereno MI (1999) Cortical surface-based analysis: I. Segmentation and Surface Reconstruction. NeuroImage 9:179–194. https://doi.org/10.1006/nimg.1998.0395

    Article  CAS  PubMed  Google Scholar 

  19. Fischl B, Sereno MI, Dale AM (1999) Cortical surface-based analysis: II: inflation, flattening, and a surface-based Coordinate System. NeuroImage 9:195–207. https://doi.org/10.1006/nimg.1998.0396

    Article  CAS  PubMed  Google Scholar 

  20. Fischl B, Sereno MI, Tootell RBH, Dale AM (1999) High-resolution intersubject averaging and a coordinate system for the cortical surface. Hum Brain Mapp 8:272–284

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

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

    Article  PubMed  Google Scholar 

  23. Avants BB, Tustison NJ, Song G et al (2011) A reproducible evaluation of ANTs similarity metric performance in brain image registration. NeuroImage 54:2033–2044. https://doi.org/10.1016/j.neuroimage.2010.09.025

    Article  PubMed  Google Scholar 

  24. Christiaens D, Reisert M, Dhollander T et al (2015) Global tractography of multi-shell diffusion-weighted imaging data using a multi-tissue model. NeuroImage 123:89–101. https://doi.org/10.1016/j.neuroimage.2015.08.008

    Article  PubMed  Google Scholar 

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

  26. Tournier JD, Calamante F, Connelly A (2010) Improved probabilistic streamlines tractography by 2nd order integration over fibre orientation distributions. In: Proceedings of the international society for magnetic resonance in medicine. John Wiley & Sons, Inc, New Jersey

  27. Tournier J, Calamante F, Connelly A (2012) MRtrix: diffusion tractography in crossing fiber regions. Int J Imaging Syst Technol 22:53–66

    Article  Google Scholar 

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

    Article  PubMed  Google Scholar 

  29. Desikan RS, Ségonne F, Fischl B et al (2006) An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. NeuroImage 31:968–980. https://doi.org/10.1016/j.neuroimage.2006.01.021

    Article  PubMed  Google Scholar 

  30. Cruces RR, Royer J, Herholz P et al (2022) Micapipe: a pipeline for multimodal neuroimaging and connectome analysis. NeuroImage 263:119612. https://doi.org/10.1016/j.neuroimage.2022.119612

    Article  PubMed  Google Scholar 

  31. Betzel RF, Griffa A, Hagmann P, Mišić B (2019) Distance-dependent consensus thresholds for generating group-representative structural brain networks. Netw Neurosci 3:475–496. https://doi.org/10.1162/netn_a_00075

    Article  PubMed  PubMed Central  Google Scholar 

  32. Coifman RR, Lafon S (2006) Diffusion maps. Appl Comput Harmon Anal 21:5–30. https://doi.org/10.1016/J.ACHA.2006.04.006

    Article  Google Scholar 

  33. Tenenbaum JB, de Silva V, Langford JC (2000) A Global Geometric Framework for Nonlinear Dimensionality Reduction. Science (1979) 290:2319–2323. https://doi.org/10.1126/science.290.5500.2319

  34. Langs G, Golland P, Ghosh SS (2015) Predicting Activation Across individuals with resting-state functional connectivity based Multi-atlas label Fusion. In: Navab N, Hornegger J, Wells WM, Frangi A (eds) Medical Image Computing and Computer-assisted intervention -- MICCAI 2015. Springer International Publishing, Cham, pp 313–320

    Google Scholar 

  35. Park B, Hong S-J, Valk SL et al (2021) Differences in subcortico-cortical interactions identified from connectome and microcircuit models in autism. Nat Commun 12:2225

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Yang Y, Park B (2023) Differences in structural connectome organization across sleep quality. Heliyon 9

  37. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc: Ser B (Methodol) 57:289–300

    Article  Google Scholar 

  38. Thomas Yeo BT, Krienen FM, Sepulcre J et al (2011) The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J Neurophysiol 106:1125–1165. https://doi.org/10.1152/jn.00338.2011

    Article  PubMed Central  Google Scholar 

  39. Mesulam MM (1998) From sensation to cognition. Brain 121:1013–1052. https://doi.org/10.1093/brain/121.6.1013

    Article  PubMed  Google Scholar 

  40. Patenaude B, Smith SM, Kennedy DN, Jenkinson M (2011) A bayesian model of shape and appearance for subcortical brain segmentation. NeuroImage 56:907–922

    Article  PubMed  Google Scholar 

  41. Pedregosa F (2011) Scikit-learn: machine learning in Python. J Mach Learn Res 12:2825

    Google Scholar 

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

    Article  PubMed  Google Scholar 

  43. Yoo S, Jang Y, Hong S-J et al (2024) Whole-brain structural connectome asymmetry in autism. NeuroImage 120534

  44. Paquola C, De Vos R, Wagstyl K et al (2019) Microstructural and functional gradients are increasingly dissociated in transmodal cortices. PLoS Biol 17:e3000284

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Becker S, Gandhi W, Pomares F et al (2017) Orbitofrontal cortex mediates pain inhibition by monetary reward. Soc Cogn Affect Neurosci 12:651–661. https://doi.org/10.1093/scan/nsw173

    Article  PubMed  PubMed Central  Google Scholar 

  46. Lev R, Granovsky Y, Yarnitsky D (2010) Orbitofrontal disinhibition of pain in migraine with aura: an interictal EEG-mapping study. Cephalalgia 30:910–918. https://doi.org/10.1177/0333102409357249

    Article  PubMed  Google Scholar 

  47. Rule RR, Shimamura AP, Knight RT (2002) Orbitofrontal cortex and dynamic filtering of emotional stimuli. Cogn Affect Behav Neurosci 2:264–270. https://doi.org/10.3758/CABN.2.3.264

    Article  PubMed  Google Scholar 

  48. Schwedt TJ (2013) Multisensory integration in migraine. Curr Opin Neurol 26

  49. Moulton EA, Becerra L, Maleki N et al (2011) Painful heat reveals hyperexcitability of the temporal Pole in Interictal and Ictal Migraine States. Cereb Cortex 21:435–448. https://doi.org/10.1093/cercor/bhq109

    Article  CAS  PubMed  Google Scholar 

  50. Demarquay G, Royet JP, Mick G, Ryvlin P (2008) Olfactory hypersensitivity in migraineurs: a H215O-PET study. Cephalalgia 28:1069–1080. https://doi.org/10.1111/j.1468-2982.2008.01672.x

    Article  CAS  PubMed  Google Scholar 

  51. Goadsby PJ, Holland PR, Martins-Oliveira M et al (2017) Pathophysiology of migraine: a disorder of sensory Processing. Physiol Rev 97:553–622. https://doi.org/10.1152/physrev.00034.2015

    Article  PubMed  PubMed Central  Google Scholar 

  52. Harriott AM, Schwedt TJ (2014) Migraine is Associated with altered Processing of sensory stimuli. Curr Pain Headache Rep 18:458. https://doi.org/10.1007/s11916-014-0458-8

    Article  PubMed  Google Scholar 

  53. Peng K-P, May A (2019) Migraine understood as a sensory threshold disease. Pain 160

  54. Meylakh N, Henderson LA (2022) Exploring alterations in sensory pathways in migraine. J Headache Pain 23:5. https://doi.org/10.1186/s10194-021-01371-y

    Article  PubMed  PubMed Central  Google Scholar 

  55. Hodkinson DJ, Veggeberg R, Kucyi A et al (2016) Cortico–cortical connections of primary sensory areas and Associated symptoms in Migraine. https://doi.org/10.1523/ENEURO.0163-16.2016. eNeuro 3:ENEURO.0163-16.2016

  56. Derbyshire SWG, Jones AKP, Gyulai F et al (1997) Pain processing during three levels of noxious stimulation produces differential patterns of central activity. Pain 73

  57. Tracey I (2008) Imaging pain. Br J Anaesth 101:32–39. https://doi.org/10.1093/bja/aen102

    Article  CAS  PubMed  Google Scholar 

  58. Goddard GV (1964) Functions of the amygdala. Psychol Bull 62:89

    Article  CAS  PubMed  Google Scholar 

  59. Kim DJ, Jassar H, Lim M et al (2021) Dopaminergic regulation of reward system connectivity underpins Pain and Emotional suffering in Migraine. J Pain Res 14:631–643. https://doi.org/10.2147/JPR.S296540

    Article  PubMed  PubMed Central  Google Scholar 

  60. Lee MJ, Park B-Y, Cho S et al (2019) Cerebrovascular reactivity as a determinant of deep white matter hyperintensities in migraine. Neurology 92:e342–e350. https://doi.org/10.1212/WNL.0000000000006822

    Article  CAS  PubMed  Google Scholar 

  61. Lee MJ, Park B, Cho S et al (2022) Cerebrovascular reactivity and deep white matter hyperintensities in migraine: a prospective CO2 targeting study. J Cereb Blood Flow Metabolism 42:1879–1889. https://doi.org/10.1177/0271678X221103006

    Article  CAS  Google Scholar 

  62. Smith RE, Tournier J-D, Calamante F, Connelly A (2013) SIFT: spherical-deconvolution informed filtering of tractograms. NeuroImage 67:298–312

    Article  PubMed  Google Scholar 

  63. Daducci A, Dal Palù A, Lemkaddem A, Thiran J-P (2014) COMMIT: convex optimization modeling for microstructure informed tractography. IEEE Trans Med Imaging 34:246–257

    Article  PubMed  Google Scholar 

  64. Schiavi S, Ocampo-Pineda M, Barakovic M et al (2020) A new method for accurate in vivo mapping of human brain connections using microstructural and anatomical information. Sci Adv 6:eaba8245

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Funding

This research was supported by the National Research Foundation of Korea (NRF-2020R1A2B5B01001826; NRF-2022R1A5A7033499), Institute for Information and Communications Technology Planning and Evaluation (IITP) funded by the Korea Government (MSIT) (No. 2022-0-00448, Deep Total Recall: Continuous Learning for Human-Like Recall of Artificial Neural Networks; No. RS-2022-00155915, Artificial Intelligence Convergence Innovation Human Resources Development (Inha University); No. RS-2021-II212068, Artificial Intelligence Innovation Hub), Institute for Basic Science (IBS-R015-D1), and the Korea Medical Device Development Fund grant funded by the Korean government (the Ministry of Science and ICT, the Ministry of Trade, Industry and Energy, the Ministry of Health & Welfare, the Ministry of Food and Drug Safety) (Project number: RS-2023-00229484).

Author information

Authors and Affiliations

Authors

Contributions

E.N. and B.P. designed the study, analyzed the data, and wrote the manuscript. M.J.L. designed the study, acquired and interpreted the data, and critically revised the manuscript to ensure important intellectual content. J.Y.N., Y.P., and Y.J. aided with the experiments and reviewed the manuscript. M.J.L. and B.P. are the corresponding authors of this study and are responsible for the integrity of the data analysis. All authors have reviewed and approved the final version of the manuscript for publication.

Corresponding authors

Correspondence to Mi Ji Lee or Bo-yong Park.

Ethics declarations

Conflict of interest

All authors declare no conflicts of interest.

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 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Noh, E., Namgung, J.Y., Park, Y. et al. Shifts in structural connectome organization in the limbic and sensory systems of patients with episodic migraine. J Headache Pain 25, 99 (2024). https://doi.org/10.1186/s10194-024-01806-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s10194-024-01806-2

Keywords