Skip to main content

Mapping brain functional networks topological characteristics in new daily persistent headache: a magnetoencephalography study

Abstract

Background

The brain functional network topology in new daily persistent headache (NDPH) is not well understood. In this study, we aim to assess the cortical functional network topological characteristics of NDPH using non-invasive neural signal recordings.

Methods

Resting-state magnetoencephalography (MEG) was used to measure power fluctuations in neuronal oscillations from distributed cortical parcels in 35 patients with NDPH and 40 healthy controls (HCs). Their structural data were collected by 3T MRI. Functional connectivity (FC) of neural networks from 1 to 80 Hz frequency ranges was analyzed with topographic patterns and calculated network topological parameters with graph theory.

Results

In the delta (1–4 Hz) and beta (13–30 Hz) bands, the lateral occipital cortex and superior frontal gyrus FC were increased in NDPH groups compared to HCs. Graph theory analysis revealed that the NDPH had significantly increased global efficiency in the delta band and decreased nodal clustering coefficient (left medial orbitofrontal cortex) in the theta (4–8 Hz) band. The clinical characteristics had a significant correlation with network topological parameters. Age at onset of patients showed a positive correlation with global efficiency in the delta band. The degree of depression of patients showed a negative correlation with the nodal clustering coefficient (left medial orbitofrontal cortex) in the theta band.

Conclusion

The FC and topology of NDPH in brain networks may be altered, potentially leading to cortical hyperexcitability. Moreover, medial orbitofrontal cortex is involved in the pathophysiological mechanism of depression in patients with NDPH. Increased FC observed in the lateral occipital cortex and superior frontal gyrus during resting-state MEG could serve as one of the imaging characteristics associated with NDPH.

Peer Review reports

Introduction

New daily persistent headache (NDPH) was initially introduced by Vanast (1986) as a benign form of headache, characterized by persistent headache with a particular temporal profile as it starts one day with a clearly remembered onset and continues in a daily pattern without remitting [1, 2]. Although epidemiological studies on the prevalence of NDPH are not abundant, the disorder is a rare condition, with the only study meeting inclusion criteria finding a prevalence of 0.03% [3]. Persistent headache seriously affects the quality of life of patients. However, the pathogenesis of NDPH is still unclear [4].

Resting-state functional networks reflect the functional organization of the brain [5]. Magnetic resonance imaging (MRI) resting-state activity exposes significant aberrant functional connectivity (FC) in multiple brain regions among patients with NDPH [6]. Even though several alterations in brain functional connections have been pinpointed in NDPH, there is still a need for a comprehensive examination of the entire brain's connectivity network and its topological characteristics [7]. Recent developments in neuroscience emphasize that brain architecture is a combination of tightly connected networks that control diverse brain functions [8]. A deeper understanding of how network organization can be altered or disrupted during NDPH is important to clarify the pathophysiology of NDPH.

Graph theory-based complex brain network analysis provides a powerful framework to examine the topological architecture of brain networks. In graph theory, the brain is a digital network of nodes and edges. In the network, nodes represent brain regions, and edges represent FC between brain regions [9]. In neuroscience, graph theory has been widely used to understand how the segregation and integration of brain regions participate in the pathophysiological processes of diseases [10].

Magnetoencephalography (MEG), which has been used in the study of neural networks in headache, is an advanced technique that can noninvasively collect neuronal magnetic activity [11]. Compared with electroencephalography (EEG), MEG has a higher spatial resolution. Moreover, neuromagnetic signals can more directly observe neuronal activity than bold signals of functional MRI. Graph theory analysis of MEG has been widely used in neuroscience to observe the intrinsic characteristics of functional networks. Characteristics in resting state can be used as markers reflecting the progression of diseases, such as Alzheimer’s disease and Parkinson’s disease [12, 13]. Therefore, we applied resting-state MEG to study whether there were abnormalities in the neural network of patients with NDPH. In addition, we evaluated whether there was a correlation between the clinical characteristics of patients with NDPH and the network topological parameters. Our study provides clues and evidence for exploring the pathophysiological mechanism of NDPH.

Methods

Study population

