Altered Resting State Functional Activity and Microstructure of the White Matter in Migraine With Aura

Introduction: Brain structure and function were reported to be altered in migraine. Importantly our earlier results showed that white matter diffusion abnormalities and resting state functional activity were affected differently in the two subtypes of the disease, migraine with and without aura. Resting fluctuation of the BOLD signal in the white matter was reported recently. The question arising whether the white matter activity, that is strongly coupled with gray matter activity is also perturbed differentially in the two subtypes of the disease and if so, is it related to the microstructural alterations of the white matter.

Methods: Resting state fMRI, 60 directional DTI images and high-resolution T1 images were obtained from 51 migraine patients and 32 healthy volunteers. The images were pre-processed and the white matter was extracted. Independent component analysis was performed to obtain white matter functional networks. The differential expression of the white matter functional networks in the two subtypes of the disease was investigated with dual-regression approach. The Fourier spectrum of the resting fMRI fluctuations were compared between groups. Voxel-wise correlation was calculated between the resting state functional activity fluctuations and white matter microstructural measures.

Results: Three white matter networks were identified that were expressed differently in migraine with and without aura. Migraineurs with aura showed increased functional connectivity and amplitude of BOLD fluctuation. Fractional anisotropy and radial diffusivity showed strong correlation with the expression of the frontal white matter network in patients with aura.

Discussion: Our study is the first to describe changes in white matter resting state functional activity in migraine with aura, showing correlation with the underlying microstructure. Functional and structural differences between disease subtypes suggest at least partially different pathomechanism, which may necessitate handling of these subtypes as separate entities in further studies.

Introduction

Migraine is the most widely known primary headache disorder, which affects nearly 10% of the population (1). The clinical appearance of attack is similar in most cases, however 20% of the patients report focal neurological symptoms (aura) before headache (migraine with aura; MWA). Different theories claimed to explain the pathomechanism of MWA and migraine without aura (MWoA) and the two forms of the headache may have distinct triggers (2). Malfunction of the trigemino-vascular complex, dysfunctional brainstem nuclei and altered metabolism of several neuropeptides may all play a part in the initiation of the migraineous headache. Aura symptoms are supposed to be originated from cortical spreading depression (CSD), a spreading wave of depolarization throughout the cortex that is thought to cause the focal symptoms associated with aura (3). CSD that is more prevalent in MWA is reportedly associated with the imbalance of excitation/inhibition in the migrainous brain (4). The altered excitability in MWA patients can cause increased amplitude of visual evoked potentials (5, 6), and also altered motion perception, where they have to judge certain properties of stimuli against a noisy background (7). Several studies also report enhanced responsiveness of the cortex to visual stimuli in a functional MRI setting (8, 9). Functional imaging placed the differences between the two subtypes on new foundations. Such methods were able to detect CSD during the aura phase (3, 10), and also revealed altered responses to sensory (11) or painful stimuli (12–14).

Resting state fMRI showed conflicting results about the alterations of brain activity during rest in interictal period. While some studies presented increased connectivity (15), some of them found the connectivity to be decreased (16). The fact, that these studies used heterogeneous populations of both MWA and MWoA patients, might explain this discrepancy. Our earlier investigation showed that the amplitude of the resting BOLD fluctuation, presumably a signature of increased excitability, was increased in MWA but not in MWoA (17).

Apart from gray matter functional networks, previous studies suggested that the white matter responds to external stimuli in a fashion similar to the gray matter, and this activation can also be measured with fMRI (18–20). Such fluctuations in activity are also present during rest, organized into distinct, consistently identifiable white matter resting state fMRI networks that bear close resemblance to gray matter functional networks (21–23). Several studies report that these networks function differently in neurological or psychiatric diseases (24–26). Since white matter microstructure is affected in migraine (27) and also differs between disease subtypes (28), a couple of questions arise. Does the functional activity of the white matter differs in disease comparing to the healthy population? Is there any difference between the various forms of migraine and if there it is, is it related to the microstructural alterations?

Resting state fMRI usually focus on the analysis of the whole frequency range of the time courses. However, it might be prudent to analyse different frequency ranges separately, since resting state networks were shown to be organize across several frequency bands (29). Several publications report resting state network abnormalities throughout several frequency bands in neuropsychiatric diseases and investigate the amplitude of low frequency fluctuations (30–33).

