Dynamic Resting-State Functional Connectivity in Migraineurs
Department of Neurology, Japanese Red Cross Shizuoka Hospital, 8-2 Ohtemachi, Aoi-ku, Shizuoka, Japan
Academic Editor: Yasushi Shibata
Special Issue: The Pathophysiology and Treatment for Migraine
Received: July 19, 2022 | Accepted: October 18, 2022 | Published: October 26, 2022
OBM Neurobiology 2022, Volume 6, Issue 4, doi:10.21926/obm.neurobiol.2204143
Recommended citation: Imai N. Dynamic Resting-State Functional Connectivity in Migraineurs. OBM Neurobiology 2022;6(4):15; doi:10.21926/obm.neurobiol.2204143.
© 2022 by the authors. This is an open access article distributed under the conditions of the Creative Commons by Attribution License, which permits unrestricted use, distribution, and reproduction in any medium or format, provided the original work is correctly cited.
Migraine is a clinical syndrome characterized by headache with specific features and accompanying symptoms such as nausea, vomiting, photophobia, and phonophobia .
Functional MRI (fMRI) investigates the mechanisms lead to sensory hypersensitivity in migraine by measuring brain’s response to sensory stimuli, while resting-state functional MRI analysis (rs-fMRI) investigates the functional connectivity of specific brain regions and networks responsible for sensory processing .
The rs-fMRI analysis methods can be classified into model-dependent and model-free methods . Model-dependent methods examine functional connectivity by correlating the resting time series of a known region of interest (ROI) with the time series of all other regions, which clearly shows which regions’ ROI is functionally connected. The model-free methods aim to search for general patterns of unique functional connectivity. Independent component analysis (ICA) is the most popular model-free method. ICA is applicable to whole brain voxel unit data, easily selects the time signals of resting-state independent components (ICs), and shows a high degree of consistency.
Conventionally, rs-fMRI analysis implicitly assumes that the statistical interdependence of signals between different brain regions is constant during the recording of a task-free experiment [4,5]. However, by dividing the time series into short time segments, it is possible to detect time-varying patterns in network strength and components on the basis of sliding time window correlation [4,5]. This method is known as dynamic connectivity analysis, and the functional connectivity for obtaining dynamic connectivity analysis is called dynamic functional connectivity (d-FC). Contrary to dynamic functional connectivity, conventional functional connectivity is called static functional connectivity (s-FC). Dynamic connectivity analysis better reflects fluctuations in neuronal activity and helps identify networks that cannot be detected by static functional analysis [4,5]. We have described previous studies using dynamic connectivity analysis in migraine patients, including our previous study [6,7,8,9,10,11]. In this review, we provided insight into resting-state dynamic functional analysis in migraine patients.
2. Overview of Previous Dynamic rs-fMRI Studies
2.1 Characteristics of Participants
All 6 studies were conducted in migraine patients [6,7,8,9,10,11]. The diagnosis of migraine differed according to the purpose of the study and included migraine , migraine without aura (MwoA) [6,10,11], episodic migraine (EM) , and chronic migraine (CM) [9,11]. The number of participants, including controls, ranged from 37 to 159 [6,7,8,9,10,11].
2.2 Migraine Phase in MRI
Migraine is divided into five phases: prodromal, aura, ictal, postictal, and interictal phases . Of these, the prodromal, aura, and postictal phases are together referred to as the peri-ictal phase. MRI images were acquired during the interictal phase in three studies [6,10,11], and one study  assumed an interictal phase. Two studies [7,8] acquired images of all phases, one of which  classified the images into the preictal phase and ictal/peri-ictal phases.
2.3 The Methods of d-FC Analysis
All studies performed d-FC analyses using the sliding window method [6,7,8,9,10,11]. However, the analyzed data using the sliding window method differed from study to study: three studies used ICs obtained from group ICA (GICA) [6,8,9], two studies used predefined ROIs [7,11] and one study used the seed obtained from regional homogeneity (ReHo) and amplitude of low frequency fluctuations (ALFF) analyses . Of these, two tools were used in GICA; two studies [6,9] used the GIFT (https://trendscenter.org/software/gift/) [12,13] and one study  used the FSL MELODIC software [14,15,16,17]. For predefined ROIs: one was 59 ROIs related to pain processing  and the other was 132 ROIs for the whole brain .
2.4 Comparison of d-FC with s-FC
Dynamic connectivity analysis can detect more fluctuations in neuronal activity than static functional analysis [4,5]. Four of the six studies compared d-FC with s-FC to determine whether the advantage of dynamic connectivity analysis can be detected in migraine patients [7,8,9,11]. Dumkrieger et al. reported that the significantly different FC between migraine and persistent posttraumatic headache (PPTH) were 17 region pairs in s-FC and 10 region pairs in d-FC . Likewise, the ROI pairs with significant differences in migraine and PPTH in s-FC and d-FC were different. Lee et al. showed a single significant different network between interictal migraine patients and controls, whereas there were no significantly different networks between ictal/peri-ictal migraine patients and controls in static connectivity analysis . In dynamic connectivity analysis, however, there were seven significant different networks each between interictal and ictal/peri-ictal migraine patients and their respective controls. Also, four ICs had significantly different connectivity in both interictal and ictal/peri-ictal migraine patients compared to controls. In addition, Zuo et al. demonstrated that there were no significantly different networks in static functional network connectivity (s-FNC) while five significantly different networks in dynamic functional network connectivity (d-FNC) were indicated . Our previous study showed that there were 18 significantly different connectivity pairs between photophobic and non-photophobic patients, 15 significantly different connectivity pairs between phonophobic and non-phonophobic patients, and one significant different connectivity pair between patients with and without osmophobia in s-FC analysis . While in d-FC analysis, there were 16 significantly different connectivity pairs between photophobic and non-photophobic patients, 8 significantly different connectivity pairs between phonophobic and non-phonophobic patients, and 14 significantly different connectivity pairs between osmophobic and non-osmophobic patients. Thus, in both s-FC and d-FC analysis, there were significantly different connectivity pairs between patients with and without photophobia, while no common significant connectivity pairs were observed between patients with and without phonophobia and between patients with and without osmophobia. The results of the above 4 studies showed that dynamic connectivity analysis can detect more significant networks and different connectivity pairs than static connectivity analysis.
3. The Details of Each Study
3.1 Abnormal Thalamocortical Network Dynamics in Migraine
Functional and morphological studies have suggested that the thalamocortical pathway plays an important role in the pathogenesis of migraine [18,19,20]. Tu et al. in 2019 reported abnormal thalamocortical network dynamics in migraine . In this study, the authors investigated the d-FNC of thalamocortical networks in 89 interictal migraine patients and 70 healthy controls (HCs), and whether clinical features were associated with abnormal connectivity.
In detail, they performed 4 major steps within the framework of characterizing dynamic d-FNC to detect atypical functional dynamics in migraine. In step 1, GICA was conducted to decompose the whole brain resting-state fMRI data into multiple ICs using the GIFT toolbox [12,13]. Following GICA, intrinsic component networks (ICNs) were selected from the ICs according to spatial activation maps . In step 2, the d-FNC in ICNs was calculated using a sliding window approach with graphical lasso (least absolute shrinkage and selection operator), the graphical lasso was a remarkably fast simple algorithm that used a coordinate descent procedure for the lasso [4,5,22]. In step 3, the d-FNC estimates were hard clustered, i.e., the data items were grouped so that each item was assigned to only one cluster to assess the recurring d-FNC patterns over time . Next, whether MwoA and HCs presented different occurrences of different functional d-FNC states during rs-fMRI was investigated. A 2-sample t-test was performed to examine the group difference in occurrence between MwoA and HCs for each d-FNC state, thereby, the presence of abnormal transient d-FNC was further investigated. In step 4, the efficiency of information transfer in dynamic brain networks is demonstrated to be variable by applying graph theory measures on different d-FNC states [24,25,26,27,28,29]. Specifically, ICNs were defined as nodes and the d-FNC between them was defined as edges, and efficiency was inversely proportional to the harmonic mean of the shortest distance (number of edges) between all possible pairs of nodes. The global efficiency was the average efficiency across all node pairs, while the local efficiency was the average of the nodal local efficiency within neighbors of the node.
The results identified 52 ICNs which were categorized into 5 d-FNC brain states to characterize and compare dynamic functional connectivity patterns. Group differences of occurrences and dFNC patterns showed migraineurs had a significantly higher occurrence rate in state 1 with strong positive connectivity within the somatomotor domain and visual domain, and a significantly lower occurrence rate in state 2 with sparsely connected. In addition, abnormalities in d-FNC with visual cortex and the precuneus of the posterior thalamus (nucleus pulvinar) were significantly correlated with migraine frequency. The posterior nucleus of the thalamus receives projections from the brainstem and transmits them to the primary and secondary somatosensory cortices, insula, primary and secondary visual cortices, primary auditory cortex, and anterior cingulate cortex . It has been suggested that abnormal connections between the posterior nucleus of the thalamus and visual cortex may be implicated in clinical features of migraine, such as photophobia and allodynia. In addition, topological measurements revealed that migraine patients had significantly reduced information transfer efficiency in the global and local d-FNC.
3.2 Differences in Static and Dynamic Functional Connectivity Between Migraine and Persistent Posttraumatic Headache
Migraine and PPTH often share common phenotypic features. The main factor distinguishing PPTH from migraine is the head injury itself . The structure of migraine and PPTH differ in regions within the right lateral orbitofrontal lobe, left caudal middle frontal lobe, left superior frontal lobe, left precuneus and right supramarginal gyrus. . Dumkrieger et al. investigated the differences in brain function between migraine and PPTH in terms of s-FC and d-FC patterns in 59 preselected regions involved in pain processing .
The study explored functional connectivity using an ROI approach. Fifty-nine ROIs (29 bilateral regions and one midline area) were predetermined according to prior literature [32,33,34,35,36,37,38,39]. For each pair of ROIs, a sliding window correlation analysis was performed and the standard deviation of the obtained time course values was calculated. Resting-state data in 10 minutes was collected, and 60 seconds window lengths were selected to increase confidence in capturing time-varying fluctuations that can be missed using shorter window lengths. The d-FC values were compared between groups using ANOVA followed by t-tests adjusting for gender and age. And the Benjamini-Hochberg procedure was used to correct for multiple comparisons 
The results showed that significant differences in s-FC between migraine and PPTH were found in 17 pairs of regions and the major connections regions were the left ventral medial prefrontal cortex and left secondary somatosensory cortex. Significant differences in d-FC between migraine and PPTH were found in 10 region pairs and the major connections regions were the left secondary somatosensory cortex and left fusiform gyrus. There was overlap in the regions of difference in s- and d-FC between PPTH and migraine, but not in functional connectivity between region pairs. In the migraine group, s-FC between the left secondary somatosensory cortex and the right cuneiform cortex was significantly correlated with headache frequency, and d-FC between the right cingulate and right amygdala was significantly correlated with pain intensity. In the PPTH group, s-FC between the left middle cingulate and the right pulvinar, right posterior insula and hypothalamus was significantly correlated with headache frequency, and d-FC between the right middle cingulate and right supramarginal gyrus was significantly correlated with headache frequency. This study showed that s-FC and d-FC coupling of pain-processing and visual-processing regions was different in migraine and PPTH. Both static and dynamic functional analysis may be useful to study brain function in migraine and PPTH and perhaps to distinguish between PPTH and migraine.
3.3 Dynamic Functional Connectivity of the Migraine Brain During Interictal and Ictal/Peri-Ictal Phases
Migraine is a clinical syndrome characterized by headache with specific features and associated symptoms, such as nausea, vomiting, photophobia and phonophobia . The period when the headache attacks is called the ictal phase. Some patients experience an aura, which is a prodromal phase occurring 24 hours before the headache and a postdromal phase following headache attacks. The aura, prodromal and postdromal phases together are called peri-ictal phase. The period outside ictal phase and peri-ictal phase is called interictal phase. Lee et al. examined the differences in dynamic functional connectivity between migraine patients and the controls . 50 migraine patients and 50 age- and sex-matched control subjects were recruited. The migraine patients were divided into interictal phase and ictal/peri-ictal phase. All subjects underwent resting-state functional magnetic resonance imaging.
First, ICs were calculated using GICA with the FSL MELODIC program [14,15,16,17,41]. Second, the spatial map and temporal/spectral properties of each IC were taken into account [42,43]. The signal ICs were employed for an investigation of functional connectivity while being regarded as brain networks. Only rs-fMRI data from interictal individuals were used for group ICA, and the ictal/peri-ictal subjects were then given access to the constructed brain networks. Based on graph theory, static and dynamic connectivity evaluations were carried out [24,44]. The group ICA-defined brain networks were regarded as graph nodes. For static connectivity analysis, the Pearson correlation coefficient of the entire time series between various nodes was used to generate graph edges [24,45,46]. Dynamic connectivity analysis was carried out in addition to static connectivity analysis using sub-time series derived from sliding windows. Next, the subject-level regularized dynamic connectivity matrices were subjected to a k-means clustering technique to characterize brain states and decrease the dimensionality of the matrices [5,22,47]. The silhouette coefficient and elbow approach were used to calculate the ideal number of clusters [5,22,48]. The group-level brain states were defined using the most common number of clusters from the subject-level brain states. The average of the dynamic connectivity matrices in the identical brain states was used to determine the EC values for each individual. Finally, EC values were determined for each subject in various brain states. The calculated EC values were utilized to evaluate the differences between patients and controls within the interictal and ictal/peri-ictal groups.
This study defined significant networks from the interictal patients and matched controls in a data-driven manner and tested patients and controls in ictal or peri-ictal data sets. Static analysis revealed no significant networks in both the interictal and interictal/preictal/period data sets. While, dynamic analysis revealed significant group differences in 7 brain networks in the interictal data set (default mode network, frontoparietal network and cerebellum network; normal controls>patients, 3 brainstem networks and pain modulatory network; normal controls
3.4 Aberrant Modulations of Static Functional Connectivity and Dynamic Functional Network Connectivity in Chronic Migraine
Default mode network (DMN), salience network (SN), and central executive network (CEN) are known as representative intrinsic functional brain networks. Previous studies have shown a significant decrease in DMN, SN and CEN connectivity for all CM participants, regardless of MOH status, when compared to controls [49,50]. To confirm previous studies and evaluate d-FNC, Zuo et al. investigated the functional characteristics of brains with chronic migraine (CM) using s-FC, s-FNC, and d-FNC analyses .
13 ICs were identified from 17 CM patients and 20 gender- and age-matched HCs, which were classified into six static state networks.
The intrinsic connection networks were first identified using ICA by Group ICA Of fMRI Toolbox (GIFT) software . The minimal description length criteria estimated 20 ICs and visual inspection identified 13 ICs as meaningful. The 13 ICs then categorized into 6 resting-state networks, such as DMN. SN, CEN, auditory network (AN), visual network (VN) and cerebellum network. The s-FNC was evaluated using MANCOVAN toolbox of GIFT software. The d-FNC was analyzed using the sliding window method with the temporal d-FNC toolbox in GIFT software [4,52].
In s-FNC, there were no significant differences between CM and HCs. However, the d-FCN analysis between CM and HCs showed a significantly decreased connectivity between the DMN (IC11) and AN (IC 5) and 4 significant increased connectivity between ECN (IC2, IC4) and DMN (IC19), between the ECN (IC 4) and AN (IC 5), and between the ECN (IC 4) and VN (IC 13).
This study suggested that the functional connectivity of the CM brain was not static but may fluctuate over time. Significantly increased connectivity was found between the ECN and DMN, the ECN and AN, and the ECN and VN. The DMN and ECN appeared to be related with cognitive functions. AN and VN probably regulated the sensory system. This study provided further evidence that the chronicity of migraine may be related to abnormalities in connectivity between sensory and cognitive brain networks.
3.5 Disrupted Dynamic Functional Connectivity of the Visual Network in Episodic Patients with Migraine Without Aura
Migraine can be classified into migraine with aura (MwA) and MwoA, depending on the presence or absence of aura . The visual aura is the most common type of aura and occurs among more than 90 % of patients with MwA . Photophobia, like visual aura, is a common visual symptom of migraine. Photosensitivity is one of the most bothersome symptoms of MwoA as well as MwA. Photophobia is not only present during the ictal phase, but also during the interictal phase. Several recent studies have addressed the altered s-FC of the visual cortex in MwoA [32,54,55]. Wei et al. investigated the relationship between d-FC of the visual cortex and clinical features in patients with MwoA in interictal phase .
Fifty-five MwoA patients and 50 sex- and age-matched HCs were included in the study. Data preprocessing was as follows. First, commonly affected brain regions were chosen as seeds, and 5 mm peak MNI coordinates were used to build them. The Dynamic Brain Connectome Toolbox (V2.1 http://restfmri.net/forum/DynamicBC) was used to conduct the d-FC analysis [22,56]. A sliding window strategy was used to convolve the rectangles with a Gaussian kernel to create temporal dynamic patterns. 64 windows per person were created with a window size of 30 TRs and a window overlap of 90%. In each sliding window, the temporal correlation coefficient (r) between each seed and the remaining brain voxels over time were measured. Multiple sliding window correlation maps were obtained for each subject. To quantify temporal variation in functional connectivity, the coefficient of variation over time for each voxel in each window were computed and the Fisher Z transform was used to obtain variables that were close to a normal distribution.
Regional changes in visual cortex were assessed using ReHo and ALFF. ReHo is one of the indicators of regional synchronization of spontaneous brain activities, and ALFF is one of the markers of the regional intensity of spontaneous brain activities. Seed-based analysis was performed to evaluate d-FC between the visual cortex and the whole brain using the Dynamic Brain Connectome Toolbox (V2.1 http://restfmri.net/forum/DynamicBC).
The results showed that patients with MwoA had reduced ReHo values in the right lingual gyrus and reduced ALFF levels in the the right lingual gyrus and angular gyrus, and increased ALFF levels in the left superior frontal gyrus, right inferior frontal gyrus, and right middle frontal gyrus compared with HCs. To evaluate the d-FC between the visual cortex and the whole brain, the right lingual gyrus was selected as the seed region because of reduced ReHo values and ALFF levels. And the d-FC analysis showed asignificant decrease in the lingual gyrus in relation to the right calcarine sulcus, bilateral cuneus, right fusiform gyrus and bilateral postcentral gyrus, and a significant increase relation to the left thalamus, right insula, right parahippocampus, right hippocampus, right angular gyrus, bilateral middle cingulate cortex, bilateral posterior cingulate cortex and bilateral precuneus. Abnormalities in the d-FC of the right lingual gyrus and bilateral cuneus were positively correlated with anxiety scores evaluated using General Anxiety Disorder scores . Previous study showed Alzheimer's disease patients with depression have decreased FC between the right lingual gyrus and dorsal anterior cingulate cortex . Anxiety, like depression, is a risk factor for migraine chronicity . A recent study reported that anxiety was more robustly associated with increased risk of migraine than depression . This study suggested that d-FC abnormalities in the visual cortex were involved in multinetwork pain integration in episodic MwoA patients and may be associated with anxiety disorders.
3.6 Functional Connectivity in Migraineurs with Photophobia, Phonophobia, or Osmophobia
Finally, differences in s-FC and d-FC between patients with and without photophobia, phonophobia, or osmophobia were explored . Sixty-two migraine patients with/without photophobia, phonophobia, or osmophobia underwent fMRI during the interictal phase. ROI-to-ROI analysis in the whole brain was performed to compare s-FC and d-FC.
This study started with the d-FC analysis using the CONN toolbox [59,60,61,62] with dynamic ICA. The dynamic properties (temporal modulation) of the ROI-to-ROI connection matrix were investigated by dynamic connectivity analysis, which also identified some connectivity circuits with similar modulations. Next, various modulating circuits were depicted through the dynamic ICA matrix and the rate of change connectivity between each pair of ROIs was measured. The rate of connectivity change was determined by the strength and sign of connectivity change conjugated to a particular component/circuit time series. The simple generalized form of the context-dependent psychophysiological interaction (gPPI) model was used to estimate the group-level modulation component gamma _l(i,j). The fast ICA with hyperbolic tangent contrast functions was then used to rotate the group-level modulation components gamma _l(i,j), and the ICA mixing matrix W was inverted to compute the dynamic iIC/circuit time series. Finally, the estimated dynamic IC/circuit time series h(t) was back-projected to a set of group-level modulation components of subject-level components gamma _nk(i,j) using a variety of conventional first-level gPPI models with gPPI psychological variables. The ROI-to-ROI connectivity matrix was used to determine the correlations between the conditions for each participant.
Analysis of s-FC showed that 18 significant different connectivity pairs patients in photophobic patients, 15 significant different connectivity pairs in phonophobic patients, while connectivity between the right cerebellar lobe and the brainstem was significantly higher in patients with osmophobia (Figure 1 A, Figure 1B, Figure 1C, static).
Figure 1 Results of s-FC and d-FC using ROI-to-ROI analysis. The Figure in (A), (B), and (C) demonstrated significantly different s-FC and d-FC between migraineurs with or without photophobia, phonophobia, and osmophobia. Patients with photophobia, phonophobia, or osmophobia have much higher connections on the red line than those who do not. Patients with photophobia, phonophobia, or osmophobia have much lower connectivity on the blue line than those without these phobias.
According to the analysis d-FC, patients with photophobia had 16, patients with phonophobia had 8, and patients with osmophobia had 14 distinct significant connectivity pairs (Figure 1A, Figure 1B, Figure 1C, dynamic). There were no appreciable overlapping linkages between the s-FC and d-FC in phonophobic and osmophobic patients.
Our study revealed that patients with photophobia, phonophobia, or osmophobia had different s-FC and d-FC values. We hypothesized that these variations in s-FC and d-FC might be due to distinct brain dynamics and that s-FC and d-FC analyses could be useful for determining the pathophysiological mechanisms behind migraine.
Univariate and multivariate approaches commonly applied to resting-state data assume that the strength of interactions between regions is not changing over time [2,3,4,5]. However, recent theories have demonstrated that the brain is a highly interconnected and dynamic system [4,5]. Therefore, conventional functional coupling analysis has limited ability to reflect the dynamic processes of neural signals and may ignore potential changes. Resting-state d-FC and d-FNC analysis may be a useful neuroimaging method to overcome this static limitation and explore temporally changing patterns. As mentioned in Section 2.4, four of six studies compared d-FC with s-FC, suggesting that dynamic connectivity analysis can detect more significant networks and different connectivity pairs compared to static connectivity analysis [7,8,9,11]. For example, the study by Dumkrieger et al. showed significantly different ROI pairs between migraine and PPTH, migraine and controls, and PPTH and controls differed between s-FC and d-FC . In our previous study, comparing migraineurs with and without photophobia, 18 significantly different functional connectivity were found in s-FC and 16 in d-FC . These significantly different functional connectivities were found mainly between the cerebellum and temporal lobe regions. And significantly different functional connectivities were found in the cerebellum and temporal lobe, but not in the occipital lobe and the visual cortex, suggesting that photophobia may be associated not only with the occipital lobe but also with a wide range of regions, including the cerebellum and temporal lobe. Meanwhile, we found 15 significantly different functional connectivity in s-FC and 8 in d-FC comparing migraineurs with and without phonophobia, one significantly different functional connectivity in s-FC and 14 in d-FC comparing migraineurs with and without osmophobia . These results suggested that different hypersensitivity symptoms such as photophobia, phonophobia and osmophobia have different time windows in which significantly different functional connectivity may be detected. The results of two other studies that may have demonstrated the advantages of d-FC are as follows [6,10]: Previous studies have suggested that thalamocortical dysrhythmia and abnormal low-frequency oscillations in thalamocortical networks are associated with clinical symptoms of migraine and strongly influence important pain processes with combining multisensory [18,19,20]. The study by Tu et al. specified taht the posterior–pulvinar thalamic complex was an abnormal functional component of migraine . Anxiety, like depression, is a risk factor for migraine chronicity . Wei et al. suggested that abnormalities in d-FC of the visual cortex may be associated with anxiety disorders . In summary, dynamic connectivity analysis can detect more significant networks and different connectivity pairs than static connectivity analysis.
Schwedt et al. mentioned several limitations of conventional fMRI studies . First, some fMRI studies of migraine involved small numbers of patients, which limited statistical power and was often suboptimal methods for determining significance, limiting the generalizability of the results. In addition, there were few replicated studies confirming fMRI results. For example, several studies have suggested that MwA undergo s-FC changes compared to MwoA and HCs, not only in the visual cortex but also in a wide range of regions involved in visual processing. However, the study by Hougaard et al. reported that among 40 MwA patients with 40 age- and sex-matched HCs, there was no difference between MwA patients and HCs in these networks, or even a trend towards no difference . Therefore, the absence of s-FC difference in a large sample of MwA and HCs questioned the findings of previous studies and emphasized the general requirement for rs-fMRI results to be replicated before making solid conclusions [63,64]. These limitations also apply to the d-FC and d-FNC analysis in the presented studies. The results of these studies have not been followed up, and have not been evaluated for reproducibility. The study by Zuo et al. was based on17 CM and 20 HCs , which was a relatively small number of patients. For heterogeneity, different methods were used to analyze d-FC, such as model-dependent methods and the model-free methods.
Dynamic connectivity analysis can detect more significant networks and different connectivity pairs compared to static connectivity analysis. Dynamic connectivity analysis might be a potential method to reveal unknown migraine pathologies.
The author did all the research work of this study.
The authors have declared that no competing interests exist.
- Headache Classification Committee of the International Headache Society. The international classification of headache disorders. 3rd ed. Cephalalgia. 2018; 38: 1-211. [CrossRef]
- Schwedt TJ, Chiang CC, Chong CD, Dodick DW. Functional MRI of migraine. Lancet Neurol. 2015; 14: 81-91. [CrossRef]
- van den Heuvel MP, Hulshoff Pol HE. Exploring the brain network: A review on resting-state fMRI functional connectivity. Eur Neuropsychopharmacol. 2010; 20: 519-534. [CrossRef]
- Hutchison RM, Womelsdorf T, Allen EA, Bandettini PA, Calhoun VD, Corbetta M, et al. Dynamic functional connectivity: Promise, issues, and interpretations. Neuroimage. 2013; 80: 360-378. [CrossRef]
- Allen EA, Damaraju E, Plis SM, Erhardt EB, Eichele T, Calhoun VD. Tracking whole-brain connectivity dynamics in the resting state. Cereb Cortex. 2014; 24: 663-676. [CrossRef]
- Tu Y, Fu Z, Zeng F, Maleki N, Lan L, Li Z, et al. Abnormal thalamocortical network dynamics in migraine. Neurology. 2019; 92: e2706-e2716. [CrossRef]
- Dumkrieger G, Chong CD, Ross K, Berisha V, Schwedt TJ. Static and dynamic functional connectivity differences between migraine and persistent post-traumatic headache: A resting-state magnetic resonance imaging study. Cephalalgia. 2019; 39: 1366-1381. [CrossRef]
- Lee MJ, Park BY, Cho S, Park H, Kim ST, Chung CS. Dynamic functional connectivity of the migraine brain: A resting-state functional magnetic resonance imaging study. Pain. 2019; 160: 2776-2786. [CrossRef]
- Zou Y, Tang W, Qiao X, Li J. Aberrant modulations of static functional connectivity and dynamic functional network connectivity in chronic migraine. Quant Imaging Med Surg. 2021; 11: 2253-2264. [CrossRef]
- Wei HL, Tian T, Zhou GP, Wang JJ, Guo X, Chen YC, et al. Disrupted dynamic functional connectivity of the visual network in episodic patients with migraine without aura. Neural Plast. 2022; 2022: 9941832. [CrossRef]
- Imai N, Moriya A, Kitamura E. Functional connectivity in migraineurs with photo‐, phono‐, or osmophobia: A static and dynamic resting‐state functional magnetic resonance imaging study. Neurol Clin Neurosci. 2020; 8: 390-398. [CrossRef]
- Erhardt EB, Rachakonda S, Bedrick EJ, Allen EA, Adali T, Calhoun VD. Comparison of multi-subject ICA methods for analysis of fMRI data. Hum Brain Mapp. 2011; 32: 2075-2095. [CrossRef]
- Du Y, Fan Y. Group information guided ICA for fMRI data analysis. Neuroimage. 2013; 69: 157-197. [CrossRef]
- Beckmann CF, Smith SM. Probabilistic independent component analysis for functional magnetic resonance imaging. IEEE Trans Med Imaging. 2004; 23: 137-152. [CrossRef]
- Beckmann CF, DeLuca M, Devlin JT, Smith SM. Investigations into resting-state connectivity using independent component analysis. Philos Trans R Soc Lond B Biol Sci. 2005; 360: 1001-1013. [CrossRef]
- Jenkinson M, Beckmann CF, Behrens TE, Woolrich MW, Smith SM. FSL. Neuroimage. 2012; 62: 782-790. [CrossRef]
- Minka TP. Automatic choice of dimensionality for PCA. Proceedings of the Advances in Neural Information Processing Systems 13 (NIPS 2000); 2000 November 28-30; Denver, CO, USA. Cambridge: MIT Press.
- Coppola G, Vandenheede M, Di Clemente L, Ambrosini A, Fumal A, De Pasqua V, et al. Somatosensory evoked high-frequency oscillations reflecting thalamo-cortical activity are decreased in migraine patients between attacks. Brain. 2005; 128: 98-103. [CrossRef]
- Pietrobon D, Moskowitz MA. Pathophysiology of migraine. Annu Rev Physiol. 2013; 75: 365-391. [CrossRef]
- Noseda R, Borsook D, Burstein R. Neuropeptides and neurotransmitters that modulate thalamo-cortical pathways relevant to migraine headache. Headache. 2017; 57 Suppl 2: 97-111. [CrossRef]
- Calhoun VD, Adali T, Pearlson GD, Pekar JJ. A method for making group inferences from functional MRI data using independent component analysis. Hum Brain Mapp. 2001; 14: 140-151. [CrossRef]
- Damaraju E, Allen EA, Belger A, Ford JM, McEwen S, Mathalon DH, et al. Dynamic functional connectivity analysis reveals transient states of dysconnectivity in schizophrenia. Neuroimage Clin. 2014; 5: 298-308. [CrossRef]
- Van de Ville D, Britz J, Michel CM. EEG microstate sequences in healthy humans at rest reveal scale-free dynamics. Proc Natl Acad Sci U S A. 2010; 107: 18179-18184. [CrossRef]
- Rubinov M, Sporns O. Complex network measures of brain connectivity: Uses and interpretations. Neuroimage. 2010; 52: 1059-1069. [CrossRef]
- Wang J, Wang X, Xia M, Liao X, Evans A, He Y. Gretna: A graph theoretical network analysis toolbox for imaging connectomics. Front Hum Neurosci. 2015; 9: 386. [CrossRef]
- Keown CL, Datko MC, Chen CP, Maximo JO, Jahedi A, Müller RA. Network organization is globally atypical in autism: A graph theory study of intrinsic functional connectivity. Biol Psychiatry Cogn Neurosci Neuroimaging. 2017; 2: 66-75. [CrossRef]
- Kim J, Criaud M, Cho SS, Díez-Cirarda M, Mihaescu A, Coakeley S, et al. Abnormal intrinsic brain functional network dynamics in Parkinson’s disease. Brain. 2017; 140: 2955-2967. [CrossRef]
- Zhang J, Wang J, Wu Q, Kuang W, Huang X, He Y, et al. Disrupted brain connectivity networks in drug-naive, first-episode major depressive disorder. Biol Psychiatry. 2011; 70: 334-342. [CrossRef]
- Hashmi JA, Loggia ML, Khan S, Gao L, Kim J, Napadow V, et al. Dexmedetomidine disrupts the local and global efficiencies of large-scale brain networks. Anesthesiology. 2017; 126: 419-430. [CrossRef]
- Brennan KC, Pietrobon D. A systems neuroscience approach to migraine. Neuron. 2018; 97: 1004-1021. [CrossRef]
- Schwedt TJ, Chong CD, Peplinski J, Ross K, Berisha V. Persistent post-traumatic headache vs. Migraine: An MRI study demonstrating differences in brain structure. J Headache Pain. 2017; 18: 87. [CrossRef]
- 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-147. [CrossRef]
- Stankewitz A, Aderjan D, Eippert F, May A. Trigeminal nociceptive transmission in migraineurs predicts migraine attacks. J Neurosci. 2011; 31: 1937-1943. [CrossRef]
- Moulton EA, Becerra L, Maleki N, Pendse G, Tully S, Hargreaves R, et al. Painful heat reveals hyperexcitability of the temporal pole in interictal and ictal migraine States. Cereb Cortex. 2011; 21: 435-448. [CrossRef]
- Mickleborough MJ, Ekstrand C, Gould L, Lorentz EJ, Ellchuk T, Babyn P, et al. Attentional network differences between migraineurs and non-migraine controls: fMRI evidence. Brain Topogr. 2016; 29: 419-428. [CrossRef]
- Maniyar FH, Sprenger T, Monteith T, Schankin C, Goadsby PJ. Brain activations in the premonitory phase of nitroglycerin-triggered migraine attacks. Brain. 2014; 137: 232-241. [CrossRef]
- Li Z, Lan L, Zeng F, Makris N, Hwang J, Guo T, et al. The altered right frontoparietal network functional connectivity in migraine and the modulation effect of treatment. Cephalalgia. 2017; 37: 161-176. [CrossRef]
- Chong CD, Gaw N, Fu Y, Li J, Wu T, Schwedt TJ. Migraine classification using magnetic resonance imaging resting-state functional connectivity data. Cephalalgia. 2017; 37: 828-844. [CrossRef]
- Amin FM, Hougaard A, Magon S, Sprenger T, Wolfram F, Rostrup E, et al. Altered thalamic connectivity during spontaneous attacks of migraine without aura: A resting-state fMRI study. Cephalalgia. 2018; 38: 1237-1244. [CrossRef]
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J Royal Stat Soc Ser B. 1995; 57: 289-300. [CrossRef]
- Smith SM, Fox PT, Miller KL, Glahn DC, Fox PM, Mackay CE, et al. Correspondence of the brain's functional architecture during activation and rest. Proc Natl Acad Sci U S A. 2009; 106: 13040-13045. [CrossRef]
- Griffanti L, Douaud G, Bijsterbosch J, Evangelisti S, Alfaro-Almagro F, Glasser MF, et al. Hand classification of fMRI ICA noise components. Neuroimage. 2017; 154: 188-205. [CrossRef]
- Kelly Jr RE, Alexopoulos GS, Wang Z, Gunning FM, Murphy CF, Morimoto SS, et al. Visual inspection of independent components: Defining a procedure for artifact removal from fMRI data. J Neurosci Methods. 2010; 189: 233-245. [CrossRef]
- Bullmore E, Sporns O. Complex brain networks: Graph theoretical analysis of structural and functional systems. Nat Rev Neurosci. 2009; 10: 186-198. [CrossRef]
- Lohmann G, Margulies DS, Horstmann A, Pleger B, Lepsien J, Goldhahn D, et al. Eigenvector centrality mapping for analyzing connectivity patterns in fMRI data of the human brain. PLoS One. 2010; 5: e10232. [CrossRef]
- Zuo XN, Ehmke R, Mennes M, Imperati D, Castellanos FX, Sporns O, et al. Network centrality in the human functional connectome. Cereb Cortex. 2012; 22: 1862-1875. [CrossRef]
- Friedman J, Hastie T, Tibshirani R. Sparse inverse covariance estimation with the graphical lasso. Biostatistics. 2008; 9: 432-441. [CrossRef]
- Kiviniemi V, Vire T, Remes J, Elseoud AA, Starck T, Tervonen O, et al. A sliding time-window ICA reveals spatial variability of the default mode network in time. Brain Connect. 2011; 1: 339-347. [CrossRef]
- Androulakis XM, Krebs K, Peterlin BL, Zhang T, Maleki N, Sen S, et al. Modulation of intrinsic resting-state fMRI networks in women with chronic migraine. Neurology. 2017; 89: 163-169. [CrossRef]
- Androulakis XM, Krebs KA, Jenkins C, Maleki N, Finkel AG, Rorden C, et al. Central executive and default mode network intranet work functional connectivity patterns in chronic migraine. J Neurol Disord. 2018; 6: 393. [CrossRef]
- Li YO, Adali T, Calhoun VD. Estimating the number of independent components for functional magnetic resonance imaging data. Hum Brain Mapp. 2007; 28: 1251-1266. [CrossRef]
- Allen EA, Erhardt EB, Damaraju E, Gruner W, Segall JM, Silva RF, et al. A baseline for the multivariate comparison of resting-state networks. Front Syst Neurosci. 2011; 5: 2. [CrossRef]
- Dodick DW. Migraine. Lancet. 2018; 391: 1315-1330. [CrossRef]
- Wei HL, Zhou X, Chen YC, Yu YS, Guo X, Zhou GP, et al. Impaired intrinsic functional connectivity between the thalamus and visual cortex in migraine without aura. J Headache Pain. 2019; 20: 116. [CrossRef]
- Qin Z, Su J, He XW, Ban S, Zhu Q, Cui Y, et al. Disrupted functional connectivity between sub-regions in the sensorimotor areas and cortex in migraine without aura. J Headache Pain. 2020; 21: 47. [CrossRef]
- Shakil S, Lee CH, Keilholz SD. Evaluation of sliding window correlation performance for characterizing dynamic functional connectivity and brain states. Neuroimage. 2016; 133: 111-128. [CrossRef]
- Liu X, Chen W, Hou H, Chen X, Zhang J, Liu J, et al. Decreased functional connectivity between the dorsal anterior cingulate cortex and lingual gyrus in Alzheimer’s disease patients with depression. Behav Brain Res. 2017; 326: 132-138. [CrossRef]
- Peres MFP, Mercante JPP, Tobo PR, Kamei H, Bigal ME. Anxiety and depression symptoms and migraine: A symptom-based approach research. J Headache Pain. 2017; 18: 37. [CrossRef]
- Whitfield-Gabrieli S, Nieto-Castanon A. Conn: A functional connectivity toolbox for correlated and anticorrelated brain networks. Brain Connect. 2012; 2: 125-141. [CrossRef]
- Nieto-Castanon A. Handbook of functional connectivity Magnetic Resonance Imaging methods in CONN. Cambridge: Hilbert Press; 2020. pp. 57-60. [CrossRef]
- Whitfield-Gabrieli S, Nieto-Castanon A. CONN – functional connectivity toolbox v17 manual. Cambridge: Massachusetts Institute of Technology; 2007. Available from: https://www.nitrc.org/frs/?group_id=279.
- Podgórski P, Waliszewska-Prosół M, Zimny A, Sąsiadek M, Bladowska J. Resting-state functional connectivity of the ageing female brain-differences between young and elderly female adults on multislice short TR rs-fMRI. Front Neurol. 2021; 12: 645974. [CrossRef]
- Hougaard A, Amin FM, Magon S, Sprenger T, Rostrup E, Ashina M. No abnormalities of intrinsic brain connectivity in the interictal phase of migraine with aura. Eur J Neurol. 2015; 22: 702-e46. [CrossRef]
- Chong CD, Schwedt TJ, Hougaard A. Brain functional connectivity in headache disorders: A narrative review of MRI investigations. J Cereb Blood Flow Metab. 2019; 39: 650-669. [CrossRef]