From May 2020 to August 2023, 40 patients with NDPH and 43 healthy controls (HCs) were recruited from the Headache Department, Neurology Centre, Beijing Tiantan Hospital, Capital Medical University. Each recruited patient needed to be diagnosed with NDPH by two specialist neurologists. The inclusion criteria of NDPH group: (1) Satisfied the diagnostic criteria of NDPH according to ICHD-3 criteria; (2) Ages 18 to 70 years; (3) None of the patients enrolled had been prophylactically treated for NDPH for at least 3 months. The exclusion criteria of NDPH group: (1) Combined with other types of primary headache or major systemic diseases; (2) Inability to complete MEG and MRI (e.g., claustrophobia or metal implants in the body); (3) Poor data quality; (4) Pregnancy or breastfeeding. The same exclusion criteria were used for the age- and gender-matched HCs, who had no history of headache and were free of physical and psychiatric disorders. Headache information (headache history, headache frequency, etc.) and clinical scales were collected before MEG acquisition. The clinical scale included Headache Impact Test -6 (HIT-6), Patient Health Questionnaire-9 (PHQ-9), Generalized Anxiety Disorder-7 (GAD-7), Pittsburgh Sleep Quality Index (PSQI), Visual Analogue Scale (VAS) and Montreal Cognitive Assessment (MoCA). The above scales assessed the intensity of headache impact, anxiety and depression symptoms, sleep quality, pain degree, and cognitive level of patients.

The study protocol was approved by the Institutional Review Committee of Beijing Tiantan Hospital of Capital Medical University (KY2022-044), which was registered on the https://www.clinicaltrials.gov (unique identifier: NCT05334927). All participants provided informed written consent according to the Declaration of Helsinki.

MRI data acquisition

All participants were imaged with a 3.0 Tesla MR scanner (GE Healthcare, Milwaukee, WI, USA) at the Nuclear Medicine of Beijing Tiantan Hospital. Examinations were performed by a neuroradiologist who was unaware of the participant's diagnosis. Participants were asked to keep their heads and neck still, stay awake, and close their eyes, with tools to reduce noise and head movements. After checking the images, exclude images with quality problems. T1-weighted volumetric images were obtained by the 3D BRAVO sequence (coronal acquisition, the field of view (FOV) = 256 mm, acquisition matrix = 256, slice number = 192, flip angle = 15°, TR = 850 ms, TE = 320 ms, voxel size = 1 × 1 × 1.5 mm3).

MEG data acquisition

The Elekta Neuromag 306-channel scanner (Elekta TRIUX ®) was used in this study to record neural activity at 2000 Hz with a low-pass filter set to 660 Hz. The Elekta Neuromag scanner with 306 channels (102 magnetometers and 204 gradiometers) was used in this study. The head position of the participants is detected by four (head position indicator) HPI coils, and re-recorded if the head movement was excessive during the scan. Then, the head position is digitally marked (Polhemus Fastrak®). About 300 points were marked on the nasion, anterior points in front of the ear points and scalp for MRI co-registration. Data acquisition was performed using a 2000 Hz sampling rate and a low-pass filter set to 660 Hz, while the participants’ electrooculogram and electrocardiogram were recorded. Resting-state MEG data were collected for five minutes for each participant. During the scan, participants were instructed to keep their heads and neck still, stay awake, and close their eyes.

Preprocessing

After checking and excluding bad channels in the original data, the data were filtered using MaxFilter. The data sampling rate was reduced to 1000 Hz and 50 Hz line noise was removed, and finally a bandpass filter of 1-80 Hz was applied. In the independent component analysis (ICA), ocular and cardiac artifacts were marked and excluded. The cleaned data were then used to construct functional connectivity networks at the source level. Neural activity was filtered into five frequency bands: delta (1-4 Hz), theta (4-8 Hz), alpha (8-13 Hz), beta (13-30 Hz), gamma (30-80 Hz). At the source level, MRI images and MEG data were registered, and calculated the surface-based source space and inverse solution. Finally, dynamic statistical parametric mapping (dSPM) was used for source estimation [14].

Compute envelope correlations of orthogonalized activity as FC using pairwise and symmetric orthogonalization in source space. The procedure for symmetric orthogonalization in is: Extract inverse label data from raw; Symmetric orthogonalization; Band-pass filter; Hilbert transform and absolute value. According to the Desikan-Killiany Atlas, the brain was divided into 68 cortical regions (nodes). This power envelope time course was then correlated between brain region for each individual. These communication links between cortical regions (nodes) correspond to “edges” in a graph theory network model, and this specific approach for connectivity estimation shows greater repeatability than a wide range of other choices (Fig. 1) [15]. Network-based statistic (NBS) analysis was used to investigate which functional connections were significantly different between NDPH and HCs.

Fig. 1
figure 1

The pipeline of neural physiological signal to the construction of the brain network