In our study, we aimed to explore differences in white matter resting state functional activity between the two migraine subtypes (MWA and MWoA) and healthy individuals. Considering the previously reported abnormalities of white matter microstructure, we also investigated whether these microstructural alterations are related to functional measures of the white matter during rest in migraine.

Materials and Methods

Participants

We recruited 51 migraine patients in our study (18 MWA, 33 MwoA). The participants belonged to the same cohort, similar to our previous investigation (17), all of them were enlisted from the Headache Outpatients Clinic of the Department of Neurology, University of Szeged. Apart from migraine, participants had no other neurologic or psychiatric disorders. All patients suffered from episodic migraine. None of the patients reported migraineous attack at in the preceeding week of the scanning. Each participant filled out a questionnaire about the prevalence and traits of the headache, medication and other diseases. None of the patient reported diabetes, or any other major cardiovascular diseases. Only one patient with aura was on interval therapy (iprazochrome). Seven of the 33 MWoA patients received interval treatment (1 topiramate, 1 amytriptilline, 5 iprazochrome).

As control group, we recruited 32 healthy volunteers, age and sex matched to the patients. Demographic data is included in Table 1.

TABLE 1

Table 1. Demographic data of the participants.

The local ethics committee of the University of Szeged approved the study (authority No.: 87/2009). All participants gave written informed consent in accordance with the Declaration of Helsinki.

Image Acquisition

MR imaging was performed on a 1.5 T GE Signa Excite HDxt MRI scanner. For every participant, we obtained high-resolution T1 weighted images (3D IR-FSPGR: TR/TE/TI: 10.3/4.2/450 ms, flip angle: 15°, ASSET: 2, FOV: 25^*25 cm, matrix: 256^*256, slice thickness: 1 mm,), a resting state fMRI protocol with echo-planar imaging technique (TE: 40 ms, TR: 3,000 ms, matrix: 64^*64 cm, FOV: 30^*30 cm, slice thickness: 6 mm, flip angle: 90°, NEX: 1, ASSET: 2,0 Ph, Phases per Loc: 128, volumes: 200) and 60 directions diffusion-weighted images with 6 non-diffusion-weighted reference volume [TE: 93.8 ms; TR: 16.000 ms; matrix: 96 * 96; FOV: 23 * 23 cm; flip angle: 90°; in-plane resolution: 2.4 * 2.4 mm; slice thickness: 2.4 mm; b: 1,000 s/mm²; number of excitations (NEX): 2; array spatial sensitivity encoding technique (ASSET) factor = 2]. Head motion was attenuated with foam pads.

Data Processing

