White matter changes in chronic and episodic migraine: a diffusion tensor imaging study
The Journal of Headache and Pain volume 21, Article number: 1 (2020)
White matter alterations have been observed in patients with migraine. However, no microstructural white matter alterations have been found particularly in episodic or chronic migraine patients, and there is limited research focused on the comparison between these two groups of migraine patients.
Fifty-one healthy controls, 55 episodic migraine patients and 57 chronic migraine patients were recruited and underwent brain T1-weighted and diffusion-weighted MRI acquisition. Using Tract-Based Spatial Statistics (TBSS), fractional anisotropy, mean diffusivity, radial diffusivity and axial diffusivity were compared between the different groups. On the one hand, all migraine patients were compared against healthy controls. On the other hand, patients from each migraine group were compared between them and also against healthy controls. Correlation analysis between clinical features (duration of migraine in years, time from onset of chronic migraine in months, where applicable, and headache and migraine frequency, where applicable) and Diffusion Tensor Imaging measures was performed.
Fifty healthy controls, 54 episodic migraine and 56 chronic migraine patients were finally included in the analysis. Significant decreased axial diffusivity (p < .05 false discovery rate and by number of contrasts corrected) was found in chronic migraine compared to episodic migraine in 38 white matter regions from the Johns Hopkins University ICBM-DTI-81 White-Matter Atlas. Significant positive correlation was found between time from onset of chronic migraine and mean fractional anisotropy in the bilateral external capsule, and negative correlation between time from onset of chronic migraine and mean radial diffusivity in the bilateral external capsule.
These findings suggest global white matter structural differences between episodic migraine and chronic migraine. Patients with chronic migraine could present axonal integrity impairment in the first months of chronic migraine with respect to episodic migraine patients. White matter changes after the onset of chronic migraine might reflect a set of maladaptive plastic changes.
According to the IIIrd edition of the International Classification of Headache Disorders (ICHD-3), patients with Chronic Migraine (CM) suffer from headache during 15 or more days per month for more than 3 months, with at least eight of these days with migrainous characteristics . Between 2 and 3% of migraine patients evolve annually from Episodic Migraine (EM) to CM . Some risk factors have been associated with progression from EM to CM, but the pathophysiological mechanisms of this conversion remain to be elucidated. Moreover, EM and CM could represent either two ranges of the same entity, or two subgroups with distinctive characteristics. Previous review studies have exposed that important distinctions exist along the continuum between EM and CM [3, 4], but it is not clear whether CM is a kind of more frequent EM or a distinct entity.
Magnetic Resonance Imaging (MRI) is one of the most powerful technologies available for the study of the migrainous brain. Among its different modalities, diffusion MRI (dMRI) is particularly well suited for the analysis of possible white matter alterations in migraine. dMRI studies have shown changes affecting the white matter in migraineurs with respect to healthy controls; these changes were mainly observed in the corpus callosum [5,6,7,8,9,10,11], thalamus [7, 8, 12, 13], thalamic radiation [7,8,9, 14, 15] and cingulate gyrus [9, 10, 15,16,17].
Some studies using dMRI have focused on migraine with and without aura. Even though significant differences have been found [10, 18], using MRI information it remains unclear whether migraine with and without aura area are actually two distinct entities, the manifestations of the same pathophysiological substrate on two different phenotypes [19, 20], or even whether aura can be better defined as a migraine phase, but not a category such as chronic and episodic migraine .
Some other studies have employed dMRI to investigate CM. To the best of our knowledge, however, only one study by Neeb et al. searched for possible differences between patients with EM and CM, not finding any significant differences . Other studies compared chronic migraine patients with healthy controls  or mixed EM and CM patients in order to compare migraine patients and controls . Previous studies alternatively used T1-weighted MRI images [23, 24] and functional MRI , respectively, to compare possible differences between patients with CM and EM. These studies described reduced grey matter volume in CM compared to EM  and more activity in CM compared to EM patients with headaches during the scanning in the anterior hypothalamus .
The present study performs a detailed comparison of the white matter in EM, CM and healthy controls over a large cohort of subjects, using dMRI data. We hypothesised that there could be white matter structural differences between CM and EM patients, although no significant differences were found in a previous study . In our study, we included a considerably higher number of patients and controls compared to that study. Our goals are to:
Investigate whether there are significant differences between CM and EM, and between these groups and healthy controls. To that end, TBSS  was employed as dMRI analysis technique.
Examine how white matter descriptors based on dMRI relate to clinical features in migraine patients. To that end, a correlation analysis between dMRI parameters and variables such as duration of migraine and time from onset of CM was performed.
Materials and methods
We conducted an observational analytic study with a case-control design. The target population included patients with migraine. Patients were firstly screened and recruited from the headache outpatient unit at the Hospital Clínico Universitario de Valladolid (Valladolid, Spain), a tertiary centre that receives patients both from specialized care and directly from primary care. Inclusion criteria were: a) Diagnosis of episodic or chronic migraine according to the ICHD-3 beta and ICHD-3 criteria [1, 27]; b) no changes in the situation of episodic or chronic migraine in the previous 3 months; c) agreeing to participate and signing the Informed Consent; d) aged between 18 and 60. We excluded patients with a) high frequency Episodic Migraine, suffering between 10 and 14 headache days per month (to avoid any confusion between high frequency EM and CM ); b) other non-craniofacial painful conditions occurring 10 or more days per month; c) known major psychiatric diseases (described in anamnesis or presence of Depression or Anxiety according to Hospital Anxiety and Depression Scale ); d) other neurological diseases; e) drug or substance abuse; f) pregnancy or childbearing; g) other headache disorders. All the patients scanned were preventive naïve. Patients were asked to keep a headache diary for 3 months (before inclusion) and were classified as EM when they had less than 10 headaches per month, or CM according to ICHD-3 criteria. Participants with EM were not allowed to have headache days with tension-type headache (TTH) phenotype. No healthy controls (HC) were included if they showed a present or past history of migraine, or if any other neurological or psychiatric condition was present, with the sole exception of infrequent TTH. We used a non-probabilistic sampling method by convenience sampling. Healthy controls balanced for age and sex were recruited through hospital and University colleagues and advertisements in these facilities by convenience sampling and snowball sampling.
For all patients, sociodemographic and clinical data were collected, including the duration of migraine disease (years), headache and migraine frequency (days per month) and time from the onset of chronic migraine (months) when applicable. The intake of symptomatic medication, i.e., combination of analgesics and triptan, was considered to verify if patients fulfilled criteria of acute medication overuse (intake at 10 or more days per month). Presence of aura was drawn.
The local Ethics Committee of Hospital Clínico Universitario de Valladolid approved the study (PI: 14–197). All participants read and signed a written consent form prior to their participation.
Images were acquired for migraine patients after at least 24 h from the last migraine attack. High-resolution 3D T1-weighted and diffusion-weighted MRI data were acquired using a Philips Achieva 3 T MRI unit (Philips Healthcare, Best, The Netherlands) with a 32-channel head coil in the MRI facility at the Universidad de Valladolid (Valladolid, Spain).
For the anatomical T1-weighted images, the following acquisition parameters were used: Turbo Field Echo (TFE) sequence, repetition time (TR) = 8.1 ms, echo time (TE) = 3.7 ms, flip angle = 8°, 256 × 256 matrix size, 1 × 1 × 1 mm3 of spatial resolution and 160 slices covering the whole brain.
Diffusion-weighted images (DWI) were acquired using the next parameters: TR = 9000 ms, TE = 86 ms, flip angle = 90°, 61 gradient directions, one baseline volume, b-value = 1000 s/mm2, 128 × 128 matrix size, 2 × 2 × 2 mm3 of spatial resolution and 66 axial slices covering the whole brain.
T1 and diffusion-weighted scans were acquired during the same session, starting with the T1 scan followed by the diffusion-weighted scan, between May 2014 and July 2018. Total acquisition time for each subject was around 18 min.
MR images were processed before carrying out the statistical analysis using TBSS . For the TBSS analysis, four Diffusion Tensor Imaging (DTI) measures were obtained: Fractional Anisotropy (FA), Mean Diffusivity (MD), Radial Diffusivity (RD) and Axial Diffusivity (AD). In a nutshell, FA reflects the degree of directionality of water diffusivity, MD is a global measure of water diffusion, RD quantifies the diffusion perpendicular to the principal direction and AD is the diffusion in the main direction of the white matter fibres .
Prior to the obtention of the four DTI measures, diverse preprocessing procedures were implemented on the DWI data. Diffusion-weighted images were denoised, using “dwidenoise” tool from MRtrix [31, 32], eddy currents and motion corrected, using “dwipreproc” tool from MRtrix , and B1 field inhomogeneity corrected, using “dwibiascorrect” tool with the “-fast” option from MRtrix [34, 35].
Once the DWI images were preprocessed, a whole brain mask for each image was generated using “dwi2mask” tool from MRtrix  and, next, diffusion tensors at each voxel were estimated using the “dtifit” tool from FSL , also obtaining FA, MD and AD maps. RD was manually calculated by obtaining the mean of the second and the third eigenvalues, which were also previously computed with “dtifit”.
For the TBSS method, all participants’ FA images were nonlinearly registered using the FNIRT tool from FSL to a template of the averaged FA images (FMRIB-58) in Montreal Neurological Institute (MNI) space; the FNIRT tool uses a b-spline representation of the registration warp field . After registration, a mean FA image was generated and thinned to create a mean FA skeleton of white matter tracts using a FA value of .2 as threshold to distinguish white from grey matter. Then, each subject’s aligned FA images were projected onto the mean FA skeleton. In a similar way, the TBSS process was repeated for MD, AD and RD, using the protocol devoted for non-FA images. To identify the white matter tracts, the Johns Hopkins University ICBM-DTI-81 White-Matter Labels Atlas [39, 40] provided in the FSL toolbox was used. However, this atlas does not cover the whole white matter across brain, so we also employed the Johns Hopkins University White-Matter Tractography Atlas , which contains a lower number of tracts but covers areas not included in the other atlas. The minimum volume to consider significant results in a region was set to 30 mm3. It must be noted that, because we use an image of 1 mm3 in the MNI space to identify the regions, the volume in mm3 is equal to the number of voxels. Moreover, when we extract the significant results in a region, we consider all the significant voxels from that region, from one or more clusters.
We estimated Sample Size according to Chong and Schwedt, 2015 . Based on the FA results from this study in three major tracts, we calculated a worst possible scenario model with an estimated effect size of difference between groups of .02 (greatest difference between groups) and a variance of .003 (greatest single group variance); a type 1-error rate of 1% and 80% power and anticipating a proportion of 10% of lost patients. The expected sample size was 167 participants.
Kolmogorov-Smirnov and Levene’s Test for equality of variances tests were used to assess normality and homogeneity of variance in age and duration of migraine in years. To test for significant differences in the age of the three groups, a one-way ANOVA was used if the null hypothesis in Kolmogorov-Smirnov and Levene tests was not rejected; otherwise, Kruskal-Wallis test was employed. To test for significant gender differences, a chi-square test was used. To compare continuous clinical features between migraine patients (i.e., duration of migraine history in years for both groups of patients and time from onset of chronic migraine in months for chronic migraine patients), a two-tailed unpaired t-test was used if the null hypothesis in Kolmogorov-Smirnov and Levene tests was not rejected; otherwise, Mann-Whitney U test was employed. To compare categorical clinical features between migraine patients, Fisher’s exact test was employed.
We executed group-wise comparisons of all migraineurs vs. healthy controls, CM vs. EM, EM vs. healthy controls, and CM vs. healthy controls. The voxel-wise TBSS differences in FA, MD, AD and RD values of white matter between the different groups were tested using a permutation-based inference tool by nonparametric statistics called “randomise”, implemented in FSL, with the threshold-free cluster enhancement (TFCE) option [42, 43]. Five thousand permutations were set to allow robust statistical inference and the significance threshold for intergroup differences was p < .05 after correcting for family wise error (FWE) applying the TFCE option. Additional clinical covariates were added to the comparisons in the cases where significant differences were found. These covariates were analysed individually to evaluate the individual effect of each covariate. In the case of presence of aura, we repeated the original TBSS analysis excluding the patients with migraine with aura. None of the design matrices included duration of migraine and time from onset of CM simultaneously as covariates due to collinearity. Time from onset of CM was also included as a covariate (only in comparisons with CM) because it may correct the results in CM patients in a more meaningful way than the total duration of migraine.
In the cases where significant differences were found, we performed a post-hoc analysis. We applied a false discovery rate (FDR) correction, using the “fdr” command from FSL, to the TFCE uncorrected p-values. The fdr command provides the uncorrected p-value which sets the level of statistical significance after the FDR correction. To correct for number of contrasts in each case, the final level of statistical significance is equal to the uncorrected p-value from the previous step divided by the number of comparisons (Bonferroni correction).
Effect size was computed using Cohen’s d value in regions with significant results from the first analysis (FWE-corrected). For every comparison, the mean value of the “most disabled” group was subtracted from the mean value of the “least disabled” or the control group. In the comparisons between both types of migraine, CM is considered the most disabled group, and EM the least disabled group.
To study the relationship between clinical parameters and DTI measures, Spearman’s rank correlation coefficient was employed in a ROI-based correlation analysis. Duration of the migraine in both types of migraine patients, time from onset of CM in chronic migraineurs, and headache and migraine frequency for both types of migraine patients were the analysed clinical parameters. It must be clarified that we obtained correlation values in CM and EM patients separately, in order to assess differences or trends within each type of migraine. Our intention was to determine the possible effect of headache or migraine frequency, in the specific range of episodic or chronic migraine, and the possible relationship with time, with special attention to the time from onset of CM. To obtain individual label maps for each subject, the inverse warp fields of the FA images to the MNI image transformation from the TBSS procedure were computed and applied to the Johns Hopkins University ICBM-DTI-81 White Matter Atlas. The ROIs that were selected for the correlation analysis were those for which significant differences were found between at least two groups in any diffusion parameter in the TBSS analysis. All DTI measures (FA, MD, RD and AD) were considered for the correlation analysis. We used all DTI measures to avoid a possible loss of complementary information given by each parameter in the assessment of differences within CM or EM. To correct for multiple comparisons, the Benjamini-Hochberg  FDR procedure was applied, and, after this correction, the level of statistical significance was set at p < .05.
During the study period, 51 healthy controls, 55 episodic migraine patients and 57 chronic migraine patients were recruited for the study after matching the inclusion and exclusion criteria. No significant structural abnormality was detected in conventional MRI studies. Due to erroneous results after applying the nonlinear registration to the FMRIB-58 image in MNI space, one healthy control, one episodic migraineur and one chronic migraineur were excluded from the study. Demographic and clinical data for the three groups with the remaining participants are summarised in Table 1. Significant differences were found in duration of migraine history in years between the two migraine groups and, as expected, in headache and migraine frequency between the migraine groups.
Considering the significant differences in duration of migraine history between episodic and chronic migraine patients (Table 1), TBSS analysis was repeated including the duration of migraine history as a covariate, as mentioned in the Statistical Analysis section. Presence of aura was added as an additional covariate to duration of migraine history in a posterior analysis. In the case of CM patients, the TBSS analysis was additionally accomplished including the time from onset of CM as a covariate. This covariate was included in a separate analysis from the one with duration of migraine history due to collinearity, as mentioned previously in the Statistical Analysis section.
Results uncorrected for covariates
No significant differences were found in any of the diffusion indices (FA, MD, RD and AD) between all migraineurs and HC. Dividing all migraineurs into EM and CM, no significant differences were found with respect to HC. With regard to the comparison between CM and EM, no significant differences were found in FA, MD or RD. However, significant lower AD values were found in CM compared to EM in widespread locations across the white matter. These locations correspond to 38 different regions from the ICBM-DTI-81 White Matter Atlas, and six regions from the White Matter Tractography Atlas, and are shown in Tables 2 and 3 and Fig. 1. The FWE-corrected results can be seen in Additional file 1: Figure S1 and Table S1.
In the case of the analysis including only patients with migraine without aura, no new significant results were observed. As in the original sample, significant lower AD values were found in CM compared to EM (FWE-corrected results), but in 10 regions from the ICBM-DTI-81 White Matter Atlas, all of them included in the 38 regions with significant differences. The FWE-corrected results excluding migraine with aura patients can be seen in Additional file 1: Table S2.
Results corrected for covariates
After the post-hoc analysis, significant lower AD values were found in CM compared to EM in six regions from the ICBM-DTI-81 White Matter Atlas, and one region (the right corticospinal tract) from the White-Matter Tractography Atlas, when including the time from onset of CM as a covariate. These results are shown in Tables 4 and 5 and Fig. 1. No significant results were observed when correcting for total duration of migraine or for the other DTI measures.
In the following subsections, the FWE-corrected results for the diverse covariates are shown.
Results corrected for duration of migraine history
Including the duration of migraine history as a covariate, no significant differences were found for FA, MD or RD. Significant decreased AD values in CM with respect to EM were found only in the middle cerebellar peduncle (675 mm3 with p < .05 FWE-corrected, minimum p-value = .028), but other 26 regions remained with p < .1 FWE-corrected (Additional file 1: Table S3). Also, significant increased AD values in EM with respect to HC were found in seven regions from the left hemisphere (Additional file 1: Table S4). These results can be seen in Additional file 1: Figure S2. No significant differences were found between all migraineurs together and HC, or between CM and HC.
Adding the presence of aura as a covariate, significant increased AD values in EM with respect to HC were found in seven regions from the left hemisphere. These seven regions were the same regions that showed significant differences in the analysis including only the duration of migraine history as a covariate. No significant differences were found for FA, MD or RD.
Results corrected for time from onset of CM
In the additional comparisons for CM patients, including time from onset of CM as a covariate, no significant differences were found for MD or RD. Significant decreased AD values in CM with respect to EM were found in 23 regions with p < .05 FWE-corrected (Additional file 1: Table S5). Significant decreased FA values in CM compared to HC were found in 15 regions, most of them from the right hemisphere (Additional file 1: Table S6). These results are shown in Additional file 1: Figure S3.
When comparing between CM and EM Axial Diffusivity values, all Cohen’s d values (except for left fornix) were negative, which means that AD values were lower in CM. For the middle cerebellar peduncle, left external capsule and right sagittal stratum, the Cohen’s d absolute values were equal or greater than .5, a medium effect size according to . A very similar trend, but with lower Cohen’s d absolute values, was found in CM with respect to HC. In the comparison between EM and HC, Cohen’s d value was positive in almost all regions, i.e., AD values were higher in EM with respect to HC. In the pontine crossing tract, a medium effect size was obtained (d = .59). These results are depicted in Additional file 1: Table S7 and Figures S4 and S5. Results for the other DTI measures can be seen in Additional file 1: Figures S6, S7, S8, S9, S10 and S11.
A summary of these results can be seen in Table 6.
After multiple comparisons correction, ROI-based significant positive correlations between time from onset of chronic migraine and mean FA in the right (ρ = .420, p = .001) and left (ρ = .439, p < .001) external capsule were found. Significant negative correlations between time from onset of chronic migraine and mean RD in the right (ρ = −.427, p = .001) and left (ρ = −.439, p < .001) external capsule were found. These results can be seen in Fig. 2.
No significant correlations were found neither for mean MD or mean AD, nor duration of the migraine, headache and migraine frequency for both migraine groups in ROI-based correlation analysis.
In a TBSS-based dMRI analysis, in relation with the first stated objective of the study, white matter structural changes in Chronic Migraine compared to Episodic Migraine patients were found in 38 regions when AD was considered as DTI measure. These findings suggest global white matter changes in CM compared to EM.
In the analysis excluding patients with migraine with aura, no new results were observed with respect to the analysis with the whole sample, and the regions with significant results were a subset of the regions of the original analysis. Using only patients without aura, higher p-values and lower number of regions with significant differences are obtained. These results suggest that the differences are caused by a loss of statistical power more than the effect of aura itself. Including duration of migraine history as a covariate, differences were found between EM and Healthy Controls in seven regions, but they did not survive the FDR and the number of contrasts corrections. The addition of the presence of aura as a covariate did not change the results with respect to using only the duration of migraine history as a covariate.
Additionally, including time from onset of chronic migraine as a covariate in the comparisons with CM patients, significant decreased values were found in AD between CM and EM. Decreased values were also found in FA in CM with respect to HC, but they did not survive the FDR and the number of contrasts corrections.
In relation with the second objective of the study, correlation analysis was executed between diverse DTI measures and clinical features. Significant correlations between time from onset of chronic migraine and mean FA (positive correlation) and mean RD (negative correlation) in the bilateral external capsule were found.
White matter differences between both groups of migraineurs were not previously found by Neeb et al. in the only study, to the best of our knowledge, which assessed patients with chronic and episodic migraine using DTI . In the present study, however, a considerably bigger cohort was included when compared to that study, and our participants were considerably younger. The influence of age in the white matter diffusion is well-known and has been extensively assessed. The most common pattern is to find decreased FA and increased RD values in older people, while the AD pattern is unclear [30, 46]. Considering the aging-effect, altered DTI measures in healthy people might reduce the diffusion differences between patients and controls in older subjects.
Differences between the two migraine groups were found using AD as a DTI measure. It is known, however, that relationships between DTI-derived parameter changes and specific microstructure alterations are difficult to establish, and therefore results must be interpreted carefully. Winklewski et al. interpreted reduced AD values, as we found in CM compared to EM, as the beginning of demyelination . This reduction, nevertheless, could be ineffective to detect prolonged demyelination . Nonetheless, based on studies in mice, a review study by Alexander et al. linked AD more to axonal damage than demyelination . In a posterior study also in mice, Sun et al. established a relationship between reduced AD values and axonal damage, and between increased RD values and myelin damage, confirmed with immunohistochemistry examinations . In a human study, Pierpaoli et al. found decreased AD in primary lesions and in regions with secondary white matter degeneration . This axonal loss hypothesis has also been exposed for Alzheimer’s disease  and migraine patients .
Previous whole-brain TBSS studies found white matter differences between migraine patients and healthy controls [6,7,8,9,10,11, 52]. Most of these studies report decreased FA in migraine patients with respect to healthy controls, but one study showed increased FA in migraineurs with respect to healthy controls; in this study, Messina et al. analysed paediatric patients , which could explain the difference. Decreased FA can be caused by factors like demyelination, lower packing density or different membrane permeability  and is modulated by characteristics such as axon diameter and packing or fibre organization .
Results for MD and RD in the literature regarding migraine using TBSS are unclear. On the one hand, increased MD and RD values in migraineurs with respect to healthy controls were obtained in [11, 52], but on the other hand, decreased values were reported in [7,8,9]. Results obtained with methods different than TBSS, such as ROI-based analysis or tractography, showed increased MD and/or RD values in migraine patients with respect to healthy controls [13, 14, 18, 55]. Increased RD could also be a biomarker of demyelination [47, 48]. There are also studies that found no differences between migraine patients and healthy controls. These studies employed methods whole-brain such as TBSS [22, 56], ROI-based analysis (ictal migraine)  and voxel-based whole brain comparison .
Interestingly, the trend when comparing AD values between EM and healthy controls is inverted when comparing CM and healthy controls or CM with respect to EM (Additional file 1: Figures S4 and S5). As previously mentioned, AD might be an indicator of axonal loss. This might indicate that the evolution from EM to CM is characterised by a loss of axonal integrity. This result also shows that, in migraine, there could be different processes of axonal behaviour involving different pathophysiological mechanisms.
In contrast to our results, some TBSS studies obtained decreased AD in migraine patients with respect to healthy controls [7,8,9]. Petrušić et al. also found decreased AD in migraine with aura patients with respect to healthy controls using a tractography approach . We detected the same trend, but only in CM patients. We obtained the opposite result in EM patients, when including duration of migraine as a covariate, but this result was not significant in the post-hoc analysis. In [7, 8], cohort differences (considering only the EM patients in our case), i.e., lower disease duration and higher attack frequency, with possibly high frequency EM patients in the sample of [7, 8], could explain the differences between the studies. In , as previously stated, paediatric patients were included in the sample by Messina et al., which could explain the difference.
The analysis of the temporal change in migraine patients adds an interesting insight to the former results.
In the case of EM patients, considering the duration of migraine as a covariate, there was a lower number of regions with significant differences with respect to CM patients. The significant difference in the duration of migraine between both groups of migraine patients could be a confounding factor in the previous results. Furthermore, significant increased AD values in EM compared to HC were observed, possibly due to more precise estimations (less variability) in the case of the values in EM patients.
These results seem to indicate a temporal evolution in CM patients that reflects an adaptation to continuous headache attacks. In the initial months with CM (“short-term” patients), axonal integrity seems to be damaged, as suggested by decreased FA values in those patients compared to HC when including time from onset of CM as a covariate. Decreased FA values in CM patients with therapy compared to HC, after 6 months follow-up, were reported previously by Gomez-Beldarrain et al. . Considering these results and the decreased AD values in CM with respect to EM, it seems that in progression from EM to CM there might be a process that causes severe white matter alterations.
Later, in CM patients, a set of plastic changes as an adaptation to the frequent headaches may happen. The white matter reorganisation is suggested by positive correlation between FA and time from onset of CM, and by the simultaneous negative correlation between RD and time from onset of CM.
In line with our CM correlation results, Szabó et al., 2017 , obtained a trend showing increased FA in migraine with aura patients with respect to healthy controls. In this study, based on increased FA values, the authors hypothesised that repeated painful conditions or increased cortical excitability might cause maladaptive plastic changes in migraine with aura. Moreover, increased FA values were found in people with repeated stimuli in learning processes [59, 60], so something similar could be happening in CM patients, who suffer repeated painful stimuli.
Regarding correlation analysis, we obtained significant correlations between time from onset of chronic migraine and DTI measures (FA and RD) in the bilateral external capsule, but no significant correlations were obtained with the duration of migraine in CM or EM patients.
The external capsule is a part of the central core, a network on top of the brainstem that includes structures like the insular surface, the extreme and internal capsules, the lentiform nucleus or the thalamus . The extreme and external capsules lie in anteroposterior disposition, and they are connected to the anteroinferior part of the insula . The lentiform nucleus is located between the external and internal capsules . The internal, external and extreme capsules connect the insular surface, basal ganglia and thalamus to the cerebral lobes . Russo et al. reported that, in the insula and lentiform nucleus, migraine patients, compared to HC, are characterised by an increased blood oxygenation level dependent response . Furthermore, Borsook et al. exposed that the insula is implicated in processes related to the clinical presentation of migraine and is a “hub of activity” in migraine . The role of the external capsule in anteroposterior connections and in connections between subcortical regions implied in migraine and cortical regions could be highly relevant in migraine pathogenesis, especially in Chronic Migraine.
In other studies, significant negative correlations between duration of migraine in years (in EM patients) and FA [5, 7], MD  or AD [7, 10] were obtained, but also significant positive correlation between duration of migraine and MD . This discrepancy could be explained with methodological and cohort differences. No significant correlation between DTI measures and the external capsule has been previously found in migraine patients, but no correlation between time from onset of chronic migraine and DTI measures has been previously assessed.
In this study, high frequency EM patients (10–14 headache days per month) were excluded. This decision was made in order to avoid misclassified patients, which could mislead the analysis . Compared to Neeb et al. , CM patients from our sample had greater headache frequency, while EM patients from our sample had lower headache frequency. This increased difference between the migraine groups, together with the larger cohort size, could be a factor explaining why we obtained significant results in the CM-EM comparison, while no significant results were obtained in . However, no significant correlations were found between headache frequency and DTI measures, which could mean that headache frequency does not have a very relevant effect on diffusion within the EM or CM groups. In any case, a deeper specific analysis, focusing on high frequency EM patients, would be needed to clarify whether this group of patients is closer to the low frequency EM group or to the CM group.
There are several strengths and limitations in this study. About the strengths, this study is, to the best of our knowledge, the white matter study with the highest number of participants simultaneously including Healthy Controls, Episodic Migraine and Chronic Migraine patients. Moreover, the selection criteria of the patients allowed us to detect significant differences not found previously, especially between EM and CM patients.
About the limitations of this study, due to time constraints in the MRI acquisition process in a clinical setting, we acquired no T2 or T2-FLAIR MRI sequences that would be helpful to assess White Matter Hyperintensities (WMHs). Migraine has been associated with an increased risk for WMHs detected on MRI ; also pain in EM patients  and an unfavourable prognosis  were found to be associated with the occurrence of WMHs. Considering our correlation results in CM patients and the state of the art, the WMHs analysis would have been interesting in this study. Medication overuse was identified in an important percentage of the CM patients in our sample (75%). This might be a confounding factor, due to possible structural differences in the white matter with respect to CM patients without overuse. The exclusion of patients with anxiety or depression implies that there was no chance to assess possible effects of these conditions on brain structure in migraine patients. Anxiety and depression are often comorbid in patients with migraine [67,68,69]. When MRI were acquired in the patients, they had no attacks in the previous 24 h, but they could be in a prodromal stage, as we did control time from past, but not to the next migraine attack. Altered brain physiology and abnormal functional connectivity have been found in prodromal stages [70, 71], so this is a possible source of bias in the results. Diagnosis of infrequent TTH in controls was done solely by history and not by using a headache diary; however, they were excluded if other headache disorders were present or the frequency of headache in the preceding year was > 1 headache day per month or > 12 headache days per year. Finally, in the analysis of the presence of aura, the number of patients with migraine with aura was too small to additionally compare the changes in migraine with aura against migraine without aura.
In summary, considering previous studies and our results, a hypothesis about the migraine process could be drafted that distinguishes three states or stages: EM, transition from EM to CM, and CM. In the EM state, there would be some white matter damage, produced mainly by a loss of axonal integrity. Then, in the transition from EM to CM, there would be a loss of axonal integrity, but probably not led by a severe additional damage in myelin. Finally, in the CM state, there would be a series of plastic changes, as an adaptation to a continuous ictal state. Considering this evolution hypothesis, EM might involve a coexistence of loss of white matter integrity and maladaptive plasticity, with more severe integrity damage in the transition to CM, and CM may show predominant maladaptive plasticity, being, in this regard, a different entity with respect to a more frequent EM. An illustration of the hypothesised temporal evolution of the three main DTI measures used here can be seen in Fig. 3.
Anyway, regarding this hypothesis, longitudinal analysis should be performed in order to confirm this possible evolution, and the interpretation of the DTI measures must be carried out cautiously, as mentioned before in this section. Additionally, about the transition from EM to CM, the high frequency EM patients could be especially interesting to investigate. Indeed, it remains to be elucidated whether high frequency EM is actually a transition phase between EM and CM, an intense EM or, considering an extreme case, a low frequency CM.
Although this can be only considered a preliminary hypothesis, it would contribute to explain the high variability in the results in the literature found when comparing migraineurs with healthy controls in terms of white matter diffusion parameters. Indeed, if different trends can be found at different stages of the disease, then the results of global comparisons would depend heavily on the internal composition of the cohort of migraine patients in each study, possibly yielding what seem to be opposite results.
The current findings suggest global white matter structural differences between Episodic Migraine and Chronic Migraine, with damaged axonal integrity in Chronic Migraine, and between both groups of migraine patients and Healthy Controls. A different temporal pathophysiological evolution with maladaptive plastic changes seems to happen in Chronic Migraine with respect to Episodic Migraine. Further research is needed for these findings to be confirmed. Also, additional clinical features should be considered, longitudinal evaluation should be performed and a possible relationship with functional changes should be assessed.
Availability of data and materials
All patients who participated in the study signed a written informed consent. This consent did not include a statement to make the individual confidential data accessible to the public. Therefore, anonymized data are available from the corresponding author on reasonable request.
Analysis of covariance
Brain extraction tool
Diffusion magnetic resonance imaging
Diffusion tensor imaging
False discovery rate
Functional magnetic resonance imaging
Family wise error
International classification of headache disorders third edition
Montreal neurological institute
Magnetic Resonance Imaging
Region of interest
Tract-based spatial statistics
Threshold-free cluster enhancement
Turbo field echo
Headache Classification Committee of the International Headache Society (2018) The international classification of headache disorders, 3rd edition. Cephalalgia 38(1):1–211
Lipton RB, Fanning KM, Buse DC, Martin VT, Reed ML, Manack Adams A et al (2018) Identifying natural subgroups of migraine based on comorbidity and concomitant condition profiles: results of the chronic migraine epidemiology and outcomes (CaMEO) study. Headache 58(7):933–947
Katsarava Z, Buse DC, Manack AN, Lipton RB (2012) Defining the differences between episodic migraine and chronic migraine. Curr Pain Headache Rep 16(1):86–92
Aurora SK, Brin MF (2017) Chronic migraine: an update on physiology, imaging, and the mechanism of action of two available pharmacologic therapies. Headache 57(1):109–125
Li XL, Fang YN, Gao QC, Lin EJ, Hu SH, Ren L et al (2011) A diffusion tensor magnetic resonance image study of Corpus callosum from adult patients with migraine complicated with depressive/anxious disorder. Headache 51(2):237–245
Yuan K, Qin W, Liu P, Zhao L, Yu D, Zhao L et al (2012) Reduced fractional anisotropy of Corpus callosum modulates inter-hemispheric resting state functional connectivity in migraine patients without Aura. PLoS One 7(9):e45476
Yu D, Yuan K, Qin W, Zhao L, Dong M, Liu P et al (2013) Axonal loss of white matter in migraine without aura: a tract-based spatial statistics study. Cephalalgia. 33(1):34–42
Yu D, Yuan K, Zhao L, Dong M, Liu P, Yang X et al (2013) White matter integrity affected by depressive symptoms in migraine without aura: a tract-based spatial statistics study. NMR Biomed 26(9):1103–1112
Messina R, Rocca MA, Colombo B, Pagani E, Falini A, Comi G et al (2015) White matter microstructure abnormalities in pediatric migraine patients. Cephalalgia. 35(14):1278–1286
Szabó N, Faragó P, Király A, Veréb A, Csete G, Tóth E et al (2017) Evidence for plastic processes in migraine with Aura: a diffusion weighted MRI study. Front Neuroanat 11:138
Shibata Y, Ishiyama S, Matsushita A (2018) White matter diffusion abnormalities in migraine and medication overuse headache: A 1.5-T tract-based spatial statistics study. Clin Neurol Neurosurg 174:167–173
DaSilva AFM, Granziera C, Tuch DS, Snyder J, Vincent M, Hadjikhani N (2007) Interictal alterations of the trigeminal somatosensory pathway and PAG in migraine. Neuroreport 18(4):301–305
Coppola G, Tinelli E, Lepre C, Iacovelli E, Di Lorenzo C, Di Lorenzo G et al (2014) Dynamic changes in thalamic microstructure of migraine without aura patients: a diffusion tensor magnetic resonance imaging study. Eur J Neurol 21(2):287–e13
Chong CD, Schwedt TJ (2015) Migraine affects white-matter tract integrity: a diffusion-tensor imaging study. Cephalalgia. 35(13):1162–1171
Petrušić I, Daković M, Kačar K, Mićić O, Zidverc-Trajković J (2018) Migraine with aura and white matter tract changes. Acta Neurol Belg 118:485–491
Gomez-Beldarrain M, Oroz I, Garcia Zapirain B, Fernandez Ruanova B, Garcia Fernandez Y, Cabrera A et al (2015) Right fronto-insular white matter tracts link cognitive reserve and pain in migraine patients [erratum in J headache pain. 2016;17:22]. J Headache Pain 17:4
Liu J, Ma S, Mu J, Chen T, Xu Q, Dun W 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(10):5250–5259
Rocca MA, Pagani E, Colombo B, Tortorella P, Falini A, Comi G et al (2008) Selective diffusion changes of the visual pathways in patients with migraine: a 3-T Tractography study. Cephalalgia. 28(10):1061–1068
Granziera C, DaSilva AFM, Snyder J, Tuch DS, Hadjikhani N (2006) Anatomical alterations of the visual motion processing network in migraine with and without Aura. PLoS Med 3(10):e402
Pietrobon D, Moskowitz MA (2013) Pathophysiology of migraine. Annu Rev Physiol 75:365–391
Dodick DW (2018) A phase-by-phase review of migraine pathophysiology. Headache 58(Suppl 1):4–16
Neeb L, Bastian K, Villringer K, Gits HC, Israel H, Reuter U et al (2015) No microstructural white matter alterations in chronic and episodic Migraineurs: a case-control diffusion tensor magnetic resonance imaging study. Headache 55(2):241–251
Valfrè W, Rainero I, Bergui M, Pinessi L (2008) Voxel-based morphometry reveals gray matter abnormalities in migraine. Headache. 48(1):109–117
Neeb L, Bastian K, Villringer K, Israel H, Reuter U, Fiebach JB (2017) Structural gray matter alterations in chronic migraine: implications for a progressive disease? Headache. 57(3):400–416
Schulte LH, Allers A, May A (2017) Hypothalamus as a mediator of chronic migraine. Neurology 88(21):2011–2016
Smith SM, Jenkinson M, Johansen-Berg H, Rueckert D, Nichols TE, Mackay CE et al (2006) Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. Neuroimage. 31(4):1487–1505
Headache Classification Committee of the International Headache Society (2013) The international classification of headache disorders, 3rd edition (beta version). Cephalalgia 33(9):629–808
Serrano D, Lipton RB, Scher AI, Reed ML, Stewart WF, Manack Adams A et al (2017) Fluctuations in episodic and chronic migraine status over the course of 1 year: implications for diagnosis, treatment and clinical trial design. J Headache Pain 18(1):101
Zigmond AS, Snaith RP (1983) The hospital anxiety and depression scale. Acta Psychiatr Scand 67(6):361–370
Pelletier A, Periot O, Dilharreguy B, Hiba B, Bordessoules M, Chanraud S et al (2015) Age-related modifications of diffusion tensor imaging parameters and white matter Hyperintensities as inter-dependent processes. Front Aging Neurosci 7:255
Veraart J, Novikov DS, Christiaens D, Ades-Aron B, Sijbers J, Fieremans E (2016) Denoising of diffusion MRI using random matrix theory. Neuroimage. 142:394–406
Tournier J-D, Smith R, Raffelt D, Tabbara R, Dhollander T, Pietsch M et al (2019) MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation. bioRxiv Available from: https://www.biorxiv.org/content/biorxiv/early/2019/02/15/551739.full.pdf. Accessed 25 July 2019
Andersson JL, Sotiropoulos SN (2015) An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging. Neuroimage. 125:1063–1078
Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H et al (2004) Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23:S208–S219
Zhang Y, Brady M, Smith S (2001) Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans Med Imaging 20:45–57
Dhollander T, Raffelt D, Connelly A (2016) Unsupervised 3-tissue response function estimation from single-shell or multi-shell diffusion MR data without a co-registered T1 image. ISMRM Work Break Barriers Diffus MRI 5
Jenkinson M, Beckmann CF, Behrens TE, Woolrich MW, Smith SM (2012) FSL. Neuroimage 62:782–790
Rueckert D, Sonoda LI, Hayes C, Hill DL, Leach MO, Hawkes DJ (1999) Nonrigid registration using free-form deformations: application to breast MR images. IEEE Trans Med Imaging 18(8):712–721
Mori S, Wakana S, Nagae-Poetscher LM, van Zijl PC (2005) MRI atlas of human white matter. Elsevier, Amsterdam
Oishi K, Zilles K, Amunts K, Faria A, Jiang H, Li X et al (2008) Human brain white matter atlas: identification and assignment of commom anatomical structures in superficial white matter. Neuroimage 43(3):447–457
Hua K, Zhang J, Wakana S, Jiang H, Li X, Reich DS et al (2008) Tract probability maps in stereotaxic spaces: analyses of white matter anatomy and tract-specific quantification. Neuroimage 39(1):336–347
Nichols TE, Holmes AP (2002) Nonparametric permutation tests for functional neuroimaging: a primer with examples. Hum Brain Mapp 15(1):1–25
Smith SM, Nichols TE (2009) Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. Neuroimage. 44(1):83–98
Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B 57(1):289–300
Cohen J (1988) Statistical power analysis for the behavioral sciences, 2nd edn. Lawrence Erlbaum Associates, Publishers, Hillsdale
Bennett IJ, Madden DJ, Vaidya CJ, Howard DV, Howard JHJ (2010) Age-related differences in multiple measures of white matter integrity: a diffusion tensor imaging study of healthy aging. Hum Brain Mapp 31(3):378–390
Winklewski PJ, Sabisz A, Naumczyk P, Jodzio K, Szurowska E, Szarmach A (2018) Understanding the physiopathology behind axial and radial diffusivity changes — what do we know? Front Neurol 9:92
Alexander AL, Lee JE, Lazar M, Field AS (2007) Diffusion tensor imaging of the brain. Neurotherapeutics 4(3):316–329
Sun S-W, Liang H-F, Cross AH, Song S-K (2008) Evolving Wallerian degeneration after transient retinal ischemia in mice characterized by diffusion tensor imaging. Neuroimage. 40(1):1–10
Pierpaoli C, Barnett A, Pajevic C, Chen R, Penix LR, Virta A et al (2001) Water diffusion changes in Wallerian degeneration and their dependence on white matter architecture. Neuroimage 13(6 Pt 1):1174–1185
Shu N, Wang Z, Qi Z, Li K, He Y (2011) Multiple diffusion indices reveals white matter degeneration in Alzheimer’s disease and mild cognitive impairment: a tract-based spatial statistics study. J Alzheimers Dis 26(Suppl 3):275–285
Szabó N, Kincses ZT, Párdutz A, Tajti J, Szok D, Tuka A et al (2012) White matter microstructural alterations in migraine: a diffusion-weighted MRI study. Pain 153(3):651–656
Kochunov P, Thompson PM, Lancaster JL, Bartzokis G, Smith S, Coyle T et al (2007) Relationship between white matter fractional anisotropy and other indices of cerebral health in normal aging: tract-based spatial statistics study of aging. Neuroimage. 35(2):478–487
Beaulieu C (2002) The basis of anisotropic water diffusion in the nervous system - a technical review. NMR Biomed 15(7–8):435–455
Ito K, Kudo M, Sasaki M, Saito A, Yamashita F, Harada T et al (2016) Detection of changes in the periaqueductal gray matter of patients with episodic migraine using quantitative diffusion kurtosis imaging: preliminary findings. Neuroradiology. 58(2):115–120
Tedeschi G, Russo A, Conte F, Corbo D, Caiazzo G, Giordano A et al (2016) Increased interictal visual network connectivity in patients with migraine with aura. Cephalalgia. 36(2):139–147
Coppola G, Di Renzo A, Tinelli E, Di Lorenzo C, Di Lorenzo G, Parisi V et al (2016) Thalamo-cortical network activity during spontaneous migraine attacks. Neurology 87(20):2154–2160
Zhang J, Wu YL, Su J, Yao Q, Wang M, Li GF et al (2017) Assessment of gray and white matter structural alterations in migraineurs without aura. J Headache Pain 18(1):74
Blumenfeld-Katzir T, Pasternak O, Dagan M, Assaf Y (2011) Diffusion MRI of structural brain plasticity induced by a learning and memory task. PLoS One 6(6):e20678
Sampaio-Baptista C, Khrapitchev AA, Foxley S, Schlagheck T, Scholz J, Jbabdi S et al (2013) Motor skill learning induces changes in white matter microstructure and myelination. J Neurosci 33(50):19499–19503
Ribas EC, Yağmurlu K, de Oliveira E, Ribas GC, Rhoton A (2018) Microsurgical anatomy of the central core of the brain. J Neurosurg 129(3):752–769
Russo A, Coppola G, Pierelli F, Parisi V, Silvestro M, Tessitore A et al (2018) Pain perception and migraine. Front Neurol 9:576
Borsook D, Veggeberg R, Erpelding N, Borra N, Linnman C, Burstein R et al (2016) The insula: a “hub of activity” in migraine. Neuroscientist 22(6):632–652
Porter A, Gladstone JP, Dodick DW (2005) Migraine and white matter hyperintensities. Curr Pain Headache Rep 9(4):289–293
Yalcin A, Ceylan M, Bayraktutan OF, Akkurt A (2018) Episodic migraine and white matter Hyperintensities: Association of Pain Lateralization. Pain Med 19(10):2051–2057
Xie H, Zhang Q, Huo K, Liu R, Jian Z-J, Bian Y-T et al (2018) Association of white matter hyperintensities with migraine features and prognosis. BMC Neurol 18:93
Lantéri-Minet M, Radat F, Chautard MH, Lucas C (2005) Anxiety and depression associated with migraine: influence on migraine subjects’ disability and quality of life, and acute migraine management. Pain 118(3):319–326
Baskin SM, Lipchik GL, Smitherman TA (2006) Mood and anxiety disorders in chronic headache. Headache 46(Suppl 3):S76–S87
Minen MT (2016) Begasse De Dhaem O, kroon Van Diest a, powers S, Schwedt TJ, Lipton R, et al. migraine and its psychiatric comorbidities. J Neurol Neurosurg Psychiatry 87(7):741–749
May A (2017) Understanding migraine as a cycling brain syndrome: reviewing the evidence from functional imaging. Neurol Sci 38(Suppl 1):125–130
Karsan N, Goadsby PJ (2018) Biological insights from the premonitory symptoms of migraine. Nat Rev Neurol 14(12):699–710
ÁPG was supported by Junta de Castilla y León (Spain) and the European Social Fund (ID: 376062, Base de Datos Nacional de Subvenciones).
This study has been partially supported by grant from Gerencia Regional de Salud CyL (GRS 943/A/14, GRS 1727/A/18). This work was also supported by Ministerio de Ciencia e Innovación of Spain with research grants RTI2018–094569-B-I00 and PRX18/00253 (Estancias de profesores e investigadores senior en centros extranjeros).
Ethics approval and consent to participate
The local Ethics Committee of Hospital Clínico Universitario de Valladolid approved the study (PI: 14–197). All participants read and signed a written consent form prior to their participation.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
White matter alterations in chronic migraine compared to episodic migraine patients. Table S1. White matter regions where decreased AD values were found in CM compared to EM (FWE-corrected). Table S2. White matter regions where decreased AD values were found in CM compared to EM considering only patients with migraine without aura (FWE-corrected). Table S3. White matter regions from the ICBM-DTI-81 White Matter Atlas for which decreased AD values were found in CM compared to EM considering the effect of duration of migraine history (FWE-corrected). Table S4. White matter regions from the ICBM-DTI-81 White Matter Atlas for which increased AD values were found in EM compared to HC considering the effect of duration of migraine history (FWE-corrected). Figure S2. White matter alterations in migraine including duration of migraine history as a covariate. Table S5. White matter regions from the ICBM-DTI-81 White Matter Atlas for which significant AD values were found in CM compared to EM considering the effect of time from onset of CM (FWE-corrected). Table S6. White matter regions from the ICBM-DTI-81 White Matter Atlas for which decreased FA values were found in CM compared to HC considering the effect of time from onset of CM (FWE-corrected). Figure S3. White matter alterations in CM including time from onset of chronic migraine as a covariate. Table S7. Cohen’s d skeleton AD values in regions with significant differences between EM and CM. Figure S4. Cohen’s AD bar plots of regions with FWE-corrected differences between EM and CM (part 1). Figure S5. Cohen’s AD bar plots of regions with FWE-corrected differences between EM and CM (part 2). Figure S6. Cohen’s d FA bar plots (part 1). Figure S7. Cohen’s d FA bar plots (part 2). Figure S8. Cohen’s d RD bar plots (part 1). Figure S9. Cohen’s d RD bar plots (part 2). Figure S10. Cohen’s d MD bar plots (part 1). Figure S11. Cohen’s d MD bar plots (part 2).
About this article
Cite this article
Planchuelo-Gómez, Á., García-Azorín, D., Guerrero, Á.L. et al. White matter changes in chronic and episodic migraine: a diffusion tensor imaging study. J Headache Pain 21, 1 (2020). https://doi.org/10.1186/s10194-019-1071-3