We computed the network topological parameters with graph theory that summarize the aspects of segregation and integration of a network. In this study, we focused on network nodal and global parameters including nodal clustering coefficient, nodal efficiency, nodal degree, global efficiency, local efficiency, shortest path length of a network. The complex network analyses were performed at a sparsity range from 0.05 to 0.39 with an interval of 0.01, and the area under curve (AUC) values under this range of sparsity were calculated for both global and nodal network parameters for statistical analyses. All network analyses were performed by Gretna (http://www.nitrc.org/projects/gretna/) and visualized by using BrainNet Viewer software (http://www.nitrc.org/projects/bnv/).

Statistical analysis

The sample size was determined based on the available data and previous literature. Assuming no negative correlation between endpoints, a sample size of 75 cases (40 HC group and 35 NDPH group) would provide 90% power to reject the null hypothesis equal means at a two-sided alpha of 0.05. IBM SPSS 26.0 was used to perform the statistical analysis. Independent sample t-test and chi-square test were used to calculate the statistical difference of clinical data between groups. Continuous variables were tested for normality using the Kruskal–Wallis test; data conforming to a normal distribution were expressed as mean ± standard deviation, and otherwise as median with interquartile range. Categorical variables were expressed as numbers (percentages). Correlations between network parameters and demographic data were calculated using Pearson’s correlation, and false discovery rate (FDR) correction was performed. In NBS analysis and graph theory analysis, age and gender were used as covariables. Finally, the patients was grouped according to the headache location (unilateral and bilateral), the network topology parameters between the groups were compared. All statistical tests were two-tailed tests, with P < 0.05 indicating statistical significance.

Results

Clinical characteristics and demographics

This study included a total of 40 patients diagnosed with NDPH and 43 HCs, who were age- and gender-matched. Three patients dropped out of the study for personal reasons. Two patients and 3 HCs were excluded due to poor quality of data, segmentation errors (Fig. 2). Finally, thirty-five patients with NDPH (40.63 ± 15.65 years; 18 females) and 40 HCs (38.00 ± 13.00 years; 24 females) were included in this study. Of these patients with NDPH, 18 presented unilateral headache (51.5%). The clinical manifestations of the NDPH group were mainly nausea and vomiting (34.3%), photophobia (42.9%), and phonophobia (62.9%). The details are summarized in Table 1. All participants were right-handed and no SARS-CoV2 infection before developing headache.

Fig. 2
figure 2

Illustration of the recruitment and exclusion process

Table 1 Participants’ demographics and clinical characteristics

Abnormal brain functional networks in patients with NDPH

NBS analysis revealed increased FC in the lateral occipital cortex and superior frontal gyrus in the patient group compared to HCs. Analysis of the NBS revealed two single connected subnetwork, one with 18 nodes and 18 connections in the delta band and one with 23 nodes and 22 connections in the beta band, with the NDPH group having lower connection strength than HCs in both bands (p < 0.05). The abnormal connectivity was mainly distributed in the bilateral frontal, temporal, and occipital regions, including both long- and short-range connections (Fig. 3).

Fig. 3
figure 3

A Increased functional connections in patients with NDPH in the delta band. B Increased functional connections in NDPH patients in the delta band. In the diagram of the brain, light blue edges represent interlobar connections, and dark blue edges represent intra-lobar connections. Nodes of the frontal lobe (light blue), temporal lobe (dark blue), parietal lobe (green), orange (occipital), and cingulate cortex (brown) are marked in the connecting network. L, left; R, right; Cu, Cuneus; LO, Lateral occipital; Lg, Lingual; PerCa, Pericalcarine; CACg, Caudal anterior cingulate; Ins, Insula; IstCg, Isthmus cingulate; PoCg, Posterior cingulate; RoACg, Rostral anterior cingulate; Fu, Fusiform; IT, Inferior temporal; MT, Middle temporal; PaH, Parahippocampal; ST, Superior temporal; TrT, Transverse temporal; CMF, Caudal middle frontal; FPol, Frontal pole; MOrF, Medial orbitofrontal; PaC, Paracentral; Op, Pars opercularis; Or, Pars orbitalis; RoMF, Rostral middle frontal; SF, Superior frontal; IP, Inferior parietal; SP, Superior parietal; SM, Supra marginal

Graph theory analysis

Group comparisons revealed that the global efficiency of the functional connectivity network in patients with NDPH were greater in the delta band compared to those of controls. Compared with HCs, the clustering coefficient of the left medial orbitofrontal cortex in patients with NDPH increased in the theta band. No significant network topological parameters difference was observed in other frequency bands. There was no difference in topological parameters between the two groups of patients with unilateral headache and bilateral headache.

Correlation analysis

Clinical measures (age at onset, disease duration, VAS, MoCA, HIT-6, PHQ-9, GAD-7, PSQI, and MIDAS) did not correlate significantly with the connectivity strengths as described by a global mean of the FC. We found that the age at onset had a positive correlation with global efficiency in the delta band (p = 0.016, r = 0.495) (Fig. 4). A significant negative correlation was found between the PHQ-9 and the clustering coefficient of the left medial orbitofrontal cortex in the theta band (p = 0.016, r = -0.520) (Fig. 5). There was no statistically significant correlation between the 2 parameters (global efficiency; clustering coefficient of the left medial orbitofrontal cortex) in the graph theory of NDPH patients and other clinical measures.

Fig. 4
figure 4

A Boxplot shows the group differences in the global efficiency of networks in the delta band. The bars and error bars represent the fitted values and the standard deviations, respectively; B Correlation between global efficiency and age at onset. HC, healthy control; NDPH, new daily persistent headache

Fig. 5
figure 5

A In the theta band, the medial orbitofrontal cortex with significantly lower nodal clustering coefficient in the NDPH group compared to HC group; B Correlation between clustering coefficient of medial orbitofrontal cortex and PHQ-9 scores. HC, healthy control; NDPH, new daily persistent headache

Discussion

In this study, we innovatively used MEG to directly observe the neural magnetic signals of the participants, and observed the correlation between various cortical brain regions through the method of oscillating envelope correlation and graph theory, and analyzed the brain topological properties of the participants. Compared with HCs, patients with NDPH showed a broad increase in FC in the delta and beta bands, especially in the frontal, temporal, and occipital lobes. The distribution of FC and graph-theory analysis indicated the differences in neural networks between the NDPH and HC group in some frequency bands. Those abnormalities may have considerable mechanistic importance given the key role of neural oscillations in the long-range synchronization of information across the brain and might contribute to the development of accessible multimodal biomarkers.

Abnormal brain functional networks of NDPH

Previous studies have found that MEG resting-state activity exposes significant cortical activity changes among patients with NDPH [16]. In the delta band, our results suggest that NDPH exhibit different connective patterns between the lateral occipital cortex and the temporal lobe compared to HCs. The occipital lobe is the visual processing center containing most of the anatomical region of the visual cortex. In our previous studies, we found that patients with NDPH had altered brain perfusion and abnormal activation of the visual system [17, 18]. Changes in the visual cortex have also been observed in other headache disorders. Studies of interictal migraine have reported that abnormal FC in the occipital lobe is involved in the pathophysiological process of headache [19]. In a magnetic resonance spectroscopy (MRS) study, cortical excitability was altered in patients with migraine, and abnormal levels of multiple neurotransmitters were observed in the occipital cortex [20,21,22]. Meanwhile, migraine-like animal models suggest that the trigeminal nucleus caudalis projects indirectly to the primary visual cortex involved in the emotional/motivational aspects of pain perception and that the primary and secondary visual cortex are involved in photophobic circuits for migraine-related pain [23]. Based on the above evidence, abnormal network patterns and neurotransmitter changes in the occipital lobe may be associated with headache. In the NDPH group, 42.9% of patients were associated with photophobia. We hypothesized that the abnormal network pattern of the lateral occipital lobe may be involved in the occurrence of photophobia.

In addition, extensive increased FC of the temporal and occipital lobe was observed in this study. The temporal lobe is involved in processing sensory input into derived meanings for the appropriate retention of visual memory, language comprehension, and emotion association [24, 25]. Imaging studies of headache showed abnormal cerebral blood flow in the bilateral frontal and temporal lobes, and the cerebral blood flow in the temporal lobe was correlated with the severity of headache [26]. Similarly, in MEG studies of migraine patients, it was observed that patients had significantly enhanced connections from the frontal to temporal cortex in response to negative emotional stimuli [27]. In our study, more than half of the patients were accompanied by phonophobia (62.9%). Some studies have shown that enhanced regional functional activity of the superior temporal gyrus is presumed to be related to auditory information processing in patients with vestibular migraine. We hypothesized that the increase of FC in temporal lobe may be involved in the pathophysiology of phonophobia in NDPH [28]. In summary, our results suggest that the increase of FC in the temporal lobe may be involved in auditory information processing and emotional memory regulation in patients with NDPH.

In the beta band, the NDPH had significantly increased FC in the frontal lobe. A recently published paper has revealed that the functional connections of frontal lobe, amygdala, and insula were changed in patients with NDPH compared with HCs [29]. MEG analysis showed that patients with NDPH had structural changes and abnormal high-frequency cortical activity in both the frontal and temporal lobes [16]. The frontal lobe involves pain regulation, cognitive control, and executive function, and thus is a primary candidate for dysfunction in many neurodevelopmental and neuropsychiatric disorders [30, 31]. Increased FC in frontal lobe may support the theory of central hyperexcitability in NDPH. These findings may facilitate to develop new therapeutic strategies for migraine. For example, in the following studies, repetitive transcranial electrical stimulation (rTMS) can be applied to patients with NDPH to inhibit frontal neuron excitation through low-frequency rTMS, and clinical efficacy can be observed to determine the role of frontal lobe in the pathophysiology of NDPH [32, 33].

Topological characteristics of NDPH

In this study, we used graph theory for the first time to analyze the functional network topological properties of NDPH patients. Compared to HCs, the global efficiency of the functional connectivity network in patients with NDPH was increased in the delta band. The global efficiency mainly characterize the global information transmission capacity of brain network, and higher global efficiency of a network reflect faster information transfer between network nodes [34]. Different levels of visual perceptual skills are associated with specific modifications in global efficiency, and higher global efficiency may imply increased hemispherical asymmetry [35]. Thus, these findings about abnormal global efficiency may support the speculation that the brain networks in patients with NDPH tend to be randomized with increased global integration to obscure the role of key nodes and reduce the modularity of the network.

Our study showed that the clustering coefficient of the left medial orbitofrontal cortex in patients with NDPH increased in the theta band. The clustering coefficient of a node represents the ratio of all existing connections between the “neighbors” of a node (nodes that are one-step away) and the maximum possible number of edges between the neighbors. The orbitofrontal cortex was defined as the part of the prefrontal cortex that receives projections from the medial dorsal nucleus of the thalamus [36,37,38]. Functional and structural prefrontal cortex abnormalities have been reported in major depressive disorder. Patients with depression have been linked to decreased activation in ventral and dorsal striatum as well as middle frontal gyrus, but hyperactivation in the prefrontal cortex [39, 40]. The mid-anterior orbitofrontal cortex has been found to consistently track subjective pleasure in neuroimaging studies [41]. We observed that patients with NDPH generally had mild to moderate depression, which may be related to decreased local integration of the medial orbitofrontal cortex.

Clinical correlation analysis of NDPH

To identify the clinical significance of aberrant MEG parameters in NDPH, we further calculated the correlations between the aberrant MEG parameters and clinical characteristics in NDPH. We found that the age at onset had a positive correlation with global efficiency in the delta band. And a significant negative correlation was found between the PHQ-9 and the clustering coefficient of the left medial orbitofrontal cortex in the theta band. By studying the developmental changes in functional brain networks from birth to adolescence, Gozdas found that the global efficiency increases with age [42]. This may explain why the greater the age of onset in patients with NDPH, the greater the global efficiency. Functional and structural changes in the medial orbitofrontal cortex, a brain region closely related to emotion, have been reported in depressed patients [43]. Our results suggest that the medial orbitofrontal cortex may play an important role in the development of depression in patients with NDPH. It is worth noting that our previous study also found that the cortical thickness of the left rostral middle frontal gyrus was reduced, so it is necessary to verify the important role of the frontal lobe, especially the prefrontal cortex, in patients with NDPH [16]. There were no correlations between the aberrant MEG parameters and other clinical characteristics in NDPH.

Study limitations and future directions

We noted that this study has some limitations. The first limitation of this study is the small number of participants, and our results still need to be validated in more data. Second, this study lacked longitudinal data of participants to observe neuromagnetic source activity during NDPH. Finally, our study found a close relationship between the degree of depression in NDPH and the medial orbitofrontal cortex. Therefore, it is necessary to distinguish various comorbidity subgroups of NDPH further and conduct clinical intervention to seek new ways to treat headache in the future.

Conclusions

This study suggests that the FC and topology of NDPH in brain networks may be altered and the cortex is overexcited. The increased FC in lateral occipital cortex and superior frontal gyrus under resting-state MEG may be one of the imaging features of NDPH. In addition, the frontal lobe (especially the medial orbitofrontal cortex) may play an important role in NDPH and is closely related to the degree of depression in patients with NDPH. These findings may contribute to developing new treatment strategies for NDPH, which can be further clinically validated by non-invasive neuroregulatory techniques (such as rTMS) in the future.

Availability of data and materials

The data that support the findings of this study are available from the corresponding author, upon reasonable request.

Abbreviations

AUC:

Area under curve

dSPM:

Dynamic statistical parametric mapping

NDPH:

New daily persistent headache

EEG:

Electroencephalography

fMRI:

Functional magnetic resonance imaging

FC:

Functional connectivity

FOV:

the field of view

FDR:

False discovery rate

GAD-7:

Generalized Anxiety Disorder-7

HCs:

Healthy controls

HIT-6:

Headache Impact Test-6

HPI:

Head position indicator

ICA:

Independent component analysis

MEG:

Magnetoencephalography

MRI:

Magnetic resonance imaging

MoCA:

Montreal Cognitive Assessment

MRS:

Magnetic resonance spectroscopy

NBS:

Network-based statistic

PHQ-9:

Patient Health Questionnaire-9

PSQI:

Pittsburgh Sleep Quality Index

rTMS:

Repetitive transcranial electrical stimulation

VAS:

Visual Analogue Scale

References

  1. Vanast W (1986) New daily persistent headache: definition of a benign syndrome. Headache 26:317

    Google Scholar 

  2. Headache Classification Committee of the International Headache Society (IHS) (2018) The international classification of headache disorders, 3rd edition. Cephalalgia 38(1). https://doi.org/10.1177/0333102417738202

  3. Grande RB, Aaseth K, Lundqvist C, Russell MB (2009) Prevalence of new daily persistent headache in the general population. The Akershus study of chronic headache. Cephalalgia 29(11):1149–55. https://doi.org/10.1111/j.1468-2982.2009.01842.x

    Article  PubMed  CAS  Google Scholar 

  4. Li D, Rozen TD (2002) The clinical characteristics of new daily persistent headache. Cephalalgia 22(1):66–69

    Article  PubMed  CAS  Google Scholar 

  5. Raichle ME (2015) The brain’s default mode network. Annu Rev Neurosci 38:433–447. https://doi.org/10.1146/annurev-neuro-071013-014030

    Article  PubMed  CAS  Google Scholar 

  6. Wang W, Yuan Z, Zhang X et al (2023) Mapping the aberrant brain functional connectivity in new daily persistent headache: a resting-state functional magnetic resonance imaging study. J Headache Pain. 24(1):46. https://doi.org/10.1186/s10194-023-01577-2

    Article  PubMed  PubMed Central  Google Scholar 

  7. Szabo E, Chang YC, Shulman J et al (2022) Alterations in the structure and function of the brain in adolescents with new daily persistent headache: a pilot MRI study. Headache 62(7):858–869. https://doi.org/10.1111/head.14360

    Article  PubMed  Google Scholar 

  8. Hallett M, de Haan W, Deco G et al (2020) Human brain connectivity: clinical applications for clinical neurophysiology. Clin Neurophysiol 131(7):1621–1651. https://doi.org/10.1016/j.clinph.2020.03.031

    Article  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  10. Proudfoot M, Colclough GL, Quinn A et al (2018) Increased cerebral functional connectivity in ALS: a resting-state magnetoencephalography study. Neurology 90(16):e1418–e1424. https://doi.org/10.1212/wnl.0000000000005333

    Article  PubMed  PubMed Central  Google Scholar 

  11. Singh SP (2014) Magnetoencephalography: Basic principles. Ann Indian Acad Neurol 17(Suppl 1):S107–S112. https://doi.org/10.4103/0972-2327.128676

    Article  PubMed  PubMed Central  Google Scholar 

  12. Stam CJ, de Haan W, Daffertshofer A et al (2009) Graph theoretical analysis of magnetoencephalographic functional connectivity in Alzheimer’s disease. Brain 132(Pt 1):213–224. https://doi.org/10.1093/brain/awn262

    Article  PubMed  CAS  Google Scholar 

  13. Olde Dubbelink KT, Hillebrand A, Stoffers D et al (2014) Disrupted brain network topology in Parkinson’s disease: a longitudinal magnetoencephalography study. Brain 137(Pt 1):197–207. https://doi.org/10.1093/brain/awt316

    Article  PubMed  Google Scholar 

  14. Dale AM, Liu AK, Fischl BR et al (2000) Dynamic statistical parametric mapping: combining fMRI and MEG for high-resolution imaging of cortical activity. Neuron. 26(1):55–67

    Article  PubMed  CAS  Google Scholar 

  15. Colclough GL, Woolrich MW, Tewarie PK, Brookes MJ, Quinn AJ, Smith SM (2016) How reliable are MEG resting-state connectivity metrics? Neuroimage 138:284–293. https://doi.org/10.1016/j.neuroimage.2016.05.070

    Article  PubMed  CAS  Google Scholar 

  16. Qiu D, Wang W, Mei Y et al (2023) Brain structure and cortical activity changes of new daily persistent headache: multimodal evidence from MEG/sMRI. J Headache Pain. 24(1):45. https://doi.org/10.1186/s10194-023-01581-6

    Article  PubMed  PubMed Central  Google Scholar 

  17. Zhang X, Wang W, Bai X et al (2023) Alterations in regional homogeneity and multiple frequency amplitudes of low-frequency fluctuation in patients with new daily persistent headache: a resting-state functional magnetic resonance imaging study. J Headache Pain 24(1):14. https://doi.org/10.1186/s10194-023-01543-y

    Article  PubMed  PubMed Central  Google Scholar 

  18. Bai X, Wang W, Zhang X et al (2022) Cerebral perfusion variance in new daily persistent headache and chronic migraine: an arterial spin-labeled MR imaging study. J Headache Pain 23(1):156. https://doi.org/10.1186/s10194-022-01532-7

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Skorobogatykh K, van Hoogstraten WS, Degan D et al (2019) Functional connectivity studies in migraine: what have we learned? J Headache Pain 20(1):108. https://doi.org/10.1186/s10194-019-1047-3

    Article  PubMed  PubMed Central  Google Scholar 

  20. Bridge H, Stagg CJ, Near J, Lau C-i, Zisner A, Cader MZ (2015) Altered neurochemical coupling in the occipital cortex in migraine with visual aura. Cephalalgia. 35(11):1025–1030. https://doi.org/10.1177/0333102414566860

    Article  PubMed  Google Scholar 

  21. Bigal ME, Hetherington H, Pan J et al (2008) Occipital levels of GABA are related to severe headaches in migraine. Neurology 70(22):2078–2080. https://doi.org/10.1212/01.wnl.0000313376.07248.28

    Article  PubMed  CAS  Google Scholar 

  22. Zhang L, Huang J, Zhang Z, Cao Z (2021) Altered metabolites in the occipital lobe in migraine without aura during the attack and the interictal period. Front Neurol 12:656349. https://doi.org/10.3389/fneur.2021.656349

    Article  PubMed  PubMed Central  Google Scholar 

  23. Brennan KC, Pietrobon D (2018) A systems neuroscience approach to migraine. Neuron 97(5):1004–1021. https://doi.org/10.1016/j.neuron.2018.01.029

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Nyberg L, McIntosh AR, Houle S, Nilsson LG, Tulving E (1996) Activation of medial temporal structures during episodic memory retrieval. Nature 380(6576):715–717

    Article  PubMed  CAS  Google Scholar 

  25. Schacter DL (1998) Memory and awareness. Science (New York, NY) 280(5360):59–60

    Article  CAS  Google Scholar 

  26. Cheng M-H, Wen S-L, Zhou HJ, Lian-Fang B, Li J-F, Xie L-J (2013) Evaluation of headache and regional cerebral flood flow in patients with migraine. Clin Nucl Med 38(11):874–877. https://doi.org/10.1097/RLU.0b013e3182a75927

    Article  PubMed  Google Scholar 

  27. Ren J, Yao Q, Tian M et al (2022) Altered effective connectivity in migraine patients during emotional stimuli: a multi-frequency magnetoencephalography study. J Headache Pain 23(1):6. https://doi.org/10.1186/s10194-021-01379-4

    Article  PubMed  PubMed Central  Google Scholar 

  28. Li Z-Y, Si L-H, Shen B, Ling X, Yang X (2023) Altered functional activity in the right superior temporal gyrus in patients with definite vestibular migraine. Neurol Sci 44(5):1719–1729. https://doi.org/10.1007/s10072-022-06570-2

    Article  PubMed  Google Scholar 

  29. Szabo E, Chang YH-C, Shulman J et al (2022) Alterations in the structure and function of the brain in adolescents with new daily persistent headache: a pilot MRI study. Headache 62(7):858–869. https://doi.org/10.1111/head.14360

    Article  PubMed  Google Scholar 

  30. Cieslik EC, Zilles K, Caspers S et al (2013) Is there “one” DLPFC in cognitive action control? Evidence for heterogeneity from co-activation-based parcellation. Cerebral Cortex (New York, NY: 1991). 23(11):2677–89. https://doi.org/10.1093/cercor/bhs256

    Article  Google Scholar 

  31. Seminowicz DA, Moayedi M (2017) The dorsolateral prefrontal cortex in acute and chronic pain. J Pain 18(9):1027–1035. https://doi.org/10.1016/j.jpain.2017.03.008

    Article  PubMed  PubMed Central  Google Scholar 

  32. Bestmann S, Walsh V (2017) Transcranial electrical stimulation. Curr Biol 27(23):R1258–R1262. https://doi.org/10.1016/j.cub.2017.11.001

    Article  PubMed  CAS  Google Scholar 

  33. Teepker M, Hötzel J, Timmesfeld N et al (2010) Low-frequency rTMS of the vertex in the prophylactic treatment of migraine. Cephalalgia 30(2):137–144. https://doi.org/10.1111/j.1468-2982.2009.01911.x

    Article  PubMed  CAS  Google Scholar 

  34. Zhang W, Guo L, Liu D, Xu G (2020) The dynamic properties of a brain network during working memory based on the algorithm of cross-frequency coupling. Cogn Neurodyn 14(2):215–228. https://doi.org/10.1007/s11571-019-09562-9

    Article  PubMed  Google Scholar 

  35. Danti S, Handjaras G, Cecchetti L, Beuzeron-Mangina H, Pietrini P, Ricciardi E (2018) Different levels of visual perceptual skills are associated with specific modifications in functional connectivity and global efficiency. Int J Psychophysiol 123:127–135. https://doi.org/10.1016/j.ijpsycho.2017.10.002

    Article  PubMed  Google Scholar 

  36. Rolls ET (2004) The functions of the orbitofrontal cortex. Brain Cogn 55(1):11–29

    Article  PubMed  Google Scholar 

  37. Kringelbach ML (2005) The human orbitofrontal cortex: linking reward to hedonic experience. Nat Rev Neurosci 6(9):691–702

    Article  PubMed  CAS  Google Scholar 

  38. Rushworth MFS, Behrens TEJ, Rudebeck PH, Walton ME (2007) Contrasting roles for cingulate and orbitofrontal cortex in decisions and social behaviour. Trends Cogn Sci 11(4):168–176

    Article  PubMed  CAS  Google Scholar 

  39. Borsini A, Wallis ASJ, Zunszain P, Pariante CM, Kempton MJ (2020) Characterizing anhedonia: a systematic review of neuroimaging across the subtypes of reward processing deficits in depression. Cogn Affect Behav Neurosci 20(4):816–841. https://doi.org/10.3758/s13415-020-00804-6

    Article  PubMed  PubMed Central  Google Scholar 

  40. Zhang W-N, Chang S-H, Guo L-Y, Zhang K-L, Wang J (2013) The neural correlates of reward-related processing in major depressive disorder: a meta-analysis of functional magnetic resonance imaging studies. J Affect Disord 151(2):531–539. https://doi.org/10.1016/j.jad.2013.06.039

    Article  PubMed  Google Scholar 

  41. Berridge KC, Kringelbach ML (2015) Pleasure systems in the brain. Neuron 86(3):646–664. https://doi.org/10.1016/j.neuron.2015.02.018

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Gozdas E, Holland SK, Altaye M (2019) Developmental changes in functional brain networks from birth through adolescence. Hum Brain Mapp 40(5):1434–1444. https://doi.org/10.1002/hbm.24457

    Article  PubMed  Google Scholar 

  43. Pizzagalli DA, Roberts AC (2022) Prefrontal cortex and depression. Neuropsychopharmacology 47(1):225–246. https://doi.org/10.1038/s41386-021-01101-7

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the participants and medical staff involved in this study, and Dr. Bingjiang Lyu for helpful discussion.

Funding

This study was supported by the National Natural Science Foundation of China (grant numbers: 32170752, 91849104, and 31770800), and the National Natural Science Foundation of Beijing (grant number: Z200024).

Author information

Authors and Affiliations

Authors

Contributions

DQ, ZLG, and YGW supported the conception and design of this project. DQ, ZLG, YGW, YLM, WW, HFT, PZ, XSL, XYY, XL and ZHX acquired data. DQ analyzed the data. ZYY, MTZ and YQZ contributed to data quality control. DQ produced the first draft. All authors contributed intellectual content to the revised manuscript and have read and approved the final manuscript.

Corresponding author

Correspondence to Yonggang Wang.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the Institutional Review Committee of Beijing Tiantan Hospital of Capital Medical University (KY2022-044), which was registered on the https://www.clinicaltrials.gov (unique identifier: NCT05334927). All participants provided informed written consent according to the Declaration of Helsinki.

Consent for publication

All authors consent for the publication.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Qiu, D., Ge, Z., Mei, Y. et al. Mapping brain functional networks topological characteristics in new daily persistent headache: a magnetoencephalography study. J Headache Pain 24, 161 (2023). https://doi.org/10.1186/s10194-023-01695-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s10194-023-01695-x

Keywords