We performed data pre-processing and statistical analyses with FSL (FMRIB's Software Library) toolkits. The Fourier transformation and statistical correlation calculation was performed by Matlab.

Pre-processing

We pre-processed resting state fMRI volumes using the FEAT toolbox, which included the removal of the first two non-steady-state volumes, removal of non-brain tissue via BET (34), motion correction (MCFLIRT), high pass filtering (cut-off sigma 100 s) and spatial smoothing with 6 mm FWHM. Because the motion parameters could alternate the BOLD signal, absolute and relative displacements were compared between groups. There were no significant differences in head movement. None of the subject's absolute displacement was larger then 0.8 mm.

All pre-processed fMRI images were registered to their own T1 images and to a standard structural image (MNI152, 2 mm isovoxel) with linear and non-linear registration respectively. To reduce calculation burden, we resampled all fMRI images to 4 mm isovoxel resolution. The pre-processed, standard space registered images were masked with a standard space white matter mask, thresholded at 0.99 probability.

We corrected diffusion data for eddy currents and movement artifacts by 12 degrees of freedom affine linear registration to the first non-diffusion-weighted reference image. Diffusion tensors at each voxel were fitted by FSL's FDT algorithm. Non-brain parts were removed with the brain extraction tool. We calculated fractional anisotropy (FA), mean diffusivity (MD), diffusivity axial (AD) and perpendicular (RD) to the principal diffusion direction for the whole brain, and aligned all subjects' FA images into a common space, using linear and non-linear registration. To match the resolution of functional volumes, we resampled the images to 4 mm isovoxel resolution.

Extraction of White Matter Resting State Networks

We identified group level, spatially independent resting state fMRI networks using independent component analysis as implemented in FSL's MELODIC (35). We used whole study population (patients and controls together) for networks' identification. Standard space white matter images entered the analysis flow after voxel wise demeaning and variance normalization. The individual participants' 4D fMRI data were concatenated into one dataset and decomposed into a set of matrices describing spatial and temporal domains, so that they showed maximally non-Gaussian distribution (fast ICA algorithm). We set the number of independent components to 30 and thresholded component spatial maps at p < 0.5. We selected components for further evaluation based on the following criteria: (I) no activation outside the brain or in the gray matter, (II) temporal fluctuation and power spectrum characteristics common to resting state gray matter networks (36, 37), (III) a spatial pattern similar to the previously described white matter networks (21).

Comparison of Networks and Power Spectrum Calculation

We used a dual regression approach for comparison between resting state white matter networks by using components' spatial maps as spatial regressors against the original subject-wise fMRI data, thus obtaining subject-specific versions of network time courses. We then regressed these time courses against the original fMRI data, which yielded individual spatial layouts for each network. To assess group differences in white matter network activity, we employed a standard GLM approach, with the group membership coded in the model. We used a non-parametric permutation test (5,000 permutations) for statistical inference, assessing cluster significance with the threshold-free cluster enhancement technique (38), and correcting for multiple comparisons via family wise error (FWE) correction (39).

To decompose network activity in the frequency domain, the time courses created in the first step underwent a Fast Fourier transform, generating a power spectrum for all individuals, for each network. We compared power spectra between groups by a two-sample t-test.

Relationship Between Resting State Fluctuation and Diffusion Parameters

Voxel-wise correlation was calculated between the diffusion parameters and the group difference functional activity of the white matter. The analysis was restricted to the voxels showing significant group difference of functional activity at a threshold of p < 0.05. We used a non-parametric permutation test (5,000 permutations) for statistical inference, assessing cluster significance with the threshold-free cluster enhancement technique (38), and correcting for multiple comparisons via family wise error (FWE) correction (39).

Correlation Between MRI Parameters and Clinical Data

Correlation was calculated between the diffusion parameters, functional activity of the white matter and clinical parameters. As dependent variable visual analog pain score (VAS), disease duration and attack frequency were used. MWA and MWoA groups were treated separately in our analysis. We used non-parametric permutation test for analysis and correcting for multiple comparisons via FWE.

Results

Demographical Data

Age and gender showed no significant differences between the three groups. The disease duration did not differ between the two migraine subgroups. MwoA patients showed significantly higher VAS of headache (p < 0.01) and attack frequency (p < 0.03) compared to MWA group.

White Matter Resting State Functional Networks

Eight white matter networks were identified, according to our inclusion criteria and publication revealed previously (21). The networks' spatial distribution are shown in the Supplementary Figure S1. Three white matter networks showed significant difference between patients and controls:

(I) IC2 included the genu (rostrum) of the corpus callosum, fiber bundles extending bilaterally into the frontal poles.

(II) IC7 contained bundles in the body of the corpus callosum, fibers extending bilaterally into the frontal and parietal lobes.

(III) IC17 contained the bilateral occipital white matter.

Different Expression of White Matter Functional Networks

MWA Compared to MWoA

The dual regression analysis identified higher expression in MWA compared to MWoA in all three networks.

For the network representing the genu of the corpus callosum (IC2), we found differences between MWA and MWoA in the right hemisphere and also in midline (Figure 1A).

FIGURE 1

Figure 1. Three white matter functional networks (IC2, A; IC7, B and IC17, C) had different expression in MWA and MWoA. Higher expression and amplitude was found in the case of MWA. The white matter functional networks are depicted in red-to-yellow thresholded to 0.5 probability. The group differences between MWA and MWoA is shown in blue-to-light blue thresholded at p < 0.05. The colorbar represents p-values. The group average Fourier spectrum presented next to the spatial maps. The red line represents MWA, the black MWoA. Green columns showed frequency values that different significantly between the two groups (p < 0.05).

For the network including the white matter fibers crossing in the body of the corpus callosum (IC7) differences appeared in the left frontal anterior area and right parietal white matter (Figure 1B).

For the network including the bilateral occipital bundles (IC17) differences were found next to the medial side of the occipital gray matter (Figure 1C).

All areas with connectivity differences showed increased amplitude of fluctuation in MWA compared to MWoA.

MWA Compared to Healthy Volunteers

IC2 showed changes in the bilateral frontal white matter (Figure 2A). IC17 revealed changes in the right occipital pathways (Figure 2B). In the area of significance, the amplitude of fluctuation was higher in MWA patients compared to healthy volunteers.

FIGURE 2

Figure 2. The two white matter functional networks showing differences between MWA and healthy volunteers (IC2, A and IC17, B). The white matter functional networks are depicted in red-to-yellow and thresholded to 0.5 probability. The group differences are shown in blue-to-light blue thresholded at p < 0.05. The colorbar represents p-values. The group average Fourier spectrum presented next to the spatial maps. The red line represents MWA, the black controls. Green columns showed frequency values that different significantly between the two groups (p < 0.05).

MWoA Compared to Healthy Volunteers

There were no significant differences between healthy volunteers and MWoA patients in any of the white matter resting state networks.

Connection Between Resting State Fluctuation and Underlying Microstructure

The resting state fMRI fluctuation of the network including the genu of the corpus callosum (IC2) was correlated with measures of the underlying microstructure. We found positive correlation with FA (r = 0.61, p < 0.05) and negative correlation with RD (r = −0.55; p < 0.05) in MWA patients (Figure 3).

FIGURE 3

Figure 3. Positive correlation was found between fractional anisotropy and frontal white matter activity fluctuations (A) (p < 0.05, corrected for multiple correlations) in MWA patients. Negative correlation was found between radial diffusivity and frontal white matter activity in the same area (B) (p < 0.05). The x axis represents the diffusion values of FA (A) and RD (B). The y axis shows the z-scores under the specified area.

No other functional network activity was associated with diffusion parameters.

Correlation Between MRI Parameters and Clinical Data

None of the resting state network, neither the diffusion parameters showed correlation with clinical parameters of the patients.

Discussion

Our investigation showed that white matter functional networks similar to the gray matter resting state functional networks show alterations in migraine with aura. Three white matter networks showed stronger expression, primarily higher amplitude of activity fluctuation in MWA. Furthermore, alterations of resting state functional fluctuation showed a strong association with parameters of the underlying white matter microstructure, a feature that has not been reported before.

Conventional fMRI measures the BOLD signal in the gray matter. However, several studies found task dependent activation in the white matter as well (22, 40). A combined DTI and fMRI study found functional activation in a white matter tract connected to the cortex (41). Moreover, similar to the gray matter, the effect of anesthesia and hypercapnia is detectable in the BOLD response of the white matter (42, 43). There are, however, significant differences that separate WM activation from GM activation: for example, the power spectrum of the WM response is reduced (22), and the hemodynamic response function (HRF) also differs from that in the gray matter, as the WM response delayed and has a lower amplitude (44, 45). Furthermore, the WM HRF was more variable in separate tracts under functional loading (46). Although these previous studies used specialized fMRI sequences, it is widely accepted that the measured fMRI activation strongly depends on the HRF model. The different response may be attributed to lower vascular density or diameter (47, 48), or the lack of post synaptic potentiation in the white matter (49). Furthermore, the gray matter's relatively greater vascular volume, higher cerebral blood flow, metabolism and connected hemodynamic response (49, 50) results in a higher BOLD signal (41) compared to that of the WM.

It was also presented that BOLD activity fluctuation can be measured in the WM too. Resting state fMRI fluctuation in the WM showed a temporal profile similar to the GM, with a comparable power spectrum. There is a detectable connection between GM and WM tract activation in rest (22, 23, 37). Symmetrical white matter resting state networks can be identified, which show strong overlap with anatomical white matter tracts (21).

Previous electrophysiological and MRI studies found functional alterations in MWA patients compared to MWoA and healthy volunteers. Visual evoked potential studies registered increased response in MWA patients (5, 6). A TMS meta-analysis also found lower phosphene threshold in MWA, but not in MWoA (51). The background of such altered excitability might be the alteration of local neurotransmitter milieu. Lower GABA level an indicator of reduced inhibition (52), and higher glutamate/glutamine ratio was found in the MWA (53).

As an indirect fMRI evidence of the altered excitability we found increased amplitude of resting BOLD fluctuation in large scale GM networks in MWA compared to MWoA (17, 52). Our analysis showed that this dysfunctions extends into the white matter, with increased connectivity and higher amplitude of BOLD fluctuation. In the background of the white matter functional alterations can stand from the connected axonal function or axonal loading transport changes, which are consequences of altered gray matter activity. The increased axonal usage generates higher energy demand that may affect the callosal much more (54). Higher activity in the occipital fibers might stem from the altered occipital neuronal activity connected to CSD phenomenon (our MWA group of patients reported mostly visual aura).

In our earlier investigations we found white matter microstructural differences in migraine (27), mostly pronounced in MWA (28). Most importantly in the current study we found connection between microstructural and functional alterations. White matter resting state activity correlated strongly with FA and RD values of the corresponding white matter pathway, which are often interpreted as measures of myelin content and microstructural integrity (55). A recent review discusses a similar connection between fractional anisotropy and fMRI fluctuation (56), however, information about structure-function relationship in the white matter is scarce. A possible explanation for this connection is that increased myelin concentration requires an increased energy supply, which could translate to chronic BOLD fluctuation changes.

The signs of altered cortical excitability perceived not just in the gray, but also in the white matter in MWA. Being either the cause of or consequence, structural alterations can be found in migraine too, that is more pronounced in MWA. These results indicate fundamental differences in the two subtypes of the disease that needs attention in clinical care and also during development of therapeutic interventions.

Of course, out investigation is not without limitations. While the unequal size of the groups and the relatively small sample size are partially accounted for by the use of non-parametric statistical approaches these factors could potentially limit the generalisability of our findings. While theoretically one should account for multiple comparisons when investigating multiple ICs or functional networks, correction for multiple comparisons are potentially over-conservative and not frequently used in the literature.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by the local ethics committee of the University of Szeged (authority No.: 87/2009). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

LV, ZK, NS, and PF planned the project and formulated the study hypothesis. ÁP, DS, and JT recruited the patients. ET, AK, BB, BK, DV, and KK organized and carried out the MRI measurements. BB, NS, and DV collected the clinical data. AK, PF, and DV analyzed the MRI data. DV, NS, ZK, PF, and JT formalized the discussion of the results. PF, ZK, NS, and AK wrote the manuscript.

Funding

The study was supported by: Neuroscience Research Group of the Hungarian Academy of Sciences and the University of Szeged, GINOP-2.3.2-15-2016-00034 grant, EFOP-3.6.1-16-2016-00008, NAP 2.0 (2017-1.2.1-NKP-2017-00002), National Brain Research Program (KTIA_13_NAP-A-II/20), UNKP-17-3 New National Excellence Program of the Ministry of Human Capacities, MTA-SZTE-Neuroscience Research Group and University of Szeged, Ministry of Human Capacities, Hungary grant (20391-3/2018/FEKUSSTRAT). BK is supported by the UNKP-18-3 New National Excellence Program of the Ministry of Human Capacities. Support was provided by the project CEITEC 2020 (LQ1601) from the Ministry of Education, Youth and Sports of the Czech Republic and from the National Center for Advancing Translational Sciences of the National Institutes of Health and European Union's Horizon 2020 research.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2019.01039/full#supplementary-material

Supplementary Figure S1. Spatial distibuion of the identified white matter networks.

References

1. Lipton RB, Stewart WF, Von Korff M. Burden of migraine: societal costs and therapeutic opportunities. Neurology. (1997) 48:S4–9. doi: 10.1212/WNL.48.3_Suppl_3.4S

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Pellegrino ABW, Davis-Martin RE, Houle TT, Turner DP, Smitherman TA. Perceived triggers of primary headache disorders: a meta-analysis. Cephalalgia. (2018) 38:1188–98. doi: 10.1177/0333102417727535

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Charles AC, Baca SM. Cortical spreading depression and migraine. Nat Rev Neurol. (2013) 9:637–44. doi: 10.1038/nrneurol.2013.192

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Pietrobon D, Moskowitz MA. Pathophysiology of migraine. Annu Rev Physiol. (2013) 75:365–91. doi: 10.1146/annurev-physiol-030212-183717

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Sand T, Zhitniy N, White LR, Stovner LJ. Visual evoked potential latency, amplitude and habituation in migraine: a longitudinal study. Clin Neurophysiol. (2008) 119:1020–7. doi: 10.1016/j.clinph.2008.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Coppola G, Bracaglia M, Di Lenola D, Di Lorenzo C, Serrao M, Parisi V, et al. Visual evoked potentials in subgroups of migraine with aura patients. J Headache Pain. (2015) 16:92. doi: 10.1186/s10194-015-0577-6

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Antal A, Temme J, Nitsche MA, Varga ET, Lang N, Paulus W. Altered motion perception in migraineurs: evidence for interictal cortical hyperexcitability. Cephalalgia. (2005) 25:788–94. doi: 10.1111/j.1468-2982.2005.00949.x

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Vincent M, Pedra E, Mourao-Miranda J, Bramati IE, Henrique AR, Moll J. Enhanced interictal responsiveness of the migraineous visual cortex to incongruent bar stimulation: a functional MRI visual activation study. Cephalalgia. (2003) 23:860–8. doi: 10.1046/j.1468-2982.2003.00609.x

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Datta R, Aguirre GK, Hu SY, Detre JA, Cucchiara B. Interictal cortical hyperresponsiveness in migraine is directly related to the presence of aura. Cephalalgia. (2013) 33:365–74. doi: 10.1177/0333102412474503

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Leao AA. Further observations on the spreading depression of activity in the cerebral cortex. J Neurophysiol. (1947) 10:409–14. doi: 10.1152/jn.1947.10.6.409

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Ren J, Xiang J, Chen Y, Li F, Wu T, Shi J. Abnormal functional connectivity under somatosensory stimulation in migraine: a multi-frequency magnetoencephalography study. J Headache Pain. (2019) 20:3. doi: 10.1186/s10194-019-0958-3

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Beese LC, Putzer D, Osada N, Evers S, Marziniak M. Contact heat evoked potentials and habituation measured interictally in migraineurs. J Headache Pain. (2015) 16:1. doi: 10.1186/1129-2377-16-1

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Russo A, Esposito F, Conte F, Fratello M, Caiazzo G, Marcuccio L, et al. Functional interictal changes of pain processing in migraine with ictal cutaneous allodynia. Cephalalgia. (2017) 37:305–14. doi: 10.1177/0333102416644969

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Kisler LB, Granovsky Y, Coghill RC, Sprecher E, Manor D, Yarnitsky D, et al. Do patients with interictal migraine modulate pain differently from healthy controls? A psychophysical and brain imaging study. Pain. (2018) 159:2667–77. doi: 10.1097/j.pain.0000000000001380

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Tedeschi G, Russo A, Conte F, Corbo D, Caiazzo G, Giordano A, et al. Increased interictal visual network connectivity in patients with migraine with aura. Cephalalgia. (2016) 36:139–47. doi: 10.1177/0333102415584360

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Niddam DM, Lai KL, Fuh JL, Chuang CYN, Chen WT, Wang SJ. Reduced functional connectivity between salience and visual networks in migraine with aura. Cephalalgia. (2016) 36:53–66. doi: 10.1177/0333102415583144

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Farago P, Tuka B, Toth E, Szabo N, Kiraly A, Csete G, et al. Interictal brain activity differs in migraine with and without aura: resting state fMRI study. J Headache Pain. (2017) 18:8. doi: 10.1186/s10194-016-0716-8

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Gawryluk JR, D'arcy RC, Mazerolle EL, Brewer KD, Beyea SD. Functional mapping in the corpus callosum: a 4T fMRI study of white matter. Neuroimage. (2011) 54:10–5. doi: 10.1016/j.neuroimage.2010.07.028

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Gawryluk JR, Mazerolle EL, D'arcy RC. Does functional MRI detect activation in white matter? A review of emerging evidence, issues, and future directions. Front Neurosci. (2014) 8:239. doi: 10.3389/fnins.2014.00239

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Marussich L, Lu KH, Wen HG, Liu ZM. Mapping white-matter functional organization at rest and during naturalistic visual perception. Neuroimage. (2017) 146:1128–41. doi: 10.1016/j.neuroimage.2016.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Peer M, Nitzan M, Bick AS, Levin N, Arzyt S. Evidence for functional networks within the human brain's white matter. J Neurosci. (2017) 37:6394–407. doi: 10.1523/JNEUROSCI.3872-16.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Ding ZH, Huang YL, Bailey SK, Gao YR, Cutting LE, Rogers BP, et al. Detection of synchronous brain activity in white matter tracts at rest and under functional loading. Proc Natl Acad Sci USA. (2018) 115:595–600. doi: 10.1073/pnas.1711567115

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Huang YL, Bailey SK, Wang PG, Cutting LE, Gore JC, Ding ZH. Voxel-wise detection of functional networks in white matter. Neuroimage. (2018) 183:544–52. doi: 10.1016/j.neuroimage.2018.08.049

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Zhang Z, Liao M, Yao ZJ, Hu B, Xie YW, Zheng WH, et al. Frequency-specific functional connectivity density as an effective biomarker for adolescent generalized anxiety disorder. Front Human Neurosci. (2017) 11:549. doi: 10.3389/fnhum.2017.00549

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Cocozza S, Pontillo G, Quarantelli M, Sacca F, Riccio E, Costabile T, et al. Default mode network modifications in Fabry disease: a resting-state fMRI study with structural correlations. Hum Brain Mapp. (2018) 39:1755–64. doi: 10.1002/hbm.23949

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Jiang Y, Luo C, Li X, Li Y, Yang H, Li J, et al. White-matter functional networks changes in patients with schizophrenia. Neuroimage. (2018) 190:172–81. doi: 10.1016/j.neuroimage.2018.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Szabo N, Kincses ZT, Pardutz A, Tajti J, Szok D, Tuka B, et al. White matter microstructural alterations in migraine: a diffusion-weighted MRI study. Pain. (2012) 153:651–6. doi: 10.1016/j.pain.2011.11.029

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Szabo N, Farago P, Kiraly A, Vereb D, Csete G, Toth E, et al. Evidence for plastic processes in migraine with aura: a diffusion weighted MRI study. Front Neuroanat. (2017) 11:138. doi: 10.3389/fnana.2017.00138

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Gohel SR, Biswal BB. Functional integration between brain regions at rest occurs in multiple-frequency bands. Brain Connect. (2015) 5:23–34. doi: 10.1089/brain.2013.0210

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Zhang Z, Liao W, Chen H, Mantini D, Ding JR, Xu Q, et al. Altered functional-structural coupling of large-scale brain networks in idiopathic generalized epilepsy. Brain. (2011) 134:2912–28. doi: 10.1093/brain/awr223

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Qi SL, Gao QJ, Shen J, Teng Y, Xie X, Sun YJ, et al. Multiple frequency bands analysis of large scale intrinsic brain networks and its application in schizotypal personality disorder. Front Comput Neurosci. (2018) 12:64. doi: 10.3389/fncom.2018.00064

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Shen YT, Li JY, Yuan YS, Wang XX, Wang M, Wang JW, et al. Disrupted amplitude of low-frequency fluctuations and causal connectivity in Parkinson's disease with apathy. Neurosci Lett. (2018) 683:75–81. doi: 10.1016/j.neulet.2018.06.043

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Xiao F, Wang T, Gao L, Fang J, Sun Z, Xu H, et al. Frequency-dependent changes of the resting BOLD signals predicts cognitive deficits in asymptomatic carotid artery stenosis. Front Neurosci. (2018) 12:416. doi: 10.3389/fnins.2018.00416

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Smith SM. Fast robust automated brain extraction. Human Brain Mapping. (2002) 17:143–55. doi: 10.1002/hbm.10062

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Beckmann CF, Deluca M, Devlin JT, Smith SM. Investigations into resting-state connectivity using independent component analysis. Philos Transac R Soc B Biol Sci. (2005) 360:1001–13. doi: 10.1098/rstb.2005.1634

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Niazy RK, Xie JY, Miller K, Beckmann CF, Smith SM. Spectral characteristics of resting state networks. Slow Brain Oscillations Sleep Resting State Vigilance. (2011) 193:259–76. doi: 10.1016/B978-0-444-53839-0.00017-X

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Ding ZH, Newton AT, Xu R, Anderson AW, Morgan VL, Gore JC. Spatio-temporal correlation tensors reveal functional structure in human brain. PLoS ONE. (2013) 8:82107. doi: 10.1371/journal.pone.0082107

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Smith SM, Nichols TE. Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. Neuroimage. (2009) 44:83–98. doi: 10.1016/j.neuroimage.2008.03.061

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Nichols TE, Holmes AP. Nonparametric permutation tests for functional neuroimaging: a primer with examples. Hum Brain Mapp. (2002) 15:1–25. doi: 10.1002/hbm.1058

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Gawryluk JR, Mazerolle EL, Beyea SD, D'arcy RC. Functional MRI activation in white matter during the Symbol Digit Modalities Test. Front Hum Neurosci. (2014) 8:589. doi: 10.3389/fnhum.2014.00589

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Mazerolle EL, Beyea SD, Gawryluk JR, Brewer KD, Bowen CV, D'arcy RCN. Confirming white matter fMRI activation in the corpus callosum: co-localization with DTI tractography. Neuroimage. (2010) 50:616–21. doi: 10.1016/j.neuroimage.2009.12.102

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Rostrup E, Law I, Blinkenberg M, Larsson HBW, Born AP, Holm S, et al. Regional differences in the CBF and BOLD responses to hypercapnia: a combined PET and fMRI study. Neuroimage. (2000) 11:87–97. doi: 10.1006/nimg.1999.0526

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Wu TL, Wang F, Anderson AW, Chen LM, Ding ZH, Gore JC. Effects of anesthesia on resting state BOLD signals in white matter of non-human primates. Magnetic Resonance Imag. (2016) 34:1235–41. doi: 10.1016/j.mri.2016.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Fraser LM, Stevens MT, Beyea SD, D'arcy RCN. White versus gray matter: fMRI hemodynamic responses show similar characteristics, but differ in peak amplitude. BMC Neurosci. (2012) 13:91. doi: 10.1186/1471-2202-13-91

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Courtemanche MJ, Sparrey CJ, Song XW, Mackay A, D'arcy RCN. Detecting white matter activity using conventional 3 Tesla fMRI: an evaluation of standard field strength and hemodynamic response function. Neuroimage. (2018) 169:145–50. doi: 10.1016/j.neuroimage.2017.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Li MW, Newton AT, Anderson AW, Ding ZH, Gore JC. Characterization of the hemodynamic response function in white matter tracts for event-related fMRI. Nat Commun. (2019) 10:1140. doi: 10.1038/s41467-019-09076-2

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Nonaka H, Akima M, Hatori T, Nagayama T, Zhang Z, Ihara F. Microvasculature of the human cerebral white matter: arteries of the deep white matter. Neuropathology. (2003) 23:111–8. doi: 10.1046/j.1440-1789.2003.00486.x

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Mozumder M, Pozo JM, Coelho S, Costantini M, Simpson J, Highley JR, et al. Quantitative histomorphometry of capillary microstructure in deep white matter. Neuroimage Clin. (2019) 23:101839. doi: 10.1016/j.nicl.2019.101839

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Logothetis NK, Pauls J, Augath M, Trinath T, Oeltermann A. Neurophysiological investigation of the basis of the fMRI signal. Nature. (2001) 412:150–7. doi: 10.1038/35084005

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Helenius J, Perkio J, Soinne L, Ostergaard L, Carano RA, Salonen O, et al. Cerebral hemodynamics in a healthy population measured by dynamic susceptibility contrast MR imaging. Acta Radiol. (2003) 44:538–46. doi: 10.1034/j.1600-0455.2003.00104.x

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Brigo F, Storti M, Tezzon F, Manganotti P, Nardone R. Primary visual cortex excitability in migraine: a systematic review with meta-analysis. Neurol Sci. (2013) 34:819–30. doi: 10.1007/s10072-012-1274-8

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Bridge H, Stagg CJ, Near J, Lau CI, Zisner A, Cader MZ. Altered neurochemical coupling in the occipital cortex in migraine with visual aura. Cephalalgia. (2015) 35:1025–30. doi: 10.1177/0333102414566860

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Gonzalez De La Aleja J, Ramos A, Mato-Abad V, Martinez-Salio A, Hernandez-Tamames JA, Molina JA, et al. Higher glutamate to glutamine ratios in occipital regions in women with migraine during the interictal state. Headache. (2013) 53:365–75. doi: 10.1111/head.12030

CrossRef Full Text | Google Scholar

54. Tracey I. Imaging pain. Br J Anaesth. (2008) 101:32–9. doi: 10.1093/bja/aen102

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Song SK, Yoshino J, Le TQ, Lin SJ, Sun SW, Cross AH, et al. Demyelination increases radial diffusivity in corpus callosum of mouse brain. Neuroimage. (2005) 26:132–40. doi: 10.1016/j.neuroimage.2005.01.028

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Warbrick T, Rosenberg J, Shah NJ. The relationship between BOLD fMRI response and the underlying white matter as measured by fractional anisotropy (FA): a systematic review. Neuroimage. (2017) 153:369–81. doi: 10.1016/j.neuroimage.2016.12.075

PubMed Abstract | CrossRef Full Text | Google Scholar